Inertial effects on the interphase drag force and rheology of dilute suspensions of buoyant droplets at low Reynolds number

本研究利用倒易定理证明，在低雷诺数下，浮力液滴稀疏悬浮液中的惯性效应会在相间拖曳力和连续相的有效应力中引入对相对速度和速度方差的二次依赖关系。

要点

想象一个拥挤的房间，每个人都试图朝同一个方向行走，但其中一些人是漂浮的气球（浮力液滴），而另一些人则是空气（连续流体）。通常，当我们研究这些气球如何在空气中移动时，我们会假设空气是完全静止的，且气球移动得非常缓慢，就像在浓稠的糖浆中移动的蜗牛一样。在那个缓慢的世界里，规则很简单：气球移动得越快，空气的反作用力就越大。

然而，在现实世界中，情况并不总是那么缓慢或简单。有时空气会带有一点点“冲劲”（惯性），而且气球可能会有些晃动，而不只是沿着直线移动。这篇由 Nicolas Fintzi 和 Jean-Lou Pierson 撰写的论文提出了一个具体的问题：当空气具有微小的速度，且气球在带着自身的能量跳动时，作用在这些漂浮气球上的力会发生什么变化？

以下是他们发现的拆解，使用了日常类比：

1. “倒易定理”就像一面魔镜

为了解决这个问题，作者并没有去模拟每一个空气液滴和每一个气球。那就像是为了了解潮汐如何移动而去数清海滩上的每一粒沙子一样。相反，他们使用了一个叫做**倒易定理（Reciprocal Theorem）**的数学工具。

你可以把它想象成一面魔镜。他们没有直接观察气球在微风中移动这一复杂且混乱的现实，而是观察了问题的一个“镜像版本”，在这个镜像中规则更简单（比如一个完全静止的房间）。通过将真实问题与这个简单的镜像进行比较，他们可以计算出复杂的力，而无需进行所有的繁重计算。这是一个捷径，让他们能够看到空气如何推拉气球的隐藏细节。

2. “抖动”很重要（速度方差）

在许多旧模型中，科学家假设所有气球的移动速度完全相同。但在现实中，有些气球可能漂移得更快，有些则更慢，有些甚至可能上下晃动。这种“抖动”或**速度方差（velocity variance）**就像一群人在走路：如果每个人都以完全相同的步调行走，那是井然有序的；但如果有人在冲刺，有人在漫步，人群就会产生一种不同的压力。

作者发现，这种“抖动”会产生额外的力：

• 阻力（Drag Force）： 空气不仅仅根据气球的平均速度产生阻力，它还会根据气球围绕平均速度进行抖动的程度产生阻力。
• 应力（“挤压”感）： 当你观察整组气球时，它们的抖动会在周围的空气中产生一种额外的“挤压”或压力。这就像一群人在紧张地挪动脚步；即使他们没有奔跑，他们的坐立不安也会在房间里营造出一种压力感。

3. “速度平方”效应

其中最重要的发现之一是，当气球移动得更快时，这些力是如何表现的。

• 在那个极慢的、类似糖浆的世界里，力与速度成正比（速度翻倍，推力也翻倍）。
• 在这个新的、稍快一点的世界里，力开始取决于速度的平方。

想象你在推一辆购物车。如果你轻轻推，很容易。如果你推得强两倍，感觉就不只是难了两倍；空气阻力和车轮与地面交互的方式会让它感觉难得多。作者表明，对于这些漂浮液滴，来自空气的“反推力”增长速度远快于速度本身，并且在很大程度上取决于液滴的抖动程度。

4. 为什么这改变了流体的“配方”

论文得出结论，如果你想描述空气和漂浮液滴组成的混合物如何表现（例如在气泡柱或浮选槽中），你不能只使用旧的、简单的配方。

• 旧配方： “根据液滴的数量增加粘度（厚度）。”
• 新配方： “增加粘度，但也要增加一个取决于液滴相对于空气移动速度的项，以及另一个取决于它们抖动程度的项。”

这意味着这种混合物表现得不像一种简单的粘稠液体（如蜂蜜），而更像是一种智能材料，它的行为会根据液滴移动的速度和混乱程度而改变。

总结

简而言之，Fintzi 和 Pierson 利用一个聪明的数学镜像证明了，当漂浮液滴在具有一定速度的流体中移动时：

1. 惯性很重要： 流体的“冲劲”改变了阻力。
2. 抖动很重要： 液滴之间随机的速度差异会产生额外的力和压力。
3. 非线性行为： 力不仅随速度增长，还随速度的平方和抖动的平方而增长。

这有助于工程师理解，要预测这些混合物（如工业分离罐中的混合物）如何流动，他们需要考虑到液滴的“坐立不安”，而不仅仅是它们的平均速度。

技术摘要

技术摘要：低雷诺数浮力液滴悬浮液中的惯性效应对界面阻力和流变性的影响

