Self-similar breakup of a liquid ligament with a solid particle

想象一下，你有一根长长的、细细的蜂蜜或浓糖浆丝悬挂在空中。如果你将两端拉开，这根丝会变得越来越细，直到最终断裂，分裂成独立的液滴。这是自然界和科技中常见的景象，从树叶上滴落的雨滴，到喷墨打印机喷射出的微小墨点，皆属此类。

通常，这种断裂发生是因为液体本身具有不稳定性；它倾向于转变为球体（液滴）以节省能量。但如果在那根粘稠的丝线内部困住了一粒微小的固体颗粒——比如一粒沙子或一颗尘埃——会发生什么呢？

本文正是研究这一情景。研究人员利用计算机模拟和数学方法，观察单个固体颗粒如何改变拉伸液体丝线的断裂方式。

以下是他们发现的简述，分解为几个简单概念：

设置：带有“结”的拉伸绳

将液体丝线想象成一根由蜂蜜制成的长而富有弹性的绳子。研究人员以恒定速度拉开绳子的两端。在绳子内部，他们放置了一个单一的固体球（即颗粒）。

起初，绳子很粗，而球体相对于绳子的宽度来说很小。这就像在一根粗花园水管里放了一颗弹珠。弹珠几乎不起作用；绳子随着拉伸变得越来越细，遵循可预测的模式。

转折点：当“水管”缩小到与“弹珠”相当时

随着绳子继续拉伸，它变得越来越窄。最终，绳子变得如此细，以至于几乎要触碰到内部弹珠的表面。

这是关键时刻。论文将这一时刻称为颗粒尺寸与绳子尺寸之比接近 1 的时候。突然间，弹珠在绳子中起到了“结”或“凸起”的作用。由于绳子非常细，这个凸起产生了一个局部扰动。

惊喜：“普适性”断裂

这是发现中最有趣的部分。研究人员用不同大小的弹珠（有些小，有些大）进行了测试。

断裂前：较大的弹珠使绳子比小弹珠更早断裂。这很合理；更大的障碍物会更早引发问题。
断裂时：一旦绳子细到足以接触弹珠，某种神奇的事情发生了。最终断裂发生的速度变得完全相同，无论弹珠是小还是大。

研究人员称此为“自相似”行为。仿佛一旦绳子细到足以接触障碍物，障碍物的具体尺寸就不再重要。液体“忘记”了颗粒的大小，并遵循一条普适的、可预测的路径走向断裂。

类比：交通堵塞

想象一条高速公路（液体丝线），车辆彼此远离行驶，使交通流分散（拉伸）。

早期阶段：如果路中间有一个小坑（小颗粒）或一块大石头（大颗粒），由于道路宽阔，这暂时影响不大。
晚期阶段：随着道路收窄至单车道，小坑和大石头都变成了巨大的障碍。
断裂时刻：当交通被挤压到撞上障碍物的那一刻，无论是小坑还是大石头，交通瘫痪并停止（即“断裂”）的方式完全相同。障碍物的尺寸不再改变最终瘫痪的时机；唯一重要的是有某种东西存在。

数学与物理

研究人员不仅观察了这一现象，还写出了一个数学公式来精确预测断裂何时发生。

他们发现，断裂时间取决于两种力之间的较量：拉伸（将绳子拉开）和粘度（液体的“厚度”或粘性）。
在粘稠的液体中（如我们类比中的蜂蜜），“粘性”占主导地位。
他们的公式成功预测了断裂时间，与计算机模拟结果完美吻合。

核心结论

该论文得出结论：虽然颗粒会改变断裂开始的时间（通过在颗粒附近使绳子更快变细），但一旦绳子细到足以接触颗粒，最终的断裂过程遵循一条普适规则。

在这种特定的“粘性”机制下，液体丝线表现得像一台拥有固定断裂程序的机器。一旦颗粒靠近到足以触发该程序，颗粒的大小就变得无关紧要，丝线每次都会以可预测的、自相似的方式断裂。

技术摘要：含固体颗粒的液体细丝自相似断裂

