[論文レビュー] Preconditioning a coupled model for reactive transport in porous media

Laila Amir, Michel Kern|arXiv (Cornell University)|Oct 4, 2017

Numerical methods for differential equations参考文献 66被引用数 4

ひとこと要約

本稿では、移動種の濃度を排除し、全固定濃度のみを解くことで、多孔質媒体内の結合反応輸送問題を解くためのプリコンディショニング付きニュートン・クリロフ法を提示する。この手法は、ブロックガウス＝ザイデルまたはブロックヤコビプリコンディショニングを用い、非線形および線形反復においてメッシュに依存しない収束性を達成し、標準的手法に比べて大幅に計算効率が向上する。

ABSTRACT

We study numerical methods for solving reactive transport problems in porous media that allow a separation of transport and chemistry at the software level, while keeping a tight numerical coupling between both subsystems. After recalling how each individual subsystem can be solved, we give a formulation that eliminates the local chemical concentrations, and keeps the total concentrations as unknowns. The block structure of the model is exploited both at the nonlinear level, by eliminating some unknowns, and at the linear level by using block Gauss-Seidel or block Jacobi preconditioning. The methods are applied to the easy 1D case of the MoMaS benchmark.

研究の動機と目的

• グローバル陰解法の反応輸送シミュレーションにおける計算コストを低減し、収束を加速すること。
• 移動種の濃度を排除し、全固定濃度のみを未知数として解く定式化を構築し、システムの簡略化を図ること。
• メッシュに依存しない収束性を保持するブロックプリコンディショニング手法の分析と適用。
• ベンチマーク1次元反応輸送問題への応用を通じて、本手法の頑健性とスケーラビリティを示すこと。
• 最小限のソフトウェア結合で、高い精度を維持したまま、複雑な反応輸送の効率的シミュレーションを可能にすること。

提案手法

• 本手法は、輸送および化学反応の完全に結合された非線形方程式系を解くためにニュートン・クリロフ法を用いる。
• 局所的な化学的濃度を排除し、全濃度のみを未知数として保持する形に問題を再定式化する。
• ニュートン法から生じる線形系に、ブロックガウス＝ザイデルおよびブロックヤコビプリコンディショニングを適用する。
• このプリコンディショニング戦略は、非線形レベルでの変数の消去と同等であり、数値的安定性を向上させる。
• ヤコビ行列はベクトル積を介して正確に計算され、輸送および化学ソルバーのヤコビ行列へのアクセスを必要とする。
• 質量作用則の直接代入を避けることで、モジュール性を保ちつつ、精度を向上させる。

実験結果

リサーチクエスチョン

• RQ1 ブロックプリコンディショニング戦略は、反応輸送シミュレーションにおいてメッシュに依存しない収束性を達成できるか？
• RQ2 移動種の濃度を排除することで、より良好な条件数を持つ線形系が得られ、収束が速くなるか？
• RQ3 プリコンディショニング付きニュートン・クリロフ法の反復回数およびCPU時間は、標準的手法と比べてどのように異なるか？
• RQ4 本手法は、現実的な反応輸送ベンチマークに適用されても、頑健性と効率性を維持できるか？
• RQ5 提案された定式化は、精度を保持しつつ、計算コストをどの程度低減できるか？

主な発見

• 非線形および線形反復回数が、さまざまなメッシュサイズに対してほぼ一定に保たれ、メッシュに依存しない収束性が確認された。
• ブロックガウス＝ザイデルおよびブロックヤコビプリコンディショナーは、プリコンディショニングなしの場合と比較して、クリロフ部分空間反復回数を顕著に削減した。
• メッシュの細分化に伴っても、ヤコビ行列の値域（field of values）がゼロから離れて保たれ、メッシュに依存しない収束性が裏付けられた。
• 移動種濃度を排除した代替定式化は、他のプリコンディショニング戦略と比較してCPU時間を削減した。
• 本手法は、大規模な時間刻み幅でも収束を達成でき、オペレータ分割法で一般的に見られる質量誤差を回避した。
• MoMaSベンチマークにおける数値結果は、既存の研究データと強く一致し、本手法の精度を検証した。