1
|
Heyes DM, Dini D, Pieprzyk S, Brańka AC, Costigliola L. Models to predict configurational adiabats of Lennard-Jones fluids and their transport coefficients. J Chem Phys 2024; 161:084502. [PMID: 39193943 DOI: 10.1063/5.0225650] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/26/2024] [Accepted: 08/08/2024] [Indexed: 08/29/2024] Open
Abstract
A comparison is made between three simple approximate formulas for the configurational adiabat (i.e., constant excess entropy, sex) lines in a Lennard-Jones (LJ) fluid, one of which is an analytic formula based on a harmonic approximation, which was derived by Heyes et al. [J. Chem. Phys. 159, 224504 (2023)] (analytic isomorph line, AIL). Another is where the density is normalized by the freezing density at that temperature (freezing isomorph line, FIL). It is found that the AIL formula and the average of the freezing density and the melting density ("FMIL") are configurational adiabats at all densities essentially down to the liquid-vapor binodal. The FIL approximation departs from a configurational adiabat in the vicinity of the liquid-vapor binodal close to the freezing line. The self-diffusion coefficient, D, shear viscosity, ηs, and thermal conductivity, λ, in macroscopic reduced units are essentially constant along the AIL and FMIL at all fluid densities and temperatures, but departures from this trend are found along the FIL at high liquid state densities near the liquid-vapor binodal. This supports growing evidence that for simple model systems with no or few internal degrees of freedom, isodynes are lines of constant excess entropy. It is shown that for the LJ fluid, ηs and D can be predicted accurately by an essentially analytic procedure from the high temperature limiting inverse power fluid values (apart from at very low densities), and this is demonstrated quite well also for the experimental argon viscosity.
Collapse
Affiliation(s)
- D M Heyes
- Department of Mechanical Engineering, Imperial College London, Exhibition Road, South Kensington, London SW7 2AZ, United Kingdom
| | - D Dini
- Department of Mechanical Engineering, Imperial College London, Exhibition Road, South Kensington, London SW7 2AZ, United Kingdom
| | - S Pieprzyk
- Institute of Molecular Physics, Polish Academy of Sciences, M. Smoluchowskiego 17, 60-179 Poznań, Poland
| | - A C Brańka
- Institute of Molecular Physics, Polish Academy of Sciences, M. Smoluchowskiego 17, 60-179 Poznań, Poland
| | - L Costigliola
- Department of Mechanical Engineering, Imperial College London, Exhibition Road, South Kensington, London SW7 2AZ, United Kingdom
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, DK-4000 Roskilde, Denmark
| |
Collapse
|
2
|
Saric D, Bell IH, Guevara-Carrion G, Vrabec J. Influence of repulsion on entropy scaling and density scaling of monatomic fluids. J Chem Phys 2024; 160:104503. [PMID: 38456532 DOI: 10.1063/5.0196592] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/08/2024] [Accepted: 02/18/2024] [Indexed: 03/09/2024] Open
Abstract
Entropy scaling is applied to the shear viscosity, self-diffusion coefficient, and thermal conductivity of simple monatomic fluids. An extensive molecular dynamics simulation series is performed to obtain these transport properties and the residual entropy of three potential model classes with variable repulsive exponents: n, 6 Mie (n = 9, 12, 15, and 18), Buckingham's exponential-six (α = 12, 14, 18, and 30), and Tang-Toennies (αT = 4.051, 4.275, and 4.600). A wide range of liquid and supercritical gas- and liquid-like states is covered with a total of 1120 state points. Comparisons to equations of state, literature data, and transport property correlations are made. Although the absolute transport property values within a given potential model class may strongly depend on the repulsive exponent, it is found that the repulsive steepness plays a negligible role when entropy scaling is applied. Hence, the plus-scaled transport properties of n, 6 Mie, exponential-six, and Tang-Toennies fluids lie basically on one master curve, which closely corresponds with entropy scaling correlations for the Lennard-Jones fluid. This trend is confirmed by literature data of n, 6 Mie, and exponential-six fluids. Furthermore, entropy scaling holds for state points where the Pearson correlation coefficient R is well below 0.9. The condition R > 0.9 for strongly correlating liquids is thus not necessary for the successful application of entropy scaling, pointing out that isomorph theory may be a part of a more general framework that is behind the success of entropy scaling. Density scaling reveals a strong influence of the repulsive exponent on this particular approach.
Collapse
Affiliation(s)
- Denis Saric
- Thermodynamics, Technical University of Berlin, Ernst-Reuter-Platz 1, 10587 Berlin, Germany
| | - Ian H Bell
- Applied Chemicals and Materials Division, National Institute of Standards and Technology, Boulder, Colorado 80305, USA
| | | | - Jadran Vrabec
- Thermodynamics, Technical University of Berlin, Ernst-Reuter-Platz 1, 10587 Berlin, Germany
| |
Collapse
|
3
|
Chaparro G, Müller EA. Simulation and Data-Driven Modeling of the Transport Properties of the Mie Fluid. J Phys Chem B 2024; 128:551-566. [PMID: 38181201 PMCID: PMC10801693 DOI: 10.1021/acs.jpcb.3c06813] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/13/2023] [Revised: 12/14/2023] [Accepted: 12/15/2023] [Indexed: 01/07/2024]
Abstract
This work reports the computation and modeling of the self-diffusivity (D*), shear viscosity (η*), and thermal conductivity (κ*) of the Mie fluid. The transport properties were computed using equilibrium molecular dynamics simulations for the Mie fluid with repulsive exponents (λr) ranging from 7 to 34 and at a fixed attractive exponent (λa) of 6 over the whole fluid density (ρ*) range and over a wide temperature (T*) range. The computed database consists of 17,212, 14,288, and 13,099 data points for self-diffusivity, shear viscosity, and thermal conductivity, respectively. The database is successfully validated against published simulation data. The above-mentioned transport properties are correlated using artificial neural networks (ANNs). Two modeling approaches were tested: a semiempirical formulation based on entropy scaling and an empirical formulation based on density and temperature as input variables. For the former, it was found that a unique formulation based on entropy scaling does not yield satisfactory results over the entire density range due to a divergent and incorrect scaling of the transport properties at low densities. For the latter empirical modeling approach, it was found that regularizing the data, e.g., modeling ρ*D* instead of D*, ln η* instead of η*, and ln κ* instead of κ*, as well as using the inverse of the temperature as an input feature, helps to ease the interpolation efforts of the artificial neural networks. The trained ANNs can model seen and unseen data over a wide range of density and temperature. Ultimately, the ANNs can be used alongside equations of state to regress effective force field parameters from volumetric and transport data.
Collapse
Affiliation(s)
- Gustavo Chaparro
- Department of Chemical Engineering,
Sargent Centre for Process Systems Engineering, Imperial College London, South Kensington Campus, London SW7 2AZ, U.K.
| | - Erich A. Müller
- Department of Chemical Engineering,
Sargent Centre for Process Systems Engineering, Imperial College London, South Kensington Campus, London SW7 2AZ, U.K.
| |
Collapse
|
4
|
McDonald J, von Spakovsky MR, Reynolds WT. Predicting Polymer Brush Behavior in Solvents Using the Steepest-Entropy-Ascent Quantum Thermodynamic Framework. J Phys Chem B 2023; 127:10370-10391. [PMID: 38006350 DOI: 10.1021/acs.jpcb.3c02713] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2023]
Abstract
The steepest-entropy-ascent quantum thermodynamic (SEAQT) framework is utilized to study the effects of temperature on polymer brushes. The brushes are represented by a discrete energy spectrum, and energy degeneracies obtained through the replica-exchange Wang-Landau algorithm. The SEAQT equation of motion is applied to the density of states to establish a unique kinetic path from an initial thermodynamic state to a stable equilibrium state. The kinetic path describes the brush's evolution in state space, as it interacts with a thermal reservoir. The predicted occupation probabilities along the kinetic path are used to determine the expected thermodynamic and structural properties. The polymer density profile of a polystyrene brush in cyclohexane solvent is predicted using the equation of motion, and it agrees qualitatively with the experimental density profiles. The Flory-Huggins parameter chosen to describe brush-solvent interactions affects the solvent distribution in the brush but has a minimal impact on the polymer density profile. Three types of nonequilibrium kinetic paths with differing amounts of entropy production are considered: a heating path, a cooling path, and a heating-cooling path. Properties such as tortuosity, radius of gyration, brush density, solvent density, and brush chain conformations are calculated for each path.
Collapse
Affiliation(s)
- Jared McDonald
- Materials Science & Engineering Department, Virginia Tech, Blacksburg, Virginia 24061, United States
| | | | - William T Reynolds
- Materials Science & Engineering Department, Virginia Tech, Blacksburg, Virginia 24061, United States
| |
Collapse
|
5
|
Kusudo H, Omori T, Joly L, Yamaguchi Y. The receding contact line cools down during dynamic wetting. J Chem Phys 2023; 159:161102. [PMID: 37877481 DOI: 10.1063/5.0171769] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2023] [Accepted: 10/04/2023] [Indexed: 10/26/2023] Open
Abstract
When a contact line (CL)-where a liquid-vapor interface meets a substrate-is put into motion, it is well known that the contact angle differs between advancing and receding CLs. Using non-equilibrium molecular dynamics simulations, we reveal another intriguing distinction between advancing and receding CLs: while temperature increases at an advancing CL-as expected from viscous dissipation, we show that temperature can drop at a receding CL. Detailed quantitative analysis based on the macroscopic energy balance around the dynamic CL showed that the internal energy change of the fluid due to the change of the potential field along the pathline out of the solid-liquid interface induced a remarkable temperature drop around the receding CL, in a manner similar to latent heat upon phase changes. This result provides new insights for modeling the dynamic CL, and the framework for heat transport analysis introduced here can be applied to a wide range of nanofluidic systems.
Collapse
Affiliation(s)
- Hiroki Kusudo
- Institute of Fluid Science, Tohoku University, 2-1-1 Katahira Aoba-ku, Sendai 980-8577, Japan
| | - Takeshi Omori
- Department of Mechanical Engineering, Osaka Metropolitan University, 3-3-138 Sugimoto, Sumiyoshi, Osaka 558-8585, Japan
| | - Laurent Joly
- Univ Lyon, Univ Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622 Villeurbanne, France
| | - Yasutaka Yamaguchi
- Institute of Fluid Science, Tohoku University, 2-1-1 Katahira Aoba-ku, Sendai 980-8577, Japan
- Water Frontier Research Center (WaTUS), Research Institute for Science & Technology, Tokyo University of Science, 1-3 Kagurazaka, Shinjuku-ku, Tokyo 162-8601, Japan
| |
Collapse
|
6
|
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
|
7
|
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]
|
8
|
Rizk F, Gelin S, Biance AL, Joly L. Microscopic Origins of the Viscosity of a Lennard-Jones Liquid. PHYSICAL REVIEW LETTERS 2022; 129:074503. [PMID: 36018701 DOI: 10.1103/physrevlett.129.074503] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/31/2022] [Accepted: 07/26/2022] [Indexed: 06/15/2023]
Abstract
Unlike crystalline solids or ideal gases, transport properties remain difficult to describe from a microscopic point of view in liquids, whose dynamics result from complex energetic and entropic contributions at the atomic scale. Two scenarios are generally proposed: one represents the dynamics in a fluid as a series of energy-barrier crossings, leading to Arrhenius-like laws, while the other assumes that atoms rearrange themselves by collisions, as exemplified by the free volume model. To assess the validity of these two views, we computed, using molecular dynamics simulations, the transport properties of the Lennard-Jones fluid and tested to what extent the Arrhenius equation and the free volume model describe the temperature dependence of the viscosity and of the diffusion coefficient at fixed pressure. Although both models reproduce the simulation results over a wide range of pressure and temperature covering the liquid and supercritical states of the Lennard-Jones fluid, we found that the parameters of the free volume model can be estimated directly from local structural parameters, also obtained in the simulations. This consistency of the results gives more credibility to the free volume description of transport properties in liquids.
Collapse
Affiliation(s)
- Farid Rizk
- Univ Lyon, Univ Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622 VILLEURBANNE, France
| | - Simon Gelin
- Univ Lyon, Univ Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622 VILLEURBANNE, France
| | - Anne-Laure Biance
- Univ Lyon, Univ Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622 VILLEURBANNE, France
| | - Laurent Joly
- Univ Lyon, Univ Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622 VILLEURBANNE, France
- Institut Universitaire de France (IUF), 1 rue Descartes, 75005 Paris, France
| |
Collapse
|
9
|
Hoang H, Galliero G. Predicting thermodiffusion in simple binary fluid mixtures. THE EUROPEAN PHYSICAL JOURNAL. E, SOFT MATTER 2022; 45:42. [PMID: 35507140 DOI: 10.1140/epje/s10189-022-00197-z] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/01/2022] [Accepted: 04/12/2022] [Indexed: 06/14/2023]
Abstract
The predictive capabilities of some existing theoretical models to quantify thermodiffusion have been investigated in this work. To do so, the tests have been performed on two model fluids, the hard-sphere and the Lennard-Jones (including spheres and dimers) ones, exploring different mixtures and thermodynamic conditions thanks to extensive molecular simulations. It has been confirmed that the thermal diffusion factor should be expressed as the sum of one term related to the isotope effect and one term related to the "chemical" effects and that a kinetic term is required to quantify thermodiffusion from the gas state to the liquid state. In addition, regarding the isotope effects, it has been obtained that none of the available theoretical models are able to yield a reasonable prediction relatively to the molecular simulations results and that the moment of inertia contribution is one order of magnitude smaller than the mass contribution in the liquid state. Finally, concerning the chemical effects, it has been shown the Shukla and Firoozabadi model, complemented with a kinetic term, is probably the most reasonable option to estimate the chemical contribution to the thermal diffusion factor, even if it fails in capturing the effect of the asymmetry in size and in shape between the species. Overall, this works confirms that there is still a lack of a generic model able to predict accurately thermal diffusion factors, or equivalently Soret coefficient, in simple binary mixtures from the gas state to the liquid state.
Collapse
Affiliation(s)
- Hai Hoang
- Institute of Fundamental and Applied Sciences, Duy Tan University, 6 Tran Nhat Duat Street, District 1, Ho Chi Minh City, 700000, Viet Nam
- Faculty of Environmental and Natural Sciences, Duy Tan University, 03 Quang Trung Street, Da Nang, Vietnam
| | - Guillaume Galliero
- Laboratoire des Fluides Complexes et leurs Réservoirs (UMR-5150 with CNRS, and TotalEnergies), Université de Pau et des Pays de l'Adour, BP 1155, 64013, Pau Cedex, France.
| |
Collapse
|
10
|
|
11
|
Molecular simulation of the viscosity of asymmetric dense mixtures. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2021.117052] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|
12
|
Rabani R, Saidi MH, Joly L, Merabia S, Rajabpour A. Enhanced local viscosity around colloidal nanoparticles probed by equilibrium molecular dynamics simulations. J Chem Phys 2021; 155:174701. [PMID: 34742212 DOI: 10.1063/5.0065050] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Nanofluids-dispersions of nanometer-sized particles in a liquid medium-have been proposed for a wide variety of thermal management applications. It is known that a solid-like nanolayer of liquid of typical thicknesses of 0.5-1 nm surrounding the colloidal nanoparticles can act as a thermal bridge between the nanoparticle and the bulk liquid. Yet, its effect on the nanofluid viscosity has not been elucidated so far. In this article, we compute the local viscosity of the nanolayer using equilibrium molecular dynamics based on the Green-Kubo formula. We first assess the validity of the method to predict the viscosity locally. We apply this methodology to the calculation of the local viscosity in the immediate vicinity of a metallic nanoparticle for a wide range of solid-liquid interaction strength, where a nanolayer of thickness 1 nm is observed as a result of the interaction with the nanoparticle. The viscosity of the nanolayer, which is found to be higher than its corresponding bulk value, is directly dependent on the solid-liquid interaction strength. We discuss the origin of this viscosity enhancement and show that the liquid density increment alone cannot explain the values of the viscosity observed. Rather, we suggest that the solid-like structure of the distribution of the liquid atoms in the vicinity of the nanoparticle contributes to the nanolayer viscosity enhancement. Finally, we observe a failure of the Stokes-Einstein relation between viscosity and diffusion close to the wall, depending on the liquid-solid interaction strength, which we rationalize in terms of the hydrodynamic slip.
Collapse
Affiliation(s)
- Reza Rabani
- Center of Excellence in Energy Conversion (CEEC), School of Mechanical Engineering, Sharif University of Technology, Tehran 11155-9567, Iran
| | - Mohammad Hassan Saidi
- Center of Excellence in Energy Conversion (CEEC), School of Mechanical Engineering, Sharif University of Technology, Tehran 11155-9567, Iran
| | - Laurent Joly
- Univ Lyon, Univ Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622 Villeurbanne, France
| | - Samy Merabia
- Univ Lyon, Univ Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622 Villeurbanne, France
| | - Ali Rajabpour
- Advanced Simulation and Computing Laboratory (ASCL), Mechanical Engineering Department, Imam Khomeini International University, Qazvin, Iran
| |
Collapse
|
13
|
Wang J, Zhang L, Li R, Chen X, Dong X, Zhao Y, Zhong Q, Sheng B, Li H, Gong M. A simple generalized viscosity model for natural gas components and their mixtures based on free-volume theory. J Mol Liq 2020. [DOI: 10.1016/j.molliq.2020.114483] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
14
|
Binti Mohd Taib M, Trusler JPM. Residual entropy model for predicting the viscosities of dense fluid mixtures. J Chem Phys 2020; 152:164104. [PMID: 32357798 DOI: 10.1063/5.0002242] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
In this work, we have investigated the mono-variant relationship between the reduced viscosity and residual entropy in pure fluids and in binary mixtures of hydrocarbons and hydrocarbons with dissolved carbon dioxide. The mixtures considered were octane + dodecane, decane + carbon dioxide, and 1,3-dimethylbenzene (m-xylene) + carbon dioxide. The reduced viscosity was calculated according to the definition of Bell, while the residual entropy was calculated from accurate multi-parameter Helmholtz-energy equations of state and, for mixtures, the multi-fluid Helmholtz energy approximation. The mono-variant dependence of reduced viscosity upon residual molar entropy was observed for the pure fluids investigated, and by incorporating two scaling factors (one for reduced viscosity and the other for residual molar entropy), the data were represented by a single universal curve. To apply this method to mixtures, the scaling factors were determined from a mole-fraction weighted sum of the pure-component values. This simple model was found to work well for the systems investigated. The average absolute relative deviation (AARD) was observed to be between 1% and 2% for pure components and a mixture of similar hydrocarbons. Larger deviations, with AARDs of up to 15%, were observed for the asymmetric mixtures, but this compares favorably with other methods for predicting the viscosity of such systems. We conclude that the residual-entropy concept can be used to estimate the viscosity of mixtures of similar molecules with high reliability and that it offers a useful engineering approximation even for asymmetric mixtures.
Collapse
Affiliation(s)
- Malyanah Binti Mohd Taib
- Department of Chemical Engineering, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom
| | - J P Martin Trusler
- Department of Chemical Engineering, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom
| |
Collapse
|
15
|
Bell IH, Messerly R, Thol M, Costigliola L, Dyre JC. Modified Entropy Scaling of the Transport Properties of the Lennard-Jones Fluid. J Phys Chem B 2019; 123:6345-6363. [PMID: 31241958 DOI: 10.1021/acs.jpcb.9b05808] [Citation(s) in RCA: 53] [Impact Index Per Article: 10.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Rosenfeld proposed two different scaling approaches to model the transport properties of fluids, separated by 22 years, one valid in the dilute gas, and another in the liquid phase. In this work, we demonstrate that these two limiting cases can be connected through the use of a novel approach to scaling transport properties and a bridging function. This approach, which is empirical and not derived from theory, is used to generate reference correlations for the transport properties of the Lennard-Jones 12-6 fluid of viscosity, thermal conductivity, and self-diffusion. This approach, with a very simple functional form, allows for the reproduction of the most accurate simulation data to within nearly their statistical uncertainty. The correlations are used to confirm that for the Lennard-Jones fluid the appropriately scaled transport properties are nearly monovariate functions of the excess entropy from low-density gases into the supercooled phase and up to extreme temperatures. This study represents the most comprehensive metastudy of the transport properties of the Lennard-Jones fluid to date.
Collapse
Affiliation(s)
- Ian H Bell
- Applied Chemicals and Materials Division , National Institute of Standards and Technology , Boulder , Colorado 80305 , United States
| | - Richard Messerly
- Applied Chemicals and Materials Division , National Institute of Standards and Technology , Boulder , Colorado 80305 , United States
| | - Monika Thol
- Thermodynamics , Ruhr-Universität Bochum , Universitätsstraße 150 , 44801 Bochum , Germany
| | - Lorenzo Costigliola
- DNRF Centre "Glass and Time," IMFUFA, Department of Science and Environment , Roskilde University , Postbox 260, DK-4000 Roskilde , Denmark
| | - Jeppe C Dyre
- DNRF Centre "Glass and Time," IMFUFA, Department of Science and Environment , Roskilde University , Postbox 260, DK-4000 Roskilde , Denmark
| |
Collapse
|
16
|
Russo A, Durán-Olivencia MA, Kalliadasis S, Hartkamp R. Macroscopic relations for microscopic properties at the interface between solid substrates and dense fluids. J Chem Phys 2019; 150:214705. [DOI: 10.1063/1.5094911] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Affiliation(s)
- Antonio Russo
- Department of Chemical Engineering, Imperial College London, London SW7 2WZ, United Kingdom
| | | | - Serafim Kalliadasis
- Department of Chemical Engineering, Imperial College London, London SW7 2WZ, United Kingdom
| | - Remco Hartkamp
- Process and Energy Department, Delft University of Technology, Leeghwaterstraat 39, 2628 CB Delft, The Netherlands
| |
Collapse
|
17
|
Giunta G, Carbone P. Cross-over in the dynamics of polymer confined between two liquids of different viscosity. Interface Focus 2019; 9:20180074. [PMID: 31065342 PMCID: PMC6501349 DOI: 10.1098/rsfs.2018.0074] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 02/27/2019] [Indexed: 11/12/2022] Open
Abstract
Using molecular dynamics simulations, we analysed the polymer dynamics of chains of different molecular weights entrapped at the interface between two immiscible liquids. We showed that on increasing the viscosity of one of the two liquids the dynamic behaviour of the chain changes from a Zimm-like dynamics typical of dilute polymer solutions to a Rouse-like dynamics where hydrodynamic interactions are screened. We observed that when the polymer is in contact with a high viscosity liquid, the number of solvent molecules close to the polymer beads is reduced and ascribed the screening effect to this reduced number of polymer-solvent contacts. For the longest chain simulated, we calculated the distribution of loop length and compared the results with the theoretical distribution developed for solid/liquid interfaces. We showed that the polymer tends to form loops (although flat against the interface) and that the theory works reasonably well also for liquid/liquid interfaces.
Collapse
Affiliation(s)
- Giuliana Giunta
- School of Chemical Engineering and Analytical Science, University of Manchester, Oxford Road, Manchester M13 9PL, UK
| | - Paola Carbone
- School of Chemical Engineering and Analytical Science, University of Manchester, Oxford Road, Manchester M13 9PL, UK
| |
Collapse
|
18
|
Song X, Qiao C, Tao J, Bao B, Han X, Zhao S. Interfacial Engineering of Thermoresponsive Microgel Capsules: Polymeric Wetting vs Colloidal Adhesion. Macromolecules 2019. [DOI: 10.1021/acs.macromol.8b02323] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
|
19
|
Probing the link between residual entropy and viscosity of molecular fluids and model potentials. Proc Natl Acad Sci U S A 2019; 116:4070-4079. [PMID: 30770449 DOI: 10.1073/pnas.1815943116] [Citation(s) in RCA: 33] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
This work investigates the link between residual entropy and viscosity based on wide-ranging, highly accurate experimental and simulation data. This link was originally postulated by Rosenfeld in 1977 [Rosenfeld Y (1977) Phys Rev A 15:2545-2549], and it is shown that this scaling results in an approximately monovariate relationship between residual entropy and reduced viscosity for a wide range of molecular fluids [argon, methane, [Formula: see text], [Formula: see text], refrigerant R-134a (1,1,1,2-tetrafluoroethane), refrigerant R-125 (pentafluoroethane), methanol, and water] and a range of model potentials (hard sphere, inverse power, Lennard-Jones, and Weeks-Chandler-Andersen). While the proposed "universal" correlation of Rosenfeld is shown to be far from universal, when used with the appropriate density scaling for molecular fluids, the viscosity of nonassociating molecular fluids can be mapped onto the model potentials. This mapping results in a length scale that is proportional to the cube root of experimentally measurable liquid volume values.
Collapse
|
20
|
Mozaffari SH, Srinivasan S, Saghir MZ. A Study on Thermodiffusion in Ternary Liquid Mixtures Using Enhanced Molecular Dynamics Algorithm with Experimental Validation. CAN J CHEM ENG 2019. [DOI: 10.1002/cjce.23199] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Affiliation(s)
- Seyedeh Hoda Mozaffari
- Department of Mech. & Indust. EngineeringRyerson University350 Victoria St.TorontoON, M5B 2K3Canada
| | - Seshasai Srinivasan
- Department of Mech. & Indust. EngineeringRyerson University350 Victoria St.TorontoON, M5B 2K3Canada
- W Booth School of Engineering Pract. & Tech., 1280 Main St. WestMcMaster UniversityHamiltonON, L8S 4L8Canada
- Department of Mech. EngineeringMcMaster UniversityHamiltonON, L8S 4L8Canada
| | - M. Ziad Saghir
- Department of Mech. & Indust. EngineeringRyerson University350 Victoria St.TorontoON, M5B 2K3Canada
| |
Collapse
|
21
|
Jamali SH, Hartkamp R, Bardas C, Söhl J, Vlugt TJH, Moultos OA. Shear Viscosity Computed from the Finite-Size Effects of Self-Diffusivity in Equilibrium Molecular Dynamics. J Chem Theory Comput 2018; 14:5959-5968. [PMID: 30296092 PMCID: PMC6236468 DOI: 10.1021/acs.jctc.8b00625] [Citation(s) in RCA: 29] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
A method is proposed for calculating
the shear viscosity of a liquid
from finite-size effects of self-diffusion coefficients in Molecular
Dynamics simulations. This method uses the difference in the self-diffusivities,
computed from at least two system sizes, and an analytic equation
to calculate the shear viscosity. To enable the efficient use of this
method, a set of guidelines is developed. The most efficient number
of system sizes is two and the large system is at least four times
the small system. The number of independent simulations for each system
size should be assigned in such a way that 50%–70% of the total
available computational resources are allocated to the large system.
We verified the method for
250 binary and 26 ternary Lennard-Jones systems, pure water, and an
ionic liquid ([Bmim][Tf2N]). The computed shear viscosities
are in good agreement with viscosities obtained from equilibrium Molecular
Dynamics simulations for all liquid systems far from the critical
point. Our results indicate that the proposed method is suitable for
multicomponent mixtures and highly viscous liquids. This may enable
the systematic screening of the viscosities of ionic liquids and deep
eutectic solvents.
Collapse
Affiliation(s)
- Seyed Hossein Jamali
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering , Delft University of Technology , Leeghwaterstraat 39 , 2628CB Delft , The Netherlands
| | - Remco Hartkamp
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering , Delft University of Technology , Leeghwaterstraat 39 , 2628CB Delft , The Netherlands
| | - Christos Bardas
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering , Delft University of Technology , Leeghwaterstraat 39 , 2628CB Delft , The Netherlands
| | - Jakob Söhl
- Delft Institute of Applied Mathematics , Delft University of Technology , van Mourik Broekmanweg 6 , 2628XE Delft , The Netherlands
| | - Thijs J H Vlugt
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering , Delft University of Technology , Leeghwaterstraat 39 , 2628CB Delft , The Netherlands
| | - Othonas A Moultos
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering , Delft University of Technology , Leeghwaterstraat 39 , 2628CB Delft , The Netherlands
| |
Collapse
|
22
|
Lautenschlaeger MP, Horsch M, Hasse H. Simultaneous determination of thermal conductivity and shear viscosity using two-gradient non-equilibrium molecular dynamics simulations. Mol Phys 2018. [DOI: 10.1080/00268976.2018.1504134] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
Affiliation(s)
- Martin P. Lautenschlaeger
- Laboratory of Engineering Thermodynamics (LTD), University of Kaiserslautern, Kaiserslautern, Germany
| | - Martin Horsch
- Laboratory of Engineering Thermodynamics (LTD), University of Kaiserslautern, Kaiserslautern, Germany
| | - Hans Hasse
- Laboratory of Engineering Thermodynamics (LTD), University of Kaiserslautern, Kaiserslautern, Germany
| |
Collapse
|
23
|
Tu KM, Kim K, Matubayasi N. Spatial-decomposition analysis of viscosity with application to Lennard-Jones fluid. J Chem Phys 2018. [DOI: 10.1063/1.5018483] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Affiliation(s)
- Kai-Min Tu
- Elements Strategy Initiative for Catalysts and Batteries, Kyoto University, Katsura, Kyoto 615-8520, Japan
| | - Kang Kim
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| | - Nobuyuki Matubayasi
- Elements Strategy Initiative for Catalysts and Batteries, Kyoto University, Katsura, Kyoto 615-8520, Japan
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| |
Collapse
|
24
|
Lötgering-Lin O, Fischer M, Hopp M, Gross J. Pure Substance and Mixture Viscosities Based on Entropy Scaling and an Analytic Equation of State. Ind Eng Chem Res 2018. [DOI: 10.1021/acs.iecr.7b04871] [Citation(s) in RCA: 55] [Impact Index Per Article: 9.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Oliver Lötgering-Lin
- Institute of Thermodynamics and Thermal Process Engineering, University of Stuttgart, Pfaffenwaldring 9, 70569 Stuttgart, Germany
| | - Matthias Fischer
- Institute of Thermodynamics and Thermal Process Engineering, University of Stuttgart, Pfaffenwaldring 9, 70569 Stuttgart, Germany
| | - Madlen Hopp
- Institute of Thermodynamics and Thermal Process Engineering, University of Stuttgart, Pfaffenwaldring 9, 70569 Stuttgart, Germany
| | - Joachim Gross
- Institute of Thermodynamics and Thermal Process Engineering, University of Stuttgart, Pfaffenwaldring 9, 70569 Stuttgart, Germany
| |
Collapse
|
25
|
Ramírez-González PV, Escobar-Barrios VA. Viscosity and normal stress forces of Lennard-Jones chains using reverse non-equilibrium molecular dynamics. Mol Phys 2017. [DOI: 10.1080/00268976.2017.1340682] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Affiliation(s)
- Patsy V. Ramírez-González
- Instituto Potosino de Investigación Científica y Tecnológica, San Luis Potosí, Mexico, Materiales Avanzados
| | | |
Collapse
|
26
|
Abstract
We report here on the computation of the microscopic flow induced by the motion of a small tagged particle in a fluid from molecular dynamic simulations. It is found that the hydrodynamical Stokes solution with slip boundary conditions is recovered at only a few diameters away from the tagged particle. However, fluctuations of the diffusing particle itself induce a renormalization of the bath viscosity and, more strikingly, an apparent violation of the non-penetrability of the particles in the laboratory frame. The expected zero normal velocity at contact is satisfied only in the particle frame, or for heavy particles. Further evidence of this generalized boundary condition is given by the evaluation of the flow in a granular gas using data from particle tracking experiments.
Collapse
Affiliation(s)
- Dominika Lesnicki
- Département de Chimie, École Normale Supérieure, PSL Research University, Sorbonne Universités-UPMC Univ Paris 06, CNRS UMR 8640 PASTEUR, 24 Rue Lhomond, 75005 Paris, France
| | - Rodolphe Vuilleumier
- Département de Chimie, École Normale Supérieure, PSL Research University, Sorbonne Universités-UPMC Univ Paris 06, CNRS UMR 8640 PASTEUR, 24 Rue Lhomond, 75005 Paris, France
| |
Collapse
|
27
|
Morciano M, Fasano M, Nold A, Braga C, Yatsyshin P, Sibley DN, Goddard BD, Chiavazzo E, Asinari P, Kalliadasis S. Nonequilibrium molecular dynamics simulations of nanoconfined fluids at solid-liquid interfaces. J Chem Phys 2017; 146:244507. [DOI: 10.1063/1.4986904] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023] Open
Affiliation(s)
- M. Morciano
- Energy Department, Politecnico di Torino, Corso Duca degli Abruzzi 24, Torino 10129, Italy
| | - M. Fasano
- Energy Department, Politecnico di Torino, Corso Duca degli Abruzzi 24, Torino 10129, Italy
| | - A. Nold
- Department of Chemical Engineering, Imperial College London, London SW7 2AZ, United Kingdom
| | - C. Braga
- Department of Chemical Engineering, Imperial College London, London SW7 2AZ, United Kingdom
| | - P. Yatsyshin
- Department of Chemical Engineering, Imperial College London, London SW7 2AZ, United Kingdom
| | - D. N. Sibley
- Department of Mathematical Sciences, Loughborough University, Loughborough LE11 3TU, United Kingdom
| | - B. D. Goddard
- The School of Mathematics and Maxwell Institute for Mathematical Sciences, The University of Edinburgh, Edinburgh EH9 3JZ, United Kingdom
| | - E. Chiavazzo
- Energy Department, Politecnico di Torino, Corso Duca degli Abruzzi 24, Torino 10129, Italy
| | - P. Asinari
- Energy Department, Politecnico di Torino, Corso Duca degli Abruzzi 24, Torino 10129, Italy
| | - S. Kalliadasis
- Department of Chemical Engineering, Imperial College London, London SW7 2AZ, United Kingdom
| |
Collapse
|
28
|
Castro-Palacio JC, Hellmann R, Vesovic V. Dilute gas viscosity of n-alkanes represented by rigid Lennard-Jones chains. Mol Phys 2016. [DOI: 10.1080/00268976.2016.1222456] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
Affiliation(s)
| | | | - Velisa Vesovic
- Department of Earth Science and Engineering, Imperial College London, London, United Kingdom
| |
Collapse
|
29
|
Bhadauria R, Aluru NR. A multiscale transport model for Lennard-Jones binary mixtures based on interfacial friction. J Chem Phys 2016; 145:074115. [DOI: 10.1063/1.4961226] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Ravi Bhadauria
- Department of Mechanical Science and Engineering, Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| | - N. R. Aluru
- Department of Mechanical Science and Engineering, Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| |
Collapse
|
30
|
Nichele J, Borges I, Oliveira AB, Alves LS. Molecular dynamics simulations of momentum and thermal diffusion properties of near-critical argon along isobars. J Supercrit Fluids 2016. [DOI: 10.1016/j.supflu.2016.04.004] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
31
|
Theoretical predictions of viscosity of methane under confined conditions. Chin J Chem Eng 2016. [DOI: 10.1016/j.cjche.2016.04.009] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
32
|
Misin M, Palmer DS, Fedorov MV. Predicting Solvation Free Energies Using Parameter-Free Solvent Models. J Phys Chem B 2016; 120:5724-31. [DOI: 10.1021/acs.jpcb.6b05352] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Maksim Misin
- Department
of Physics, SUPA, University of Strathclyde, 107 Rottenrow, Glasgow, G4 0NG, U.K
| | - David S. Palmer
- Department
of Pure and Applied Chemistry, University of Strathclyde, 295 Cathedral
Street, Glasgow, G1 1XL, U.K
| | - Maxim V. Fedorov
- Department
of Physics, SUPA, University of Strathclyde, 107 Rottenrow, Glasgow, G4 0NG, U.K
- Skolkovo Institute of Science and Technology, 3 Nobel Street, Moscow 143026, Russian Federation
| |
Collapse
|
33
|
|
34
|
Delage-Santacreu S, Galliero G, Hoang H, Bazile JP, Boned C, Fernandez J. Thermodynamic scaling of the shear viscosity of Mie n-6 fluids and their binary mixtures. J Chem Phys 2015; 142:174501. [DOI: 10.1063/1.4919296] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Stephanie Delage-Santacreu
- Laboratoire de Mathématiques et leurs Applications (UMR-5142 with CNRS), Université de Pau et des Pays de l’Adour, BP 1155, F-64013 PAU Cedex, France
| | - Guillaume Galliero
- Laboratoire des Fluides Complexes et leurs Reservoirs (UMR-5150 with CNRS and TOTAL), Université de Pau et des Pays de l’Adour, BP 1155, F-64013 PAU Cedex, France
| | - Hai Hoang
- Laboratoire des Fluides Complexes et leurs Reservoirs (UMR-5150 with CNRS and TOTAL), Université de Pau et des Pays de l’Adour, BP 1155, F-64013 PAU Cedex, France
| | - Jean-Patrick Bazile
- Laboratoire des Fluides Complexes et leurs Reservoirs (UMR-5150 with CNRS and TOTAL), Université de Pau et des Pays de l’Adour, BP 1155, F-64013 PAU Cedex, France
| | - Christian Boned
- Laboratoire des Fluides Complexes et leurs Reservoirs (UMR-5150 with CNRS and TOTAL), Université de Pau et des Pays de l’Adour, BP 1155, F-64013 PAU Cedex, France
| | - Josefa Fernandez
- Laboratorio de Propiedades Termofisicas, Universidade Santiago de Compostela, Campus Vida, E-15782 Santiago de Compostela, Spain
| |
Collapse
|
35
|
Angélil R, Diemand J, Tanaka KK, Tanaka H. Bubble evolution and properties in homogeneous nucleation simulations. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2014; 90:063301. [PMID: 25615216 DOI: 10.1103/physreve.90.063301] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/01/2014] [Indexed: 06/04/2023]
Abstract
We analyze the properties of naturally formed nanobubbles in Lennard-Jones molecular dynamics simulations of liquid-to-vapor nucleation in the boiling and the cavitation regimes. The large computational volumes provide a realistic environment at unchanging average temperature and liquid pressure, which allows us to accurately measure properties of bubbles from their inception as stable, critically sized bubbles, to their continued growth into the constant speed regime. Bubble gas densities are up to 50% lower than the equilibrium vapor densities at the liquid temperature, yet quite close to the gas equilibrium density at the lower gas temperatures measured in the simulations: The latent heat of transformation results in bubble gas temperatures up to 25% below those of the surrounding bulk liquid. In the case of rapid bubble growth-typical for the cavitation regime-compression of the liquid outside the bubble leads to local temperature increases of up to 5%, likely significant enough to alter the surface tension as well as the local viscosity. The liquid-vapor bubble interface is thinner than expected from planar coexistence simulations by up to 50%. Bubbles near the critical size are extremely nonspherical, yet they quickly become spherical as they grow. The Rayleigh-Plesset description of bubble-growth gives good agreement in the cavitation regime.
Collapse
Affiliation(s)
- Raymond Angélil
- Institute for Computational Science, University of Zurich, 8057 Zurich, Switzerland
| | - Jürg Diemand
- Institute for Computational Science, University of Zurich, 8057 Zurich, Switzerland
| | - Kyoko K Tanaka
- Institute of Low Temperature Science, Hokkaido University, Sapporo 060-0819, Japan
| | - Hidekazu Tanaka
- Institute of Low Temperature Science, Hokkaido University, Sapporo 060-0819, Japan
| |
Collapse
|
36
|
Galliero G. Equilibrium, interfacial and transport properties of n-alkanes: Towards the simplest coarse grained molecular model. Chem Eng Res Des 2014. [DOI: 10.1016/j.cherd.2014.05.028] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
37
|
Hoang H, Galliero G. Local shear viscosity of strongly inhomogeneous dense fluids: from the hard-sphere to the Lennard-Jones fluids. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2013; 25:485001. [PMID: 24132101 DOI: 10.1088/0953-8984/25/48/485001] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/02/2023]
Abstract
This work aims at providing a tractable approach to model the local shear viscosity of strongly inhomogeneous dense fluids composed of spherical molecules, in which the density variations occur on molecular distance. The proposed scheme, which relies on the local density average model, has been applied to the quasi-hard-sphere, the Week-Chandler-Andersen and the Lennard-Jones fluids. A weight function has been developed to deal with the hard-sphere fluid given the specificities of momentum exchange. To extend the approach to the smoothly repulsive potential, we have taken into account that the non-local contributions to the viscosity due to the interactions of particles separated by a given distance are temperature dependent. Then, using a simple perturbation scheme, the approach is extended to the Lennard-Jones fluids. It is shown that the viscosity profiles of inhomogeneous dense fluids deduced from this approach are consistent with those directly computed by non-equilibrium molecular dynamics simulations.
Collapse
Affiliation(s)
- Hai Hoang
- Laboratoire des Fluides Complexes et leurs Réservoirs (UMR-5150 with CNRS and TOTAL), Université de Pau et des Pays de l'Adour, BP 1155, F-64013 PAU Cedex, France
| | | |
Collapse
|
38
|
Contribution to the modeling of the shear viscosity of sulfur hexafluoride (SF6): Comparative study of some representative models. Chem Phys 2013. [DOI: 10.1016/j.chemphys.2013.06.017] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
39
|
Llovell F, Marcos RM, Vega LF. Free-Volume Theory Coupled with Soft-SAFT for Viscosity Calculations: Comparison with Molecular Simulation and Experimental Data. J Phys Chem B 2013; 117:8159-71. [DOI: 10.1021/jp401307t] [Citation(s) in RCA: 86] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- F. Llovell
- MATGAS Research Center, Campus de la UAB, 08193 Bellaterra, Barcelona, Spain
| | - R. M. Marcos
- Departament d’Enginyeria
Mecànica, ETSE, Universitat Rovira i Virgili, 43007 Tarragona, Spain
| | - L. F. Vega
- MATGAS Research Center, Campus de la UAB, 08193 Bellaterra, Barcelona, Spain
- Carburos Metálicos/Air Products Group, C/Aragón, 300, 08009 Barcelona,
Spain
| |
Collapse
|
40
|
Detcheverry F, Bocquet L. Thermal fluctuations of hydrodynamic flows in nanochannels. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2013; 88:012106. [PMID: 23944413 DOI: 10.1103/physreve.88.012106] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/22/2013] [Indexed: 06/02/2023]
Abstract
Flows at the nanoscale are subject to thermal fluctuations. In this work, we explore the consequences for a fluid confined within a channel of nanometric size. First, the phenomenon is illustrated on the basis of molecular dynamics simulations. The center of mass of the confined fluid is shown to perform a stochastic, non-Markovian motion, whose diffusion coefficient satisfies Einstein's relation, and which can be further characterized by the fluctuation relation. Next, we develop an analytical description of the thermally induced fluid motion. We compute the time- and space-dependent velocity correlation function, and characterize its dependence on the nanopore shape, size, and boundary slip at the surface. The experimental implications for mass and charge transports are discussed for two situations. For a particle confined within the nanopore, we show that the fluid fluctuating motion results in an enhanced diffusion. The second situation involves a charged nanopore in which fluid motion within the double layer induces a fluctuating electric current. We compute the corresponding contribution to the current power spectrum.
Collapse
Affiliation(s)
- François Detcheverry
- Institut Lumière Matière, UMR5306 Université Lyon 1-CNRS, Université de Lyon 69622 Villeurbanne, France
| | | |
Collapse
|
41
|
Llovell F, Marcos RM, Vega LF. Transport Properties of Mixtures by the Soft-SAFT + Free-Volume Theory: Application to Mixtures of n-Alkanes and Hydrofluorocarbons. J Phys Chem B 2013; 117:5195-205. [DOI: 10.1021/jp401754r] [Citation(s) in RCA: 54] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Affiliation(s)
- F. Llovell
- MATGAS Research Center (Carburos Metálicos/Air Products, CSIC, UAB), Campus de la UAB, 08193 Bellaterra, Barcelona, Spain
| | - R. M. Marcos
- Departament d’Enginyeria
Mecànica, ETSE, Universitat Rovira i Virgili, 43007 Tarragona, Spain
| | - L. F. Vega
- MATGAS Research Center (Carburos Metálicos/Air Products, CSIC, UAB), Campus de la UAB, 08193 Bellaterra, Barcelona, Spain
- Carburos Metálicos/Air Products Group, C/Aragón, 300, 08069 Barcelona,
Spain
| |
Collapse
|
42
|
Hellmann R, Riesco N, Vesovic V. Calculation of the transport properties of a dilute gas consisting of Lennard-Jones chains. J Chem Phys 2013; 138:084309. [PMID: 23464153 DOI: 10.1063/1.4793221] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The transport properties in the dilute gas limit have been calculated by the classical-trajectory method for a gas consisting of chain-like molecules. The molecules were modelled as rigid chains consisting of spherical segments that interact through a combination of site-site Lennard-Jones 12-6 potentials. Results are reported for shear viscosity, self-diffusion, and thermal conductivity for chains consisting of 1, 2, 3, 4, 5, 6, 7, 8, 10, 13, and 16 segments in the reduced temperature range of 0.3-50. The results indicate that the transport properties increase with temperature and decrease with chain length. At high temperatures the dependence of the transport properties is governed effectively by the repulsive part of the potential. No simple scaling with chain length has been observed. The higher order correction factors are larger than observed for real molecules so far, reaching asymptotic values of 1.019-1.033 and 1.060-1.072 for viscosity and thermal conductivity, respectively. The dominant contribution comes from the angular momentum coupling. The agreement with molecular dynamics calculations for viscosity is within the estimated accuracy of the two methods for shorter chains. However, for longer chains differences of up to 7% are observed.
Collapse
Affiliation(s)
- Robert Hellmann
- Institut für Chemie, Universität Rostock, 18059 Rostock, Germany
| | | | | |
Collapse
|
43
|
Hoang H, Galliero G. Shear behavior of a confined thin film: Influence of the molecular dynamics scheme employed. J Chem Phys 2013; 138:054707. [DOI: 10.1063/1.4789582] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
44
|
Meng X, Wang Y, Zhao Y, Huang J. Distinct transport properties of O2and CH4across a carbon nanotube. Mol Phys 2013. [DOI: 10.1080/00268976.2012.762126] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|
45
|
Delage Santacreu S, Galliero G, Odunlami M, Boned C. Low density shear viscosity of Lennard-Jones chains of variable rigidities. J Chem Phys 2012. [DOI: 10.1063/1.4767528] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023] Open
|
46
|
Baidakov VG, Protsenko SP, Kozlova ZR. Metastable Lennard-Jones fluids. I. Shear viscosity. J Chem Phys 2012; 137:164507. [DOI: 10.1063/1.4758806] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
47
|
Hoang H, Galliero G. Local viscosity of a fluid confined in a narrow pore. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2012; 86:021202. [PMID: 23005753 DOI: 10.1103/physreve.86.021202] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/09/2012] [Indexed: 06/01/2023]
Abstract
In this paper, molecular dynamics simulations of a simple Lennard-Jones fluid confined in narrow slit pores and undergoing shear have been performed. The aim is to investigate the effects of density inhomogeneities at the fluid-solid interfaces on the shear viscosity profiles. It has been found that the local viscosity was varying strongly with the distance from the solid walls for both dilute and dense fluid states with oscillations correlated to the density ones. To describe the computed viscosity profiles, we propose a scheme that uses the local average density model, combined with an adequate weight function, for the configurational viscosity and a semiempirical model for the translational viscosity. It is shown that the proposed approach is able to provide viscosity profiles in good agreement with those coming from simulations for different pore widths and for different fluid states (dilute to dense).
Collapse
Affiliation(s)
- Hai Hoang
- Laboratoire des Fluides Complexes et leurs Réservoirs, UMR-5150 with CNRS and TOTAL, Université de Pau et des Pays de l'Adour, Boîte Postale 1155, PAU Cedex 64013, France
| | | |
Collapse
|
48
|
Yong X, Zhang LT. Nanoscale simple-fluid behavior under steady shear. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2012; 85:051202. [PMID: 23004740 DOI: 10.1103/physreve.85.051202] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/17/2012] [Revised: 04/18/2012] [Indexed: 06/01/2023]
Abstract
In this study, we use two nonequilibrium molecular dynamics algorithms, boundary-driven shear and homogeneous shear, to explore the rheology and flow properties of a simple fluid undergoing steady simple shear. The two distinct algorithms are designed to elucidate the influences of nanoscale confinement. The results of rheological material functions, i.e., viscosity and normal pressure differences, show consistent Newtonian behaviors at low shear rates from both systems. The comparison validates that confinements of the order of 10 nm are not strong enough to deviate the simple fluid behaviors from the continuum hydrodynamics. The non-Newtonian phenomena of the simple fluid are further investigated by the homogeneous shear simulations with much higher shear rates. We observe the "string phase" at high shear rates by applying both profile-biased and profile-unbiased thermostats. Contrary to other findings where the string phase is found to be an artifact of the thermostats, we perform a thorough analysis of the fluid microstructures formed due to shear, which shows that it is possible to have a string phase and second shear thinning for dense simple fluids.
Collapse
Affiliation(s)
- Xin Yong
- Department of Mechanical, Aerospace & Nuclear Engineering, Rensselaer Polytechnic Institute, Troy, New York 12180, USA
| | | |
Collapse
|
49
|
|
50
|
de Wijn AS, Riesco N, Jackson G, Martin Trusler JP, Vesovic V. Viscosity of liquid mixtures: The Vesovic-Wakeham method for chain molecules. J Chem Phys 2012; 136:074514. [DOI: 10.1063/1.3685605] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
|