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

by   Bert Mortier, et al.

We develop a novel Monte Carlo strategy for the simulation of the Boltzmann-BGK model with both low-collisional and high-collisional regimes present. The presented solution to maintain accuracy in low-collisional regimes and remove exploding simulation costs in high-collisional regimes uses hybridized particles that exhibit both kinetic behaviour and diffusive behaviour depending on the local collisionality. In this work, we develop such a method that maintains the correct mean, variance, and correlation of the positional increments over multiple time steps of fixed step size for all values of the collisionality, under the condition of spatial homogeneity during the time step. In the low-collisional regime, the method reverts to the standard velocity-jump process. In the high-collisional regime, the method collapses to a standard random walk process. We analyze the error of the presented scheme in the low-collisional regime for which we obtain the order of convergence in the time step size. We furthermore provide an analysis in the high-collisional regime that demonstrates the asymptotic-preserving property.



There are no comments yet.


page 1

page 2

page 3

page 4


Multilevel Asymptotic-Preserving Monte Carlo for Particle Simulations

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

A multilevel Monte Carlo method for asymptotic-preserving particle schemes

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

Balanced data assimilation with a blended numerical model

A challenge arising from the local Bayesian assimilation of data in an a...

Uniformly accurate integrators for Klein-Gordon-Schrödinger systems from the classical to non-relativistic limit regime

In this paper we present a novel class of asymptotic consistent exponent...

Geometrically Convergent Simulation of the Extrema of Lévy Processes

We develop a novel Monte Carlo algorithm for the simulation from the joi...

A New Simulation Algorithm for Absorbing Receiver in Molecular Communication

The simulation of diffusion-based molecular communication systems with a...

Simulations of Dynamical Gas-Dust Circumstellar Disks: Going Beyond the Epstein Regime

In circumstellar disks, the size of dust particles varies from submicron...
This week in AI

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