6533b82bfe1ef96bd128df4b

RESEARCH PRODUCT

Asymptotic preserving IMEX finite volume schemes for low Mach number Euler equations with gravitation

Leonid YelashGeorgij BispenMria Lukov-medvid'ov

subject

Numerical AnalysisFinite volume methodPhysics and Astronomy (miscellaneous)DiscretizationApplied MathematicsMathematical analysis010103 numerical & computational mathematics01 natural sciencesComputer Science ApplicationsEuler equations010101 applied mathematicsLinear mapComputational Mathematicssymbols.namesakeNonlinear systemMach numberModeling and SimulationStability theorysymbolsCompressibility0101 mathematicsMathematics

description

In this paper we will present and analyze a new class of the IMEX finite volume schemes for the Euler equations with a gravity source term. We will in particular concentrate on a singular limit of weakly compressible flows when the Mach number M1. In order to efficiently resolve slow dynamics we split the whole nonlinear system in a stiff linear part governing the acoustic and gravity waves and a non-stiff nonlinear part that models nonlinear advection effects. For time discretization we use a special class of the so-called globally stiffly accurate IMEX schemes and approximate the stiff linear operator implicitly and the non-stiff nonlinear operator explicitly. For spatial discretization the finite volume approximation is used with the central and Rusanov/LaxFriedrichs numerical fluxes for the linear and nonlinear subsystem, respectively. In the case of a constant background potential temperature we prove theoretically that the method is asymptotically consistent and asymptotically stable uniformly with respect to small Mach number. We also analyze experimentally convergence rates in the singular limit when the Mach number tends to zero.

https://doi.org/10.1016/j.jcp.2017.01.020