库仑碰撞的多层蒙特卡洛模拟

我们提出了一种针对等离子体物理的全新、高效的多层蒙特卡洛数值方法，用于模拟库仑碰撞。该方法分离并最优地最小化了 Landau-Fokker-Planck 方程 Langevin 表述中固有的有限时间步长和有限采样误差。其实现方式是结合底层方程具有不同时间步长的多个解。对于所需的精度水平 epsilon，该方法的计算成本为 order(epsilon^{-2}) 或 order(epsilon^{-2} (\ln epsilon)^2)，具体取决于底层离散化方案分别为 Milstein 或 Euler-Maruyama。相比之下，直接模拟蒙特卡洛或二元碰撞方法的成本为 order(epsilon^{-3})。我们成功地在 2D 经典束流扩散测试案例中展示了该方法，利用 Levy 面积近似处理相关的 Milstein 交叉项，并在 epsilon = 10^{-5} 时实现了 100 倍的计算节省。我们讨论了该方法在碰撞构成计算速率限制步骤的问题中的重要性及其局限性。