Total-variation-diminishing implicit-explicit Runge-Kutta methods for the simulation of double-diffusive convection in astrophysics
Abstract
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.
- Publication:
-
Journal of Computational Physics
- Pub Date:
- May 2012
- DOI:
- 10.1016/j.jcp.2011.12.031
- arXiv:
- arXiv:1106.5672
- Bibcode:
- 2012JCoPh.231.3561K
- Keywords:
-
- Mathematics - Numerical Analysis;
- 65M06;
- 65M08;
- 65M20;
- 65L05;
- 76F65
- E-Print:
- 45 pages, 17 figures, 12 tables, accepted for publication in Journal of Computational Physics