Computer Simulation of Steady Shear Flow of Liquid Crystals with the Gay-Berne Potential

Abstract

Computer simulations of steady shear flow of nematic phases are performed using a non-equilibrium molecular dynamics. The SLLOD-type algorithm is applied to the shear flows of Gay-Berne fluids which are nematic in the equilibrium states. Two Gay-Berne potentials with μ = 1, v = 2 and μ = 2, v = 1 are used and both molecular aspect ratios are set equal to 3. In the initial stage of shearing various patterns of the transient behavior are found in the order parameter and orientation angle, according to the initial conditions and shear rate. In the steady state the director is inclined to the flow at a constant angle regardless of the shear rate and the system becomes more highly oriented state. The fluids exhibit shear-thinning in viscosity and non-zero first and second normal stress differences.