热辐射输运问题的隐式蒙特卡罗方法求解

An implicit Monte Carlo method for thermal radiation transport

热辐射与物质相互作用及辐射光子在物质中的传输是惯性约束聚变研究中的重要课题.介绍了基于隐式蒙特卡罗方法的辐射输运方程,在该方程的积分-微分形式基础上,推导了利于蒙特卡罗方法模拟的等价的积分输运方程;基于积分方程设计数值模拟流程,编写三维蒙特卡罗数值模拟程序;针对热辐射输运典型问题及benchmark问题开展了数值实验,计算结果验证了方法的适应性及程序的正确性.

The equations of thermal radiative transfer are mathematical models that describe the physical processes of photons scattering through and being absorbed in and emitted by a high-energy background material. These processes are important research objects of inertial confinement fusion（ICF）. Numerical simulation is an indispensable method for this transportation equation. As an important method in the field of particle transportation, the Monte Carlo method has been widely employed in the fields of linear transportation of neutrons and high energy photons. However, traditional Monte Carlo method does not hold true when it is applied to the simulation of thermal radiation. In this research, an implicit Monte Carlo method based on calculable modeling and numerical simulation is studied. An integral transport equation that is suitable to Monte Carlo simulation is derived. A three-dimensional simulation code is developed, by which some thermal radiative transportation problems are simulated. The results of numerical experiments support that the implicit Monte Carlo method is applicable for thermal radiation transfer simulating. The work is expected to provide an important calculation method and tool for the thermal radiative transfer simulations of ICF.