用于蒙特卡洛算法中稀有事件采样的保矩与网格自适应重加权方法
计算物理
2021-07-07 v2 数值分析
数值分析
等离子体物理
摘要
我们提出了一种用于稀有事件采样的新型轮盘赌方案，该方案既保持结构又无偏。蒙特卡洛标记分裂和删除的边界在运行期间自动放置并自适应。扩展现有代码以使用新方案无需进行重大更改，因为标记的运动方程没有被改变。这些方案也可以在标记之间存在非线性和非局部耦合的情况下应用。作为一个说明性应用，我们在 ASCOT-RFOF 代码中实现了该方法，该代码用于模拟聚变等离子体中射频波产生的快离子。在此应用中，使用该方法可以在不增加计算量的情况下，将典型快离子能量的蒙特卡洛噪声水平至少降低一个数量级。
引用
@article{arxiv.2003.14048,
title = {Moment-Preserving and Mesh-Adaptive Reweighting Method for Rare-Event Sampling in Monte-Carlo Algorithms},
author = {C. U. Schuster and T. Johnson and G. Papp and R. Bilato and S. Sipilä and J. Varje and M. Hasenöhrl},
journal= {arXiv preprint arXiv:2003.14048},
year = {2021}
}