Dense-fluid shear viscosity via nonequilibrium molecular dynamics

Abstract

A novel fluid-transport calculation by computer simulation, via nonequilibrium molecular dynamics, of laboratory methods of transport measurement is described. Shear viscosity of soft-sphere ($r^{-12\mathrm{}}$ potential) and Lennard-Jones particles ($r^{-12\mathrm{}}-r^{-6\mathrm{}}$ potential) has been obtained from molecular dynamic modeling of Couette flow. Soft-sphere deviations from Enskog theory are similar to those found for hard spheres by Alder, Gass, and Wainwright, using time-correlations of equilibrium molecular dynamic system fluctuations. For the Lennard-Jones shear viscosity near the triple-point region, there is agreement between the equilibrium calculation of Levesque, Verlet, and Kurkijarvi and the nonequilibrium results using 108 atoms in a cube. However, systems two and three cubes wide give lower results, which, when extrapolated with inverse width, yield close agreement with the experimental argon shear viscosity. Comparison of the Lennard-Jones shear viscosity with experimental argon data along the saturated vapor-pressure line of argon confirms our successful simulation of macroscopic viscous flow with few-particle nonequilibrium molecular dynamic systems. A new result of the nonequilibrium molecular dynamics is the characterization of nonequilibrium distribution functions, which might provide the basis for a perturbation theory of transport. Since momentum transport is primarily accomplished by the repulsive potential core for high temperatures, the Lennard-Jones shear viscosity must behave like the soft-sphere system for high temperatures [viscosity divided by ${\left(\mathrm{t}\mathrm{e}\mathrm{m}\mathrm{p}\mathrm{e}\mathrm{r}\mathrm{a}\mathrm{t}\mathrm{u}\mathrm{r}\mathrm{e}\right)}^{\frac{2}{3}}$ is a function of density divided by ${\mathrm{}\left(\mathrm{temperature}\right)\mathrm{}}^{\mathrm{¼}}$]. In fact, the calculated excess shear viscosity (that part above the zero-density temperature dependence) has been successfully correlated in terms of the 12th-power scaling variables for temperatures as low as the critical value (along the freezing line). The utilization of soft-sphere scaling variables yields relatively simple functions for describing both the excess shear viscosity and the thermal-conductivity behavior throughout the fluid phase. The introduction of these scaling variables also clearly reveals two features: (i) weak temperature dependence, and (ii) the sign of the temperature derivative at constant density (negative for shear viscosity and positive for thermal conductivity). While both of these features have been experimentally observed in simple fluid experimental data, their cause has not been previously traced to the dominance of the core potential. Thus, the soft-sphere scaling variables should be useful for correlating experimental data.