Shear viscosity of inhomogeneous fluids

Statistical physics of inhomogeneous transport: Unification of diffusion laws and inference from first-passage statistics

Characterization of composite materials, whose properties vary in space over microscopic scales, has become a problem of broad interdisciplinary interest. In particular, estimation of the inhomogeneous transport coefficients, e.g., the diffusion coefficient or the heat conductivity, which shape important processes in biology and engineering, is a challenging task. The analysis of such systems is further complicated because two alternative formulations of the inhomogeneous transport equations exist in the literature-the Smoluchowski and Fokker-Planck equations, which are also related to the so-called Ito-Stratonovich dilemma. Using the theory of statistical physics, we show that the two formulations, usually regarded as distinct models, are physically equivalent. From this result we develop efficient estimates for the transverse space-dependent diffusion coefficient in fluids near a phase boundary. Our method requires only measurements of escape probabilities and mean exit times of molecules leaving a narrow spatial region. We test our estimates in three case studies: (i) a Langevin model of a Büttikker-Landauer ratchet; atomistic molecular-dynamics simulations of liquid-water molecules in contact with (ii) vapor, and (iii) soap (surfactant) film which has promising applications in physical chemistry. Our analysis reveals that near the surfactant monolayer the mobility of water molecules is slowed down almost twice with respect to the bulk liquid. Moreover, the diffusion coefficient of water correlates with the transition from hydrophilic to hydrophobic parts of the film.

New Framework for Computing a General Local Self-Diffusion Coefficient Using Statistical Mechanics

Widely applicable, modified Green-Kubo expressions for the local diffusion coefficient (D_l) are obtained using linear response theory. In contrast to past definitions in use, these expressions are statistical mechanical results. Molecular simulations of systems with anisotropic diffusion and an inhomogeneous density profile confirm the validity of the results. Diffusion coefficients determined from different expressions in terms of currents and velocity correlations agree in the limit of large systems. Furthermore, they apply to arbitrarily small local regions, making them readily applicable to nanoscale and inhomogeneous systems where knowledge of D_l is important.

Enhanced local viscosity around colloidal nanoparticles probed by equilibrium molecular dynamics simulations

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.

Wall friction should be decoupled from fluid viscosity for the prediction of nanoscale flow

The accurate determination of fluid viscosity based on the microscopic information of molecules is very crucial for the prediction of nanoscale flow. Despite the challenge of this problem, researchers have done a lot of meaningful work and developed several distinctive methods. However, one of the common approaches to calculate the fluid viscosity is using the Green-Kubo formula by considering all the fluid molecules in nanospace, inevitably causing the involvement of the frictional interaction between fluid and the wall into the fluid viscosity. This practice is certainly not appropriate because viscosity is essentially related only to the interactions among fluid molecules. Here, we clarify that the wall friction should be decoupled from fluid viscosity by distinguishing the frictional region and the viscous region for the accurate prediction of nanoscale flow. By comparing the fluid viscosities calculated from the Green-Kubo formula in the whole region and viscous region and the viscosity obtained from the velocity profile through the Hagen-Poiseuille equation, it is found that only the calculated viscosity in the viscous region agrees well with the viscosity from the velocity profile. To demonstrate the applicability of this clarification, the Lennard-Jones fluid and water confined between Lennard-Jones, graphene, and silica walls, even with different fluid-wall interactions, are extensively tested. This work clearly defines the viscosity of fluids at nanoscales from the inherent nature of physics, aiming at the accurate prediction of nanoscale flow from the classical continuum hydrodynamic theory.

Non-local viscosity from the Green-Kubo formula

