Shear viscosity of simple fluids in porous media: molecular dynamic simulations and correlation models

Equilibrium molecular dynamic simulations have been used to calculate the shear viscosity of liquid argon in macrovolume system and in porous media at different temperatures, densities and pore widths. On the other hand, based on the Chapman–Enskog theory and Heyes relationships, two correlation models which can describe the viscosity of simple liquids in porous media are proposed as a function of the reduced temperature, density and pore width. The validity of the models is evaluated by comparing the calculated viscosity to simulation data.     

Viscosity of confined inhomogeneous nonequilibrium fluids

We use the nonlocal linear hydrodynamic constitutive model, proposed by Evans and Morriss [Statistical Mechanics of Nonequilibrium Liquids (Academic, London, 1990)], for computing an effective spatially dependent shear viscosity of inhomogeneous nonequilibrium fluids. The model is applied to a simple atomic fluid undergoing planar Poiseuille flow in a confined channel of several atomic diameters width. We compare the spatially dependent viscosity with a local generalization of Newton's law of viscosity and the Navier-Stokes viscosity, both of which are known to suffer extreme inaccuracies for highly inhomogeneous systems. The nonlocal constitutive model calculates effective position dependent viscosities that are free from the notorious singularities experienced by applying the commonly used local constitutive model. It is simple, general, and has widespread applicability in nanofluidics where experimental measurement of position dependent transport coefficients is currently inaccessible. In principle the method can be used to predict approximate flow profiles of any arbitrary inhomogeneous system. We demonstrate this by predicting the flow profile for a simple fluid undergoing planar Couette flow in a confined channel of several atomic diameters width.     

Shear-dependent viscosity in simple fluids

Computer simulation results are used to examine the shear-dependent viscosity of simple fluids. The prime purpose of the investigation was to examine the possible non-newtonian effects when the particles in a fluid interact through central forces. A range of spherical interactions was examined, including both potentials with a barrier and simple colloidal models, and in all cases the fluids were shear-thinning. The non-equilibrium radial distribution function for a soft potential is considerably more sensitive to shear rate than is that for the strongly repulsive model. Calculations of the hydrostatic pressure show that it increases with shear rate, demonstrating that the phenomena of shear thickening and shear dilatency — used interchangeably in the literature — are separate phenomena.     

First-order mean spherical approximation for inhomogeneous fluids

The first-order mean-spherical approximation (FMSA) [Y. Tang, J. Chem. Phys., 118, 4140 (2003)] is extended to the studies of inhomogeneous fluids by combining with Rosenfeld's perturbative method [Y. Rosenfeld, J. Chem. Phys. 98, 8126 (1993)]. In the extension, the key input-direct correlation function of FMSA-is applied to constructing the free energy density functional. Preserving its high fidelity at the bulk limit, the FMSA shows satisfactory performance for Yukawa fluids near hard and attractive walls. The results are better than or comparable to several other theories reported before for the geometry. The FMSA is found, in particular, more satisfactory than the traditional mean-field theory for predicting density profiles around hard walls. The FMSA is also compared with the full MSA for inhomogeneous fluids, showing no appreciable differences. The inhomogeneous FMSA goes successfully through the self-consistency test for reproducing the radial distribution function of the bulk Yukawa fluid. As far as the computation is concerned, the FMSA can be executed much faster than any nonmean-field theories, and the speed is virtually identical to that of the mean-field theory.     

Model for dynamics of inhomogeneous and bulk fluids

