[論文レビュー] One-Dimensional Navier-Stokes Finite Element Flow Model

Taha Sochi|arXiv (Cornell University)|Apr 8, 2013

Advanced Numerical Methods in Computational Mathematics参考文献 31被引用数 15

ひとこと要約

本技術報告では、弱形式とガレルキン離散化を用いて、弾性血管および接続されたネットワーク内のニュートン流体の流れをシミュレートする1次元ナビエ＝ストークス有限要素モデルを提示する。本モデルは、さまざまな境界条件および初期条件下でも妥当な収束性と頑健性を示し、効率的で低コストの血液流れのシミュレーションを可能にする。

ABSTRACT

This technical report documents the theoretical, computational, and practical aspects of the one-dimensional Navier-Stokes finite element flow model. The document is particularly useful to those who are interested in implementing, validating and utilizing this relatively-simple and widely-used model.

研究の動機と目的

• 大血管におけるニュートン流体の流れを対象とした、計算的に効率的で数学的に頑健な1次元ナビエ＝ストークス方程式の有限要素実装を開発すること。
• 血管壁の力学に基づく圧力-面積の構成関係を組み合わせることで、可撓性のある血管へのモデルの拡張を図ること。
• 単一血管および分岐ネットワークにおける時間に依存しないおよび時間に依存する流れをシミュレートする包括的なフレームワークを提供すること。
• 解析的解と収束解析を通じて、数値実装の妥当性を検証すること。
• ハエモダイナミクスおよびバイオフラUIDメカニクス分野の研究者に対して、実装、境界条件、無次元化、解の収束に関する実用的ガイダンスを提供すること。

提案手法

• 空間的離散化にガレルキン有限要素法を用いて、1次元ナビエ＝ストークス方程式を弱形式で定式化する。
• 単一血管および分岐ネットワークにおける時間に依存しないおよび時間に依存する流れの両方に、有限要素法を適用する。
• 分岐ノードにおける質量保存およびエネルギー保存を確保するため、適合性および一致条件を課す。
• 数値的安定性の向上およびパラメータ比較の容易化を目的に、無次元化を実施する。
• 有限要素上の弱形式の残差の数値積分にガウス求積法を用いる。
• 血管壁の弾性（ヤング率、ポisson比）に基づく線形化された形で、圧力-面積関係を統合する。

実験結果

リサーチクエスチョン

• RQ1 時間に依存しないおよび時間に依存する両流れに対して、1次元ナビエ＝ストークス方程式を有限要素法で効果的に離散化する方法は何か？
• RQ2 収束性および物理的整合性を確保するための、最適な境界条件および初期条件は何か？
• RQ3 無次元化は、有限要素解の安定性および精度にどのように影響するか？
• RQ4 本モデルは、弾性血管内の流れに対して解析的解をどの程度正確に再現できるか？
• RQ5 本1次元手法を用いた分岐血管ネットワークのシミュレーションにおいて、主な実装上の課題は何か？

主な発見

• 時間に依存しないおよび時間に依存する両流れに対して、有限要素実装はメッシュの細分化に伴い解析的解に収束し、精度が向上することが確認された。
• 時間に依存する定式化は、時間に依存しない場合よりもより信頼性が高く、漸近的に定常状態解を近似可能である。
• 無次元化により数値的安定性が向上し、異なる血管幾何形状および材料特性間での一貫した比較が可能になった。
• 非一様なネットワークにおいてもモデルは精度を維持するが、血管間の急激な変化は収束を妨げる可能性がある。
• ヤング率およびポisson比を用いた血管壁弾性に基づく圧力-面積関係の組み込みにより、可撓性血管の正確なシミュレーションが可能になった。
• 本手法は計算的に効率的であり、モデルの有効範囲内で信頼性の高い結果を提供するため、対称的で大血管の流れに対して2次元／3次元モデルの代替手段として実用的である。