扩散势存在下圆柱孔内带电介电液滴的扩散泳动

《Colloids and Interfaces》：Diffusiophoresis of a Charged Dielectric Fluid Droplet in a Cylindrical Pore in the Presence of Diffusion Potential

　　

研究人员对诱导扩散势（diffusion potential）存在下（例如NaCl电解质溶液中）圆柱孔内介电液滴（dielectric droplet）的扩散泳动（diffusiophoretic motion）进行了理论分析。采用基于切比雪夫多项式（Chebyshev polynomials）的拼接伪谱法（patched pseudo-spectral method）结合几何映射方案处理不规则求解域，求解了基本电动力学控制方程。详细考察了边界约束效应（boundary confinement effect）对液滴迁移率（mobility）的影响。研究中发现了有趣的电动力学现象，例如窄圆柱孔内的迁移率反转（mobility reversal），由于强烈的边界约束效应，液滴运动方向与基于经典库仑静电力定律的预期方向相反。此外，在特定孔半径处还发现了“固化现象”（solidification phenomenon），此时无论液滴粘度如何，液滴均像刚性粒子一样运动，内部无循环涡流。找到了对应的临界点λ（液滴半径与圆柱半径之比），每当通过临界点时，液滴表面的旋转取向都会改变。静电和流体动力学方面深远的边界约束效应是这些奇特现象的原因。所呈现的结果在微流控和纳流控操作以及药物递送应用中具有直接应用前景。

**论文解读文章**

**研究背景与问题**

液滴因其微小尺寸及内部循环涡流带来的优异传质与传热特性，在药物递送等实际应用中备受关注。例如，载药脂质体（liposome）在人体血液中迁移时，目标病变区域（如骨折组织）释放特定化学物质（如Ca²⁺、OH^-、H₂PO₄^-），形成溶质浓度梯度，从而触发扩散泳动（diffusiophoresis）——一种无需外加电场、无焦耳热效应的自导引迁移机制。圆柱孔是模拟人体微血管（直径1000-2500 nm）或生物组织细胞外基质（extracellular matrix, ECM）孔隙的理想几何模型。已有研究（Lee等）针对KCl溶液中介电液滴（dielectric droplet）在圆柱孔内的扩散泳动进行了分析，但由于K⁺与Cl^-扩散系数几乎相等，扩散势（diffusion potential）可忽略。然而，人体血液等生物流体中Na⁺和Cl^-是主要离子，二者扩散系数差异显著（Na⁺：1.334×10^-9 m²/s，Cl^-：2.032×10^-9 m²/s），会诱导出扩散势，从而产生额外的电泳分量（electrophoresis component），与原有的化学泳分量（chemiphoresis component）耦合。因此，研究人员开展了此项理论分析，旨在揭示扩散势存在时边界约束效应对液滴迁移行为的影响，并对比弱荷与强荷液滴的情况，为微流控、纳流控及药物递送提供基础理论指导。该论文发表在《Colloids and Interfaces》。

**主要关键技术方法**

研究人员采用基于切比雪夫多项式（Chebyshev polynomials）的拼接伪谱法（patched pseudo-spectral method），结合几何映射方案（geometric mapping scheme）处理液滴与圆柱孔之间的不规则求解域，求解了耦合的泊松方程（Poisson equation）、离子输运方程（Nernst-Planck方程）和流体动量方程（Navier-Stokes方程）。通过无量纲化处理，利用边界条件（包括恒定表面电荷密度条件、扩散势相关参数β）计算液滴表面的水动力拖曳力（hydrodynamic drag force）和电驱动力（electric driving force），进而获得液滴的无量纲迁移率（mobility）。数值方案经网格收敛性检验（51点θ方向×100点r方向）及单液滴（无限大介质）渐近行为验证，确保可靠。

**研究结果**

**迁移率反转与边界约束效应**
研究人员通过数值计算获得了不同粘度比σ_H（0.01、0.5、100）及不同表面电荷σ^*（2.03、10.53、-2.03）下液滴迁移率μ^*随孔半径比λ（液滴半径a与圆柱孔半径之比）的变化曲线。结果表明：在NaCl溶液中，所有带电液滴在λ较大时（如>5）均正向迁移（朝低浓度区域），且迁移率接近无限大介质中的值；但λ较小时（窄孔），出现迁移率反转（mobility reversal），即液滴运动方向与经典库仑静电力预期相反。例如，正电荷液滴在窄孔中可能朝高浓度区域运动，这是因为电泳分量不再是唯一主导因素，强烈的边界约束效应（静电与流体动力学双重作用）改变了电动力学场景。与KCl溶液对比，NaCl中扩散势的存在显著增大了迁移率幅值。

