脑动脉瘤血流的计算流体动力学模拟

Overview

资料来源:约瑟夫·穆斯凯特、维塔利·雷兹和克雷格·戈尔根，珀杜大学韦尔登生物医学工程学院，西拉斐特，印第安纳州

本视频的目的是描述基于患者或动物特异性血管的计算流体动态（CFD）模拟的最新进展。在这里，创建了基于主题的容器分割，并使用开源和商业工具的组合，在流程模型中确定了高分辨率数值解决方案。大量研究表明，血管内的血管内血液动力学条件影响动脉粥样硬化、动脉瘤和其他周围动脉疾病的发展和进展;同时，在体内很难直接测量宫内压力、壁剪切应力（WSS）和粒子停留时间（PRT）。

差价合约允许非侵入性评估此类变量。此外，CFD还用于模拟手术技术，为医生提供术后流动状况的更好预见性。磁共振成像（MRI）中的两种方法，即具有飞行时间（TOF-MRA）或对比度增强 MRA （CE-MRA）和相位对比（PC-MRI）的磁共振血管造影（MRA），使我们能够获得容器几何和时间解析的三维速度场分别。TOF-MRA 基于通过应用于成像体积的重复 RF 脉冲抑制来自静态组织的信号。信号来自不饱和旋转，随流动血液进入体积。CE-MRA 是一种更好的技术，用于成像具有复杂循环流的容器，因为它使用造影剂（如高岭土）来增加信号。

另外，PC-MRI 利用双极梯度生成与流体速度成正比的相移，从而提供时间解析的速度分布。虽然 PC-MRI 能够提供血流速度，但此方法的准确性受有限的时空分辨率和速度动态范围的影响。CFD 提供卓越的分辨率，并可以评估从高速喷射到在病变血管中观察到的缓慢循环涡流的速度范围。因此，即使CFD的可靠性取决于建模假设，它为高质量、全面地描述患者特定的流动场提供了可能，从而可以指导诊断和治疗。

Procedure

本教程的前身是创建患者特定的血管模型。在本演示中，使用模拟工具、3D 系统地理魔法设计 X 和 Altair HyperMesh 从 MRA 数据生成四面体体积网格。

1. 为模型生成容器中心线

