颈动脉分叉个体实例模型数值模拟.pdfVIP

第二十一届全国水动力学研讨会暨第八届全国水动力学学术会议 暨两岸船舶与海洋工程水动力学研讨会文集 颈动脉分叉个体实例模型的数值模拟 1 1 2 1 1 王青 王炜哲 费智敏 刘应征 曹兆敏 1 上海交通大学机械与动力工程学院 2 上海交通大学医学院附属仁济医院 摘要：颈动脉分叉处是动脉粥样硬化的多发部位。大量研究表明，动脉粥样硬化的发生与该 处的血液动力学特性密切相关。而研究个体化的颈动脉分叉模型可以帮助获得个体更为真实 的颈动脉血液流动细节。本研究通过对CT 医学影像的三维重建，获得人体个体化真实颈动脉 分叉模型，并且利用计算流体力学（CFD）的研究方法进行数值模拟，分析个体化模型的血液 动力学特性。结果表明：在一个心动周期内，此个体化颈动脉分叉模型在颈动脉窦和外颈动 脉分叉入口外壁附近血液流动发生流动分离，产生旋涡等复杂流动。这些区域的外壁面的壁 面切应力较低，并随入口速度的脉动而振荡，这些特征的产生与颈动脉几何形态密切相关。 关键词：颈动脉分叉；血液动力学；计算流体力学；动脉粥样硬化 1 引言 动脉粥样硬化是严重威胁人类健康的心血管疾病，其发生发展的机理一直是生物医学工 程领域的研究热点之一，但至今仍没有彻底弄清。从形态学上看，动脉粥样硬化总是发生在 大动脉的弯曲部、分叉口附近这些局部部位。这些位置上的血液流动的动力学特性在动脉粥 样硬化的发生和发展过程起着决定性影响。Nerem[1]等的研究认为动脉粥样硬化多发的局部部 位血液流动复杂，流动具有多向性。根据Caro[2] 的理论，在低切应力区域，血脂质分子滞留时 间长，容易在该区域长生动脉粥样硬化斑块。目前普遍认为 Ku[3]等人提出的低切应力或者切 应力的周期性振荡是其主要动力学因素。颈动脉分叉处是动脉粥样硬化的多发部位，由动脉 粥样硬化引起的颈动脉狭窄是脑缺血和中风产生的主要原因。颈动脉分叉由于同时存在弯曲 和分叉，还有颈动脉窦，流动情况复杂，对研究这些局部位置的血液动力学特性有重要意义。 解剖学报告表指出，在颈动脉分叉血管中粥样硬化斑块优选区位于分支血管的外壁和侧壁， 内壁很少发生，如图 1 所示。不少学者通过建立颈动脉分叉模型来研究颈动脉形状和流动参 数对颈动脉流动特性的影响，分析其与动脉粥样硬化的关系[4] 。人体血管结构相当复杂并且存 在很大的个体差异，研究个体化颈动脉分叉模型可以获得相对个体而言更为真实的流动细节， 并可进一步探讨由于这些个体差异会引起怎样的流动差异, 以及这些差异在心血管疾病方面的 影响。 目前，还没有很好的实验手段能观察到真实血管内部流动细节，而数值模拟在研究流动 细节方面有很大优点。本研究通过医学影像技术三维重建个体化的颈动脉分叉模型，用计算 - 353 - 第二十一届全国水动力学研讨会暨第八届全国水动力学学术会议 暨两岸船舶与海洋工程水动力学研讨会文集 流体力学（CFD ）的方法进行数值模拟，研究真实的个体化颈动脉分叉处的血液动力学特性及 其与动脉粥样硬化的关系。从长远来看，这一方法可以为动脉粥样硬化的预防，早期干预和 制定个体化治疗方案提供可靠的血液动力学依据。 图1 动脉粥样硬化引起颈动脉狭窄的示意图 2 数值模拟方法 2.1 数学模型 在正常生理条件下，颈动脉内血液流动为层流。Xu[6]等人的研究表明把血液视为牛顿流 体或非牛顿流体对流动形态影响不大。Zhao[7]等人研究表明将血管壁视为刚性不影响整体的流 动特性及壁面切应力的分布。因此为了简化计算模型，将颈动脉血管视为无渗透、刚性血管, 并假设血液是在恒定温度下具有相同的密度和黏度的牛顿流体，内颈动脉（ICA ）和外颈动脉 （ECA ）流量分配比恒定为正常人体的颈动脉分叉流量比例7:3[5] 。血液流动为不可压缩非定 常层流，遵守质量守恒和动量守恒方程，即连续方程和N-S 方程。 连续方程： ∂u ∂v ∂w + + 0