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 robust, discrete-gradient descent procedure for optimisation with time-dependent PDE and norm constraints

An Alternating Structure-Adapted Bregman Proximal Gradient Descent Algorithm for Constrained Nonconvex Nonsmooth Optimization Problems and Its Inertial Variant.

Adjoint-Based Enforcement of State Constraints in PDE Optimization Problems

Optimization on manifolds: A symplectic approach

A geometric integration approach to smooth optimisation: Foundations of the discrete gradient method

Goal-oriented adaptivity for a conforming residual minimization method in a dual discontinuous Galerkin norm

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

An SDE Perspective on Stochastic Inertial Gradient Dynamics with Time-Dependent Viscosity and Geometric Damping

An Eulerian–Lagrangian Method for Optimization Problems Governed by Multidimensional Nonlinear Hyperbolic PDEs

Global Optimization with Orthogonality Constraints Via Stochastic Diffusion on Manifold

An Adjoint Method for a High-Order Discretization of Deforming Domain Conservation Laws for Optimization of Flow Problems

A relaxation-based probabilistic approach for PDE-constrained optimization under uncertainty with pointwise state constraints

Adjoint-Based Design Optimization of Stability Constrained Systems

The Stochastic Steepest Descent Method for Robust Optimization in Banach Spaces

A Novel Convergence Enhancement Method Based on Online Dimension Reduction Optimization.

A narrow-stencil framework for convergent numerical approximations of fully nonlinear second order PDEs

Numerical Methods for Optimal Control Problems with SPDEs

Nonsmooth Optimization over the Stiefel Manifold and Beyond: Proximal Gradient Method and Recent Variants

Accelerated optimization algorithms and ordinary differential equations: the convex non Euclidean case

Multigoal-oriented optimal control problems with nonlinear PDE constraints

The Modified MSA, a Gradient Flow and Convergence