Numerical simulation of viscoelastic flows through a planar contraction

Abstract In this study, three nonlinear rheological models consisting of the Giesekus, the FENE-P, and the Phan-Thien-Tanner model are used to simulate the flow of a viscoelastic fluid through a planar 4:1 contraction. Both stress and velocity fields are examined at different sections of the flow and the predictions of the numerical simulations are compared with the experimental results of L.M. Quinzani, R.C. Armstrong adn R.A. Brown, J. Non-Newtonian Fluid Mech., 52 (1994) 1–36. Overall, the numerical simulations allow a description of the essential features of the flow, and reproduce much of the experimental results with good accuracy. Excellent qualitative agreement between the numerical results and the experimental observations is reported. However, the agreement remains semi-quantitative especially for the first normal stress difference around the entry section of the flow. This is to be expected in view that the simulations were limited to only one-mode models.

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

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

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

[4]  Robert C. Armstrong,et al.  Calculation of steady-state viscoelastic flow through axisymmetric contractions with the EEME formulation , 1992 .

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

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

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

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

[9]  Michel Fortin,et al.  A preconditioned generalized minimal residual algorithm for the numerical solution of viscoelastic fluid flows , 1990 .

[10]  Gilberto Schleiniger,et al.  A remark on the Giesekus viscoelastic fluid , 1991 .

[11]  D. V. Boger,et al.  Further observations of elastic effects in tubular entry flows , 1986 .

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

[13]  E. Mitsoulis.  Numerical simulation of planar entry flow for a polyisobutylene solution using an integral constitutive equation , 1993 .

[14]  J. M. Rallison,et al.  Creeping flow of dilute polymer solutions past cylinders and spheres , 1988 .

[15]  R. Keiller.  Entry-flow calculations for the Oldroyd-B and FENE equations , 1993 .

[16]  H. Giesekus.  Stressing behaviour in simple shear flow as predicted by a new constitutive model for polymer fluids , 1983 .

[17]  G. Homsy,et al.  Numerical simulation of non-newtonian free shear flows at high reynolds numbers , 1994 .

[18]  R. Armstrong,et al.  Nonlinear dynamics of viscoelastic flow in axisymmetric abrupt contractions , 1991, Journal of Fluid Mechanics.

[19]  M. Crochet,et al.  Further Results On the Flow of a Viscoelastic Fluid Through An Abrupt Contraction , 1986 .

[20]  S. White,et al.  Review of the entry flow problem: Experimental and numerical , 1987 .

[21]  D. V. Boger.  Viscoelastic Flows Through Contractions , 1987 .

[22]  R. Armstrong,et al.  Modeling the rheology of polyisobutylene solutions , 1990 .

[23]  A. Beris,et al.  Finite element calculation of viscoelastic flow in a journal bearing: II. Moderate eccentricity , 1986 .

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

[25]  Brian J. Edwards,et al.  Thermodynamics of flowing systems : with internal microstructure , 1994 .

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