🔬 physics

A generalized, quasi-universal subgrid model for macroscopic moving contact line flows

本文提出了一种广义的、准普适的亚网格模型，该模型将考克斯（Cox）的匹配渐近关系与经验接触角模型相结合，旨在实现对宏观移动接触线流动的预测性、网格无关性模拟，且无需进行针对特定案例的校准或引入唯象参数。

原作者： Vyankatesh Manoj Mundhada, Manoj Kumar Tripathi

发布于 2026-09-29✓ Author reviewed ⓘ

📖 1 分钟阅读☕ 轻松阅读

CC BY 4.0

原作者： Vyankatesh Manoj Mundhada, Manoj Kumar Tripathi

原始论文采用 CC BY 4.0 许可（http://creativecommons.org/licenses/by/4.0/）。 ✨ 这是对下方论文的AI生成解释。它不是由作者撰写的。如需技术准确性，请参阅原始论文。阅读完整免责声明

想象一滴水落在平面上。随着它的铺展，液体、空气和固体交汇的边缘——即接触线（contact line）——并不会像汽车轮胎在柏油路上行驶那样简单地滑动。相反，这个边界的行为挑战了用于描述大多数流动液体的标准流体力学规则。几十年来，科学家们一直难以在计算机上模拟这种运动，因为液滴边缘的物理过程涉及的尺度极其微小——甚至达到了单个分子的尺寸——试图计算每一个细节所需的计算能力将超过现有的总和。为了绕过这个问题，研究人员使用了“子网格模型”（subgrid models）作为捷径，通过猜测微观边缘发生的情况来处理问题。然而，这些捷径存在一个重大缺陷：它们产生的结果往往取决于计算机将空间划分为多细的网格，这使得模拟结果对于现实世界的预测而言是不可靠的。

印度科学教育研究学院（Indian Institute of Science Education and Research, Bhopal）的一支研究团队开发了一种处理此问题的新方法。他们创建了一个模型，使计算机能够在不需要了解那些无法测量的微观细节的情况下，准确地模拟移动接触线。他们的方法弥合了液滴边缘微观且不可见的物理现象与宏观可见流体行为之间的鸿沟。通过这种方式，他们产生的模拟结果无论计算机的网格大小如何，都能保持一致性和准确性，这是实现从喷墨打印到石油开采等各种领域中液体行为预测的关键一步，且无需针对每种新情况进行调整。

问题的核心在于流体运动方式与它们接触固体表面方式之间的冲突。在标准流体力学中，一个被称为“无滑移条件”（no-slip condition）的规则规定，接触固体表面的流体在边界处必须完全停止运动。虽然这在大多数流动中表现良好，但在移动液滴的接触线处，它会失效，从而产生一个力变为无穷大的数学奇点。为了解决这个问题，科学家们在历史上一直假设流体在表面会发生极微小的滑动，但这种滑动的距离如此之小，以至于无法直接测量。因此，以往的计算机模型必须猜测这些微观数值，或依赖于拟合参数来匹配实验数据。这意味着，如果研究人员想要模拟不同类型的液体或不同的速度，他们通常必须重新开始并重新校准模型，这剥夺了模型预测新场景的能力。

这项研究中的研究人员意识到，以往模拟中的不一致性源于对液滴角度随远离接触线而变化方式的误解。他们借鉴了20世纪80年代一项成熟的理论，该理论描述了由于表面张力与流体内部摩擦力之间的竞争，液体界面形状是如何弯曲的。该理论表明，你观察到的角度取决于你距离边缘有多远。研究人员将这种理论理解与现有的经验公式（即通过观察真实液滴得出的规则）相结合，创造了一个新的统一模型。该模型不再强迫计算机去猜测微观滑动长度或墙壁处的精确角度，而是根据计算机网格单元的大小来计算应使用的正确角度。

为了测试他们的想法，该团队运行了一系列关于液滴在平板上铺展的计算机模拟。他们将新方法与旧技术进行了对比，网格大小从粗糙到非常精细不等。在旧方法中，改变网格大小会显著改变结果，导致模拟的液滴以不同的速率铺展或呈现出不同的形状。而在新模型下，结果则保持几乎完全相同。无论计算机使用的是粗网格还是细网格，模拟的液滴行为都保持一致。这证明了该模型已成功消除了对网格的依赖，即所谓的“网格依赖性”现象。

