[論文レビュー] 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

ひとこと要約

本稿では、材料粒子を用いてエネルギーを運搬させ、時間ステップをimplicitに処理することで、従来のimplicitモンテカルロ（IMC）スキームに内在するトランスポートエラー（テレポートエラー）を解消する2次元Multi-Frequency Implicit Semi-analog Monte-Carlo（ISMC）放射移動ソルバーを提示する。ISMC法は、1次元および2次元幾何における灰色および多周波数問題において、正確でテレポートエラーのない解を得られ、特に光学的に厚い領域や高勾配の吸収率領域で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.

研究の動機と目的

• 光学的に厚い材料において物理的に不自然な熱波の伝播を引き起こす、implicitモンテカルロ（IMC）スキームに長年存在するテレポートエラー問題に対処すること。
• 従来は明示的スキームに限られていた半アナログモンテカルロ（SMC）フレームワークをimplicit形式に拡張し、精度を損なわずに大きな時間ステップを許容すること。
• 2次元直交座標および円柱座標幾何における多数のベンチマーク問題（灰色および多周波数（多群）ケースを含む）においてISMC法の妥当性を検証すること。
• ISMCが古典的IMCよりも高速に収束することを示し、特に高吸収率勾配領域や光学的に厚い材料において顕著であることを明らかにすること。
• IMCがテレポートおよび収束問題により失敗するような、大きな時間ステップおよび粗い空間分解能下でも安定性と精度を維持できることを保証すること。

提案手法

• 材料粒子がエネルギーを運搬し、光子は既存の材料粒子からのみ発生する半アナログモンテカルロ（SMC）フレームワークを採用し、人工的な光子伝搬を排除する。
• 時間ステップの開始時の材料温度を用いて、材料粒子と光子間のエネルギー交換を解くことでimplicit時間ステッピングを実装する。
• 多周波数問題には周波数依存の吸収率モデルを用い、エネルギー依存の断面積および温度依存の放射・吸収率を設定する。
• 各時間ステップに一定数の粒子を割り当て、閉じた系において粒子数制御の必要性を排除し、発生と吸収の追跡により粒子保存則を維持する。
• 離散的な空間メッシュを用いて2次元での放射移動方程式を解き、境界条件として反射境界条件を適用し、t=0でブラックボディ源項を設定する。
• 多周波数ベンチマークではΔt = 10⁻¹³ sの時間ステップを用い、全粒子数を2.5×10⁷に制限することで数値的安定性および収束性を確保する。

実験結果

リサーチクエスチョン

• RQ1 ISMC法は、光学的に厚い材料および高吸収率勾配を有する2次元放射移動問題において、テレポートエラーを完全に解消できるか？
• RQ2 1次元および2次元の灰色および多周波数ベンチマークにおいて、ISMCは古典的IMCと比較して収束速度および精度の点で優れているか？
• RQ3 IMCがテレポートおよび最大原理違反により失敗するような、大きな時間ステップおよび粗い空間分解能下でも、ISMCは安定性と精度を維持できるか？
• RQ4 光学的に薄い領域におけるISMCの統計的ノイズはIMCと比較してどの程度か、粒子数の増加によって低減可能か？
• RQ5 不連続な材料特性を有する複雑な2次元幾何（例えばホールラウムや格子構造）において、ISMC法はスケーラビリティおよびロバスト性を有するか？

主な発見

• ISMC法は、光学的に厚いか高勾配吸収率領域を含むすべての2次元ベンチマークで、粗い空間分解能および小さな時間ステップ下でもテレポートエラーを完全に解消した。
• 2次元マーシャック波およびホールラウム問題において、ISMCは特に光学的に厚い材料において、IMCよりもはるかに速く正しい温度および放射エネルギー密度プロファイルに収束した。
• アルミブロックとフォームを含む2次元多周波数問題において、ISMCの結果はPN参照解および元のIMC実装と非常に良好に一致し、顕著なテレポート効果は観察されなかった。
• 離散的な発生・吸収に起因する光学的に薄い領域における統計的ノイズは高めであるが、ISMCは安定かつ収束的であり、粒子数の増加によってノイズを低減可能であった。
• 同じ粒子数を用いた場合、ISMCはIMCと同等の実行時間を持ち、大きな時間ステップ下でも最大原理違反に対する同程度の耐性を示した。
• すべてのベンチマークにおいて、ISMCはIMCよりも顕著に低い空間分解能要件を満たし、不透明または高勾配領域において優れた効率性を示した。