Numerical aspects in modeling high Deborah number flow and elastic instability

Investigating highly nonlinear viscoelastic flow in 2D domain, we explore problem as well as property possibly inherent in the streamline upwinding technique (SUPG) and then present various results of elastic instability. The mathematically stable Leonov model written in tensor-logarithmic formulation is employed in the framework of finite element method for spatial discretization of several representative problem domains. For enhancement of computation speed, decoupled integration scheme is applied for shear thinning and Boger-type fluids. From the analysis of 4:1 contraction flow at low and moderate values of the Deborah number (De) the solution with SUPG method does not show noticeable difference from the one by the computation without upwinding. On the other hand, in the flow regime of high De, especially in the state of elastic instability the SUPG significantly distorts the flow field and the result differs considerably from the solution acquired straightforwardly. When the strength of elastic flow and thus the nonlinearity further increase, the computational scheme with upwinding fails to converge and evolutionary solution does not become available any more. All this result suggests that extreme care has to be taken on occasions where upwinding is applied, and one has to first of all prove validity of this algorithm in the case of high nonlinearity. On the contrary, the straightforward computation with no upwinding can efficiently model representative phenomena of elastic instability in such benchmark problems as 4:1 contraction flow, flow over a circular cylinder and flow over asymmetric array of cylinders. Asymmetry of the flow field occurring in the symmetric domain, enhanced spatial and temporal fluctuation of dynamic variables and flow effects caused by extension hardening are properly described in this study.

[1]  Fernando T. Pinho,et al.  Dynamics of high-Deborah-number entry flows: a numerical study , 2011 .

[2]  Youngdon Kwon,et al.  Finite element analysis of planar 4:1 contraction flow with the tensor-logarithmic formulation of differential constitutive equations , 2004 .

[3]  Youngdon Kwon,et al.  Mathematical characteristics of the pom-pom model , 2002 .

[4]  Norio Takeuchi,et al.  Mixed finite element method for analysis of viscoelastic fluid flow , 1977 .

[5]  A. I. Leonov,et al.  On instabilities of single-integral constitutive equations for viscoelastic liquids , 1994 .

[6]  A. I. Leonov.  On a class of constitutive equations for viscoelastic liquids , 1987 .

[7]  Raanan Fattal,et al.  Flow of viscoelastic fluids past a cylinder at high Weissenberg number : stabilized simulations using matrix logarithms , 2005 .

[8]  Daniel D. Joseph,et al.  Fluid Dynamics Of Viscoelastic Liquids , 1990 .

[9]  Mahesh Gupta,et al.  VISCOELASTIC MODELLING OF ENTRANCE FLOW USING MULTIMODE LEONOV MODEL , 1997 .

[10]  Raanan Fattal,et al.  Constitutive laws for the matrix-logarithm of the conformation tensor , 2004 .

[11]  A. Beris,et al.  Finite element calculation of viscoelastic flow in a journal bearing: I. small eccentricities , 1984 .

[12]  Youngdon Kwon.  Numerical description of elastic flow instability and its dependence on liquid viscoelasticity in planar contraction , 2012 .

[13]  P. Raviart,et al.  On a Finite Element Method for Solving the Neutron Transport Equation , 1974 .

[14]  Eric S. G. Shaqfeh,et al.  Purely elastic instabilities in viscometric flows , 1996 .

[15]  P. Townsend,et al.  Expansion flows of non-Newtonian liquids , 1994 .

[16]  Marcel Crochet,et al.  Flows of polymer solutions through contractions .1. Flows of polyacrylamide solutions through planar contractions , 1996 .

[17]  Gareth H. McKinley,et al.  Investigating the stability of viscoelastic stagnation flows in T-shaped microchannels , 2009 .

[18]  R. Keunings.  MICRO-MACRO METHODS FOR THE MULTISCALE SIMULATION OF VISCOELASTIC FLOW USING MOLECULAR MODELS OF KINETIC THEORY , 2004 .

[19]  R. Bird.  Dynamics of Polymeric Liquids , 1977 .

