Modeling and simulation of three-dimensional planar contraction flow of viscoelastic fluids with PTT, Giesekus and FENE-P constitutive models

Abstract The planar contraction flow is a representative benchmark problem that is often used as a stringent test for the robustness of the numerical method to predict real flow patterns of viscoelastic fluids. A numerical methodology based on the penalty finite element method with a decoupled algorithm is presented in the study to simulate three-dimensional flow of viscoelastic fluids through planar contraction. The viscoelastic rheological responses are described by using three kind differential constitutive models including the Phan-Thien–Tanner (PTT) model, the Giesekus model and the finite extensible nonlinear elastic dumbbell with a peterlin closure approximation (FENE-P) model. The discrete elastic viscous split stress (DEVSS) formulation in cooperating with the inconsistent streamline upwind (SU) scheme is employed to improve the computational stability. The simulation results of flow velocity and stress on different sections of the flow field are compared with Quinzani’s experimental results that detected by laser-doppler velocimetry and flow-induced birefringence technologies. It is found that the simulation results predicted with three kind differential constitutive models agree well with the experimental results. The numerical methodology proposed in the study can be used successfully to predict complex flow patterns of viscoelastic fluids.

[1]  N. Phan-Thien,et al.  Three dimensional numerical simulations of viscoelastic flows through planar contractions , 1998 .

[2]  K. Funatsu,et al.  Numerical simulation of three‐dimensional viscoelastic flow within dies , 1999 .

[3]  David O. Olagunju Local similarity solutions for boundary layer flow of a FENE-P fluid , 2006, Appl. Math. Comput..

[4]  Constantin Fetecau,et al.  Flow of a generalized Oldroyd-B fluid due to a constantly accelerating plate , 2008, Appl. Math. Comput..

[5]  Youbing Yin,et al.  Oscillating flow of a viscoelastic fluid in a pipe with the fractional Maxwell model , 2006, Appl. Math. Comput..

[6]  Antony N. Beris,et al.  Pseudospectral simulation of turbulent viscoelastic channel flow , 1999 .

[7]  Michel Deville,et al.  Lattice Boltzmann method for the simulation of viscoelastic fluid flows , 2010 .

[8]  R. Armstrong,et al.  Birefringence and laser-Doppler velocimetry (LDV) studies of viscoelastic flow through a planar contraction , 1994 .

[9]  J. Ouyang,et al.  Numerical Simulation of the Non-isothermal Viscoelastic Flow Past a Confined Cylinder , 2010 .

[10]  M. R. Apelian,et al.  Numerically stable finite element techniques for viscoelastic calculations in smooth and singular geometries , 1988 .

[11]  N. Phan-Thien,et al.  An adaptive viscoelastic stress splitting scheme and its applications: AVSS/SI and AVSS/SUPG , 1996 .

[12]  Argimiro Resende Secchi,et al.  Viscoelastic flow analysis using the software OpenFOAM and differential constitutive equations , 2010 .

[13]  J. Guillet,et al.  Three-dimensional numerical simulation of viscoelastic contraction flows using the Pom-Pom differential constitutive model , 2005 .

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

[15]  H. Demir,et al.  Numerical modelling of viscoelastic cavity driven flow using finite difference simulations , 2005, Appl. Math. Comput..

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

[17]  Shangyou Zhang,et al.  Axisymmetric finite element solution of non-isothermal parallel-plate flow , 2005, Appl. Math. Comput..

[18]  Jason S. Howell,et al.  A two-parameter defect-correction method for computation of steady-state viscoelastic fluid flow , 2008, Appl. Math. Comput..

[19]  Leo G. Rebholz,et al.  Enforcing energy, helicity and strong mass conservation in finite element computations for incompressible Navier-Stokes simulations , 2011, Appl. Math. Comput..

[20]  F. Baaijens,et al.  Numerical simulations of the planar contraction flow for a polyethylene melt using the XPP model , 2004 .

[21]  Robert C. Armstrong,et al.  Finite element method for viscoelastic flows based on the discrete adaptive viscoelastic stress splitting and the discontinuous Galerkin method : DAVSS-G/DG , 1999 .

[22]  Murilo F. Tomé,et al.  Numerical simulation of viscoelastic flows using integral constitutive equations: A finite difference approach , 2008, J. Comput. Phys..

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

[24]  Xiaowei Gao A promising boundary element formulation for three‐dimensional viscous flow , 2005 .

[25]  S. Ramazani,et al.  Predictions of some internal microstructural models for polymer melts and solutions in shear and elongational flows , 2004 .

[26]  H. Giesekus A simple constitutive equation for polymer fluids based on the concept of deformation-dependent tensorial mobility , 1982 .

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

[28]  J. Azaiez,et al.  Numerical simulation of viscoelastic flows through a planar contraction , 1996 .

[29]  R. Armstrong,et al.  Finite element methdos for calculation of steady, viscoelastic flow using constitutive equations with a Newtonian viscosity , 1990 .

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

[31]  Guoqun Zhao,et al.  Measurement and simulation of low‐density polyethylene extrudate swell through a circular die , 2009 .