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

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的可靠性取决于建模假设，它为高质量、全面地描述患者特定的流动场提供了可能，从而可以指导诊断和治疗。

Principles

TOF-MRA、CE-MRA 和 PC-MRI 通常用作 CFD 仿真的输入几何和流边界条件。如上所述，每个主体测量容器几何体和流入边界条件（通过横截面的速度剖面）。对于本研究中包含的数据，TOF-MRA 分辨率为 0.26 x 0.26 x 0.50 mm，而 PC-MRI 分辨率为 1.00 x 1.00 x 1.20 mm。4D流MRI序列用于通过心脏周期获得三维速度分布。TOF 数据使用各种工具自动进行伪分段。图像分辨率（即体素的大小）直接影响几何体生成的模型的质量。4D Flow MRI Equation 1 根据以下方程使用相移 Equation 2 确定每个体素的血液速度:

Equation 3 (1)

Equation 4 (2)

测量的相位和速度取决于梯度 Equation 5 场、陀螺比 Equation 6 、自旋 Equation 7 的初始位置、自旋速度和 Equation 8 自旋加速度。 Equation 9 在初始化 MRI 扫描时定义磁场和材料常量。4D 流 MRI 编码在三个正交方向，以获得三维速度场。然后，可以为每个患者或动物特定病例生成 3D 模型。过程部分中详细介绍的方法将通过数值求解 Navier-Stokes 方程来引入 CFD 模拟，这些方程被概括为:

Equation 10 ₍₃₎

哪里 Equation 11 是密度， Equation 12 是流速，p是压力，而亩是流动的动态粘度。

Procedure

本教程的前身是创建患者特定的血管模型。在本演示中，使用模拟工具、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软件。

Results

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

图1:A）容器几何体中4D流MRI数据的可视化。B） CFD模拟结果的可视化。一般来说，CFD流线型使对脑动脉瘤内血流模式的理解更加全面。

图1显示CFD结果与体内4D流MRI一致。图 1 （A）显示了动脉瘤区域内复杂的循环流动模式，即动脉的气球状扩张，通过 4D Flow MRI 检测到。然而，病变顶部和底部的停滞流动区域没有填充流线。这是因为这些区域的信号噪声比较低。CFD模拟流量如图 1 （B）所示，可提供更高分辨率的速度场，尤其是在容器壁附近。因此，CFD 模型能够提供临床相关、流派生指标（如压力、WSS 和 PRT）的更高精度估计值，可用于预测动脉瘤疾病进展。

此外，CFD 模拟可用于对替代治疗方案产生的术后流动条件进行建模。例如，图 2 （A）和（B）比较流经同一容器，流入率不同。通过规定不同的边界条件，如模拟无流动的血管闭塞，显示各种手术治疗后的流量。

图2:A）右前脑动脉（ACA）手术剪裁的模拟。B）模拟左ACA的手术剪裁。为简单起见，此图在未修改的入口处保持术前流入率;在现实中，流速会增加在开放容器，以补偿。C）正常血流量规定此模型的进气条件。来自 4D Flow MRI 的患者数据为流程模式的真实可视化提供了入口条件。

模拟各种手术治疗产生的术后流动场的能力是CFD模型的一个重要优势。通过应用逼真的、特定于患者的几何形状和流入数据，可以演示不同的治疗方案，为医生提供有关计划程序对流程模式的影响的信息。例如，图 2 （A）和（B）显示了如果一条或另一条近端动脉被剪切时会发生的循环流。可以模拟血管削流或部署分流器等治疗方法，使医生和患者能够决定在每个特定情况下最起作用的。

Application and Summary

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

PC-MRI 的新方法能够显示更大的动态速度范围。然而，这受到患者扫描时间的严格限制。通常，患者数据以较低的分辨率采集，以减少在扫描仪中花费的时间。遗憾的是，这可能导致锯齿数据或信号丢弃，当速度编码梯度（VENC）设置过高时，问题会加剧。这可能错过缓慢和循环的流数据。将患者特定的流量和几何体与 CFD 配对，为捕获高分辨率血流动力学提供了一种有效的方法。

基于患者的建模具有固有的有用性，是它能够提供详细信息，而无需在通常具有非常不同特征的患者、疾病或治疗中进行概括。模拟允许医生和工程师在执行实际程序之前对替代治疗方案进行建模。模拟血流动力学可用于模拟流量分流支架、动脉旁路移植和导管造影注射等应用。当临床医生和患者希望获得最佳结果时，CFD 提供了一种观察术后流动的方法，从而提供更好的预见性。除了在引入设备或处理后描绘流量外，CFD 还允许估计墙壁的剪切应力。这与低 WSS 通常与动脉疾病进展相关的知识相配合，允许预测或概率建模。使用计算工具识别动脉瘤生长、血栓形成或出血的前体，为更早地识别高危患者开辟了可能性。总之，将患者特定的图像数据与CFD模拟相结合是疾病评估和手术预测的有力工具。

确认

作者要感谢西北大学的苏珊娜·施奈尔博士和迈克尔·马克尔博士向我们提供了数据中使用的4D患者数据。

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软件。

跳至...

0:07

Overview

1:20

Principles of Computational Fluid Dynamics

3:15

Generating Vessel Centerlines

4:24

Remapping 4D Flow MRI and Determining the Boundary Conditions

8:09

CFD Simulations

10:22

Results

11:17

Applications

12:17

Summary

此集合中的视频:

article

Now Playing

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

Biomedical Engineering

11.8K Views

article

使用光学和共聚焦显微镜成像生物样本

Biomedical Engineering

35.8K Views

article

生物样品的SEM成像

Biomedical Engineering

23.8K Views

article

纳米药物载体的生物分布:SEM的应用

Biomedical Engineering

9.3K Views

article

腹大塔的高频超声成像

Biomedical Engineering

14.5K Views

article

腹部主动脉瘤的定量应变图谱

Biomedical Engineering

4.6K Views

article

红外主塔中图像血液和脂质的光声断层扫描

Biomedical Engineering

5.7K Views

article

心脏磁共振成像

Biomedical Engineering

14.8K Views

article

腹主动脉瘤的近红外荧光成像

Biomedical Engineering

8.3K Views

article

非侵入性血压测量技术

Biomedical Engineering

11.9K Views

article

心电图（心电图）信号的采集与分析

Biomedical Engineering

105.9K Views

article

可吸收生物材料的拉伸强度

Biomedical Engineering

7.5K Views

article

小鼠脊髓的微CT成像

Biomedical Engineering

8.0K Views

article

大鼠非侵入性ACL损伤后膝关节退化的可视化

Biomedical Engineering

8.2K Views

article

结合 SPECT 和 CT 成像，使心脏功能可视化

Biomedical Engineering

11.0K Views

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