[论文解读] 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.

研究动机与目标

• 建立一个可靠、可复现且高保真的患者特异性血流模拟数学与计算框架。
• 利用FDA喷嘴基准模型的实验数据，验证开源Feel++库中高阶有限元方法（FEM）求解器的性能。
• 评估网格细化、多项式阶次（P2, P3）以及几何逼近（G1, G2）对复杂流态下模拟精度的影响。
• 在不同雷诺数下评估关键验证指标$E_z$和$E_Q$，以量化与实验数据的一致性。
• 推动医疗设备开发中CFD的可靠性与可复现性提升。

提出的方法

• 采用稳定化Galerkin有限元法求解守恒形式的不可压缩非定常Navier-Stokes方程。
• 采用广义Taylor-Hood有限元对偶$P_{N+1}P_NG_{k_{\text{geo}}}$对速度、压力及几何逼近进行处理。
• 使用基于PETSc预条件的开源Feel++库进行高效并行求解，采用域分解法。
• 利用Gmsh生成结构化、各向异性的四面体网格，关键区域（如喉部和壁面边界）进行局部加密。
• 初始条件设定为$\mathbf{u} = \mathbf{0}$，入口边界$\Gamma_{\text{inlet}}$采用Poiseuille速度分布，壁面$\Gamma_{\text{wall}}$采用无滑移边界条件，出口$\Gamma_{\text{outlet}}$采用自由流出边界条件。
• 时间推进采用$\Delta t = 10^{-3}$ s（Re = 500）和$\Delta t = 10^{-4}$ s（Re = 2000和3500），以准确捕捉瞬态行为。

实验结果

研究问题

• 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网格在两处位置达到约10%，表明其对网格分辨率较为敏感。
• $P_2P_1G_2$逼近在最粗网格上即可实现$E_Q < 2\%$，而$P_3P_2G_1$在M0网格上未能提升精度，表明在粗网格中几何逼近的重要性超过速度阶次。
• Re = 2000模拟中，数值射流破裂点位于实验结果的下游，可能归因于数值积分精度不足或边界条件调优问题。
• 本研究证实，结合合理网格划分与几何逼近的高阶FEM可实现复杂血流动力学中高度准确且可复现的CFD模拟。