摘要：用蒙特卡罗（MC）方法模拟高温、高压、多介质、大变形辐射输运问题时,由于网格体积悬殊,导致各网格通量的统计误差涨落很大,随着时间步的增加,误差积累甚至会导致计算结果失真.为此,发展了针对全局网格计算的源偏倚抽样技巧.用于源偏倚抽样的价值函数基于上个时间步各网格通量及误差,通过加权构造产生,它比传统MC通过解伴随方程获取价值的性价比要高得多.数值试验表明,全局源偏倚抽样通过自适应分配当前时间步各网格的粒子数,有效地降低了当前步重要网格通量误差.

关键词：非定常 辐射输运 蒙特卡罗 源偏倚抽样 

单位：北京应用物理与计算数学研究所计算物理实验室 北京100088 中国工程物理研究院北京研究生部 北京100088