Abstract:New calculations have been made of the self-diffusion coefficient D, the shear viscosity ηs, the bulk viscosity ηb and thermal conductivity λ of the hard sphere fluid, using molecular dynamics (MD) computer simulation. A newly developed hard sphere MD scheme was used to model the hard sphere fluid over a wide range up to the glass transition (∼0.57 packing fraction). System sizes of up to 32 000 hard spheres were considered. This set of transport coefficient data was combined with others taken from the literature to test a number of previously proposed analytical formulae for these quantities together with some new ones given here. Only the self-diffusion coefficient showed any substantial N dependence for N < 500 at equilibrium fluid densities (ε 0.494). D increased with N, especially at intermediate densities in the range ε ∼ 0.3–0.35. The expression for the packing fraction dependence of D proposed by Speedy, R. J., 1987, Molec. Phys., 62, 509 was shown to fit these data well for N ∼ 500 particle systems. We found that the packing fraction ε dependence of the two viscosities and thermal conductivity, generically denoted by X, were represented well by the simple formula X/X 0 = 1/[1 − (ε/ε1)]m within the equilibrium fluid range 0 > ε > 0.493. This formula has two disposable parameters, ε and m, and X 0 is the value of the property X in the limit of zero density. This expression has the same form as the Krieger-Dougherty formula (Kreiger, I. M., 1972, Adv. Colloid. Interface Sci., 3, 111) which is used widely in the colloid literature to represent the packing fraction dependence of the Newtonian shear viscosity of monodisperse colloidal near-hard spheres. Of course, in the present case, X o was the dilute gas transport coefficient of the pure liquid rather than the solvent viscosity. It was not possible to fit the transport coefficient normalized by their Enskog values with such a simple expression because these ratios are typically of order unity until quite high packing fractions and then diverge rapidly at higher values over a relatively narrow density range. At the maximum equilibrium fluid packing fraction ε = 0.494 for both the hard sphere fluid and the corresponding colloidal case a very similar value was found for ηs/ηo −30–40, suggesting that the ‘crowding’ effects and their consequences for the dynamics in this region of the phase diagram in the two types of liquid have much in common. For the hard sphere by MD, Do/D ∼ 11 at the same packing fraction, possibly indicating the contribution from ‘hydrodynamic enhancement’ of this transport coefficient, which is largely absent for the shear viscosity. Interestingly the comparable ratio for hard sphere colloids is the same.

Simulation and Data-Driven Modeling of the Transport Properties of the Mie Fluid

A Novel Investigation of a Micropolar Fluid Characterized by Nonlinear Constitutive Diffusion Model in Boundary Layer Flow and Heat Transfer.

A Comparative Study of TIP4P-2005, SPC/E, SPC, and TIP3P-Ew Models for Predicting Water Transport Coefficients Using EMD and NEMD Simulations

A New Single Equation of State to describe the Dynamic Viscosity and Self-Diffusion Coefficient for all Fluid Phases of Water from 200 K to 1800 K based on a New Original Microscopic Model

Shear stress relaxation and diffusion in simple liquids by molecular dynamics simulations: Analytic expressions and paths to viscosity

Molecular Dynamics Simulation of Self-Diffusion Coefficient and Its Relation with Temperature Using Simple Lennard-Jones Potential

Transport Properties of Anisotropic Polar Fluids:2. Dipolar Interaction

Viscosity and self-diffusion of supercooled and stretched water from molecular dynamics simulations

The shear viscosities of common water models by non-equilibrium molecular dynamics simulations

Entropy scaling framework for transport properties using molecular-based equations of state

Understanding coupled mass-heat transport in fluids by approach-to-equilibrium molecular dynamics

Simulation of high Schmidt number fluids with dissipative particle dynamics: Parameter identification and robust viscosity evaluation

Some universal trends of the Mie(n,m) fluid thermodynamics

Transport Properties of Anisotropic Polar Fluids: 1. Quadrupolar Interaction

Molecular kinetic modelling of non-equilibrium transport of confined van der Waals fluids

Transport coefficients of hard sphere fluids

Molecular dynamics study of transport coefficients for hard‐chain fluids

Molecular diffusion of mass and energy predicted by ab initio potential energy surfaces for air components at high temperatures

Self-diffusion coefficient and viscosity of methane and carbon dioxide via molecular dynamics simulations based on new ab initio-derived force fields

Diffusion of methane in a mica slit pore: Molecular dynamics simulations and correlation models

Transport properties of carbon dioxide and methane from molecular dynamics simulations