[20]  F. Baaijens,et al.  Numerical analysis of start-up planar and axisymmetric contraction flows using multi-mode differential constitutive models , 1993 .

[21]  R. Poole,et al.  Purely elastic flow asymmetries. , 2007, Physical review letters.

[22]  Arkadii I. Leonov,et al.  On the rheological modeling of viscoelastic polymer liquids with stable constitutive equations , 1995 .

[23]  R. Courant,et al.  On the Partial Difference Equations, of Mathematical Physics , 2015 .

[24]  M. Fortin,et al.  A new approach for the FEM simulation of viscoelastic flows , 1989 .

[25]  R. Keunings,et al.  Die Swell of a Maxwell Fluid - Numerical Prediction , 1980 .

[26]  Mehmet Sahin.  Parallel large-scale numerical simulations of purely-elastic instabilities behind a confined circular cylinder in a rectangular channel , 2013 .

[27]  Ma Martien Hulsen,et al.  Simulation of the Doi-Edwards model in complex flow , 1999 .

[28]  Kwang Sun Park,et al.  Decoupled algorithm for transient viscoelastic flow modeling , 2012, Korea-Australia Rheology Journal.

[29]  A. I. Leonov.  Nonequilibrium thermodynamics and rheology of viscoelastic polymer media , 1976 .

[30]  R. Larson.  The Structure and Rheology of Complex Fluids , 1998 .

[31]  R. Larson.  Constitutive equations for polymer melts and solutions , 1988 .

[32]  Ma Martien Hulsen,et al.  Decoupled second-order transient schemes for the flow of viscoelastic fluids without a viscous solvent contribution , 2010 .

[33]  K. Walters,et al.  Long range memory effects in flows involving abrupt changes in geometry: Part 2: the expansion/contraction/expansion problem , 1977 .

[34]  R. Courant,et al.  Über die partiellen Differenzengleichungen der mathematischen Physik , 1928 .

[35]  T. Hughes,et al.  Streamline upwind/Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations , 1990 .

[36]  K. Walters,et al.  Long-range memory effects in flows involving abrupt changes in geometry: Part I: flows associated with I-shaped and T-shaped geometries , 1977 .

[37]  J. S. Lee,et al.  Comparison of numerical simulations and birefringence measurements in viscoelastic flow between eccentric rotating cylinders , 1992 .

[38]  R. Keunings.  On the high weissenberg number problem , 1986 .

[39]  L. Collins,et al.  Numerical approach to simulating turbulent flow of a viscoelastic polymer solution , 2003 .

[40]  Andreas G. Boudouvis,et al.  On Hadamard stability and dissipative stability of the molecular stress function model of non-linear viscoelasticity , 2009 .

[41]  F. Baaijens.  Mixed finite element methods for viscoelastic flow analysis : a review , 1998 .

[42]  Robert C. Armstrong,et al.  Approximation error in finite element calculation of viscoelastic fluid flows , 1982 .

[43]  F. Brezzi.  On the existence, uniqueness and approximation of saddle-point problems arising from lagrangian multipliers , 1974 .

[44]  M. Fortin,et al.  A new mixed finite element method for computing viscoelastic flows , 1995 .

[45]  K. Walters,et al.  The behaviour of polymer solutions in extension-dominated flows, with applications to Enhanced Oil Recovery , 1989 .

[46]  Youngdon Kwon,et al.  On instability of the Doi–Edwards model in simple flows , 1999 .

[47]  R. Keunings.  FINITE ELEMENT METHODS FOR INTEGRAL VISCOELASTIC FLUIDS , 2003 .

[48]  A. I. Leonov,et al.  Stability constraints in the formulation of viscoelastic constitutive equations , 1995 .

[49]  M. Crochet,et al.  A new mixed finite element for calculating viscoelastic flow , 1987 .

[50]  I. Babuska.  The finite element method with Lagrangian multipliers , 1973 .

[51]  R. A. Silverman,et al.  The Mathematical Theory of Viscous Incompressible Flow , 1972 .