Total-variation-diminishing implicit–explicit Runge–Kutta methods for the simulation of double-diffusive convection in astrophysics

• Published in: Journal of computational physics Vol. 231; no. 9; pp. 3561 - 3586
• Main Authors: Kupka, Friedrich, Happenhofer, Natalie, Higueras, Inmaculada, Koch, Othmar
• Format: Journal Article
• Language: English
• Published: Kidlington Elsevier Inc 01-05-2012; Elsevier
• Subjects: Computational techniques; Double-diffusive convection; Exact sciences and technology; Hydrodynamics; Mathematical methods in physics; Numerical methods; Physics; SSP; Stellar convection and pulsation; Strong stability preserving; Total-variation-diminishing; TVD; Total-variation-diminishing; SSP; TVD; Stellar convection and pulsation; Numerical methods; Strong stability preserving; Hydrodynamics; Double-diffusive convection; Total energy; Numerical method; Equations of motion; Stellar pulsations; Euler scheme; Convection; Calculation methods; Thermodynamics; Equations of state; Discretization; Energy equation; Calculation; Runge-Kutta methods; Navier-Stokes equations; Finite difference method
• ISSN: 0021-9991; 1090-2716
• DOI: 10.1016/j.jcp.2011.12.031

We put forward the use of total-variation-diminishing (or more generally, strong stability preserving) implicit–explicit Runge–Kutta methods for the time integration of the equations of motion associated with the semiconvection problem in the simulation of stellar convection. The fully compressible Navier–Stokes equation, augmented by continuity and total energy equations, and an equation of state describing the relation between the thermodynamic quantities, is semi-discretized in space by essentially non-oscillatory schemes and dissipative finite difference methods. It is subsequently integrated in time by Runge–Kutta methods which are constructed such as to preserve the total variation diminishing (or strong stability) property satisfied by the spatial discretization coupled with the forward Euler method. We analyse the stability, accuracy and dissipativity of the time integrators and demonstrate that the most successful methods yield a substantial gain in computational efficiency as compared to classical explicit Runge–Kutta methods.