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.

Learning data-driven discretizations for partial differential equations

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

Linking Machine Learning with Multiscale Numerics: Data-Driven Discovery of Homogenized Equations

DL-PDE: Deep-learning based data-driven discovery of partial differential equations from discrete and noisy data

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

Robust data-driven discovery of governing physical laws with error bars

Promising directions of machine learning for partial differential equations

Solving Non-linear Kolmogorov Equations in Large Dimensions by Using Deep Learning: A Numerical Comparison of Discretization Schemes

Physics-informed learning of governing equations from scarce data

A step towards holistic discretisation of stochastic partial differential equations

Machine Discovery of Partial Differential Equations from Spatiotemporal Data: A Sparse Bayesian Learning Framework.

Machine Learning for Partial Differential Equations

Learning to Discretize: Solving 1D Scalar Conservation Laws via Deep Reinforcement Learning

A Decision-Making Machine Learning Approach in Hermite Spectral Approximations of Partial Differential Equations

Data-Driven Finite Elements Methods: Machine Learning Acceleration of Goal-Oriented Computations

Deep Hidden Physics Models: Deep Learning of Nonlinear Partial Differential Equations

Solving Partial Differential Equations Using Deep Learning and Physical Constraints

Machine Learning of Partial Differential Equations from Noise Data

Recovering missing CFD data for high-order discretizations using deep neural networks and dynamics learning

Physics-Guided Discovery of Highly Nonlinear Parametric Partial Differential Equations

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