Implicit shock tracking using an optimization-based high-order discontinuous Galerkin method

Abstract

A novel framework for resolving discontinuous solutions of conservation laws, e.g., contact lines, shock waves, and interfaces, using implicit tracking and a high-order discontinuous Galerkin (DG) discretization was introduced in [38]. Central to the framework is an optimization problem whose solution is a discontinuity-aligned mesh and the corresponding high-order approximation to the flow that does not require explicit meshing of the unknown discontinuity surface. The method was shown to deliver highly accurate solutions on coarse, high-order discretizations without nonlinear stabilization and recover optimal convergence rates even for problems with discontinuous solutions. This work extends the implicit tracking framework such that robustness is improved and convergence accelerated. In particular, we introduce an improved formulation of the central optimization problem and an associated sequential quadratic programming (SQP) solver. The new error-based objective function penalizes violation of the DG residual in an enriched test space and is shown to have excellent tracking properties. The SQP solver simultaneously converges the nodal coordinates of the mesh and DG solution to their optimal values and is equipped with a number of features to ensure robust, fast convergence: Levenberg-Marquardt approximation of the Hessian with weighted elliptic regularization, backtracking line search, and rigorous convergence criteria. We use the proposed method to solve a range of inviscid conservation laws of varying difficulty. We show the method is able to deliver accurate solutions on coarse, high-order meshes and the SQP solver is robust and usually able to drive the first-order optimality system to tight tolerances.

arXiv comment35 pages, 20 figures

References in corpus (1)

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

Cited by in corpus (18)

• A robust, high-order implicit shock tracking method for simulation of complex, high-speed flows
• The Moving Discontinuous Galerkin Finite Element Method with Interface Condition Enforcement for Compressible Viscous Flows
• Registration-based model reduction of parameterized two-dimensional conservation laws
• Extrapolated Shock Tracking: Bridging shock-fitting and embedded boundary methods
• A Moving Discontinuous Galerkin Finite Element Method with Interface Conservation Enforcement for Compressible Flows
• A Least-Squares Formulation of the Moving Discontinuous Galerkin Finite Element Method with Interface Condition Enforcement
• Implicit shock tracking for unsteady flows by the method of lines
• Preserving general physical properties in model reduction of dynamical systems via constrained-optimization projection
• Extrapolated DIscontinuity Tracking for complex 2D shock interactions
• High-Order Mesh Morphing for Boundary and Interface Fitting to Implicit Geometries
• A globally convergent method to accelerate large-scale optimization using on-the-fly model hyperreduction: application to shape optimization
• Combining high-order metric interpolation and geometry implicitization for curved $r$-adaption
• Space-time registration-based model reduction of parameterized one-dimensional hyperbolic PDEs
• Defining metric-aware size-shape measures to validate and optimize curved high-order meshes
• Registration-based model reduction in complex two-dimensional geometries
• A globalized and preconditioned Newton-CG solver for metric-aware curved high-order mesh optimization
• A coupled discontinuous Galerkin-Finite Volume framework for solving gas dynamics over embedded geometries
• UnDiFi-2D: an Unstructured Discontinuity Fitting code for 2D grids