A High-Order Asymptotic-Preserving Scheme for Kinetic Equations Using Projective Integration
Author(s) -
Pauline Lafitte,
Annelies Lejon,
Giovanni Samaey
Publication year - 2016
Publication title -
siam journal on numerical analysis
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 2.78
H-Index - 134
eISSN - 1095-7170
pISSN - 0036-1429
DOI - 10.1137/140966708
Subject(s) - mathematics , relaxation (psychology) , mathematical analysis , stability (learning theory) , euler equations , time derivative , euler method , scaling , simple (philosophy) , geometry , psychology , social psychology , philosophy , epistemology , machine learning , computer science
We investigate a high-order, fully explicit, asymptotic-preserving scheme for a kinetic equation with linear relaxation, both in the hydrodynamic and diffusive scalings in which a hyperbolic, resp., parabolic, limiting equation exists. The scheme first takes a few small (inner) steps with a simple, explicit method (such as direct forward Euler) to damp out the stiff components of the solution and estimate the time derivative of the slow components. These estimated time derivatives are then used in an (outer) Runge--Kutta method of arbitrary order. We show that, with an appropriate choice of inner step size, the time-step restriction on the outer time step is similar to the stability condition for the limiting macroscopic equation. Moreover, the number of inner time steps is also independent of the scaling parameter. We analyze stability and consistency, and illustrate with numerical results.
Accelerating Research
Robert Robinson Avenue,
Oxford Science Park, Oxford
OX4 4GP, United Kingdom
Address
John Eccles HouseRobert Robinson Avenue,
Oxford Science Park, Oxford
OX4 4GP, United Kingdom