反向蒙特卡罗法模拟参与性介质中热辐射传递

Computation of radiative heat transfer in participating media using Backward Monte Carlo method

为解决传统蒙特卡罗法结果收敛时间长的问题,当仅需要了解某个指定位置和方向的信息时,提出了反向蒙特卡罗法替代方法.该反向运算法则与传统的蒙特卡罗法类似,除了能束射线反向跟踪外.能束的传播是从探测位置到参与性介质中,而不是从像传统的蒙特卡罗法从源到探测位置处.利用反向蒙特卡罗法模拟计算一维和二维的热辐射传递问题,并且与相关文献进行比较,结果吻合较好.

The Monte Carlo (MC) method has been widely used to solve radiative transfer problems due to its flexibility and simplicity in simulating the energy transport process in arbitrary geometries with complex boundary conditions. However, the major drawback of the conventional (or forward) Monte Carlo method is the long computational time for converged solution. Backward Monte Carlo was considered as an alternative approach when solutions were needed at certain locations and directions. The backward algorithm was similar to the conventional method, except that the energy bundle was tracked in a direction-reversal manner. Its migration was recorded from the detection location into the participating medium, rather than from the source to the detection location as in the conventional MC. ID and 2D radiative heat transfer problems were calculated by backward MC method. The calculation results agreed well with the relative documents.