**固化现象与临界点**
在λ≈4.577处，所有粘度比下的液滴迁移率曲线相交于同一点，即无论液滴内部流体粘度如何，均以相同速度运动，表面剪切应力为零，内部无循环涡流。研究人员将该现象称为“固化现象”（solidification phenomenon），源于静电马克斯韦牵引力（Maxwell traction）与流体动力学扩散渗流（diffusioosmosis）驱动产生的旋转作用相互抵消。该临界点对脂质体设计具有重要意义，因为在血管运输中零剪切应力可避免膜破裂及药物提前释放。此外，当λ经过该临界点时，液滴表面的旋转取向发生改变：λ<4.577时，水动力旋转力占优，液滴向下运动（负迁移率），外部无轴涡流；λ>4.577时，静电牵引力逆转并产生外部轴涡流（内凹涡流环），形成流动再协调。

**强荷液滴与负电荷液滴**
对于高表面电荷密度（σ^*=10.53）的液滴，迁移率反转和固化现象依然存在，但迁移率幅值不同。对于负电荷液滴（σ^*=-2.03），在窄孔中也观察到与直观预测相反的极性依赖行为，进一步证实了边界约束效应的普遍性。研究人员指出，不能将无限大介质中单液滴的电动力学知识直接应用于圆柱孔内的液滴运动。

**总结讨论**
论文结论部分指出：在扩散势存在下（如NaCl溶液），圆柱孔内带电介电液滴的扩散泳动展现出奇特现象。窄孔中出现迁移率反转，且移动方向违背直觉预期。电泳分量不再是确定运动方向的唯一因素，而是由深远的静电与流体动力学边界约束效应共同决定。发现了两个临界点（λ值），液滴表面旋转取向在每次经过临界点时改变。此外，在特定条件下观察到“固化现象”，液滴如同刚性粒子运动，表面剪切速率为零。这些结果对微流控、纳流控及药物递送具有潜在应用价值。但研究局限于球形非变形液滴、固定表面电荷和理想圆柱约束，直接应用于实际时需考虑膜变形、电荷调节、复杂电解质、背景流及非均匀几何等。未来可扩展至不等式价电解质（如BaCl₂、LaCl₃）及多物种体系。

**研究结论翻译（来自论文第4节）**
扩散泳动（diffusiophoresis）中带电介电液滴（dielectric droplet）在圆柱孔（cylindrical pore）内的行为得到了理论分析，重点关注扩散势（diffusion potential）的影响，例如在NaCl电解质溶液中。发现了有趣的现象。在窄通道中观察到迁移率反转（mobility reversal）。此外，液滴运动方向违背了直观预测。由扩散势产生的电泳分量（electrophoresis component）并非决定液滴运动方向的唯一主导因素，这与单液滴在无限大电解质溶液介质中的通常情况不同。静电和流体动力学方面深远的边界约束效应（boundary confinement effect）是其原因。此外，找到了两个临界点，液滴表面旋转取向在每次经过这些临界点时发生变化。此外，在特定情况下还观察到“固化现象”（solidification phenomenon），液滴像刚性粒子一样运动，内部无循环涡流。所有液滴无论粘度如何都以相同速度运动，表明液滴表面剪切速率为零。当发生“固化现象”时，液滴表面的特征旋转运动完全消失。所呈现的结果在微流控（microfluidic）和纳流控（nanofluidic）操作以及药物递送（drug delivery）应用中具有潜在应用价值。然而，本研究侧重于球形非变形介电液滴的基本电动力学响应，具有固定表面电荷和理想圆柱约束。直接应用于药物递送等实际操作应谨慎，因为真实脂质体和生物孔道可能涉及膜变形、电荷调节、复杂电解质、表面吸附、背景流及非均匀几何。此外，作为可能的未来方向，本框架可自然地扩展到超出二元、电荷对称电解质（如本文处理的NaCl）的体系。由于扩散势是本研究的核心，具有不等化合价的电解质（如BaCl₂或LaCl₃）以及多物种电解质溶液可能导致更丰富的液滴迁移行为，正如近年来相关研究文献所表明的。