385 research outputs found
Towards an ultra efficient kinetic scheme. Part I: basics on the BGK equation
In this paper we present a new ultra efficient numerical method for solving
kinetic equations. In this preliminary work, we present the scheme in the case
of the BGK relaxation operator. The scheme, being based on a splitting
technique between transport and collision, can be easily extended to other
collisional operators as the Boltzmann collision integral or to other kinetic
equations such as the Vlasov equation. The key idea, on which the method
relies, is to solve the collision part on a grid and then to solve exactly the
transport linear part by following the characteristics backward in time. The
main difference between the method proposed and semi-Lagrangian methods is that
here we do not need to reconstruct the distribution function at each time step.
This allows to tremendously reduce the computational cost of the method and it
permits for the first time, to the author's knowledge, to compute solutions of
full six dimensional kinetic equations on a single processor laptop machine.
Numerical examples, up to the full three dimensional case, are presented which
validate the method and assess its efficiency in 1D, 2D and 3D
Cauchy problem for the Boltzmann-BGK model near a global Maxwellian
In this paper, we are interested in the Cauchy problem for the Boltzmann-BGK
model for a general class of collision frequencies. We prove that the
Boltzmann-BGK model linearized around a global Maxwellian admits a unique
global smooth solution if the initial perturbation is sufficiently small in a
high order energy norm. We also establish an asymptotic decay estimate and
uniform -stability for nonlinear perturbations.Comment: 26 page
An efficient numerical method for solving the Boltzmann equation in multidimensions
International audienceIn this paper we deal with the extension of the Fast Kinetic Scheme (FKS) [J. Comput. Phys., Vol. 255, 2013, pp 680-698] originally constructed for solving the BGK equation, to the more challenging case of the Boltzmann equation. The scheme combines a robust and fast method for treating the transport part based on an innovative Lagrangian technique supplemented with fast spectral schemes to treat the collisional operator by means of an operator splitting approach. This approach along with several implementation features related to the parallelization of the algorithm permits to construct an efficient simulation tool which is numerically tested against exact and reference solutions on classical problems arising in rarefied gas dynamic. We present results up to the 3DĂ—3D case for unsteady flows for the Variable Hard Sphere model which may serve as benchmark for future comparisons between different numerical methods for solving the multidimensional Boltzmann equation. For this reason, we also provide for each problem studied details on the computational cost and memory consumption as well as comparisons with the BGK model or the limit model of compressible Euler equations
A particle micro-macro decomposition based numerical scheme for collisional kinetic equations in the diffusion scaling
In this work, we derive particle schemes, based on micro-macro decomposition,
for linear kinetic equations in the diffusion limit. Due to the particle
approximation of the micro part, a splitting between the transport and the
collision part has to be performed, and the stiffness of both these two parts
prevent from uniform stability. To overcome this difficulty, the micro-macro
system is reformulated into a continuous PDE whose coefficients are no longer
stiff, and depend on the time step in a consistent way. This
non-stiff reformulation of the micro-macro system allows the use of standard
particle approximations for the transport part, and extends the work in [5]
where a particle approximation has been applied using a micro-macro
decomposition on kinetic equations in the fluid scaling. Beyond the so-called
asymptotic-preserving property which is satisfied by our schemes, they
significantly reduce the inherent noise of traditional particle methods, and
they have a computational cost which decreases as the system approaches the
diffusion limit
- …