1
|
Yang Y, Wan J, Li J, Zhao G, Shang X. Molecular modeling of interfacial properties of the hydrogen + water + decane mixture in three-phase equilibrium. Phys Chem Chem Phys 2023; 25:29641-29655. [PMID: 37881893 DOI: 10.1039/d3cp04406h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2023]
Abstract
The understanding of interfacial phenomena between H2 and geofluids is of great importance for underground H2 storage, but requires further study. We report the first investigation on the three-phase fluid mixture containing H2, H2O, and n-C10H22. Molecular dynamics simulation and PC-SAFT density gradient theory are employed to estimate the interfacial properties under various conditions (temperature ranges from 298 to 373 K and pressure is up to around 100 MPa). Our results demonstrate that interfacial tensions (IFTs) of the H2-H2O interface in the H2 + H2O + C10H22 three-phase mixture are smaller than IFTs in the H2 + H2O two-phase mixture. This decrement of IFT can be attributed to C10H22 adsorption in the interface. Importantly, H2 accumulates in the H2O-C10H22 interface in the three-phase systems, which leads to weaker increments of IFT with increasing pressure compared to IFTs in the water + C10H22 two-phase mixture. In addition, the IFTs of the H2-C10H22 interface are hardly influenced by H2O due to the limited amount of H2O dissolved in nonaqueous phases. Nevertheless, positive surface excesses of H2O are seen in the H2-C10H22 interfacial region. Furthermore, the values of the spreading coefficient are mostly negative revealing the presence of the three-phase contact for the H2 + H2O + C10H22 mixture under studied conditions.
Collapse
Affiliation(s)
- Yafan Yang
- State Key Laboratory for Geomechanics and Deep Underground Engineering, China University of Mining and Technology, Xuzhou, China.
| | - Jingyu Wan
- State Key Laboratory for Geomechanics and Deep Underground Engineering, China University of Mining and Technology, Xuzhou, China.
| | - Jingfa Li
- School of Mechanical Engineering and Hydrogen Energy Research Center, Beijing Institute of Petrochemical Technology, Beijing, China
| | - Guangsi Zhao
- State Key Laboratory for Geomechanics and Deep Underground Engineering, China University of Mining and Technology, Xuzhou, China.
| | - Xiangyu Shang
- State Key Laboratory for Geomechanics and Deep Underground Engineering, China University of Mining and Technology, Xuzhou, China.
| |
Collapse
|
2
|
Bråten V, Schaefer D, Stephan S, Hasse H. Molecular dynamics simulation study on the mass transfer across vapor-liquid interfaces in azeotropic mixtures. J Chem Phys 2023; 159:084503. [PMID: 37622596 DOI: 10.1063/5.0165421] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2023] [Accepted: 07/24/2023] [Indexed: 08/26/2023] Open
Abstract
Mass transfer through fluid interfaces is an important phenomenon in industrial applications as well as in naturally occurring processes. In this work, we investigate the mass transfer across vapor-liquid interfaces in binary mixtures using molecular dynamics simulations. We investigate the influence of interfacial properties on mass transfer by studying three binary azeotropic mixtures known to have different interfacial behaviors. Emphasis is placed on the effect of the intermolecular interactions by choosing mixtures with the same pure components but different cross-interactions such that different azeotropic behaviors are obtained. The molar flux is created by utilizing a non-stationary molecular dynamics simulation approach, where particles of one component are inserted into the vapor phase over a short period of time before the system's response to this insertion is monitored. From a direct comparison of the density profiles and the flux profiles in close proximity to the interface, we analyze the particles' tendency to accumulate in the interfacial region throughout the different stages of the simulation. We find that for mixtures with strong attractive cross-interactions, the inserted particles are efficiently transported into the liquid phase. For systems with weak attractive cross-interactions, the inserted particles show a tendency to accumulate in the interfacial region, and the flux through the system is lower. The results from this work indicate that the accumulation of particles at the interface can act as a hindrance to mass transfer, which has practical relevance in technical processes.
Collapse
Affiliation(s)
- Vilde Bråten
- Laboratory of Engineering Thermodynamics (LTD), RPTU Kaiserslautern, Kaiserslautern, Germany
| | - Dominik Schaefer
- Laboratory of Engineering Thermodynamics (LTD), RPTU Kaiserslautern, Kaiserslautern, Germany
| | - Simon Stephan
- Laboratory of Engineering Thermodynamics (LTD), RPTU Kaiserslautern, Kaiserslautern, Germany
| | - Hans Hasse
- Laboratory of Engineering Thermodynamics (LTD), RPTU Kaiserslautern, Kaiserslautern, Germany
| |
Collapse
|
3
|
Schaefer D, Stephan S, Langenbach K, Horsch MT, Hasse H. Mass Transfer through Vapor-Liquid Interfaces Studied by Non-Stationary Molecular Dynamics Simulations. J Phys Chem B 2023; 127:2521-2533. [PMID: 36896991 DOI: 10.1021/acs.jpcb.2c08752] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/11/2023]
Abstract
Molecular dynamics (MD) simulations are highly attractive for studying the influence of interfacial effects, such as the enrichment of components, on the mass transfer through the interface. In a recent work, we have presented a steady-state MD simulation method for investigating this phenomenon and tested it using model mixtures with and without interfacial enrichment. The present study extends this work by introducing a non-stationary MD simulation method. A rectangular simulation box that contains a mixture of two components 1 + 2 with a vapor phase in the middle and two liquid phases on both sides is used. Starting from a vapor-liquid equilibrium state, a non-stationary molar flux of component 2 is induced by inserting particles of component 2 into the center of the vapor phase in a pulse-like manner. During the isothermal relaxation process, particles of component 2 pass through the vapor phase, cross the vapor-liquid interface, and enter the liquid phase. The system thereby relaxes into a new vapor-liquid equilibrium state. During the relaxation process, spatially resolved responses for the component densities, fluxes, and pressure are sampled. To reduce the noise and provide measures for the uncertainty of the observables, a set of replicas of simulations is carried out. The new simulation method was applied to study mass transfer in two binary Lennard-Jones mixtures: one that exhibits a strong enrichment of the low-boiling component 2 at the vapor-liquid interface and one that shows no enrichment. Even though both mixtures have similar transport coefficients in the bulk phases, the results for mass transfer differ significantly, indicating that the interfacial enrichment influences the mass transfer.
Collapse
Affiliation(s)
- Dominik Schaefer
- Laboratory of Engineering Thermodynamics (LTD), TU Kaiserslautern, 67663 Kaiserslautern, Germany
| | - Simon Stephan
- Laboratory of Engineering Thermodynamics (LTD), TU Kaiserslautern, 67663 Kaiserslautern, Germany
| | - Kai Langenbach
- Institute of Chemical Engineering, University of Innsbruck, 6020 Innsbruck, Austria
| | - Martin T Horsch
- Norwegian University of Life Sciences, Faculty of Science and Technology, Department of Data Science, Drøbakveien 31, 1430 Ås, Norway
| | - Hans Hasse
- Laboratory of Engineering Thermodynamics (LTD), TU Kaiserslautern, 67663 Kaiserslautern, Germany
| |
Collapse
|
4
|
Aasen A, Wilhelmsen Ø, Hammer M, Reguera D. Free energy of critical droplets-from the binodal to the spinodal. J Chem Phys 2023; 158:114108. [PMID: 36948791 DOI: 10.1063/5.0142533] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/03/2023] Open
Abstract
Arguably, the main challenge of nucleation theory is to accurately evaluate the work of formation of a critical embryo in the new phase, which governs the nucleation rate. In Classical Nucleation Theory (CNT), this work of formation is estimated using the capillarity approximation, which relies on the value of the planar surface tension. This approximation has been blamed for the large discrepancies between predictions from CNT and experiments. In this work, we present a study of the free energy of formation of critical clusters of the Lennard-Jones fluid truncated and shifted at 2.5σ using Monte Carlo simulations, density gradient theory, and density functional theory. We find that density gradient theory and density functional theory accurately reproduce molecular simulation results for critical droplet sizes and their free energies. The capillarity approximation grossly overestimates the free energy of small droplets. The incorporation of curvature corrections up to the second order with the Helfrich expansion greatly remedies this and performs very well for most of the experimentally accessible regions. However, it is imprecise for the smallest droplets and largest metastabilities since it does not account for a vanishing nucleation barrier at the spinodal. To remedy this, we propose a scaling function that uses all relevant ingredients without adding fitting parameters. The scaling function reproduces accurately the free energy of the formation of critical droplets for the entire metastability range and all temperatures examined and deviates from density gradient theory by less than one kBT.
Collapse
Affiliation(s)
- Ailo Aasen
- SINTEF Energy Research, NO-7465 Trondheim, Norway
| | | | | | - David Reguera
- Departament de Física de la Matèria Condensada, Universitat de Barcelona, Martí i Franquès 1, 08028 Barcelona, Spain
| |
Collapse
|
5
|
Urschel M, Stephan S. Determining Brown's Characteristic Curves Using Molecular Simulation. J Chem Theory Comput 2023; 19:1537-1552. [PMID: 36802608 DOI: 10.1021/acs.jctc.2c01102] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/23/2023]
Abstract
Brown's characteristic curves define lines on the thermodynamic surface where special thermodynamic conditions hold. These curves are an important tool for the development of thermodynamic models of fluids. Yet, practically no experimental data for Brown's characteristic curves is available. In this work, a rigorous and generalized method for determining Brown's characteristic curves based on molecular simulation was developed. As multiple thermodynamic equivalent definitions apply for the characteristic curves, different simulation routes were compared. Based on this systematic approach, the most favorable route for determining each characteristic curve was identified. The computational procedure developed in this work combines molecular simulation, molecular-based equation of state, and the evaluation of the second virial coefficient. The new method was tested on a simple model system (the classical Lennard-Jones fluid) and different types of real substances (toluene, methane, ethane, propane, and ethanol). It is thereby shown that the method is robust and yields accurate results. Moreover, a computer code implementation of the method is presented.
Collapse
Affiliation(s)
- Maximilian Urschel
- Laboratory of Engineering Thermodynamics (LTD), TU Kaiserslautern, Kaiserslautern 67663, Germany
| | - Simon Stephan
- Laboratory of Engineering Thermodynamics (LTD), TU Kaiserslautern, Kaiserslautern 67663, Germany
| |
Collapse
|
6
|
Rehner P, Bauer G, Gross J. FeO s: An Open-Source Framework for Equations of State and Classical Density Functional Theory. Ind Eng Chem Res 2023. [DOI: 10.1021/acs.iecr.2c04561] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/11/2023]
Affiliation(s)
- Philipp Rehner
- Energy and Process Systems Engineering, Department of Mechanical and Process Engineering, ETH Zurich, Tannenstrasse 3, Zurich 8092, Switzerland
| | - Gernot Bauer
- Institute of Thermodynamics and Thermal Process Engineering, University of Stuttgart, Pfaffenwaldring 9, Stuttgart 70569, Germany
| | - Joachim Gross
- Institute of Thermodynamics and Thermal Process Engineering, University of Stuttgart, Pfaffenwaldring 9, Stuttgart 70569, Germany
| |
Collapse
|
7
|
Fertig D, Stephan S. Influence of dispersive long-range interactions on transport and excess properties of simple mixtures. Mol Phys 2023. [DOI: 10.1080/00268976.2022.2162993] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Affiliation(s)
- David Fertig
- Laboratory of Engineering Thermodynamics (LTD), TU Kaiserslautern, Kaiserslautern, Germany
| | - Simon Stephan
- Laboratory of Engineering Thermodynamics (LTD), TU Kaiserslautern, Kaiserslautern, Germany
| |
Collapse
|
8
|
Staubach J, Stephan S. Interfacial properties of binary azeotropic mixtures of simple fluids: Molecular dynamics simulation and density gradient theory. J Chem Phys 2022; 157:124702. [PMID: 36182407 DOI: 10.1063/5.0100728] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Interfacial properties of binary azeotropic mixtures of Lennard-Jones truncated and shifted fluids were studied by molecular dynamics (MD) simulation and density gradient theory (DGT) in combination with an equation of state. Three binary mixtures were investigated, which differ in the energetic cross interaction parameter that yields different types of azeotropic behavior. This study covers a wide temperature and composition range. Mixture A exhibits a heteroazeotrope at low temperatures, which changes to a low-boiling azeotrope at high temperatures, mixture B exhibits a low-boiling azeotrope, and mixture C exhibits a high-boiling azeotrope. The phase behavior and fluid interfacial properties as well as their relation were studied. Vapor-liquid, liquid-liquid, and vapor-liquid-liquid equilibria and interfaces were considered. Density profiles, the surface tension, the interfacial thickness, as well as the relative adsorption and enrichment of the components at the interface were studied. The results obtained from the two independent methods (MD and DGT) are overall in good agreement. The results provide insights into the relation of the phase behavior, particularly the azeotropic behavior, of simple fluid mixtures and the corresponding interfacial properties. Strong enrichment was found for the mixture with a heteroazeotrope in the vicinity of the three-phase equilibrium, which is related to a wetting transition.
Collapse
Affiliation(s)
- Jens Staubach
- Laboratory of Engineering Thermodynamics (LTD), TU Kaiserslautern, 67663 Kaiserslautern, Germany
| | - Simon Stephan
- Laboratory of Engineering Thermodynamics (LTD), TU Kaiserslautern, 67663 Kaiserslautern, Germany
| |
Collapse
|
9
|
Transport Properties of Binary Lennard-Jones Mixtures: Insights from Entropy Scaling and Conformal Solution Theory. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.120401] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
10
|
Schmitt S, Vo T, Lautenschlaeger MP, Stephan S, Hasse H. Molecular dynamics simulation study of heat transfer across solid–fluid interfaces in a simple model system. Mol Phys 2022. [DOI: 10.1080/00268976.2022.2057364] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Affiliation(s)
- Sebastian Schmitt
- Laboratory of Engineering Thermodynamics (LTD), TU Kaiserslautern, Kaiserslautern, Germany
| | - Truong Vo
- Laboratory of Engineering Thermodynamics (LTD), TU Kaiserslautern, Kaiserslautern, Germany
| | | | - Simon Stephan
- Laboratory of Engineering Thermodynamics (LTD), TU Kaiserslautern, Kaiserslautern, Germany
| | - Hans Hasse
- Laboratory of Engineering Thermodynamics (LTD), TU Kaiserslautern, Kaiserslautern, Germany
| |
Collapse
|
11
|
van Westen T, Hammer M, Hafskjold B, Aasen A, Gross J, Wilhelmsen Ø. Perturbation theories for fluids with short-ranged attractive forces: A case study of the Lennard-Jones spline fluid. J Chem Phys 2022; 156:104504. [DOI: 10.1063/5.0082690] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
It is generally not straightforward to apply molecular-thermodynamic theories to fluids with short-ranged attractive forces between their constituent molecules (or particles). This especially applies to perturbation theories, which, for short-ranged attractive fluids, typically must be extended to high order or may not converge at all. Here, we show that a recent first-order perturbation theory, the uv-theory, holds promise for describing such fluids. As a case study, we apply the uv-theory to a fluid with pair interactions defined by the Lennard-Jones spline potential, which is a short-ranged version of the LJ potential that is known to provide a challenge for equation-of-state development. The results of the uv-theory are compared to those of third-order Barker–Henderson and fourth-order Weeks–Chandler–Andersen perturbation theories, which are implemented using Monte Carlo simulation results for the respective perturbation terms. Theoretical predictions are compared to an extensive dataset of molecular simulation results from this (and previous) work, including vapor–liquid equilibria, first- and second-order derivative properties, the critical region, and metastable states. The uv-theory proves superior for all properties examined. An especially accurate description of metastable vapor and liquid states is obtained, which might prove valuable for future applications of the equation-of-state model to inhomogeneous phases or nucleation processes. Although the uv-theory is analytic, it accurately describes molecular simulation results for both the critical point and the binodal up to at least 99% of the critical temperature. This suggests that the difficulties typically encountered in describing the vapor–liquid critical region are only to a small extent caused by non-analyticity.
Collapse
Affiliation(s)
- Thijs van Westen
- Institute of Thermodynamics and Thermal Process Engineering, University of Stuttgart, Pfaffenwaldring 9, D-70569 Stuttgart, Germany
| | - Morten Hammer
- Department of Gas Technology, SINTEF Energy Research, NO-7465 Trondheim, Norway
| | - Bjørn Hafskjold
- Porelab, Department of Chemistry, Norwegian University of Science and Technology, NO-7491 Trondheim, Norway
| | - Ailo Aasen
- Department of Gas Technology, SINTEF Energy Research, NO-7465 Trondheim, Norway
| | - Joachim Gross
- Institute of Thermodynamics and Thermal Process Engineering, University of Stuttgart, Pfaffenwaldring 9, D-70569 Stuttgart, Germany
| | - Øivind Wilhelmsen
- Porelab, Department of Chemistry, Norwegian University of Science and Technology, NO-7491 Trondheim, Norway
| |
Collapse
|
12
|
Heier M, Stephan S, Diewald F, Müller R, Langenbach K, Hasse H. Molecular Dynamics Study of Wetting and Adsorption of Binary Mixtures of the Lennard-Jones Truncated and Shifted Fluid on a Planar Wall. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2021; 37:7405-7419. [PMID: 34097830 PMCID: PMC8280723 DOI: 10.1021/acs.langmuir.1c00780] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 03/19/2021] [Revised: 05/14/2021] [Indexed: 06/12/2023]
Abstract
The wetting of surfaces is strongly influenced by adsorbate layers. Therefore, in this work, sessile drops and their interaction with adsorbate layers on surfaces were investigated by molecular dynamics simulations. Binary fluid model mixtures were considered. The two components of the fluid mixture have the same pure component parameters, but one component has a stronger and the other a weaker affinity to the surface. Furthermore, the unlike interactions between both components were varied. All interactions were described by the Lennard-Jones truncated and shifted potential with a cutoff radius of 2.5σ. The simulations were carried out at constant temperature for mixtures of different compositions. The parameters were varied systematically and chosen such that cases with partial wetting as well as cases with total wetting were obtained and the relation between the varied molecular parameters and the phenomenological behavior was elucidated. Data on the contact angle as well as on the mole fraction and thickness of the adsorbate layer were obtained, accompanied by information on liquid and gaseous bulk phases and the corresponding phase equilibrium. Also, the influence of the adsorbate layer on the wetting was studied: for a sufficiently thick adsorbate layer, the wall's influence on the wetting vanishes, which is then only determined by the adsorbate layer.
Collapse
Affiliation(s)
- Michaela Heier
- Laboratory
of Engineering Thermodynamics, Technische
Universität Kaiserslautern, 67663 Kaiserslautern, Germany
| | - Simon Stephan
- Laboratory
of Engineering Thermodynamics, Technische
Universität Kaiserslautern, 67663 Kaiserslautern, Germany
| | - Felix Diewald
- Institute
of Applied Mechanics, Technische Universität
Kaiserslautern, 67663 Kaiserslautern, Germany
| | - Ralf Müller
- Institute
of Applied Mechanics, Technische Universität
Kaiserslautern, 67663 Kaiserslautern, Germany
| | - Kai Langenbach
- Laboratory
of Engineering Thermodynamics, Technische
Universität Kaiserslautern, 67663 Kaiserslautern, Germany
- Thermal
Separation Science (Endowed Professorship of the State Tyrol), University of Innsbruck, 6020 Innsbruck, Austria
| | - Hans Hasse
- Laboratory
of Engineering Thermodynamics, Technische
Universität Kaiserslautern, 67663 Kaiserslautern, Germany
| |
Collapse
|
13
|
Affiliation(s)
- Simon Homes
- Thermodynamik und Thermische Verfahrenstechnik, Technische Universität Berlin, Berlin, Germany
| | - Matthias Heinen
- Thermodynamik und Thermische Verfahrenstechnik, Technische Universität Berlin, Berlin, Germany
| | - Jadran Vrabec
- Thermodynamik und Thermische Verfahrenstechnik, Technische Universität Berlin, Berlin, Germany
| | - Johann Fischer
- Institut für Verfahrens- und Energietechnik, Universität für Bodenkultur, Wien, Austria
| |
Collapse
|
14
|
Ebrahimi Viand R, Höfling F, Klein R, Delle Site L. Theory and simulation of open systems out of equilibrium. J Chem Phys 2020; 153:101102. [PMID: 32933284 DOI: 10.1063/5.0014065] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023] Open
Abstract
We consider the theoretical model of Bergmann and Lebowitz for open systems out of equilibrium and translate its principles in the adaptive resolution simulation molecular dynamics technique. We simulate Lennard-Jones fluids with open boundaries in a thermal gradient and find excellent agreement of the stationary responses with the results obtained from the simulation of a larger locally forced closed system. The encouraging results pave the way for a computational treatment of open systems far from equilibrium framed in a well-established theoretical model that avoids possible numerical artifacts and physical misinterpretations.
Collapse
Affiliation(s)
- R Ebrahimi Viand
- Freie Universität Berlin, Institute of Mathematics, Arnimallee 6, 14195 Berlin, Germany
| | - F Höfling
- Freie Universität Berlin, Institute of Mathematics, Arnimallee 6, 14195 Berlin, Germany
| | - R Klein
- Freie Universität Berlin, Institute of Mathematics, Arnimallee 6, 14195 Berlin, Germany
| | - L Delle Site
- Freie Universität Berlin, Institute of Mathematics, Arnimallee 6, 14195 Berlin, Germany
| |
Collapse
|
15
|
Stephan S, Schaefer D, Langenbach K, Hasse H. Mass transfer through vapour–liquid interfaces: a molecular dynamics simulation study. Mol Phys 2020. [DOI: 10.1080/00268976.2020.1810798] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
Affiliation(s)
- Simon Stephan
- Laboratory of Engineering Thermodynamics (LTD), TU Kaiserslautern, Kaiserslautern, Germany
| | - Dominik Schaefer
- Laboratory of Engineering Thermodynamics (LTD), TU Kaiserslautern, Kaiserslautern, Germany
| | - Kai Langenbach
- Laboratory of Engineering Thermodynamics (LTD), TU Kaiserslautern, Kaiserslautern, Germany
| | - Hans Hasse
- Laboratory of Engineering Thermodynamics (LTD), TU Kaiserslautern, Kaiserslautern, Germany
| |
Collapse
|
16
|
Stephan S, Deiters UK. Characteristic Curves of the Lennard-Jones Fluid. INTERNATIONAL JOURNAL OF THERMOPHYSICS 2020; 41:147. [PMID: 32863513 PMCID: PMC7441092 DOI: 10.1007/s10765-020-02721-9] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/21/2020] [Accepted: 07/25/2020] [Indexed: 05/25/2023]
Abstract
Equations of state based on intermolecular potentials are often developed about the Lennard-Jones (LJ) potential. Many of such EOS have been proposed in the past. In this work, 20 LJ EOS were examined regarding their performance on Brown's characteristic curves and characteristic state points. Brown's characteristic curves are directly related to the virial coefficients at specific state points, which can be computed exactly from the intermolecular potential. Therefore, also the second and third virial coefficient of the LJ fluid were investigated. This approach allows a comparison of available LJ EOS at extreme conditions. Physically based, empirical, and semi-theoretical LJ EOS were examined. Most investigated LJ EOS exhibit some unphysical artifacts.
Collapse
Affiliation(s)
- Simon Stephan
- Laboratory of Engineering Thermodynamics (LTD), TU Kaiserslautern, Erwin-Schrödinger-Straße 44, 67663 Kaiserslautern, Germany
| | - Ulrich K. Deiters
- Institute of Physical Chemistry, University of Cologne, Greinstraße 4-6, 50939 Cologne, Germany
| |
Collapse
|
17
|
Cea-Klapp E, Quinteros-Lama H, Polishuk I, Garrido JM. Effect of size disparity on the gas-liquid interfacial properties of Lennard-Jones monomer plus dimer mixtures. J Mol Liq 2020. [DOI: 10.1016/j.molliq.2020.113280] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
18
|
Stephan S, Hasse H. Enrichment at vapour–liquid interfaces of mixtures: establishing a link between nanoscopic and macroscopic properties. INT REV PHYS CHEM 2020. [DOI: 10.1080/0144235x.2020.1777705] [Citation(s) in RCA: 25] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Simon Stephan
- Laboratory of Engineering Thermodynamics (LTD), TU Kaiserslautern, Kaiserslautern, Germany
| | - Hans Hasse
- Laboratory of Engineering Thermodynamics (LTD), TU Kaiserslautern, Kaiserslautern, Germany
| |
Collapse
|
19
|
Zhang Y, Xi S, Parambathu AV, Chapman WG. Density functional study of one- and two-component bottlebrush molecules in solvents of varying quality. Mol Phys 2020. [DOI: 10.1080/00268976.2020.1767812] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Affiliation(s)
- Yuchong Zhang
- Department of Chemical and Biomolecular Engineering, Rice University, Houston, TX, USA
| | - Shun Xi
- Department of Chemical and Biomolecular Engineering, Rice University, Houston, TX, USA
| | | | - Walter G. Chapman
- Department of Chemical and Biomolecular Engineering, Rice University, Houston, TX, USA
| |
Collapse
|
20
|
Stephan S, Hasse H. Interfacial properties of binary mixtures of simple fluids and their relation to the phase diagram. Phys Chem Chem Phys 2020; 22:12544-12564. [DOI: 10.1039/d0cp01411g] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
Enrichment at vapour–liquid interfaces can be interpreted as a wetting transition in the vicinity of a three phase equilibrium.
Collapse
Affiliation(s)
- Simon Stephan
- Laboratory of Engineering Thermodynamics (LTD)
- TU Kaiserslautern
- Kaiserslautern 67663
- Germany
| | - Hans Hasse
- Laboratory of Engineering Thermodynamics (LTD)
- TU Kaiserslautern
- Kaiserslautern 67663
- Germany
| |
Collapse
|
21
|
Stephan S, Hasse H. Molecular interactions at vapor-liquid interfaces: Binary mixtures of simple fluids. Phys Rev E 2020; 101:012802. [PMID: 32069593 DOI: 10.1103/physreve.101.012802] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/28/2019] [Indexed: 06/10/2023]
Abstract
Properties of vapor-liquid equilibria and planar interfaces of binary Lennard-Jones truncated and shifted mixtures were investigated with molecular dynamics simulations, density gradient theory, and conformal solution theory at constant liquid phase composition and temperature. The results elucidate the influence of the liquid phase interactions on the interfacial properties (surface tension, surface excess, interfacial thickness, and enrichment). The studied mixtures differ in the ratios of the dispersion energies of the two components ɛ_{2}/ɛ_{1} and the binary interaction parameter ξ. By varying ξ and ɛ_{2}/ɛ_{1}, a variety of types of phase behavior is covered by this paper. The dependence of the interfacial properties on the variables ξ and ɛ_{2}/ɛ_{1} reveals regularities that can be explained by conformal solution theory of the liquid phase. It is thereby shown that the interfacial properties of the mixtures are dominated by the mean liquid phase interactions whereas the vapor phase has only a minor influence.
Collapse
Affiliation(s)
- Simon Stephan
- Laboratory of Engineering Thermodynamics (LTD), TU Kaiserslautern, 67663 Kaiserslautern, Germany
| | - Hans Hasse
- Laboratory of Engineering Thermodynamics (LTD), TU Kaiserslautern, 67663 Kaiserslautern, Germany
| |
Collapse
|
22
|
Stephan S, Dyga M, Urbassek HM, Hasse H. The Influence of Lubrication and the Solid-Fluid Interaction on Thermodynamic Properties in a Nanoscopic Scratching Process. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2019; 35:16948-16960. [PMID: 31815481 DOI: 10.1021/acs.langmuir.9b01033] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Liquid lubricants play an important role in contact processes; for example, they reduce friction and cool the contact zone. To gain better understanding of the influence of lubrication on the nanoscale, both dry and lubricated scratching processes in a model system are compared in the present work using molecular dynamics simulations. The entire range between total dewetting and total wetting is investigated by tuning the solid-fluid interaction energy. The investigated scratching process consists of three sequential movements: A cylindrical indenter penetrates an initially flat substrate, then scratches in the lateral direction, and is finally retracted out of the contact with the substrate. The indenter is fully submersed in the fluid in the lubricated cases. The substrate, the indenter, and the fluid are described by suitably parametrized Lennard-Jones model potentials. The presence of the lubricant is found to have a significant influence on the friction and on the energy balance of the process. The thermodynamic properties of the lubricant are evaluated in detail. A correlation of the simulation results for the profiles of the temperature, density, and pressure of the fluid in the vicinity of the chip is developed. The work done by the indenter is found to mainly dissipate and thereby heat up the substrate and eventually the fluid. Only a minor part of the work causes plastic deformation of the substrate.
Collapse
|
23
|
Stephan S, Hasse H. Influence of dispersive long-range interactions on properties of vapour–liquid equilibria and interfaces of binary Lennard-Jones mixtures. Mol Phys 2019. [DOI: 10.1080/00268976.2019.1699185] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
Affiliation(s)
- Simon Stephan
- Laboratory of Engineering Thermodynamics (LTD), TU Kaiserslautern, Kaiserslautern, Germany
| | - Hans Hasse
- Laboratory of Engineering Thermodynamics (LTD), TU Kaiserslautern, Kaiserslautern, Germany
| |
Collapse
|
24
|
Liu J, Xi S, Chapman WG. Competitive Sorption of CO 2 with Gas Mixtures in Nanoporous Shale for Enhanced Gas Recovery from Density Functional Theory. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2019; 35:8144-8158. [PMID: 31030516 DOI: 10.1021/acs.langmuir.9b00410] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/24/2023]
Abstract
CO2 competitive sorption with shale gas under various conditions from simple to complex pore characteristics is studied using a molecular density functional theory (DFT) that reduces to perturbed chain-statistical associating fluid theory in the bulk fluid region. The DFT model is first verified by grand canonical Monte Carlo simulation in graphite slit pores for pure and binary component systems at different temperatures, pressures, pore sizes, and bulk gas compositions for methane/ethane with CO2. Then, the model is utilized in multicomponent systems that include CH4, C2H6, and C3+ components of different compositions. It is shown that the selectivity of CO2 decreases with increases in temperature, pressure, nanopore size, and average molecular weight of shale gas. Extending the model to more realistic situations, we consider the impact of water present in the pore and consider the effect of permeation of fluid molecules into the kerogen that forms the pore walls. The water-graphite interaction is calibrated with contact angle from molecular simulation data from the literature. The kerogen pore model prediction of gas absolute sorption is compared with experimental and molecular simulation values in the literature. It is shown that the presence of water reduces the CO2 adsorption but improves the CO2 selectivity. The dissolution of gases into the kerogen matrix also leads to the increase in CO2 selectivity. The effect of kerogen type and maturity on the gas sorption amount and CO2 selectivity is also studied. The associated mechanisms are discussed to provide fundamental understanding for gas recovery by CO2.
Collapse
Affiliation(s)
- Jinlu Liu
- Department of Chemical and Biomolecular Engineering , Rice University , 6100 Main Street , Houston , Texas 77005 , United States
| | - Shun Xi
- Department of Chemical and Biomolecular Engineering , Rice University , 6100 Main Street , Houston , Texas 77005 , United States
| | - Walter G Chapman
- Department of Chemical and Biomolecular Engineering , Rice University , 6100 Main Street , Houston , Texas 77005 , United States
| |
Collapse
|
25
|
Stephan S, Langenbach K, Hasse H. Interfacial properties of binary Lennard-Jones mixtures by molecular simulation and density gradient theory. J Chem Phys 2019; 150:174704. [DOI: 10.1063/1.5093603] [Citation(s) in RCA: 35] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023] Open
Affiliation(s)
- Simon Stephan
- Laboratory of Engineering Thermodynamics (LTD), TU Kaiserslautern, Kaiserslautern, Germany
| | - Kai Langenbach
- Laboratory of Engineering Thermodynamics (LTD), TU Kaiserslautern, Kaiserslautern, Germany
| | - Hans Hasse
- Laboratory of Engineering Thermodynamics (LTD), TU Kaiserslautern, Kaiserslautern, Germany
| |
Collapse
|
26
|
Diewald F, Heier M, Horsch M, Kuhn C, Langenbach K, Hasse H, Müller R. Three-dimensional phase field modeling of inhomogeneous gas-liquid systems using the PeTS equation of state. J Chem Phys 2018; 149:064701. [PMID: 30111148 DOI: 10.1063/1.5035495] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Recently, an equation of state (EoS) for the Lennard-Jones truncated and shifted (LJTS) fluid has become available. As it describes metastable and unstable states well, it is suited for predicting density profiles in vapor-liquid interfaces in combination with density gradient theory (DGT). DGT is usually applied to describe interfaces in Cartesian one-dimensional scenarios. In the present work, the perturbed LJ truncated and shifted (PeTS) EoS is implemented into a three-dimensional phase field (PF) model which can be used for studying inhomogeneous gas-liquid systems in a more general way. The results are compared with the results from molecular dynamics simulations for the LJTS fluid that are carried out in the present work and good agreement is observed. The PF model can therefore be used to overcome the scale limit of molecular simulations. A finite element approach is applied for the implementation of the PF model. This requires the first and second derivatives of the PeTS EoS which are calculated using hyper-dual numbers. Several tests and examples of applications of the new PeTS PF model are discussed.
Collapse
Affiliation(s)
- Felix Diewald
- Institute of Applied Mechanics, Technische Universität Kaiserslautern, 67653 Kaiserslautern, Germany
| | - Michaela Heier
- Laboratory of Engineering Thermodynamics, Technische Universität Kaiserslautern, 67653 Kaiserslautern, Germany
| | - Martin Horsch
- Laboratory of Engineering Thermodynamics, Technische Universität Kaiserslautern, 67653 Kaiserslautern, Germany
| | - Charlotte Kuhn
- Computational Mechanics, Technische Universität Kaiserslautern, 67653 Kaiserslautern, Germany
| | - Kai Langenbach
- Laboratory of Engineering Thermodynamics, Technische Universität Kaiserslautern, 67653 Kaiserslautern, Germany
| | - Hans Hasse
- Laboratory of Engineering Thermodynamics, Technische Universität Kaiserslautern, 67653 Kaiserslautern, Germany
| | - Ralf Müller
- Institute of Applied Mechanics, Technische Universität Kaiserslautern, 67653 Kaiserslautern, Germany
| |
Collapse
|