1
|
Malgaretti P, Bafile U, Vallauri R, Jedlovszky P, Sega M. Surface viscosity in simple liquids. J Chem Phys 2023; 158:114705. [PMID: 36948818 DOI: 10.1063/5.0141971] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/18/2023] Open
Abstract
The response of Newtonian liquids to small perturbations is usually considered to be fully described by homogeneous transport coefficients like shear and dilatational viscosity. However, the presence of strong density gradients at the liquid/vapor boundary of fluids hints at the possible existence of an inhomogeneous viscosity. Here, we show that a surface viscosity emerges from the collective dynamics of interfacial layers in molecular simulations of simple liquids. We estimate the surface viscosity to be 8-16 times smaller than that of the bulk fluid at the thermodynamic point considered. This result can have important implications for reactions at liquid surfaces in atmospheric chemistry and catalysis.
Collapse
Affiliation(s)
- Paolo Malgaretti
- Forschungszentrum Jülich GmbH, Helmholtz Institute Erlangen-Nürnberg for Renewable Energy (IEK-11), Cauerstr.1, D-91058 Erlangen, Germany
| | - Ubaldo Bafile
- Consiglio Nazionale delle Ricerche, Istituto di Fisica Applicata "Nello Carrara," I-50019 Sesto Fiorentino, Italy
| | - Renzo Vallauri
- Consiglio Nazionale delle Ricerche, Istituto di Fisica Applicata "Nello Carrara," I-50019 Sesto Fiorentino, Italy
| | - Pál Jedlovszky
- Department of Chemistry, Eszterházy Károly University, Leányka u. 6, H-3300 Eger, Hungary
| | - Marcello Sega
- Department of Chemical Engineering, University College London, London WC1E 7JE, United Kingdom
| |
Collapse
|
2
|
Maffioli L, Smith ER, Ewen JP, Daivis PJ, Dini D, Todd BD. Slip and stress from low shear rate nonequilibrium molecular dynamics: The transient-time correlation function technique. J Chem Phys 2022; 156:184111. [PMID: 35568555 DOI: 10.1063/5.0088127] [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/15/2022] Open
Abstract
We derive the transient-time correlation function (TTCF) expression for the computation of phase variables of inhomogenous confined atomistic fluids undergoing boundary-driven planar shear (Couette) flow at constant pressure. Using nonequilibrium molecular dynamics simulations, we then apply the TTCF formalism to the computation of the shear stress and the slip velocity for atomistic fluids at realistic low shear rates, in systems under constant pressure and constant volume. We show that, compared to direct averaging of multiple trajectories, the TTCF method dramatically improves the accuracy of the results at low shear rates and that it is suitable to investigate the tribology and rheology of atomistically detailed confined fluids at realistic flow rates.
Collapse
Affiliation(s)
- Luca Maffioli
- Department of Mathematics, School of Science, Computing and Engineering Technologies, Swinburne University of Technology, P.O. Box 218, Hawthorn, Victoria 3122, Australia
| | - Edward R Smith
- Mechanical and Aerospace Engineering, Brunel University London, Kingston Lane, Uxbridge UB8 3PH, United Kingdom
| | - James P Ewen
- Department of Mechanical Engineering, Imperial College London, South Kensington Campus, Exhibition Road, London SW7 2AZ, United Kingdom
| | - Peter J Daivis
- School of Science, RMIT University, GPO Box 2476, Victoria 3001, Australia
| | - Daniele Dini
- Department of Mechanical Engineering, Imperial College London, South Kensington Campus, Exhibition Road, London SW7 2AZ, United Kingdom
| | - B D Todd
- Department of Mathematics, School of Science, Computing and Engineering Technologies, Swinburne University of Technology, P.O. Box 218, Hawthorn, Victoria 3122, Australia
| |
Collapse
|
3
|
Affiliation(s)
- J. S. Hansen
- “Glass and Time”, IMFUFA, Department of Science and Environment, Roskilde University, Roskilde, Denmark
| |
Collapse
|
4
|
A Review of Recent Developments in Molecular Dynamics Simulations of the Photoelectrochemical Water Splitting Process. Catalysts 2021. [DOI: 10.3390/catal11070807] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023] Open
Abstract
In this review, we provide a short overview of the Molecular Dynamics (MD) method and how it can be used to model the water splitting process in photoelectrochemical hydrogen production. We cover classical non-reactive and reactive MD techniques as well as multiscale extensions combining classical MD with quantum chemical and continuum methods. Selected examples of MD investigations of various aqueous semiconductor interfaces with a special focus on TiO2 are discussed. Finally, we identify gaps in the current state-of-the-art where further developments will be needed for better utilization of MD techniques in the field of water splitting.
Collapse
|
5
|
Sgouros AP, Tsagkalakis DS, Theodorou DN. Effect of Surface Nanopatterning on Slip: The Case of Couette Flow of Long-Chain Polyethylene Melt Flowing Past Gold Surfaces. J Phys Chem B 2021; 125:6681-6696. [PMID: 34126736 DOI: 10.1021/acs.jpcb.1c02546] [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/30/2022]
Abstract
The manifestation of slip during flow of a polymer melt past a solid surface depends on several parameters, such as film thickness, the strength of polymer-solid interactions compared to the cohesive energy of the polymer, and the roughness of the surface. Understanding the role of these molecular aspects for slip is crucial in microfluidics, friction-tuning, polymer extrusion, and nanocomposites applications. The present article investigates the effect of surface nanopatterning on slip, via Couette-flow simulations of long chain polyethylene melts past nanopatterned gold surfaces. Slip is quantified in terms of the true and effective slip velocity, and the slip length. When polymer chains are adsorbed to surfaces with periodic features (e.g., crystal planes), they develop preferential ordering in a way that enables them to minimize their free energy. The orientation of a chain is affected by that of its neighbors; thus, when several chains come together, they are prone to form regions with crystal-like orientation. We show that, in some cases, the introduction of nanopatterns on the surface can perturb and induce reorganization of these regions, and in turn affect slip. The nanopatterns are realized as periodic defect stripes of variable width, depth, areal density, and orientation angle. In situations in which the width of the defects becomes comparable to the diameter of individual chain backbones, slip is minimized (stick conditions). Cutting the nanopatterns in low symmetry directions can affect the quality of their edges and lead to enhanced friction. To characterize these edges we have devised a scheme for the quantification of the mean square roughness and mean position of the surface, which is general and applicable in 2 and 3 dimensions for any kind of material, either crystalline of amorphous. Applying the patterns on the opposing solid surfaces in a symmetric or antisymmetric manner has a profound effect on flow. We show that the application of nanopatterns in symmetric configurations generates zero net flow and induces additional shear along directions normal to the direction of the flow. The application of symmetry-breaking configurations can guide flow toward preferential directions, a result with possible applications in microfluidic devices.
Collapse
Affiliation(s)
- A P Sgouros
- School of Chemical Engineering, National Technical University of Athens (NTUA), Athens, GR-15780, Greece
| | - D S Tsagkalakis
- School of Chemical Engineering, National Technical University of Athens (NTUA), Athens, GR-15780, Greece
| | - D N Theodorou
- School of Chemical Engineering, National Technical University of Athens (NTUA), Athens, GR-15780, Greece
| |
Collapse
|
6
|
Alosious S, Kannam SK, Sathian SP, Todd BD. Nanoconfinement Effects on the Kapitza Resistance at Water-CNT Interfaces. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2021; 37:2355-2361. [PMID: 33570421 DOI: 10.1021/acs.langmuir.0c03298] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
The Kapitza resistance (Rk) at the water-carbon nanotube (CNT) interface, with water on the inside of the nanotube, was investigated using molecular dynamics simulations. We propose a new equilibrium molecular dynamics (EMD) method, also valid in the weak flow regime, to determine the Kapitza resistance in a cylindrical nanoconfinement system where nonequilibrium molecular dynamics (NEMD) methods are not suitable. The proposed method is independent of the correlation time compared to Green-Kubo-based methods, which only work in short correlation time intervals. Rk between the CNT and the confined water strongly depends on the diameter of the nanotube and is found to decrease with an increase in the CNT diameter, the opposite to what is reported in the literature when water is on the outside of the nanotube. Rk is furthermore found to converge to the planar graphene surface value as the number of water molecules per unit surface area approaches the value in the graphene surface and a higher overlap of the vibrational spectrum. A slight increase in Rk with the addition of the number of CNT walls was observed, whereas the chirality and flow do not have any impact.
Collapse
Affiliation(s)
- Sobin Alosious
- Department of Applied Mechanics, Indian Institute of Technology Madras, Chennai 600036, India
- Department of Mathematics, School of Science, Swinburne University of Technology, Melbourne, Victoria 3122, Australia
| | - Sridhar Kumar Kannam
- Department of Mathematics, School of Science, Swinburne University of Technology, Melbourne, Victoria 3122, Australia
| | - Sarith P Sathian
- Department of Applied Mechanics, Indian Institute of Technology Madras, Chennai 600036, India
| | - B D Todd
- Department of Mathematics, School of Science, Swinburne University of Technology, Melbourne, Victoria 3122, Australia
| |
Collapse
|
7
|
Herrero C, Tocci G, Merabia S, Joly L. Fast increase of nanofluidic slip in supercooled water: the key role of dynamics. NANOSCALE 2020; 12:20396-20403. [PMID: 33021296 DOI: 10.1039/d0nr06399a] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
Nanofluidics is an emerging field offering innovative solutions for energy harvesting and desalination. The efficiency of these applications depends strongly on liquid-solid slip, arising from a favorable ratio between viscosity and interfacial friction. Using molecular dynamics simulations, we show that wall slip increases strongly when water is cooled below its melting point. For water on graphene, the slip length is multiplied by up to a factor of five and reaches 230 nm at the lowest simulated temperature, T ∼ 225 K; experiments in nanopores can reach much lower temperatures and could reveal even more drastic changes. The predicted fast increase in water slip can also be detected at supercoolings reached experimentally in bulk water, as well as in droplets flowing on anti-icing surfaces. We explain the anomalous slip behavior in the supercooled regime by a decoupling between viscosity and bulk density relaxation dynamics, and we rationalize the wall-type dependence of the enhancement in terms of interfacial density relaxation dynamics. While providing fundamental insights on the molecular mechanisms of hydrodynamic transport in both interfacial and bulk water in the supercooled regime, this study is relevant to the design of anti-icing surfaces, could help explain the subtle phase and dynamical behaviors of supercooled confined water, and paves the way to explore new behaviors in supercooled nanofluidic systems.
Collapse
Affiliation(s)
- Cecilia Herrero
- Univ Lyon, Univ Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622, Villeurbanne, France.
| | - Gabriele Tocci
- Department of Chemistry, Universität Zürich, 8057 Zürich, Switzerland
| | - Samy Merabia
- 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. and Institut Universitaire de France (IUF), 1 rue Descartes, 75005 Paris, France
| |
Collapse
|
8
|
Sam A, Hartkamp R, Kumar Kannam S, Babu JS, Sathian SP, Daivis PJ, Todd BD. Fast transport of water in carbon nanotubes: a review of current accomplishments and challenges. MOLECULAR SIMULATION 2020. [DOI: 10.1080/08927022.2020.1782401] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Affiliation(s)
- Alan Sam
- Department of Applied Mechanics, Indian Institute of Technology Madras, Chennai, India
| | - Remco Hartkamp
- Process and Energy Department, Delft University of Technology, Delft, The Netherlands
| | - Sridhar Kumar Kannam
- Faculty of Science, Engineering and Technology, Swinburne University of Technology, Melbourne, Australia
| | - Jeetu S. Babu
- Department of Mechanical Engineering, Amrita Vishwa Vidyapeetham, Amritapuri, India
| | - Sarith P. Sathian
- Department of Applied Mechanics, Indian Institute of Technology Madras, Chennai, India
| | - Peter J. Daivis
- School of Applied Sciences, RMIT University, Melbourne, Australia
| | - B. D. Todd
- Department of Mathematics, Swinburne University of Technology, Melbourne, Australia
| |
Collapse
|
9
|
Alosious S, Kannam SK, Sathian SP, Todd BD. Kapitza resistance at water-graphene interfaces. J Chem Phys 2020; 152:224703. [PMID: 32534537 DOI: 10.1063/5.0009001] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Heat transfer across fluid-solid interfaces in nanoconfinement has received significant attention due to its relevance in nanoscale systems. In this study, we investigate the Kapitza resistance at the water-graphene interface with the help of classical molecular dynamics simulation techniques in conjunction with our recently proposed equilibrium molecular dynamics (EMD) method [S. Alosious et al., J. Chem. Phys. 151, 194502 (2019)]. The size effect of the Kapitza resistance on different factors such as the number of graphene layers, the cross-sectional area, and the width of the water block was studied. The Kapitza resistance decreases slightly with an increase in the number of layers, while the influence of the cross-sectional area and the width of the water block is negligible. The variation in the Kapitza resistance as a function of the number of graphene layers is attributed to the large phonon mean free path along the graphene cross-plane. An optimum water-graphene system, which is independent of size effects, was selected, and the same was used to determine the Kapitza resistance using the predicted EMD method. The values obtained from both the EMD and the non-equilibrium molecular dynamics (NEMD) methods were compared for different potentials and water models, and the results are shown to be in good agreement. Our method allows us to compute the Kapitza resistance using EMD simulations, which obviates the need to create a large temperature gradient required for the NEMD method.
Collapse
Affiliation(s)
- Sobin Alosious
- Department of Applied Mechanics, Indian Institute of Technology Madras, Chennai 600036, India
| | - Sridhar Kumar Kannam
- Department of Mathematics, Faculty of Science, Engineering and Technology, Swinburne University of Technology, Melbourne, Victoria 3122, Australia
| | - Sarith P Sathian
- Department of Applied Mechanics, Indian Institute of Technology Madras, Chennai 600036, India
| | - B D Todd
- Department of Mathematics, Faculty of Science, Engineering and Technology, Swinburne University of Technology, Melbourne, Victoria 3122, Australia
| |
Collapse
|
10
|
Xiong H, Devegowda D, Huang L. Oil–water
transport in
clay‐hosted
nanopores: Effects of
long‐range
electrostatic forces. AIChE J 2020. [DOI: 10.1002/aic.16276] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Affiliation(s)
- Hao Xiong
- Mewbourne School of Petroleum and Geological EngineeringThe University of Oklahoma Norman Oklahoma USA
| | - Deepak Devegowda
- Mewbourne School of Petroleum and Geological EngineeringThe University of Oklahoma Norman Oklahoma USA
| | - Liangliang Huang
- Chemical, Biological & Materials EngineeringThe University of Oklahoma Norman Oklahoma USA
| |
Collapse
|
11
|
Deng X, Wei X, Wang X, Sheng P. Decomposing thermal fluctuations with hydrodynamic modes. Phys Rev E 2020; 101:063104. [PMID: 32688559 DOI: 10.1103/physreve.101.063104] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2019] [Accepted: 05/29/2020] [Indexed: 06/11/2023]
Abstract
We have obtained analytically the complete set of hydrodynamic modes (HMs) for a two-dimensional (2D) fluid confined within a channel with the Navier slip boundary condition at the hydrodynamic boundary. The HMs are orthogonal to each other and hence each represents an independent degree of freedom. We show that the HMs can be used to recursively generate a time series of random thermal fluctuations of displacement velocity, with identical statistical distributions as those obtained from MD simulations. By projecting the HMs onto molecular dynamics (MD) configurations and evaluating the resulting decay time from the autocorrelation function, we obtain from MD the eigenvalues of the HMs. Multiplying two different HMs and integrating as a function of z from center of the channel towards the fluid-solid interface, the position of the hydrodynamic boundary (HDB) is unambiguously identified as the point at which the integral vanishes. Invariably the HDB is located inside the fluid domain and not on the liquid-solid interface. With the knowledge of the HDB position, the value of the slip length can be obtained directly from HM's dispersion relation. We show that in terms of the complete set of HMs, the fluctuation-dissipation theorem may be expressed in a simple expression involving the average of the inverse of the eigenvalues. Besides offering an alternative perspective on thermal fluctuations and hydrodynamic boundary, the present work opens the possibility of using modulated boundary conditions to manipulate thermal fluctuations in mesoscopic channels, which can lead to interesting statistical mechanical consequences.
Collapse
Affiliation(s)
- Xiaohui Deng
- Department of Physics, Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong, China
| | - Xiaoyu Wei
- Department of Mathematics, Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong, China
| | - Xiaoping Wang
- Department of Mathematics, Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong, China
| | - Ping Sheng
- Department of Physics, Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong, China
| |
Collapse
|
12
|
|
13
|
Nakano H, Sasa SI. Equilibrium measurement method of slip length based on fluctuating hydrodynamics. Phys Rev E 2020; 101:033109. [PMID: 32289937 DOI: 10.1103/physreve.101.033109] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2019] [Accepted: 03/02/2020] [Indexed: 11/07/2022]
Abstract
We perform equilibrium molecular dynamics simulations for nanoscale fluids confined between two parallel walls and investigate how the autocorrelation function of force acting on one wall is related to the slip length. We demonstrate that for atomically smooth surfaces, the autocorrelation function is accurately described by linearized fluctuating hydrodynamics (LFH). Excellent agreement between the simulation and the LFH solution is found over a wide range of scales, specifically, from the timescale of fluid relaxation even to that of molecular motion. Fitting the simulation data yields a reasonable estimation of the slip length. We show that LFH provides a starting point for examining the relationship between the slip length and the force fluctuations.
Collapse
Affiliation(s)
| | - Shin-Ichi Sasa
- Department of Physics, Kyoto University, Kyoto 606-8502, Japan
| |
Collapse
|
14
|
de la Torre JA, Duque-Zumajo D, Camargo D, Español P. Microscopic Slip Boundary Conditions in Unsteady Fluid Flows. PHYSICAL REVIEW LETTERS 2019; 123:264501. [PMID: 31951457 DOI: 10.1103/physrevlett.123.264501] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/14/2019] [Indexed: 06/10/2023]
Abstract
An algebraic tail in the Green-Kubo integral for the solid-fluid friction coefficient hampers its use in the determination of the slip length. A simple theory for discrete nonlocal hydrodynamics near parallel solid walls with extended friction forces is given. We explain the origin of the algebraic tail and give a solution of the plateau problem in the Green-Kubo expressions. We derive the slip boundary condition with a microscopic expression for the slip length and the hydrodynamic wall position, and assess it through simulations of an unsteady plug flow.
Collapse
Affiliation(s)
- J A de la Torre
- Dept. Física Fundamental, Universidad Nacional de Educación a Distancia, Aptdo. 60141 E-28080, Madrid, Spain
| | - D Duque-Zumajo
- Dept. Física Fundamental, Universidad Nacional de Educación a Distancia, Aptdo. 60141 E-28080, Madrid, Spain
| | - D Camargo
- Facultad de Ingeniería y Arquitectura, Universidad Pontificia Bolivariana, CO-230003 Montería, Colombia
| | - Pep Español
- Dept. Física Fundamental, Universidad Nacional de Educación a Distancia, Aptdo. 60141 E-28080, Madrid, Spain
| |
Collapse
|
15
|
Duque-Zumajo D, de la Torre JA, Camargo D, Español P. Discrete hydrodynamics near solid walls: Non-Markovian effects and the slip boundary condition. Phys Rev E 2019; 100:062133. [PMID: 31962479 DOI: 10.1103/physreve.100.062133] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/14/2019] [Indexed: 06/10/2023]
Abstract
A simple Markovian theory for the prediction of averages and correlations of discrete hydrodynamics near parallel solid walls is presented. The discrete momentum of bins is defined through a finite element basis function. The effect of the walls on the fluid is through irreversible extended friction forces appearing in the very equations of hydrodynamics. The Markovian assumption is critically assessed from the exponential decay of the eigenvalues of the correlation matrix of the discrete transverse momentum. We observe that for bins smaller than molecular dimensions, allowing one to resolve the density layering near the wall, the dynamics near the wall is non-Markovian. Bins larger than the molecular size do behave in a Markovian way. We measure the nonlocal viscosity and frictions kernels that appear in the discrete hydrodynamic equations, which are given in terms of Green-Kubo formulas. They suffer dramatically from the plateau problem. We use a recent procedure for reliably extracting the transport kernels out of the plateau-problematic Green-Kubo formula. With the so-measured transport kernels the nonlocal theory predicts very well the decay of the average of the transverse momentum when the initial velocity profile is a plug flow. The theory allows us to derive the slip boundary condition with microscopic expressions for the slip length and the hydrodynamic position of the wall. The slip boundary condition is not satisfied at the initial stages of the discontinous plug flow, but good agreement is obtained at later stages.
Collapse
Affiliation(s)
- D Duque-Zumajo
- Dept. Física Fundamental, Universidad Nacional de Educación a Distancia, Aptdo. 60141 E-28080, Madrid, Spain
| | - J A de la Torre
- Dept. Física Fundamental, Universidad Nacional de Educación a Distancia, Aptdo. 60141 E-28080, Madrid, Spain
| | - Diego Camargo
- Facultad de Ingeniería y Arquitectura, Universidad Pontificia Bolivariana, Montería 230002, Colombia
| | - Pep Español
- Dept. Física Fundamental, Universidad Nacional de Educación a Distancia, Aptdo. 60141 E-28080, Madrid, Spain
| |
Collapse
|
16
|
Alosious S, Kannam SK, Sathian SP, Todd BD. Prediction of Kapitza resistance at fluid-solid interfaces. J Chem Phys 2019; 151:194502. [PMID: 31757152 DOI: 10.1063/1.5126887] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Understanding the interfacial heat transfer and thermal resistance at an interface between two dissimilar materials is of great importance in the development of nanoscale systems. This paper introduces a new and reliable linear response method for calculating the interfacial thermal resistance or Kapitza resistance in fluid-solid interfaces with the use of equilibrium molecular dynamics (EMD) simulations. The theoretical predictions are validated against classical molecular dynamics (MD) simulations. MD simulations are carried out in a Lennard-Jones (L-J) system with fluid confined between two solid slabs. Different types of interfaces are tested by varying the fluid-solid interactions (wetting coefficient) at the interface. It is observed that the Kapitza length decreases monotonically with an increasing wetting coefficient as expected. The theory is further validated by simulating under different conditions such as channel width, density, and temperature. Our method allows us to directly determine the Kapitza length from EMD simulations by considering the temperature fluctuation and heat flux fluctuations at the interface. The predicted Kapitza length shows an excellent agreement with the results obtained from both EMD and non-equilibrium MD simulations.
Collapse
Affiliation(s)
- Sobin Alosious
- Department of Applied Mechanics, Indian Institute of Technology Madras, Chennai 600036, India
| | - Sridhar Kumar Kannam
- Department of Mathematics, Faculty of Science, Engineering and Technology, Swinburne University of Technology, Melbourne, Victoria 3122, Australia
| | - Sarith P Sathian
- Department of Applied Mechanics, Indian Institute of Technology Madras, Chennai 600036, India
| | - B D Todd
- Department of Mathematics, Faculty of Science, Engineering and Technology, Swinburne University of Technology, Melbourne, Victoria 3122, Australia
| |
Collapse
|
17
|
Celebi AT, Nguyen CT, Hartkamp R, Beskok A. The role of water models on the prediction of slip length of water in graphene nanochannels. J Chem Phys 2019; 151:174705. [DOI: 10.1063/1.5123713] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Alper Tunga Celebi
- Process and Energy Department, Delft University of Technology, Leeghwaterstraat 39, 2628 CB Delft, The Netherlands
| | - Chinh Thanh Nguyen
- Lyle School of Engineering, Southern Methodist University, 3101 Dyer Street, Dallas, Texas 75205, USA
| | - Remco Hartkamp
- Process and Energy Department, Delft University of Technology, Leeghwaterstraat 39, 2628 CB Delft, The Netherlands
| | - Ali Beskok
- Lyle School of Engineering, Southern Methodist University, 3101 Dyer Street, Dallas, Texas 75205, USA
| |
Collapse
|
18
|
Ogawa K, Oga H, Kusudo H, Yamaguchi Y, Omori T, Merabia S, Joly L. Large effect of lateral box size in molecular dynamics simulations of liquid-solid friction. Phys Rev E 2019; 100:023101. [PMID: 31574745 DOI: 10.1103/physreve.100.023101] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2019] [Indexed: 11/07/2022]
Abstract
Molecular dynamics simulations are a powerful tool to characterize liquid-solid friction. A slab configuration with periodic boundary conditions in the lateral dimensions is commonly used, where the measured friction coefficient could be affected by the finite lateral size of the simulation box. Here we show that for a very wetting liquid close to its melting temperature, strong finite size effects can persist up to large box sizes along the flow direction, typically ∼30 particle diameters. We relate the observed decrease of friction in small boxes to changes in the structure of the first adsorbed layer, which becomes less commensurable with the wall structure. Although these effects disappear for lower wetting cases or at higher temperatures, we suggest that the possible effect of the finite lateral box size on the friction coefficient should not be automatically set aside when exploring unknown systems.
Collapse
Affiliation(s)
- Koshun Ogawa
- Department of Mechanical Engineering, Osaka University, 2-1 Yamadaoka, Suita 565-0871, Japan
| | - Haruki Oga
- Department of Mechanical Engineering, Osaka University, 2-1 Yamadaoka, Suita 565-0871, Japan
| | - Hiroki Kusudo
- Department of Mechanical Engineering, Osaka University, 2-1 Yamadaoka, Suita 565-0871, Japan
| | - Yasutaka Yamaguchi
- Department of Mechanical Engineering, Osaka University, 2-1 Yamadaoka, Suita 565-0871, Japan.,Water Frontier Science & Technology Research Center, Research Institute for Science & Technology, Tokyo University of Science, 1-3 Kagurazaka, Shinjuku-ku, Tokyo 162-8601, Japan
| | - Takeshi Omori
- Department of Mechanical Engineering, Osaka University, 2-1 Yamadaoka, Suita 565-0871, Japan
| | - Samy Merabia
- 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
| |
Collapse
|
19
|
Oga H, Yamaguchi Y, Omori T, Merabia S, Joly L. Green-Kubo measurement of liquid-solid friction in finite-size systems. J Chem Phys 2019. [DOI: 10.1063/1.5104335] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Haruki Oga
- Department of Mechanical Engineering, Osaka University, 2-1 Yamadaoka, Suita 565-0871, Japan
| | - Yasutaka Yamaguchi
- Department of Mechanical Engineering, Osaka University, 2-1 Yamadaoka, Suita 565-0871, Japan
- Water Frontier Science and Technology Research Center (W-FST), Research Institute for Science and Technology, Tokyo University of Science, 1-3 Kagurazaka, Shinjuku-ku, Tokyo, 162-8601, Japan
| | - Takeshi Omori
- Department of Mechanical Engineering, Osaka University, 2-1 Yamadaoka, Suita 565-0871, Japan
| | - Samy Merabia
- 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
| |
Collapse
|
20
|
Herrero C, Omori T, Yamaguchi Y, Joly L. Shear force measurement of the hydrodynamic wall position in molecular dynamics. J Chem Phys 2019; 151:041103. [DOI: 10.1063/1.5111966] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Cecilia Herrero
- Univ Lyon, Univ Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622 Villeurbanne, France
| | - Takeshi Omori
- Department of Mechanical Engineering, Osaka University, 2-1 Yamadaoka, Suita 565-0871, Japan
| | - Yasutaka Yamaguchi
- Department of Mechanical Engineering, Osaka University, 2-1 Yamadaoka, Suita 565-0871, Japan
- Water Frontier Science and Technology Research Center (W-FST), Research Institute for Science and Technology, Tokyo University of Science, 1-3 Kagurazaka, Shinjuku-ku, Tokyo 162-8601, Japan
| | - Laurent Joly
- Univ Lyon, Univ Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622 Villeurbanne, France
| |
Collapse
|
21
|
Duque-Zumajo D, Camargo D, de la Torre JA, Chejne F, Español P. Discrete hydrodynamics near solid planar walls. Phys Rev E 2019; 99:052130. [PMID: 31212438 DOI: 10.1103/physreve.99.052130] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/25/2019] [Indexed: 06/09/2023]
Abstract
We derive, with the projection operator technique, the equations of motion for the time-dependent average of the discrete mass and momentum densities of a fluid confined by planar walls under the assumption that the flow field is translationally invariant along the directions tangent to the walls. Shear flow and sound propagation perpendicular to the walls can be described with the discrete hydrodynamic equations. The interaction with the walls is not given through boundary conditions but rather in terms of impenetrability and friction forces appearing in the discrete hydrodynamic equations. Microscopic expressions for the transport coefficients entering the discrete equations are provided. We further show that the obtained discrete equations can be interpreted as a Petrov-Galerkin finite-element discretization of the continuum equations presented by Camargo et al. [J. Chem. Phys. 148, 064107 (2018)JCPSA60021-960610.1063/1.5010401] when restricted to planar geometries and flows.
Collapse
Affiliation(s)
- D Duque-Zumajo
- Departamento Física Fundamental, Universidad Nacional de Educación a Distancia, Apartado 60141, 28080 Madrid, Spain
| | - Diego Camargo
- Escuela de Ingeniería y Arquitectura, Universidad Pontificia Bolivariana, Montería, Colombia
| | - J A de la Torre
- Departamento Física Fundamental, Universidad Nacional de Educación a Distancia, Apartado 60141, 28080 Madrid, Spain
| | - Farid Chejne
- Facultad de Minas, Universidad Nacional de Colombia, Medellín, Colombia
| | - Pep Español
- Departamento Física Fundamental, Universidad Nacional de Educación a Distancia, Apartado 60141, 28080 Madrid, Spain
| |
Collapse
|
22
|
Camargo D, de la Torre JA, Delgado-Buscalioni R, Chejne F, Español P. Boundary conditions derived from a microscopic theory of hydrodynamics near solids. J Chem Phys 2019; 150:144104. [DOI: 10.1063/1.5088354] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023] Open
Affiliation(s)
- Diego Camargo
- Facultad de Ingeniería y Arquitectura, Universidad Pontificía Bolivariana, Montería, Colombia
| | - J. A. de la Torre
- Dept. Física Fundamental, Universidad Nacional de Educación a Distancia, Aptdo. 60141, E-28080 Madrid, Spain
| | - Rafael Delgado-Buscalioni
- Departamento de Física Teórica de la Materia Condensada Universidad Autónoma de Madrid, and Condensed Matter Physics Center (IFIMAC), Madrid 28049, Spain
| | - Farid Chejne
- Universidad Nacional de Colombia, Bogotá, Colombia
| | - Pep Español
- Dept. Física Fundamental, Universidad Nacional de Educación a Distancia, Aptdo. 60141, E-28080 Madrid, Spain
| |
Collapse
|
23
|
Sam A, Hartkamp R, Kannam SK, Sathian SP. Prediction of fluid slip in cylindrical nanopores using equilibrium molecular simulations. NANOTECHNOLOGY 2018; 29:485404. [PMID: 30207542 DOI: 10.1088/1361-6528/aae0bd] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
We introduce an analytical method to predict the slip length (L s) in cylindrical nanopores using equilibrium molecular dynamics (EMD) simulations, following the approach proposed by Sokhan and Quirke for planar channels [39]. Using this approach, we determined the slip length of water in carbon nanotubes (CNTs) of various diameters. The slip length predicted from our method shows excellent agreement with the results obtained from nonequilibrium molecular dynamics (NEMD) simulations. The data show a monotonically decreasing slip length with an increasing nanotube diameter. The proposed EMD method can be used to precisely estimate slip length in high slip cylindrical systems, whereas, L s calculated from NEMD is highly sensitive to the velocity profile and may cause large statistical errors due to large velocity slip at the channel surface. We also demonstrated the validity of the EMD method in a BNNT-water system, where the slip length is very small compared to that in a CNT pore of similar diameter. The developed method enables us to calculate the interfacial friction coefficient directly from EMD simulations, while friction can be estimated using NEMD by performing simulations at various external driving forces, thereby increasing the overall computational time. The EMD analysis revealed a curvature dependence in the friction coefficient, which induces the slip length dependency on the tube diameter. Conversely, in flat graphene nanopores, both L s and friction coefficient show no strong dependency on the channel width.
Collapse
Affiliation(s)
- Alan Sam
- Department of Applied Mechanics, Indian Institute of Technology Madras, Chennai, India
| | | | | | | |
Collapse
|
24
|
Smith ER, Theodorakis PE, Craster RV, Matar OK. Moving Contact Lines: Linking Molecular Dynamics and Continuum-Scale Modeling. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2018; 34:12501-12518. [PMID: 29727189 DOI: 10.1021/acs.langmuir.8b00466] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
Despite decades of research, the modeling of moving contact lines has remained a formidable challenge in fluid dynamics whose resolution will impact numerous industrial, biological, and daily life applications. On the one hand, molecular dynamics (MD) simulation has the ability to provide unique insight into the microscopic details that determine the dynamic behavior of the contact line, which is not possible with either continuum-scale simulations or experiments. On the other hand, continuum-based models provide a link to the macroscopic description of the system. In this Feature Article, we explore the complex range of physical factors, including the presence of surfactants, which governs the contact line motion through MD simulations. We also discuss links between continuum- and molecular-scale modeling and highlight the opportunities for future developments in this area.
Collapse
|
25
|
Abstract
The fluid dynamics of macroscopic and microscopic systems is well developed and has been extensively validated. Its extraordinary success makes it tempting to apply Navier–Stokes fluid dynamics without modification to systems of ever decreasing dimensions as studies of nanofluidics become more prevalent. However, this can result in serious error. In this paper, we discuss several ways in which nanoconfined fluid flow differs from macroscopic flow. We give particular attention to several topics that have recently received attention in the literature: slip, spin angular momentum coupling, nonlocal stress response and density inhomogeneity. In principle, all of these effects can now be accurately modelled using validated theories. Although the basic principles are now fairly well understood, much work remains to be done in their application.
Collapse
|
26
|
Suga K, Mori Y, Moritani R, Kaneda M. Combined effects of molecular geometry and nanoconfinement on liquid flows through carbon nanotubes. Phys Rev E 2018; 97:053109. [PMID: 29906844 DOI: 10.1103/physreve.97.053109] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2018] [Indexed: 06/08/2023]
Abstract
Molecular dynamics simulations are carried out to investigate the geometry effects of diatomic molecules on liquid flows in carbon nanotubes (CNTs). Oxygen molecules are considered as the fluid inside armchair (n,n) (n=6-20) CNTs. The simulated fluid temperature and bulk pressure for the liquid state are T=133 K and ρ_{b}=1346kg/m^{3}, respectively. In the agglomerated molecular cluster, nanoconfinement-induced structural changes are observed. As the CNT diameter decreases, it is confirmed that the flow rate significantly increases with irregular trends (discontinuity points in the profiles). From the discussion of the structure of the agglomerated fluid molecules, it is found that those trends are not simply caused by the structural changes. The main factor to induce the irregularity is confirmed to be the interlayer molecular movement affected by the combination of the molecular geometry and the arrangement of the multilayered structure.
Collapse
Affiliation(s)
- Kazuhiko Suga
- Department of Mechanical Engineering, Osaka Prefecture University, 1-1 Gakuen-cho, Naka-ku, Sakai, Osaka 599-8531, Japan
| | - Yuki Mori
- Department of Mechanical Engineering, Osaka Prefecture University, 1-1 Gakuen-cho, Naka-ku, Sakai, Osaka 599-8531, Japan
| | - Rintaro Moritani
- Department of Mechanical Engineering, Osaka Prefecture University, 1-1 Gakuen-cho, Naka-ku, Sakai, Osaka 599-8531, Japan
| | - Masayuki Kaneda
- Department of Mechanical Engineering, Osaka Prefecture University, 1-1 Gakuen-cho, Naka-ku, Sakai, Osaka 599-8531, Japan
| |
Collapse
|
27
|
Camargo D, de la Torre JA, Duque-Zumajo D, Español P, Delgado-Buscalioni R, Chejne F. Nanoscale hydrodynamics near solids. J Chem Phys 2018; 148:064107. [DOI: 10.1063/1.5010401] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023] Open
Affiliation(s)
- Diego Camargo
- Facultad de Minas, Universidad Nacional de Colombia, Medellin, Colombia
- Facultad Mecánica, Universidad Pontificia Bolivariana, Montería, Colombia
| | - J. A. de la Torre
- Departamento Física Fundamental, Universidad Nacional de Educación a Distancia, Aptdo. 60141, E-28080 Madrid, Spain
| | - D. Duque-Zumajo
- Departamento Física Fundamental, Universidad Nacional de Educación a Distancia, Aptdo. 60141, E-28080 Madrid, Spain
| | - Pep Español
- Departamento Física Fundamental, Universidad Nacional de Educación a Distancia, Aptdo. 60141, E-28080 Madrid, Spain
| | - Rafael Delgado-Buscalioni
- Departamento de Física Teórica de la Materia Condensada, Universidad Autónoma de Madrid, Madrid 28049, Spain
| | - Farid Chejne
- Facultad de Minas, Universidad Nacional de Colombia, Medellin, Colombia
| |
Collapse
|
28
|
Ang EYM, Ng TY, Yeo J, Lin R, Liu Z, Geethalakshmi KR. Effects of CNT size on the desalination performance of an outer-wall CNT slit membrane. Phys Chem Chem Phys 2018; 20:13896-13902. [DOI: 10.1039/c8cp01191e] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We investigate the effect of varying carbon nanotube (CNT) size on the desalination performance through slit confinements formed by horizontally aligned CNTs stacked on top of one another.
Collapse
Affiliation(s)
- Elisa Y. M. Ang
- School of Mechanical and Aerospace Engineering
- Nanyang Technological University
- Singapore 639798
- Singapore
| | - Teng Yong Ng
- School of Mechanical and Aerospace Engineering
- Nanyang Technological University
- Singapore 639798
- Singapore
| | - Jingjie Yeo
- Department of Civil and Environmental Engineering
- Massachusetts Institute of Technology
- Cambridge
- Massachusetts 02139
- USA
| | - Rongming Lin
- School of Mechanical and Aerospace Engineering
- Nanyang Technological University
- Singapore 639798
- Singapore
| | - Zishun Liu
- International Center for Applied Mechanics
- State Key Laboratory for Strength and Vibration of Mechanical Structures
- Xi’an Jiaotong University
- Xi'an 710049
- P. R. China
| | - K. R. Geethalakshmi
- School of Mechanical and Aerospace Engineering
- Nanyang Technological University
- Singapore 639798
- Singapore
| |
Collapse
|
29
|
Bhadauria R, Aluru NR. A multiscale transport model for non-classical nanochannel electroosmosis. J Chem Phys 2017; 147:214105. [DOI: 10.1063/1.5005127] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [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
|
Bhadauria R, Aluru NR. Multiscale modeling of electroosmotic flow: Effects of discrete ion, enhanced viscosity, and surface friction. J Chem Phys 2017. [DOI: 10.1063/1.4982731] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/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
|
31
|
Shu JJ, Teo JBM, Chan WK. A new model for fluid velocity slip on a solid surface. SOFT MATTER 2016; 12:8388-8397. [PMID: 27714378 DOI: 10.1039/c6sm01178k] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
A general adsorption model is developed to describe the interactions between near-wall fluid molecules and solid surfaces. This model serves as a framework for the theoretical modelling of boundary slip phenomena. Based on this adsorption model, a new general model for the slip velocity of fluids on solid surfaces is introduced. The slip boundary condition at a fluid-solid interface has hitherto been considered separately for gases and liquids. In this paper, we show that the slip velocity in both gases and liquids may originate from dynamical adsorption processes at the interface. A unified analytical model that is valid for both gas-solid and liquid-solid slip boundary conditions is proposed based on surface science theory. The corroboration with the experimental data extracted from the literature shows that the proposed model provides an improved prediction compared to existing analytical models for gases at higher shear rates and close agreement for liquid-solid interfaces in general.
Collapse
Affiliation(s)
- Jian-Jun Shu
- School of Mechanical & Aerospace Engineering, Nanyang Technological University, 50 Nanyang Avenue, Singapore 639798.
| | - Ji Bin Melvin Teo
- School of Mechanical & Aerospace Engineering, Nanyang Technological University, 50 Nanyang Avenue, Singapore 639798.
| | - Weng Kong Chan
- School of Mechanical & Aerospace Engineering, Nanyang Technological University, 50 Nanyang Avenue, Singapore 639798.
| |
Collapse
|
32
|
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
|
33
|
Neek-Amal M, Peeters FM, Grigorieva IV, Geim AK. Commensurability Effects in Viscosity of Nanoconfined Water. ACS NANO 2016; 10:3685-3692. [PMID: 26882095 DOI: 10.1021/acsnano.6b00187] [Citation(s) in RCA: 100] [Impact Index Per Article: 12.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
The rate of water flow through hydrophobic nanocapillaries is greatly enhanced as compared to that expected from macroscopic hydrodynamics. This phenomenon is usually described in terms of a relatively large slip length, which is in turn defined by such microscopic properties as the friction between water and capillary surfaces and the viscosity of water. We show that the viscosity of water and, therefore, its flow rate are profoundly affected by the layered structure of confined water if the capillary size becomes less than 2 nm. To this end, we study the structure and dynamics of water confined between two parallel graphene layers using equilibrium molecular dynamics simulations. We find that the shear viscosity is not only greatly enhanced for subnanometer capillaries, but also exhibits large oscillations that originate from commensurability between the capillary size and the size of water molecules. Such oscillating behavior of viscosity and, consequently, the slip length should be taken into account in designing and studying graphene-based and similar membranes for desalination and filtration.
Collapse
Affiliation(s)
- Mehdi Neek-Amal
- Department of Physics, Shahid Rajaee Teacher Training University , Lavizan, 16785-136 Tehran, Iran
| | - Francois M Peeters
- Departement Fysica, Universiteit Antwerpen , Groenenborgerlaan 171, B-2020 Antwerpen, Belgium
| | - Irina V Grigorieva
- School of Physics and Astronomy, University of Manchester , Manchester M13 9PL, United Kingdom
| | - Andre K Geim
- School of Physics and Astronomy, University of Manchester , Manchester M13 9PL, United Kingdom
| |
Collapse
|
34
|
Ramos-Alvarado B, Kumar S, Peterson GP. Hydrodynamic slip length as a surface property. Phys Rev E 2016; 93:023101. [PMID: 26986407 DOI: 10.1103/physreve.93.023101] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/09/2015] [Indexed: 06/05/2023]
Abstract
Equilibrium and nonequilibrium molecular dynamics simulations were conducted in order to evaluate the hypothesis that the hydrodynamic slip length is a surface property. The system under investigation was water confined between two graphite layers to form nanochannels of different sizes (3-8 nm). The water-carbon interaction potential was calibrated by matching wettability experiments of graphitic-carbon surfaces free of airborne hydrocarbon contamination. Three equilibrium theories were used to calculate the hydrodynamic slip length. It was found that one of the recently reported equilibrium theories for the calculation of the slip length featured confinement effects, while the others resulted in calculations significantly hindered by the large margin of error observed between independent simulations. The hydrodynamic slip length was found to be channel-size independent using equilibrium calculations, i.e., suggesting a consistency with the definition of a surface property, for 5-nm channels and larger. The analysis of the individual trajectories of liquid particles revealed that the reason for observing confinement effects in 3-nm nanochannels is the high mobility of the bulk particles. Nonequilibrium calculations were not consistently affected by size but by noisiness in the smallest systems.
Collapse
Affiliation(s)
- Bladimir Ramos-Alvarado
- The George W. Woodruff School of Mechanical Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| | - Satish Kumar
- The George W. Woodruff School of Mechanical Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| | - G P Peterson
- The George W. Woodruff School of Mechanical Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| |
Collapse
|
35
|
Hansen JS, Dyre JC, Daivis P, Todd BD, Bruus H. Continuum Nanofluidics. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2015; 31:13275-89. [PMID: 26457405 DOI: 10.1021/acs.langmuir.5b02237] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/13/2023]
Abstract
This paper introduces the fundamental continuum theory governing momentum transport in isotropic nanofluidic systems. The theory is an extension of the classical Navier-Stokes equation, and includes coupling between translational and rotational degrees of freedom as well as nonlocal response functions that incorporate spatial correlations. The continuum theory is compared with molecular dynamics simulation data for both relaxation processes and fluid flows, showing excellent agreement on the nanometer length scale. We also present practical tools to estimate when the extended theory should be used. It is shown that in the wall-fluid region the fluid molecules align with the wall, and in this region the isotropic model may fail and a full anisotropic description is necessary.
Collapse
Affiliation(s)
- Jesper S Hansen
- DNRF Centre "Glass and Time", IMFUFA, Department of Sciences, Roskilde University , Postbox 260, DK-4000 Roskilde, Denmark
| | - Jeppe C Dyre
- DNRF Centre "Glass and Time", IMFUFA, Department of Sciences, Roskilde University , Postbox 260, DK-4000 Roskilde, Denmark
| | - Peter Daivis
- Applied Physics, School of Applied Sciences, RMIT University , GPO Box 2476, Melbourne, Victoria 3001, Australia
| | - Billy D Todd
- Department of Mathematics, Faculty of Science, Engineering and Technology, and Center for Molecular Simulation, Swinburne University of Technology , P.O. Box 218, Hawthorn, Victoria 3122, Australia
| | - Henrik Bruus
- Department of Physics, Technical University of Denmark , DTU Physics Building 309, DK-2800 Kongens Lyngby, Denmark
| |
Collapse
|
36
|
Yasuoka H, Takahama R, Kaneda M, Suga K. Confinement effects on liquid-flow characteristics in carbon nanotubes. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2015; 92:063001. [PMID: 26764798 DOI: 10.1103/physreve.92.063001] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/21/2015] [Indexed: 06/05/2023]
Abstract
Liquid flow dynamics through the armchair (6,6)-(160,160) carbon nanotubes (CNTs) is elucidated by molecular dynamics simulations. The liquid is modeled by nonpolar argon atoms to understand the fundamental flow physics. The velocity profiles and slip lengths are discussed considering the radial distributions of the fluid density by the presently proposed finite difference-based velocity fitting method. It is found that as the CNT diameter D increases, the slip length and the flow rate enhancement show three-step transitional profiles in the region of D≤2.3 nm. The slip length and the flow rate stepwise increase at the first transition while they drop at the second and third transitions. The first transition corresponds to the structural change from the single-file chain to single-ring structures of the molecule cluster. The second and third transitions take place when the ring structure starts to develop another inner layer.
Collapse
Affiliation(s)
- Haruka Yasuoka
- Department of Mechanical Engineering, Osaka Prefecture University, 1-1 Gakuen-cho, Naka-ku, Sakai, Osaka 599-8531 Japan
| | - Ryo Takahama
- Department of Mechanical Engineering, Osaka Prefecture University, 1-1 Gakuen-cho, Naka-ku, Sakai, Osaka 599-8531 Japan
| | - Masayuki Kaneda
- Department of Mechanical Engineering, Osaka Prefecture University, 1-1 Gakuen-cho, Naka-ku, Sakai, Osaka 599-8531 Japan
| | - Kazuhiko Suga
- Department of Mechanical Engineering, Osaka Prefecture University, 1-1 Gakuen-cho, Naka-ku, Sakai, Osaka 599-8531 Japan
| |
Collapse
|
37
|
Bhadauria R, Sanghi T, Aluru NR. Interfacial friction based quasi-continuum hydrodynamical model for nanofluidic transport of water. J Chem Phys 2015; 143:174702. [DOI: 10.1063/1.4934678] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.2] [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
| | - Tarun Sanghi
- 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
|
38
|
Chen S, Wang H, Qian T, Sheng P. Determining hydrodynamic boundary conditions from equilibrium fluctuations. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2015; 92:043007. [PMID: 26565332 DOI: 10.1103/physreve.92.043007] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/24/2014] [Indexed: 06/05/2023]
Abstract
The lack of a first-principles derivation has made the hydrodynamic boundary condition a classical issue for the past century. The fact that the fluid can have interfacial structures adds additional complications and ambiguities to the problem. Here we report the use of molecular dynamics to identify from equilibrium thermal fluctuations the hydrodynamic modes in a fluid confined by solid walls, thereby extending the application of the fluctuation-dissipation theorem to yield not only the accurate location of the hydrodynamic boundary at the molecular scale, but also the relevant parameter value(s) for the description of the macroscopic boundary condition. We present molecular dynamics results on two examples to illustrate the application of this approach-one on the hydrophilic case and one on the hydrophobic case. It is shown that the use of the orthogonality condition of the modes can uniquely locate the hydrodynamic boundary to be inside the fluid in both cases, separated from the molecular solid-liquid interface by a small distance Δ that is a few molecules in size. The eigenvalue equation of the hydrodynamic modes directly yields the slip length, which is about equal to Δ in the hydrophilic case but is larger than Δ in the hydrophobic case. From the decay time we also obtain the bulk viscosity which is in good agreement with the value obtained from dynamic simulations. To complete the picture, we derive the Green-Kubo relation for a finite fluid system and show that the boundary fluctuations decouple from the bulk only in the infinite-fluid-channel limit; and in that limit we recover the interfacial fluctuation-dissipation theorem first presented by Bocquet and Barrat. The coupling between the bulk and the boundary fluctuations provides both the justification and the reason for the effectiveness of the present approach, which promises broad utility for probing the hydrodynamic boundary conditions relevant to structured or elastic interfaces, as well as two-phase immiscible flows.
Collapse
Affiliation(s)
- Shuyu Chen
- Department of Physics, Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong, China
| | - Han Wang
- CAEP Software Center for High Performance Numerical Simulation, Beijing, China
- Zuse Institute Berlin, Takustrasse 7, 14195 Berlin, Germany
| | - Tiezheng Qian
- Department of Mathematics, Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong, China
| | - Ping Sheng
- Department of Physics, Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong, China
- Institute for Advanced Study, Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong, China
| |
Collapse
|
39
|
Babu JS, Uday S, Sekhar S, Sathian SP. Physicochemical analysis of slip flow phenomena in liquids under nanoscale confinement. THE EUROPEAN PHYSICAL JOURNAL. E, SOFT MATTER 2015; 38:109. [PMID: 26490250 DOI: 10.1140/epje/i2015-15109-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/11/2015] [Revised: 09/12/2015] [Accepted: 09/15/2015] [Indexed: 06/05/2023]
Abstract
Eyring theory employs the statistical mechanical theory of absolute reaction rates to analyse the transport mechanisms in fluids. A physicochemical methodology combining molecular dynamics (MD) and Eyring theory of reaction rates is proposed for investigating the liquid slip on a solid wall in the nanoscale domain. The method involves the determination of activation energy required for the flow process directly from the MD trajectory information and then calculate the important transport properties of the confined fluid from the activation energy. In order to demonstrate the universal applicability of the proposed methodology in nanofluidics, the slip flow behavior of argon, water and ionic liquid confined in various nanostructures has been investigated. The slip length is found to be size dependent in all the cases. The novelty of this method is that the variations in slip length are explained on the basis of molecular interactions and the subsequent changes in the activation energy.
Collapse
Affiliation(s)
- Jeetu S Babu
- Computational Nanotechnology Laboratory, School of Nano Science and Technology, National Institute of Technology Calicut, 673601, Kozhikode, India
| | - Swathi Uday
- Computational Nanotechnology Laboratory, School of Nano Science and Technology, National Institute of Technology Calicut, 673601, Kozhikode, India
| | - Suneeth Sekhar
- Computational Nanotechnology Laboratory, School of Nano Science and Technology, National Institute of Technology Calicut, 673601, Kozhikode, India
| | - Sarith P Sathian
- Department of Applied Mechanics, Indian Institute of Technology Madras, 600036, Chennai, India.
| |
Collapse
|
40
|
Zhang Y. Novel nano bearings constructed by physical adsorption. Sci Rep 2015; 5:14539. [PMID: 26412488 PMCID: PMC4585955 DOI: 10.1038/srep14539] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2015] [Accepted: 08/25/2015] [Indexed: 12/02/2022] Open
Abstract
The paper proposes a novel nano bearing formed by the physical adsorption of the confined fluid to the solid wall. The bearing is formed between two parallel smooth solid plane walls sliding against one another, where conventional hydrodynamic lubrication theory predicted no lubricating effect. In this bearing, the stationary solid wall is divided into two subzones which respectively have different interaction strengths with the lubricating fluid. It leads to different physical adsorption and slip properties of the lubricating fluid at the stationary solid wall respectively in these two subzones. It was found that a significant load-carrying capacity of the bearing can be generated for low lubricating film thicknesses, because of the strong physical adsorption and non-continuum effects of the lubricating film.
Collapse
Affiliation(s)
- Yongbin Zhang
- College of Mechanical Engineering, Changzhou University, Changzhou, 213016, Jiangsu Province, China
| |
Collapse
|
41
|
|
42
|
Li L, Mo J, Li Z. Flow and slip transition in nanochannels. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2014; 90:033003. [PMID: 25314525 DOI: 10.1103/physreve.90.033003] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/11/2014] [Indexed: 06/04/2023]
Abstract
We experimentally investigate the Poiseuille flows in nanochannels. It is found that the flow rate undergoes a transition between two linear regimes as the shear rate is varied. The transition indicates that the nonslip boundary condition is valid at low shear rate. When the shear rate is larger than a critical value, slip takes place and the slip length increases linearly with increasing shear rate before approaching a constant value. The results reported in this work can help advance the understanding of flow slip in nanochannels.
Collapse
Affiliation(s)
- Long Li
- Department of Mechanical and Aerospace Engineering, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
| | - Jingwen Mo
- Department of Mechanical and Aerospace Engineering, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
| | - Zhigang Li
- Department of Mechanical and Aerospace Engineering, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong
| |
Collapse
|
43
|
Agnihotri MV, Chen SH, Beck C, Singer SJ. Displacements, Mean-Squared Displacements, and Codisplacements for the Calculation of Nonequilibrium Properties. J Phys Chem B 2014; 118:8170-8. [DOI: 10.1021/jp5012523] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Affiliation(s)
- Mithila V. Agnihotri
- Biophysics Program and ‡Department of
Chemistry and Biochemistry, Ohio State University, Columbus, Ohio 43210, United States
| | - Si-Han Chen
- Biophysics Program and ‡Department of
Chemistry and Biochemistry, Ohio State University, Columbus, Ohio 43210, United States
| | - Corey Beck
- Biophysics Program and ‡Department of
Chemistry and Biochemistry, Ohio State University, Columbus, Ohio 43210, United States
| | - Sherwin J. Singer
- Biophysics Program and ‡Department of
Chemistry and Biochemistry, Ohio State University, Columbus, Ohio 43210, United States
| |
Collapse
|
44
|
De Luca S, Todd BD, Hansen JS, Daivis PJ. Molecular dynamics study of nanoconfined water flow driven by rotating electric fields under realistic experimental conditions. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2014; 30:3095-3109. [PMID: 24575940 DOI: 10.1021/la404805s] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/03/2023]
Abstract
In our recent work, J. Chem. Phys. 2013, 138, 154712, we demonstrated the feasibility of unidirectional pumping of water, exploiting translational-rotational momentum coupling using nonequilibrium molecular dynamics simulations. Flow can be sustained when the fluid is driven out of equilibrium by an external spatially uniform rotating electric field and confined between two planar surfaces exposing different degrees of hydrophobicity. The permanent dipole moment of water follows the rotating field, thus inducing the molecules to spin, and the torque exerted by the field is continuously injected into the fluid, enabling a steady conversion of spin angular momentum into linear momentum. The translational-rotational coupling is a sensitive function of the rotating electric field parameters. In this work, we have found that there exists a small energy dissipation region attainable when the frequency of the rotating electric field matches the inverse of the dielectric relaxation time of water and when its amplitude lies in a range just before dielectric saturation effects take place. In this region, that is, when the frequency lies in a small window of the microwave region around ∼20 GHz and amplitude ∼0.03 V Å(-1), the translational-rotational coupling is most effective, yielding fluid velocities of magnitudes of ∼2 ms(-1) with only moderate fluid heating. In this work, we also confine water to a realistic nanochannel made of graphene giving a hydrophobic surface on one side and β-cristobalite giving a hydrophilic surface on the other, reproducing slip-and-stick velocity boundary conditions, respectively. This enables us to demonstrate that in a realistic environment, the coupling can be effectively exploited to achieve noncontact pumping of water at the nanoscale. A quantitative comparison between nonequilibrium molecular dynamics and analytical solutions of the extended Navier-Stokes equations, including an external rotating electric field has been performed, showing excellent agreement when the electric field parameters match the aforementioned small energy dissipation region.
Collapse
Affiliation(s)
- Sergio De Luca
- Department of Mathematics, Faculty of Science, Engineering and Technology, and Centre for Molecular Simulation, Swinburne University of Technology , Melbourne, Victoria 3122, Australia
| | | | | | | |
Collapse
|
45
|
Huang K, Szlufarska I. Green-Kubo relation for friction at liquid-solid interfaces. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2014; 89:032119. [PMID: 24730802 DOI: 10.1103/physreve.89.032119] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/19/2013] [Indexed: 06/03/2023]
Abstract
We have developed a Green-Kubo relation that enables accurate calculations of friction at solid-liquid interfaces directly from equilibrium molecular dynamics (MD) simulations and that provides a pathway to bypass the time-scale limitations of typical nonequilibrium MD simulations. The theory has been validated for a number of different interfaces and it is demonstrated that the liquid-solid slip is an intrinsic property of an interface. Because of the high numerical efficiency of our method, it can be used in the design of interfaces for applications in aqueous environments, such as nano- and microfluidics.
Collapse
Affiliation(s)
- Kai Huang
- Materials Science Program, University of Wisconsin, Madison, Wisconsin 53706-1595, USA
| | - Izabela Szlufarska
- Materials Science Program, University of Wisconsin, Madison, Wisconsin 53706-1595, USA and Department of Materials Science and Engineering, University of Wisconsin, Madison, Wisconsin 53706-1595, USA
| |
Collapse
|
46
|
Ritos K, Mattia D, Calabrò F, Reese JM. Flow enhancement in nanotubes of different materials and lengths. J Chem Phys 2014; 140:014702. [DOI: 10.1063/1.4846300] [Citation(s) in RCA: 75] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022] Open
|
47
|
Bhatia SK, Nicholson D. Friction between solids and adsorbed fluids is spatially distributed at the nanoscale. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2013; 29:14519-14526. [PMID: 24168469 DOI: 10.1021/la403445j] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/02/2023]
Abstract
The widespread developments in the use of nanomaterials in catalysis, adsorption, and nanofluidics present significant new challenges in achieving optimal adsorbed fluid flow characteristics. Here we demonstrate, using molecular dynamics simulations of nanoconfined fluids, that at nanoscales, fluid-solid friction is not restricted to a sharp interface as is commonly assumed; instead it is distributed over the whole adsorbed fluid phase, and is strongest in an interfacial region that is not negligible in comparison to the system size. Our simulations yield position-dependent dynamical fluid-solid friction coefficients, and lead to a modification of conventional hydrodynamics, incorporating distributed momentum loss in the fluid due to fluid-solid interaction. The results demonstrate that the usual concepts of slip length or interfacial friction coefficient are meaningful only for uniform fluids, and lose their significance for adsorbates in nanospaces, which are intrinsically inhomogeneous. We show that static friction coefficients, based on equilibrium density distributions, follow the same spatial dependence as the dynamical coefficients. These results open up possibilities for tailoring nanomaterials and surfaces to engineer low friction pathways for adsorbed fluid flow by tuning the potential energy landscape.
Collapse
Affiliation(s)
- Suresh K Bhatia
- School of Chemical Engineering, The University of Queensland Brisbane QLD 4072, Australia
| | | |
Collapse
|
48
|
De Luca S, Todd BD, Hansen JS, Daivis PJ. Electropumping of water with rotating electric fields. J Chem Phys 2013; 138:154712. [PMID: 23614441 DOI: 10.1063/1.4801033] [Citation(s) in RCA: 38] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Pumping of fluids confined to nanometer dimension spaces is a technically challenging yet vitally important technological application with far reaching consequences for lab-on-a-chip devices, biomimetic nanoscale reactors, nanoscale filtration devices and the like. All current pumping mechanisms require some sort of direct intrusion into the nanofluidic system, and involve mechanical or electronic components. In this paper, we present the first nonequilibrium molecular dynamics results to demonstrate that non-intrusive electropumping of liquid water on the nanoscale can be performed by subtly exploiting the coupling of spin angular momentum to linear streaming momentum. A spatially uniform rotating electric field is applied to water molecules, which couples to their permanent electric dipole moments. The resulting molecular rotational momentum is converted into linear streaming momentum of the fluid. By selectively tuning the degree of hydrophobicity of the solid walls one can generate a net unidirectional flow. Our results for the linear streaming and angular velocities of the confined water are in general agreement with the extended hydrodynamical theory for this process, though also suggest refinements to the theory are required. These numerical experiments confirm that this new concept for pumping of polar nanofluids can be employed under laboratory conditions, opening up significant new technological possibilities.
Collapse
Affiliation(s)
- Sergio De Luca
- Mathematics, Faculty of Engineering and Industrial Sciences, and Centre for Molecular Simulation, Swinburne University of Technology, Melbourne, Victoria 3122, Australia
| | | | | | | |
Collapse
|
49
|
Kannam SK, Todd BD, Hansen JS, Daivis PJ. How fast does water flow in carbon nanotubes? J Chem Phys 2013; 138:094701. [PMID: 23485316 DOI: 10.1063/1.4793396] [Citation(s) in RCA: 206] [Impact Index Per Article: 18.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023] Open
Abstract
The purpose of this paper is threefold. First, we review the existing literature on flow rates of water in carbon nanotubes. Data for the slip length which characterizes the flow rate are scattered over 5 orders of magnitude for nanotubes of diameter 0.81-10 nm. Second, we precisely compute the slip length using equilibrium molecular dynamics (EMD) simulations, from which the interfacial friction between water and carbon nanotubes can be found, and also via external field driven non-equilibrium molecular dynamics simulations (NEMD). We discuss some of the issues in simulation studies which may be reasons for the large disagreements reported. By using the EMD method friction coefficient to determine the slip length, we overcome the limitations of NEMD simulations. In NEMD simulations, for each tube we apply a range of external fields to check the linear response of the fluid to the field and reliably extrapolate the results for the slip length to values of the field corresponding to experimentally accessible pressure gradients. Finally, we comment on several issues concerning water flow rates in carbon nanotubes which may lead to some future research directions in this area.
Collapse
Affiliation(s)
- Sridhar Kumar Kannam
- Mathematics Discipline, Faculty of Engineering and Industrial Science, and Centre for Molecular Simulation, Swinburne University of Technology, Melbourne, Victoria 3122, Australia.
| | | | | | | |
Collapse
|
50
|
Sun J, Wang W, Wang HS. Dependence between velocity slip and temperature jump in shear flows. J Chem Phys 2013; 138:234703. [DOI: 10.1063/1.4810810] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|