2025/12/10 by G.J. Parker, M. Umansky, Parker, Gregory J. +3
Mathematics · Physics and Astronomy · #Computational Physics (physics.comp-ph) #Dust and Plasma Wave Phenomena #FOS: Mathematics #FOS: Physical sciences #Gas Dynamics and Kinetic Theory #Magnetic confinement fusion research #Numerical Analysis (math.NA) #Plasma Physics (physics.plasm-ph)
paper · pdf · doi:10.48550/arxiv.2512.09334
openalex publication_date 2025/12/10 · openalex created_date 2025/12/12 · openalex updated_date 2026/07/28
We propose and investigate a new multi-level Monte Carlo scheme for numerical solutions of the kinetic Boltzmann equation for neutral species in edge plasmas. In particular, this method explicitly exploits a key structural property of neutral particle dynamics: the prevalence of frequent collisions for which the outgoing velocity is determined by local plasma parameters. Using this property, we derive a multi-level algorithm based on collision event propagator and show, both analytically and through numerical experiments, that it reproduces the results of standard Monte Carlo methods. We further demonstrate that, in the context of coupled plasma-neutral edge simulations employing correlated Monte Carlo, the proposed scheme retains trajectory correlation to machine precision as the system evolves, whereas conventional methods exhibit rapid decorrelation. These results indicate that the propagator-based multi-level Monte Carlo scheme is a promising candidate for use in fully implicit Jacobian-free Newton-Krylov (JFNK) solvers for coupled plasma-neutral systems.