Andrea Thomann

Andrea Thomann, Markus Zenk, Christian Klingenberg

We present a well-balanced finite volume solver for the compressible Euler equations with gravity where the approximate Riemann solver is derived using a relaxation approach. Besides the well-balanced property, the scheme is robust with respect to the physical admissible states. Another feature of the method is that it can maintain general stationary solutions of the hydrostatic equilibrium up to machine precision. For the first order scheme we present a well-balanced and positivity preserving second order extension using a modified minmod slope limiter. To maintain the well-balanced property, we reconstruct in equilibrium variables. Numerical examples are performed to demonstrate the accuracy, well-balanced and positivity preserving property of the presented scheme for up to 3 space dimensions.

Michael Dumbser, Andrea Thomann, Maurizio Tavelli et al.

We introduce a novel structure-preserving vertex-staggered semi-implicit four-split discretization of a unified first order hyperbolic formulation of continuum mechanics that is able to describe at the same time fluid and solid materials within the same mathematical model. The governing PDE system goes back to pioneering work of Godunov, Romenski, Peshkov and collaborators. Previous structure-preserving discretizations of this system allowed to respect the curl-free properties of the distortion field and the specific thermal impulse in the absence of source terms and were consistent with the low Mach number limit with respect to the adiabatic sound speed. However, the evolution of the thermal impulse and the distortion field were still discretized explicitly, thus requiring a rather severe CFL stability restriction on the time step based on the shear sound speed and the finite, but potentially large, speed of heat waves. Instead, the new four-split semi-implicit scheme presented in this paper has a material time step restriction only. For this purpose, the governing PDE system is split into four subsystems: i) a convective subsystem, which is the only one that is treated explicitly; ii) a heat subsystem, iii) a subsystem containing momentum, distortion field and specific thermal impulse; iv) a pressure subsystem. The three subsystems ii)-iv) are all discretized implicitly, hence a rather mild CFL restriction based on the velocity of the continuum is imposed. The method is asymptotically consistent with the low Mach number limit and the stiff relaxation limits. Moreover, it maintains an exactly curl-free distortion field and thermal impulse in the case of linear source terms or in their absence. The scheme is benchmarked against classical test cases verifying its theoretical properties.