An accurate model for the density of states (DOS) for strongly inhomogeneous and bulk fluids has been proposed based on gamma distributions. The contribution to the density of states from the collective dynamics is modeled as an incomplete gamma distribution and the high frequency region is obtained from the solution of the memory equation using a sech memory kernel. Using only the frequency moments as input, the model parameters for the collective dynamics are obtained by matching moments of the resulting distribution. The model results in an analytical expression for the self-diffusivity of the fluid. We present results for soft sphere fluids confined in slit-shaped pores as well as bulk soft sphere liquids. Comparisons of the DOS, velocity autocorrelation functions, and memory kernels with molecular dynamics simulations reveal that the model predicts features in the DOS over the entire frequency range and is able to capture changes in the DOS as a function of fluid density and temperature. As a result the predicted VACFs, memory kernels, and self-diffusivities are accurately predicted over a wide range of conditions. Since the frequency moments for bulk liquids can be obtained from pair correlation functions, our method provides a direct route from fluid structure to dynamics. For fluids confined in slit-shaped pores, where the frequency moments are obtained from molecular dynamics simulations, the predicted self-diffusivities capture the resulting oscillations due to variations in the solvation pressure, and in the case of smooth walled pores, the predictions are superior to those obtained using kinetic theory.     

Density functional theory for inhomogeneous associating chain fluids

We propose a nonlocal density functional theory for associating chain molecules. The chains are modeled as tangent spheres, which interact via Lennard-Jones (12,6) attractive interactions. A selected segment contains additional, short-ranged, highly directional interaction sites. The theory incorporates an accurate treatment of the chain molecules via the intramolecular potential formalism and should accurately describe systems with strongly varying external fields, e.g., attractive walls. Within our approach we investigate the structure of the liquid-vapor interface and capillary condensation of a simple model of associating chains with only one associating site placed on the first segment. In general, the properties of inhomogeneous associating chains depend on the association energy. Similar to the bulk systems we find the behavior of associating chains of a given length to be in between that for the nonassociating chains of the same length and that for the nonassociating chains twice as large.     

Shear viscosity of dispersions of particles with liquid shells

Shear viscosity of the dispersions of particles with a liquid shell is calculated on the basis of cell method. The problem of choosing an optimal radius of the cell is considered. Exact expression for the shear viscosity of dispersion, which is suitable at the arbitrary volume fraction of dispersed phase, is derived.     

Theoretical scheme for the shear viscosity of Lennard-Jones fluids

We use the shear viscosity expression from the Enskog theory of dense gases in a perturbative scheme for the Lennard-Jones (LJ) fluid. This perturbative scheme is formulated by combining the analytic rational function approximation method of Bravo Yuste and Santos [Phys. Rev. A 43, 5418 (1991)] for the radial distribution function of hard-sphere fluids and the well known Mansoori-Canfield/Rasaiah-Stell perturbation theory to determine an effective diameter for the LJ fluid. The scheme is reliable on a wide range of temperatures and densities, and is very accurate around the critical point. Using this information, we build an accurate empirical formula for the shear viscosity in the liquid phase, which fits the recent data [K. Meier et al., J. Chem. Phys. 121, 3671 (2004)] in the whole simulation range.     

Transient molecular dynamics simulations of viscosity for simple fluids

A transient molecular dynamics (TMD) method has been developed for simulation of fluid viscosity. In this method a sinusoidal velocity profile is instantaneously overlaid onto equilibrated molecular velocities, and the subsequent decay of that velocity profile is observed. The viscosity is obtained by matching in a least-squares sense the analytical solution of the corresponding momentum transport boundary-value problem to the simulated decay of the initial velocity profile. The method was benchmarked by comparing results obtained from the TMD method for a Lennard-Jones fluid with those previously obtained using equilibrium molecular dynamics (EMD) simulations. Two different constitutive models were used in the macroscopic equations to relate the shear rate to the stress. Results using a Newtonian fluid model agree with EMD results at moderate densities but exhibit an increasingly positive error with increasing density at high densities. With the initial velocity profiles used in this study, simulated transient velocities displayed clear viscoelastic behavior at dimensionless densities above 0.7. However, the use of a linear viscoelastic model reproduces the simulated transient velocity behavior well and removes the high-density bias observed in the results obtained under the assumption of Newtonian behavior. The viscosity values obtained using the viscoelastic model are in excellent agreement with the EMD results over virtually the entire fluid domain. For simplicity, the Newtonian fluid model can be used at lower densities and the viscoelastic model at higher densities; the two models give equivalent results at intermediate densities.     

