Unraveling the pressure-viscosity behavior and shear thinning in glycerol using atomic scale molecular dynamics simulations

In order to increase the usage and explore new applications of glycerol as a replacement for fossil-based lubricants its properties needs to be known at the fundamental level. In this study, the viscosity of pure glycerol at high pressures and strain rates has been investigated using of molecular dynamics (MD) simulations, utilizing both the Green-Kubo (GK) formalism and the SLLOD algorithm. Although the viscosity acquired by the GK method is in agreement with the corresponding experimental values at low pressure, a significant distinction was identified between the viscosity obtained by the GK method and the experimental values at higher pressures (P > 0.5 GPa). This results in a clear difference between the viscosity-pressure coefficient attained by the GK method and the corresponding experimental value. The SLLOD method using a non-equilibrium MD (NEMD) platform was exploited to take into account the simultaneous effects of strain rate and pressure on viscosity. As a result, the pressure-viscosity coefficient acquired by the SLLOD algorithm approaches the experimental value. By combining the experimental outputs for viscosity at low strain rates (gamma < 10(4 )s(-1)) with the SLLOD outputs at higher rates (gamma > 10(5 )s(-1)), the evolutions of glycerol viscosity with pressure and strain rate were ultimately achieved. Implementing this computational platform depicts the shear thinning process in pure glycerol in a wide range of pressures and strain rates.