Arbitrary-Lagrangian-Eulerian discontinuous Galerkin schemes with a posteriori subcell finite volume limiting on moving unstructured meshes

TL;DR

This paper introduces high-order accurate Arbitrary-Lagrangian-Eulerian discontinuous Galerkin schemes with a novel a posteriori subcell finite volume limiter on moving unstructured meshes, suitable for complex hyperbolic PDEs including dissipative processes.

Contribution

It presents a new fully discrete one-step DG scheme with a posteriori limiting on moving meshes, combining high-order accuracy and adaptive troubled-cell detection.

Findings

Achieves up to fourth-order convergence in space and time.

Effectively handles complex hyperbolic PDEs with moving meshes.

Demonstrates applicability to magnetohydrodynamics and inertial confinement fusion flows.

Abstract

We present a new family of high order accurate fully discrete one-step Discontinuous Galerkin (DG) finite element schemes on moving unstructured meshes for the solution of nonlinear hyperbolic PDE in multiple space dimensions, which may also include parabolic terms in order to model dissipative transport processes. High order piecewise polynomials are adopted to represent the discrete solution at each time level and within each spatial control volume of the computational grid, while high order of accuracy in time is achieved by the ADER approach. In our algorithm the spatial mesh configuration can be defined in two different ways: either by an isoparametric approach that generates curved control volumes, or by a piecewise linear decomposition of each spatial control volume into simplex sub-elements. Our numerical method belongs to the category of direct Arbitrary-Lagrangian-Eulerian (ALE) schemes, where a space-time conservation formulation of the governing PDE system is considered and which already takes into account the new grid geometry directly during the computation of the numerical fluxes. Our new Lagrangian-type DG scheme adopts the novel a posteriori sub-cell finite volume limiter method, in which the validity of the candidate solution produced in each cell by an unlimited ADER-DG scheme is verified against a set of physical and numerical detection criteria. Those cells which do not satisfy all of the above criteria are flagged as troubled cells and are recomputed with a second order TVD finite volume scheme. The numerical convergence rates of the new ALE ADER-DG schemes are studied up to fourth order in space and time and several test problems are simulated. Finally, an application inspired by Inertial Confinement Fusion (ICF) type flows is considered by solving the Euler equations and the PDE of viscous and resistive magnetohydrodynamics (VRMHD).