Predicting the shear viscosities of liquids using LAMMPS
The ability to predict the liquids properties is of great industrial interest, since they are components of oils and fuels. In this work, we apply classical molecular dynamics methods using LAMMPS to calculate the equation of state and self-diffusion and viscosity coefficients of various compounds [1-8]. The interatomic force fields are compared by their ability to predict the experimental data. The molecular dynamics methods are used to predict the shear viscosities in the pressure range from 0.1 MPa to 1 GPa using both equilibrium and non-equilibrium methods. Blindly predicted viscosity values for model lubricants match the experimental data measured later. We demonstrate the universal scaling between the viscosity and diffusivity at pressures up to 800 MPa.
1. V. Pisarev and N. Kondratyuk. Fluid Phase Equilib. 2019. 501. 112273. 2. N.D. Kondratyuk and V.V. Pisarev. Fluid Phase Equilib. 2019. 498. 151-159. 3. N. Kondratyuk. J. Chem. Phys. 2019. 151(7). 074502. 4. N. Kondratyuk, D. Lenev and V. Pisarev. J. Chem. Phys. 2020. 152. 191104. 5. N.D. Kondratyuk, V.V. Pisarev, J.P. Ewen. J. Chem. Phys. 2020. 153(15). 154502. 6. I. Bakulin, N. Kondratyuk, A. Lankin, G. Norman. J. Chem. Phys. 2021. 155, 154501. 7. V. Deshchenya, N. Kondratyuk, A. Lankin, G. Norman. J. Mol. Liq. 2022. 367. 120456. 8. B. Nikitiuk, D. Salikova, N. Kondratyuk, V. Pisarev. J. Mol. Liq. 2022. 368. 120714.