[論文レビュー] High order finite element simulations for fluid dynamics validated by experimental data from the fda benchmark nozzle model

Vincent Chabannes, Christophe Prud'Homme+2 more|arXiv (Cornell University)|Jan 9, 2017

Cardiovascular Function and Risk Factors被引用数 3

ひとこと要約

本論文は、オープンソースのFeel++ FEMライブラリを用いて、FDAベンチマークノズルにおける血液流れの高次有限要素シミュレーションを提示し、レイノルズ数Re = 500、2000、3500の層流、遷移流、乱流領域において実験データと照合する。軸方向速度および圧力差に関して、実験データと優れた一致を示し、$E_z$および$E_Q$の妥当性指標は、特に高次幾何と速度-圧力離散化を用いることで、高い精度と頑健性を示している。

ABSTRACT

The objective of the present work is to construct a sound mathematical, numerical and computational framework relevant to blood flow simulations and to assess it through a careful validation against experimental data. We perform simulations of a benchmark proposed by the FDA for fluid flow in an idealized medical device, under different flow regimes. The results are evaluated using metrics proposed in the literature and the findings are in very good agreement with the validation experiment.

研究の動機と目的

• 患者特異的血液流れのシミュレーションのための信頼性が高く、再現可能で高精度な数学的・計算的フレームワークを確立すること。
• オープンソースのFeel++ライブラリに実装された高次有限要素法（FEM）ソルバを、FDAのベンチマークノズルモデルからの実験データと照合すること。
• メッシュの細分化、多項式次数（P2、P3）、幾何近似（G1、G2）が、複雑な流れ領域におけるシミュレーション精度に与える影響を評価すること。
• 異なるレイノルズ数における主要な妥当性指標$E_z$および$E_Q$を評価し、実験データとの一致度を定量化すること。
• 医療機器開発におけるCFDの信頼性と再現可能性を向上させることに貢献すること。

提案手法

• 安定化ガレルキン有限要素法を用いて、保存形式の非圧縮性非定常ナビエ＝ストークス方程式を解く。
• 速度、圧力、幾何近似に一般化されたテイラー＝フッド有限要素ペア$P_{N+1}P_NG_{k_{\text{geo}}}$を採用する。
• PETScベースのプリコンディショニング（例：PCD）および領域分割法を用いて、効率的な並列解法を実装したオープンソースのFeel++ライブラリを使用する。
• Gmshを用いて構造的で非等方的な四面体メッシュを生成し、喉部や壁境界などの重要な領域を細分化する。
• 初期条件として$\mathbf{u} = \mathbf{0}$、$\Gamma_{\text{inlet}}$でポアウイユ入力プロファイル、$\Gamma_{\text{wall}}$でスリップなし境界条件、$\Gamma_{\text{outlet}}$で自由流出境界条件を適用する。
• 一時的挙動を解像するため、Re = 500では$\Delta t = 10^{-3}$ s、Re = 2000および3500では$\Delta t = 10^{-4}$ sの時間離散化を実装する。

実験結果

リサーチクエスチョン

• RQ1 高次有限要素法は、FDAノズルベンチマークにおいて、実験的測定値の軸方向速度および圧力プロファイルをどの程度正確に再現できるか？
• RQ2 メッシュの細分化および多項式次数（P2、P3）が、シミュレーション指標$E_z$および$E_Q$の収束性と精度に与える影響は何か？
• RQ3 高次幾何近似（G1対G2）は、特に遷移流および乱流領域において、実験データとの一致度にどのように影響するか？
• RQ4 Re = 2000のシミュレーションにおいて、数値的ジェット崩壊位置が実験よりもより下流に位置する理由は何か？
• RQ5 オープンソースのFeel++ FEMフレームワークは、血行動態シミュレーションにおいて、実験ベンチマークと同等の精度・頑健性・再現性を達成できるか？

主な発見

• Re = 500の条件下では、すべてのメッシュ細分化および多項式次数において、正規化された軸方向速度および圧力差が実験データと優れた一致を示した。
• 妥当性指標$E_z$は、すべてのシミュレーションで一貫して低く（< 0.05）維持され、計算された軸方向速度と実験値との強い一致を示している。
• 質量保存誤差$E_Q$は、粗いメッシュ（M0）を除きすべてのケースで2%未満であり、メッシュ解像度への感受性が示された。M0では2か所で約10%に達した。
• $P_2P_1G_2$近似では粗いメッシュでも$E_Q < 2\%$を達成したが、$P_3P_2G_1$ではM0で精度の向上が見られず、粗いグリッドでは幾何近似の方が速度次数よりも重要であることが示唆された。
• Re = 2000のシミュレーションにおける数値的ジェット崩壊位置は、実験よりも下流に位置しており、数値積分の精度や境界条件の調整の問題が関係している可能性がある。
• 本研究は、適切なメッシュ生成と幾何近似を伴う高次FEMが、複雑な血行動態流れにおける高精度で再現可能なCFDシミュレーションを可能にすることを確認した。