Poiseuille flow to measure the viscosity of particle model fluids

The most important property of a fluid is its viscosity, it determines the flow properties. If one simulates a fluid using a particle model, calculating the viscosity accurately is difficult because it is a collective property. In this article we describe a new method that has a better signal to noise ratio than existing methods. It is based on using periodic boundary conditions to simulate counter-flowing Poiseuille flows without the use of explicit boundaries. The viscosity is then related to the mean flow velocity of the two flows. We apply the method to two quite different systems. First, a simple generic fluid model, dissipative particle dynamics, for which accurate values of the viscosity are needed to characterize the model fluid. Second, the more realistic Lennard-Jones fluid. In both cases the values we calculated are consistent with previous work but, for a given simulation time, they are more accurate than those obtained with other methods.     

Dependence of the viscosity coefficients of magnetic fluids on parameters of state

Dependences of the viscosity coefficients of magnetic fluids on parameters of state were investigated numerically using previously derived dynamic equations. It was shown that the volume viscosity and shear viscosity coefficients of a magnetic fluid based on kerosene increase with increasing density and concentration and decrease with increasing temperature; the coefficients increase with an increase in the magnetic field gradient. The results obtained are in satisfactory agreement with the experimental data.     

A viscosity equation for polyatomic fluids under normal and high pressures

In this work, we relate the self-diffusion coefficient to the residual entropy of the system according to the free volume theory and scaling principle. The viscosity equation for a freely jointed Lennard-Jones chain fluid is then obtained from the expression of self-diffusion coefficient by applying the Stokes–Einstein equation. The real polyatomic compounds are modeled as chains of tangent Lennard-Jones segments. The segment size and energy parameter as well as chain length (expressed by the number of segments) are obtained from the experimental viscosity data. The proposed viscosity equation reproduces the experimental viscosity data with an average absolute deviation of 5.12% for 18 polyatomic compounds (1600 data points) over wide ranges of temperature and pressure. For engineering applications, the generalized model parameters for normal alkanes with the number of carbon atoms n > 3 are proposed. The segment energy parameter is suggested to be evaluated from the critical temperature, and the segment size parameter and chain length are correlated with the number of carbon atoms in an alkane molecule.     

Shear and longitudinal viscosity of non-ionic C12E8 aqueous solutions

