[论文解读] Multiscale modeling of vascularized tissues via non-matching immersed methods

Luca Heltai, Alfonso Caiazzo|arXiv (Cornell University)|Aug 29, 2019

Advanced Mathematical Modeling in Engineering参考文献 31被引用 14

一句话总结

本文提出了一种非匹配浸入有限元方法，用于高效模拟血管化组织的多尺度特性，其中受压的1D血管被表示为共维数为二的流形（如血管中心线）上的奇异或超奇异源项，从而显著降低计算成本。该方法通过统计模拟实现了对双相组织有效力学性能的体外表征，结果表明血管压力会诱导可测量的宏观剪切效应和各向异性，且与血管取向相关。

ABSTRACT

We consider a multiscale approach based on immersed methods for the efficient computational modeling of tissues composed of an elastic matrix (in two or three-dimensions) and a thin vascular structure (treated as a co-dimension two manifold) at a given pressure. We derive different variational formulations of the coupled problem, in which the effect of the vasculature can be surrogated in the elasticity equations via singular or hyper-singular forcing terms. These terms only depend on information defined on co-dimension two manifolds (such as vessel center line, cross sectional area, and mean pressure over cross section), thus drastically reducing the complexity of the computational model. We perform several numerical tests, ranging from simple cases with known exact solutions to the modeling of materials with random distributions of vessels. In the latter case, we use our immersed method to perform an in silico characterization of the mechanical properties of the effective biphasic material tissue via statistical simulations.

研究动机与目标

• 开发一种高效的计算模型，用于模拟血管化组织，能够在不于网格层面解析血管结构的情况下捕捉受压微血管的力学影响。
• 解决磁共振弹性成像（MRE）中的多尺度组织力学挑战，其中临床分辨率限制了对完整微血管结构的解析。
• 实现对具有随机血管分布的双相组织中有效力学性能（如剪切模量和压缩模量）的统计表征。
• 提供一个将微尺度血管几何结构与压力关联至宏观组织行为的框架，以改进体内组织表征。
• 通过在固定网格上快速模拟多种血管构型，为弹性成像中的未来反问题提供支持。

提出的方法

• 采用变分浸入有限元公式，将1D血管作为奇异或超奇异源项嵌入3D（或2D）弹性基质问题中。
• 通过共维数为二的流形（如中心线）表示血管，源项仅依赖于血管中心线、横截面积和平均压力。
• 采用非匹配有限元方法，使同一背景网格可在不同血管构型下重复使用，无需重新划分网格。
• 应用超奇异公式，将源项力集中于血管中心线，从而在保持精度的同时最小化计算成本。
• 推导均匀血管分布的均质化模型，以解析地关联血管压力与体积分数至有效组织刚度。
• 通过随机生成血管结构的统计模拟，估算在不同几何与压力条件下组织的有效力学性能。

实验结果

研究问题

• RQ1 如何在不解析血管几何结构的前提下，准确表示受压1D血管对3D弹性基质的力学影响？
• RQ2 血管取向与压力对双相组织有效剪切模量和压缩模量有何影响？
• RQ3 非匹配浸入方法是否能高效地在固定背景网格上模拟多种血管构型，以实现统计表征？
• RQ4 在复杂血管网络中，超奇异公式与奇异公式在精度与计算成本方面有何比较？
• RQ5 血管压力在多大程度上诱导组织的宏观剪切效应？该效应与血管网络拓扑结构有何关联？

主要发现

• 该浸入方法即使在复杂、弯曲且相交的血管网络中，也能以高计算效率成功捕捉受压血管对周围弹性基质的力学影响。
• 统计模拟表明，组织的有效剪切模量与下层血管的取向和压力密切相关，表明存在显著的各向异性。
• 超奇异公式将力局部化于血管中心线，在保持与奇异公式相当精度的同时，显著降低了计算成本。
• 该方法通过在单一固定网格上模拟多种随机血管分布构型，实现了对血管化组织力学参数的体外估算。
• 结果表明，血管压力可诱导可测量的宏观剪切效应，且剪切量的大小与血管的空间排布高度相关。
• 该框架通过允许在不重新组装线性系统的情况下快速调整血管压力与几何结构，为弹性成像中的未来反问题提供了支持。