血流动力学仿真
理论与物理
血流动力学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）来处理，取决于所分析的血管部位和流速范围。
各项的物理意义
- 结构-热耦合项：温度变化引起的热膨胀诱发结构变形，变形又影响温度场。$\sigma = D(\varepsilon - \alpha \Delta T)$。【日常示例】夏天铁轨伸长导致缝隙变窄——温度上升→热膨胀→产生应力的典型例子。电子电路板在焊接后翘曲，也是不同材料热膨胀系数差异所致。发动机气缸体因高温部和低温部的温差产生热应力，最坏情况下会导致裂纹。
- 流体-结构耦合（FSI）项：流体压力·剪切力使结构变形，结构变形又改变流体区域，是双向的相互作用。【日常示例】强风使悬索桥的缆索振动（涡激振动）——风力使结构摇晃，摇晃的结构改变风流，进而放大振动。心脏血流与血管壁的弹性变形、飞机机翼的颤振（气动弹性不稳定性）也是典型的FSI问题。有时仅单向耦合即可，但变形较大时双向耦合是必须的。
- 电磁-热耦合项：焦耳发热 $Q = J^2/\sigma$ 引起温度上升，温度变化又改变电阻，形成反馈回路。【日常示例】电暖炉的镍铬丝通电后发热（焦耳热）变红——温度升高则电阻改变，电流分布也随之变化。IH电磁炉的涡流发热、输电线路因温度升高导致的垂度增加也是此耦合的例子。
- 数据传递项：通过插值解决不同物理场间网格不匹配的问题。【日常示例】天气预报中，将“气温数据”和“风速数据”结合计算体感温度时，若观测地点不同就需要插值——CAE的耦合分析中，结构网格与CFD网格通常也不一致，因此界面处的数据传递（插值）精度直接关系到结果的可靠性。
假设条件与适用范围
- 弱耦合假设（单向耦合）：当一方物理场影响另一方而反向影响可忽略时有效
- 需要强耦合的情况：FSI中的大变形、电磁-热耦合中温度依赖性较强时
- 时间尺度的分离：各物理场的特征时间差异较大时，可通过子循环提高效率
- 界面条件的协调性：需确认耦合界面处的能量·动量守恒在数值上得到满足
- 不适用的案例：三个以上物理场同时强耦合时，有时需要整体式方法
量纲分析与单位制
| 变量 | SI单位 | 注意事项·换算备忘 |
|---|---|---|
| 热膨胀系数 $\alpha$ | 1/K | 钢：约12×10⁻⁶、铝：约23×10⁻⁶ |
| 耦合界面力 | N/m²（压力）或N（集中力） | 确认流体侧与结构侧的力平衡 |
| 数据传递误差 | 无量纲（%） | 插值精度取决于网格密度比。5%以下为目标 |
数值解法与实现
稳定化有限元法
听说血流CFD需要SUPG/PSPG稳定化，为什么？
等阶插值（P1-P1）不满足inf-sup条件，会导致压力出现伪振荡。添加SUPG/PSPG稳定化后可获得稳定解。
在血管分叉部或支架周边，局部Re数较高，稳定化参数的设计直接关系到计算稳定性。
耦合算法
血管的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（Vascular Modeling Toolkit）是血管专用开源软件，可从中心线提取到边界层网格进行一体化处理。
格子玻尔兹曼法（LBM）适合血流分析的意外理由
近年来，格子玻尔兹曼法（LBM）作为一种血流分析的数值方法备受关注。它采用与常规CFD（Navier-Stokes方程的有限体积解法）完全不同的方法，将流体视为“虚拟粒子的统计分布”。之所以适合血流分析，是因为它擅长计算像红细胞这样复杂形状粒子混杂的悬浮流体，并且网格生成比传统方法简单得多。在毛细血管级别的计算中，可以逐个追踪红细胞进行粒子分解计算，东北大学的研究小组在超级计算机“富岳”上首次实现了毛细血管网络（直径5~10μm）的LBM计算，在超过5亿单元的模型中再现了红细胞的流动。
整体式方法
将所有物理场作为一个联立方程组同时求解。对强耦合问题稳定，但实现复杂，内存消耗大。
分区法（分离迭代法）
各物理场独立求解，在界面交换数据。实现容易，可利用现有求解器。适用于弱耦合。
界面数据传递
最近邻法（最简单但精度低）、投影法（具有守恒性）、RBF插值（对网格不一致鲁棒性强）。守恒性与精度的平衡很重要。
子迭代
在每个耦合步内进行充分迭代，确保界面条件的协调性。残差基准基于各物理场的典型值进行缩放。
Aitken松弛
自动调整耦合迭代的松弛系数。防止因过度松弛导致发散，是加速收敛的自适应方法。
稳定性条件
注意附加质量效应（流体-结构耦合中结构密度≈流体密度时）。不稳定时可应用Robin型界面条件或IQN-ILS法。
Aitken松弛的比喻
Aitken松弛类似于“平衡跷跷板”。一方推得太用力，另一方就会弹起，反弹又导致推得更用力——为了抑制这种振荡，自动调整推力大小的就是Aitken松弛。当耦合迭代振荡不收敛时，根据上次的修正量自动调整下次修正量的自适应方法。
实践指南
边界条件设置
入口和出口的边界条件如何设置？
入口设置从4D Flow MRI或多普勒超声获得的流量波形。常用Womersley解作为初始速度剖面。
Womersley数 $\alpha = R\sqrt{\omega/\nu}$ 越大（主动脉中 $\alpha \approx 15$），剖面越平坦。
出口的Windkessel模型如何设置？
使用三元件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分析。
初学者易犯的错误
“单向耦合足够了吧？”——这个判断失误在耦合分析中最危险。如果结构变形微小，单向确实足够。但像心脏瓣膜开闭这样变形会大幅改变流路的情况，单向则完全不行。判断标准是“变形量是否超过特征长度的1%”。超过的话，双向耦合是必须的。如果错误地只用单向耦合，结果会“看似合理实则大错特错”——这是最可怕的模式。
边界条件的思考方式
耦合界面的数据交换如同“国境的出入境管理”。各国（物理场）有各自的法律（控制方程），但若不在国境（界面）准确管理人员与物资（力、温度、位移）的往来，两国的经济（能量平衡）就会崩溃。网格不一致时的插值如同“翻译”——误译（插值误差）越小，结果越好。
软件比较
血流动力学分析工具比较
有哪些可用于血流动力学FSI分析的工具？
从专用工具到通用求解器，范围很广。
| 工具 | 类型 | FSI支持 | 特点 |
|---|---|---|---|
| SimVascular | OSS（斯坦福大学） | CMM法 | 血管专用。从图像→网格→FSI的一体化流程 |
| CRIMSON | OSS（密歇根大学） | 支持 | SimVascular衍生。GUI丰富 |
| ANSYS Fluent |
なった
詳しく
報告