Abstract:This work presents a comprehensive numerical implementation of a viscoelastic thixotropic model known as the modified-Bautista–Manero (MBM) model (a model widely used to characterize the rheological behavior of wormlike micellar solutions). This implementation is integrated into the OpenFOAM computational fluid dynamics software, specifically using the rheoTool library. In contrast to previous studies, our approach includes simulations of inertial flows and we account for convective terms within all the governing equations. Moreover, we use two numerical techniques that have been proven to be efficient methods in computational simulations of viscoelastic fluids through a planar contraction: (1) the high-resolution scheme called convergent and universally bounded interpolation scheme for the treatment of advection to discretize the convective terms of the momentum, viscoelastic, and thixotropic equations; and (2) the logarithm-conformation reformulation to tackle the high-Weissenberg number problem. This expansion of the model's scope results in a more realistic representation of the fluid's behavior and, consequently, more accurate computational results. To validate the implementation of the MBM model, we performed numerical simulations in two different industrial-relevant geometries: the planar 4:1 contraction and the 4:1:4 contraction–expansion configurations. Two fluids with different structural responses are studied: one with quick structural recovery and the other with mildly slow recovery. The numerical method developed in this work allowed us to run stable simulations of these fluids up to Reynolds and Deborah number values of Re = 1.11 and De = 14.0, respectively, and the influence of inertial, viscoelastic, and thixotropic effects on various flow field variables is investigated. These variables include velocity, viscosity, normal stresses, and corner vortex size. Transient and steady solutions of corner vortexes across a range of Deborah and Reynolds numbers are also presented. In addition, the results obtained in this work are also directly compared with numerical simulations obtained using the non-thixotropic rubber network-based exponential Phan-Thien–Tanner model. For the planar 4:1 contraction simulations, the analysis of the results is based on three dimensionless numbers (the viscoelastic Mach, the elasticity, and the thixoelastic numbers), which allows to conclude that vortex enhancement is mainly seen when high elasticity is coupled with quick structural reformation and very low inertial effects. We further report from the planar 4:1:4 contraction–expansion results that an increase in inertia leads both to vortex inhibition in the upstream channel and to slight vortex enhancement in the downstream channel. Finally, we demonstrate the strong influence of fluid convection on fluidity profiles and on the size of the upstream/downstream corner vortex.

A model for micro-front dynamics using a φ 4 equation

Numerical Analysis of the Micro-Couette Flow Using a Non-Newton–fourier Model with Enhanced Wall Boundary Conditions

A Novel Investigation of a Micropolar Fluid Characterized by Nonlinear Constitutive Diffusion Model in Boundary Layer Flow and Heat Transfer.

Modelling Slip Flow in Micro/nano Gaps with Moving Boundary

Understanding the dynamics of fluid–structure interaction with an Air Deflected Microfluidic Chip (ADMC)

A deterministic-particle-based scheme for micro-macro viscoelastic flows

Strong pinning of propagation fronts in adverse flow

Generalized Maxwell Velocity Slip Boundary Model

A well-balanced positivity preserving two-dimensional shallow flow model with wetting and drying fronts over irregular topography

Front pinning in capillary filling of chemically coated channels

Molecular kinetic modelling of nanoscale slip flow using a continuum approach

Modeling and Simulation of Dynamics of Three-Component Flows on Solid Surface

Bridging the scales in capillary rise dynamics with complexity-reduced models

Numerical analysis of the entry flow for continuous shear thickening fluids in circular microtubes

A Squeeze Flow Problem with A Navier Slip Condition

A fluid particle model

Continuum simulation for regularized non-local μ(I) model of dense granular flows

Numerical simulation of a thixotropic-viscoelastic model in contraction geometries

A Semi-Empirical Force Balance-Based Model to Capture Sessile Droplet Spread on Smooth Surfaces: A Moving Front Kinetic Monte Carlo Study

A mesoscopic model for microscale hydrodynamics and interfacial phenomena: Slip, films, and contact angle hysteresis

A modeling approach for capillary-driven flow of a wetting fluid in a rectangular open microchannel of arbitrary axisymmetric width profile: Application to heat pipes