Accelerated simulation of Boltzmann-BGK equations near the diffusive limit with asymptotic-preserving multilevel Monte Carlo

by   Emil Løvbak, et al.

Kinetic equations model the position-velocity distribution of particles subject to transport and collision effects. Under a diffusive scaling, these combined effects converge to a diffusion equation for the position density in the limit of an infinite collision rate. Despite this well-defined limit, numerical simulation is expensive when the collision rate is high but finite, as small time steps are then required. In this work, we present an asymptotic-preserving multilevel Monte Carlo particle scheme that makes use of this diffusive limit to accelerate computations. In this scheme, we first sample the diffusive limiting model to compute a biased initial estimate of a Quantity of Interest, using large time steps. We then perform a limited number of finer simulations with transport and collision dynamics to correct the bias. The efficiency of the multilevel method depends on being able to perform correlated simulations of particles on a hierarchy of discretization levels. We present a method for correlating particle trajectories and present both an analysis and numerical experiments. We demonstrate that our approach significantly reduces the cost of particle simulations in high-collisional regimes, compared with prior work, indicating significant potential for adopting these schemes in various areas of active research.


page 1

page 2

page 3

page 4


A multilevel Monte Carlo method for asymptotic-preserving particle schemes

Kinetic equations model distributions of particles in position-velocity ...

Multilevel Asymptotic-Preserving Monte Carlo for Particle Simulations

We develop a novel Multilevel Asymptotic-Preserving Monte Carlo (ML-APMC...

Particle simulation methods for the Landau-Fokker-Planck equation with uncertain data

The design of particle simulation methods for collisional plasma physics...

A fast Chebyshev method for the Bingham closure with application to active nematic suspensions

Continuum kinetic theories provide an important tool for the analysis an...

Performance of the BGSDC integrator for computing fast ion trajectories in nuclear fusion reactors

Modelling neutral beam injection (NBI) in fusion reactors requires compu...

A cell-based population control of Monte Carlo particles for the global variance reduction for transport equations

We present a population control method with sampling and regulation step...

Monte Carlo simulations of interacting particles with fast and accurate electrostatics

To minimise systematic errors in Monte Carlo simulations of charged part...

Please sign up or login with your details

Forgot password? Click here to reset