计算机模拟揭示主动脉下膜(SAM)病变的血流动力学扰动与瓣叶扑动机制

《Computers in Biology and Medicine》：In silico analysis of the haemodynamic disturbances caused by the subaortic membrane pathology

为探究主动脉瓣下狭窄(SAS)中瓣叶异常扑动（fluttering）现象这一尚未阐明的临床难题，Alessandra Monteleone等人利用无网格光滑粒子流体动力学(SPH)方法，开展了首次针对SAM病理的三维流固耦合(FSI)模拟研究。该研究比较了柔性薄层与刚性厚纤维肌性嵴两种常见SAM类型，揭示了其通过诱导涡旋周期性脱落触发瓣叶自持振荡的机制，并定量分析了膜刚度对振荡幅度和频率的显著调控作用，为理解SAM相关的瓣膜损伤和血流动力学并发症提供了关键的物理基础。

　　

心脏，这个不知疲倦的“泵”，其精妙的结构是生命得以延续的基石。主动脉瓣作为血液从左心室泵入全身动脉系统的最后一道“单向阀门”，其正常开合至关重要。然而，有一种被称为主动脉瓣下狭窄(SAS)的心脏疾病，会在主动脉瓣的入口下方形成一道异常的纤维组织“屏障”，即主动脉瓣下膜(SAM)。这道膜就像水管中的一个狭窄环，阻碍了血液的顺畅流出，不仅增加了心脏负荷，还常常伴有一个令人费解的现象——主动脉瓣瓣叶的快速、不规则抖动，即“瓣叶扑动”(leaflet fluttering)。

瓣叶扑动绝非小事。临床观察和体外实验均已证实，这种现象在人工心脏瓣膜中会加速材料疲劳、钙化、血液破坏(溶血)和瓣膜关闭不全(反流)，导致瓣膜过早失效。然而，在SAM这种天然瓣膜的病理条件下，扑动究竟如何产生？它的物理起源是什么？不同类型的SAM（如柔软的薄纤维膜与坚厚的纤维肌性嵴）对扑动有何不同影响？这些扑动又如何进一步影响局部的血流环境， potentially 损伤瓣膜组织本身？由于在活体内直接观测和量化这一动态过程极为困难，其背后的机制长期以来笼罩在迷雾之中，限制了我们对疾病进展的深入理解和治疗策略的优化。

为了拨开这层迷雾，揭开瓣叶扑动的神秘面纱，研究人员Alessandra Monteleone和Gaetano Burriesci开展了一项开创性的计算机模拟研究。他们运用前沿的计算力学手段，首次在三维空间中“重建”并深入分析了SAM存在下的心脏血流与瓣膜结构的复杂相互作用。这项研究成功揭示了瓣叶扑动源于SAM下游周期性涡旋的生成与脱落，并量化了膜刚度对这一现象的关键调控作用。相关成果发表于生物医学与计算机交叉领域的著名期刊《Computers in Biology and Medicine》，为从流体物理学角度理解这一临床现象提供了全新的视角和坚实的理论基础。

研究人员为开展此项研究，主要采用了以下几个关键技术方法：首先是基于无网格光滑粒子流体动力学(SPH)的“单物理场流固耦合(mono-physics FSI)”方法。该方法用粒子统一表示血液和固体结构（瓣叶、SAM），无需预设流固交界面，能高效处理大变形和接触问题。其次是建立了理想化的三维几何模型，包括健康主动脉瓣、柔性SAM（模拟离散薄层，杨氏模量50 kPa）和刚性SAM（模拟厚纤维肌性嵴，杨氏模量500 kPa）三种场景，其中SAM将流入腔横截面积阻塞了50%。最后，研究设定了简化的边界条件，在射血期对心室侧和主动脉侧施加恒定的压差以产生稳定的平均血流速度，从而专注于分离和揭示扑动现象的基本物理机制，排除了心脏搏动周期等其他因素的干扰。

研究结果

1. 投影动态瓣膜面积(PDVA)量化扑动

研究人员采用投影动态瓣膜面积(PDVA)作为量化瓣叶开合振荡的关键参数。结果显示：

• •
  健康瓣膜：PDVA在经历初始两次平滑振荡后，迅速稳定在约283 mm²，不再波动。
• •
  病理瓣膜（存在SAM）：两种SAM模型均诱发了持续的自持周期性振荡。与健康组相比，柔性SAM的平均PDVA降低了约3%，而刚性SAM降低了约9%。更重要的是，刚性SAM使振荡幅度增加了一倍以上（从3.12 mm²增至6.77 mm²），振荡频率也增加了8%。这表明SAM不仅减少了瓣膜的有效开口面积，而且其刚度显著调制了扑动的剧烈程度。

