Analysis of kinetic-diffusion Monte Carlo simulation and source term estimation scheme in nuclear fusion applications

In plasma edge simulations, the behavior of neutral particles is often described by a Boltzmann–BGK equation. Solving this kinetic equation and estimating the moments of its solution are essential tasks, typically carried out using Monte Carlo (MC) methods. However, for large-scale reactors, like ITER and DEMO, high collision rates lead to a substantial computational cost. To accelerate the computational process, an asymptotic-preserving kinetic-diffusion Monte Carlo (KDMC) simulation method and a corresponding fluid estimation technique have recently been proposed. In this work, we present a numerical analysis of the convergence of KDMC combined with the associated fluid estimation in a one-dimensional setting under the diffusive scaling assumption. The analysis consists of proving theoretical upper bounds for both KDMC and the fluid estimation, and numerical verifications of these bounds. In addition, we compare the analyzed algorithm with a purely fluid-based method using the fully kinetic MC method as a reference. The algorithm consistently achieves lower error than the fluid-based method, and even one order of magnitude lower in a fusion-relevant test case. ... mehrMoreover, the algorithm exhibits a significant speedup compared to the reference kinetic MC method. Overall, our analysis confirms the effectiveness of KDMC with the associated fluid estimation in nuclear fusion applications.