Yamamoto Takehiro, Suga Takanori, Mori Noriyasu
Department of Mechanophysics Engineering, Osaka University, 2-1, Yamadaoka, Suita, Osaka 565-0871, Japan.
Phys Rev E Stat Nonlin Soft Matter Phys. 2005 Aug;72(2 Pt 1):021509. doi: 10.1103/PhysRevE.72.021509. Epub 2005 Aug 24.
Brownian dynamics (BD) simulations were carried out for suspensions of oblate spheroid particles interacting via the Gay-Berne (GB) potential. The oblate spheroid particles were applied as a model of disc-like particles and the system of suspension of the particles was considered. Numerically analyzed were both the change in phase with the number density of the particles at equilibrium state and the behavior of the particles in simple shear flows. The system changed from isotropic phase to nematic one with increasing the particle concentration. In the simulation of shear flows, the shear was imposed upon the systems in nematic phase at equilibrium. The systems exhibited various motions of the director depending on the shear rate, e.g. the continuous rotation of director at low shear rates, the wagging at moderate shear rates, and the flow aligning at high shear rates. Temporal change in inner structure of suspensions was also analyzed and collapse of initial particle configurations due to shear was found. Moreover, rheological properties of the suspension were investigated. The numerical simulation predicted the shear-thinning in viscosity, negative first normal stress difference, and positive second normal stress difference, and these results qualitatively agreed with the predictions using a constitutive equation for discotic nematics. The present study proved that the BD simulation using spheroid particles interacting via the GB potential is an effective approach for investigating the flow behavior and flow-induced structure of suspensions of disklike particles at a particulate level.
对通过盖伊 - 伯恩(GB)势相互作用的扁球形颗粒悬浮液进行了布朗动力学(BD)模拟。将扁球形颗粒用作盘状颗粒的模型,并考虑了颗粒悬浮液系统。对平衡状态下颗粒数密度引起的相变化以及颗粒在简单剪切流中的行为进行了数值分析。随着颗粒浓度的增加,系统从各向同性相转变为向列相。在剪切流模拟中,在平衡的向列相系统上施加剪切力。系统根据剪切速率表现出指向矢的各种运动,例如在低剪切速率下指向矢的连续旋转、在中等剪切速率下的摆动以及在高剪切速率下的流动排列。还分析了悬浮液内部结构的时间变化,发现由于剪切导致初始颗粒构型的崩塌。此外,研究了悬浮液的流变性质。数值模拟预测了粘度的剪切变稀、负的第一法向应力差和正的第二法向应力差,这些结果与使用盘状向列相本构方程的预测在定性上一致。本研究证明,使用通过GB势相互作用的球形颗粒进行的BD模拟是在颗粒水平上研究盘状颗粒悬浮液的流动行为和流动诱导结构的有效方法。