A multilevel Monte Carlo method for asymptotic-preserving particle schemes

by   Emil Løvbak, et al.

Kinetic equations model distributions of particles in position-velocity phase space. Often, one is interested in studying the long-time behavior of particles in the diffusive limit, in which the collision rate tends to infinity. Classical particle-based techniques suffer from a strict time-step restriction in this limit, to maintain stability. Asymptotic-preserving schemes avoid this problem, but introduce an additional time discretization error, possibly resulting in an unacceptably large bias for larger time steps. Here, we present and analyze a multilevel Monte Carlo scheme that reduces this bias by combining estimates using a hierarchy of different time step sizes. We demonstrate how to correlate trajectories from this scheme, using different time steps. We also present a strategy for selecting the levels in the multilevel scheme. Our approach significantly reduces the computation required to perform accurate simulations of the considered kinetic equations, compared to classical Monte Carlo approaches.



page 26


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

Kinetic equations model the position-velocity distribution of particles ...

Multilevel Asymptotic-Preserving Monte Carlo for Particle Simulations

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

Kinetic-diffusion asymptotic-preserving Monte Carlo algorithm for Boltzmann-BGK in the diffusive scaling

We develop a novel Monte Carlo strategy for the simulation of the Boltzm...

Accelerating Surface Tension Calculation in SPH via Particle Classification Monte Carlo Integration

Surface tension has a strong influence on the shape of fluid interfaces....

Monte Carlo simulations of interacting particles with fast and accurate electrostatics

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

Accelerating the estimation of energetic particle confinement statistics in stellarators using multifidelity Monte Carlo

In the design of stellarators, energetic particle confinement is a criti...

The Seven-League Scheme: Deep learning for large time step Monte Carlo simulations of stochastic differential equations

We propose an accurate data-driven numerical scheme to solve Stochastic ...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.