Numerical Simulation of Rheological Models for Complex Fluids Using Hierarchical Grids

Abstract

In this work, we implement models that are able to describe complex rheological behaviour (such as shear-banding and elastoviscoplasticity) in the HiGTree/HiGFlow system, which is a recently developed Computational Fluid Dynamics (CFD) software that can simulate Newtonian, Generalised-Newtonian and viscoelastic flows using finite differences in hierarchical grids. The system uses a moving least squares (MLS) meshless interpolation technique, allowing for more complex mesh configurations while still keeping the overall order of accuracy. The selected models are the Vasquez-Cook-McKinley (VCM) model for shear-banding micellar solutions and the Saramito model for viscoelastic fluids with yield stress. Development of solvers and numerical simulations of inertial flows of these models in 2D channels and planar-contraction 4:1 are carried out in the HiGTree/HiGFlow system. Our results are compared with those predicted by two other methodologies: the OpenFOAM-based software RheoTool that uses the Finite-Volume-Method and an in-house code that uses the Vorticity-Velocity-Formulation (VVF). We found an excellent agreement between the numerical results obtained by these three different methods. A mesh convergence analysis using uniform and refined meshes is also carried out, where we show that great convergence results in tree-based grids are obtained thanks to the finite difference method and the meshless interpolation scheme used by the HiGFlow software. More importantly, we show that our methodology implemented in the HiGTreee/HiGFlow system can successfully reproduce rheological behaviour of high interest by the rheology community, such as non-monotonic flow curves of micellar solutions and plug-flow velocity profiles of yield-stress viscoelastic fluids.

Fluctuating Flexoelectric Membranes in Asymmetric Viscoelastic Media: Power Spectrum through Mechanical Network and Transfer Function Models

“…The computed norm of the complex wave vector as a function of dimensionless frequency is shown in Figure 7. It clearly shows the three-mode response: flexo-viscous with linear frequency scaling, intermediate flexo-viscoelastic with square-root frequency scaling, and flexo-bulk elastic linear scaling [29,30]. The role of the elastic number is a vertical shift [31,32].…”

“…The dimensionless elastic number E (Equation ( 26)) is the ratio of two fundamental numbers: (i) the Weissenberg (We) and (ii) the Reynolds number (Re) [28]. The elastic number E is given by [29][30][31]:…”

SPH Simulation of Molten Metal Flow Modeling Lava Flow Phenomena with Solidification

“…However, its The complex viscosity, the solidification of lava, and the deformation of the solidliquid interfaces make numerical simulations difficult. Some previous studies [14][15][16][17] proposed methods to simulate the behavior of non-Newtonian fluids. Some previous studies [5][6][7] took a simplified approach to the solidification of lava and the deformation of the solid-liquid interfaces.…”