颈动脉模型中非定常流动.pdfVIP

应 用 力 学 进展 颈动脉模型中的非定常流动 于晓红，) 章梓雄 (香港大学机械工程系，香港) 摘要 心血管疾病成为危害世界人民健康的严重疾病之一，作为高发部位之一的颈动脉成为近几 十年来研究人员们竞相研究的对象利用数值模拟方法对颈动脉模型中的非定常流进行研究，通过 描述颈动脉中特定位置的壁面剪切应力随时间变化的规律，特征时刻内流线在颈动脉中的分布、及 颈动脉窦截面内的速度分布等，展现了血液在颈动脉中的流动特点，说明动脉粥样硬化斑块易在内 颈动脉出现的血液动力学机制． 关键词 数值模拟，颈动脉，非定常流动，动脉粥样硬化壁面剪切应力 1研究背景及目的 无论是发达国家还是发展中国家，心血管疾病已经成为威胁世界人民健康的严重疾病之一． 经过多年的研究表明心血管疾病的发生部位并不是随机的，而是经常出现在几个特定的部位， 如，颈动脉，冠状动脉，腹动脉，股动脉等，这种位置上惊人的确定性一定与某种引发疾病产生 的力学机制有关． 颈动脉，作为向脑部供血的主要通道，如果出现狭窄或拴塞将直接影响到对脑部的供血，与 缺血性中风的出现紧密相关．深入理解颈动脉内血液流动的特性，将有助于预测动脉粥样硬化 及栓塞在颈动脉中的产生与发展．颈动脉是具有分岔结构的血管，为实验测量带来了一定的难 度．并且通常的实验只能测量血管模型内的定常流动，对于脉动特性为显著特性的血液流动， 实验手段已显出了其局限性．随着计算机与计算流体力学(CFD)的发展，利用计算的方法模拟 复杂血管内部的脉动流成为可能，从而为深入的了解复杂血管内部流动提供了一种有效的非介 入性的手段． 本文将利用计算流体力学的方法模拟血液在颈动脉模型中的脉动流，在深入理解血液在人 体颈动脉中流动的同时阐明动脉粥样硬化在颈动脉中产生的血液动力学机制． 2计算模型及方法 2．1计算模型 真实颈动脉是具有分岔，两分支不处于同一平面，且他们与主支之间的夹角也因人而异的 血管，这些特点给研究工作带来了一定的难度，尤其是对颈动脉内流的实验研究，更是一个考 验．近年来已经有人通过磁共振的方法得到颈动脉的真实图像，并且利用计算流体力学的方法 模拟了其中的流动n但真实的颈动脉也是千差万别的，所以即使使用真实模型进行的模拟也 只能给整个问题一个定性化的说明．基于上述原因，本文采用了Perktold等口J给出的简化的颈 动脉模型来进行数值模拟(如图1)，对颈动脉的真实流动做定性化的研究．这里采用的颈主动脉 11 hkuhk Email：yuxh@hkusua 于晓红等：颈动脉模型中的非定常流动 的直径d=6．2 mm，内颈动脉的直径为o．7d，外颈动脉的直径为0．59d，颈动脉窦的直径为1．06d 内颈动脉与颈主动脉之间的夹角为15。，外颈动脉与颈主动脉之间的夹角为25。． 内颈动脉 颁主动脉 图1文中计算使用的颈动脉模型 2．2基本假设、边界条件及求僻方法 假设管内流动的是各向同性、不可压的牛顿流体，这些假设对于直径大于2mm的血管 都是适用的，Perktold等【2】也证明了非牛顿流性质对流动的影响只是局部的．同时，文中 还假设管壁是刚性的，因为管壁的弹性对 管内流动的影响也是有限的”j．控制方 程使用的是不可压非定常形式的N—S方程 组，并使用虚拟可压性的非定常公式求解 该方程组，在这个方法中，虚拟时间项作 为独立的变量出现，在每一个真实的时间 步长上方程的解是由这个虚拟时间迭代产 生． 壁面边界考虑为无滑移和无渗透的边 界条件，颈主动脉的入口处输入生理上较 为真实的脉动的流量[4l(图2)，分配到内 颈动脉的流量为总流量的o65，分配到外 2mm，平均流量为5．6ml／s， 动脉的直径，，是心跳的频率．对于文中使用的流量波形，在d=6 15文中使用结构网格对计算空 血液的运动粘