Symmetric structure-preserving discretization of N-phase incompressible fluid mixtures with arbitrary density ratios

本文提出了一种对称全离散数值方法，用于处理具有任意密度比的不可压缩 N 相 Navier-Stokes-Cahn-Hilliard 混合模型，该方法严格保持了关键物理特性，包括精确的质量与体积守恒、能量耗散以及相饱和约束。

要点

想象一下你正在观察一锅正在搅拌的汤，里面有油、水和醋在旋转。在现实世界中，这些液体并不会完美混合；它们会形成清晰的分层或液滴，并根据它们的重量和“粘性”相互推挤和拉扯。在计算机上模拟这一切是非常困难的，尤其是当你加入第三种液体（增加成分）或者当这些成分的重量差异巨大时（比如将沉重的蜂蜜与轻盈的空气混合）。

这篇论文介绍了一种用于模拟这些复杂流体混合物的全新“配方”。以下是作者所做工作的拆解，使用了简单的类比：

问题所在：“失衡的秤”

当科学家尝试模拟这些流体时，经常会遇到一个叫做“漂移”的问题。想象一下，一个本该保持完美平衡的秤，由于计算机微小的舍入误差，随着时间的推移，秤可能会慢慢倾斜，导致看起来质量在凭空消失或出现。

在具有不同密度的复杂混合物中，情况甚至更糟。如果计算机将一种液体视为“主角”，而将其他液体视为“配角”，模拟过程可能会产生偏差。它可能会在无意中偏袒某种液体而非另一种，从而破坏了现实世界的对称性。作者想要一种能够对每种液体都一视同仁的方法，就像一场每种相位都有平等投票权的民主选举，确保总体的“物质”（质量和体积）永远不会发生神奇的变化。

解决方案：一种“对称且能量诚实”的方法

作者创建了一个新的数学框架（一套计算机规则），它表现得像一个完美平衡的账本。

1. “平等地位”规则：
   大多数旧方法会选择一种液体作为“参考对象”（就像为一支队伍挑选队长）。本文的方法将所有 $N$ 种液体视为平等的伙伴。无论你有 3 种还是 10 种液体，数学处理它们的方式都是对称的。这防止了计算机在无意中偏袒某种液体而非另一种。

2. “无漂移”保证：
   作者证明了他们的方法可以保证三件事无论模拟运行多久都不会改变：

   • 总体积： 汤既不会膨胀也不会收缩。
   • 总质量： 没有液体会凭空消失或出现。
   • 个体质量： 油、水和醋的量始终保持完全一致（它们可以到处移动，但总量被锁定了）。

3. “能量银行”类比：
   把流体系统想象成一个银行账户。“能量”就是这个系统里的钱。在现实世界中，摩擦和混合总是需要消耗金钱（能量转化为热能而损失）。作者的方法确保计算机模拟的行为表现得像一家严格的银行：能量收支表始终是减少或持平的，绝不会无故增加。这被称为“能量耗散”，它能让模拟过程保持稳定且符合物理现实。

他们是如何做到的

为了实现这一点，作者必须重写计算机使用的方程组。

• “饱和”约束： 他们确保在模拟的每一个点，液体都填满了 100% 的空间（没有空隙）。如果液体开始完美地填充空间，数学逻辑保证它们在每一步都会持续完美地填充空间。
• “任意密度”特性： 以前的方法在处理重量差异极大的液体（例如重金属液体与轻气体的混合）时会遇到困难。这种新方法即使在密度比例极端的情况下也能完美运行。

证明：进行测试

作者不仅编写了数学公式，还通过三个场景进行了测试：

1. 收敛性测试： 他们检查了随着计算机“网格”变得更细，数学计算是否变得更精确。结果正如预测的那样，确实如此。
2. 相分离： 他们模拟了一个混乱的混合物分离成清晰液滴的过程。计算机正确地展示了液滴的形成过程以及能量如何逐渐降低，且没有任何“幽灵质量”出现。
3. 上升的气泡： 他们模拟了一个气泡在液体中上升的过程。他们将结果与已知基准进行了对比，发现其方法完美符合物理规律，精确地保持了气泡的体积。他们甚至模拟了一个气泡穿过两种不同液体层的过程，展示了该方法处理复杂多层交互的能力。

核心结论

这篇论文提供了一个稳健的、“对称的”工具，用于模拟复杂的流体混合物。它确保计算机模拟在每一步都遵循基本的物理定律（质量和能量守恒），即使在处理许多重量差异巨大的不同液体时也是如此。这就像是从一个漏水的桶升级到了一个密封且完美平衡的容器，用于你的流体模拟。

技术摘要

技术摘要：对称结构保持型 $N$ 相不可压缩流体混合物离散化方法

问题陈述
弥散界面模型，特别是 Navier–Stokes–Cahn–Hilliard (NSCH) 系列，通过将界面表示为光滑过渡层，被广泛用于模拟复杂流体的界面动力学。虽然针对二相（两相）流的数值方法已经非常成熟，但将其扩展到具有任意密度比的 $N$ 相不可压缩混合物时面临着重大挑战。现有的结构保持方法通常依赖于二相公式或区分一个参考相，这破坏了物理系统固有的对称性。