问题陈述
液体细丝的变薄断裂是涵盖喷雾形成、喷墨打印至微流体液滴生成等诸多应用中的基本过程。在没有外部拉伸的情况下，这种断裂由瑞利 - 普拉托（Rayleigh–Plateau）不稳定性主导，其中表面扰动因毛细力而增长。其动力学特征通常由奥内佐格数（$Oh$）表征，该数代表粘性力与惯性力和表面张力力之比。虽然洁净细丝的行为已得到充分理解，但许多实际系统涉及含有固体颗粒或杂质的悬浮液。先前的研究考察了含多颗粒或浓悬浮液的射流，证明了断裂模式的改变。然而，单个固体颗粒如何作为局部扰动来改变拉伸细丝的断裂动力学和 pinch-off（颈缩断裂）时间的基本机制，仍知之甚少。具体而言，尚不清楚细丝拉伸、粘度与局部颗粒诱导扰动之间的相互作用如何决定最终的断裂机制。

方法论
作者采用数值模拟与解析建模相结合的方法来研究该问题。

数值设置：研究利用基于彩色梯度法（color-gradient approach）的两相不混溶流体格子玻尔兹曼方法（LBM）。系统模拟了一个初始半径为 $R_i$ 、长度为 $L$ （长径比 $L/R_i = 28.8$ ）的圆柱形液体细丝，周围为空气。一个半径为 $r_p$ 的固体球形颗粒被放置在细丝中心，并施加无滑移边界条件。
拉伸机制：通过在细丝末端施加线性速度剖面（ $u_x = \frac{2}{L}(x - L/2)u_{max}$ ）来施加轴向拉伸，导致质量流出计算域，细丝半径随时间减小。拉伸速率由无量纲流出率 $\tilde{\dot{\epsilon}}$ 表征。
参数：模拟聚焦于粘性机制区域，其中 $Oh \sim 1.0$ 。颗粒尺寸相对于初始细丝半径（ $r_p/R_i$ ）进行变化，并调整流出率以控制拉伸时间尺度相对于毛细 - 粘性时间尺度（ $\tau_v = \mu R_i / \sigma$ ）的关系。
解析方法：作者通过建立解析表达式来推导 pinch-off 时间，将细丝半径建模为无扰动平均半径（由质量守恒控制）与颗粒诱导的局部扰动振幅之和。他们假设扰动根据瑞利 - 普拉托不稳定性动力学增长，在 $Oh \sim 1$ 机制下由粘性效应主导。

关键结果

两阶段动力学：断裂过程表现出两个截然不同的阶段。最初，当细丝半径 $R(t)$ 显著大于颗粒半径（ $r_p/R(t) < 1$ ）时，细丝根据质量守恒单调变薄，颗粒的影响可忽略不计。
自相似性的 onset：随着拉伸进行，细丝表面接近颗粒（ $\beta(t) = r_p/R(t) \sim 1$ ），触发局部扰动。一旦达到这种邻近状态，随后的扰动增长和最终的 pinch-off 动力学变得自相似且与颗粒尺寸无关。
通用 pinch-off 时间：虽然断裂的总时间（ $t_p$ ）取决于颗粒尺寸（较大的颗粒触发更早的断裂），但自相似 pinch-off 时间（ $t_{ps} = t_p - t_{exp}$ ），定义为从扰动增长 onset 到断裂的持续时间，在给定的拉伸速率下，对于不同的颗粒尺寸会坍缩为一个单一值。
解析验证：作者基于拉伸速率与粘性毛细不稳定性增长率（ $\lambda = \sigma / \mu R(t)$ ）之间的竞争，推导出了 pinch-off 时间（ $t_{p,a}$ ）的解析表达式。该表达式 $t_{p,a} = \frac{1}{\lambda_0 e^{\dot{\epsilon} t_{p,a}} + \dot{\epsilon} \ln(R_i/a_0)}$ 显示，在不同颗粒尺寸和流出率下，与数值模拟数据具有出色的定量一致性。

意义与主张
本文声称揭示了一种通用机制，即局部扰动如何控制粘性机制中含固体颗粒的细丝断裂。其主要贡献在于证明：一旦细丝表面接近固体颗粒，颗粒的具体尺寸便与最终的断裂动力学无关；系统进入一个通用的自相似机制。

作者强调，这种通用性源于粘性时间尺度（ $t_\mu = \mu R / \sigma$ ）允许颗粒诱导的扰动持续并增长，这与惯性机制（ $Oh \ll 1$ ）不同，在后者中断裂发生得太快，以至于此类扰动无法显著改变结果。所提出的解析框架成功捕捉了 $Oh \sim 1$ 时的这种通用标度机制，提供了一种预测 pinch-off 时间的工具，该时间取决于拉伸速率和流体属性，但在扰动被触发后与具体的颗粒尺寸无关。该研究局限于对称的、静止的颗粒构型，作为隔离这些效应的最小化设置，更复杂的、被平流的颗粒相互作用留待未来研究。