Multilevel Asymptotic-Preserving Monte Carlo for Particle Simulations

by   Bert Mortier, et al.

We develop a novel Multilevel Asymptotic-Preserving Monte Carlo (ML-APMC) method for simulating the kinetic Boltzmann transport equation with Bhatnagar-Gross-Krook (BGK) collision operator. This equation occurs, for instance, in mathematical models of the neutral particles in the plasma edge of nuclear fusion reactors. The main features of our method are a new and improved recipe for correlating particle trajectories with different time step sizes, and a new and more general level selection strategy. We illustrate the efficiency of our ML-APMC method by applying it to a one-dimensional fusion test case with nonhomogeneous and anisotropic plasma background. Our method yields significant speedups, both in the low and high collisional regime. In the high-collisional case, our ML-APMC outperforms the single-level APMC method by several orders of magnitude.



There are no comments yet.


page 17

page 21


A multilevel Monte Carlo method for asymptotic-preserving particle schemes

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

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...

Multilevel Adaptive Sparse Grid Quadrature for Monte Carlo models

Many problems require to approximate an expected value by some kind of M...

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

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

Multilevel Monte Carlo for Smoothing via Transport Methods

In this article we consider recursive approximations of the smoothing di...

A higher order moment preserving reduction scheme for the Stochastic Weighted Particle Method

The Stochastic Weighted Particle Method (SWPM) is a Monte Carlo techniqu...

Monte-Carlo Methods for the Neutron Transport Equation

This paper continues our treatment of the Neutron Transport Equation (NT...
This week in AI

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