On a Class of Implicit–Explicit Runge–Kutta Schemes for Stiff Kinetic Equations Preserving the Navier–Stokes Limit

Abstract

Implicit–explicit (IMEX) Runge–Kutta (RK) schemes are popular high order time discretization methods for solving stiff kinetic equations. As opposed to the compressible Euler limit (leading order asymptotics of the Boltzmann equation as the Knudsen number $$\varepsilon $$ goes to zero), their asymptotic behavior at the Navier–Stokes (NS) level (next order asymptotics) was rarely studied. In this paper, we analyze a class of existing IMEX RK schemes and show that, under suitable initial conditions, they can capture the NS limit without resolving the small parameter $$\varepsilon $$ , i.e., $$\varepsilon =o(\Delta t)$$ , $$\Delta t^m=o(\varepsilon )$$ , where m is the order of the explicit RK part in the IMEX scheme. Extensive numerical tests for BGK and ES-BGK models are performed to verify our theoretical results.