Asymptotic-preserving Projective Integration Schemes for Kinetic Equations in the Diffusion Limit

Abstract

We investigate a projective integration scheme for a kinetic equation in the limit of vanishing mean free path, in which the kinetic description approaches a diffusion phenomenon. The scheme first takes a few small steps with a simple, explicit method, such as a spatial centered flux/forward Euler time integration, and subsequently projects the results forward in time over a large time step on the diffusion time scale. We show that, with an appropriate choice of the inner step size, the time-step restriction on the outer time step is similar to the stability condition for the diffusion equation, whereas the required number of inner steps does not depend on the mean free path. We also provide a consistency result. The presented method is asymptotic-preserving, in the sense that the method converges to a standard finite volume scheme for the diffusion equation in the limit of vanishing mean free path. The analysis is illustrated with numerical results, and we present an application to the Su-Olson test.

Projective Integration for Nonlinear BGK Kinetic Equations

“…However, as ε goes to 0, we obtain the limiting system (1.5) for which a standard finite volume/forward Euler method only needs to satisfy a stability restriction of the form ∆t ≤ C∆ x, with C a constant that depends on the specific choice of the scheme. In [8], it was proposed to use a projective integration method to accelerate such a brute-force integration; the idea, originating from [5], is the following. Starting from a computed numerical solution f n at time t n = n∆t, one first takes K + 1 inner steps of size δt using (1.7), denoted as f n,k+1 , in which the superscripts (n, k) denote the numerical solution at time t n,k = n∆t + kδt.…”

Erratum: A High-Order Asymptotic-Preserving Scheme for Kinetic Equations Using Projective Integration

“…see the bottom of page 22 in [1], and compare with the consistency derivation in [2]. However, equation (5) is valid only if one has…”

A Collision-Based Hybrid Method for Time-Dependent, Linear, Kinetic Transport Equations

“…From (24), is is clear that Integrating the equations in (26) over S 2 , adding them together, and applying the results of (27) and (28) gives…”