Embedded discontinuous Galerkin transport schemes with localised limiters
arXiv:1509.04431 · doi:10.1016/j.jcp.2016.02.021
Abstract
Motivated by finite element spaces used for representation of temperature in the compatible finite element approach for numerical weather prediction, we introduce locally bounded transport schemes for (partially-)continuous finite element spaces. The underlying high-order transport scheme is constructed by injecting the partially-continuous field into an embedding discontinuous finite element space, applying a stable upwind discontinuous Galerkin (DG) scheme, and projecting back into the partially-continuous space; we call this an embedded DG scheme. We prove that this scheme is stable in L2 provided that the underlying upwind DG scheme is. We then provide a framework for applying limiters for embedded DG transport schemes. Standard DG limiters are applied during the underlying DG scheme. We introduce a new localised form of element-based flux-correction which we apply to limiting the projection back into the partially-continuous space, so that the whole transport scheme is bounded. We provide details in the specific case of tensor-product finite element spaces on wedge elements that are discontinuous P1/Q1 in the horizontal and continuous P2 in the vertical. The framework is illustrated with numerical tests.
References in corpus (1)
Cited by in corpus (10)
- Algebraic entropy fixes and convex limiting for continuous finite element discretizations of scalar hyperbolic conservation laws
- Compatible finite element spaces for geophysical fluid dynamics
- Subcell flux limiting for high-order Bernstein finite element discretizations of scalar hyperbolic conservation laws
- A Compatible Finite Element Discretisation for the Moist Compressible Euler Equations
- A new perspective on flux and slope limiting in discontinuous Galerkin methods for hyperbolic conservation laws
- Vertical slice modelling of nonlinear Eady waves using a compatible finite element method
- The Recovered Space Advection Scheme for Lowest-Order Compatible Finite Element Methods
- Active flux methods for hyperbolic conservation laws -- flux vector splitting and bound-preservation
- Second Order Unconditional Positive Preserving Schemes for Non-equilibrium Reactive Flows with Mass and Mole Balance
- Improving the accuracy of discretisations of the vector transport equation on the lowest-order quadrilateral Raviart-Thomas finite elements