A two-way coupling scheme to model the effects of particle rotation on the rheological properties of a semidilute suspension. (15th September 2018)

• Record Type: Journal Article
• Authors: Fukui, Tomohiro; Kawaguchi, Misa; Morinishi, Koji
• Is Part Of: Computers & fluids. Volume 173(2018)
• Page Start: 6
• Page End: 16
• Publication Date: 2018-09-15
• Subjects: Effective viscosity -- Spatial variation -- Rheology -- Thixotropy -- Einstein's viscosity formula -- Two-way coupling; Fluid dynamics -- Data processing -- Periodicals; 532.050285
• DOI: 10.1016/j.compfluid.2018.04.038
• Languages: English
• ISSNs: 0045-7930

Abstract

Highlights: The effective viscosity is overestimated when the particle's rotation is neglected. The effective viscosity varies in time and space with the particles' position. The viscosity variation becomes more significant with increasing area fraction. The coefficient β of the term quadratic in ϕ is found to be 0.773 for 2D. The value of β may increase due to fluid resistance against the particle's rotation. Abstract: According to Einstein's viscosity formula, the effective viscosity of a suspension is related to the volume fraction ϕ of suspended particles. Higher-order terms, which include contributions from hydrodynamic interactions, become more important when the concentration of the suspension reaches around 25%. The coefficient β of the term quadratic in ϕ, however, has not been incontrovertibly validated. In order to investigate the effects on the macroscopic rheology of a suspension due to the rotational motions of the microscopic suspended particles—which are usually neglected in analytical and numerical studies—we have performed simulations of a pressure-driven flow using a two-way coupling method (defined in text). We used the regularized lattice Boltzmann equation for the fluid part of the calculation, and we employed Newton's second law of motion and the equation of angular motion for the translational and rotational motions of the rigid particles. As a result, we found that the particles rotate to achieve kinetic balance with the surrounding hydrodynamic forces, which results in a decrease in fluid resistance. It is therefore important to take into account the particles' rotational motions and the accompanying changes in the flow field in order to evaluate the effective viscosity of a suspension. The spatial variation of the relative viscosity determined from the wall shear stress becomes more significant with increasing area fraction. We found that the value of the coefficient β may increase due to hydrodynamic resistance against the rotational motions of the particles.