ADER-ZIST - a mixed explicit ADER-DG and zonally implicit space-time DG approach for Cartesian cut cell meshes
Abstract
We present high order accurate mixed explicit-implicit discontinuous Galerkin (DG) schemes for the solution of hyperbolic systems of conservation laws on Cartesian cut cell meshes. Such meshes may contain arbitrarily small cut cells, which leads to the well-known small cell problem: the CFL condition of explicit schemes would severely restrict the admissible time step. To overcome this difficulty, we use high order explicit fully-discrete one-step ADER-DG schemes only on Cartesian cells away from the embedded boundaries. The zone close to the embedded boundary, which consists of cut cells and their direct face neighbors, is discretized with an implicit space-time DG scheme that is coupled to the explicit scheme directly via the numerical flux. The resulting scheme, called ADER-ZIST (ADER-DG with Zonally Implicit Space-Time DG), is conservative by construction, and its time step is determined only by the size of the Cartesian background cells. Curved embedded boundaries are represented by high order isoparametric sub-triangles inside the cut cells, and the nonlinear system of the implicit cells is solved via an inexact matrix-free Newton-Krylov method. For problems with shock waves, we add an artificial viscosity term whose strength is controlled by an oscillation-eliminating smoothness sensor. ADER-ZIST is assessed on a comprehensive set of test problems for the compressible Euler equations, ranging from smooth flows around straight and curved embedded boundaries to acoustic scattering and shock-dominated flows. The scheme essentially reaches its nominal order of accuracy N+1 for smooth solutions and for polynomial degrees up to N=5, on meshes with minimal volume fractions of the small cut cells down to the order of 10^-7.
Related
- A compact high-order and positivity-preserving active flux method for compressible Navier--Stokes equations on Cartesian meshes
- Wavenumber-explicit stability of three-dimensional Cartesian UPMLs with piecewise-constant profiles for Helmholtz and Maxwell equations
- Explicit Cartesian Coordinates for Three-Dimensional Clothoids
- A structure-preserving staggered semi-implicit four-split finite volume scheme for continuum mechanics on unstructured meshes
- Average Energy Dissipation Rates of Implicit-Explicit Runge-Kutta Methods for Gradient Flow Problems
BibTeX
@misc{dumbser2026aderzist,
title = {{ADER-ZIST - a mixed explicit ADER-DG and zonally implicit space-time DG approach for Cartesian cut cell meshes}},
author = {Michael Dumbser and Christian Engwer and Sandra May and Peter Munch},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.13245},
eprint = {2610.13245},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.13245},
note = {CC BY-NC-ND 4.0},
}
RIS
TY - UNPB
TI - ADER-ZIST - a mixed explicit ADER-DG and zonally implicit space-time DG approach for Cartesian cut cell meshes
AU - Michael Dumbser
AU - Christian Engwer
AU - Sandra May
AU - Peter Munch
PY - 2026
DA - 2026/10/09
PB - arXiv
DO - 10.48550/arxiv.2610.13245
UR - https://arxiv.org/abs/2610.13245
AB - We present high order accurate mixed explicit-implicit discontinuous Galerkin (DG) schemes for the solution of hyperbolic systems of conservation laws on Cartesian cut cell meshes. Such meshes may contain arbitrarily small cut cells, which leads to the well-known small cell problem: the CFL condition of explicit schemes would severely restrict the admissible time step. To overcome this difficulty, we use high order explicit fully-discrete one-step ADER-DG schemes only on Cartesian cells away from the embedded boundaries. The zone close to the embedded boundary, which consists of cut cells and their direct face neighbors, is discretized with an implicit space-time DG scheme that is coupled to the explicit scheme directly via the numerical flux. The resulting scheme, called ADER-ZIST (ADER-DG with Zonally Implicit Space-Time DG), is conservative by construction, and its time step is determined only by the size of the Cartesian background cells. Curved embedded boundaries are represented by high order isoparametric sub-triangles inside the cut cells, and the nonlinear system of the implicit cells is solved via an inexact matrix-free Newton-Krylov method. For problems with shock waves, we add an artificial viscosity term whose strength is controlled by an oscillation-eliminating smoothness sensor. ADER-ZIST is assessed on a comprehensive set of test problems for the compressible Euler equations, ranging from smooth flows around straight and curved embedded boundaries to acoustic scattering and shock-dominated flows. The scheme essentially reaches its nominal order of accuracy N+1 for smooth solutions and for polynomial degrees up to N=5, on meshes with minimal volume fractions of the small cut cells down to the order of 10^-7.
ER -