基于特异速度的玻尔兹曼格式能否准确分辨接触间断？

原始PVU格式为何抹平静止接触

无粘气体动力学的动理学或玻尔兹曼格式之所以有吸引力，是因为当速度分布为麦克斯韦分布时，欧拉方程可作为玻尔兹曼方程的矩出现；而且分子速度玻尔兹曼方程中的对流项是线性的，因此迎风处理更简单，格式也不依赖宏观方程组的特征结构。传统动理学格式具有动态一致性、正性（L1稳定）、熵稳定，并在高马赫数下稳健，但它们在1D中产生过度扩散的接触间断、在2D中产生过度扩散的滑移面，人为加宽层流边界层，并在孤立接触间断附近出现虚假的压力和速度振荡。

论文从特异速度迎风（PVU）玻尔兹曼格式出发，该格式将流体速度与特异速度分离，并分别对二者施加迎风处理，在宏观层面得到Zha-Bilgen型的对流-压力分裂。原始PVU格式稳健但过度扩散，无法识别任何静止间断。作者对半离散对流步的修正偏微分方程（MPDE）分析揭示了原因：与特异速度部分相关的数值扩散系数在马赫数趋于零时不消失，这正是该格式扩散静止接触间断的主要原因。

重新设计特异速度扩散使其在接触处关闭

用特异速度重写的玻尔兹曼方程允许将流体速度通量与特异速度通量独立处理。在玻尔兹曼层面，对流通量采用迎风、压力通量采用中心处理，本可精确分辨孤立的静止接触间断，但这种离散在强非线性间断处不稳定。因此论文总体上保留稳定的完全迎风处理，转而在数值通量中引入扩散控制参数，并要求它们满足动态一致性条件，使得两侧守恒状态相等时参数也相等。

目标是让特异速度数值扩散在马赫数趋于零时消失，在其他情形下保持低但非零，从而锐化线性间断，同时不损失跨激波和膨胀扇的稳健性。变体1将对流控制设为1，将特异速度控制设为局部声速与界面平均声速之比。变体2使用速度比的对流控制，并将特异速度控制设为exp(-c^2/u^2)，当整体速度为零时该值为零，因此数值通量再现了通过孤立静止接触间断连接左右状态的精确黎曼解。

变体3消除接触处的残余能量通量

在静止接触间断处，变体1仍会产生一个非零的数值能量通量项，其正比于p a^2之差除以两倍平均声速，正是这一项扩散了1D静止接触。变体3将该非零项中的密度近似为p除以平均声速的平方，使该项变为正比于(p_j - p_{j+1})乘以平均声速。只要压力相等，该项就自动关闭，因此格式能够识别一维静止接触间断。

同样的思路扩展到二维，将数值能量通量中的密度近似为p_K除以平均声速的平方。论文指出，这种二维扩散在静止滑移面处不会完全关闭，因此变体3对滑移面的分辨不如变体2。所有变体的稳定性条件来自对紧支撑PVU格式的正性分析，呈CFL型形式，涉及法向速度加三倍RT的平方根，并在单元面上求和。

基准算例中精确捕捉1D接触与2D滑移面

论文报告了若干1D和2D可压缩流基准算例。在水平滑移流算例中，马赫3.0气流滑过马赫2.0气流，计算域为单位正方形，离散为20×20笛卡尔网格，全场密度和压力相同，左右分别为超声速入流和出流。变体2精确再现了网格对齐、静止、孤立的滑移面，且无数值扩散。作者用该算例说明一个对比较格式的读者很重要的观点：能准确捕捉一维静止接触间断的格式，未必能以同样精度捕捉二维滑移面。

奇偶解耦算例让马赫6正激波通过划分为800×20方形单元的矩形通道，中心线受到扰动，演化至55个时间单位。两种变体都未出现这种激波不稳定性，论文称该不稳定性在HLLC等格式中会导致解完全崩溃。正性算例使用大初始马赫数，Mx = 10、My = -3或2，压力0.1 bar，温度35 K，单元尺寸0.02×0.02，CFL为0.25，迭代100次；两种变体都保持压力和密度为正且不崩溃。

激波算例显示稳健性得以保留

在激波绕射算例中，马赫5.09正激波绕过90度角，计算域为单位正方形，划分为200×200笛卡尔网格，角点位于x = 0.05、y = 0.45，解演化至0.1561个时间单位。数值结果没有违反熵的解和振荡，而一阶HLLC解在上边界附近的激波前方出现不稳定性，且在高阶结果中更明显。

双马赫反射算例让与底壁成60度的斜马赫10激波通过四单位长、一单位宽、划分为480×160单元的域，演化至0.2个时间单位。两种变体都未在马赫杆中产生非物理扭结，而这是HLLC格式中已知的不稳定性。在1920×640网格上使用Barth-Jespersen限制器和原始变量重构重复该算例，变体2干净地捕捉了主滑移线上的复杂涡结构。在马赫20、160×20网格的超声速半圆柱算例中，格式也避免了碳核不稳定性。

该格式尚未解决的问题

证据是数值和基准算例层面的。论文通过算例展示了1D静止接触间断和2D滑移面的精确捕捉，并证明了紧支撑PVU格式在CFL型条件下的正性，但基于麦克斯韦分布的PVU格式并不继承这一正性证明，因为麦克斯韦分布在速度空间具有无限支撑。作者为稳定性条件采用3 sqrt(RT)速度截断，指出对麦克斯韦分布一阶矩在该区间积分可捕获其总质量的约99.73%，而在Mandal和Deshpande的数值实验中，第二至第五阶矩分别覆盖总面积的约99.73、97.92、97.43和93.2%。

论文还声明该研究未使用任何数据，且不存在利益冲突。变体以一阶和二阶格式呈现，二阶扩展在时间上使用三阶段强稳定性保持Runge-Kutta。报告的比较对象是原始PVU格式，并在若干激波算例中与论文描述的HLLC行为进行比较。

何时优先选择哪个变体

论文自身的优先级从其比较中清晰可见。变体2是在水平滑移流算例中精确再现网格对齐静止滑移面的那个，因为其特异速度扩散在整体速度为零时消失。变体3是消除导致1D静止接触扩散的残余数值能量通量的那个，但其二维扩散在静止滑移面处不会完全关闭，因此对滑移面的分辨不如变体2。作者得出结论：能精确捕捉稳态二维滑移面的格式，在激波-边界层干扰算例中更有效地捕捉负表面摩擦系数。

对于决定从这篇论文中获取什么的读者，实际的分野在于静止线性间断处的精确性与强非线性间断处的行为之间。变体2以整体速度为零为代价换取精确的滑移面捕捉；变体3通过在压力相等时关闭能量通量项换取更干净的1D接触。两者在奇偶解耦、正性、激波绕射、双马赫反射和超声速半圆柱算例中都保留了原格式的稳健性，因此选择取决于应用必须最锐利分辨哪种间断，而非一般性的精度排名。

本文引用的文献