We study through MD simulations the correlation matrix of the discrete transverse momentum density field in real space for an unconfined Lennard-Jones fluid at equilibrium. Mori theory predicts this correlation under the Markovian approximation from the knowledge of the non-local shear viscosity matrix, which is given in terms of a Green-Kubo formula. However, the running Green-Kubo integral for the non-local shear viscosity does not have a plateau. By using a recently proposed correction for the Green-Kubo formula that eliminates the plateau problem [Español et al., Phys. Rev. E 99, 022126 (2019)], we unambiguously obtain the actual non-local shear viscosity. The resulting Markovian equation, being local in time, is not valid for very short times. We observe that the Markovian equation with non-local viscosity gives excellent predictions for the correlation matrix from a time at which the correlation is around 80% of its initial value. A local in space approximation for the viscosity gives accurate results only after the correlation has decayed to 40% of its initial value.

Challenges in Nanofluidics—Beyond Navier–Stokes at the Molecular Scale

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.

Pressure gradients fail to predict diffusio-osmosis

We present numerical simulations of diffusio-osmotic flow, i.e. the fluid flow generated by a concentration gradient along a solid-fluid interface. In our study, we compare a number of distinct approaches that have been proposed for computing such flows and compare them with a reference calculation based on direct, non-equilibrium molecular dynamics simulations. As alternatives, we consider schemes that compute diffusio-osmotic flow from the gradient of the chemical potentials of the constituent species and from the gradient of the component of the pressure tensor parallel to the interface. We find that the approach based on treating chemical potential gradients as external forces acting on various species agrees with the direct simulations, thereby supporting the approach of Marbach et al (2017 J. Chem. Phys. 146 194701). In contrast, an approach based on computing the gradients of the microscopic pressure tensor does not reproduce the direct non-equilibrium results.

Quantitative relationship between fluid inhomogeneities and flow enhancement in nanotubes

Flow enhancement in nanotubes is of great potential to achieve ultra-fast fluidic transport. However, the mechanism of such a fast transport and the reduction as the tube enlarges to bulk scale is still unclear. In this study, we establish a model to quantitatively correlate the flow enhancement and the fluid inhomogeneity to describe the enhanced transport and its evolution with the tube dimension. We found the fluid inhomogeneity at the solid-liquid interface in nanotubes and its independence with tube size by dissipative particle dynamics (DPD) simulation. Based on that, we establish novel theoretical models for the penetration rate in nanotubes for the first time with parameters related to the fluid inhomogeneity which can achieve quantitative prediction of nanoflow enhancement and are valid through all scales.

Continuum Nanofluidics

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.

How Local and Average Particle Diffusivities of Inhomogeneous Fluids Depend on Microscopic Dynamics

Computer simulations and a stochastic Fokker-Planck equation based approach are used to compare the single-particle diffusion coefficients of equilibrium hard-sphere fluids exhibiting identical inhomogeneous static structure and governed by either Brownian (i.e., overdamped Langevin) or Newtonian microscopic dynamics. The physics of inhomogeneity is explored via the imposition of one-dimensional sinusoidal density profiles of different wavelengths and amplitudes. When imposed density variations are small in magnitude for distances on the scale of a particle diameter, bulk-like average correlations between local structure and mobility are observed. In contrast, when density variations are significant on that length scale, qualitatively different structure-mobility correlations emerge that are sensitive to the governing microscopic dynamics. Correspondingly, a previously proposed scaling between long-time diffusivities for bulk isotropic fluids of particles exhibiting Brownian versus Newtonian dynamics [Pond et al. Soft Matter 2011, 7, 9859-9862] cannot be generalized to describe the position-dependent behaviors of strongly inhomogeneous fluids. While average diffusivities in the inhomogeneous and homogeneous directions are coupled, their qualitative dependencies on inhomogeneity wavelength are sensitive to the details of the microscopic dynamics. Nonetheless, average diffusivities of the inhomogeneous fluids can be approximately predicted for either type of dynamics based on knowledge of bulk isotropic fluid behavior and how inhomogeneity modifies the distribution of available volume. Analogous predictions for average diffusivities of experimental, inhomogeneous colloidal dispersions (based on known bulk behavior) suggest that they will exhibit qualitatively different trends than those predicted by models governed by overdamped Langevin dynamics that do not account for hydrodynamic interactions.

Nonlocal response functions for predicting shear flow of strongly inhomogeneous fluids. I. Sinusoidally driven shear and sinusoidally driven inhomogeneity

