[論文レビュー] Overview of the Incompressible Navier-Stokes simulation capabilities in the MOOSE Framework

John W. Peterson, Alexander Lindsay+1 more|arXiv (Cornell University)|Oct 24, 2017

Advanced Numerical Methods in Computational Mathematics参考文献 41被引用数 3

ひとこと要約

本論文では、MOOSE有限要素フレームワーク内での非圧縮性ナビエ＝ストークス（INS）ソルバの実装と検証を提示している。安定化有限要素法（SUPG、PSPG、GLSを含む）を用いて、複雑な2次元および3次元幾何形状に対しても、正確で安定した解を得ている。主な貢献は、MOOSEに統合された堅牢で拡張可能かつ十分に検証されたINSモジュールであり、マルチフィジックス結合を可能にし、広範なレイノルズ数をカバーする。

ABSTRACT

The Multiphysics Object Oriented Simulation Environment (MOOSE) framework is a high-performance, open source, C++ finite element toolkit developed at Idaho National Laboratory. MOOSE was created with the aim of assisting domain scientists and engineers in creating customizable, high-quality tools for multiphysics simulations. While the core MOOSE framework itself does not contain code for simulating any particular physical application, it is distributed with a number of physics "modules" which are tailored to solving e.g. heat conduction, phase field, and solid/fluid mechanics problems. In this report, we describe the basic equations, finite element formulations, software implementation, and regression/verification tests currently available in MOOSE's navier_stokes module for solving the Incompressible Navier-Stokes (INS) equations.

研究の動機と目的

• MOOSE有限要素フレームワーク内に、安定的で正確かつ拡張可能な非圧縮性ナビエ＝ストークス（INS）ソルバを実装・検証すること。
• INSシミュレーションを、熱伝導、固体力学、相場など、他のMOOSE物理モジュールとシームレスに結合できるようにすること。
• 確立された有限要素安定化技法を用いて、広範なレイノルズ数および複雑な2次元／3次元幾何形状をカバーすること。
• コミュニティ開発、オープンソース、保守可能なINSモジュールをMOOSEフレームワークに直接統合し、コード変更を最小限に抑えること。
• 将来の高性能シミュレーションに向け、スケーラビリティの主な課題（プリコンディショニング、多重グリッド法など）を特定・解決すること。

提案手法

• INS方程式は、ストリームライン・アップウインド・ガレルキン（SUPG）、圧力安定化ガレルキン（PSPG）、ガレルキン最小二乗法（GLS）に基づく安定化有限要素法を用いて解かれる。
• メッシュ処理、有限要素空間、線形代数にはlibMeshライブラリが活用され、MOOSEはアプリケーション開発のための一貫性があり高水準なインターフェースを提供する。
• navier_stokesモジュールは、連続および非連続ガレルキン（DG）離散化をサポートしており、インフラ・サブスクリプション条件を満たす速度-圧力混合型式に重点を置いている。
• 検証には、製造された解決法（MMS）と、Re = 100およびRe = 200での球体周囲流れのようなベンチマーク問題が用いられている。
• コード品質と保守性を確保するため、継続的インテグレーションワークフローに基づいてソフトウェアが開発されている。
• フレームワークは、流体-構造連成のためのtensor_mechanicsや界面追跡のためのlevel_setなど、他のMOOSEモジュールとの結合をサポートしている。

実験結果

リサーチクエスチョン

• RQ1 確立された安定化技法を用いて、MOOSE有限要素フレームワーク内に安定的で正確な非圧縮性ナビエ＝ストークスソルバを実装する方法は何か？
• RQ2 MOOSEのnavier_stokesモジュールは、広範なレイノルズ数および複雑な幾何形状において、性能と正確性の特性をどのように示すか？
• RQ3 INSモジュールは、流体-構造連成や熱対流のようなマルチフィジックスシミュレーションにおいて、他のMOOSE物理モジュールと効果的に結合できるか？
• RQ4 特にプリコンディショニングと線形方程式系の解法効率に注目した場合、INSソルバを大規模問題にスケーリングする上で直面する主な課題は何か？
• RQ5 将来的に、乱流モデル、非ニュートン流体、低マッハ数モデルなどの拡張は可能で、科学的インパクトを広げるために必要とされるか？

主な発見

• navier_stokesモジュールは、Re = 100およびRe = 200での球体周囲流れを良好にシミュレートでき、再循環領域や境界層剥離といったベンチマークデータとよく一致している。
• 製造された解決法（MMS）による検証により、速度場および圧力場の最適収束率が確認され、実装の正確性が裏付けられている。
• ソフトウェアは、小規模から中規模の問題に対して本番運用可能であり、MOOSEと共に配布されており、新しいアプリケーションに統合するにはMakefileの変更のみで可能である。
• 実装は、適合および非適合有限要素離散化をサポートしており、乱流モデルや非ニュートン流体の今後の拡張に対しても拡張可能である。
• フィールド分割プリコンディショニングと幾何的多重グリッド法が、スケーラビリティ向上の最も有望な道筋と特定されたが、シュール補完の近似における現在の制限が広範な利用を妨げている。
• 共有解フィールドと一貫性のあるインターフェースを介して、熱伝導や固体力学との結合を含む複雑なマルチフィジックスシミュレーションが可能になっている。