[論文レビュー] A Unified Gas-kinetic Particle Method for Multiscale Photon Transport

Weiming Li, Chang Liu+3 more|arXiv (Cornell University)|Oct 14, 2018

Gas Dynamics and Kinetic Theory参考文献 30被引用数 5

ひとこと要約

本稿では、マルチスケール光子輸送を一元的かつ滑らかに扱える統一ガス動力学的粒子（UGKP）法を提案する。時間に依存するフラックス関数を運動方程式から導出することで、粒子輸送と衝突を結合し、漸近的保存性を達成。光線的厚さの極限では拡散方程式を正確に回復し、自由輸送状態では正確な粒子追跡が可能であり、ベンチマーク問題を用いた全スケールでの検証が行われた。

ABSTRACT

In this work, we present a unified gas-kinetic particle (UGKP) method for the simulation of multiscale photon transport. The multiscale nature of the particle method mainly comes from the recovery of the time evolution flux function in the unified gas-kinetic scheme (UGKS) through a coupled dynamic process of particle transport and collision. This practice improves the original operator splitting approach in the Monte Carlo method, such as the separated treatment of particle transport and collision. As a result, with the variation of the ratio between numerical time step and local photon's collision time, different transport physics can be fully captured in a single computation. In the diffusive limit, the UGKP method could recover the solution of the diffusion equation with the cell size and time step being much larger than the photon's mean free path and the mean collision time. In the free transport limit, it presents an exact particle tracking process as the original Monte Carlo method. In the transition regime, the weights of particle free transport and collision are determined by the ratio of local numerical time step to the photon's collision time. Several one-dimensional numerical examples covering all transport regimes from the optically thin to optically thick are computed to validate the accuracy and efficiency of the current scheme. In comparison with the $S_N$ discrete ordinate method, the UGKP method is based on particles and avoids the discretization of particle velocity space, which does not suffer from the ray effect.

研究の動機と目的

• 自由輸送から拡散まで、あらゆるスケールの光子輸送を、領域分割や別々のソルバーを必要とせずに正確にシミュレートできる統一数値法の開発を目的とする。
• 従来のモンテカルロ法が拡散状態で過剰な衝突追跡を要し、計算コストが著しく増加するという限界を克服することを目的とする。
• モンテカルロ法の利点（光線効果なし）を保ちつつ、漸近的保存性を達成する粒子ベースの定式化として、統一ガス動力学的スキーム（UGKS）を拡張することを目的とする。
• 放射体と物質の結合系に対する検証を通じて、全スケールの輸送状態において、エネルギー交換と平衡状態が正確に再現されることを確認することを目的とする。

提案手法

• UGKP法は、放射輸送方程式の方向空間を粒子ベースで離散化し、決定的法で必要な速度空間の離散化を回避する。
• 時間積分により運動方程式を統合することで、マルチスケールのフラックス関数を構築し、粒子輸送と衝突プロセスを動的に結合する。
• フラックスの時間発展は、数値時間ステップと局所的な光子平均自由行程の比に依存し、スケール間の自動的遷移を可能にする。
• 粒子の重みは、自由輸送と衝突の相対的寄与を表しており、局所的な衝突時間スケールに従って決定される。
• 演算子分割を置き換える統一された粒子時間発展スキームを採用し、全スケールで輸送と衝突を一貫して取り扱う。
• 放射体と物質の結合系への拡張には、放射と物質間のエネルギー交換項を組み込み、物質エネルギーは熱的緩和モデルに従って進化させる。

実験結果

リサーチクエスチョン

• RQ1 粒子ベースの手法が、時間ステップの制限なしに、光線的厚さの極限で拡散方程式を正確に回復できるか、すなわち漸近的保存性を達成できるか？
• RQ2 全スケールの輸送状態で精度を保つために、粒子輸送と衝突を1つのアルゴリズムでどのように結合できるか？
• RQ3 自由輸送状態ではモンテカルロ法と同等の正確な粒子追跡を再現し、同時に拡散状態でも安定性を保てるか？
• RQ4 放射体と物質の結合系を、全スケール、特に拡散極限での平衡状態を含めて正確にシミュレートできるか？
• RQ5 多スケール輸送問題において、UGKP法は離散座標法（DOM）やインクリメンタルモンテカルロ法と比べて、精度と効率の面で優れているか？

主な発見

• 光線的厚さの極限（ε = 10⁻⁴）において、UGKP法は、光子平均自由行程や衝突時間よりもはるかに大きな時間ステップとセルサイズでも、正確に拡散方程式を回復した。
• 自由輸送状態（ε = 1）では、UGKP解は離散座標法で得た基準解と高い精度で一致し、正確な粒子追跡が実現していることが確認された。
• 中間スケール（ε = 10⁻²）でも、UGKP法は基準解と優れた一致を示し、遷移スケールでの頑健性を示した。
• 拡散状態（ε = 10⁻⁴）では、UGKP解は平衡拡散解に収束し、漸近的保存性が妥当であることが検証された。
• 放射体と物質の結合系では、拡散状態において放射と物質のエネルギーが平衡に達し、t = 0.03におけるEとcu_rがよく一致した。
• 決定的法に内在する光線効果を回避し、スケール間の領域分割やハイブリッドソルバーの必要性も排除した。