Locally Structure-Preserving div-curl operators for high order Discontinuous Galerkin schemes
arXiv:2206.08609 · doi:10.1016/j.jcp.2023.112130
Abstract
We propose a novel Structure-Preserving Discontinuous Galerkin (SPDG) operator that recovers at the discrete level the algebraic property related to the divergence of the curl of a vector field, which is typically referred to as div-curl problem. A staggered Cartesian grid is adopted in 3D, where the vector field is naturally defined at the corners of the control volume, while its curl is evaluated as a cell-centered quantity. Firstly, the curl operator is rewritten as the divergence of a tensor, hence allowing compatible finite difference schemes to be devised and to be proven to mimic the algebraic div-curl property. Successively, a high order DG divergence operator is built upon integration by parts, so that the structure-preserving finite difference div-curl operator is exactly retrieved for first order discretizations. We further demonstrate that the novel SPDG schemes are capable of obtaining a zero div-curl identity with machine precision from second up to sixth order accuracy. In a second part, we show the applicability of these SPDG methods by solving the incompressible Navier-Stokes equations written in vortex-stream formulation. This hyperbolic system deals with divergence-free involutions related to the velocity and vorticity field as well as to the stream function, thus it provides an ideal setting for the validation of the novel schemes. A compatible discretization of the numerical viscosity is also proposed in order to maintain the structure-preserving property of the div-curl DG operators even in the presence of artificial or physical dissipative terms. Finally, to overcome the time step restriction dictated by the viscous sub-system, Implicit-Explicit (IMEX) Runge-Kutta time stepping techniques are tailored to handle the SPDG framework.
References in corpus (14)
- An Unsplit Godunov Method for Ideal MHD via Constrained Transport
- Structure and structure-preserving algorithms for plasma physics
- A staggered space-time discontinuous Galerkin method for the three-dimensional incompressible Navier-Stokes equations on unstructured tetrahedral meshes
- A mass, energy, enstrophy and vorticity conserving (MEEVC) mimetic spectral element discretization for the 2D incompressible Navier-Stokes equations
- High order pressure-based semi-implicit IMEX schemes for the 3D Navier-Stokes equations at all Mach numbers
- A divergence-free semi-implicit finite volume scheme for ideal, viscous and resistive magnetohydrodynamics
- A structure-preserving staggered semi-implicit finite volume scheme for continuum mechanics
- An embedded variable step IMEX scheme for the incompressible Navier-Stokes equations
- Reinterpretation and Extension of Entropy Correction Terms for Residual Distribution and Discontinuous Galerkin Schemes: Application to Structure Preserving Discretization
- On GLM curl cleaning for a first order reduction of the CCZ4 formulation of the Einstein field equations
- High order ADER schemes and GLM curl cleaning for a first order hyperbolic formulation of compressible flow with surface tension
- Continuous finite element subgrid basis functions for Discontinuous Galerkin schemes on unstructured polygonal Voronoi meshes
- A Mixed Mimetic Spectral Element Model of the 3D Compressible Euler Equations on the Cubed Sphere
- Div-Curl Problems and -regular Stream Functions in 3D Lipschitz Domains