Multispecies Bhatnagar-Gross-Krook models and the Onsager reciprocal relations

本文表明，大多数现有的多组分系统速度无关的Bhatnagar-Gross-Krook模型未能满足昂萨格倒易关系，从而使得将其输运性质校准以匹配特定流体变得复杂。

要点

想象一下，你试图预测两种不同颜色的弹珠（比如红色和蓝色）在摇晃盒子时如何混合。科学家拥有一套描述这些弹珠行为的“规则手册”，称为玻尔兹曼方程。它极其精确，但也极其复杂，求解它就像在飓风肆虐时试图数清海滩上的每一粒沙子。

为了简化问题，科学家们创建了这个规则手册的简化版本，称为BGK 模型（以 Bhatnagar、Gross 和 Krook 的名字命名）。你可以将 BGK 模型视为一张“作弊纸”或“捷径”，它在不进行繁重数学计算的情况下，近似模拟弹珠的复杂行为。几十年来，它一直被用于模拟各种现象，从机翼上的气流到恒星中的等离子体。

问题所在：指南针失灵

E. S. Benilov 撰写的这篇论文指出，在处理不同气体的混合物（如氧气和氮气，或水蒸气与空气）时，这种“作弊纸”最常见的版本存在一个重大缺陷。

作者发现，这些流行的 BGK 模型违反了一条基本的物理定律，即昂萨格倒易关系。

以下是昂萨格关系的一个简单类比：
想象一条连接两个城镇的双向街道，交通在两地之间流动。

• 规则 A： 如果你在城镇 A 修了一座山丘，它会影响城镇 B 的车流速度。
• 规则 B： 如果你在城镇 B 修了一座山丘，它会影响城镇 A 的车流速度。

昂萨格关系指出，这两种效应必须完美地“相互协调”。如果城镇 A 的山丘使城镇 B 的车流速度减慢了 10%，那么城镇 B 的山丘必须使城镇 A 的车流速度减慢一个数学上相关联的量。这是一条对称规则；宇宙要求这些相互作用达到平衡。

论文发现了什么

Benilov 用这条规则检验了标准的 BGK“作弊纸”。他发现：

1. 模型是不对称的： 在 BGK 模型中，城镇 A 的“山丘”（温度变化）对城镇 B 的“交通”（质量流）完全没有影响。然而，城镇 B 的“山丘”（密度变化）确实会影响城镇 A 的“交通”（热流）。
2. 不匹配： 由于方程的一侧为零而另一侧不为零，对称性被破坏了。这个模型就像是一个永远向一侧倾斜的天平。
3. 后果： 由于该模型破坏了这一基本规则，因此不可能对其进行“校准”。校准就像调谐收音机以获得清晰的信号。如果你试图调整 BGK 模型的旋钮（参数），使其与特定流体的真实世界数据相匹配，这是行不通的。该模型在根本上存在缺陷，无论你怎么调整，它都无法达到完美的准确性。

“水蒸气”例外（以及为何它无法力挽狂澜）

你可能会想：“好吧，也许这仅对奇怪的气体很重要。那像水蒸气和空气这样的常见混合物呢？”

论文也对此进行了检查。即使温度对质量流的影响微乎其微（对于水蒸气和空气确实如此），该模型仍然失效。要使该模型在这种特定情况下生效，你必须将某个旋钮调至无穷大，这实际上通过使所有运动停止而彻底破坏了模型。因此，该模型在复杂混合物和简单混合物中均告失败。

有没有好的模型？

论文指出，确实存在少数其他更复杂的 BGK 模型遵循这些规则，但它们也有自己的问题（例如违反其他物理定律，如确保熵总是增加的"H 定理”）。

作者总结道，目前没有任何现有的 BGK 模型是完美的。一个完美的模型需要：

• 守恒质量、动量和能量。
• 遵循热力学定律（熵）。
• 公平地对待相同粒子。
• 保持温度和浓度为正值。
• 允许科学家将其调谐以匹配任何真实流体。
• 并且遵守昂萨格倒易关系（对称规则）。

目前，我们拥有的每一个模型都至少未能通过其中一项测试。

核心结论

这篇论文是对使用这些模型的科学家们的一个警告。如果你使用标准的 BGK 模型来模拟气体混合物，你正在使用一个在根本上与物理定律“失谐”的工具。它可能给你一个关于正在发生什么的粗略概念，但你不能信任它能为真实世界的流体提供精确的、经过校准的结果，因为它违反了自然界的一个核心对称规则。作者希望未来有人能构建一个“完美”的模型，解决所有这些问题。

技术摘要

技术摘要：多组分 Bhatnagar–Gross–Krook 模型与昂萨格倒易关系

