A unified gas-kinetic wave-particle method for multiscale binary-species gas mixtures

本文提出了一种统一气体动理学波粒（UGKWP）方法，用于模拟多尺度双组分气体混合物，该方法通过整合修正的平衡模型、基于 Shakhov 的普朗特数修正以及改进的粒子输运机制，准确捕捉了从连续流到稀薄流各流态下各组分特有的速度差与温度差，同时展示了与高超声速流动的直接模拟蒙特卡洛（DSMC）结果的高度一致性。

要点

想象一下，你试图预测两种不同气体（比如氩气和氖气）在混合并以极高速度运动时的行为，就像航天器重返大气层时掠过其表面的空气那样。

这是一个棘手的问题，因为气体的行为取决于其拥挤程度。在拥挤的房间（高密度）中，气体表现得像一种平滑的流体，如同河流中流动的水。在稀疏的房间（低密度）中，气体则表现得像个体随机相互碰撞的人，如同人群在空旷的大公园中穿行。

大多数计算机程序难以同时处理这两种情况。它们通常必须做出选择：要么模拟平滑流动（这在空旷空间中会失效），要么模拟单个粒子（这在拥挤区域中过于缓慢且昂贵）。

解决方案：“波 - 粒”混合方法

本文介绍了一种名为**统一气体动理学波 - 粒（UGKWP）**的新方法。可以将这种方法想象成一个智能交通控制器，能够瞬间在两种观察气体的方式之间切换：

1. 波视角（群体）：当气体密度高时，该方法将其视为平滑、连续的波。它不追踪每一个分子，而是计算“平均”行为，就像预测河流的流动一样。这既快速又高效。
2. 粒子视角（个体）：当气体稀疏或运动极快（例如在激波附近）时，该方法切换为追踪单个粒子。它将粒子模拟为在周围弹跳的微小台球。这捕捉到了波所忽略的混乱、非平滑行为。

这种新方法的奇妙之处在于，它并非简单地来回切换，而是同时进行两者。它能够自动决定有多少气体表现得像波，有多少表现得像粒子，直至最细微的层面。

“双组分”挑战

本文的具体突破在于处理两种不同气体混合在一起（双组分混合物）的情况。

想象一个舞池，里面有两大群舞者：重舞者（氩气）和轻舞者（氖气）。

• 问题：当它们混合时，轻舞者可能比重舞者移动得更快。它们也可能具有不同的温度。标准方法通常将它们视为完全相同，或者在它们如何交换能量和动量方面感到困惑。
• 修正：作者建立了一套新的“规则手册”（数学模型），用于描述这两组人如何相互作用。他们精确计算出了这两组人应该达到的“目标”状态。
  • 他们修正了“摩擦”（粘度），使重舞者和轻舞者不会不切实际地相互滑过。
  • 他们修正了“热传递”（普朗特数），使冷热区域能够正确混合。
  • 他们甚至改进了对“最快舞者”（高速粒子）的处理方式，意识到高速粒子比低速粒子碰撞更频繁，从而改变了它们的运动方式。

他们的测试

为了证明其方法有效，他们进行了多项模拟：

1. 激波：他们模拟了一堵气体墙撞击另一股气体（类似于音爆）。他们的方法比旧方法更准确地预测了温度和密度的变化，特别是在撞击前极快运动的气体区域。
2. 气体混合：他们观察了氩气和氖气在管道中的混合。他们的方法正确预测了两种气体如何分离和移动，即使在气体非常稀薄的情况下，其结果也与“金标准”模拟方法（DSMC）的结果相符。
3. 滑动平板：他们模拟了两块移动平板之间的气体（库埃特流动）。他们的方法捕捉到了气体在边缘处的滑移现象，这是一个很难准确获取的细节。
4. 高超声速圆柱体：最后，他们模拟了气体以超音速绕圆柱体飞行的情况。其表面上的压力、摩擦和热量结果与金标准粒子模拟的结果几乎完美吻合。

结论

本文提出了一种模拟气体混合物的新且更智能的方法。它将流体方程的速度与粒子追踪的准确性相结合。通过专门修正两种不同气体相互作用的数学模型，它提供了一种可靠的工具，用于理解复杂的流动，特别是那些涉及高速航空航天飞行器的流动，其中不同气体混合、升温并以极端方式表现。

技术摘要

技术摘要：面向多尺度二元气体混合物的统一气体动理学波粒方法

问题陈述
模拟多尺度气体混合物，特别是在高超声速条件下，由于热非平衡、化学反应与复杂组分间相互作用的紧密耦合，带来了巨大挑战。现有数值方法往往难以在连续流与稀薄流区间之间高效衔接，同时准确捕捉扩散、粘性及热传导等组分特异性现象。具体而言，标准动理学模型在气体混合物的连续流极限下，常无法恢复正确的输运系数（粘度、扩散及热传导）；而直接模拟蒙特卡洛（DSMC）等粒子方法在处理近连续流时计算成本高昂。此外，在高速流动中捕捉不同组分的独特速度与温度分布，需要能够考虑跨组分碰撞与非平衡效应的稳健物理模型。

