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

Full description

Saved in:
Bibliographic Details
Published in:The Journal of chemical physics Vol. 134; no. 2; p. 024526
Main Authors: Thomas, Jason C, Rowley, Richard L
Format: Journal Article
Language:English
Published: United States 14-01-2011
Subjects:
Online Access:Get more information
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary: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.
ISSN:1089-7690
DOI:10.1063/1.3543798