Implicit shock tracking for unsteady flows by the method of lines
From MaRDI portal
Publication:6358779
DOI10.1016/J.JCP.2021.110906arXiv2101.08913MaRDI QIDQ6358779FDOQ6358779
Authors: Andrew Shi, Per-Olof Persson, Matthew J. Zahr 
Publication date: 21 January 2021
Abstract: 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 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.
Hyperbolic equations and hyperbolic systems (35Lxx) Numerical methods for partial differential equations, initial value and time-dependent initial-boundary value problems (65Mxx) Basic methods in fluid mechanics (76Mxx)
This page was built for publication: Implicit shock tracking for unsteady flows by the method of lines
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6358779)
