Compressible and immiscible fluids with arbitrary density ratio

本文提出了一种基于热力学的新理论，用于对具有任意高密度比的可压缩且不混溶流体进行建模，通过推导能够正确解释真实动量演化的密度演化方程，解决了传统纳维-斯托克斯方程和欧拉方程的局限性。

要点

想象一下，你正在试图描述一场由两个截然不同的舞伴进行的舞蹈：一个沉重、移动缓慢的大象（水）和一个轻盈、移动快速的羽毛（空气）。在流体物理学的世界里，它们是“不互溶”的，这意味着它们不像油和水那样混合，但它们相遇的地方有一个模糊的边界。

长期以来，科学家们一直难以为这种密度差异巨大的系统（例如大象比羽毛重1000倍）编写“舞蹈规则”（即方程）。旧的规则有一个重大缺陷：它们将大象的重量视为一个常数、不变的数值，甚至在那个大象转化为羽毛的模糊边界处也是如此。这就像是在试图描述一个人变成幽灵的过程时说：“他们全程仍然是一个实实在在的人”，这显然是不合理的。

以下是本文旨在解决这一问题的简单分解。

1. 旧规则的问题

计算流体如何运动的传统方法（使用纳维-斯托克斯方程和欧拉方程）依赖于一种被称为 Boussinesq 近似 的捷径。

• 类比： 想象你正在推一辆购物车。如果车里装满了砖头，它就很重；如果车里只有空气，它就很轻。旧规则假设，如果你推的是一辆部分装有砖头、部分装有空气的购物车，那么在你推动的过程中，车的重量永远不会改变。它仅仅假设重量是一个固定的平均值。
• 缺陷： 在现实中，随着小车的移动和砖头的位移（扩散），重量会发生变化。旧规则忽略了“动量”（质量 $\times$ 速度）会随质量的变化而改变这一事实。它们还假设空气和水占据的空间是完全可预测的，忽略了当它们在边界处混合时，由于压力的影响，它们所占据的空间实际上会发生摆动和变化。

2. 新方法：能量最小化

作者并没有去猜测密度是如何变化的，而是从一个基本原理出发：自然界总是倾向于使用尽可能少的能量。

• 类比： 想象一个球从山上滚下。它并不关心自己具体迈了多少步，它只想到达山脚（能量最低点）。作者利用这种“能量山丘”的概念，推导出了水和空气如何相互作用的新规则。
• 核心创新： 作者引入了一个概念，称为 “过剩体积”（Excess Volume）。
  • 想象你有一桶水和一桶空气。如果你把它们倒在一起，你可能会预期总容量正好等于两桶容量之和。但在微观层面，当它们相遇时，分子可能会排列得更紧密或更松散，从而产生一点点“额外”或“缺失”的空间。
  • 旧规则假设这种额外的空间在任何地方都为零。而这篇论文指出：“不，这种额外的空间确实存在，它在不同位置发生变化，并且会影响密度。”

3. 新的“舞蹈规则”（研究结果）

通过考虑这种变化的“额外空间”并使用能量最小化原则，作者推导出了一套新的方程，这些方程主要实现了三件事：

A. 一种新的听声方式（声速）
论文表明，声速不仅仅是一个随机数字；它直接源于流体在被挤压时能量的变化。

• 隐喻： 把声音想象成人群中的涟漪。涟漪传播的速度取决于人们（分子）排列得有多紧密以及他们拥有多少能量。新的公式能够自然地计算出这种速度，而不需要预先告知其数值。它甚至表明，在气体中，声速大致等于气体分子跳动运动的平均速度。

B. 关于压力与速度的新规则（伯努利定律）
你可能听说过伯努利原理：“流体运动越快，压力越小。”

