Effects of surface viscosities on the motion of a droplet enclosing a translating particle

本研究提出了精确解析解与数值分析，表明表面扩张粘度虽可调控包裹平移粒子的同心液滴的运动，但表面剪切粘度却因对称性破缺而独特地增强了偏心构型中的运动，共同揭示了界面流变学、受限效应与几何形状如何支配复合粒子的动力学行为。

要点

想象一个微小而不可见的微观世界：一颗实心大理石被囚禁在一团漂浮的油泡内部。现在，设想有人从内部推动这颗大理石，试图让它游向外部的水中。那么，油泡会发生什么变化？它是静止不动，还是会被拖拽着一起移动？

本文正是探讨这一情景，但有一个转折：油泡的表面并非一层简单、滑溜的薄膜，而是覆盖着一种“分子蜂蜜”或“粘性薄膜”，能够抵抗拉伸和滑动。研究人员希望了解，这种粘性薄膜如何改变当内部大理石推动油泡时，油泡的运动方式。

以下是他们研究发现的通俗解读：

设置：大理石与油泡

将整个系统想象成一个俄罗斯套娃，但由流体构成。

• 内层套娃：一颗被以恒定速度推动的实心刚性大理石（即颗粒）。
• 外层套娃：包裹着大理石的一滴液体（即油泡）。
• 表皮：油泡的表面具有特殊性质。它拥有表面剪切粘度（抵抗侧向滑动，就像试图在地板上拖动一块厚重的地毯）和表面膨胀粘度（抵抗拉伸或收缩，就像试图吹胀一个非常厚实、坚硬的气球）。

完美居中的情况（“同心”设置）

首先，研究人员考察了大理石恰好位于油泡正中心的场景。

• “滑动”阻力（剪切粘度）：令人惊讶的是，如果大理石完美居中，油泡表皮的“滑动”粘性完全不起作用。对于这种特定设置而言，表皮仿佛完美光滑。无论油泡表皮抵抗滑动的程度如何，油泡的移动速度都保持不变。

• “拉伸”阻力（膨胀粘度）：这里情况变得复杂。“拉伸”粘性确实会改变结果，但它的作用如同拔河，涉及两种相互对抗的力：

  1. 刹车：粘性表皮使油泡更难移动，如同刹车片。
  2. 引擎：由于内部大理石被以固定速度推动，表皮越粘，大理石为了保持运动就必须施加更大的推力。这种额外的推力实际上有助于拖拽油泡前行。

  结果：取决于大理石在油泡内的贴合紧密程度以及流体的粘度，“刹车”可能会占上风（使油泡减速），或者“引擎”可能会占上风（使油泡加速）。这是一种微妙的平衡。

偏心情况（“偏心”设置）

接下来，研究人员将大理石移开，使其不再位于中心，而是更靠近油泡的一侧。

• “滑动”阻力回归：突然间，“滑动”粘性（剪切粘度）变得重要了！当大理石偏离中心时，油泡表皮开始以某种方式滑动，从而产生新的效应。
• 助推：在这种偏心位置下，滑动粘性实际上帮助油泡移动得更快。这就像摩擦力现在转而为你所用，给油泡额外的推力。大理石偏离中心的程度越大，这种助推效应就越显著。
• 主导力量：然而，如果你同时存在两种粘性（滑动和拉伸），“拉伸”效应通常占据主导地位。它决定了速度，而“滑动”助推则成为一个较小、次要的细节。

全局视角

研究人员利用高级数学和计算机模拟证明了这些观点。他们发现：

1. 对称性是关键：当系统完美平衡（居中）时，一种类型的粘性会从方程中消失。
2. 失衡产生新力：当系统失衡（偏心）时，那种“缺失”的粘性会重新出现，并实际上有助于运动。
3. “粘性”表皮是一把双刃剑：它既可以像刹车一样使系统减速，也可以通过迫使内部大理石施加更大推力来加速系统。

简而言之，这篇论文揭示出，液滴的“表皮”并非仅仅是被动的包裹物。根据内部物体所处位置的不同，该表皮可以充当刹车、引擎或助推器，从根本上改变整个系统在流体中的运动方式。

技术摘要

技术摘要：表面粘度对包裹平移粒子的液滴运动的影响

