🧬 biology

Microstructure-Induced Variability in Residence Time Distributions of Fibrous Porous Media

本研究利用计算流体力学证明，多孔介质中纤维的空间排列方式（而非仅由体积分数决定）是决定停留时间分布变异性的主要因素，其中有序的六角结构有助于促进塞流行为，而聚集构型则会诱发显著的弥散并延长平均停留时间。

原作者： E.G. Karvelas, A. Tsiantis, P. Chatzisaroglou, T.D. Papathanasiou

发布于 2026-09-07

📖 1 分钟阅读☕ 轻松阅读

原作者： E.G. Karvelas, A. Tsiantis, P. Chatzisaroglou, T.D. Papathanasiou

原始论文采用 CC BY 4.0 许可（https://creativecommons.org/licenses/by/4.0/）。 ⚕️ 这是一篇未经同行评审的预印本的AI生成解释。这不是医疗建议。请勿根据此内容做出健康决定。阅读完整免责声明

在生物体内，组织并非由细胞构成的实心块状，而是由嵌入在纤维网格中的细胞所组成的复杂、海绵状的网络。这种被称为细胞外基质的网格，既是提供结构的支架，也是物质运输的高速公路。必需的营养物质、氧气和化学信号必须通过这个网络才能到达每个细胞，同时代谢废物也必须通过它被带走。这些物质的移动是两种力量之间微妙的平衡：一种是扩散，即分子从高浓度区域向低浓度区域随机漂移；另一种是对流，即流体流动将它们推向前行。在健康的组织中，这种运输是高效且可预测的。然而，当网格的结构发生变化时——可能是由于疾病、损伤或医疗植入物的设计——流体的流动方式可能会变得混乱，导致生命维持物质的输送不均，或导致药物的危险堆积。

一组研究人员致力于理解纤维在这一网格中的排列方式是如何决定流体运输的速度和可预测性的。他们专注于一个被称为“停留时间分布”的概念，这本质上是衡量一滴流体从材料的一端移动到另一端所需的时间。如果材料是均匀的，每一滴流体所需的时间大致相同。如果材料是无序的，一些流体会通过快速通道疾驰而过，而另一些则会卡在缓慢的死角中，从而产生广泛且不可预测的旅行时间分布。研究人员想要知道，仅仅是纤维的密度是否足以解释这些差异，还是纤维排列的具体模式才是真正的关键所在。

为了进行调查，科学家们构建了详细的纤维材料计算机模型，模拟水流通过其中的过程。他们创建了三种截然不同的纤维排列类型进行对比。第一种是有序的六角形图案，其中纤维按照完美的重复网格排列，类似于蜂窝。第二种是随机排列，其中纤维在没有任何固定模式的情况下散布，模仿了许多生物组织中存在的自然无序性。第三种是集群排列，其中纤维聚集在一起形成致密的团块，在团块之间留下巨大的开放空间。随后，他们运行了模拟实验，观察在不同流动条件下，一种被动示踪剂（一种用于追踪运动的无害物质）如何在每种结构中移动。

结果显示，纤维的空间排列远比纤维的密度本身更为重要。在完美的有序六角形结构中，流体移动得非常均匀。示踪分子几乎在同一时间到达出口，形成了一个尖锐且可预测的脉冲。这种行为被称为“塞流”（plug flow），即流体作为一个凝聚的整体移动，且扩散极小。即使研究人员增加纤维的数量使网格变得更密，流动依然保持得异常一致且可预测。有序的网格迫使流体进入均匀的通道，防止了缓慢区域或快速河流的形成。

相比之下，随机排列展现了不同的情况。虽然流体仍在移动，但旅行时间更加分散。一些分子找到了宽阔、开放的路径并快速穿过，而另一些则在狭窄、蜿蜒的间隙中艰难前行，速度变慢。这产生了一个更宽的到达时间分布，表明纤维放置中的无序性为运输过程引入了一定程度的不确定性。流体并非作为一个整体移动，而是作为一个速度各异的流体进行。

在集群排列中，观察到的效应最为剧烈。在这里，流体运输变得高度不稳定。致密的纤维团块创造了流体几乎不移动的停滞区，而团块之间的开放空间则充当了高速通道。这导致了一种情况：虽然部分示踪分子到达得非常快，但很大一部分分子被困在缓慢区域，需要更长的时间才能排出。结果是旅行时间数据呈现出一个长而沉重的“尾部”，这意味着虽然平均时间看起来尚可，但对于任何单个分子而言，其实际体验是高度多变的。团块的存在极大地降低了运输的可预测性，导致流体的分散程度远高于其他两种情况。

研究人员发现，这种变异性是由纤维的拓扑结构而非仅仅是纤维的数量驱动的。即使在不同模型中使用的总纤维量相同，纤维的组合方式也会完全改变结果。有序网格维持了稳定的流动，随机散射引入了一些混合，而集群则创造了一个带有长时间延迟的混沌环境。这表明，在生物系统中（例如肿瘤或工程化组织支架），纤维集群的具体模式可能是影响营养到达细胞或药物分布程度的主要因素。具有集群纤维的组织可能会遭受药物暴露不均的问题，即由于携带药物的流体卡在了缓慢区域，导致某些细胞接收到高剂量，而另一些细胞则完全没有接收到。

这些基于理想化二维模型计算机模拟得出的发现，为如何设计用于医疗用途的纤维材料提供了新的视角。如果目标是确保药物或营养物质均匀地到达组织的每个部分，则必须精确控制纤维的排列，以避免集群现象。相反，如果目标是减缓某种物质的释放，那么创造特定的集群可能是一种有用的策略。虽然这项研究受限于简化的模型，且未考虑真实人体组织的复杂三维特性，但它确立了一个明确的原则：网格的架构与材料本身的密度同样重要，甚至更为重要。通过理解纤维模式如何影响流动，科学家和工程师可以更好地预测生物组织中的运输过程，并设计出更有效的医疗方案和植入物。

