从稳态到搏动流：分子通信中纳米颗粒在中等尺寸动脉传播的新范式

本文聚焦分子通信(MC)和磁性药物靶向(MDT)中纳米颗粒在血流中的传播问题，针对现有模型多忽略血流脉动特性和血管壁弹性的局限，研究人员开展了“搏动流与弹性血管壁下SPIONs在桡动脉中传播”的主题研究。通过对比经典层流与搏动流模型，发现搏动流下的信道冲激响应(CIR)呈现与心动周期同步的尖锐峰值序列，70%的颗粒运动发生在每个心动周期的前30%，显著影响码间干扰(ISI)和MDT的导向策略设计，对MC解调算法和MDT治疗优化具有重要指导意义。

在生物医学工程与纳米技术交叉的前沿领域，分子通信(Molecular Communication, MC)作为一种模仿生物系统信息传递的新范式，尤其在人体循环系统内实现靶向给药和信息传输方面展现出巨大潜力。其中，超顺磁性氧化铁纳米颗粒(Superparamagnetic Iron-Oxide Nanoparticles, SPIONs)因其良好的生物相容性和磁响应特性，被视为理想的药物载体和信息载体，在磁性药物靶向(Magnetic Drug Targeting, MDT)和体内通信系统中扮演着关键角色。然而，构建精确的颗粒传播模型面临巨大挑战：人体血液循环并非稳定不变的层流，而是由心脏搏动驱动的、具有强烈时变特性的搏动流(Pulsatile Flow)，同时血管壁也并非刚性，而是具有弹性的。现有的大多数模型为了简化分析，往往将血流假设为稳态的层流，并忽略血管壁的弹性变形，这可能导致模型预测与真实生理环境存在显著偏差，进而影响MDT的靶向精度和MC通信系统的可靠性。

为了弥补这一研究空白，来自德国埃尔朗根-纽伦堡大学智能电子与系统研究所的Luiz C. P. Wille等研究人员在《IEEE Transactions on Molecular, Biological, and Multi-Scale Communications》上发表了他们的最新研究成果。他们以临床中常用于药物注射的中等尺寸动脉——桡动脉(Radial Artery)为研究对象，创新性地将搏动流和血管壁弹性纳入考量，通过系统的仿真分析，对比了SPIONs在经典层流和更接近生理状态的搏动流下的传播特性，深入探究了信道冲激响应(Channel Impulse Response, CIR)的变化规律。

为了开展这项研究，研究人员主要采用了以下几种关键技术方法：首先，他们建立了代表左桡动脉的直管几何模型，参数基于解剖学数据（半径R = 1.575 mm，长度L = 23.5 cm）。其次，他们构建并比较了三种流体动力学仿真模型：（1）稳态层流（泊肃叶流）一维模型，作为经典对比基线；（2）基于广义达西弹性模型(Generalized Darcy's Elastic Model, GDEM)的一维搏动流模型，该模型能同时考虑血流脉动和血管壁弹性（弹性模量E = 1413 kPa）；（3）使用COMSOL®软件求解纳维-斯托克斯方程的三维搏动流刚性壁模型，用于验证一维简化模型的准确性。最后，在所有流场仿真的基础上，进行了大规模的粒子轨迹模拟（10⁵个颗粒），通过计算颗粒通过固定传感器位置（x_sensor= 0.8L）的速率和累积数量，得到了关键的CIR和累积粒子到达曲线。

II. 流模型和仿真设置

研究人员设定了统一的仿真环境，包括血管几何、血液物性参数（粘度η = 4×10^-3Pa·s，密度ρ = 1060 kg/m³）、颗粒属性（直径d_p= 200 nm）以及边界条件。通过计算雷诺数(Re = 75)和佩克莱特数(Pe = 2.5×10⁸)，确认在该桡动脉场景下可忽略湍流和扩散的影响，颗粒运动主要由流体的对流主导。入射流量Q_in(t)采用真实的心跳周期数据（T = 0.8 s）。

III. 粒子仿真结果

仿真结果清晰地揭示了搏动流与稳态流之间的根本差异。在层流模型（红色曲线）中，CIR呈现为单个的长尾峰，颗粒到达传感器的时间分布较广。而所有搏动流模型（黑色-1D刚性壁，绿色-1D弹性壁，蓝色-3D刚性壁）的CIR则表现出与心动周期同步的一系列尖锐峰值，其高度约为层流峰的5倍，主峰间隔为0.8 s（一个心动周期），并伴有间隔0.13 s的次峰，对应心脏的收缩期和舒张期。甚至出现了由血液短暂反向流动导致的负峰。值得注意的是，70%的颗粒在每個心动周期的前30%时间内通过传感器，表明颗粒运动高度集中在血流速度较高的阶段。

在累积粒子到达方面，尽管瞬时CIR差异巨大，但搏动流和层流模型在每个心动周期结束时（t = T, 2T, ...）的净颗粒传输量是相同的。这是因为搏动流中由高频分量引起的振荡速度场在一个周期内的积分为零。对比不同模型，1D和3D刚性壁模型的结果高度一致，表明在直血管中简单的1D模型足以准确描述颗粒传播。弹性壁模型则引起了CIR峰值的幅度和相位微小变化，但其净传输效应与刚性壁模型在周期末仍保持一致。

IV. 稳态和搏动累积CIR的解析比较

研究人员从理论上证明了累积CIR周期性匹配的必然性。他们通过傅里叶分析将搏动流速度场分解为稳态分量（对应层流）和振荡分量。振荡分量在一个周期内的积分贡献为零，因此净位移只由稳态分量决定，这解释了为何尽管瞬时行为不同，但周期平均传输效果与层流模型一致。

V. 结论

该研究通过精细的仿真分析，明确指出了在分子通信和磁性药物靶向应用中考虑血流搏动性和血管弹性的重要性。主要结论包括：1）搏动流下的信道冲激响应(CIR)与稳态层流有本质不同，表现为周期性尖锐峰值，这将对码间干扰(Inter-Symbol Interference, ISI)和解调算法设计产生重大影响，需要优化符号速率与心率的关系以及注射时机。2）对于直血管，计算效率高的一维模型（如GDEM）可以替代复杂的三维模型，为大规模网络仿真提供了可行性。3）血管壁弹性会引入CIR的幅值和相位调制，在追求更高精度的模型时应予以考虑。

这项研究为未来在更复杂血管网络（如分叉处）中研究颗粒传播奠定了基础，并指出了将传感器模型扩展为有限尺寸、以及进行实验验证是下一步的重要研究方向。该工作显著提升了对生理真实环境下纳米颗粒传播动力学的理解，对优化基于纳米颗粒的生物医学应用具有关键的指导意义。