[論文レビュー] Radiative Transfer For Variable 3D Atmospheres

F. Golse, F. Hecht+3 more|arXiv (Cornell University)|Aug 11, 2022

Thermal Radiation and Cooling Technologies被引用数 1

ひとこと要約

本稿では、H行列圧縮と有限要素離散化を用いて、空間的・スペクトル的に変化する吸収・散乱を有する大気における3次元放射移動を、O(N log N)の高速な数値的手法として提示する。この手法により、複雑な地形における温室効果ガス（CO2、CH4）の吸収に起因する温度差の高精度なシミュレーションが可能となり、50,000メッシュポイントと680周波数帯域を用いた計算をM1ラップトップ上で5分未塔で実行可能である。この手法は積分方程式形式とナビエ＝ストークス方程式およびエネルギー収支方程式を組み合わせることで実現されている。

ABSTRACT

To study the temperature in a gas subjected to electromagnetic radiations, one may use the Radiative Transfer equations coupled with the Navier-Stokes equations. The problem has 7 dimensions; however with minimal simplifications it is equivalent to a small number of integro-differential equations in 3 dimensions. We present the method and a numerical implementation using an H-matrix compression scheme. The result is a very fast: 50K physical points, all directions of radiation and 680 frequencies require less than 5 minutes on an Apple M1 Laptop. The method is capable of handling variable absorptioN and scattering functionS of spatial positions and frequencies. The implementation is done using htool, a matrix compression library interfaced with the PDE solver freefem++. Applications to the temperature in the French Chamonix valley is presented at different hours of the day with and without snow / clouds and with a variable absorption taken from the Gemini measurements. The result is precise enough to assert temperature differences due to increased absorption in the vibrational frequency subrange of greenhouse gasses.

研究の動機と目的

• 空間的・スペクトル的に変化する吸収・散乱を有する大気における3次元放射移動を解く数値的に効率的な手法の開発。
• SHDOM やモンテカルロ法などの従来手法に比べて高い計算コストがかかる問題、特に細かい周波数分解能および方向分解能を必要とする場合の課題を克服すること。
• CO2 や CH4 の振動帯域における吸収係数の微小な変化に起因する温度変化の正確なシミュレーションを可能とすること。
• 実際の地形（フランス・シャモニーアルプス）を用いた高解像度メッシュと現実的な大気データを用いて、数値ツールの妥当性を検証すること。
• 標準的なハードウェア上でも気候関連のシミュレーションに実用的かつスケーラブルであることを示すこと。

提案手法

• 放射移動方程式を立体角方向に平均化した積分方程式形式を用いることで、7次元の放射移動問題を3次元系に低次元化する。
• 積分方程式に対して固定点反復スキームを適用し、数学的に単調かつ収束することが保証されており、安定性と精度が確保される。
• H行列圧縮を用いて、密行列演算の計算コストをO(N²)からO(N log N)に削減し、高速な行列・ベクトル積を実現する。
• 複雑な3次元地形を扱うために、自由にメッシュを生成可能な非構造化三角形メッシュを用いた有限要素法をFreeFem++で実装し、実地形データから高精度なメッシュ生成を実現する。
• ナビエ＝ストークス方程式は定常状態で、平均化されたレイノルズ形式を用い、乱流は有効粘性係数でモデル化する。この方程式は、放射加熱を含むエネルギー方程式と結合される。
• 放射源項は圧縮されたH行列を用いて計算され、境界・内部行列要素の混合成分に対しても特別な処理が施されており、積分方程式の効率的解法が可能となる。

実験結果

リサーチクエスチョン

• RQ1 空間的・スペクトル的に変化する吸収・散乱を有する複雑な3次元大気における、高速かつスケーラブルな数値的手法を開発可能か？
• RQ2 H行列圧縮は、精度を損なわずに放射移動シミュレーションの計算コストをどの程度低減できるか？
• RQ3 CO2 や CH4 の振動周波数帯域における吸収に起因する微小な温度差を、本手法が解像可能か？
• RQ4 空間的吸収係数が変化する現実的な3次元大気条件下でも、積分方程式に対する固定点反復スキームは安定的かつ収束的か？
• RQ5 非構造化メッシュを用いて、垂直勾配が急で複雑な地形（例：フランス・アルプス）にも適用可能か？

主な発見

• Apple M1ラップトップ上での計算で、50,000物理点、全方向の放射、680周波数帯域をカバーする完全な3次元放射移動シミュレーションが5分未塔で実行可能である。
• H行列圧縮により、計算コストがO(N²)からO(N log N)に削減され、大規模問題の効率的取り扱いが可能となった。
• 数値解法は、Gemini衛星観測データを用いてモデル化されたCO2およびCH4の振動帯域における吸収増加に起因する温度差を解像できるほど高精度である。
• 本手法はシャモニーグループの実際の温度プロファイルを正確に再現できており、雪の有無、雲の有無、吸収係数の変動に起因する明確な熱的応答の違いを示している。
• 積分方程式に対する反復スキームは、単調かつ収束することが数学的に証明されており、上界解と下界解から数値誤差の上限を評価可能である。
• コードはGitHubで公開されており、より大きな領域やスーパーコンピュータへの拡張が可能であり、数百万メッシュノードを数分で処理できる。