Convoluted models and high-Weissenberg predictions for micellar thixotropic fluids in contraction–expansion flows

Abstract This study is concerned with finite element/volume modelling of contraction–expansion axisymmetric pipe flows for thixotropic and non-thixotropic viscoelastic models. To obtain solutions at high Weissenberg numbers (Wi) under a general differential form W i τ p ∇ = 2 ( 1 − β ) D − f τ p , both thixotropic Bautista–Manero micellar and non-thixotropic EPTT f-functionals have been investigated. Here, three key modifications have been implemented: first, that of convoluting EPTT and micellar Bautista–Manero f-functionals, either in a multiplicative (Conv*) or additive (Conv+) form; second, by adopting f-functionals in absolute form (ABS-f-correction); and third, by imposing pure uniaxial-extension velocity-gradient components at the pure-stretch flow-centreline (VGR-correction). With this combination of strategies, highly non-linear solutions have been obtained to impressively high Wi [=O(5000+)]. This capability permits analysis of industrial applications, typically displaying non-linear features such as thixotropy, yield stress and shear banding. The scope of applications covers enhanced oil-recovery, industrial processing of plastics and foods, as well as in biological and microfluidic flows. The impact of rheological properties across convoluted models (moderate-hardening, shear-thinning) has been observed through steady-state solutions and their excess pressure-drop (epd) production, stress, f-functional field structure, and vortex dynamics. Three phases of vortex-behaviour have been observed with rise in elasticity, along with upstream–downstream Moffatt vortices and plateauing epd-behaviour at high-Wi levels. Moreover, enhancement of positive-definiteness in stress has improved high-Wi solution attenuation.

[1]  Johnson–Segalman model with a diffusion term in cylindrical Couette flow , 1998, cond-mat/9805286.

[2]  Gareth H. McKinley,et al.  A network scission model for wormlike micellar solutions: I. Model formulation and viscometric flow predictions , 2007 .

[3]  J. Rothstein.  Strong Flows of Viscoelastic Wormlike Micelle Solutions , 2009 .

[4]  P. Olmsted.  Perspectives on shear banding in complex fluids , 2008 .

[5]  S. Manneville.  Recent experimental probes of shear banding , 2008, 0903.5389.

[6]  M. F. Webster,et al.  Transient viscoelastic flows in planar contractions , 2004 .

[7]  O. Manero,et al.  On the shear banding flow of elongated micellar solutions , 2000 .

[8]  Jt Johan Padding,et al.  Constitutive equations for extensional flow of wormlike micelles: stability analysis of the Bautista–Manero model , 2005 .

[9]  R. Tanner,et al.  A new constitutive equation derived from network theory , 1977 .

[10]  M. Webster,et al.  A second-order hybrid finite-element/volume method for viscoelastic flows 1 Dedicated to Professor M , 1998 .

[11]  S. Manneville,et al.  Shear-banding in surfactant wormlike micelles: elastic instabilities and wall slip , 2012 .

[12]  Jiang Yang.  Viscoelastic wormlike micelles and their applications , 2002 .

[13]  J. Berret.  Rheology of Wormlike Micelles: Equilibrium Properties and Shear Banding Transitions , 2004, cond-mat/0406681.

[14]  M.A. Hulsen.  A sufficient condition for a positive definite configuration tensor in differential models , 1990 .

[15]  G. McKinley,et al.  Wormlike Micellar Solutions: III. VCM Model Predictions in Steady and Transient Shearing Flows , 2014 .

[16]  M. F. Webster,et al.  Recovery and stress-splitting schemes for viscoelastic flows , 1998 .

[17]  M. Aboubacar,et al.  Time‐dependent algorithms for viscoelastic flow: Finite element/volume schemes , 2005 .

[18]  M. F. Webster,et al.  A new constitutive model for worm-like micellar systems – Numerical simulation of confined contraction–expansion flows , 2014 .

[19]  M. F. Webster,et al.  Excess pressure drop and drag calculations for strain-hardening fluids with mild shear-thinning: Con , 2011 .

[20]  M. Webster,et al.  Alternative subcell discretisations for viscoelastic flow: Velocity-gradient approximation , 2008 .

[21]  O. Manero,et al.  Dynamics of worm-like micelles: the Cox–Merz rule , 2002 .

[22]  H. R. Tamaddon-Jahromi,et al.  High-Weissenberg predictions for micellar fluids in contraction-expansion flows , 2015 .

[23]  N. Wagner,et al.  Large amplitude oscillatory shear (LAOS) measurements to obtain constitutive equation model parameters: Giesekus model of banding and nonbanding wormlike micelles , 2012 .

[24]  O. Manero,et al.  On the yield stress of complex materials , 2013, Korea-Australia Rheology Journal.

[25]  H. R. Tamaddon-Jahromi,et al.  Numerical modelling of thixotropic and viscoelastoplastic materials in complex flows , 2014, Rheologica Acta.

[26]  M. S. Chandio,et al.  ON CONSISTENCY OF CELL-VERTEX FINITE-VOLUME FORMULATIONS FOR VISCOELASTIC FLOW , 2000 .

[27]  J. P. Aguayo,et al.  Numerical investigation of transient contraction flows for worm-like micellar systems using Bautista , 2011 .

[28]  O. Manero,et al.  Understanding thixotropic and antithixotropic behavior of viscoelastic micellar solutions and liquid crystalline dispersions. I. The model , 1999 .

[29]  J. Rothstein,et al.  Flow of wormlike micelle solutions past a confined circular cylinder , 2010 .

[30]  Gareth H. McKinley,et al.  Wormlike micellar solutions: II. Comparison between experimental data and scission model predictions , 2010 .