问题陈述
Bhatnagar–Gross–Krook (BGK) 模型是玻尔兹曼动力学方程的一种广泛使用的唯象近似，尤其适用于气体混合物。虽然原始的单组分模型及其向混合物的推广（例如 Gross 和 Krook [2]）因其定性正确性和适应性而备受推崇，但它们尚未接受关于是否符合昂萨格倒易关系的严格检验。这些关系对输运系数施加了基本的热力学约束，具体将质量通量中的温度梯度系数（热扩散或索雷效应）与热通量中的密度梯度系数（杜福尔效应）联系起来。从第一性原理推导出的模型会自动满足这些关系，而唯象模型则未必。模型若不符合这些关系，不仅会引发对其物理相关性的质疑，而且关键的是，会阻碍将模型参数校准以匹配特定真实流体的输运特性。

方法论
本文构建了适用于单原子气体二元混合物的最一般多组分 BGK 模型。该公式化包含了碰撞频率 ($\nu_{ij}$) 以及参数 ($\alpha_i, \beta_i, \gamma_i$)，这些参数定义了向具有组分依赖性速度和温度的混合麦克斯韦分布 ($M_{ij}$) 的弛豫过程。作者通过对这些参数施加特定约束，确保模型守恒质量、动量和能量。

为了分析输运性质，本文采用了扩散近似而非标准的流体动力学近似。这涉及以下步骤：

1. 标度化：引入小参数 $\epsilon$ 以拉伸空间坐标 ($\nabla \to \epsilon \nabla$) 和时间 ($\partial_t \to \epsilon^2 \partial_t$)，这适用于缓慢、微弱的扩散流动。
2. 展开：将分布函数、宏观速度、温度和密度按 $\epsilon$ 的幂次进行展开。
3. 推导：将这些展开式代入重标度的 BGK 方程，以求解分布函数的一阶修正项 ($f^{(1)}_i$)。
4. 通量计算：对 $f^{(1)}_i$ 进行积分，推导出用密度梯度和温度梯度表示的质量通量 ($J_i$) 和热通量 ($Q$) 的显式表达式。
5. 比较：将推导出的 BGK 通量与满足昂萨格关系的标准流体动力学表达式（恩斯科格–查普曼理论）进行比较。

主要贡献与结果
本文的核心结果是证明最常见的多组分 BGK 模型版本——其中碰撞频率和模型系数与分子速度无关——不满足昂萨格倒易关系。

具体而言，分析揭示了一个根本性的不一致：

• 根据昂萨格关系，如果连接密度梯度与热通量的系数（杜福尔效应，$C$）非零，那么连接温度梯度与质量通量的系数（热扩散率，$B$）必须非零，且两者通过 $C = pB$ 相关联。
• BGK 模型的推导得出的质量通量中，热扩散项严格为零 ($B=0$)，除非取某些特定的、非物理的极限。
• 然而，由此产生的热通量中包含一项与密度梯度成正比（即杜福尔效应）的项，该项非零且与组分的质量差成正比。
• 因此，BGK 模型预测了非零的杜福尔效应，同时索雷效应为零，从而违反了昂萨格关系 $C = pB$。

文章进一步表明，这种差异无法通过校准模型参数（例如控制相对速度弛豫的参数 $\beta$）来解决。不存在任何 $\beta$ 值能够同时满足一般混合物的标准输运方程和昂萨格关系。即使在热扩散率可忽略的流体（即 $B \approx 0$）这一特定情况下，BGK 模型也无法在不迫使整个扩散通量消失的情况下重现正确的热通量行为。

意义与主张
本文主张，这种不合规性给多组分 BGK 模型的实际应用带来了重大问题。由于从该模型推导出的输运系数不遵守必要的热力学倒易性，因此不可能“校准”该模型——即选择参数值以使模型的输运特性与特定真实流体（如水蒸气和空气）相匹配。

作者指出，虽然该模型未能精确满足昂萨格关系，但在两个特定极限下，它能渐近地满足这些关系：

1. 当质量比 $m_1/m_2 \to 0$ 时（适用于电子质量相对于离子可忽略不计的电离等离子体）。
2. 当质量比 $m_1/m_2 \to 1$ 时（此时热通量表达式中的第一项消失）。

文章最后列出了“完美”多组分动力学模型的要求，包括守恒定律、H 定理、不可区分性原理、温度/浓度的正定性、能够表示任意普朗特数/施密特数，以及符合昂萨格关系。作者指出，迄今为止，尚无现有的 BGK 型模型被证明能同时满足所有这些要求。研究结果表明，尽管当前的 BGK 模型是有用的近似，但它们无法满足昂萨格关系，这限制了其在一般流体动力学应用中精确校准时的物理保真度和实用性。