血管内的血流动力学分析

血流动力学是研究血液在血管系统中运动规律及其与血管壁相互作用的学科。随着医学影像技术和计算流体力学（CFD）的快速发展，血管内血流动力学分析已成为评估动脉粥样硬化、动脉瘤破裂风险以及介入手术效果的重要手段。本文系统阐述血管内血流动力学的基本理论、数值建模流程以及典型案例，帮助读者快速上手并在实际工程或临床研究中开展可靠的仿真工作。

基本理论

血液的流体特性
- 在大多数动脉中，血液可近似视为非牛顿流体，但在剪切率较高的区域（>100 s⁻¹）常采用牛顿流体假设，黏度取 3.5 mPa·s 左右。
- 血液密度约为 1060 kg/m³。
血管壁的力学行为
- 刚性壁模型适用于短期稳态分析或血管壁相对刚硬的部位（如主动脉弓）。
- 弹性壁（FSI）模型能够捕捉脉动血流对血管扩张的反馈，常用于动脉瘤、瓣膜等高变形区域。
流动特征
- 层流 vs. 湍流：雷诺数 (Re = \frac{\rho U D}{\mu})。在大动脉 (Re) 常在 500–2000 之间，属于过渡流，需要采用 k‑ω SST 或 LES 湍流模型。
- 脉动性：血流呈周期性脉冲，常用 Womersley 参数 (\alpha = D\sqrt{\frac{\omega\rho}{\mu}}) 描述频率效应。

关键方程

连续方程（质量守恒）
[
\nabla \cdot \mathbf{u} = 0
]
动量方程（Navier‑Stokes）
[
\rho\left(\frac{\partial \mathbf{u}}{\partial t} + \mathbf{u}\cdot\nabla\mathbf{u}\right)
= -\nabla p + \mu \nabla^{2}\mathbf{u} + \mathbf{f}
]
壁面本构关系（弹性壁）
- 常用 线性弹性模型：(\sigma = E \varepsilon)
- 对于大变形可采用 超弹性模型（如 Mooney‑Rivlin）。

在实际 CFD 软件中，这些方程通过离散化（有限体积、有限元或谱元）求解，求解器会自动处理耦合项。

建模方法

方法	适用范围	优点	缺点
三维稳态 CFD	大血管、局部几何复杂	计算成本低，易收敛	不能捕捉脉动效应
三维瞬态 CFD	脉动血流、血管壁耦合	高保真度	计算资源需求大
一维血管网络模型	全身血流分布、参数敏感性分析	速度快，可与全身循环耦合	空间分辨率低
流固耦合（FSI）	动脉瘤、瓣膜、血管支架	同时考虑血流与壁变形	需要专门的耦合平台，收敛困难

边界条件与网格划分要点

入口流速：

inlet
{
    type            timeVaryingMappedFixedValue;
    fileName        "inletVelocity.dat";
    value           uniform (0 0 0);
}

出口阻抗（RCR）示例（OpenFOAM 自定义函数）
- R：近端阻抗，C：血管顺应性，Rd：远端阻抗。
- 通过 codedFixedValue 实现时间卷积。
网格质量：
- 非结构化六面体优于纯四面体。
- 最大非正交角 < 15°，**最小单元体积比** > 0.1，确保数值稳定。

数值求解与收敛

时间步长：满足 CFL 条件 (\text{CFL}= \frac{U\Delta t}{\Delta x} \le 0.5)。
残差收敛：对瞬态仿真，每个时间步残差 < 1e‑5，且 周期性监测（如血流量、压力）达到稳态或统计收敛。
加速技术：使用 PISO 或 PIMPLE 组合，提高瞬态耦合效率。

结果分析与临床意义

壁面剪切应力（WSS）
- 低 WSS（<0.4 Pa）与动脉粥样硬化斑块形成高度相关。
- 高 WSS 峰值可提示血管壁受力集中，易于动脉瘤破裂。
振荡剪切指数（OSI）
[
\text{OSI} = \frac{1}{2}\left(1-\frac{|\int_{0}^{T}\mathbf{\tau}w dt|}{\int{0}^{T}|\mathbf{\tau}_w| dt}\right)
]
- OSI 越大，血流方向波动越显著，提示血管壁受扰动。
压力梯度
- 在狭窄段，压力下降可用于 血流动力学阻力 评估，帮助决定是否进行支架植入。

示例案例：颈动脉内膜中层瘤（CIMT）血流分析

几何获取：从患者的颈部 CTA 获得左颈总动脉（CCA）至颈内动脉（ICA）分叉的 3D 重建，分辨率 0.5 mm。
网格：体积单元约 2.5 M，壁面层 10 层，最小 (y^{+}=0.5)。
入口条件：采用 1 Hz 心率的 Womersley 入口速度曲线，峰值 0.6 m/s。
出口阻抗：RCR 参数 R=1.2e⁹ Pa·s/m³，C=1.5e⁻⁹ m³/Pa，Rd=0.8e⁹ Pa·s/m³。
求解：使用 pimpleFoam，时间步长 0.001 s，模拟 3 个心动周期。
关键结果：
- 峰值 WSS 达 2.8 Pa，局部低 WSS 区域（<0.3 Pa）位于斑块下游 5 mm 处。
- OSI 最大值 0.42，提示血流分离区可能导致斑块进展。
- 通过比较支架植入前后 WSS 分布，发现支架显著提升了低 WSS 区域的剪切应力，理论上有助于抑制斑块增长。

该案例展示了血流动力学分析在 个体化治疗决策 中的价值。

小结

血管内血流动力学分析融合了流体力学、血管壁力学以及医学影像技术，是评估血管疾病、优化介入手术的重要工具。掌握 基本方程、合理的建模假设、精细的网格划分以及稳健的数值求解策略，即可在科研或临床环境中实现高质量的血流仿真。后续工作可进一步结合 机器学习 对仿真结果进行快速预测，或将 多尺度模型（1D‑3D‑FSI）用于全身循环的系统性分析。

血管内的血流动力学分析

血管内的血流动力学分析

基本理论

关键方程

建模方法

推荐工作流（以 OpenFOAM 为例）

边界条件与网格划分要点

数值求解与收敛

结果分析与临床意义

示例案例：颈动脉内膜中层瘤（CIMT）血流分析

小结