• 转折： 旧规则在水在管道中流动时表现得很好，但当出现巨大的密度跳变（如水撞击空气）时，它就会失效。作者创建了一个 广义伯努利定律 (Generalized Bernoulli's Law)。
• 隐喻： 想象一条河流流入瀑布。旧规则说能量保持不变。新规则则说：“等等，当水变成水雾（空气）时，由于水的‘稠密程度’在发生变化，一部分能量损失或转化了。”新方程解释了这种能量转移，使其即使在流体性质发生剧烈变化时依然准确。

C. 密度的“凸起”
这可能是最直观的结果。

• 旧观点： 如果观察水和空气之间的边界，旧模型认为密度会像一个向下的斜坡一样，从“重的水”平滑过渡到“轻的空气”。
• 新观点： 作者的数学模型预测了一个 凸起。当你穿过边界时，密度实际上会在降至空气水平之前先略微上升。
• 隐喻： 想象一群人（水）试图挤过一扇门进入一个空旷的大厅（空气）。当他们挤过门口时，在散开之前，可能会在瞬间挤压得更紧凑。新理论预测了这种“包装凸起”，这与先进的计算机模拟（称为密度泛函理论）所观察到的现象一致，而旧的简单模型却遗漏了这一点。

总结

这篇论文提出了一种新的方法，用于编写针对密度差异巨大的流体（如水和空气）的物理定律。

1. 它不再假装边界处的重量是恒定的。
2. 它承认分子占据的空间会发生变化（过剩体积）。
3. 它利用“最低能量”原则推导出了新的规则，解释了声音如何传播、压力如何变化，以及为什么水-空气边界处的密度实际上存在一个小“凸起”。

作者声称，这个新的框架适用于 任何 混合流体，无论它们的重量差异有多大，这为更精确地模拟诸如降雨、海浪破碎或气泡上升等现象打开了大门。

技术摘要

技术摘要：具有任意密度比的可压缩与不混溶流体：弥散界面理论

问题陈述
目前的计算流体力学（CFD）模型在模拟具有高密度比（如水-空气，比例 $\approx 1000$）的不混溶流体系统时面临困难。现有方法主要依赖于两种手段，两者均存在显著局限性：

1. Boussinesq 近似法： 该方法假设界面内密度（$\rho$）恒定，利用方程 $\rho(du/dt) = -\nabla p + \dots$。这种表述未能考虑到“真实”的动量演化 $d(\rho u)/dt$，忽略了密度变化及由此产生的动量传递。此外，它无法在不添加人为项的情况下自然地产生重力和浮力，因为这些力源于密度差异。
2. 插值法： 将密度定义为体积浓度的线性函数（$\rho = \rho_w \phi + \rho_a(1-\phi)$）。这种方法在不可压缩性（$d\rho/dt = 0$）与扩散性（$d\phi/dt \neq 0$）之间造成了根本性的冲突。此外，它假设整个区域内的过剩体积（通常为零）是均匀的，这与物理现实相悖，因为局部过剩体积取决于局部压力。因此，标准的连续性方程和线性插值无法捕捉到在水-空气界面处观察到的非单调密度剖面。

方法论
作者提出了一种基于能量最小化热力学原理的新型理论框架，这与传统的 Navier-Stokes 和 Euler 方程推导不同。该方法论包括：

• 体积与密度的重新表述： 该模型不再假设过剩体积为零，而是引入了界面内随空间变化的过剩体积（$v_e$）。这种体积在流体组分之间进行重新分配，导致偏密度（$\rho'_w, \rho'_a$）不再是常数，而是取决于局部压力和混合状态。这需要两个耦合的演化方程：一个用于体积浓度（$\phi$），另一个用于密度（$\rho$）。
• 耗散-守恒原理： 作者建立了一个将热力学第一定律（守恒）与第二定律（耗散）联系起来的定理。通过将总系统能量（$E$）定义为自由能（$F$）、压力能（$P$）、动能（$K$）和壁面自由能（$W$）之和，他们通过最小化能量耗散来导出动力学方程。
• 修正的 Gibbs-Duhem 关系式： 为了解决现有模型中空间梯度与泛函导数之间的数学不兼容问题，推导出了一个修正的等温 Gibbs-Duhem 关系式。该关系式明确考虑了与速度相关的熵和惯性功，从而导出了一个包含自由能和速度熵贡献的严格压力定义。
• 动能定义： 动能被表述为 $K = \int \rho u \cdot u \, d\Omega$（线性动量密度 $\zeta = \rho u$），而非传统的 $\frac{1}{2}\rho u^2$。这种区分允许将动量平衡（同时改变 $u$ 和 $\rho$）与速度耗散（仅改变 $u$）区分开来。

核心贡献
该推导产生了一个广义密度演化方程以及三项主要的理论进展：

1. 广义伯努利定律： 标准的伯努利原理（$p + \frac{1}{2}\rho u^2 = \text{const}$）被证明是一个特例。广义形式考虑了界面处的速度熵变化和密度演化：$p + \rho u^2 = \text{constant}$（在特定的稳态无粘条件下），承认由于密度差异，能量在界面处会发生跳跃。
2. 通用声速公式： 基于自由能对密度的导数推导出了一个通用的声速表达式：$c = \sqrt{d\Psi/d\rho}$。对于理想气体，这近似于分子的平均热运动速度。
3. 广义状态方程 (EOS)： 该理论表明，密度演化方程能够自然地生成流体的状态方程（例如，对于理想气体，$p \sim \rho \ln \rho$），而无需预设 EOS，这与标准的弥散界面方法形成对比。

结果与验证
论文通过三个案例研究验证了该理论：

• 声波传播： 推导出的声速方程与摄动分析一致，并能正确预测空气中声波的色散关系。
• 理想气体 EOS： 在稳态下，该模型重现了理想气体状态方程，证实了密度演化方程本身就捕捉到了热力学状态关系。
• 水-空气界面动力学： 通过管内移动平坦水-空气界面的解析解，证明该模型预测了非单调密度剖面。具体而言，密度从体相水值增加到界面中心的峰值，然后下降至体相气体值。这一结果与密度泛函理论（DFT）的发现 [11, 12] 相符，并反驳了线性插值法所预测的单调递减趋势。该模型也能正确重现弯曲界面的 Young-Laplace 方程。

意义
该论文声称提供了一个适用于任意高密度比的不混溶流体的通用框架，且同时适用于可压缩和不可压缩机制。通过将密度演化视为与过剩体积和粘性耗散耦合的动态变量，该模型解决了 Boussinesq 近似和线性插值法中的矛盾。这项工作表明，粘性耗散本质上会改变密度，挑战了在耗散流中严格的不可压缩假设。作者认为，该理论为 CFD 开辟了新的窗口，为流体-流体界面处（尤其是密度比极高时）的动量和能量传递提供了物理上一致的描述。