The behavior of lubricants at operational conditions, such as at high pressures, is a topic of great industrial interest. In particular, viscosity and the viscosity-pressure relation are especially important for applications and their determination by computational simulations is very desirable. In this study we evaluate methods to compute these quantities based on fully atomistic molecular dynamics simulations which are computationally demanding but also have the potential to be most accurate. We used the 9,10dimethyloctadecane molecule, main component of PAO-2 base oil as the lubricant for our tests. The methods used for the viscosity simulations are the Green-Kubo equilibrium molecular dynamics (EMD-GK) and non-equilibrium molecular dynamics (NEMD), at pressures of up to 1.0 GPa and various temperatures (40-150 degrees Celsius). We present the theory behind these methods and investigate how the simulation parameters affect the results obtained, to ensure viscosity convergence with respect to the A c c e p t e d M a n u s c r i p t simulation intervals and all other parameters. We show that by using each method in its regime of applicability, we can achieve good agreement between simulated and measured values. NEMD simulations at high pressures captured zero shear viscosity successfully, while at 40 degrees Celsius EMD-GK is only applicable to pressures up to 0.3 GPa, where the viscosity is lower. In NEMD, longer and multiply repeated simulations improve the confidence interval of viscosity, which is essential at lower pressures. Another aspect of these methods is the choice of the utilized force field for the atomic interactions. This was investigated by selecting two different commonly used force fields.
The continuum theory of elasticity and/or homogeneously discretized finite element models have been commonly used to investigate and analyze subsurface stresses in Hertzian contacts. These approaches, however, do not effectively capture the influence of the random microstructure topology on subsurface stress distributions in Hertzian contacts. In this paper, a finite element model for analyzing subsurface stresses in an elastic half-space subjected to a general Hertzian contact load with explicit consideration of the material microstructure topology is presented. The random internal geometry of polycrystalline microstructures is modeled using a 3D Voronoi tessellation, where each Voronoi cell represents a distinct material grain. The grains are then meshed using finite elements, and an algorithm was developed to eliminate poorly shaped elements resulting from “near degeneracy” in the Voronoi tessellations. Hertzian point and line contacts loads are applied as distributed surface loads, and the model’s response is evaluated with commercial finite element software ABAQUS. Internal stress results obtained from the current model compare well with analytical solutions from theory of elasticity. The influence of the internal microstructure topology on the subsurface stresses is demonstrated by analyzing the model’s response to an over rolling element using a critical plane approach.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.