|
Special
Applied Math Seminar |
|
A new family of Monte Carlo schemes is introduced for the numerical solution of the Boltzmann equation of rarefied gas dynamics. The schemes are inspired by the Wild sum expansion of the solution of the Boltzmann equation for Maxwellian molecules, and consist in a novel time discretization of the equation. In particular, high order terms in the expansion are replaced by the equilibrium Maxwellian distribution. The two main features of the schemes are high order accuracy in time and asymptotic preservation. The first property allows to recover accurate solutions with time steps larger than those required by Direct Simulation Monte Carlo (more precisely by the Nanbu-Babowsky scheme) while the latter guarantees that for vanishing Knudsen number, the numerical solution relaxes to the local Maxwellian. Conservation of mass, momentum and energy are preserved by the scheme. Numerical results on several space homogeneous problems show the improvement of the new schemes over standard DSMC. Applications to a one-dimensional shock wave problem are also presented. |