Simulation of lubricant rheology in thin film lubrication part 2: Simulation of couette flow

Abstract

In the second part of this two-part study, molecular dynamics simulations are performed for a fluid of spherical molecule in Couette flow. The simulation uses the same system as that in Part 1, but with the top wall translating in the cursive Greek chi direction to generate a Couette flow. The shear equivalent viscosity ηes of a fluid is found to increase as film thickness decreases. The dependence is similar to that of the flow equivalent viscosity ηef obtained in Part 1, but the shear viscosity ηes exhibits a smaller value and slower increase rate. A further comparison between ηes and ηef shows that at the film thickness where the viscosity ηef diverges, the corresponding shear viscosity keeps a relatively small value, which is attributed to the larger shear rate applied in simulation of Couette flow. In a region where shear rate is low, the shear viscosity remains almost constant until a critical shear rate γ̇c is reached, then the 'shear-thinning' follows, i.e. the viscosity declines in a power-law, ηes∼γ̇-2/3, and the mean shear stress approaches a constant value - the 'limiting shear stress'. If the film becomes molecularly thin, the lubricant behaves like a viscoelastic material, indicated by the considerable value of shear stress existing at the zero shear rate. The mean velocity profile of molecular flow shows a linear distribution, but with inflexions on the profile near the wall-fluid interface. When a very high shear rate is applied, however, the flow in thin films seems to be divided into two parts, half stays almost at rest and half is rapidly sheared. Once the stress exceeds the limiting shear stress, a slip in velocity appears at the solid-fluid interface.