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 note on a NSFD discretization of a Sobolev type PDE

A note on an NSFD scheme for a mathematical model of respiratory virus transmission

Fully-Discretely Nonlinearly-Stable Flux Reconstruction Methods for Compressible Flows

A discretization scheme for path-dependent FBSDEs

Compact finite-difference scheme for some Sobolev type equations with Dirichlet boundary conditions

Parametric Finite Element Discretization of the Surface Stokes Equations

Stable semi-implicit SDC methods for conservation laws

Discretization of fractional fully nonlinear equations by powers of discrete Laplacians

Efficient finite element method for 2D singularly perturbed parabolic convection diffusion problems with discontinuous source term

Numerical Dynamics of a Nonstandard Finite Difference Method for a Class of Delay Differential Equations

Nyström discretization of integrodifference equations: numerical continuation of periodic solutions and Floquet multipliers

A finite difference scheme for two-dimensional singularly perturbed convection-diffusion problem with discontinuous source term

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

Scaling-robust built-in a posteriori error estimation for discontinuous least-squares finite element methods

A dual-pairing summation-by-parts finite difference framework for nonlinear conservation laws

High-order nonstandard finite difference methods preserving dynamical properties of one-dimensional dynamical systems

Finite Spectral Essentially Non-Oscillatory Scheme

Extending nonstandard finite difference scheme rules to systems of nonlinear ODEs with constant coefficients

A PDE approach to fractional diffusion: a space-fractional wave equation

A data-driven learned discretization approach in finite volume schemes for hyperbolic conservation laws and varying boundary conditions

Constraint-satisfying Krylov solvers for structure-preserving discretisations