Multilevel Monte Carlo simulation of Coulomb collisions
From MaRDI portal
Publication:728652
Abstract: We present a new, for plasma physics, highly efficient multilevel Monte Carlo numerical method for simulating Coulomb collisions. The method separates and optimally minimizes the finite-timestep and finite-sampling errors inherent in the Langevin representation of the Landau-Fokker-Planck equation. It does so by combining multiple solutions to the underlying equations with varying numbers of timesteps. For a desired level of accuracy epsilon, the computational cost of the method is order(epsilon^{-2}) or order(epsilon^{-2} (ln epsilon)^2), depending on the underlying discretization, Milstein or Euler-Maruyama respectively. This is to be contrasted with a cost of order(epsilon^{-3}) for direct simulation Monte Carlo or binary collision methods. We successfully demonstrate the method with a classic beam diffusion test case in 2D, making use of the Levy area approximation for the correlated Milstein cross terms, and generating a computational saving of a factor of 100 for epsilon = 10^{-5}. We discuss the importance of the method for problems in which collisions constitute the computational rate limiting step, and its limitations.
Recommendations
- Monte Carlo methods and their analysis for Coulomb collisions in multicomponent plasmas
- Higher-order time integration of Coulomb collisions in a plasma using Langevin equations
- A Monte Carlo method with negative particles for Coulomb collisions
- A Hybrid Method for Accelerated Simulation of Coulomb Collisions in a Plasma
- On the accuracy of direct simulation the Landau collision integral by the Boltzmann integral
Cites work
- scientific article; zbMATH DE number 1215244 (Why is no real title available?)
- scientific article; zbMATH DE number 2000348 (Why is no real title available?)
- scientific article; zbMATH DE number 3806623 (Why is no real title available?)
- scientific article; zbMATH DE number 933353 (Why is no real title available?)
- scientific article; zbMATH DE number 3024704 (Why is no real title available?)
- scientific article; zbMATH DE number 3068091 (Why is no real title available?)
- A Monte Carlo method for Coulomb collisions in hybrid plasma models
- A binary collision model for plasma simulation with a particle code
- A grid-based Coulomb collision model for PIC codes
- A review of Vlasov-Fokker-Planck numerical modeling of inertial confinement fusion plasma
- Adaptive multilevel Monte Carlo simulation
- An entropy based thermalization scheme for hybrid simulations of Coulomb collisions
- Analyzing multi-level Monte Carlo for options with non-globally Lipschitz payoff
- Antithetic Multilevel Monte Carlo Estimation for Multidimensional SDEs
- Antithetic multilevel Monte Carlo estimation for multi-dimensional SDEs without Lévy area simulation
- Direct simulation Monte Carlo schemes for Coulomb interactions in plasmas
- Fokker-Planck Equation for an Inverse-Square Force
- Higher-order implicit strong numerical schemes for stochastic differential equations
- Higher-order time integration of Coulomb collisions in a plasma using Langevin equations
- How to generate random matrices from the classical compact groups
- ISDEP: integrator of stochastic differential equations for plasmas
- Improved multilevel Monte Carlo convergence using the Milstein scheme
- Infinite-dimensional quadrature and approximation of distributions
- Langevin representation of Coulomb collisions in PIC simulations
- Monte Carlo complexity of global solution of integral equations
- Monte Carlo methods and their analysis for Coulomb collisions in multicomponent plasmas
- Multilevel Monte Carlo Path Simulation
- Multilevel Monte Carlo approximation of distribution functions and densities
- Multilevel quasi-Monte Carlo path simulation
- On irregular functionals of SDEs and the Euler scheme
- Particle simulation of Coulomb collisions: Comparing the methods of Takizuka \& Abe and Nanbu
- Random Generation of Stochastic Area Integrals
- Randomized multilevel quasi-Monte Carlo path simulation
- Small-angle Coulomb collision model for particle-in-cell simulations
- The Boltzmann equation and its applications
- The semi-Lagrangian method for the numerical resolution of the Vlasov equation
- Variable Step Size Control in the Numerical Solution of Stochastic Differential Equations
Cited in
(26)- Effect of collision-partner selection schemes on the accuracy and efficiency of the direct simulation Monte Carlo method
- Weakly convergent stochastic simulation of electron collisions in plasmas
- Multilevel importance sampling for rare events associated with the McKean-Vlasov equation
- Direct simulation Monte Carlo schemes for Coulomb interactions in plasmas
- An explicitly solvable energy-conserving algorithm for pitch-angle scattering in magnetized plasmas
- Higher-order time integration of Coulomb collisions in a plasma using Langevin equations
- Meta variance reduction for Monte Carlo estimation of energetic particle confinement during stellarator optimization
- Small-angle Coulomb collision model for particle-in-cell simulations
- Accelerated Simulation of Boltzmann-BGK Equations near the Diffusive Limit with Asymptotic-Preserving Multilevel Monte Carlo
- Monte Carlo methods and their analysis for Coulomb collisions in multicomponent plasmas
- Particle simulation of Coulomb collisions: Comparing the methods of Takizuka \& Abe and Nanbu
- Collisional \(N\)-body numerical integrator with applications to charged particle dynamics
- An entropy based thermalization scheme for hybrid simulations of Coulomb collisions
- Multiscale modelling and splitting approaches for fluids composed of Coulomb-interacting particles
- Adaptive time-stepping Monte Carlo integration of Coulomb collisions
- Monte Carlo simulation of excitation and ionization collisions with complexity reduction
- Variance reduced particle solution of the Fokker-Planck equation with application to rarefied gas and plasma dynamics
- Understanding the accuracy of Nanbu's numerical Coulomb collision operator
- Composition schemes for the stochastic differential equation describing collisional pitch-angle diffusion
- Simulating Coulomb and log-gases with hybrid Monte Carlo algorithms
- Multilevel and multi-index Monte Carlo methods for the McKean-Vlasov equation
- MLMC for nested expectations
- Data-driven stochastic particle scheme for collisional plasma simulations
- Accurate particle time integration for solving Vlasov-Fokker-Planck equations with specified electromagnetic fields
- A Monte Carlo method with negative particles for Coulomb collisions
- On the accuracy of direct simulation the Landau collision integral by the Boltzmann integral
This page was built for publication: Multilevel Monte Carlo simulation of Coulomb collisions
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q728652)
