血流动力学仿真
血流动力学的理论基础
血流动力学FSI概述
血管内血流仿真中什么时候需要FSI？
动脉瘤破裂风险评估、支架置入后再狭窄预测、冠状动脉旁路移植的设计等。血管壁是弹性体，随脉动径向变形5～10%。这种壁面变形会影响血流模式，因此需要FSI。
控制方程
血管壁的力学特性如何建模？
血管壁作为非线性超弹性体处理。广泛采用Holzapfel-Gasser-Ogden模型。
其中$\mu$是基质刚度，$k_1, k_2$是胶原蛋白纤维的刚度参数，$I_4^{(i)}$是纤维方向的伪不变量。
流体侧是非压缩性Navier-Stokes方程吗？
是的。在ALE框架中描述。
大血管内通常将血液近似为牛顿流体（$\mu \approx 3.5$ mPa·s），但在低剪切速率区域需要采用Carreau-Yasuda模型等非牛顿模型。
血液是"非牛顿流体"——粘度随流速变化的神奇液体
工程学中最初学到的流体是水和空气这样的"牛顿流体"，但血液完全不同。血液的粘度会随剪切速率（流动速度的梯度）变化，这是一种非牛顿流体。在低流速下粘度很高，高流速下粘度会降低。这是因为红细胞在低流速时会相互重叠（形成罗氏）堆积，在高流速下会分散和排列，从而阻力下降。在毛细血管中（直径～8μm），红细胞（直径～7μm）会一个接一个地通过，血液的行为几乎更像是粒子行进而不是粘性流体。用CFD计算血流时，应根据分析的血管部位和流速范围决定是采用Carreau-Yasuda模型等非牛顿粘度模型，还是简单的牛顿近似（μ≒3.5mPa·s）。
血流动力学的数值计算方法
稳定化有限元法
听说血流CFD需要SUPG/PSPG稳定化，为什么呢？
等次插值（P1-P1）不满足inf-sup条件，压力会出现虚假振荡。加入SUPG/PSPG稳定化后可得到稳定解。
在血管分岔部和支架周围，局部雷诺数很高，稳定化参数的设计直接影响计算稳定性。
耦合算法
血管的FSI可以用弱耦合吗？
血管壁很薄，附加质量效应明显，弱耦合容易发散。推荐采用Robin-Neumann分割或单一结构法。基于广义Robin条件的半隐式耦合法引起关注。
基于图像的建模工作流程
从临床图像到网格生成的流程是怎样的？
1. 获取CT血管成像或MRA影像
2. 分割（Mimics, ITK-SNAP, 3D Slicer）
3. 表面平滑化和缺陷修补
4. 体积网格生成（TetGen, VMTK, ICEM CFD）
5. 插入边界层网格（VMTK的boundary layer功能）
6. 确定入口、出口面并添加延伸管
VMTK（血管建模工具箱）是血管专用开源软件，可一次性处理从中心线提取到边界层网格的所有步骤。
格子Boltzmann法（LBM）为何适合血流分析
近年来格子Boltzmann法（LBM）在血流分析中备受关注。LBM采用完全不同于传统CFD（Navier-Stokes有限体积法）的方法，将流体视为"虚拟粒子的统计分布"。为什么LBM适合血流？因为它善于计算红细胞等复杂形状粒子混杂的悬浊液流体，网格生成的难度比传统方法大幅降低。在毛细血管水平的计算中，LBM可以逐个追踪红细胞，进行粒子分辨计算。东北大学的团队在日本超级计算机"富岳"上首次实现了毛细血管网络（直径5～10μm）的LBM计算，建立了超过5亿单元的模型，重现了红细胞的流动。
血流动力学的实际应用
边界条件设置
入口和出口的边界条件如何设置？
入口处设置从4D Flow MRI或多普勒超声获得的流量波形。通常以Womersley解作为初始速度分布。
Womersley数 $\alpha = R\sqrt{\omega/\nu}$ 越大（大动脉约$\alpha \approx 15$），速度分布越平坦。
出口的Windkessel模型如何设置？
采用3单元Windkessel，设置近端阻力 $R_p$、远端阻力 $R_d$ 和顺应性 $C$。有多个出口时，根据Murray规律按直径的立方比分配流量。
壁面剪应力评估指标
结果中应该看哪些指标？
时间平均壁面剪应力（TAWSS）和振荡剪应力指标（OSI）是代表性指标。
低TAWSS（< 0.4 Pa）且高OSI（> 0.3）的区域被认为动脉硬化风险高。FDA指南也推荐这些评估指标。
用CFD识别脑动脉瘤的"危险区域"——临床应用前沿
利用血流仿真评估脑动脉瘤（蛛网膜下腔出血的原因）破裂风险的研究正在进行。通过患者特异性的MRI/CT血管造影建立3D血管模型，用CFD计算壁面剪应力（WSS）。在WSS极低的区域（以0.4Pa为目标），内皮细胞容易变性，瘤壁变薄的风险会升高。东京大学和庆应义塾大学医院的联合研究表明，通过CFD计算"低WSS区域的面积比"，可以以约80%的准确度预判破裂动脉瘤与未破裂动脉瘤。"是否应该手术还是观察随访"这个困难的判断，正在被血流分析提供定量依据。
血流动力学软件对比
血流动力学分析工具对比
血流动力学FSI分析有哪些可用工具？
从专用工具到通用求解器有很多选择。
| 工具 | 类别 | FSI支持 | 特点 |
|---|---|---|---|
| SimVascular | 开源（斯坦福大学） | CMM法 | 血管专用。图像→网格→FSI一体化工作流程 |
| CRIMSON | 开源（密歇根大学） | 有 | SimVascular衍生版。GUI功能完善 |
| Ansys Fluent + Mechanical | 商用 | System Coupling | 通用性强。Windkessel需自定义实现 |
| STAR-CCM+ | 商用 | 强耦合FSI | co-simulation接口 |
| COMSOL Multiphysics | 商用 | 单一结构法 | 小规模模型适用。内置超弹性模型 |
| Alya | 开源（巴塞罗那超计算中心） | 有 | 针对HPC优化。支持数亿单元 |
FDA规制要求用哪个工具更合适？
FDA强调V&V（验证与验证），不指定特定工具。需要按照ASME V&V 40进行验证。在FDA基准nozzle模型上有验证实绩的Ansys Fluent和STAR-CCM+在规制申报中使用最多。
研究用途的话开源软件更好吗？
SimVascular和CRIMSON源代码公开，可验证性强，方法自由度也高。在Nature Medicine等顶级期刊的使用实绩丰富。商用工具的优点是软件QMS体系完善，规制申报文件的准备容易。
Simendo vs. Fluent——血流分析专用工具诞生的背景
通用CFD软件（Fluent、CFX、OpenFOAM）也可用于血流分析，但医疗器械制造商使用的专用工具如"Simendo（现已成为HeartFlow）"和"Mimics+3-matic（Materialise）"是为血流分析专门开发的工作流程。这些工具产生的背景是解决汎用CFD血流分析工作流程的瓶颈——从DICOM图像制作血管网格需要专业知识，边界条件设置没有明确指南。HeartFlow开发的"FFRCT（非侵入性冠状动脉血流储备分数计测）"基于血流CFD分析结果定量评估心肌梗死风险，2019年获得FDA批准。这是CFD直接成为"临床诊断依据"的历史性事件。
血流动力学前沿研究
0D-3D耦合模型
全身血管系统用3D求解是不可能的吧。怎么处理的？
只用3D-FSI求解感兴趣的区域，其余血管系统用阶跃参数模型（0D模型）表示。这称为几何多尺度耦合。SimVascular将这个0D-3D耦合作为标准配置。
机器学习加速
实时血流动力学预测可能吗？
用POD-Galerkin ROM曾经实现了血管FSI计算的3位数加速。图神经网络直接将患者特异形状转换为图结构来预测血流动力学参数的方法也在研究中。
血细胞溶血评估
机械瓣膜中的红细胞损伤如何预测？
传统的幂律溶血模型被张量基溶血模型取代。
还有用DEM逐个追踪红细胞的方法和格子Boltzmann法微观仿真也在研究。从宏观CFD通过拉格朗日粒子追踪获得代表性血球轨迹，然后在该处运行微观仿真的多尺度策略最为现实。
药物洗脱支架的多物理耦合
药物洗脱支架（DES）的设计也用FSI吗？
除了血流-结构耦合外，还需要解决药物运输和血管壁渗透的多物理耦合问题。
COMSOL Multiphysics在这种多物理耦合的单一环境设置上更容易。
血流的数字孪生——患者个体动脉硬化风险预测的未来
血流力学的最前沿是开发"患者血管形状+血压波形+血液粘度"作为输入、预测动脉硬化进展速度的数字孪生。通常动脉硬化进展需要10～20年，但CFD仿真数小时内就能显示"10年后这个部位的斑块不稳定化的风险"。欧洲HEART研究项目正在开发结合400多名患者数据与血流仿真的风险模型，目标是与AI结合解读CFD结果，比既有临床指标提高20～30%的心血管事件预测准确度。血流分析正在从"过去状态诊断"进化到"未来风险预测"。
血流动力学故障排除
非生理压力振荡
计算结果的压力波形出现非物理振荡。
整理一下典型原因和对策。
| 原因 | 症状 | 对策 |
|---|---|---|
| 出口边界条件不当 | 压力波反射 | 引入Windkessel模型，添加出口延伸管 |
| 初始条件不匹配 | 初期数步内出现大幅振荡 | 从定常解重启，流量渐增 |
| 时间步长过大 | 高频振荡 | 将 $\Delta t$ 设为0.1 ms以下 |
| 弱耦合不稳定 | 振幅随时间增长 | 切换到强耦合 |
Windkessel参数调整
参数设置有困难。有诀窍吗？
推荐值如下。
分割影响评估
图像分割精度对结果的影响有多大？
已有报告称内径误差±10%会导致WSS变化±40%。对策是多个操作员对比、执行UQ、用统计形状模型（SSM）系统评估形状变差。ASME V&V 40要求作为Credibility Evidence的一部分评估输入不确定性。
"血压波形异常"——入口边界条件设置上90%都栽跟头
血流CFD分析中最容易犯错的是入口边界条件设置。心脏搏动产生的脉动流的"速度分布随时间变化"，因此简单的"恒定流速"完全不合适。理想做法是从患者的多普勒超声测定的特异性入口流量波形进行Fourier展开后输入CFD，但如果不将入口设为"Womersley流动分布"（振荡流的解析解），入口附近就会产生非物理流动。常见症状是"入口起向下游5～10D（直径）出现异常压力振荡"，这几乎100%是因为采用了简单的"均匀流入"而没有正确实现Womersley边界条件。根本对策是将入口部分延长至20D以上，或正确实现Womersley边界条件。
细节
错误