[論文レビュー] Numerical simulation of parabolic moving and growing interface problems using small mesh deformation

Ulrich Langer, Huidong Yang|arXiv (Cornell University)|Jul 31, 2015

Advanced Numerical Methods in Computational Mathematics参考文献 21被引用数 4

ひとこと要約

本稿では、3次元の放物型移動・成長界面問題を解くための、新しいALEベースの数値法を提案する。この手法は、四面体および八面体要素から成るマクロ要素に局所的な小規模なメッシュ変形を適用する。界面捕捉と代数的多重グリッド（AMG）ソルバーを組み合わせることで、大規模な変形や成長に対しても安定した性能を達成し、数値実験ではすべてのスパンでAMGプリコンディショニング付き共役勾配法（CG）が最適なCPU効率を示した。

ABSTRACT

In this work, we develop a cutting method for solving problems with moving and growing interfaces in 3D. This new method is able to resolve large displacement or deformation of immersed objects by combining the Arbitrary Lagrangian-Eulerian method with only small local mesh deformation defined on the reference domain, that is decomposed into the macro-elements. The linear system of algebraic equations arising after the temporal and spatial discretizations of a model parabolic interface heat-conduction-like problem with vector-valued functions is solved by either an all-at-once or a segregated algebraic multigrid method.

研究の動機と目的

• 3次元において大規模な変位や変形を伴う界面問題を安定的かつ効率的にシミュレートするための数値法の開発。
• 大規模な変形によって生じるメッシュの著しい歪みに対処できない従来のALE法の限界を克服すること。
• 参考マクロ要素の分割に基づく、局所的な小規模なメッシュ変形のみを用いることで、高い精度と効率を維持すること。
• 得られる線形方程式系に対する、すべてのスパンで処理する（all-at-once）と分離処理する（segregated）代数的多重グリッド（AMG）ソルバーの性能を比較すること。
• ベクトル値の熱伝導に類似した問題を用いて、剛体移動と界面成長の両方を、複雑な3次元領域内でも安定してシミュレートできることを示すこと。

提案手法

• 本手法は、マクロ要素に分解された参考領域を用い、各マクロ要素は4つの四面体と1つの八面体から構成され、局所的なメッシュ変形が可能となる。
• 固定メッシュのALEアプローチと、移動・成長対象とバックグラウンドの四面体要素との間の事前計算済み交差を用いた界面捕捉を組み合わせる。
• モデル問題は、領域Ω₁とΩ₂で不連続な拡散係数a₁とa₂を有する、ベクトル値の解uが満たす放物型界面方程式である。
• 時間的および空間的離散化により、大規模な線形方程式系が得られ、これはすべてのスパンで処理するか、分離処理する代数的多重グリッド（AMG）法で解かれる。
• CGおよびGMRESソルバーにAMGプリコンディショナーを適用し、モノリシックおよびシュール補行列系の両方で使用する。
• マクロ要素上での一次多項式有限要素基底関数を活用することで、追加の自由度を導入せずに済ませる。

実験結果

リサーチクエスチョン

• RQ1 マクロ要素上での小規模な局所的メッシュ変形が、大規模な変位を伴う3次元放物型界面問題において、メッシュの品質と精度を維持できるか。
• RQ2 すべてのスパンで処理するAMGプリコンディショニング付きCG法と、分離処理するAMGプリコンディショニング付きCGおよびGMRES法の間で、収束性および計算コストの観点から性能に差があるか。
• RQ3 提案された界面捕捉法が、追加の自由度を必要とせずに、成長界面を効果的に扱えるか。
• RQ4 時間依存の界面問題から生じる線形方程式系を解く際、CG、GMRES、および分離処理アプローチの相対的な効率はいかほどか。
• RQ5 本手法が、複雑な3次元領域内での浸漬物体の剛体移動と等方的成長を、両方とも安定してシミュレートできるか。

主な発見

• すべてのスパンで処理するAMGプリコンディショニング付きCG法が、移動および成長界面シミュレーションの両方において、CPU時間の観点で最良の計算効率を達成した。
• AMGプリコンディショニング付きGMRES法が、反復回数の観点で最良の性能を示し、1反復あたりの収束が速いことを示した。
• 分離処理法は反復回数の観点ですべてのスパンで処理するアプローチと同等の性能を示したが、CPU時間がやや高かった。
• 移動物体のケースでは、GMRES法が最も少ない反復回数（最良の収束）を達成した一方、CG法が最も低いCPU時間を記録した。
• 成長する球体のケースでも、同じ性能順序が成立：GMRESが最良の反復回数、CGが最良のCPU時間、分離処理法は中程度の効率を示した。
• すべてのシミュレーションにおいて相対残差誤差が1e-9未満に保たれ、線形ソルバーの収束が安定していることを確認した。