An analysis of the time integration algorithms for the finite element solutions of incompressible Navier-Stokes equations based on a stabilised formulation

Abstract This work is concerned with the analysis of time integration procedures for the stabilised finite element formulation of unsteady incompressible fluid flows governed by the Navier–Stokes equations. The stabilisation technique is combined with several different implicit time integration procedures including both finite difference and finite element schemes. Particular attention is given to the generalised-α method and the linear discontinuous in time finite element scheme. The time integration schemes are first applied to two model problems, represented by a first order differential equation in time and the one dimensional advection–diffusion equation, and subjected to a detailed mathematical analysis based on the Fourier series expansion. In order to establish the accuracy and efficiency of the time integration schemes for the Navier–Stokes equations, a detailed computational study is performed of two standard numerical examples: unsteady flow around a cylinder and flow across a backward facing step. It is concluded that the semi-discrete generalised-α method provides a viable alternative to the more sophisticated and expensive space–time methods for simulations of unsteady flows of incompressible fluids governed by the Navier–Stokes equations.

[1]  G. Hulbert,et al.  A generalized-α method for integrating the filtered Navier–Stokes equations with a stabilized finite element method , 2000 .

[2]  O. Zienkiewicz,et al.  A general algorithm for compressible and incompressible flow—Part I. the split, characteristic‐based scheme , 1995 .

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

[4]  Ramon Codina,et al.  A general algorithm for compressible and incompressible flow—Part II. tests on the explicit form , 1995 .

[5]  Eugenio Oñate,et al.  A stabilized finite element method for incompressible viscous flows using a finite increment calculus formulation , 2000 .

[6]  Claes Johnson.  Numerical solution of partial differential equations by the finite element method , 1988 .

[7]  T. Hughes,et al.  A new finite element formulation for computational fluid dynamics: VI. Convergence analysis of the generalized SUPG formulation for linear time-dependent multi-dimensional advective-diffusive systems , 1987 .

[8]  Thomas J. R. Hughes,et al.  A new finite element formulation for computational fluid dynamics: IX. Fourier analysis of space-time Galerkin/least-squares algorithms , 1991 .

[9]  Djordje Peric,et al.  Computational modelling of viscoplastic fluids based on a stabilised finite element method , 2001 .

[10]  J. C. Simo,et al.  Unconditional stability and long-term behavior of transient algorithms for the incompressible Navier-Stokes and Euler equations , 1994 .

[11]  Thomas J. R. Hughes,et al.  A boundary integral modification of the Galerkin least squares formulation for the Stokes problem , 1994 .

[12]  Miguel R. Visbal,et al.  On spurious behavior of CFD simulations , 1999 .

[13]  Tayfun E. Tezduyar,et al.  Finite element stabilization parameters computed from element matrices and vectors , 2000 .

[14]  T. Tezduyar,et al.  A new strategy for finite element computations involving moving boundaries and interfaces—the deforming-spatial-domain/space-time procedure. I: The concept and the preliminary numerical tests , 1992 .

[15]  Claes Johnson,et al.  Computational Differential Equations , 1996 .

[16]  J. Z. Zhu,et al.  The finite element method , 1977 .

[17]  Kenneth E. Jansen,et al.  A better consistency for low-order stabilized finite element methods , 1999 .

[18]  H. Kreiss,et al.  Time-Dependent Problems and Difference Methods , 1996 .

[19]  T. Hughes,et al.  A new finite element formulation for computational fluid dynamics. X - The compressible Euler and Navier-Stokes equations , 1991 .

[20]  J. Szmelter.  Incompressible flow and the finite element method , 2001 .

[21]  R. Temam,et al.  Navier-stokes equations: Theory and approximation , 1998 .

[22]  Franco Brezzi,et al.  $b=\int g$ , 1997 .

[23]  Michel Fortin,et al.  Mixed and Hybrid Finite Element Methods , 2011, Springer Series in Computational Mathematics.

[24]  S. Mittal,et al.  A new strategy for finite element computations involving moving boundaries and interfaces—the deforming-spatial-domain/space-time procedure. II: Computation of free-surface flows, two-liquid flows, and flows with drifting cylinders , 1992 .

[25]  R. Codina.  Comparison of some finite element methods for solving the diffusion-convection-reaction equation , 1998 .

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

[27]  T. Hughes,et al.  A new finite element formulation for computational fluid dynamics: II. Beyond SUPG , 1986 .

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

[29]  A. Spence,et al.  Is the steady viscous incompressible two‐dimensional flow over a backward‐facing step at Re = 800 stable? , 1993 .

[30]  Jintai Chung,et al.  A Time Integration Algorithm for Structural Dynamics With Improved Numerical Dissipation: The Generalized-α Method , 1993 .

[31]  O. Pironneau.  Finite Element Methods for Fluids , 1990 .

[32]  Thomas J. R. Hughes,et al.  A space-time Galerkin/least-squares finite element formulation of the Navier-Stokes equations for moving domain problems , 1997 .

[33]  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 .

[34]  Claes Johnson,et al.  Finite element methods for linear hyperbolic problems , 1984 .