1
|
Hulikal Chakrapani T, Hajibeygi H, Moultos OA, Vlugt TJH. Mutual Diffusivities of Mixtures of Carbon Dioxide and Hydrogen and Their Solubilities in Brine: Insight from Molecular Simulations. Ind Eng Chem Res 2024; 63:10456-10481. [PMID: 38882502 PMCID: PMC11177264 DOI: 10.1021/acs.iecr.4c01078] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2024] [Revised: 05/14/2024] [Accepted: 05/20/2024] [Indexed: 06/18/2024]
Abstract
H2-CO2 mixtures find wide-ranging applications, including their growing significance as synthetic fuels in the transportation industry, relevance in capture technologies for carbon capture and storage, occurrence in subsurface storage of hydrogen, and hydrogenation of carbon dioxide to form hydrocarbons and alcohols. Here, we focus on the thermodynamic properties of H2-CO2 mixtures pertinent to underground hydrogen storage in depleted gas reservoirs. Molecular dynamics simulations are used to compute mutual (Fick) diffusivities for a wide range of pressures (5 to 50 MPa), temperatures (323.15 to 423.15 K), and mixture compositions (hydrogen mole fraction from 0 to 1). At 5 MPa, the computed mutual diffusivities agree within 5% with the kinetic theory of Chapman and Enskog at 423.15 K, albeit exhibiting deviations of up to 25% between 323.15 and 373.15 K. Even at 50 MPa, kinetic theory predictions match computed diffusivities within 15% for mixtures comprising over 80% H2 due to the ideal-gas-like behavior. In mixtures with higher concentrations of CO2, the Moggridge correlation emerges as a dependable substitute for the kinetic theory. Specifically, when the CO2 content reaches 50%, the Moggridge correlation achieves predictions within 10% of the computed Fick diffusivities. Phase equilibria of ternary mixtures involving CO2-H2-NaCl were explored using Gibbs Ensemble (GE) simulations with the Continuous Fractional Component Monte Carlo (CFCMC) technique. The computed solubilities of CO2 and H2 in NaCl brine increased with the fugacity of the respective component but decreased with NaCl concentration (salting out effect). While the solubility of CO2 in NaCl brine decreased in the ternary system compared to the binary CO2-NaCl brine system, the solubility of H2 in NaCl brine increased less in the ternary system compared to the binary H2-NaCl brine system. The cooperative effect of H2-CO2 enhances the H2 solubility while suppressing the CO2 solubility. The water content in the gas phase was found to be intermediate between H2-NaCl brine and CO2-NaCl brine systems. Our findings have implications for hydrogen storage and chemical technologies dealing with CO2-H2 mixtures, particularly where experimental data are lacking, emphasizing the need for reliable thermodynamic data on H2-CO2 mixtures.
Collapse
Affiliation(s)
- Thejas Hulikal Chakrapani
- Reservoir Engineering, Geoscience and Engineering Department, Faculty of Civil Engineering and Geosciences, Delft University of Technology, Delft 2628 CN, The Netherlands
| | - Hadi Hajibeygi
- Reservoir Engineering, Geoscience and Engineering Department, Faculty of Civil Engineering and Geosciences, Delft University of Technology, Delft 2628 CN, The Netherlands
| | - Othonas A Moultos
- Engineering Thermodynamics, Process and Energy Department, Faculty of Mechanical Engineering, Delft University of Technology, Delft 2628 CB, The Netherlands
| | - Thijs J H Vlugt
- Engineering Thermodynamics, Process and Energy Department, Faculty of Mechanical Engineering, Delft University of Technology, Delft 2628 CB, The Netherlands
| |
Collapse
|
2
|
van Rooijen WA, Habibi P, Xu K, Dey P, Vlugt TJH, Hajibeygi H, Moultos OA. Interfacial Tensions, Solubilities, and Transport Properties of the H 2/H 2O/NaCl System: A Molecular Simulation Study. JOURNAL OF CHEMICAL AND ENGINEERING DATA 2024; 69:307-319. [PMID: 38352074 PMCID: PMC10859954 DOI: 10.1021/acs.jced.2c00707] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/12/2022] [Accepted: 12/23/2022] [Indexed: 02/16/2024]
Abstract
Data for several key thermodynamic and transport properties needed for technologies using hydrogen (H2), such as underground H2 storage and H2O electrolysis are scarce or completely missing. Force field-based Molecular Dynamics (MD) and Continuous Fractional Component Monte Carlo (CFCMC) simulations are carried out in this work to cover this gap. Extensive new data sets are provided for (a) interfacial tensions of H2 gas in contact with aqueous NaCl solutions for temperatures of (298 to 523) K, pressures of (1 to 600) bar, and molalities of (0 to 6) mol NaCl/kg H2O, (b) self-diffusivities of infinitely diluted H2 in aqueous NaCl solutions for temperatures of (298 to 723) K, pressures of (1 to 1000) bar, and molalities of (0 to 6) mol NaCl/kg H2O, and (c) solubilities of H2 in aqueous NaCl solutions for temperatures of (298 to 363) K, pressures of (1 to 1000) bar, and molalities of (0 to 6) mol NaCl/kg H2O. The force fields used are the TIP4P/2005 for H2O, the Madrid-2019 and the Madrid-Transport for NaCl, and the Vrabec and Marx for H2. Excellent agreement between the simulation results and available experimental data is found with average deviations lower than 10%.
Collapse
Affiliation(s)
- W. A. van Rooijen
- Reservoir
Engineering, Geoscience and Engineering Department, Faculty of Civil
Engineering and Geosciences, Delft University
of Technology, Stevinweg 1, 2628CN, Delft, The Netherlands
| | - P. Habibi
- Engineering
Thermodynamics, Process and Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628CB, Delft, The Netherlands
- Department
of Materials Science and Engineering, Faculty of Mechanical, Maritime
and Materials Engineering, Delft University
of Technology, Mekelweg
2, 2628CD, Delft, The Netherlands
| | - K. Xu
- Department
of Materials Science and Engineering, Faculty of Mechanical, Maritime
and Materials Engineering, Delft University
of Technology, Mekelweg
2, 2628CD, Delft, The Netherlands
| | - P. Dey
- Department
of Materials Science and Engineering, Faculty of Mechanical, Maritime
and Materials Engineering, Delft University
of Technology, Mekelweg
2, 2628CD, Delft, The Netherlands
| | - T. J. H. Vlugt
- Engineering
Thermodynamics, Process and Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628CB, Delft, The Netherlands
| | - H. Hajibeygi
- Reservoir
Engineering, Geoscience and Engineering Department, Faculty of Civil
Engineering and Geosciences, Delft University
of Technology, Stevinweg 1, 2628CN, Delft, The Netherlands
| | - O. A. Moultos
- Engineering
Thermodynamics, Process and Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628CB, Delft, The Netherlands
| |
Collapse
|
3
|
Shi R, Cooper AJ, Tanaka H. Impact of hierarchical water dipole orderings on the dynamics of aqueous salt solutions. Nat Commun 2023; 14:4616. [PMID: 37550299 PMCID: PMC10406952 DOI: 10.1038/s41467-023-40278-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2022] [Accepted: 07/13/2023] [Indexed: 08/09/2023] Open
Abstract
Ions exhibit highly ion-specific complex behaviours when solvated in water, which remains a mystery despite the fundamental importance of ion solvation in nature, science, and technology. Here we explain these ion-specific properties by the ion-induced hierarchical dipolar, translational, and bond-orientational orderings of ion hydration shell under the competition between ion-water electrostatic interactions and inter-water hydrogen bonding. We first characterise this competition by a new length λHB(q), explaining the ion-specific effects on solution dynamics. Then, by continuously tuning ion size and charge, we find that the bond-orientational order of the ion hydration shell highly develops for specific ion size and charge combinations. This ordering drastically stabilises the hydration shell; its degree changes the water residence time around ions by 11 orders of magnitude for main-group ions. These findings are fundamental to ionic processes in aqueous solutions, providing a physical principle for electrolyte design and application.
Collapse
Affiliation(s)
- Rui Shi
- Zhejiang Province Key Laboratory of Quantum Technology and Device, School of Physics, Zhejiang University, Hangzhou, 310027, China.
- Department of Fundamental Engineering, Institute of Industrial Science, The University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo, 153-8505, Japan.
| | - Anthony J Cooper
- Department of Fundamental Engineering, Institute of Industrial Science, The University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo, 153-8505, Japan
- Department of Physics, University of California, Santa Barbara, CA, 93106-9530, USA
| | - Hajime Tanaka
- Department of Fundamental Engineering, Institute of Industrial Science, The University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo, 153-8505, Japan.
- Research Center for Advanced Science and Technology, The University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo, 153-8904, Japan.
| |
Collapse
|
4
|
Yang Y, Che Ruslan MFA, Narayanan Nair AK, Qiao R, Sun S. Interfacial properties of the hexane + carbon dioxide + water system in the presence of hydrophilic silica. J Chem Phys 2022; 157:234704. [PMID: 36550045 DOI: 10.1063/5.0130986] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022] Open
Abstract
Molecular dynamics simulations were conducted to study the interfacial behavior of the CO2 + H2O and hexane + CO2 + H2O systems in the presence of hydrophilic silica at geological conditions. Simulation results for the CO2 + H2O and hexane + CO2 + H2O systems are in reasonable agreement with the theoretical predictions based on the density functional theory. In general, the interfacial tension (IFT) of the CO2 + H2O system exponentially (linearly) decreased with increasing pressure (temperature). The IFTs of the hexane + CO2 + H2O (two-phase) system decreased with the increasing mole fraction of CO2 in the hexane/CO2-rich phase xCO2 . Here, the negative surface excesses of hexane lead to a general increase in the IFTs with increasing pressure. The effect of pressure on these IFTs decreased with increasing xCO2 due to the positive surface excesses of carbon dioxide. The simulated water contact angles of the CO2 + H2O + silica system fall in the range from 43.8° to 76.0°, which is in reasonable agreement with the experimental results. These contact angles increased with pressure and decreased with temperature. Here, the adhesion tensions are influenced by the variations in fluid-fluid IFT and contact angle. The simulated water contact angles of the hexane + H2O + silica system fall in the range from 58.0° to 77.0° and are not much affected by the addition of CO2. These contact angles increased with pressure, and the pressure effect was less pronounced at lower temperatures. Here, the adhesion tensions are mostly influenced by variations in the fluid-fluid IFTs. In all studied cases, CO2 molecules could penetrate into the interfacial region between the water droplet and the silica surface.
Collapse
Affiliation(s)
- Yafan Yang
- State Key Laboratory for Geomechanics and Deep Underground Engineering, China University of Mining and Technology, Xuzhou, Jiangsu 221116, China
| | - Mohd Fuad Anwari Che Ruslan
- Physical Science and Engineering Division (PSE), Computational Transport Phenomena Laboratory, King Abdullah University of Science and Technology (KAUST), Thuwal 23955-6900, Saudi Arabia
| | - Arun Kumar Narayanan Nair
- Physical Science and Engineering Division (PSE), Computational Transport Phenomena Laboratory, King Abdullah University of Science and Technology (KAUST), Thuwal 23955-6900, Saudi Arabia
| | - Rui Qiao
- Department of Mechanical Engineering, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Shuyu Sun
- Physical Science and Engineering Division (PSE), Computational Transport Phenomena Laboratory, King Abdullah University of Science and Technology (KAUST), Thuwal 23955-6900, Saudi Arabia
| |
Collapse
|
5
|
Hosseinzadeh Dehaghani Y, Assareh M, Feyzi F. Simultaneous Prediction of Equilibrium, Interfacial, and Transport Properties of CO 2-Brine Systems Using Molecular Dynamics Simulation: Applications to CO 2 Storage. Ind Eng Chem Res 2022. [DOI: 10.1021/acs.iecr.2c02249] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Yasaman Hosseinzadeh Dehaghani
- Thermodynamics Research Laboratory, School of Chemical, Petroleum and Gas Engineering, Iran University of Science and Technology (IUST), Tehran 16846-13114, Iran
| | - Mehdi Assareh
- Thermodynamics Research Laboratory, School of Chemical, Petroleum and Gas Engineering, Iran University of Science and Technology (IUST), Tehran 16846-13114, Iran
| | - Farzaneh Feyzi
- Thermodynamics Research Laboratory, School of Chemical, Petroleum and Gas Engineering, Iran University of Science and Technology (IUST), Tehran 16846-13114, Iran
| |
Collapse
|
6
|
Narayanan Nair AK, Anwari Che Ruslan MF, Ramirez Hincapie ML, Sun S. Bulk and Interfacial Properties of Brine or Alkane in the Presence of Carbon Dioxide, Methane, and Their Mixture. Ind Eng Chem Res 2022. [DOI: 10.1021/acs.iecr.2c00249] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Arun Kumar Narayanan Nair
- Physical Science and Engineering Division (PSE), Computational Transport Phenomena Laboratory, King Abdullah University of Science and Technology (KAUST), Thuwal, 23955-6900, Saudi Arabia
| | - Mohd Fuad Anwari Che Ruslan
- Physical Science and Engineering Division (PSE), Computational Transport Phenomena Laboratory, King Abdullah University of Science and Technology (KAUST), Thuwal, 23955-6900, Saudi Arabia
| | - Marcia Luna Ramirez Hincapie
- Physical Science and Engineering Division (PSE), Computational Transport Phenomena Laboratory, King Abdullah University of Science and Technology (KAUST), Thuwal, 23955-6900, Saudi Arabia
| | - Shuyu Sun
- Physical Science and Engineering Division (PSE), Computational Transport Phenomena Laboratory, King Abdullah University of Science and Technology (KAUST), Thuwal, 23955-6900, Saudi Arabia
| |
Collapse
|
7
|
Ofori K, Phan CM, Barifcani A, Iglauer S. An investigation of some H2S thermodynamical properties at the water interface under pressurised conditions through molecular dynamics. Mol Phys 2021. [DOI: 10.1080/00268976.2021.2011972] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Affiliation(s)
- Kofi Ofori
- Department of Petroleum Engineering, Curtin University, Perth, Australia
| | - Chi M. Phan
- Department of Chemical Engineering, Curtin University, Perth, Australia
| | - Ahmed Barifcani
- Department of Petroleum Engineering, Curtin University, Perth, Australia
| | - Stefan Iglauer
- Department of Petroleum Engineering, Curtin University, Perth, Australia
- School of Engineering, Edith Cowan University, Joondalup, Australia
| |
Collapse
|
8
|
Muniz MC, Gartner TE, Riera M, Knight C, Yue S, Paesani F, Panagiotopoulos AZ. Vapor-liquid equilibrium of water with the MB-pol many-body potential. J Chem Phys 2021; 154:211103. [PMID: 34240989 DOI: 10.1063/5.0050068] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Abstract
Among the many existing molecular models of water, the MB-pol many-body potential has emerged as a remarkably accurate model, capable of reproducing thermodynamic, structural, and dynamic properties across water's solid, liquid, and vapor phases. In this work, we assessed the performance of MB-pol with respect to an important set of properties related to vapor-liquid coexistence and interfacial behavior. Through direct coexistence classical molecular dynamics simulations at temperatures of 400 K < T < 600 K, we calculated properties such as equilibrium coexistence densities, vapor-liquid interfacial tension, vapor pressure, and enthalpy of vaporization and compared the MB-pol results to experimental data. We also compared rigid vs fully flexible variants of the MB-pol model and evaluated system size effects for the properties studied. We found that the MB-pol model predictions are in good agreement with experimental data, even for temperatures approaching the vapor-liquid critical point; this agreement was largely insensitive to system sizes or the rigid vs flexible treatment of the intramolecular degrees of freedom. These results attest to the chemical accuracy of MB-pol and its high degree of transferability, thus enabling MB-pol's application across a large swath of water's phase diagram.
Collapse
Affiliation(s)
- Maria Carolina Muniz
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA
| | - Thomas E Gartner
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, USA
| | - Marc Riera
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, USA
| | - Christopher Knight
- Computational Science Division, Argonne National Laboratory, Argonne, Illinois 60439, USA
| | - Shuwen Yue
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, USA
| | | |
Collapse
|
9
|
Interaction of low salinity surfactant nanofluids with carbonate surfaces and molecular level dynamics at fluid-fluid interface at ScCO 2 loading. J Colloid Interface Sci 2020; 586:315-325. [PMID: 33148450 DOI: 10.1016/j.jcis.2020.10.095] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/20/2020] [Revised: 10/12/2020] [Accepted: 10/23/2020] [Indexed: 10/23/2022]
Abstract
HYPOTHESIS The advanced low salinity aqueous formulations are yet to be validated as an injection fluid for enhanced oil recovery (EOR) from the carbonate reservoirs and CO2 geosequestration. Interaction of various ionic species present in the novel low salinity surfactant nanofluids with scCO2/CO2 saturated aqueous phase interface and at the interface of CO2 saturated aqueous phase/mixed wet (with CO2 and Decane) limestone surface at the conditions of low salinity at reservoir conditions are to yet to be understood. EXPERIMENTS This study, carried out for the first time in low salinity at scCO2 loading conditions at 20 MPa pressure and 343 K temperature, comprises of wettability study of the limestone surface by aqueous phase contact angle measurements using ZrO2 nanoparticles (in the concentration range of 100-2000 mg/L) and 0.82 mM Hexadecyltrimethylammonium bromide (CTAB) surfactant. Molecular dynamics simulations results were used to understand the underlying mechanism of wettability alteration and interfacial tension (IFT) change. FINDINGS This study reveals that a low dosage (100 mg/L) of ZrO2 nanoparticles forming ZrO2-CTAB nano-complexes helps in wettability alteration of the rock surface to more water-wetting state; certain ionic species augment this effect when used in appropriate concentration. Also, these nano-complexes helps in scCO2/CO2 saturated aqueous phase IFT reduction. This study can be used to design advanced low salinity injection fluids for water alternating gas injection for EOR and CO2 geosequestration projects.
Collapse
|
10
|
Ansari N, Karmakar T, Parrinello M. Molecular Mechanism of Gas Solubility in Liquid: Constant Chemical Potential Molecular Dynamics Simulations. J Chem Theory Comput 2020; 16:5279-5286. [PMID: 32551636 DOI: 10.1021/acs.jctc.0c00450] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Accurate prediction of gas solubility in a liquid is crucial in many areas of chemistry, and a detailed understanding of the molecular mechanism of the gas solvation continues to be an active area of research. Here, we extend the idea of the constant chemical potential molecular dynamics (CμMD) approach to the calculation of the gas solubility in the liquid under constant gas chemical potential conditions. As a representative example, we utilize this method to calculate the isothermal solubility of carbon dioxide in water. Additionally, we provide microscopic insight into the mechanism of solvation that preferentially occurs in areas of the surface where the hydrogen network is broken.
Collapse
Affiliation(s)
- Narjes Ansari
- Department of Chemistry and Applied Biosciences, ETH Zurich, 8092 Zurich, Switzerland.,Facoltà di informatica, Istituto di Scienze Computazionali, Università della Svizzera Italiana, CH-6900 Lugano, Switzerland
| | - Tarak Karmakar
- Department of Chemistry and Applied Biosciences, ETH Zurich, 8092 Zurich, Switzerland.,Facoltà di informatica, Istituto di Scienze Computazionali, Università della Svizzera Italiana, CH-6900 Lugano, Switzerland
| | - Michele Parrinello
- Department of Chemistry and Applied Biosciences, ETH Zurich, 8092 Zurich, Switzerland.,Facoltà di informatica, Istituto di Scienze Computazionali, Università della Svizzera Italiana, CH-6900 Lugano, Switzerland.,Italian Institute of Technology, Via Morego 30, 16163 Genova, Italy
| |
Collapse
|
11
|
Naeiji P, Woo TK, Alavi S, Ohmura R. Molecular dynamics simulations of interfacial properties of the CO 2-water and CO 2-CH 4-water systems. J Chem Phys 2020; 153:044701. [PMID: 32752701 DOI: 10.1063/5.0008114] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022] Open
Abstract
Molecular dynamics simulations were performed to study the interfacial behavior of the pure carbon dioxide-water system and a binary 40:60 mol. % gas mixture of (carbon dioxide + methane)-water at the temperatures of 275.15 K and 298.15 K and pressures near 4 MPa for CO2 and up to 10 MPa for methane. The simulations are used to study the dynamic equilibrium of the gases at the water-gas interface, to determine the z-density profiles for the gases and water, and calculate the interfacial tension γ under the different temperature/pressure conditions close to those of the formation of clathrate hydrates of these gases. At the same hydrostatic gas phase pressure, the CO2-water interface has a lower interfacial tension than the CH4-water interface. A greater number of CO2 molecules, as much as three times more than methane at the same pressure, were adsorbed at the interfacial layer, which reflects the stronger electrostatic quadrupolar and van der Waals interactions between CO2 and water molecules at the interface. The water surfaces are covered by less than a monolayer of gas even when the pressure of the system goes near the saturation pressure of CO2. The surface adsorbed molecules are in dynamic equilibrium with the bulk gas and with exchange between the gas and interface regions occurring repeatedly within the timescale of the simulations. The effects of the changes in the CO2-water interfacial tension with external temperature and pressure conditions on the formation of the clathrate hydrates and other CO2 capture and sequestration processes are discussed.
Collapse
Affiliation(s)
- Parisa Naeiji
- Department of Chemistry and Biomolecular Sciences, University of Ottawa, Ottawa, Ontario K1A 0R6, Canada
| | - Tom K Woo
- Department of Chemistry and Biomolecular Sciences, University of Ottawa, Ottawa, Ontario K1A 0R6, Canada
| | - Saman Alavi
- Department of Chemistry and Biomolecular Sciences, University of Ottawa, Ottawa, Ontario K1A 0R6, Canada
| | - Ryo Ohmura
- Department of Mechanical Engineering, Keio University, 3-14-1 Hiyoshi, Kohoku-Ku, Yokohama 223-8522, Japan
| |
Collapse
|
12
|
Anwar J, Leitold C, Peters B. Solid–solid phase equilibria in the NaCl–KCl system. J Chem Phys 2020; 152:144109. [DOI: 10.1063/5.0003224] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023] Open
Affiliation(s)
- Jamshed Anwar
- Department of Chemistry, Lancaster University, Lancaster LA1 4YW, United Kingdom
| | - Christian Leitold
- Department of Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| | - Baron Peters
- Department of Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
- Department of Chemistry and Biochemistry, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| |
Collapse
|
13
|
Silvestri A, Ataman E, Budi A, Stipp SLS, Gale JD, Raiteri P. Wetting Properties of the CO 2-Water-Calcite System via Molecular Simulations: Shape and Size Effects. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2019; 35:16669-16678. [PMID: 31714788 DOI: 10.1021/acs.langmuir.9b02881] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Assessment of the risks and environmental impacts of carbon geosequestration requires knowledge about the wetting behavior of mineral surfaces in the presence of CO2 and the pore fluids. In this context, the interfacial tension (IFT) between CO2 and the aqueous fluid and the contact angle, θ, with the pore mineral surfaces are the two key parameters that control the capillary pressure in the pores of the candidate host rock. Knowledge of these two parameters and their dependence on the local conditions of pressure, temperature, and salinity is essential for the correct prediction of structural and residual trapping. We have performed classical molecular dynamics simulations to predict the CO2-water IFT and the CO2-water-calcite contact angle. The IFT results are consistent with previous simulations, where simple point charge water models have been shown to underestimate the water surface tension, thus affecting the simulated IFT values. When combined with the EPM2 CO2 model, the SPC/Fw water model indeed underestimates the IFT in the low-pressure region at all temperatures studied. On the other hand, at high pressure and low temperature, the IFT is overestimated by ∼5 mN/m. Literature data regarding the CO2/water/calcite contact angle on calcite are contradictory. Using our new set of force field parameters, we performed NVT simulations at 323 K and 20 MPa to calculate the contact angle of a water droplet on the calcite {10.4} surface in a CO2 atmosphere. We performed simulations for both spherical and cylindrical droplet configurations for different initial radii to study the size dependence of the water contact angle on calcite in the presence of CO2. Our results suggest that the contact angle of a cylindrical droplet, is independent of droplet size, for droplets with a radius of 50 Å or more. On the contrary, spherical droplets make a contact angle that is strongly influenced by their size. At the largest size explored in this study, both spherical and cylindrical droplets converge to the same contact angle, 38°, indicating that calcite is strongly wetted by water.
Collapse
Affiliation(s)
- A Silvestri
- Curtin Institute for Computation, The Institute for Geoscience Research (TIGeR), School of Molecular and Life Sciences , Curtin University , PO Box U1987, Perth , WA 6845 , Australia
| | - E Ataman
- Nano-Science Center, Department of Chemistry , University of Copenhagen , Universitetsparken 5 , København Ø DK-2100 , Denmark
| | - A Budi
- Institute for Frontier Materials , Deakin University , Geelong , VIC 3216 , Australia
| | - S L S Stipp
- Department of Physics , Technical University of Denmark , Fysikvej , DK-2800 Kongens Lyngby , Denmark
| | - J D Gale
- Curtin Institute for Computation, The Institute for Geoscience Research (TIGeR), School of Molecular and Life Sciences , Curtin University , PO Box U1987, Perth , WA 6845 , Australia
| | - P Raiteri
- Curtin Institute for Computation, The Institute for Geoscience Research (TIGeR), School of Molecular and Life Sciences , Curtin University , PO Box U1987, Perth , WA 6845 , Australia
| |
Collapse
|
14
|
Lopez-Lazaro C, Bachaud P, Moretti I, Ferrando N. Predicting the phase behavior of hydrogen in NaCl brines by molecular simulation for geological applications. ACTA ACUST UNITED AC 2019. [DOI: 10.1051/bsgf/2019008] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022]
Abstract
Hydrogen is targeted to have a significant influence on the energy mix in the upcoming years. Its underground injection is an efficient solution for large-scale and long-term storage. Furthermore, natural hydrogen emissions have been proven in several locations of the world, and the potential underground accumulations constitute exciting carbon-free energy sources. In this context, comprehensive models are necessary to better constrain hydrogen behavior in geological formations. In particular, solubility in brines is a key-parameter, as it directly impacts hydrogen reactivity and migration in porous media. In this work, Monte Carlo simulations have been carried out to generate new simulated data of hydrogen solubility in aqueous NaCl solutions in temperature and salinity ranges of interest for geological applications, and for which no experimental data are currently available. For these simulations, molecular models have been selected for hydrogen, water and Na+ and Cl− to reproduce phase properties of pure components and brine densities. To model solvent-solutes and solutes-solutes interactions, it was shown that the Lorentz-Berthelot mixing rules with a constant interaction binary parameter are the most appropriate to reproduce the experimental hydrogen Henry constants in salted water. With this force field, simulation results match measured solubilities with an average deviation of 6%. Additionally, simulation reproduced the expected behaviors of the H2O + H2 + NaCl system, such as the salting-out effect, a minimum hydrogen solubility close to 57 °C, and a decrease of the Henry constant with increasing temperature. The force field was then used in extrapolation to determine hydrogen Henry constants for temperatures up to 300 °C and salinities up to 2 mol/kgH2O. Using the experimental measures and these new simulated data generated by molecular simulation, a binary interaction parameter of the Soreide and Whiston equation of state has been fitted. The obtained model allows fast and reliable phase equilibrium calculations, and it was applied to illustrative cases relevant for hydrogen geological storage or H2 natural emissions.
Collapse
|
15
|
Yang Y, Che Ruslan MFA, Narayanan Nair AK, Sun S. Effect of Ion Valency on the Properties of the Carbon Dioxide-Methane-Brine System. J Phys Chem B 2019; 123:2719-2727. [PMID: 30843699 DOI: 10.1021/acs.jpcb.8b12033] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Abstract
Molecular dynamics simulations and theoretical analysis were carried out to study the bulk and interfacial properties of carbon dioxide-methane-water and carbon dioxide-methane-brine systems under geological conditions. The density gradient theory with the bulk phase properties estimated using the cubic-plus-association (CPA) equation of state (EoS) can well describe the increase in the interfacial tension (IFT) of the CO2-water system in the presence of methane. The theoretical estimates of species mole fractions in the carbon dioxide-methane-water system are in good quantitative agreement with the experimental results. Furthermore, simulations of carbon dioxide-methane-brine system show that the IFT of the CaCl2 case is generally higher than that of the NaCl case. This is probably due to the stronger hydration of Ca2+ ions and their stronger repulsion from the interface compared to Na+. While the overall shape of the ionic profiles is not much affected by the ion type, the water profiles show a local enrichment at the interface in the system with CaCl2. In contrast to the case of NaCl, the slopes of the plots of IFT vs CaCl2 concentration are dependent on temperature. Species mole fractions in the carbon dioxide-methane-brine system predicted by combining the CPA EoS with the Debye-Hückel electrostatic term are in good agreement with simulation results.
Collapse
Affiliation(s)
- Yafan Yang
- Physical Science and Engineering Division (PSE), Computational Transport Phenomena Laboratory , King Abdullah University of Science and Technology (KAUST) , Thuwal 23955-6900 , Saudi Arabia
| | - Mohd Fuad Anwari Che Ruslan
- Physical Science and Engineering Division (PSE), Computational Transport Phenomena Laboratory , King Abdullah University of Science and Technology (KAUST) , Thuwal 23955-6900 , Saudi Arabia
| | - Arun Kumar Narayanan Nair
- Physical Science and Engineering Division (PSE), Computational Transport Phenomena Laboratory , King Abdullah University of Science and Technology (KAUST) , Thuwal 23955-6900 , Saudi Arabia
| | - Shuyu Sun
- Physical Science and Engineering Division (PSE), Computational Transport Phenomena Laboratory , King Abdullah University of Science and Technology (KAUST) , Thuwal 23955-6900 , Saudi Arabia
| |
Collapse
|
16
|
Yang Y, Narayanan Nair AK, Sun S. Molecular Dynamics Simulation Study of Carbon Dioxide, Methane, and Their Mixture in the Presence of Brine. J Phys Chem B 2017; 121:9688-9698. [DOI: 10.1021/acs.jpcb.7b08118] [Citation(s) in RCA: 38] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Affiliation(s)
- Yafan Yang
- Physical Science and Engineering Division
(PSE), Computational Transport Phenomena Laboratory, King Abdullah University of Science and Technology (KAUST), Thuwal, Jeddah 23955-6900, Saudi Arabia
| | - Arun Kumar Narayanan Nair
- Physical Science and Engineering Division
(PSE), Computational Transport Phenomena Laboratory, King Abdullah University of Science and Technology (KAUST), Thuwal, Jeddah 23955-6900, Saudi Arabia
| | - Shuyu Sun
- Physical Science and Engineering Division
(PSE), Computational Transport Phenomena Laboratory, King Abdullah University of Science and Technology (KAUST), Thuwal, Jeddah 23955-6900, Saudi Arabia
| |
Collapse
|
17
|
Jiang H, Economou IG, Panagiotopoulos AZ. Molecular Modeling of Thermodynamic and Transport Properties for CO 2 and Aqueous Brines. Acc Chem Res 2017; 50:751-758. [PMID: 28234455 DOI: 10.1021/acs.accounts.6b00632] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Molecular simulation techniques using classical force-fields occupy the space between ab initio quantum mechanical methods and phenomenological correlations. In particular, Monte Carlo and molecular dynamics algorithms can be used to provide quantitative predictions of thermodynamic and transport properties of fluids relevant for geologic carbon sequestration at conditions for which experimental data are uncertain or not available. These methods can cover time and length scales far exceeding those of quantum chemical methods, while maintaining transferability and predictive power lacking from phenomenological correlations. The accuracy of predictions depends sensitively on the quality of the molecular models used. Many existing fixed-point-charge models for water and aqueous mixtures fail to represent accurately these fluid properties, especially when descriptions covering broad ranges of thermodynamic conditions are needed. Recent work on development of accurate models for water, CO2, and dissolved salts, as well as their mixtures, is summarized in this Account. Polarizable models that can respond to the different dielectric environments in aqueous versus nonaqueous phases are necessary for predictions of properties over extended ranges of temperatures and pressures. Phase compositions and densities, activity coefficients of the dissolved salts, interfacial tensions, viscosities and diffusivities can be obtained in near-quantitative agreement to available experimental data, using relatively modest computational resources. In some cases, for example, for the composition of the CO2-rich phase in coexistence with an aqueous phase, recent results from molecular simulations have helped discriminate among conflicting experimental data sets. The sensitivity of properties on the quality of the intermolecular interaction model varies significantly. Properties such as the phase compositions or electrolyte activity coefficients are much more sensitive than phase densities, viscosities, or component diffusivities. Strong confinement effects on physical properties in nanoscale media can also be directly obtained from molecular simulations. Future work on molecular modeling for CO2 and aqueous brines is likely to be focused on more systematic generation of interaction models by utilizing quantum chemical as well as direct experimental measurements. New ion models need to be developed for use with the current generation of polarizable water models, including ion-ion interactions that will allow for accurate description of dense, mixed brines. Methods will need to be devised that go beyond the use of effective potentials for incorporation of quantum effects known to be important for water, and reactive force fields developed that can handle bond creation and breaking in systems with carbonate and silicate minerals. Another area of potential future work is the integration of molecular simulation methods in multiscale models for the chemical reactions leading to mineral dissolution and flow within the porous media in underground formations.
Collapse
Affiliation(s)
- Hao Jiang
- Department
of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, United States
| | - Ioannis G. Economou
- Chemical Engineering Program, Texas A&M University at Qatar, P.O. Box 23874, Doha, Qatar
| | | |
Collapse
|
18
|
Silvestri A, Stipp SLS, Andersson MP. Predicting CO2–H2O Interfacial Tension Using COSMO-RS. J Chem Theory Comput 2017; 13:804-810. [DOI: 10.1021/acs.jctc.6b00818] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- A. Silvestri
- Nano-Science Center, Department
of Chemistry, University of Copenhagen, Universitetsparken 5, DK-2100 København Ø, Denmark
| | - S. L. S. Stipp
- Nano-Science Center, Department
of Chemistry, University of Copenhagen, Universitetsparken 5, DK-2100 København Ø, Denmark
| | - M. P. Andersson
- Nano-Science Center, Department
of Chemistry, University of Copenhagen, Universitetsparken 5, DK-2100 København Ø, Denmark
| |
Collapse
|
19
|
Zhao L, Ji J, Tao L, Lin S. Ionic Effects on Supercritical CO2-Brine Interfacial Tensions: Molecular Dynamics Simulations and a Universal Correlation with Ionic Strength, Temperature, and Pressure. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2016; 32:9188-9196. [PMID: 27564433 DOI: 10.1021/acs.langmuir.6b02485] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
For geological CO2 storage in deep saline aquifers, the interfacial tension (IFT) between supercritical CO2 and brine is critical for the storage security and design of the storage capacitance. However, currently, no predictive model exists to determine the IFT of supercritical CO2 against complex electrolyte solutions involving various mixed salt species at different concentrations and compositions. In this paper, we use molecular dynamics (MD) simulations to investigate the effect of salt ions on the incremental IFT at the supercritical CO2-brine interface with respect to that at the reference supercritical CO2-water interface. Supercritical CO2-NaCl solution, CO2-CaCl2 solution and CO2-(NaCl+CaCl2) mixed solution systems are simulated at 343 K and 20 MPa under different salinities and salt compositions. We find that the valence of the cations is the primary contributor to the variation in IFT, while the Lennard-Jones potentials for the cations pose a smaller impact on the IFT. Interestingly, the incremental IFT exhibits a general linear correlation with the ionic strength in the above three electrolyte systems, and the slopes are almost identical and independent of the solution types. Based on this finding, a universal predictive formula for IFTs of CO2-complex electrolyte solution systems is established, as a function of ionic strength, temperature, and pressure. The predicted IFTs using the established formula agree perfectly (with a high statistical confidence level of ∼96%) with a wide range of experimental data for CO2 interfacing with different electrolyte solutions, such as those involving MgCl2 and Na2SO4. This work provides an efficient and accurate route to directly predict IFTs in supercritical CO2-complex electrolyte solution systems for practical engineering applications, such as geological CO2 sequestration in deep saline aquifers and other interfacial systems involving complex electrolyte solutions.
Collapse
Affiliation(s)
- Lingling Zhao
- Key Laboratory of Energy Thermal Conversion and Control of Ministry of Education, School of Energy & Environment, Southeast University , Nanjing, Jiangsu 210096, China
| | - Jiayuan Ji
- Key Laboratory of Energy Thermal Conversion and Control of Ministry of Education, School of Energy & Environment, Southeast University , Nanjing, Jiangsu 210096, China
| | - Lu Tao
- Key Laboratory of Energy Thermal Conversion and Control of Ministry of Education, School of Energy & Environment, Southeast University , Nanjing, Jiangsu 210096, China
| | - Shangchao Lin
- Department of Mechanical Engineering, Materials Science & Engineering Program, FAMU-FSU College of Engineering, Florida State University , Tallahassee, Florida 32310, United States
| |
Collapse
|
20
|
Bourg IC, Beckingham LE, DePaolo DJ. The Nanoscale Basis of CO2 Trapping for Geologic Storage. ENVIRONMENTAL SCIENCE & TECHNOLOGY 2015; 49:10265-10284. [PMID: 26266820 DOI: 10.1021/acs.est.5b03003] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
Carbon capture and storage (CCS) is likely to be a critical technology to achieve large reductions in global carbon emissions over the next century. Research on the subsurface storage of CO2 is aimed at reducing uncertainties in the efficacy of CO2 storage in sedimentary rock formations. Three key parameters that have a nanoscale basis and that contribute uncertainty to predictions of CO2 trapping are the vertical permeability kv of seals, the residual CO2 saturation Sg,r in reservoir rocks, and the reactive surface area ar of silicate minerals. This review summarizes recent progress and identifies outstanding research needs in these areas. Available data suggest that the permeability of shale and mudstone seals is heavily dependent on clay fraction and can be extremely low even in the presence of fractures. Investigations of residual CO2 trapping indicate that CO2-induced alteration in the wettability of mineral surfaces may significantly influence Sg,r. Ultimately, the rate and extent of CO2 conversion to mineral phases are uncertain due to a poor understanding of the kinetics of slow reactions between minerals and fluids. Rapidly improving characterization techniques using X-rays and neutrons, and computing capability for simulating chemical interactions, provide promise for important advances.
Collapse
Affiliation(s)
- Ian C Bourg
- Department of Civil and Environmental Engineering and Princeton Environmental Institute, Princeton University , E-208 E-Quad, Princeton, New Jersey 08544, United States
- Earth Sciences Division, Lawrence Berkeley National Laboratory , 1 Cyclotron Road, Berkeley, California 94720, United States
| | - Lauren E Beckingham
- Earth Sciences Division, Lawrence Berkeley National Laboratory , 1 Cyclotron Road, Berkeley, California 94720, United States
| | - Donald J DePaolo
- Earth Sciences Division, Lawrence Berkeley National Laboratory , 1 Cyclotron Road, Berkeley, California 94720, United States
| |
Collapse
|
21
|
Jiang H, Mester Z, Moultos OA, Economou IG, Panagiotopoulos AZ. Thermodynamic and Transport Properties of H2O + NaCl from Polarizable Force Fields. J Chem Theory Comput 2015; 11:3802-10. [DOI: 10.1021/acs.jctc.5b00421] [Citation(s) in RCA: 53] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Hao Jiang
- Department
of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, United States
| | - Zoltan Mester
- Department
of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, United States
| | - Othonas A. Moultos
- Chemical Engineering Program, Texas A&M University at Qatar, P.O. Box 23874, Doha, Qatar
| | - Ioannis G. Economou
- Chemical Engineering Program, Texas A&M University at Qatar, P.O. Box 23874, Doha, Qatar
| | | |
Collapse
|
22
|
Javanbakht G, Sedghi M, Welch W, Goual L. Molecular Dynamics Simulations of CO2/Water/Quartz Interfacial Properties: Impact of CO2 Dissolution in Water. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2015; 31:5812-5819. [PMID: 25965772 DOI: 10.1021/acs.langmuir.5b00445] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
The safe trapping of carbon dioxide (CO2) in deep saline aquifers is one of the major concerns of CO2 sequestration. The amount of capillary trapping is dominated by the capillary pressure of water and CO2 inside the reservoir, which in turn is controlled by the interfacial tension (IFT) and the contact angle (CA) of CO2/water/rock systems. The measurement of IFT and CA could be very challenging at reservoir conditions, especially in the presence of toxic cocontaminants. Thus, the ability to accurately predict these interfacial properties at reservoir conditions is very advantageous. Although the majority of existing molecular dynamics (MD) studies of CO2/water/mineral systems were able to capture the trends in IFT and CA variations with pressure and temperature, their predictions often deviated from experimental data, possibly due to erroneous models and/or overlooked chemical reactions. The objective of this study was to improve the MD predictions of IFT and CA of CO2/water/quartz systems at various pressure and temperature conditions by (i) considering the chemical reactions between CO2 and water and (ii) using a new molecular model for α-quartz surface. The results showed that the presence of carbonic acid at the CO2/water interface improved the predictions of IFT, especially at low temperature and high pressure where more CO2 dissolution occurs. On the other hand, the effect on CA was minor. The slight decrease in CA observed across the pressure range investigated could be attributed to an increase in the total number of H-bonds between fluid molecules and quartz surface.
Collapse
Affiliation(s)
- Gina Javanbakht
- Department of Chemical and Petroleum Engineering, University of Wyoming, 1000 E. University Avenue, Laramie, Wyoming 82071, United States
| | - Mohammad Sedghi
- Department of Chemical and Petroleum Engineering, University of Wyoming, 1000 E. University Avenue, Laramie, Wyoming 82071, United States
| | - William Welch
- Department of Chemical and Petroleum Engineering, University of Wyoming, 1000 E. University Avenue, Laramie, Wyoming 82071, United States
| | - Lamia Goual
- Department of Chemical and Petroleum Engineering, University of Wyoming, 1000 E. University Avenue, Laramie, Wyoming 82071, United States
| |
Collapse
|
23
|
Moultos OA, Orozco GA, Tsimpanogiannis IN, Panagiotopoulos AZ, Economou IG. Atomistic molecular dynamics simulations of H2O diffusivity in liquid and supercritical CO2. Mol Phys 2015. [DOI: 10.1080/00268976.2015.1023224] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Affiliation(s)
- Othonas A. Moultos
- Chemical Engineering Program , Texas A&M University at Qatar , Doha, Qatar
| | - Gustavo A. Orozco
- Department of Chemical and Biological Engineering, Princeton University , Princeton, United States
| | - Ioannis N. Tsimpanogiannis
- Chemical Engineering Program , Texas A&M University at Qatar , Doha, Qatar
- Environmental Research Laboratory, Institute of Nuclear and Radiology Sciences and Technology, Energy and Safety, National Center for Scientific Research “Demokritos” , Aghia Paraskevi, Greece
| | | | | |
Collapse
|
24
|
Orozco GA, Moultos OA, Jiang H, Economou IG, Panagiotopoulos AZ. Molecular simulation of thermodynamic and transport properties for the H2O+NaCl system. J Chem Phys 2014; 141:234507. [DOI: 10.1063/1.4903928] [Citation(s) in RCA: 51] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022] Open
Affiliation(s)
- Gustavo A. Orozco
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA
| | - Othonas A. Moultos
- Chemical Engineering Program, Texas A&M University at Qatar, P.O. Box 23874, Doha, Qatar
| | - Hao Jiang
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA
| | - Ioannis G. Economou
- Chemical Engineering Program, Texas A&M University at Qatar, P.O. Box 23874, Doha, Qatar
| | | |
Collapse
|
25
|
Moultos OA, Tsimpanogiannis IN, Panagiotopoulos AZ, Economou IG. Atomistic molecular dynamics simulations of CO₂ diffusivity in H₂O for a wide range of temperatures and pressures. J Phys Chem B 2014; 118:5532-41. [PMID: 24749622 DOI: 10.1021/jp502380r] [Citation(s) in RCA: 39] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
Molecular dynamics simulations were employed for the calculation of diffusion coefficients of CO2 in H2O. Various combinations of existing force fields for H2O (SPC, SPC/E, and TIP4P/2005) and CO2 (EPM2 and TraPPE) were tested over a wide range of temperatures (283.15 K < T < 623.15 K) and pressures (0.1 MPa < P < 100.0 MPa). All force-field combinations qualitatively reproduce the trends of the experimental data; however, two specific combinations were found to be more accurate. In particular, at atmospheric pressure, the TIP4P/2005-EPM2 combination was found to perform better for temperatures lower than 323.15 K, while the SPC/E-TraPPE combination was found to perform better at higher temperatures. The pressure dependence of the diffusion coefficient of CO2 in H2O at constant temperature is shown to be negligible at temperatures lower than 473.15 K, in good agreement with experiments. As temperature increases, the pressure effect becomes substantial. The phenomenon is driven primarily by the higher compressibility of liquid H2O at near-critical conditions. Finally, a simple power-law-type phenomenological equation is proposed to correlate the simulation values; the proposed correlation should be useful for engineering calculations.
Collapse
Affiliation(s)
- Othonas A Moultos
- Chemical Engineering Program, Texas A&M University at Qatar , P.O. Box 23874, Doha, Qatar
| | | | | | | |
Collapse
|
26
|
Nikoubashman A, Register RA, Panagiotopoulos AZ. Sequential Domain Realignment Driven by Conformational Asymmetry in Block Copolymer Thin Films. Macromolecules 2014. [DOI: 10.1021/ma402526q] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Affiliation(s)
- Arash Nikoubashman
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, United Sates
| | - Richard A. Register
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, United Sates
| | | |
Collapse
|