Parallelized FVM algorithm for three-dimensional viscoelastic flows

Abstract A parallel implementation for the finite volume method (FVM) for three-dimensional (3D) viscoelastic flows is developed on a distributed computing environment through Parallel Virtual Machine (PVM). The numerical procedure is based on the SIMPLEST algorithm using a staggered FVM discretization in Cartesian coordinates. The final discretized algebraic equations are solved with the TDMA method. The parallelisation of the program is implemented by a domain decomposition strategy, with a master/slave style programming paradigm, and a message passing through PVM. A load balancing strategy is proposed to reduce the communications between processors. The three-dimensional viscoelastic flow in a rectangular duct is computed with this program. The modified Phan-Thien–Tanner (MPTT) constitutive model is employed for the equation system closure. Computing results are validated on the secondary flow problem due to non-zero second normal stress difference N2. Three sets of meshes are used, and the effect of domain decomposition strategies on the performance is discussed. It is found that parallel efficiency is strongly dependent on the grid size and the number of processors for a given block number. The convergence rate as well as the total efficiency of domain decomposition depends upon the flow problem and the boundary conditions. The parallel efficiency increases with increasing problem size for given block number. Comparing to two-dimensional flow problems, 3D parallelized algorithm has a lower efficiency owing to largely overlapped block interfaces, but the parallel algorithm is indeed a powerful means for large scale flow simulations.

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

[2]  P. S. Larsen,et al.  Secondary flows in straight ducts of rectangular cross section , 1991 .

[3]  D. B. Spalding,et al.  A general purpose computer program for multi-dimensional one- and two-phase flow , 1981 .

[4]  G. P Sasmal,et al.  A finite volume approach for calculation of viscoelastic flow through an abrupt axisymmetric contraction , 1995 .

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

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

[7]  Jack Dongarra,et al.  PVM: Parallel virtual machine: a users' guide and tutorial for networked parallel computing , 1995 .

[8]  Joel H. Ferziger,et al.  Computational methods for fluid dynamics , 1996 .

[9]  S. Patankar.  Numerical Heat Transfer and Fluid Flow , 2018, Lecture Notes in Mechanical Engineering.

[10]  Nhan Phan-Thien,et al.  VISCOELASTIC FLOW BETWEEN ECCENTRIC ROTATING CYLINDERS: UNSTRUCTURED CONTROL VOLUME METHOD , 1996 .

[11]  N. Phan-Thien,et al.  A domain decomposition implementation of the SIMPLE method with PVM , 1997 .

[12]  John R. Whiteman,et al.  Numerical modelling of viscoelastic liquids using a finite-volume method , 1992 .

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

[14]  A. Chorin.  A Numerical Method for Solving Incompressible Viscous Flow Problems , 1997 .

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

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

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

[18]  J. Yoo,et al.  A finite volume technique to simulate the flow of a viscoelastic fluid , 1991 .

[19]  M. Perić,et al.  Computation of fluid flow with a parallel multigrid solver , 1993 .

[20]  Dionysis Assimacopoulos,et al.  A finite volume approach in the simulation of viscoelastic expansion flows , 1998 .

[21]  F. Pinho,et al.  Numerical simulation of non-linear elastic flows with a general collocated finite-volume method , 1998 .

[22]  Jeff Nichols,et al.  Implementation of the direct SCF and RPA methods on loosely coupled networks of workstations , 1993, J. Comput. Chem..

[23]  Squeezing a viscoelastic liquid from a wedge: an exact solution , 1984 .

[24]  B. A. Tanyi,et al.  ITERATIVE SOLUTION OF INCOMPRESSIBLE NAVIER–STOKES EQUATIONS ON THE MEIKO COMPUTING SURFACE , 1996 .

[25]  Xiaolin Luo.  A control volume approach for integral viscoelastic models and its application to contraction flow of polymer melts , 1996 .

[26]  Roland Glowinski,et al.  A distributed Lagrange multiplier/fictitious domain method for viscoelastic particulate flows , 2000 .

[27]  Daniel D. Joseph,et al.  Numerical simulation of viscoelastic flow past a cylinder , 1990 .

[28]  B. Lakshminarayana,et al.  An Assessment of Computational Fluid Dynamic Techniques in the Analysis and Design of Turbomachinery—The 1990 Freeman Scholar Lecture , 1991 .

[29]  N. Phan-Thien,et al.  The flow of an Oldroyd-B fluid past a cylinder in a channel: adaptive viscosity vorticity (DAVSS-ω) formulation , 1999 .

[30]  Tayfun E. Tezduyar,et al.  Flow simulation and high performance computing , 1996 .

[31]  J. Dooley,et al.  On the development of secondary motions in straight channels induced by the second normal stress difference: experiments and simulations , 1997 .

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