[論文レビュー] A fast algorithm for radiative transport in isotropic media

Kui Ren, Rongting Zhang+1 more|arXiv (Cornell University)|Oct 4, 2016

Numerical methods in inverse problems参考文献 63被引用数 4

ひとこと要約

本稿では、等方的媒体における定常状態放射移動方程式（ERT）を解く高速なアルゴリズムを提示する。この手法は、方向平均解についてのフリードホルム積分方程式に再定式化し、その後、核に依存しない高速多重極法（FMM）とGMRESを組み合わせて効率的に解く。この方法により、計算量が線形に増加し、均質媒体および不均質媒体の両方を変更なしに処理でき、主な計算フェーズで方向離散化を回避するため、従来の手法に比べて顕著な高速化を実現する。

ABSTRACT

We propose in this work a fast numerical algorithm for solving the equation of radiative transfer (ERT) in isotropic media. The algorithm has two steps. In the first step, we derive an integral equation for the angularly averaged ERT solution by taking advantage of the isotropy of the scattering kernel, and solve the integral equation with a fast multipole method (FMM). In the second step, we solve a scattering-free transport equation to recover the original ERT solution. Numerical simulations are presented to demonstrate the performance of the algorithm for both homogeneous and inhomogeneous media.

研究の動機と目的

• 高次元性のため計算的に困難な、等方的媒体における放射移動方程式（ERT）を効率的に解く数値的手法の開発。
• 方向離散化に依存する従来の決定的スケーラーの高い計算コストとメモリ要件を克服すること。
• 統一されたフレームワークを用いて、均質および不均質媒体の両方でERTを高速に解けるようにすること。
• 積分方程式の定式化と高速多重極法を活用し、正確性を保ちながら解法プロセスを加速すること。

提案手法

• ERTを方向平均解 ∫Φ(x,v)d𝑣 についての第二種フリードホルム積分方程式に再定式化する。
• 得られた体積積分方程式を解くために、GMRESなどの反復解法を適用する。
• Chebyshev補間を用いた核に依存しない高速多重極法（FMM）により、GMRESにおける行列・ベクトル積を高速化する。
• 低ランク近似技術を用いてFMM内の相互作用項を圧縮し、計算コストを削減する。
• 方向平均解を取得した後、散乱なしの輸送方程式を解くことで、解の全方向的依存性を回復する。
• 不均質媒体における繰り返し計算を削減するため、全吸収係数の線積分評価をキャッシュ化する。

実験結果

リサーチクエスチョン

• RQ1 等方的媒体におけるERTは、主な計算フェーズで方向離散化を回避する体積積分方程式に再定式化可能か？
• RQ2 ERTの積分定式化に高速多重極法を効果的に適用することで、問題サイズに比例して線形スケーリングが達成可能か？
• RQ3 構造的変更なしに、均質および不均質媒体の両方で正確性と効率性を維持できるか？
• RQ4 Chebyshev補間ノードを少数に制限した場合、FMM加速解法の精度と収束性にどのような影響があるか？
• RQ5 この手法は、複雑または異方的媒体におけるより一般的な輸送スケーラーの前処理に用いることができるか？

主な発見

• 提案手法は、空間的問題サイズに対して線形の計算複雑性を達成し、大規模問題においてもスケーラビリティを示した。
• 少数のChebyshev補間ノードを用いたFMMは、真の解に対する高精度な近似をもたらし、前処理としての利用可能性を示唆している。
• 本手法は、均質および不均質媒体の両方で頑健であり、媒体タイプの切り替えに対してもアルゴリズム構造の変更が不要である。
• 均質および不均質媒体の両方において、積分方程式のカーネルは構造的に同一のままであり、実装を単純化している。
• 数値シミュレーションにより、本手法の効率性と正確性が確認され、複雑な散乱および吸収プロファイルを含む多様なシナリオで妥当性が検証された。
• 計算量が最も高いフェーズで方向離散化を回避することで、メモリ使用量が顕著に削減され、より大きな問題の解法が可能になった。