Time Relaxed Monte Carlo Methods for the Boltzmann Equation

Authors:
Lorenzo Pareschi;Giovanni Russo
Affiliations:
-;-
Venue:
SIAM Journal on Scientific Computing
Year:
2001

Citing 0
Cited 10

High order numerical methods for the space non-homogeneous Boltzmann equation

Journal of Computational Physics
Uniformly stable numerical schemes for the Boltzmann equation preserving the compressible Navier-Stokes asymptotics

Journal of Computational Physics
Gridless DSMC

Journal of Computational Physics
A class of asymptotic-preserving schemes for kinetic equations and related problems with stiff sources

Journal of Computational Physics
Fluid Solver Independent Hybrid Methods for Multiscale Kinetic Equations

SIAM Journal on Scientific Computing
On DSMC Calculations of Rarefied Gas Flows with Small Number of Particles in Cells

SIAM Journal on Scientific Computing
Exponential Runge-Kutta Methods for Stiff Kinetic Equations

SIAM Journal on Numerical Analysis
A high-order Monte Carlo algorithm for the direct simulation of Boltzmann equation

Journal of Computational Physics
Locally refined discrete velocity grids for stationary rarefied flow simulations

Journal of Computational Physics
An efficient particle Fokker-Planck algorithm for rarefied gas flows

Journal of Computational Physics

Quantified Score

Hi-index	0.04

Visualization

Abstract

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 of 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 (DSMC), while the latter guarantees that for the 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.