问题阐述
本研究旨在对由浮力液滴组成的稀薄两相分散流进行建模，其处于低但有限的雷诺数（$Re = aU\rho_f/\mu_f$）和稀薄体积分数（$\phi \ll 1$）状态。尽管在斯托克斯流（粘性主导）和势流（无粘）极限领域已有大量文献，但在低但有限的 $Re$状态下，针对浮力液滴的高阶流体力学力矩的闭合模型仍然匮乏。具体而言，本研究旨在量化流体惯性如何修正界面阻力以及力的第一和第二阶矩。这些力矩对于闭合两流体模型中的平均动量方程至关重要，因为它们贡献了连续相的有效应力以及两相间的动量交换。作者重点研究了平移球形液滴的具体构型，其中相对速度诱发了惯性，且表面张力足以维持球形形状（$Ca \ll 1$）。

研究方法
作者结合倒易定理（reciprocal theorem），采用严谨的统计平均程序来推导闭合关系。

1. 平均方程： 研究首先定义了分散相和连续相的系综平均质量和动量方程。闭合问题涉及将系综平均项（如界面力 $F$ 和有效应力 $\Sigma^{\text{eff}}$）表示为宏观变量（数密度、体积分数、平均速度）的形式。
2. 条件平均： 为了解决闭合问题，作者利用了条件平均场。他们考虑了一个具有特定质心速度 $w$ 的“测试液滴”浸没在背景流中。围绕该液滴的扰动场的控制方程是通过对纳维-斯托克斯方程进行条件平均而得到的。
3. 广义倒易定理： 本研究的一个核心方法论贡献是推导出了适用于液滴的广义倒易定理。该定理将作用在任意流场中液滴上的力矩与已知的辅助解（斯托克斯流解及奇点解，如点源和点力）联系起来。不同于以往仅关注力或应力张量（stresslet）的工作，该公式系统地计算了任意阶的力矩。
4. 摄动分析： 作者将此定理应用于均匀稳态流中的液滴。他们计算了阻力、第一矩（应力张量）和第二矩的 $O(Re)$惯性修正。这涉及匹配内层（斯托克斯）和外层（奥森）解，以处理有限$Re$ 下扰动场的不均匀性。
5. 系综平均： 最后，通过对液滴速度分布进行系综平均，得到导出的力矩。这一步明确地将依赖于分散相速度方差（$\langle \delta_p u'_\alpha u'_\alpha \rangle$）的项引入到宏观方程中。

主要结果
研究给出了流体力学力和力矩惯性修正的显式解析表达式：

• 阻力： $O(Re)$阶阻力修正的量级为$O(\rho_f \phi U^2)$。结果恢复了 Taylor 和 Acrivos (1964) 关于平移液滴阻力的经典解，但其是在倒易定理和系综平均的框架下推导出来的。
• 第一矩（应力张量）： 力的第一矩惯性修正的量级为 $O(\rho_f \phi U^2)$。值得注意的是，作者发现只要液滴保持球形，液滴密度比（$\zeta = \rho_p/\rho_f$）不会出现在第一阶和第二阶矩的领先阶惯性修正中。这与阻力的情况不同，在其他语境下，密度会通过变形进入阻力计算。
• 第二矩： 第二矩的 $O(Re)$修正量级为$O(a \rho_f \phi U^2)$。作者推导出了完整的张量表达式，包括其迹（trace）和无迹部分。
• 速度方差贡献： 一个关键发现是，通过液滴速度分布对力矩进行系综平均，会引入与分散相速度方差成正比的附加项。具体而言，平均阻力和连续相的有效应力取决于相对速度的平方以及速度方差。
• 有效应力： 得到的连续相有效应力张量包含了正比于 $\rho_rho_f \phi U_r U_r$ 和 $\rho_f \phi \langle u' u' \rangle$ 的项。这些项表明该悬浮液表现出非牛顿行为，其有效应力取决于相对速度的平方以及体积分数和相对速度的空间梯度。

意义与主张
作者将这项工作定位为实现更精确的两流体建模的必要步骤，特别是在浮选和液液萃取等存在浮力和湍流共存的工业过程中。

• 框架的新颖性： 本文声称是首个在平均两相流方程框架内，研究低但有限 $Re$ 数下浮力液滴，并明确包含来自第一和第二力矩的流变学贡献的研究。
• 方法论的严谨性： 通过推导广义倒易关系，作者提供了一个系统性的工具，用于计算任意阶的力矩，扩展了以往用于固体颗粒或势流的方法。
• 闭合意义： 研究强调，从孤立颗粒推导出的标准阻力定律在应用于平均框架时必须进行修正。对液滴速度分布的平均过程自然地产生了高阶矩（速度方差项），这些项在以往的研究中常被忽略或仅进行经验建模。作者证明了这些项是基于平均过程的物理基础。
• 流变行为： 该工作表明，在低 $Re$ 数下，浮力液滴的稀薄悬浮液表现为二阶梯度流体（或在某些极限下为 Reiner-Rivlin 流体），其有效应力取决于速度梯度和方差。然而，作者也提醒道，虽然第一矩的惯性贡献是显著的，但由于尺度分离假设（$a/L \ll 1$），第二矩的惯性贡献在实际模型中可能可以忽略不计。
• 局限性： 作者谦虚地指出，该方程组并未完全闭合，因为尚未针对速度方差项本身提出明确的模型，尽管他们引用了现有的动力学理论和湍流建模方法作为潜在的闭合手段。此外，结果假设液滴为球形且界面洁净，未考虑马兰戈尼（Marangoni）应力和显著的液滴变形。