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

page 17

page 21

research
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 ...
research
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...
research
04/22/2023

High-order implicit shock tracking boundary conditions for flows with parametrized shocks

High-order implicit shock tracking (fitting) is a class of high-order, o...
research
05/18/2021

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

Solving nonlinear systems is an important problem. Numerical continuatio...
research
08/08/2023

A space-time high-order implicit shock tracking method for shock-dominated unsteady flows

High-order implicit shock tracking (fitting) is a class of high-order, o...
research
02/28/2022

Arbitrarily high order implicit ODE integration by correcting a neural network approximation with Newton's method

As a method of universal approximation deep neural networks (DNNs) are c...
research
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...

Please sign up or login with your details

Forgot password? Click here to reset