技术摘要：纤维状多孔介质中微观结构诱导的停留时间分布变异性

问题陈述
生物组织作为层次化结构发挥作用，其中细胞居住在纤维状细胞外基质（ECM）中，该基质负责提供机械支撑和溶质传输。尽管以往的研究已基于纤维密度和排列方式表征了纤维状多孔介质的稳态水力渗透率，但在理解微观结构组织如何影响溶质传输时间的瞬态统计分布方面仍存在空白。传统的停留时间分布（RTD）模型通常假设均匀的塞流（plug flow）和空间均匀的扩散，而这些假设在真实的纤维状介质中是被违背的，因为不规则的纤维排列会产生优先流通道（“河流”）和停滞区。本研究调查了在理想化的纤维几何结构中，纤维排列（有序 vs. 无序）、体积分数 ( $\phi$ ) 和 Peclet 数 ( $Pe$ ) 如何影响水力渗透率和 RTD 统计特性，旨在将微观结构组织的影响从仅由孔隙率决定的影响中分离出来。

方法论
研究采用计算流体力学（CFD）工具 OpenFOAM，在包含固定半径（ $10\ \mu m$ ）圆柱形纤维的二维单元格（ $L \times H$ ）中进行模拟。研究生成了三种不同的微观结构配置：

六方晶格（有序）： 纤维放置在固定的、相等的纤维间距处。
随机： 通过随机顺序添加（RSA）算法生成。
聚集： 使用分子动力学（HOOMD-blue）通过 NVT/Langevin 动力学在降低温度时诱导聚集而成。

模拟过程分为两个步骤：

稳态流动： 在恒定压降下求解不可压缩 Navier-Stokes 方程，以获得稳态速度场，使用约 $4 \times 10^6$ 个单元格的网格，并在纤维周围进行局部加密。
瞬态传输： 通过在收敛的速度场上求解非稳态对流-扩散方程来模拟被动示踪剂。

关键参数包括纤维体积分数 ( $\phi = 0.1$ 至 $0.35$) 和 Peclet 数 ( $Pe = 0.01$ 至 $5$)。研究利用基于矩（moment-based）的 RTD 框架来分析突破曲线，计算了平均无量纲传输时间 ( $\langle t^* \rangle$ ) 和有效 Peclet 数 ( $Pe_{eff} = 2/\sigma^2_\theta$ )，其中 $\sigma^2_\theta$ 是归一化 RTD 的方差。模型验证是针对 Gebart 的横向渗透率解析模型进行的。

主要结果
研究揭示了微观结构组织对传输变异性的层次化影响：

有序（六方）排列： 这些配置产生分布较窄的 RTD，其峰值在 $\theta \approx 1$ 附近表现得非常尖锐，表明接近塞流行为。即使随着 $\phi$ 的增加， $Pe_{eff}$ 仍保持在高水平，且传输时间的方差极小。系统表现出准塞流行为，纵向色散有限。
随机排列： 这些配置产生了较宽的 RTD，具有扁平化的峰值，并且与有序结构相比向更高的传输时间偏移。在匹配条件下， $Pe_{eff}$ 比六方晶格低一个数量级，反映了由于速度场异质性和优先流通道导致的色散增加。
聚集排列： 这些微观结构表现出最显著的传输变异性。RTD 呈现强烈的右偏分布和长尾特征，表明在致密纤维聚集体内部出现了缓慢路径和局部滞留区。这导致了升高的平均停留时间和 $Pe_{eff}$ 的急剧下降。与有序或随机阵列不同，聚集介质中的传输强烈依赖于 $\phi$ ，其突破时间随密度显著变化。

意义与主张
本文声称，微观结构组织而非仅靠体积分数，是纤维状多孔介质中传输变异性的主要决定因素。研究结果表明，纤维排列中的空间相关性（聚集）会诱导宏观色散，而这种色散无法仅通过孔隙率指标来预测。

作者将这些结果定位为关于理想化二维纤维介质中微观结构-传输耦合的基础见解。他们认为，这些发现对于纤维材料的设计具有相关性，例如工程支架和药物递送载体，在这些应用中传输的可预测性至关重要。具体而言，研究暗示通过控制纤维排列和密度，可以实现对传输可预测性的调控，从而潜在地影响营养物质的均匀化、氧气渗透以及药物停留时间。

局限性与范围
作者明确指出，这是一个基于理想化二维几何结构的假设生成型研究。他们并未声称定量地模拟特定的生物组织。承认的主要局限性包括：

使用二维模拟，未能捕捉到诸如纤维波纹度、面外连通性或轴向曲折度等三维效应。
假设分子扩散系数在空间上是均匀的，这隔离了速度场异质性的影响，但也排除了真实组织中存在的空间变化的表观扩散率。
缺乏将 RTD 统计数据与纤维状组织传输数据进行直接实验验证的过程。
将 $Pe_{eff}$ 解释为一种总结性的宏观色散指数，而非针对聚集介质中高度偏态 RTD 的字面高斯等效 Peclet 数。

因此，虽然该研究提供了一个将微观结构与停留时间统计联系起来的定量框架，但作者警告说，将这些结果推演至三维病理组织（如肿瘤 ECM）需要经过具有生理代表性几何形状和实验数据的进一步验证。

技术摘要：纤维状多孔介质中微观结构诱导的停留时间分布变异性

类似论文