多发性颅内动脉瘤的血液动力学数值模拟的研究.pdfVIP

多发性颅内动脉瘤的血液动力学数值模拟研究1) *,2) * + 陈家亮 ，丁光宏 ，杨新健 * + （复旦大学力学与工程科学系，上海 200433 ） （北京市神经外科研究所，北京 100050) 摘要：血液动力学因素被认为与颅内动脉瘤的整个自然病程有着密切的关系。本研究利用血液动力学数值模拟 方法分析多发性动脉瘤患者的局部血液流动情况，着重针对每一例多发性动脉瘤患者,从血液动力学的角度分析 不同位置动脉瘤之间的相互作用。我们选取了四例多发性动脉瘤患者作为研究对象，每例病人在同一根血管上 有两个动脉瘤邻近生长。结果显示对于每一例多发性动脉瘤患者,不同位置上的动脉瘤具有各自不同的血液流动 形态和壁面切应力分布,而这主要是受到相邻载瘤动脉弯曲引起的二次流动所影响。另外我们还发现对于其中三 例多发性动脉瘤患者，虽然两个动脉瘤在同一根血管上邻近生长，但是从血液动力学的角度分析这两个动脉瘤 之间是没有明显的相互作用，即下游动脉瘤的血液流动形态和壁面切应力分布不会受到上游动脉瘤存在的影 响。这一分析结果也就从血液动力学的角度说明在临床上有选择地保留一些危险系数比较小的伴生动脉瘤不会 影响到栓塞动脉瘤的治疗效果,因为这些小的伴生动脉瘤不会明显改变载瘤动脉中的血液流动情况，从而也就不 会影响栓塞动脉瘤残颈区域的血液流动。 关键词：多发性颅内动脉瘤；血液动力学数值模拟；壁面切应力；血液流动形态 引 言 颅内动脉瘤是指颅内动脉的异常性扩张，主要发生在大脑动脉的分叉或者弯曲较为剧烈处及其邻近区 [1] 域 。颅内动脉瘤破裂导致蛛网膜下腔出血占脑血管意外病人中的第三位，其中多发性颅内动脉瘤患者 [2] 占整个颅内动脉瘤破裂患者人数的35%左右 。对于此类病人,通常在临床介入手术过程中会根据病人的 实际情况有选择地栓塞脉瘤，而保留一些临床认为危险系数比较小的伴生动脉瘤。这些没有栓塞的伴生 动脉瘤是否是合理的,是否会影响栓塞动脉瘤术后治疗效果这一问题在临床上还没有定论。虽然目前颅内 动脉瘤长期生长过程以及破裂机制仍然没有完全清楚，但在已知的各种影响因素中，血液动力学因素被 认为与颅内动脉瘤的整个自然病程有着密切的关系，并且在临床优化动脉瘤手术方案以及后期评估动脉 [3-6] 瘤治疗效果的过程中也起到了相当重要的作用 。本研究利用血液动力学数值模拟方法分析多发性动脉 瘤患者的局部血液流动情况，着重针对每一例多发性动脉瘤患者,从血液动力学的角度分析不同位置动脉 瘤之间的相互作用，从而为临床有选择性地治疗动脉瘤提供一定参考。 1) 国家自然基金(和上海市重点学科基金 （B 112）资助项目 2) E-mail ：marshalfdu@ 1 1 多发性颅内动脉瘤血液动力学计算模型建立方法 颅内动脉瘤的三维几何形态数据来自临床病人的三维螺旋血管造影原始数据的图像重建,保存为 STL 格式文件。然后将已获得的颅内动脉瘤三维几何模型数据STL 文件导入商业化计算流体力学网格划 分软件ICEM CFD 中进行流体计算所需网格的建立。在计算流体网格划分这一步中最重要的就是确定网 格类型和网格大小，在本研究中所有计算模型都采用混合网格类型 （hybrid mesh ）,在ICEM CFD 中具体 称为四面体-三棱柱混合网格，即在大部分区域用四面体单元进行划分而在靠近血管壁和瘤壁的区域用三 棱柱网格进行划分。对于四面体单元区域，我们统一定义最大网格尺寸约为血管平均直径的 10%.对于三 棱柱单元，需要定义这个区域的大小和首层网格距血管壁的距离，对于所有动脉瘤模型我们统一定义模 型三棱柱单元区域的尺寸为该模型中最小血管直径的25%，首层网格距血管壁的距离，我们统一定义为 0.02mm，网格类型如图1 所示。然后将ICEM CFD 产生的四面体-三棱柱混合网格导入到Ansys CFX 中 进行流体计算文件的配置，这其中包括血液模