Implicit shock tracking for unsteady flows by the method of lines

01/22/2021
by   Andrew Shi, et al.
0

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.

READ FULL TEXT
POST COMMENT

Comments

There are no comments yet.

Authors

page 17

page 21

12/24/2019

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

A novel framework for resolving discontinuous solutions of conservation ...
05/01/2021

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...
05/18/2021

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

Solving nonlinear systems is an important problem. Numerical continuatio...
11/04/2019

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...
07/02/2020

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

We present a new family of high-order shock-capturing finite difference ...
04/14/2021

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...
06/28/2021

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.