Journal
SIAM JOURNAL ON SCIENTIFIC COMPUTING
Volume 31, Issue 1, Pages 334-368Publisher
SIAM PUBLICATIONS
DOI: 10.1137/07069479X
Keywords
transport equations; diffusion limit; asymptotic preserving schemes; stiff terms
Categories
Ask authors/readers for more resources
We propose a new numerical scheme for linear transport equations. It is based on a decomposition of the distribution function into equilibrium and nonequilibrium parts. We also use a projection technique that allows us to reformulate the kinetic equation into a coupled system of an evolution equation for the macroscopic density and a kinetic equation for the nonequilibrium part. By using a suitable time semi-implicit discretization, our scheme is able to accurately approximate the solution in both kinetic and diffusion regimes. It is asymptotic preserving in the following sense: when the mean free path of the particles is small, our scheme is asymptotically equivalent to a standard numerical scheme for the limit diffusion model. A uniform stability property is proved for the simple telegraph model. Various boundary conditions are studied. Our method is validated in one-dimensional cases by several numerical tests and comparisons with previous asymptotic preserving schemes.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available