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

by   Matthew J. Zahr, et al.

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.



There are no comments yet.


page 18

page 20

page 22

page 24

page 27

page 29

page 32


Implicit shock tracking for unsteady flows by the method of lines

A recently developed high-order implicit shock tracking (HOIST) framewor...

A robust, high-order implicit shock tracking method for simulation of complex, high-speed flows

High-order implicit shock tracking is a new class of numerical methods t...

A Discontinuous Galerkin method for Shock Capturing using a mixed high-order and sub-grid low-order approximation space

This article considers a new discretization scheme for conservation laws...

SANM: A Symbolic Asymptotic Numerical Solver with Applications in Mesh Deformation

Solving nonlinear systems is an important problem. Numerical continuatio...

A comparison of interpolation techniques for non-conformal high-order discontinuous Galerkin methods

The capability to incorporate moving geometric features within models fo...

Towards a Scalable Hierarchical High-order CFD Solver

Development of highly scalable and robust algorithms for large-scale CFD...

Naturally curved quadrilateral mesh generation using an adaptive spectral element solver

We describe an adaptive version of a method for generating valid natural...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.