[论文解读] One-Dimensional Navier-Stokes Finite Element Flow Model

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

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

一句话总结

本技术报告提出了一种用于模拟弹性血管及互联网络中牛顿流体流动的一维纳维-斯托克斯有限元模型，采用弱形式和伽辽金空间离散化。该方法可实现高效、低成本的大血管血流模拟，并在各种边界和初始条件下验证了其收敛性和鲁棒性。

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.

研究动机与目标

• 开发一种计算高效且数学上稳健的一维纳维-斯托克斯方程有限元实现方法，用于大血管中牛顿流体的流动模拟。
• 通过基于血管壁力学的压力-面积本构关系，将模型扩展至可扩张血管。
• 提供一个全面的框架，用于模拟单个血管及分支血管网络中的定常与非定常流动。
• 通过解析解和收敛性分析验证数值实现的正确性。
• 为血流动力学与生物流体力学领域的研究人员提供关于实现方法、边界条件、无量纲化及解收敛性的实用指导。

提出的方法

• 采用伽辽金有限元法对一维纳维-斯托克斯方程进行弱形式化，实现空间离散化。
• 将有限元方法应用于单个血管及分支血管网络中的定常与非定常流动。
• 在分支节点处施加相容性与匹配条件，以保证质量与能量守恒。
• 采用无量纲化以提高数值稳定性，并便于不同血管几何与材料属性之间的参数比较。
• 采用高斯积分法对有限元上的弱形式残差进行数值积分。
• 基于血管壁弹性（杨氏模量、泊松比）的线性化形式，实现压力-面积关系的积分。

实验结果

研究问题

• RQ1 如何通过有限元法对一维纳维-斯托克斯方程进行有效离散化，以处理定常与非定常流动？
• RQ2 为确保一维流动模拟中的收敛性与物理一致性，应采用何种最优边界与初始条件？
• RQ3 无量纲化如何影响有限元解的稳定性与精度？
• RQ4 该模型在多大程度上能准确再现弹性血管中流动的解析解？
• RQ5 使用该一维方法模拟分支血管网络时，面临的主要实现挑战是什么？

主要发现

• 有限元实现对定常与非定常流动均能成功收敛至解析解，且随着网格加密，精度进一步提升。
• 非定常格式的收敛性优于定常格式，可实现对稳态解的渐近逼近。
• 无量纲化显著提升了数值稳定性，并实现了不同血管几何与材料属性间的稳定比较。
• 即使在非均匀网络中，模型仍能保持较高精度，但血管间突变过渡可能影响收敛性。
• 通过引入基于血管壁弹性（杨氏模量与泊松比）的压力-面积关系，模型可准确模拟可扩张血管。
• 该方法计算效率高，在模型有效范围内提供可靠结果，是二维/三维模型在对称大血管流动模拟中的可行替代方案。