Abstract:This work introduces a novel discontinuity-tracking framework for resolving discontinuous solutions of conservation laws with high-order numerical discretizations that support inter-element solution discontinuities, such as discontinuous Galerkin methods. The proposed method aims to align inter-element boundaries with discontinuities in the solution by deforming the computational mesh. A discontinuity-aligned mesh ensures the discontinuity is represented through inter-element jumps while smooth basis functions interior to elements are only used to approximate smooth regions of the solution, thereby avoiding Gibbs' phenomena that create well-known stability issues. Therefore, very coarse high-order discretizations accurately resolve the piecewise smooth solution throughout the domain, provided the discontinuity is tracked. Central to the proposed discontinuity-tracking framework is a discrete PDE-constrained optimization formulation that simultaneously aligns the computational mesh with discontinuities in the solution and solves the discretized conservation law on this mesh. The optimization objective is taken as a combination of the the deviation of the finite-dimensional solution from its element-wise average and a mesh distortion metric to simultaneously penalize Gibbs' phenomena and distorted meshes. We advocate a gradient-based, full space solver where the mesh and conservation law solution converge to their optimal values simultaneously and therefore never require the solution of the discrete conservation law on a non-aligned mesh. The merit of the proposed method is demonstrated on a number of one- and two-dimensional model problems including 2D supersonic flow around a bluff body. We demonstrate optimal $\mathcal{O}(h^{p+1})$ convergence rates in the $L^1$ norm for up to polynomial order $p=6$ and show that accurate solutions can be obtained on extremely coarse meshes.

H -Adaptive Mesh Method with Double Tolerance Adaptive Strategy for Hyperbolic Conservation Laws

Mesh Adaptation for Curing the Pathological Behaviors of an Upwind Scheme

Mathematical Mesh Adaptation of Numerical Manifold Element Method

Adaptive Mesh Methods for One- and Two-Dimensional Hyperbolic Conservation Laws.

Adjoint-based Error Estimation and Mesh Adaptation for the Correction Procedure Via Reconstruction Method

An H-Adaptive Mesh Method for Boltzmann-BGK/hydrodynamics Coupling

The Multi-Dimensional Limiters for Solving Hyperbolic Conservation Laws on Unstructured Grids II: Extension to High Order Finite Volume Schemes

On Multi-Mesh H -Adaptive Methods

A compact simple HWENO scheme with ADER time discretization for hyperbolic conservation laws I: structured meshes

Adaptive Distribute Points Technique for Convection Dominated Problems Combined with Meshless Method

A quasi-Lagrangian moving mesh discontinuous Galerkin method for hyperbolic conservation laws

Finite difference alternative unequal-sized weighted essentially non-oscillatory schemes for hyperbolic conservation laws

A Self-Adaptive Numerical Method to Solve Convection-Dominated Diffusion Problems

Moving Mesh Methods in Multiple Dimensions Based on Harmonic Maps

A Path-Conservative ADER Discontinuous Galerkin Method for Non-Conservative Hyperbolic Systems: Applications to Shallow Water Equations

Robust High-Order Space-Time Conservative Schemes for Solving Conservation Laws on Hybrid Meshes

A fifth-order finite difference HWENO scheme combined with limiter for hyperbolic conservation laws

An optimization-based approach for high-order accurate discretization of conservation laws with discontinuous solutions

An Anisotropic hp-mesh Adaptation Method for Time-Dependent Problems Based on Interpolation Error Control

Towards the efficient calculation of quantity of interest from steady Euler equations I: a dual-consistent DWR-based h-adaptive Newton-GMG solver