WIAS Preprint No. 1320, (2008)

Deviational particle Monte Carlo for the Boltzmann equation



Authors

  • Wagner, Wolfgang

2010 Mathematics Subject Classification

  • 65R20 82C40

2008 Physics and Astronomy Classification Scheme

  • 05.10.Ln

Keywords

  • Monte Carlo method, Boltzmann equation, reference Maxwellian, deviational particles

Abstract

The paper describes the deviational particle Monte Carlo method for the Boltzmann equation. The approach is an application of the general ``control variates'' variance reduction technique to the problem of solving a nonlinear equation. The deviation of the solution from a reference Maxwellian is approximated by a system of positive and negative particles. Previous results from the literature are modified and extended. New algorithms are proposed that cover the nonlinear Boltzmann equation (instead of a linearized version) with a general interaction model (instead of hard spheres). The algorithms are obtained as procedures for generating trajectories of Markov jump processes. This provides the framework for deriving the limiting equations, when the number of particles tends to infinity. These equations reflect the influence of various numerical approximation parameters. Detailed simulation schemes are provided for the variable hard sphere interaction model.

Appeared in

  • Monte Carlo Methods Appl., 14 (2008) pp. 191--268.

Download Documents