Parallel finite element computation of unsteady incompressible flows

A parallel semi-explicit iterative finite element computational procedure for modelling unsteady incompressible fluid flows is presented. During the procedure, element flux vectors are calculated in parallel and then assembled into global flux vectors. Equilibrium iterations which introduce some ‘local implicitness’ are performed at each time step. The number of equilibrium iterations is governed by an implicitness parameter. The present technique retains the advantages of purely explicit schemes, namely (i) the parallel speed-up is equal to the number of parallel processors if the small communication overhead associated with purely explicit schemes is ignored and (ii) the computation time as well as the core memory required is linearly proportional to the number of elements. The incompressibility condition is imposed by using the artificial compressibility technique. A pressure-averaging technique which allows the use of equal-order interpolations for both velocity and pressure, this simplifying the formulation, is employed. Using a standard Galerkin approximation, three benchmark steady and unsteady problems are solved to demonstrate the accuracy of the procedure. In all calculations the Reynolds number is less than 500. At these Reynolds numbers it was found that the physical dissipation is sufficient to stabilize the convective term with no need for additional upwind-type dissipation. © 1998 John Wiley & Sons, Ltd.

[1]  S. Mittal,et al.  Incompressible flow computations with stabilized bilinear and linear equal-order-interpolation velocity-pressure elements , 1992 .

[2]  Miguel Cervera,et al.  The intrinsic time for the streamline upwind/Petrov-Galerkin formulation using quadratic elements , 1992 .

[3]  H. Saunders,et al.  Finite element procedures in engineering analysis , 1982 .

[4]  T. Hughes,et al.  A new finite element formulation for computational fluid dynamics: I. Symmetric forms of the compressible Euler and Navier—Stokes equations and the second law of thermodynamics , 1986 .

[5]  Charbel Farhat,et al.  Large Scale FE Parallel Nonlinear Computations Using a Homotopy Method , 1987, PPSC.

[6]  I. Babuska Error-bounds for finite element method , 1971 .

[7]  J. Hyun,et al.  Numerical Simulation of Flows Driven by a Torsionally Oscillating Lid in a Square Cavity , 1992 .

[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]  B. Armaly,et al.  Experimental and theoretical investigation of backward-facing step flow , 1983, Journal of Fluid Mechanics.

[10]  O. C. Zienkiewicz,et al.  Quadratic finite element schemes for two-dimensional convective-transport problems , 1977 .

[11]  T. Hughes,et al.  A new finite element formulation for computational fluid dynamics: V. Circumventing the Babuscka-Brezzi condition: A stable Petrov-Galerkin formulation of , 1986 .

[12]  A. Roshko On the development of turbulent wakes from vortex streets , 1953 .

[13]  T. Belytschko,et al.  A generalized Galerkin method for steady convection-diffusion problems with application to quadratic shape function elements , 1985 .

[14]  Wing Kam Liu,et al.  Finite Element Analysis of Incompressible Viscous Flows by the Penalty Function Formulation , 1979 .

[15]  M. Bercovier,et al.  A finite element for the numerical solution of viscous incompressible flows , 1979 .

[16]  O. C. Zienkiewicz,et al.  A general explicit or semi-explicit algorithm for compressible and incompressible flows , 1992 .

[17]  P. Fischer,et al.  PARALLEL SIMULATION OF VISCOUS INCOMPRESSIBLE FLOWS , 1994 .

[18]  Thomas J. R. Hughes,et al.  Scalability of finite element applications on distributed-memory parallel computers , 1994 .

[19]  J. Tinsley Oden,et al.  Optimal h-p finite element methods , 1994 .

[20]  T. Hughes,et al.  Stabilized finite element methods. I: Application to the advective-diffusive model , 1992 .

[21]  On quadratic elements in finite element solutions of steady‐state convection—diffusion equation , 1980 .

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

[23]  Peter Hansbo,et al.  A velocity pressure streamline diffusion finite element method for Navier-Stokes equations , 1990 .

[24]  A. J. Baker,et al.  Finite element computational fluid mechanics , 1983 .

[25]  T. Hughes,et al.  The Galerkin/least-squares method for advective-diffusive equations , 1988 .

[26]  M. F. Webster,et al.  Numerical solution of unsteady viscous flows , 1992 .

[27]  J. Donea A Taylor–Galerkin method for convective transport problems , 1983 .

[28]  Thomas J. R. Hughes,et al.  Finite element methods for first-order hyperbolic systems with particular emphasis on the compressible Euler equations , 1984 .

[29]  Peter Hansbo,et al.  A velocity-pressure streamline diffusion finite element method for the incompressible Navier-Stokes equation , 1990 .

[30]  Charbel Farhat,et al.  A general approach to nonlinear FE computations on shared-memory multiprocessors , 1989 .

[31]  Andrew R. Mitchell,et al.  Upwinding of high order Galerkin methods in conduction−convection problems , 1978 .