方法论
作者提出了一种专为多尺度二元气体混合物设计的统一气体动理学波粒（UGKWP）方法。该方法在有限体积框架下运行，采用直接建模策略，将气体分布函数分解为解析流体动力学波（近平衡态）与离散粒子（非平衡态）。

关键方法论组成部分包括：

1. 动理学模型：本研究采用基于 Andries-Aoki-Perthame (AAP) 框架的单弛豫 Bhatnagar–Gross–Krook (BGK) 型模型。为恢复正确的输运系数，目标平衡分布函数 ($g_\alpha$) 依据 Groppi 等人提出的模型构建。这涉及引入第二个弛豫参数，以同时恢复剪切粘度和扩散系数。
2. 普朗特数修正：为解决热传导系数问题，将 Shakhov 模型扩展至气体混合物。该修正调整分布函数内的普朗特数 ($Pr$)，确保在连续流区域正确恢复热通量。
3. 通量与源项公式：
   • 波演化：近平衡部分被视为确定性流体动力学波。利用沿特征线的积分解推导通量，该通量通过特征积分解内在地计及扩散效应，而非仅通过算子分裂将扩散处理为源项。
   • 粒子演化：离散粒子捕捉非平衡自由输运部分。针对高速粒子引入物理修正的碰撞时间 ($\tau^*$)。这一从单组分工作扩展至混合物的修正，增加了相对于组分漂移速度具有高热速度的粒子的碰撞频率，从而提高了激波前温度分布和驻点热通量的计算精度。
   • 一致性：该方法通过在通量和源项计算中使用相同的目标宏观变量 ($\tilde{W}^*_\alpha$) 和时间层级系数，确保了波描述与粒子描述之间的严格一致性。这消除了对动量和能量交换源项隐式求解器的需求。
4. 算法：算法流程包括初始状态采样、自由输运（将粒子分为无碰撞和碰撞两类）、碰撞（更新宏观变量并采样新粒子）以及迭代。

主要贡献

• 混合物的统一框架：本文建立了一个统一的气体动理学框架，实现了二元气体混合物在纳维 - 斯托克斯求解器（连续流）与玻尔兹曼求解器（稀薄流）之间的无缝过渡。
• 准确的输运系数：通过整合 Groppi 模型和扩展的 Shakhov 修正，该方法在连续流极限下成功恢复了正确的粘度、扩散和热传导系数。
• 改进的高速粒子模型：将物理修正的碰撞时间 ($\tau^*$) 扩展至气体混合物，解决了高马赫数流动建模中的局限性，特别是提高了激波前区域温度分布的预测精度。
• 一致的波粒耦合：该公式实现了波描述与粒子描述之间的一致性，无需额外的源项隐式模块，从而促进了向多组分系统的扩展。

结果
所提出的方法通过一系列涵盖从连续流到稀薄流区域的数值测试进行了验证：

• 激波结构：在 $Ma=1.5$和$Ma=3.0$ 条件下对二元气体混合物激波结构的模拟表明，UGKWP 方法准确捕捉了密度和温度分布。与基于标准 AAP 的 UGKS 相比，该结果与参考玻尔兹曼解的吻合度更高，特别是在激后温度分布（得益于普朗特数修正）和激前区域（得益于 $\tau^*$ 修正）。
• 质量扩散：在氩 - 氖混合物中，该方法在连续流区域 ($Kn=0.0061$) 成功恢复了纳维 - 斯托克斯扩散系数，并在一系列克努森数 ($Kn=0.061, 0.61$) 下与 DSMC 结果表现出良好的一致性。
• 库埃特流动：对四种稀薄参数 ($\hat{\delta}_0 = 100$ 至 $0.1$) 下的 Ar-Ne 混合物库埃特流动的模拟，展示了该方法捕捉速度滑移和组分分离的能力，其结果介于离散 UGKS 与离散速度法 (DVM) 基准结果之间。
• 高超声速流动：针对马赫数 ($Ma=3, 9, 18$) 和克努森数 ($Kn=0.01$至$1$) 各异的圆柱绕高超声速流动，该方法准确预测了壁面压力、剪切应力和热通量系数。结果与 DSMC 模拟表现出极好的一致性，捕捉到了驻点线上速度与温度的组分特异性差异。

意义与主张
本文主张，所提出的 UGKWP 方法为模拟涉及气体混合物的复杂多尺度流动奠定了坚实基础。通过准确捕捉组分间的速度与温度差异并恢复正确的输运系数，该方法满足了航空航天工程领域，特别是临近空间飞行器对先进数值模拟日益增长的需求。作者指出，该方法非常适合未来扩展至更复杂的物理模型，包括热非平衡、化学反应和等离子体输运，且未编造具体的未经验证的应用案例。这项工作证明，统一的方法在保持计算效率和物理保真度的同时，能够有效处理高超声速流动的多尺度特性。