Implicit shock tracking for unsteady flows by the method of lines

Journal Article · Mon Jan 10 00:00:00 EST 2022 · Journal of Computational Physics

DOI:https://doi.org/10.1016/j.jcp.2021.110906· OSTI ID:1883333

Shi, A. ^[1]; Persson, P.-O. ^[2]; Zahr, M. J. ^[3]

---

1. Univ. of California, Berkeley, CA (United States)
2. Univ. of California, Berkeley, CA (United States); Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States)
3. Univ. of Notre Dame, IN (United States)

Here, a recently developed high-order implicit shock tracking (HOIST) framework for resolving discontinuous solutions of inviscid, steady conservation laws 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 rapid 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.

• Research Organization: Lawrence Berkeley National Laboratory (LBNL), Berkeley, CA (United States)
• Sponsoring Organization: US Air Force Office of Scientific Research (AFOSR); USDOE Office of Science (SC), Advanced Scientific Computing Research (ASCR); USDOE Office of Science (SC), High Energy Physics (HEP)
• Grant/Contract Number: AC02-05CH11231
• OSTI ID: 1883333
• Alternate ID(s): OSTI ID: 1841461
• Journal Information: Journal of Computational Physics, Journal Name: Journal of Computational Physics Vol. 454; ISSN 0021-9991
• Publisher: Elsevier
• Country of Publication: United States
• Language: English

Related Subjects

71 CLASSICAL AND QUANTUM MECHANICS
GENERAL PHYSICS
Discontinuous Galerkin
High-order methods
High-speed flows
Method of lines
Shock fitting
Shock tracking