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.

Solving Minimal Residual Methods in $W^{-1,p'}$ with large Exponents $p$

Solving Minimal Residual Methods in with Large Exponents p

$W^{m,p}$-Solution ($p\geq2$) of Linear Degenerate Backward Stochastic Partial Differential Equations in the Whole Space

Minimal residual discretization of a class of fully nonlinear elliptic PDE

The Dirichlet problem for the $p(x)$-Laplacian with unbounded exponent $p(x)$

High-order exponential time differencing multi-resolution alternative finite difference WENO methods for nonlinear degenerate parabolic equations

Dissipative WENO stabilization of high-order discontinuous Galerkin methods for hyperbolic problems

The limiting behavior of solutions to p-Laplacian problems with convection and exponential terms

Super-localized orthogonal decomposition for convection-dominated diffusion problems

Wm,P-Solution (P >= 2) Of Linear Degenerate Backward Stochastic Partial Differential Equations In The Whole Space

A convenient inclusion of inhomogeneous boundary conditions in minimal residual methods

Strongly consistent low-dissipation WENO schemes for finite elements

Interior Sobolev regularity for fully nonlinear parabolic equations

Localized Nodal Solutions for P-Laplacian Equations with Critical Exponents

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

Discrete-to-Continuum Rates of Convergence for $p$-Laplacian Regularization

Positive solutions to nonlinear elliptic problems involving Sobolev exponent

$W^{1,p}$ Regularity for the Conormal Derivative Problem with Parabolic BMO Nonlinearity in Reifenberg Domains

A Novel Low-Dimensional Method for Analytically Solving Partial Differential Equations

Stable solution and extremal solution for the fractional $p$-Laplacian equation

Normalized solutions for p-Laplacian equations with a $$L^{2}$$-supercritical growth