We present an extensive set of measurements of steady shear viscosity (eta degrees(s)), longitudinal elastic modulus (M'), and ultrasonic absorption (alpha) in the one-phase isotropic liquid region of the non-ionic surfactant C12E8 aqueous solutions. Within a given temperature interval, this phase extends along the entire surfactant concentration range that could be fully covered in the experiments. In agreement with previous studies, the overall results support the presence of two separated intervals of concentration corresponding to different structural properties. In the surfactant-rich region the temperature dependence of eta degrees(s) follows an equation characteristic of glass-like systems. The ultrasonic absorption spectra show unambiguous evidence of viscoelastic behavior that can be described by a Cole-Cole relaxation formula. In this region, when both the absorption and the frequency are scaled by the static shear viscosity (eta degrees(s)), the scaled attenuation reduces to a single universal curve for all temperatures and concentrations. In the water-rich region the behavior of eta degrees(s), M', and alpha are more complex and reflect the presence of dispersed aggregates whose size increases with temperature and concentration. At these concentrations the ultrasonic spectra are characterized by a multiple decay rate. The high-frequency tail falls in the same frequency range seen at high surfactant content and exhibits similar behaviors. This contribution is ascribed to the mixture of hydrophilic terminations and water present at the micellar interfaces that resembles the condition of a concentrated polymer solution. An additional low-frequency contribution is also observed, which is ascribed to the exchange of water molecules and/or surfactant monomers between the aggregates and the bulk solvent region.     

Shear viscosity of molten alkali halides from equilibrium and nonequilibrium molecular-dynamics simulations

The shear viscosity of molten NaCl and KCl was calculated through equilibrium (EMD) and nonequilibrium molecular-dynamics (NEMD) simulations in the canonical (N,V,T) ensemble. Two rigid-ion potentials were investigated, namely, the Born-Mayer-Huggins-Tosi-Fumi potential and the Michielsen-Woerlee-Graaf-Ketelaar potential with the parameters proposed by Ladd. The NEMD simulations were performed using the SLLOD equations of motion [D. J. Evans and G. P. Morriss, Phys. Rev. A 30, 1528 (1984)] with a Gaussian isokinetic thermostat and the results are compared with those obtained from Green-Kubo EMD (N,V,T) simulations and experimental shear viscosity data. The NEMD zero strain rate shear viscosity, eta(0), was obtained by fitting a simplified Carreau-type equation and by application of mode-coupling theory, i.e., a eta-gamma(1/2) linear relationship. The values obtained from the first method are found to be significantly lower than those predicted by the second. The agreement between the EMD and NEMD results with experimental data is satisfactory for the two potentials investigated. The ion-ion radial distribution functions obtained with the two rigid-ion potentials for both molten salts are discussed in terms of the differences between the two models.     

Transient molecular dynamics simulations of liquid viscosity for nonpolar and polar fluids

A transient molecular dynamics (TMD) method for obtaining fluid viscosity is extended to multisite, force-field models of both nonpolar and polar liquids. The method overlays a sinusoidal velocity profile over the peculiar particle velocities and then records the transient decay of the velocity profile. The viscosity is obtained by regression of the solution of the momentum equation with an appropriate constitutive equation and initial and boundary conditions corresponding to those used in the simulation. The transient velocity decays observed appeared to include both relaxation and retardation effects. The Jeffreys viscoelastic model was found to model accurately the transient responses obtained for multisite models for n-butane, isobutane, n-hexane, water, methanol, and 1-hexanol. TMD viscosities obtained for saturated liquids over a wide range of densities agreed well for the polar fluids, both with nonequilibrium molecular dynamics (NEMD) results using the same force-field models and with correlations based on experimental data. Viscosities obtained for the nonpolar fluids agreed well with the experimental and NEMD results at low to moderate densities, but underpredicted experimental values at higher densities where shear-thinning effects and viscous heating may impact the TMD simulations.     

Correlation analysis of the power law parameters for viscosity of some engineering fluids

Knowledge and estimation of transport properties of fluids which are sensitive to temperature variation like viscosity are necessary in mass flow and heat transfer computation. In the present work, based on the use of econometric and statistical techniques for regression analysis and correlation tests, we propose an original equation modelling the relationship between the two parameters of viscosity power law equation. Empirical validation using data set of some pure fluids provided from the literature gives excellent statistical results which allow us to redefine the power law equation using only a single parameter. This result is important in fluids engineering since the validation of the proposed equation simplifies the estimation of viscous behaviour and the ensuing calculations by reducing the number of viscosity equation parameters and facilitating manipulations for fluids with low and moderate viscosity in engineering data which permits to estimate one unavailable parameter when the second one is available.     

Rotational viscosity of fluids composed of linear molecules: an equilibrium molecular dynamics study

In this paper, we investigate the rotational viscosity for a chlorine fluid and for a fluid composed of small linear molecules by using equilibrium molecular dynamics simulations. The rotational viscosity is calculated over a large range of state points. It is found that the rotational viscosity is almost independent of temperature in the range studied here but exhibits a power-law dependency on density. The rotational viscosity also shows a power-law relationship with the molecular length, and the ratio between the shear and rotational viscosities approaches 0.5 for the longest molecule studied here. By changing the number of atoms or united atomic units per molecule and by keeping the molecule length fixed, we show that fluids composed of molecules which have a rodlike shape have a lower rotational viscosity. We argue that this phenomenon is due to the reduction in intermolecular connectivity, which leads to larger fluctuations around the values possessed by the fluid on average. The conclusions here can be extended to fluids composed of uniaxial molecules of arbitrary length.