研究还探讨了该模型在各种条件下的表现，例如不同的流速和不同类型的液体。他们测试了液体向外铺展（前进）和向内收缩（后退）的情景。在每种情况下，新模型都产生了稳定的、与网格无关的结果。一个关键发现是，该模型不需要源自特定实验的“拟合参数”即可工作。一旦研究人员确定了原始经验公式所依据的尺度，该模型就能在无需进一步调整的情况下预测流动的行为。这与以往的方法有着显著区别，后者通常需要为每一组新条件进行独特的校准。

研究人员还调查了为什么一些旧的模拟在过去看起来效果不错。他们发现，在某些特定情况下，特别是气液系统中，网格大小引入的误差本身就足够小，以至于未被察觉。此外，如果计算机网格恰好与原始实验数据收集时的特定尺度相匹配，结果就会在偶然间与现实吻合。然而，这纯属运气，而非稳健的解决方案。新模型消除了这种偶然性，提供了一个在广泛的物理条件下均能可靠工作的框架，包括不同的密度和粘度。

这项工作最具有前景的方面之一是其预测能力。由于该模型不依赖于已知微观滑动长度或精确壁面接触角（这些量是极难测量的），因此它可以应用于已知这些数值缺失的新问题。研究人员展示了，只要知道原始实验观察所依据的尺度，该模型就可以被普遍应用。这意味着工程师和科学家可以利用这些模拟来设计更好的涂层系统、管理多孔岩石中的流体流动或理解肺部的生物过程，并确信其结果并非来自计算机网格的人为产物。

研究结论指出，虽然该模型的阶近似法对于大多数角度都适用，但涉及极大接触角的更复杂场景可能需要稍细致的数学修正。然而，将微观理论与宏观观察相结合的基本方法已被证明是成功的。通过将接触线视为一个响应观察尺度的动态特征，而非一个具有神秘属性的固定点，研究人员提供了一个既准确又具有广泛适用性的工具。这项工作代表了从仅仅拟合数据向理解底层物理学的转变，为实现能够真正预测现实世界中液体行为的模拟提供了路径。

技术摘要：一种用于宏观移动接触线流动的广义、准通用亚格点模型

问题陈述
数值模拟移动接触线（MCL）流动面临着根本性的挑战，即在固体表面上由于无滑移边界条件导致的应力奇异性，以及该现象固有的多尺度特性。为了消除这种奇异性而解析所需的微观长度尺度（纳米级滑移长度），宏观模拟在计算上是难以实现的。因此，必须使用亚格点模型来弥合微观物理学与宏观分辨率之间的差距。然而，现有方法通常存在两个关键局限性：

网格依赖性： 当直接应用经验接触角模型而不考虑其推导时的尺度时，数值结果（特别是接触线位置和界面形状）会随网格分辨率的变化而显著变化。
缺乏预见性： 许多尺度不变模型需要给定现象学参数，例如微观滑移长度（ $\lambda$ 或 $s_m$ ）和微观壁面接触角（ $\theta_w$ ），这些参数在实验中或先验地都难以确定。这种对特定案例校准的依赖限制了模拟的预见能力。

方法论
作者提出了一种将 Cox 的匹配渐近理论与经验动态接触角模型相结合的准通用亚格点模型。核心方法包括以下步骤：

