[论文解读] Multi-Frequency Implicit Semi-analog Monte-Carlo (ISMC) Radiative Transfer Solver in Two-Dimensions (without Teleportation)

Elad Steinberg, Shay I. Heizler|arXiv (Cornell University)|Aug 5, 2021

Radiative Heat Transfer Studies参考文献 31被引用 16

一句话总结

本论文提出了一种二维多频隐式半随机蒙特卡罗（ISMC）辐射转移求解器，通过使用物质粒子携带能量并实现隐式时间推进，消除了传统隐式蒙特卡罗（IMC）方法固有的跃迁误差。ISMC方法在1D和2D几何中对灰体和多频问题均实现了准确、无跃迁误差的解，且在光学厚和高梯度消光系数区域收敛速度快于IMC。

ABSTRACT

We study the multi-dimensional radiative transfer phenomena using the ISMC scheme, in both gray and multi-frequency problems. Implicit Monte-Carlo (IMC) schemes have been in use for five decades. The basic algorithm yields teleportation errors, where photons propagate faster than the correct heat front velocity. Recently [Po\"ette and Valentin, J. Comp. Phys., 412, 109405 (2020)], a new implicit scheme based on the semi-analog scheme was presented and tested in several one-dimensional gray problems. In this scheme, the material energy of the cell is carried by material-particles, and the photons are produced only from existing material particles. As a result, the teleportation errors vanish, due to the infinite discrete spatial accuracy of the scheme. We examine the validity of the new scheme in two-dimensional problems, both in Cartesian and Cylindrical geometries. Additionally, we introduce an expansion of the new scheme for multi-frequency problems. We show that the ISMC scheme presents excellent results without teleportation errors in a large number of benchmarks, especially against the slow classic IMC convergence.

研究动机与目标

• 解决隐式蒙特卡罗（IMC）方法中长期存在的跃迁误差问题，该问题在光学厚材料中导致非物理解的热波传播。
• 将此前仅限于显式格式的半随机蒙特卡罗（SMC）框架扩展为隐式形式，以在不损失精度的前提下实现大时间步长。
• 在二维笛卡尔和柱对称几何中对一系列基准问题（包括灰体和多频（多群）情况）验证ISMC方法。
• 证明ISMC在高消光系数梯度或光学厚材料区域的收敛速度优于经典IMC。
• 确保在大时间步长和低空间分辨率条件下保持稳定性和精度，而IMC因跃迁误差和收敛性问题在此类条件下失效。

提出的方法

• 采用半随机蒙特卡罗（SMC）框架，其中物质粒子携带能量，光子仅从现有物质粒子发射，从而消除人工光子传播。
• 通过求解时间步长内物质粒子与光子之间的能量交换，实现隐式时间推进，使用时间步长开始时的物质温度。
• 在多频问题中采用频率相关的消光系数模型，包含能量依赖的截面以及温度相关的发射和吸收速率。
• 每个时间步采用固定数量的粒子，避免在封闭系统中使用种群控制，并通过发射和吸收过程的追踪保持粒子守恒。
• 在二维空间上使用离散网格求解辐射转移方程，施加反射边界条件，并在t=0时应用黑体源项。
• 在多频基准测试中使用Δt = 10⁻¹³ s的时间步长，并将总粒子数限制在2.5×10⁷以内，以确保数值稳定性和收敛性。

实验结果

研究问题

• RQ1 ISMC方法是否能在光学厚材料和高消光系数梯度的二维辐射转移问题中消除跃迁误差？
• RQ2 在1D和2D灰体与多频基准测试中，ISMC与经典IMC相比在收敛速度和精度方面表现如何？
• RQ3 在大时间步长和粗空间分辨率条件下，ISMC是否能保持稳定性和精度，而IMC因跃迁误差和最大值原理违反而失效？
• RQ4 ISMC在光学薄区域的统计噪声与IMC相比如何，是否可通过增加粒子数来缓解？
• RQ5 在具有不连续材料属性的复杂2D几何（如空腔和晶格结构）中，ISMC方法是否具备可扩展性和鲁棒性？

主要发现

• 在所有2D基准测试中，包括光学厚和高梯度消光系数区域，ISMC方法成功消除了跃迁误差，即使在粗空间分辨率和小时间步长下也表现良好。
• 在2D Marshak波和空腔问题中，ISMC比IMC更快收敛到正确的温度和辐射能量密度分布，尤其在光学厚材料中表现显著。
• 在含铝块和泡沫的2D多频问题中，ISMC结果与PN参考解以及原始IMC实现高度一致，未观察到任何跃迁效应。
• 尽管由于离散发射和吸收导致光学薄区域统计噪声较高，ISMC仍保持稳定和收敛，且可通过增加粒子数降低噪声。
• 在使用相同粒子数时，ISMC的运行时间与IMC相当，并在大时间步长下对最大值原理违反具有相似的容忍度。
• 在所有基准测试中，ISMC在不透明或高梯度区域的收敛解所需的空间分辨率显著低于IMC，证明了其在这些区域的优越效率。