[論文レビュー] Higher order minimum entropy approximations in radiative transfer

Philipp Monreal, Martin Frank|ArXiv.org|Dec 16, 2008

Gas Dynamics and Kinetic Theory参考文献 13被引用数 14

ひとこと要約

本稿は、放射移動における高次最小エントロピー閉じ込めのための数値フレームワークを構築し、$χ_n$ 法を一次（$ olimits\mathcal{M}_1$）から$ olimits\mathcal{M}_2$ およびそれ以上の次数へと拡張する。モーメント系のエントロピー最小化問題を数値的に解くことで、$ olimits\mathcal{M}_1$ モデルがゼロネットフラックス状況で失敗する問題を克服し、2本のビームまたは内部源を伴う1次元シミュレーションにおいて、より高い精度と物理的整合性を示す。

ABSTRACT

In this paper we approximate the radiative transfer equations by the method of moments, constructing mesoscopic approximations of arbitrary order of the otherwise microscopic system. To define the necessary closure a minimum entropy approach is utilized. While in radiative transfer, the minimum entropy closure for moment systems up to the first-order moment is well known, higher-order minimum entropy closures have not been implemented. This is probably due to the fact that the closure cannot be expressed in analytical form. Our focus thus lies in developing some general results about the minimum entropy system and in deriving a numerical closure. By extending to higher order, among increasing the precision, we are able to overcome difficulties that arise for the first order minimum entropy method. Numerical experiments in a 1-dimensional domain irradiated by two beams or with internal source show the accuracy of this approach.

研究の動機と目的

• 高次最小エントロピー閉じ込めを放射移動に実装するという長年の課題に取り組むこと。これは解析的形が存在せず、ほとんど未開拓であった。
• 一次を超えるモーメント系の閉じ込めに不可欠なエディングトン因子を計算するための数値アルゴリズムを開発すること。
• $\mathcal{M}_2$ やそれ以上の次数のモデルが、特にゼロネットフラックスにおける放射非平衡状態の取り扱いにおいて、$\mathcal{M}_1$ の根本的限界を克服できることを示すこと。
• 2本のビームや1次元領域内の有限時間内部源を伴うベンチマークテストケースを用いて、数値的シミュレーションにより、$\mathcal{M}_1$ や球面調和法と比較して物理的整合性と精度が向上することを検証すること。
• 将来的な多次元問題への拡張および放射線治療やプラズマ物理学における幅広い応用の基盤を築くこと。

提案手法

• 角度変数についての積分を用いて、放射移動方程式からモーメント法を適用し、マクロな偏微分方程式の階層を導出する。
• 最小エントロピー原理を適用してモーメント系を閉じ、制約下でのエントロピー最大化により物理的整合性を保証する。
• 離散グリッド上でラグランジュ乗数を用いて双対最適化問題を解くことで、閉じ込めの鍵たるエディングトン因子を数値的に計算するアルゴリズムを開発する。
• 空間的および角度的離散化後に得られる常微分方程式系を、適応的ルンゲ・クッタ法で解く。
• モーメントの対称性を活用し、計算コストを削減する。半分のモーメント空間のみを計算し、結果を補間する。
• 2本の反対方向ビームと有限時間内部源を伴う1次元領域におけるベンチマークテストケースを用いて、手法の妥当性を検証する。

実験結果

リサーチクエスチョン

• RQ1 解析的表現が存在しないにもかかわらず、$\mathcal{M}_n$（$n \geq 2$）の高次最小エントロピー閉じ込めを数値的に実装可能か？
• RQ2 $\mathcal{M}_1$ が失敗する状況、特にネットフラックスがゼロのとき、$\mathcal{M}_2$ モデルはどのように性能を発揮するか？
• RQ3 $\mathcal{M}_2$ モデルは、$\mathcal{M}_1$ や球面調和近似（$P_1$, $P_3$）と比較して、どの程度精度と物理的整合性を向上させるか？
• RQ4 数値的に計算されたエディングトン因子は、モーメント系の安定化、双曲型性、エントロピー散逸の実現にどのような役割を果たすか？
• RQ5 $\mathcal{M}_2$ モデルは、有限時間持続の内部源のような、一時的で不連続な源項の時間発展を正確に捉えることができるか？

主な発見

• $\mathcal{M}_1$ モデルは、ゼロネットフラックスを伴う2本のビームテストケースで物理的に非現実な衝撃波を生成し、放射非平衡状態の取り扱いに根本的な欠陥を示している。
• $\mathcal{M}_2$ モデルは、物理的に整合的で正のエネルギー分布を維持しながら、非現実な衝撃波を完全に排除した。
• 有限時間持続の内部源を伴う無限領域において、$\mathcal{M}_2$ は$\mathcal{M}_1$ よりも顕著に高い精度を示し、ベンチマーク解からのずれが小さくなった。
• $\mathcal{M}_2$ モデルは光速よりも遅い情報伝播を示し、因果律を満たしている。これは拡散モデルとは対照的である。
• 数値的に計算されたエディングトン因子は安定しており、時間依存シミュレーションに適した双曲型的かつエントロピー散逸型の系を構築可能にした。
• この手法は解の正値性を維持し、特に不連続性付近で顕著に現れる球面調和近似の振動を回避した。