在 $N$ 相设定下的核心困难在于如何在离散层面同时满足几个相互竞争的要求：

1. 对称性： 数值方案必须对所有相进行对等处理，不得通过消除变量（通过饱和约束）或引入特权参考相来破坏对称性。
2. 守恒性： 方法必须保持相体积和质量，以及总体积和总质量。
3. 饱和约束： 体积份额之和在每个空间点和每个时间步都必须等于 1 ($\sum \phi_\alpha = 1$)，以防止出现非物理性的空隙或重叠。
4. 能量耗散： 离散解必须满足与连续介质热力学一致的离散能量耗散律。
5. 任意密度比： 即使在各相密度差异显著的情况下，该方法也必须保持稳定性和准确性，因为这种情况会引入动量平衡方程与质量平衡方程之间的强耦合。

方法论
作者提出了一种针对 $N$ 相不可压缩 NSCH 混合物的全离散有限元方法。该方法构建于统一的混合物理论框架之上，该系统由单个混合物动量方程和 $N$ 个相质量平衡方程控制。

• 等效重构： 为了便于结构保持，作者首先推导了连续方程的等效强形式。该重构利用特定的测试函数（权重）来导出能量演化。至关重要的是，他们修改了动量方程，消除了质量平衡权重中的非线性动能项，从而确保所有测试函数都属于适用于有限元离散的标准 Sobolev 空间。
• 对称处理： 该方案直接基于全套相变量 ($\phi_\alpha$) 和单个混合物速度 ($v$) 进行构建，无需消除任何相或引入参考相。这保持了各相之间的置换对称性。
• 离散策略：
  • 时间离散： 采用时间步进法构建全离散方案。对于自由能的无梯度部分，使用时间平均法（受 [37] 启发）处理，以确保凸性和稳定性。引入正则化项来处理密度和迁移率矩阵，以应对体积份额可能暂时脱离物理区间 $[0, 1]$ 的情况。
  • 空间离散： 采用使用 $H^1$ 包含空间的有限元近似。速度空间的选择需与压力空间满足 inf-sup 稳定性。
  • 关键结构特征： 该方案利用对流项的反对称形式以及质量和动量方程的特定离散权重。这使得推导精确的离散守恒律和离散能量耗散不等式成为可能。

主要贡献

1. 对称 $N$ 相离散化： 本文提出了一种针对 $N$ 相不可压缩 NSCH 模型的全离散方法，该方法在各相之间具有真正的对称性，避免了参考相公式引入的不对称性。
2. 精确的结构保持： 作者证明了每一个离散解都满足：
   • 精确的相体积守恒。
   • 精确的相质量守恒。
   • 精确的总体积和总质量守恒。
   • 离散能量耗散律。
   • 只要初始数据满足，则逐点保持体积饱和约束 ($\sum \phi_\alpha = 1$)。
3. 任意密度比： 该方法对于具有任意密度对比度的混合物具有鲁棒性，而在这种情形下，许多现有方案会因为密度与质量平衡之间的耦合而难以维持稳定性。
4. 混合物感知一致性： 该公式遵循连续介质模型的“混合物感知”特性，确保在相合并或某相不存在时，系统能够正确地进行简化（在单纯形面上进行归约）。

数值结果
作者通过以下几项数值实验验证了其理论特性：

• 收敛性测试 ($N=3$)： 一个具有三相和任意密度的测试表明，该方案达到了所选有限元空间预期的收敛阶，误差随网格细化而减小。
• 相分离 ($N=3$)： 三相分离模拟展示了相互连接的结构和三相交界处。结果证实了总能量是非增的，且质量/体积守恒误差处于机器精度水平。
• 上升气泡 ($N=2$)： 使用 Hysing 等人的标准上升气泡基准测试对方法进行了验证。气泡质心和上升速度的结果与尖锐界面代码及其他弥散界面公式的参考数据吻合良好。该方案精确保持了体积守恒。
• 三流体上升气泡 ($N=3$)： 模拟了一个复杂的场景，即气泡在两个分层液体层中上升。该方法能够正确捕捉物理行为：根据气泡半径相对于临界毛细半径的大小，气泡要么停留在界面处，要么穿透上层液体。

意义与主张
本文声称为具有复杂界面动力学和任意密度对比度的不可压缩 $N$ 相混合物流提供了一个鲁棒的计算框架。其主要意义在于能够在不牺牲对称性的情况下，在全离散层面保留连续方程定义的物理热力学和守恒结构。

作者强调，以往的方法往往会破坏对称性，或者在存在任意密度比时无法精确保持饱和约束。通过构建一种能对称处理所有相并逐点保持饱和约束的方案，这项工作填补了多相流数值模拟中的一个关键空白。结果表明，该方法对于代表性的多相流问题是稳定且准确的，为需要长时间积分和强耦合的系统模拟提供了可靠工具。

文章最后指出，虽然目前的工作建立了对称结构保持离散化的基础，但未来的工作可以将其思想扩展到高阶时间离散、更通用的混合物感知闭合以及可压缩 $N$ 相流。