We present theoretical expressions for the density, strain rate, and shear pressure profiles in strongly inhomogeneous fluids undergoing steady shear flow with periodic boundary conditions. The expressions that we obtain take the form of truncated functional expansions. In these functional expansions, the independent variables are the spatially sinusoidal longitudinal and transverse forces that we apply in nonequilibrium molecular-dynamics simulations. The longitudinal force produces strong density inhomogeneity, and the transverse force produces sinusoidal shear. The functional expansions define new material properties, the response functions, which characterize the system's nonlocal response to the longitudinal force and the transverse force. We find that the sinusoidal longitudinal force, which is mainly responsible for the generation of density inhomogeneity, also modulates the strain rate and shear pressure profiles. Likewise, we find that the sinusoidal transverse force, which is mainly responsible for the generation of sinusoidal shear flow, can also modify the density. These cross couplings between density inhomogeneity and shear flow are also characterized by nonlocal response functions. We conduct nonequilibrium molecular-dynamics simulations to calculate all of the response functions needed to describe the response of the system for weak shear flow in the presence of strong density inhomogeneity up to the third order in the functional expansion. The response functions are then substituted directly into the truncated functional expansions and used to predict the density, velocity, and shear pressure profiles. The results are compared to the directly evaluated profiles from molecular-dynamics simulations, and we find that the predicted profiles from the truncated functional expansions are in excellent agreement with the directly computed density, velocity, and shear pressure profiles.

Structure, thermodynamics, and position-dependent diffusivity in fluids with sinusoidal density variations

Molecular dynamics simulations and a stochastic method based on the Fokker-Planck equation are used to explore the consequences of inhomogeneous density profiles on the thermodynamic and dynamic properties of the hard-sphere fluid and supercooled liquid water. Effects of the inhomogeneity length scale are systematically considered via the imposition of sinusoidal density profiles of various wavelengths. For long-wavelength density profiles, bulk-like relationships between local structure, thermodynamics, and diffusivity are observed as expected. However, for both systems, a crossover in behavior occurs as a function of wavelength, with qualitatively different correlations between the local static and dynamic quantities emerging as density variations approach the scale of a particle diameter. Irrespective of the density variation wavelength, average diffusivities of hard-sphere fluids in the inhomogeneous and homogeneous directions are coupled and approximately correlate with the volume available for insertion of another particle. Unfortunately, a quantitatively reliable static predictor of position-dependent dynamics has yet to be identified for even the simplest of inhomogeneous fluids.

Green-Kubo relation for friction at liquid-solid interfaces

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.

Local shear viscosity of strongly inhomogeneous dense fluids: from the hard-sphere to the Lennard-Jones fluids

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.

Effects of nanoscale density inhomogeneities on shearing fluids

It is well known that density inhomogeneities at the solid-liquid interface can have a strong effect on the velocity profile of a nanoconfined fluid in planar Poiseuille flow. However, it is difficult to control the density inhomogeneities induced by solid walls, making this type of system unsuitable for a comprehensive study of the effect on density inhomogeneity on nanofluidic flow. In this paper, we employ an external force compatible with periodic boundary conditions to induce the density variation, which greatly simplifies the problem when compared to flow in nonperiodic nanoconfined systems. Using the sinusoidal transverse force method to produce shearing velocity profiles and the sinusoidal longitudinal force method to produce inhomogeneous density profiles, we are able to observe the interactions between the two property inhomogeneities at the level of individual Fourier components. This gives us a method for direct measurement of the coupling between the density and velocity fields and allows us to introduce various feedback control mechanisms which customize fluid behavior in individual Fourier components. We briefly discuss the role of temperature inhomogeneity and consider whether local thermal expansion due to nonuniform viscous heating is sufficient to account for shear-induced density inhomogeneities. We also consider the local Newtonian constitutive relation relating the shear stress to the velocity gradient and show that the local model breaks down for sufficiently large density inhomogeneities over atomic length scales.

Local viscosity of a fluid confined in a narrow pore

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).