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.

A structure preserving discretization for the Derrida-Lebowitz-Speer-Spohn equation based on diffusive transport

Convergence of spectral discretization for the flow of diffeomorphisms

Existence of weak solutions to a class of fourth order partial differential equations with Wasserstein gradient structure

A second-order structure-preserving discretization for the Cahn-Hilliard/Allen-Cahn system with cross-kinetic coupling

Structure Preserving Primal Dual Methods for Gradient Flows with Nonlinear Mobility Transport Distances

Structure-preserving Local Discontinuous Galerkin method for nonlinear cross-diffusion systems

Augmented saddle point formulation of the steady-state Stefan--Maxwell diffusion problem

Discrete approximation of the minimizing movement scheme for evolution equations of Wasserstein gradient flow type with nonlinear mobility

A discrete de Rham discretization of interface diffusion problems with application to the Leaky Dielectric Model

A second-order fully-balanced structure-preserving variational discretization scheme for the Cahn-Hilliard Navier-Stokes system

Unconditional convergence for discretizations of dynamical optimal transport

Symmetry-preserving finite-difference discretizations of arbitrary order on structured curvilinear staggered grids

Structure-Preserving Discretization of Incompressible Fluids

Fully-discrete provably Lyapunov consistent discretizations for convection-diffusion-reaction PDE systems

High Order Asymptotic Preserving DG-IMEX Schemes for Discrete-Velocity Kinetic Equations in a Diffusive Scaling

Decoupled structure-preserving discretization of incompressible MHD equations with general boundary conditions

On discretization in time in simulations of particulate flows

A structure-preserving scheme for computing effective diffusivity and anomalous diffusion phenomena of random flows

Structure preserving schemes for a class of Wasserstein gradient flows

Asymptotic Preserving Semi-Lagrangian Discontinuous Galerkin Methods for Multiscale Kinetic Transport Equations

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