Marco Artiano

Marco Artiano, Hendrik Ranocha, Saurav Samantaray

We consider the compressible Euler system with anelastic scaling, modeling isentropic flows under the influence of gravity. In the zero-Mach-number limit, the solution of the compressible Euler system converges to a variable density anelastic incompressible limit system. In this work, we present the design and analysis of a class of higher-order linearly implicit IMEX Runge-Kutta schemes that are asymptotic preserving, i.e., they respect the transitory nature of the governing equations in the limit. The presence of gravitational potential warrants the incorporation of the well-balancing property. The scheme is developed as a novel combination of a penalization of a linear steady state, a finite-volume balance-preserving reconstruction, and a source term discretization preserving steady states. The penalization plays a crucial role in obtaining a linearly implicit scheme, and well-balanced flux-source discretization ensures accuracy in very low Mach number regimes. Some results of numerical case studies are presented to corroborate the theoretical assertions.

Marco Artiano, Arpit Babbar, Michael Schlottke-Lakemper et al.

Jin-Xin relaxation is a method for approximating non-linear hyperbolic conservation laws by a linear system of hyperbolic equations with an $\varepsilon$ dependent stiff source term. The system formally relaxes to the original conservation law as $\varepsilon \to 0$. An asymptotic analysis of the Jin-Xin relaxation system shows that it can be seen as a convection-diffusion equation with a diffusion coefficient that depends on the relaxation parameter $\varepsilon$. This work makes use of this property to use the Jin-Xin relaxation system as a shock-capturing method for high-order discontinuous Galerkin (DG) or flux reconstruction (FR) schemes. The idea is to use a smoothness indicator to choose the $\varepsilon$ value in each cell, so that we can use larger $\varepsilon$ values in non-smooth regions to add extra numerical dissipation. We show how this can be done by using a single stage method by using the compact Runge-Kutta FR method that handles the stiff source term by using IMplicit-EXplicit Runge-Kutta (IMEX-RK) schemes. Numerical results involving Burgers' equation and the compressible Euler equations are shown to demonstrate the effectiveness of the proposed method.

Marco Artiano, Hendrik Ranocha

Many entropy-conservative and entropy-stable (summarized as entropy-preserving) methods for hyperbolic conservation laws rely on Tadmor's theory for two-point entropy-preserving numerical fluxes and its higher-order extension via flux differencing using summation-by-parts (SBP) operators, e.g., in discontinuous Galerkin spectral element methods (DGSEMs). The underlying two-point formulations have been extended to nonconservative systems using fluctuations by Castro et al. (2013, doi:10.1137/110845379) with follow-up generalizations to SBP methods. We propose specific forms of entropy-preserving fluctuations for nonconservative hyperbolic systems that are simple to interpret and allow an algorithmic construction of entropy-preserving methods. We analyze necessary and sufficient conditions, and obtain a full characterization of entropy-preserving three-point methods within the finite volume framework. This formulation is extended to SBP methods in multiple space dimensions on Cartesian and curvilinear meshes. Additional properties such as well-balancedness extend naturally from the underlying finite volume method to the SBP framework. We use the algorithmic construction enabled by the chosen formulation to derive several new entropy-preserving schemes for nonconservative hyperbolic systems, e.g., the compressible Euler equations of an ideal gas using the internal energy equation and a dispersive shallow-water model. Numerical experiments show the robustness and accuracy of the proposed schemes.