Bubble detachment from circular cavities and flat surfaces

技术摘要：圆形空腔与平坦表面上的气泡脱离

问题陈述
本文研究了一个基础流体力学问题，即确定附着在水平表面上的气体气泡在脱离前所能达到的最大稳定体积。这一现象对于两个主要领域至关重要：一是提高电解水的效率（其中附着的气泡会阻碍离子传输并增加电阻）；二是理解自然发生气体的环境传输（例如海底甲烷渗漏）。虽然先前的解析模型（如 Tate 和 Fritz）及实验研究已经探讨了特定机制，但它们往往无法解释非球形界面、非对称不稳定性以及在全量程空腔尺寸和接触角下的不同附着模式之间的转换。作者试图将这些不同的案例统一到一个框架内，以回答三个具体问题：（Q1）脱离时的体积是多少？（Q2）气泡会长多宽？（Q3）它是如何脱离的？

研究方法
作者将该系统建模为一个简化、轴对称且附着在水平固体表面上的气泡，这在物理上等同于一个悬滴（pendant drop）。研究考虑了两种截然不同的附着模式：

1. 钉扎型气泡（Pinned bubbles）： 固定接触半径 ($r_0$)，接触角 ($\phi_0$) 可变。
2. 铺展型气泡（Spreading bubbles）： 固定接触角 ($\phi_0$)，接触半径 ($r_0$) 可变。

其控制物理过程由 Young–Laplace 方程描述，该方程平衡了表面张力、静水压力和重力，并以毛细长度 ($\lambda$) 进行缩放。作者使用打靶法（shooting method）对该方程进行数值求解。积分从气泡顶点（$r=0, z=h, \phi=\pi$）开始，并沿弧长采用 Adams–Bashforth 法进行。通过对广泛的初始顶点曲率半径 ($R_h$) 进行采样，作者生成了一系列平衡剖面。脱离判定基于两个标准：

• 轴对称不稳定性： 当 $\partial V/\partial h = 0$ 时发生，表明体积在给定高度下达到最大值。
• 非轴对称不稳定性： 当接触角达到 $\phi_0 = \pi$ 时发生，导致气泡沿空腔边缘发生侧向滑动。

主要贡献与结果
本研究提供了一个统一的数值框架，绘制了跨越整个接触半径（$0 < r_0/\lambda < \infty$）和接触角（$0 < \phi_0 < \pi$）参数空间的脱离行为图谱。

• 脱离体积：

  • 对于钉扎型气泡，脱离体积最初遵循经典的 Tate 预测（$V = 2\pi\lambda^2 r_0$）对于小半径的情况。然而，随着浮力使界面发生畸变，会出现偏差。钉扎型气泡的最大稳定体积在接触半径 $r_0 = 3.219\lambda$ 时为 $V = 18.96\lambda^3$。当 $r_0 > 3.219\lambda$ 时，气泡对非轴对称扰动变得不稳定，会发生侧向脱离，且脱离时的体积随半径增大而减小。当 $r_0 > 3.832\lambda$ 时，不存在稳定的解。
  • 对于铺展型气泡，体积随接触角单调增加，并在 $\phi_0 \to \pi$ 时趋于与钉扎情况相同的最大体积（$V = 18.96\lambda^3$）。这一结果低于 Fritz 在大接触角下的预测。
  • 研究结果表明，最大的钉扎型气泡和最大的铺展型气泡对应于 Young–Laplace 方程中完全相同的物理解。

• 脱离模式：
  本文识别了三种不同的脱离模式：

  1. 轴对称脱离（向上）： 发生在 $r_0 < 3.219\lambda$ 的钉扎型气泡和所有铺展型气泡中。
  2. 侧向脱离（向侧方）： 发生在 $3.129\lambda < r_0 < 3.832\lambda$ 的钉扎型气泡中，此时 $\phi_0$ 达到 $\pi$。
  3. 立即侧向脱离： 发生在 $r_0 > 3.832\lambda$ 时，此时无法形成稳定的气泡。

• 气泡尺寸与合并：
  作者推导了气泡最大水平半径 ($r_{max}$) 的经验拟合公式。对于钉扎型气泡，$r_{max}$ 在小半径时按 $r_0^{1/3}\lambda^{2/3}$ 比例缩放，在大半径时等于 $r_0$。对于铺展型气泡，$r_{max}$ 在小角度时随 $\phi_0$ 线性缩放，在大角度时呈二次方关系缩放。这些发现使得预测成核位点之间的合并距离成为可能。

• 脱钉转换（Unpinning Transition）：
  研究确立了在钉扎模式下成核的气泡（在半径为 $r_c$ 的空腔内）何时向铺展模式转换的判据。如果润湿角 $\phi_w$ 在接触角达到增长过程中的最小值之前就被达到，则会发生转换，具体条件为 $r_c < 3.5\lambda (\phi_w/\pi)^2$。

意义与主张
作者声称，其结果通过提供完整的数值精确解，考虑了非球形界面和稳定性极限，从而整合并扩展了先前的解析估计（Tate 和 Fritz）。研究结果已通过涵盖电解、沸腾、碳酸气泡以及悬滴的独立实验数据集进行了验证。

论文断言，其统一框架提供了：

• 精确的脱离判据，用于完整的电解槽模拟，有助于设计具有优化润湿性和空腔几何形状的电极，以促进气泡早期脱离。
• 定量指导，针对空间应用中的产气电极，指出降低重力（如在火星或月球上）会增加毛细长度，从而允许产生显著更大的稳定气泡。
• 环境气体传输模型的物理极限，特别指出从深海沉积物中脱离的单个甲烷气泡无法长到足以到达大气层的规模，因为在水中最大脱离直径约为 9.24 mm，低于深水逃逸所需的约 12.4 mm。

作者明确指出，由于对称性，其结果直接适用于悬滴，并且数值代码是免费提供的以支持进一步研究。他们也谦虚地提到，未来的工作可以扩展到涉及接触角滞后的非均匀表面模型。