Implicit shock tracking for unsteady flows by the method of lines

by   Andrew Shi, et al.

A recently developed high-order implicit shock tracking (HOIST) framework for resolving discontinuous solutions of inviscid, steady conservation laws [41, 43] is extended to the unsteady case. Central to the framework is an optimization problem which simultaneously computes a discontinuity-aligned mesh and the corresponding high-order approximation to the flow, which provides nonlinear stabilization and a high-order approximation to the solution. This work extends the implicit shock tracking framework to the case of unsteady conservation laws using a method of lines discretization via a diagonally implicit Runge-Kutta method by "solving a steady problem at each timestep". We formulate and solve an optimization problem that produces a feature-aligned mesh and solution at each Runge-Kutta stage of each timestep, and advance this solution in time by standard Runge-Kutta update formulas. A Rankine-Hugoniot based prediction of the shock location together with a high-order, untangling mesh smoothing procedure provides a high-quality initial guess for the optimization problem at each time, which results in Newton-like convergence of the sequential quadratic programing (SQP) optimization solver. This method is shown to deliver highly accurate solutions on coarse, high-order discretizations without nonlinear stabilization and recover the design accuracy of the Runge-Kutta scheme. We demonstrate this framework on a series of inviscid, unsteady conservation laws in both one- and two- dimensions. We also verify that our method is able to recover the design order of accuracy of our time integrator in the presence of a strong discontinuity.



There are no comments yet.


page 17

page 21


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

A novel framework for resolving discontinuous solutions of conservation ...

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...

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

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

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...

An order-adaptive compact approximation Taylor method for systems of conservation laws

We present a new family of high-order shock-capturing finite difference ...

High Order Residual Distribution Conservative Finite Difference HWENO Scheme for Steady State Problems

In this paper, we develop a high order residual distribution (RD) method...

A Hybrid High-Order method for incompressible flows of non-Newtonian fluids with power-like convective behaviour

In this work, we design and analyze a Hybrid High-Order (HHO) discretiza...
This week in AI

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