An improved suspension balance model applied to shear-induced phase segregation

• Published in: International journal of multiphase flow Vol. 184; p. 105120
• Main Authors: Schlatter, Lauren, da Silva Ferreira, Gabriel Gonçalves, da Cunha Lage, Paulo Laranjeira
• Format: Journal Article
• Language: English
• Published: Elsevier Ltd 01.03.2025
• Subjects: Non-Newtonian fluids; OpenFOAM; Shear-induced migration; Suspension balance model; Suspension rheology; OpenFOAM; Shear-induced migration; Non-Newtonian fluids; Suspension balance model; Suspension rheology
• ISSN: 0301-9322
• DOI: 10.1016/j.ijmultiphaseflow.2024.105120

When suspensions are subject to non-homogeneous shear, particles migrate towards lower-shear rate regions of the flow due to the anisotropy on the particle phase’s normal stress. This phenomenon leads to phase segregation, complicating the interpretation of rheological experimental data for characterising those suspensions. This work presents an improved version of the well-known Suspension Balance Model (SBM), featuring a frame-independent formulation of the particles’ normal stress with an improved momentum interpolation scheme that prevents numerical oscillations. The particle’s stress model also includes a local formulation for the microscopically generated extra stress, ensuring grid convergence. Our model, implemented in OpenFOAM-v7®, successfully simulated various shear-dominated flows. The simulated Couette rheometric data showed that the characterisation of suspensions from such data could not capture their non-Newtonian behaviour. Besides, the obtained rheological model for the suspensions depended on the geometry of the rheometer, being unable to predict the flow of the same suspension in different conditions and geometries. [Display omitted] •We developed a frame-independent Suspension Balance Model solver using OpenFOAM.•A modified stress tensor model allows the presence of velocity plateaus.•An improved momentum interpolation avoids oscillations in the phase fraction fields.•A cyclic boundary condition simulated fully developed flows up to 500 times faster.•The suspension’s apparent viscosity strongly depends on rheometer flow geometry.