Quiet direct simulation Monte-Carlo with random timesteps

  • Authors:
  • William Peter

  • Affiliations:
  • Johns Hopkins University, Applied Physics Laboratory, 11100 Johns Hopkins Road, Laurel, MD 20723-6099, United States

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2007

Quantified Score

Hi-index 31.46

Visualization

Abstract

Use of a high-order deterministic sampling technique in direct simulation Monte-Carlo (DSMC) simulations eliminates statistical noise and improves computational performance by orders of magnitude. In this paper it is also shown that if a random timestep is used in place of a fixed timestep, there is an additional improvement in performance. This performance can be increased by using a timestep that samples a random variable with a high-kurtosis probability density function. As a simple example of the method, the one-dimensional diffusion equation with an exponentially-distributed timestep is simulated, and a performance gain of approximately two is obtained. Applications to numerical simulations of fluids and plasmas are indicated.