问题表述
本研究考察了斯托克斯流 regime 下，由刚性球形粒子（半径 $a$）包裹在粘性球形液滴（半径 $b$）内构成的复合粒子系统的流体动力学。该系统浸没在无界的牛顿流体中。刚性粒子以预设速度 $U_P$ 沿 $z$ 轴驱动，诱导包裹液滴以未知速度 $U_D$ 发生平移运动。研究重点在于界面流变学，特别是表面剪切粘度（$\eta$）和表面膨胀粘度（$\kappa$）如何影响这一运动。液滴界面采用 Boussinesq–Scriven 本构定律进行建模。研究考虑了同心构型（粒子居中，偏心率 $\epsilon = 0$）和偏心构型（粒子发生位移，$\epsilon > 0$）。分析假设毛细数较小（$Ca \ll 1$），允许液滴保持球形并忽略形状变形。

方法论
作者采用结合精确解析解和数值模拟的双重方法：

1. 解析解（同心情况）： 对于同心构型（$\epsilon = 0$），作者利用球坐标系下的流函数公式推导出了精确解析解。斯托克斯方程被简化为双调和方程，未知系数通过满足速度连续性、应力平衡（包含 Boussinesq–Scriven 模型）以及全局力平衡的边界条件来确定。
2. 谱边界积分法（偏心情况）： 为了解决解析解难以处理的偏心几何构型，作者利用三维谱边界积分方法。该方法将控制斯托克斯方程重写为仅定义在粒子表面（$\Gamma_P$）和液滴界面（$\Gamma_D$）上的边界积分方程。表面使用实球谐函数进行参数化，生成的线性方程组采用广义最小残差法（GMRES）求解。
3. 验证： 数值实现针对同心情况的精确解析解进行了验证，并利用源自洛伦兹互易定理的一般积分关系进行了核实，该关系将系统上的总流体动力与面积平均表面速度联系起来。

主要贡献与结果

• 同心构型：

  • 与剪切粘度的无关性： 解析解表明，诱导的液滴速度 $U_D$ 完全独立于表面剪切粘度（$Bq_\eta$）。这一发现与之前关于平移液滴和具有表面粘度的同心自泳体（squirmers）的观察结果一致。
  • 膨胀粘度的作用： 液滴速度取决于表面膨胀粘度（$Bq_\kappa$），但其影响是非单调的，且取决于限制比（$\alpha = b/a$）和粘度对比度（$\lambda = \mu_i/\mu_e$）。
    • 在紧密限制系统中（例如 $\alpha = 1.5$），增加 $Bq_\kappa$ 会增强 $U_D$。
    • 在弱限制系统中（例如 $\alpha = 10$），增加 $Bq_\kappa$ 通常会抑制 $U_D$。
    • 中等限制表现出临界粘度比阈值，在此阈值处趋势发生逆转。
  • 机制： 作者将这些趋势合理化解释为两种效应之间的竞争：(i) 界面阻力增加（迁移率降低）阻碍运动；(ii) 为了克服增加的界面阻力以维持固定的粒子速度 $U_P$，需要更大的驱动力。在紧密限制下，增加的驱动力占主导地位，从而增强 $U_D$。在弱限制下，迁移率降低占主导地位，从而抑制 $U_D$。
  • 迁移率解释： 结果表明，表面膨胀粘度主要起到有效增加内部流体粘度的作用，特别是在弱限制系统中。

• 偏心构型：

  • 剪切粘度依赖性的出现： 与同心情况不同，偏心构型表现出对表面剪切粘度的依赖性。表面剪切粘度的存在始终增强了诱导的液滴速度 $U_D$。
  • 偏心度与不对称性： 随着偏心度 $\epsilon$ 的增加，由剪切粘度引起的增强效应变得更加显著。偏心情况下的流场表现出前后不对称性，在粒子与界面之间的狭窄间隙中速度幅值更高。
  • 综合效应： 当同时存在剪切粘度和膨胀粘度时，液滴速度的修正主要由膨胀贡献决定，使得剪切粘度的附加效应相比之下较弱。

意义与主张
本文建立了具有复杂界面的活性复合粒子动力学的定量基准。作者声称，他们的发现提供了关于界面流变学、几何限制和对称性破缺如何共同支配系统动力学的机理见解。具体而言，该工作阐明了表面粘度如何根据具体的几何和流变学参数增强或抑制液滴运动。精确的解析解和经过验证的数值工具被呈现为未来涉及更复杂构型（如非轴对称运动、旋转夹杂物或可变形液滴）研究的重要参考。作者明确指出，这项工作代表了探索具有复杂界面的活性复合粒子动力学的第一步，并未提出具体的新实验应用，也未将范围扩展到粘弹性效应（他们指出这是未来的工作方向）。