Numerical simulation of three-dimensional viscoelastic planar contraction flow using the software OpenFOAM

Abstract In this study, we have used a highly accurate and novel approach to solve numerically for true three-dimensional (3D) viscoelastic flows into sudden contractions. Motivation for this development has been for the advancement of 3D viscoelastic flows in complex geometries, where this new methodology is available as a general solver, written for the open source code OpenFOAM ( Weller et al., 1998 ). The proposed approach is able to include the multiple relaxation times of differential constitutive equations, and has been performed using the Finite Volume Method (FVM), based on a Discrete Elastic Viscous Split Stress (DEVSS) technique ( Fortin et al., 1997 ). In this work, both the Giesekus and Phan-Thien–Tanner (PTT) shear-thinning models were implemented to reproduce flow through a planar 4:1 contraction, where numerical convergence was achieved for a Weissenberg number ( We ) of 2.9. Direct comparison with experimental data and literature involving 2D and 3D numerical simulations shows this method to be both stable and effective.

[1]  S. White,et al.  The importance of extensional flow properties on planar entry flow patterns of polymer melts , 1986 .

[2]  Lars Davidson,et al.  A PRESSURE CORRECTION METHOD FOR UNSTRUCTURED MESHES WITH ARBITRARY CONTROL VOLUMES , 1996 .

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

[4]  Tim A. Osswald,et al.  Polymer Processing: Modeling and Simulation , 2006 .

[5]  H. Jasak,et al.  Preconditioned Linear Solvers for Large Eddy Simulation , 2007 .

[6]  A. Gosman,et al.  Solution of the implicitly discretised reacting flow equations by operator-splitting , 1986 .

[7]  Robert C. Armstrong,et al.  Dynamics of polymeric liquids: Fluid mechanics , 1987 .

[8]  D. Spalding,et al.  A calculation procedure for heat, mass and momentum transfer in three-dimensional parabolic flows , 1972 .

[9]  Gerrit W. M. Peters,et al.  3D viscoelastic analysis of a polymer solution in a complex flow , 1999 .

[10]  Weeratunge Malalasekera,et al.  An introduction to computational fluid dynamics - the finite volume method , 2007 .

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

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

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

[14]  F. Pinho,et al.  Visualizations of Boger fluid flows in a 4:1 square/square contraction , 2005 .

[15]  F. Pinho,et al.  Viscoelastic flow in a 3D square/square contraction: Visualizations and simulations , 2008 .

[16]  J. Song,et al.  Numerical simulation of viscoelastic flow through a sudden contraction using a type dependent difference method , 1987 .

[17]  Michel Fortin,et al.  Numerical analysis of the modified EVSS method , 1997 .

[18]  J. Dooley,et al.  Secondary motions in straight and tapered channels: Experiments and three-dimensional finite element simulation with a multimode differential viscoelastic model , 1999 .

[19]  N. Phan-Thien,et al.  Numerical study of secondary flows of viscoelastic fluid in straight pipes by an implicit finite volume method , 1995 .

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

[21]  Henk A. van der Vorst,et al.  Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear Systems , 1992, SIAM J. Sci. Comput..

[22]  D. V. Boger,et al.  Experimental dilemmas in non-Newtonian fluid mechanics and their theoretical resolution , 2000 .

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

[24]  M. Deville,et al.  Unsteady finite volume simulation of Oldroyd-B fluid through a three-dimensional planar contraction , 1997 .

[25]  Hrvoje Jasak,et al.  A tensorial approach to computational continuum mechanics using object-oriented techniques , 1998 .

[26]  T. Phillips,et al.  Viscoelastic flow through a planar contraction using a semi-Lagrangian finite volume method , 1999 .

[27]  Achi Brandt,et al.  Multilevel Adaptive Computations in Fluid Dynamics , 1980 .

[28]  D. M. Binding,et al.  On the use of flow through a contraction in estimating the extensional viscosity of mobile polymer solutions , 1988 .

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

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