The locally conservative Galerkin (LCG) method for solving the incompressible Navier–Stokes equations

In this paper, the locally conservative Galerkin (LCG) method (Numer. Heat Transfer B Fundam. 2004; 46:357–370; Int. J. Numer. Methods Eng. 2007) has been extended to solve the incompressible Navier–Stokes equations. A new correction term is also incorporated to make the formulation to give identical results to that of the continuous Galerkin (CG) method. In addition to ensuring element-by-element conservation, the method also allows solution of the governing equations over individual elements, independent of the neighbouring elements. This is achieved within the CG framework by breaking the domain into elemental sub-domains. Although this allows discontinuous trial function field, we have carried out the formulation using the continuous trial function space as the basis. Thus, the changes in the existing CFD codes are kept to a minimum. The edge fluxes, establishing the continuity between neighbouring elements, are calculated via a post-processing step during the time-stepping operation. Therefore, the employed formulation needs to be carried out using either a time-stepping or an equivalent iterative scheme that allows post-processing of fluxes. The time-stepping algorithm employed in this paper is based on the characteristic-based split (CBS) scheme. Both steady- and unsteady-state examples presented show that the element-by-element formulation employed is accurate and robust. Copyright © 2007 John Wiley & Sons, Ltd.

[1]  K. Morgan,et al.  Three‐dimensional incompressible flow calculations using the characteristic based split (CBS) scheme , 2004 .

[2]  Robert L. Lee,et al.  THE CONSISTENT GALERKIN FEM FOR COMPUTING DERIVED BOUNDARY QUANTITIES IN THERMAL AND/OR FLUIDS PROBLEMS , 1987 .

[3]  P. Nithiarasu,et al.  An element‐wise, locally conservative Galerkin (LCG) method for solving diffusion and convection–diffusion problems , 2008 .

[4]  P. Nithiarasu A Matrix Free Fractional Step Method for Static and Dynamic Incompressible Solid Mechanics , 2006 .

[5]  P. Raviart,et al.  On a Finite Element Method for Solving the Neutron Transport Equation , 1974 .

[6]  U. Ghia,et al.  High-Re solutions for incompressible flow using the Navier-Stokes equations and a multigrid method , 1982 .

[7]  G. Chavent,et al.  A finite-element method for the 1-D water flooding problem with gravity , 1982 .

[8]  Thomas J. R. Hughes,et al.  Consistent finite element calculations of boundary and internal fluxes , 1998 .

[9]  Arnaud G. Malan,et al.  An improved unsteady, unstructured, artificial compressibility, finite volume scheme for viscous incompressible flows: Part I. Theory and implementation , 2002 .

[10]  Ramon Codina,et al.  Numerical comparison of CBS and SGS as stabilization techniques for the incompressible Navier–Stokes equations , 2006 .

[11]  A mixed finite element method for boundary flux computation , 1986 .

[12]  Thomas J. R. Hughes,et al.  Conservation properties for the Galerkin and stabilised forms of the advection–diffusion and incompressible Navier–Stokes equations , 2005 .

[13]  Chi-Wang Shu,et al.  TVB Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws III: one-dimensional systems , 1989 .

[14]  G. Chavent,et al.  The Local Projection P0P1-Discontinuous-Galerkin-Finite Element Method for Scalar Conservation Laws , 1987 .

[15]  O. C. Zienkiewicz,et al.  Natural convection in porous medium‐fluid interface problems ‐ A finite element analysis by using the CBS procedure , 2001 .

[16]  Thomas J. R. Hughes,et al.  The Continuous Galerkin Method Is Locally Conservative , 2000 .

[17]  Spencer J. Sherwin,et al.  On discontinuous Galerkin methods , 2003 .

[18]  Ann L Gaitonde,et al.  A dual-time method for two-dimensional unsteady incompressible flow calculations , 1998 .

[19]  P. Nithiarasu,et al.  Laminar and turbulent flow calculations through a model human upper airway using unstructured meshes , 2006 .

[20]  O. C. Zienkiewicz,et al.  Analysis of an explicit and matrix free fractional step method for incompressible flows , 2006 .

[21]  Perumal Nithiarasu,et al.  A unified fractional step method for compressible and incompressible flows, heat transfer and incompressible solid mechanics , 2008 .

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

[23]  Arnaud G. Malan,et al.  An improved unsteady, unstructured, artificial compressibility, finite volume scheme for viscous incompressible flows: Part II. Application , 2002 .

[24]  K. Morgan,et al.  Petrov-Galerkin solutions of the incompressible Navier-Stokes equations in primitive variables with adaptive remeshing , 1993 .

[25]  O. C. Zienkiewicz,et al.  The solution of non‐linear hyperbolic equation systems by the finite element method , 1984 .

[26]  Chi-Wang Shu,et al.  The Runge-Kutta Discontinuous Galerkin Method for Conservation Laws V , 1998 .

[27]  Chi-Wang Shu,et al.  The Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws. IV. The multidimensional case , 1990 .

[28]  Chi-Wang Shu,et al.  TVB Runge-Kutta local projection discontinuous galerkin finite element method for conservation laws. II: General framework , 1989 .

[29]  Dimitris E. Papantonis,et al.  A characteristic-based method for incompressible flows , 1994 .

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

[31]  P. Nithiarasu A SIMPLE LOCALLY CONSERVATIVE GALERKIN (LCG) FINITE-ELEMENT METHOD FOR TRANSIENT CONSERVATION EQUATIONS , 2004 .

[32]  Chi-Wang Shu,et al.  The Runge-Kutta local projection $P^1$-discontinuous-Galerkin finite element method for scalar conservation laws , 1988, ESAIM: Mathematical Modelling and Numerical Analysis.

[33]  W. H. Reed,et al.  Triangular mesh methods for the neutron transport equation , 1973 .

[34]  T. Hughes,et al.  Finite element method for high-speed flows - Consistent calculation of boundary flux , 1987 .

[35]  P. Nithiarasu An efficient artificial compressibility (AC) scheme based on the characteristic based split (CBS) method for incompressible flows , 2003 .

[36]  Guido Kanschat,et al.  The local discontinuous Galerkin method for linearized incompressible fluid flow: a review , 2005 .

[37]  O. C. Zienkiewicz,et al.  The Characteristic‐Based Split (CBS) scheme—a unified approach to fluid dynamics , 2006 .