血流动力学数值计算和模拟.pdf

摘要 摘要 动脉粥样硬化等血管疾病成为现代社会人类健康最主要威胁。研究动脉 血流运动是血流动力学主要任务之一，这对人们更好理解并分析血管疾病的 发生与发展十分重要。已有研究发现，局部血流动力学特性的改变和血管复杂 流场与血管病变之间关系十分密切。血管病变常常发生在分叉、弯曲等具有 特殊几何形状的位置附近。血液是由多种血细胞和血浆组成的悬浮液，其流 变性质为非牛顿流体。另外，血管壁具有弹性，动脉血流在心脏搏动影响下， 呈现明显脉动流特点。这些特征使得动脉血流过程是一个复杂的瞬态流．固耦 合系统。因此，要深入研究血流过程，单纯的理论或实验方法都会有较大局限 性。 随着计算流体动力学的不断发展和成熟，基于流体计算的数值模拟技术 与CAD技术或三维数字成像技术相结合，成为当前研究复杂血流过程的主要 手段。要进行数值模拟，数值求解不可压牛顿或非牛顿流体流动是中心问题 之一。不可压牛顿流体流动控制方程是不可压Naiver-Stokes方程。与之相比， 不可压非牛顿流动控制方程只是动量方程中的扩散项替换为粘性应力项。从 数值计算角度看，粘性应力项的出现并不会对数值离散构成本质困难，因此 影方法是应用最多的两种不可压流动算法。由于动脉血流脉动性的特点，更 适合稳态流动计算的SIMPLE方法没有投影方法更有效率。在与有限体积或 有限元方法结合求解血流动力学问题时，投影方法更有优势。另外，(病变)血 管复杂几何要求数值求解方法应建立在具有良好区域适应性的无结构网格上， 有限体积和有限元成为两种主要选择。在数值求解不可压NaiVer．Stokes方程 时，若用有限体积离散动量方法，可以自然地应用更多优良无振荡对流项数 值离散格式，这是其优点。若用Galerkin有限元方法离散动量方程，则数值解 有可能出现数值振荡，这需要结合流线扩散有限元方法(SUPG)、最小二乘有 限元(GLS)等方法来克服，但如此会使得数值实现过程变得更加复杂。离散动 量方程时，有限体积方法是更好的选择。然而，无结构网格有限体积投影方法 在离散压力Piosson方程时，仍然是使用有限体积方法，这并不是其优势所在， I 摘要 而恰恰是有限元方法的长处。发展一种结合两种方法优势的不可压流动数值 方法，并将其推广到血流动力学计算方面，具有现实意义。另外，在离散动量 方程对流项时，使用无振荡高阶格式会避免数值解非物理振荡产生。在各种 无振荡格式中，M币高分辨率格式基于三节点模板来重构对流项界面速度，构 造思想明确易行，很适合推广实施于无结构网格。本文在前人研究的基础上， 给出一种M汴格式的入实施方案，在此方案下，M厚格式具有形式一致性，易 于分析与代码实现。我们定义了新的无结构网格下正规变量形式，使得M汗格 式无结构网格上的计算更具效率。 第一章简要总结血流动力学所研究的基本问题，对血流动力学数值计算 方法，主要是不可压牛顿月}牛顿流动的数值计算进行了回顾。 第二章给出血流动力学中牛顿和非牛顿流动的控制方程，包括连续方程、 动量方程和能量方程，列出常见的非牛顿流体的本构方程和控制方程相应的 初始和边界条件。 第三章介绍对流扩散方程M，F格式及构造NVF高分辨率格式的CBC准则。 提出一种NvF格式的A型实施方案；定义一种无结构网格下正规变量定义新形 式，并给出无结构网格下NVF格式的实现方案。 第四章首先回顾了不可压流动计算中各种速度压力耦合算法。详细总结 了各种投影方法，并将旋度型压力增量投影方法推广至不可压非牛顿流计算。 详细讨论了无结构网格下有限体积方法控制元构成方式及相应变量配置方 式。在此基础上，提出一种无结构网格顶点型有限体积有限元投影方法。利用 本方法数值计算了腔体流动、平板间流动和后台阶流动(包括牛顿和非牛顿流 体流动，稳态流动和脉动流动)等典型验证算例。 第五章利用无结构投影方法，对二维血流动力学典型模型，如(多重)狭 窄、分叉等，进行数值计算。根据数值结果分析了与心血管疾病紧密相关血流 动力学因素，如壁面剪切应力，压力降等。 第六章总结本文主要内容，展望血流动力学数值计算和模拟的发展，并提 出进一步要解决的一些问题。 格式，有限体积方法，有限元方法，投影方法，无结构网格，非交错型变量配 置