On the design of general-purpose flux limiters for implicit FEM with a consistent mass matrix

The algebraic flux correction (AFC) paradigm is extended to finite element discretizations with a consistent mass matrix. It is shown how to render an implicit Galerkin scheme positivity-preserving and remove excessive artificial diffusion in regions where the solution is sufficiently smooth. To this end, the original discrete operators are modified in a mass-conserving fashion so as to enforce the algebraic constraints to be satisfied by the numerical solution. A node-oriented limiting strategy is employed to control the raw antidiffusive fluxes which consist of a convective part and a contribution of the consistent mass matrix. The former offsets the artificial diffusion due to ‘upwinding’ of the spatial differential operator and lends itself to an upwind-biased flux limiting. The latter eliminates the error induced by mass lumping and calls for the use of a symmetric flux limiter. The concept of a target flux and a new definition of upper/lower bounds make it possible to combine the advantages of algebraic FCT and TVD schemes introduced previously by the author and his coworkers. Unlike other high-resolution schemes for unstructured meshes, the new algorithm reduces to a consistent (high-order) Galerkin scheme in smooth regions and is designed to provide an optimal treatment of both stationary and time-dependent problems. Its performance is illustrated by application to the linear advection equation for a number of 1D and 2D configurations.

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

[2]  Stefan Turek,et al.  Flux correction tools for finite elements , 2002 .

[3]  Elaine S. Oran,et al.  Numerical Simulation of Reactive Flow , 1987 .

[4]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

[5]  Robert J. MacKinnon,et al.  Positivity‐preserving, flux‐limited finite‐difference and finite‐element methods for reactive transport , 2003 .

[6]  Philip L. Roe,et al.  Multidimensional upwind schemes based on fluctuation-splitting for systems of conservation laws , 1993 .

[7]  R. LeVeque.  High-resolution conservative algorithms for advection in incompressible flow , 1996 .

[8]  Stefan Turek,et al.  Efficient Solvers for Incompressible Flow Problems - An Algorithmic and Computational Approach , 1999, Lecture Notes in Computational Science and Engineering.

[9]  Dmitri Kuzmin,et al.  Positive finite element schemes based on the flux-corrected transport procedure , 2001 .

[10]  Martin Berzins,et al.  Variable-order finite elements and positivity preservation for hyperbolic PDEs , 2004 .

[11]  Alexander Sokolichin.  Mathematische Modellbildung und numerische Simulation von Gas-Flüssigkeits-Blasenströmungen , 2004 .

[12]  Martin Berzins,et al.  Modified mass matrices and positivity preservation for hyperbolic and parabolic PDEs , 2001 .

[13]  Dmitri Kuzmin,et al.  Algebraic Flux Correction I. Scalar Conservation Laws , 2005 .

[14]  J. Peraire,et al.  TVD ALGORITHMS FOR THE SOLUTION OF THE COMPRESSIBLE EULER EQUATIONS ON UNSTRUCTURED MESHES , 1994 .

[15]  L. Quartapelle,et al.  An analysis of time discretization in the finite element solution of hyperbolic problems , 1987 .

[16]  D. Kuzmin,et al.  High-resolution FEM-TVD schemes based on a fully multidimensional flux limiter , 2004 .

[17]  Antony Jameson,et al.  Positive schemes and shock modelling for compressible flows , 1995 .

[18]  Alain Dervieux,et al.  Construction of TVD-like Artificial Viscosities on Two-Dimensional Arbitrary FEM Grids , 1993 .

[19]  D. Kuzmin,et al.  Algebraic Flux Correction III. Incompressible Flow Problems , 2005 .

[20]  K. Morgan,et al.  FEM-FCT - Combining unstructured grids with high resolution. [Flux Corrected Transport , 1988 .

[21]  Dmitri Kuzmin,et al.  Algebraic Flux Correction II. Compressible Euler Equations , 2005 .

[22]  Stefan Turek,et al.  High-resolution FEM?FCT schemes for multidimensional conservation laws , 2004 .

[23]  T. Jongen,et al.  Design of an unconditionally stable, positive scheme for the K-ϵ and two-layer turbulence models , 1997 .

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

[25]  S. Zalesak.  Fully multidimensional flux-corrected transport algorithms for fluids , 1979 .

[26]  Philip L. Roe,et al.  Multidimensional upwinding: Its relation to finite elements , 1995 .

[27]  J. Peraire,et al.  Finite Element Flux-Corrected Transport (FEM-FCT) for the Euler and Navier-Stokes equations , 1987 .

[28]  Rainald Löhner,et al.  The design of flux-corrected transport (fct) algorithms on structured grids , 2005 .

[29]  Dmitri Kuzmin,et al.  Adaptive mesh refinement for high‐resolution finite element schemes , 2006 .

[30]  Jozef Brilla.  Error analysis for Laplace transform —Finite element solution of hyperbolic equations , 1983 .

[31]  A. Jameson.  Computational algorithms for aerodynamic analysis and design , 1993 .

[32]  M. Lukáčová-Medvid'ová,et al.  Combined finite element-finite volume solution of compressible flow , 1995 .

[33]  Philip L. Roe,et al.  A Well-Behaved TVD Limiter for High-Resolution Calculations of Unsteady Flow , 1997 .

[34]  Barry Koren,et al.  A robust upwind discretization method for advection, diffusion and source terms , 1993 .

[35]  K. Hain.  The partial donor cell method , 1978 .

[36]  George Em Karniadakis,et al.  The Development of Discontinuous Galerkin Methods , 2000 .

[37]  Vittorio Selmin.  Finite element solution of hyperbolic equations II. Two dimensional case , 1987 .