理论框架： 本研究利用 Cox 的流体动力学理论，该理论将表观接触角 $\theta(s)$ （距离接触线为 $s$ 处）与微观壁面角 $\theta_w$ 以及接触线毛细数（ $Ca_{CL}$ ）联系起来。该理论引入了对长度尺度比（ $s/s_m$ ）的对数依赖关系。
耦合策略： 为了避免使用未知的微观参数（ $\theta_w$ 和 $s_m$ ），作者推导出了一个消除这些变量的关系式。通过认识到经验模型（如 Kistler, Jiang, Tanner）提供的是在特定观测尺度 $s_0$ （通常是用于制定该模型的实验数据的分辨率）下的接触角 $\theta_0$ ，作者在经验模型与网格尺度之间建立了桥梁。
修正模型的推导： 作者在经验关系中的毛细数中引入了一个尺度相关的修正因子 $\chi(Ca_{CL}, s)$ 。利用 Cox 函数 $G(\theta)$ 的一阶泰勒展开，作者推导出了修正后的毛细数：
$\chi(s) = Ca_{CL} \left[ 1 + \zeta \ln\left(\frac{s}{s_0}\right) \right]$
其中 $\zeta$ 是从特定经验模型和粘度比中解析得出的参数。这使得能够施加一个与网格尺寸相关的接触角 $\theta(s)$ ，从而在不需要了解 $s_m$ 或 $\theta_w$ 的情况下满足 Cox 理论。
数值实现： 该模型在 Volume-of-Fluid (VoF) 框架内实现，使用的是 Basilisk 代码。接触线速度在每个时间步内与接触角和位置同时计算，以确保一致性，避免了与时间交错量相关的误差。该方法在广泛的雷诺数（ $Re$ ）和韦伯数（ $We$ ）范围内，针对轴对称液滴铺展和撤回板场景进行了测试。

主要结果

网格无关性： 所提出的“修正型”接触角模型成功消除了网格依赖性。在从 $\Delta_5$ 到 $\Delta_9$ （其中 $\Delta_i = D_0/2^i$ ）的网格上进行的模拟表明，接触线位置和液滴界面形状收敛于单一解。相比之下，使用原始（未经修正的）经验模型会导致显著误差（界面形状误差高达 ~1.5%，接触线位置误差更高），且这些误差不会随网格细化而收敛。
误差降低： 对于液-液系统（其网格依赖性最为严重），修正后的模型在接触线位置（ $E_r$ ）和界面形状（ $E_s$ ）的平均相对误差较原始模型降低了约一个数量级。例如，在 $Re=100, We=1, \theta_e=60^\circ$ 时，接触线位置误差从 ~0.69%（原始 Kistler）降至 ~0.05%（修正 Kistler）。
跨机制的鲁棒性： 该模型在变化的 $Re$ (1–200)、 $We$ (0.5–100) 和平衡接触角 ( $30^\circ$ – $135^\circ$ ) 下均表现出有效性。
高阶修正： 研究发现，对于较大的钝角（ $\ge 120^\circ$ ），修正因子 $\zeta$ 的一阶线性近似会引入更高的误差。研究推导出了二阶修正，并证明其在这些特定机制下能进一步降低误差，尽管一阶近似对于锐角和中等角度仍然足够。
气-液系统： 作者指出，对于气-液系统，原始数值方法通常表现出可以忽略不计的网格依赖性，无需修改，这可能是因为接触角随距离的变化很小，且网格尺寸偶然间与实验观测尺度相匹配。

意义与主张
本文声称提出了第一个能够在无需特定案例校准或先验了解微观参数（ $\theta_w$ 和 $s_m$ ）的情况下，进行移动接触线流动的预见性模拟的框架。

准通用性： 该模型被称为“准通用”的，是因为它仅依赖于 Cox 理论和所选经验接触角模型的有效性，而与特定的流体参数（如 $Re, We$ 或密度/粘度比）无关。
观测尺度依赖性： 该方法仅需要一次性了解与经验模型相关的观测长度尺度（ $s_0$ ）。一旦为特定模型固定了 $s_0$ ，即可将该亚格点模型应用于任何宏观流动问题。
消除拟合参数： 不同于以往需要通过拟合参数来匹配实验数据的尺度不变方法，本方法通过解析推导必要的修正，从而实现了真正的预见能力。
泛用性： 虽然已在二维轴对称几何中得到验证，但该公式是基于局部接触线描述推导的，并非仅限于轴对称几何，这表明其具有扩展到三维接触线的潜力。

总之，这项工作通过将流体动力学理论与经验数据严格耦合，解决了 MCL 模拟中网格依赖性的根本瓶颈，为可靠的、无需参数化的复杂多相流数值预测提供了一条路径。

类似论文