打开 vmtk 启动器 python GUI。在 PypePad 中，键入:vmtkcenterlines -ifile [STL 文件保存到桌面。
选择"运行，全部运行"将数据加载到程序中。将打开一个新窗口，显示输入模型的说明和呈现。旋转模型并将诅咒器放置在每个入口位置。按空格键放置种子。
在所有入口放置种子后，按"Q"继续。对所有插座重复相同的种子放置。放置出口种子后，再次按"Q"，让程序运行。这将将中心线文件保存到桌面。

2. 可视化软件中的数据集设置

启动开源可视化工具 ParaView（此过程中使用的版本 5.4.1）。
选择"文件""
Log in or to access full content. Learn more about your institution’s access to JoVE content here

Results

在本演示中，生成了脑动脉瘤的受试者特定模型，并使用 CFD 模拟流场。通过提供详细的流量特征和量化从成像数据中无法获得的造声能量，CFD 可用于增加低分辨率的 4D Flow MRI 数据。 图 1显示了 CFD 如何更完整地描述近壁、重新循环区域中的流量。

Log in or to access full content. Learn more about your institution’s access to JoVE content here

Application and Summary

此处描述的框架可用于执行特定于患者的 CFD 仿真。高分辨率网格用于插值低分辨率 4D 流 MRI 数据;这将隔离流量数据，并最大限度地减少与容器壁外部噪声相关的误差。通过使用基于患者的边界条件为入口和出口流量，仿真能够匹配与 MRI 成像的血动力学条件。

PC-MRI 的新方法能够显示更大的动态速度范围。然而，这受到患者扫描时间的严格限制。通常，患者数据以较低的分?...

Log in or to access full content. Learn more about your institution’s access to JoVE content here

Tags

本教程的前身是创建患者特定的血管模型。在本演示中，使用模拟工具、3D 系统地理魔法设计 X 和 Altair HyperMesh 从 MRA 数据生成四面体体积网格。

1. 为模型生成容器中心线

打开 vmtk 启动器 python GUI。在 PypePad 中，键入:vmtkcenterlines -ifile [STL 文件保存到桌面。
选择"运行，全部运行"将数据加载到程序中。将打开一个新窗口，显示输入模型的说明和呈现。旋转模型并将诅咒器放置在每个入口位置。按空格键放置种子。
在所有入口放置种子后，按"Q"继续。对所有插座重复相同的种子放置。放置出口种子后，再次按"Q"，让程序运行。这将将中心线文件保存到桌面。

2. 可视化软件中的数据集设置

启动开源可视化工具 ParaView（此过程中使用的版本 5.4.1）。
选择"文件""打开..."并查找以前创建的文件:患者特定的卷网格、中心线文件和 EnSight.case 文件。单击"确定"后，所有数据都应加载到界面中。
在左下角"属性"表中，选择"应用"。此命令将加载和读取用户在 ParaView 中加载或更改的所有信息。通过单击管道浏览器中的名称来突出显示卷网格以激活此选择。
同样，在"属性"表中，滚动并将不可变值更改为0.2 - 0.5之间的不全值。现在，中心线和几何渲染应可见。

3. 使用体积网格重新映射 4D 流 MRI 数据并删除杂色

从顶部菜单，选择过滤器，字母，重新采样与数据集。将打开一个新窗口。将源设置为体积网格，将"输入"设置为 EnSight.case 文件。设置这些后，请选择"确定"。
在"属性"表中，选择"应用"以应用筛选器。
与之前一样，突出显示新的重新采样与数据集* 名称以激活它。减少此新呈现的不恰量，如前面所述。此外，在顶部菜单中将中心线从"曲面"更改为"点"。

4. 确定入口和出口流量边界条件

在界面的右侧，旁边要最大化和最小化渲染选项，请选择"创建视图"工具（垂直线为正方形）。选择"铺张视图"选项。
从"显示"下拉框中，选择中心线文件。一次只能选择一种类型的文件。循环访问数据，选择各种点以识别每个入口和出口内的位置。
现在，使用点下方的"展图视图"计算靠近（4.2）中相同位置的两个点之间的法线矢量。
找到每个入口和出口位置的法向量后，选择"过滤器"、"字母表"和"切片"。请务必事先激活"重新采样与数据集*"。
切片筛选器需要显示在来自Resample 与Dataset*的重新采样数据集中的新分支下方。在"属性"表中，将平面原点设置为用于计算法向量的两个点之一的相同 XYZ 点位置。使用（4.3）中的法线矢量填充"正态"值。选择"应用"。
突出显示/激活刚刚创建的Slice+筛选器，然后选择"筛选器、字母表"和"曲面流"，然后应用。在管道浏览器中激活新的SurfaceFlow+项，并应用"筛选器、字母表、组时间步骤"应用。
在展图视图中，打开组时间步*数据。通过复制和粘贴或使用导出电子表格将这些数据导出到 Microsoft Excel。
在 Excel 中，计算对应于每个出口处的流速与总进气流率的比率的加权值。由于 4D Flow MRI 数据固有的噪声和误差，请识别最小（通常数据不太可靠）的容器，以"打开"以确保质量的保持。
在CFD仿真软件中，使用读瞬态表命令导入瞬态流波形;因此，以在线教程中描述的兼容.txt 格式保存入口流数据。

5. 设置差价合约模拟

打开差价合约模拟软件。在这里，我们使用ANSYS Fluent（本过程中描述的18.1版本作为默认版本）。选择"文件、读取和大小写..."，然后打开以前在 ParaView 中使用的体积网格 .cas 文件。显示网格（此过程使用使用 Altair HyperMesh 生成的 .cas 文件），通过选择"显示..."显示。
缩放几何体以确保模型的正确物理尺寸非常重要。选择"缩放..."，并应用特定情况所需的任何单位转换，然后关闭。
选择"材质，创建/编辑"以输入血液的材料属性。本教程分别使用1060 kg/s和0.0035 kg/ms的生理相关值来表示密度和粘度。
通过规定质量流或速度流率作为每个入口的时间函数来设置瞬态流边界条件。使用从 4D Flow MRI 测量中获得的波形来规定入口边界条件。出口点在（4.8）中找到的加权值。
在"解，方法"下，设置用于纳维埃-斯托克斯方程空间和时间离散化的数值方案。对于此过程，使用偶联，它启用全压速度耦合、最小平方单元（梯度）、压力二阶方案、动量方程的三阶 MUSCL方案以及用于及时离散的二阶隐式方案。确保左上角的时间参数已设置为瞬态。
在"解决方案、初始化"下，选择"标准初始化"。将所有初始值设置为0后，选择"初始化"。现在程序设置为运行。指定一个解决方案文件夹，以在计算活动下每个自动保存每个（时间步骤）保存结果。
在最后的步骤中，在"运行计算"下设置时间步长大小。使用（4.7）中的 Excel 边界条件数据来确定此值。缩短时间步长有助于收敛并提高数值解的准确性，同时增加解数时间。最好至少运行三个完整的心脏周期的模拟，以消除初始瞬变的影响。
最后，设置300 - 500之间的每个时间步长的最大迭代次数。一旦达到收敛，软件将在每个时间步长自动停止迭代，然后进入以下时间步长。通过运行具有平均速度值的稳定流模拟，然后将结果用作脉动流模拟的初始条件，可以改进收敛。选择准备运行解算器时的"计算"。
软件将运行每次迭代，直到实现收敛或最大迭代导致迭代继续。这些文件将自动保存在位置从（5.5），解决方案数据可以可视化在ANSYSCFD-Post或ParaView软件。

PLAYLIST

Loading...

版权所属 © 2025 MyJoVE 公司版权所有，本公司不涉及任何医疗业务和医疗服务。