2. 血流模式与涡旋动力学

通过分析速度场、涡量(vorticity)图和流线，研究揭示了SAM引起的血流剧烈改变：

• •
  健康瓣膜：流线平滑有序，涡旋近乎消失。
• •
  病理瓣膜：SAM的存在形成了一个聚焦的中央射流，导致峰值流速大幅增加（柔性SAM+41%，刚性SAM+74%）。这个射流在SAM下游和瓣叶附近促进了强烈的涡旋结构形成。涡量图显示，在病理情况下，中性涡量区被压缩，高涡量尾迹从SAM游离缘脱落，主导了下游流动。

3. 瓣叶扑动的物理机制：涡旋的生成、迁移与脱落

通过分析一个完整PDVA振荡周期内四个特征时刻（波谷、上升过中线、波峰、下降过中线）的流线，研究清晰描绘了扑动的驱动机制：

• •
  过程：SAM的游离缘会周期性地产生并释放涡旋（研究中标记为A⁺, B⁺等）。这些涡旋在迁移过程中，会进入开放瓣叶的流入侧凹面。
• •
  驱动：位于瓣叶凹面的涡旋（如B⁺）对瓣叶施加一个向外的（扩张的）力，促使瓣叶进一步打开，PDVA增加。
• •
  释放：当涡旋增长到一定程度，会从瓣叶游离缘逃逸，向下游移动。此时瓣叶负荷减轻，PDVA开始下降。
• •
  循环：涡旋的周期性生成、迁移、对瓣叶施力及最终脱落，构成了一个自持的反馈循环，直接导致了瓣叶的周期性扑动。刚性SAM产生的涡旋更大、更持久，因此其驱动的扑动幅度和频率也更高。

4. 瓣叶表面的动力学与壁面剪应力(WSS)变化

• •
  瓣叶动力学：对瓣叶表面法向速度的分析显示，在扑动周期中，交替的膨胀（正速度）和收缩（负速度）“波前”有规律地从瓣叶基部向游离缘传播，直观反映了瓣叶的异常运动模式。
• •
  壁面剪应力(WSS)：WSS是作用在血管或瓣膜内壁上的摩擦力，其大小和方向与内皮细胞功能、组织重塑密切相关。研究发现：
  • •
    健康瓣膜：WSS向量方向一致，主要以轴向为主，模式规整。
  • •
    病理瓣膜：WSS分布变得高度不规则，方向杂乱，空间异质性显著。尤其在瓣叶的主动脉侧面，WSS幅值在所有分析时刻均显著低于健康情况。这种异常的WSS环境可能破坏瓣膜结构的完整性，并影响其钙化潜力。

研究结论与意义

本研究通过创新的计算模拟，首次在三维尺度上阐明了SAM导致主动脉瓣瓣叶扑动的流体动力学机制。核心结论指出：瓣叶的振荡性运动是SAM边缘涡旋形成和脱落的结果；SAM的刚度显著调节了扑动行为，刚性膜会导致更大振幅和更高频率的振荡。 尽管模拟未直接计算结构损伤和血液学并发症，但所识别的振荡与文献中报道的这些病理相关的血流动力学条件相符。

研究的重要意义在于多个层面：

1. 1.
   机制突破：它超越了此前临床观察中“高速射流冲击”的简单解释，从涡旋脱落和流固耦合的角度，为瓣叶扑动这一普遍而重要的现象提供了首个全面、定量的物理学解释。
2. 2.
   临床启示：研究明确区分了柔性膜和刚性膜对血流动力学影响的差异，这有助于解释为何不同类型的SAM患者临床表现和预后可能不同，为临床风险分层和手术时机选择提供了潜在的力学依据。
3. 3.
   技术验证：研究所采用的“单物理场FSI”方法成功应对了心脏瓣膜模拟中的大变形、复杂接触等挑战，证明了其在该领域的强大适用性，为后续更复杂的患者特异性模拟奠定了基础。
4. 4.
   病理链接：研究详细描绘了SAM存在下异常的WSS环境，这为理解SAM如何通过改变瓣叶表面的力学微环境，进而可能引发内皮功能障碍、促进瓣膜退行性变、钙化甚至血栓形成，提供了物理关联和思考方向。

总之，这项研究如同一台高精度的“数值显微镜”，让我们得以窥见SAM病理内部那些肉眼难以捕捉的瞬态流动与结构变形的细节。它不仅解答了一个长期存在的生理物理学问题，更架起了一座连接基础流体力学与临床心脏瓣膜病学之间的桥梁，为未来设计更优的诊疗方案和人工瓣膜提供了宝贵的理论工具和新颖视角。