Numerical Simulation of a Thixotropic-Viscoelastic Model in Contraction Geometries

In this work, we incorporate a thixotropic-viscoelastic model into the widely used Computational Fluid Dynamics (CFD) software OpenFOAM, along with the rheoTool library. The model we implement is known as the Modified-Bautista-Manero (MBM), and effectively describes the rheological behavior of worm-like micellar solutions in extensional flows. We provide a detailed explanation of the numerical implementation of the model, specifically using the log-conformation tensor approach. Unlike previous works focused on this kind of fluids, we simulate inertial flows while considering convective terms in the governing equations, thus obtaining a more realistic behavior on the calculated results. The MBM model implementation is validated through numerical simulations on two different industrial-relevant geometries: the planar 4:1 contraction and the 4:1:4 contraction-expansion configurations. Furthermore, we investigate the influence of inertial, viscoelastic, and thixotropic effects on various flow field variables. These variables include velocity, viscosity, normal stresses, and corner vortex size. Our analysis encompasses both transient and steady solutions of corner vortexes across a range of Deborah and Reynolds numbers. Our results are also directly compared with simulations obtained using the non-thixotropic rubber network-based exponential Phan-Thien-Tanner (EPTT) model. From our planar 4:1 contraction results, we found that vortex-enhancement is seen when high elasticity is coupled with quick structural reformation and very low inertial effects. From our planar 4:1:4 contraction-expansion simulations, we show that an increase in inertia leads both to vortex-inhibition in the upstream channel and slight vortex-enhancement in the downstream channel. Lastly, we show the strong effect of the convection of fluidity into the fluidity profiles and into the upstream/downstream corner vortex sizes.