Journal
JOURNAL OF COMPUTATIONAL PHYSICS
Volume 255, Issue -, Pages 699-719Publisher
ACADEMIC PRESS INC ELSEVIER SCIENCE
DOI: 10.1016/j.jcp.2013.07.017
Keywords
Kinetic equations; Discrete velocity models; Semi-Lagrangian schemes; Boltzmann-BGK equation; Euler solver; High order scheme
Ask authors/readers for more resources
In a recent paper we presented a new ultra efficient numerical method for solving kinetic equations of the Boltzmann type (Dimarco and Loubere, 2013) [17]. The key idea, on which the method relies, is to solve the collision part on a grid and then to solve exactly the transport part by following the characteristics backward in time. On the contrary to classical semi-Lagrangian methods one does not need to reconstruct the distribution function at each time step. This allows to tremendously reduce the computational cost and to perform efficient numerical simulations of kinetic equations up to the six-dimensional case without parallelization. However, the main drawback of the method developed was the loss of spatial accuracy close to the fluid limit. In the present work, we modify the scheme in such a way that it is able to preserve the high order spatial accuracy for extremely rarefied and fluid regimes. In particular, in the fluid limit, the method automatically degenerates into a high order method for the compressible Euler equations. Numerical examples are presented which validate the method, show the higher accuracy with respect to the previous approach and measure its efficiency with respect to well-known schemes (Direct Simulation Monte Carlo, Finite Volume, MUSCL, WENO). (C) 2013 Elsevier Inc. All rights reserved.
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