[論文レビュー] A robust conservative mixed finite element method for compressible flow on pipe networks

Herbert Egger|arXiv (Cornell University)|Sep 16, 2016

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

ひとこと要約

本稿では、圧縮流体のパイプネットワークにおける流れをシミュレートするための、頑健で保存的な混合有限要素法を提示する。変分式により全エネルギーと質量の保存が保証され、空間方向には適合ガラーキン離散化を、時間方向には質量を正確に保存し、エネルギーをわずかに散逸させる陰解法を採用することで、多様なネットワーク構造やパrameter設定に対しても適切な定式化と安定性が保証される。

ABSTRACT

We consider the numerical approximation of compressible flow in a pipe network. Appropriate coupling conditions are formulated that allow us to derive a variational characterization of solutions and to prove global balance laws for the conservation of mass and energy on the whole network. This variational principle, which is the basis of our further investigations, is amenable to a conforming Galerkin approximation by mixed finite elements. The resulting semi-discrete problems are well-posed and automatically inherit the global conservation laws for mass and energy from the continuous level. We also consider the subsequent discretization in time by a problem adapted implicit time stepping scheme which leads to conservation of mass and a slight dissipation of energy of the full discretization. The well-posedness of the fully discrete scheme is established and a fixed-point iteration is proposed for the solution of the nonlinear systems arising in every single time step. Some computational results are presented for illustration of our theoretical findings and for demonstration of the robustness and accuracy of the new method.

研究の動機と目的

• パイプネットワークにおける圧縮流体の流れにおいて、質量とエネルギーを全領域で保存する数値法の開発。
• 接続部における物理的整合性と質量・エネルギー保存を保証する結合条件の導出。
• 保存則を自然に組み込んだ変分原理の構築と、適合ガラーキン離散化を可能にする形式の定式化。
• 半離散および完全離散スキームの適切な定式化と保存性の確立。
• リーマン型の流れ問題と定常状態収束に関する数値実験を通じた、ロバスト性と精度の検証。

提案手法

• 圧縮流体の支配方程式に基づく変分式を導出し、質量と運動量のバランスに加え、圧力則 $ p = c\rho^\gamma $ を組み込む。
• パイプ接続部における結合条件を、質量フラックスと滞留エンタルピーの保存を満たす形で定式化し、粘性効果を含む。
• 空間離散化に混合有限要素法を適用し、適合ガラーキン近似を実現するとともに、半離散レベルでの正確な質量保存を保証する。
• 質量を正確に保存し、エネルギーに制御された数値的散逸を導入する陰時間ステッピングスキームを採用する。
• 各時間ステップで生じる非線形方程式系は、固定点反復法により解き、ロバスト性と収束性を確保する。
• エネルギー推定と安定性解析を用いて、完全離散スキームの適切な定式化を確立する。

実験結果

リサーチクエスチョン

• RQ1 パイプネットワークにおける圧縮流体の流れにおいて、質量とエネルギーの全領域保存を内蔵する変分式を導出できるか？
• RQ2 接続部における結合条件をどのように定式化すれば、物理的整合性と質量・エネルギー保存を保証できるか？
• RQ3 変分式に適合する混合有限要素法を適用した場合、離散レベルでも質量とエネルギー保存が保たれるか？
• RQ4 質量を保存し、エネルギー散逸を制御的に導入する陰時間離散化を設計できるか？
• RQ5 完全離散スキームは、モデルパrameter、メッシュの細分化、初期条件に対してどれほどロバストか？

主な発見

• 提案手法は離散レベルで正確に質量を保存するため、正しい定常状態解を得る上で極めて重要である。
• 陰時間離散化によりエネルギーがわずかに散逸され、物理的粘性散逸を模倣し、安定性が保証される。
• 半離散および完全離散スキームともに適切に定式化されており、メッシュと時間刻みの細分化に伴い収束が観察された。
• 数値実験では、初期条件が異なっても（不連続密度プロファイルを含む）ロバストな挙動を示し、数値的減衰により不連続性が迅速に滑らかに処理された。
• 微小な粘性係数や複雑なネットワーク構造に対しても、保存性と安定性が維持された。
• 固定点反復法はすべてのテストケースで信頼性高く収束し、非線形ソルバーのロバスト性が裏付けられた。