《Journal of Fluid Mechanics》：Analysis of thermally induced rarefied gas flow in a square cavity: from near-free-molecular to near-continuum conditions
编辑推荐：
中文翻译（摘要部分）
研究人员研究了一个封闭的二维方腔内由壁面温度差所诱导的气体流动。该问题在全部气体稀薄化率范围内得到了研究。基于已知的自由分子描述，研究人员推导了在高但有限努森数（Knudsen number, Kn）下的渐近解。后者基于Bhatnagar
中文翻译（摘要部分）
研究人员研究了一个封闭的二维方腔内由壁面温度差所诱导的气体流动。该问题在全部气体稀薄化率范围内得到了研究。基于已知的自由分子描述，研究人员推导了在高但有限努森数（Knudsen number, Kn）下的渐近解。后者基于Bhatnagar–Gross–Krook（BGK）、Shakhov和完整硬球（hard-sphere, HS）玻尔兹曼碰撞算子推导得出。这些结果辅以任意稀薄化率下的直接模拟蒙特卡洛（Direct Simulation Monte Carlo, DSMC）和粒子方案计算，以及在Kn趋于零时的热力学“纯传导”连续介质解。研究人员发现，气体在两个极限状态下均处于预期中的宏观静止状态：在无碰撞的自由分子极限中，热量以弹道方式输运；在零阶Navier–Stokes–Fourier（NSF）连续介质极限中，热量以纯传导方式输运。然而，在高但有限的Kn下，即使少量的分子碰撞也足以产生旋涡流动，其强度随接近自由分子极限而趋于零。对BGK渐近修正项的分析确定了该环流的微观起源和旋转方向。渐近结果与数值计算在直至O(1)量级的Kn范围内高度吻合。研究人员分析了气体总动能随气体稀薄化程度的变化，发现其在Kn约为O(10-1)量级时出现极大值。极大值的位置和大小均被量化。文中还讨论了HS、BGK和Shakhov模型预测之间的定量差异。
论文解读
一、研究背景与问题
气体在稀薄化条件下的流动行为与经典连续介质方程所预测的结果存在根本性差异。当流动的特征长度和时间尺度分别与分子平均自由程和平均碰撞时间相当时，Navier–Stokes–Fourier（NSF）描述所依赖的基本假设不再成立，非平衡效应变得显著。在此类条件下，分子间碰撞、分子与表面的相互作用以及弹道输运均发挥关键作用，产生连续介质理论无法捕捉的输运现象。因此，气体动理学理论在微纳系统、真空技术和低密度空气动力学环境的分析中具有核心重要性。
在众多非连续介质条件下出现的非平衡现象中，热致流动尤为引人关注。这类运动由温度不均匀性驱动，即使不存在外部施加的压力梯度或运动边界也可发生。根据气体稀薄化程度和热构型的不同，多种机制可能发挥作用，包括热蠕滑（thermal creep）、热应力（thermal stress）、辐射计力载荷（radiometric force loading）以及热边缘流（thermal edge flow）等。这些效应本质上源于动理学效应，在均匀压力和静止壁面条件下的连续介质极限中，NSF描述中并不存在。热致流动因其在努森泵（Knudsen pump）、微皮拉尼真空计（micro-Pirani gauge）及其他真空应用中的基础重要性和实际意义而受到广泛关注。
在热致流动的框架内，具有非均匀加热壁面的封闭方腔构型——常见于真空封装器件中——已被广泛研究。在此类封闭空间中，边界上的温度差可产生内部环流运动。从理论角度看，该构型为研究弹道分子输运、分子间碰撞与壁面边界条件之间的相互作用提供了一个典型的模型问题；从实际角度看，理解此类热致腔体流动与各种微系统和低压环境中受限稀薄气体的热管理和气体动力学行为密切相关。
然而，现有研究存在三个主要空白。第一，除Cai（2008）的自由分子分析和Sone等人的一般自由分子解外，现有处理方法几乎完全依赖计算密集的数值模拟（主要为DSMC方法），对诱导运动的动理学机制提供的洞察有限。第二，现有解析工作集中于滑移和早期过渡区域，其中热蠕滑和热应力占主导，而相反的近自由分子极限相对未被探索。第三，达到高稀薄化程度的分析大多局限于线性化动理学模型或小壁面温差情形，强非线性、有限温差区域仍未得到探索。因此，关于弱分子间相互作用如何修正无碰撞状态并产生有限Kn环流的解析描述仍然缺失。
二、研究内容与结论
本研究直接针对上述空白展开。研究采用Mousivand和Roohi（2022b）的双热壁—双冷壁方腔构型，推导了完整非线性硬球（HS）玻尔兹曼碰撞算子以及BGK和Shakhov模型的近自由分子渐近解。与早期高Kn分析不同，这些解适用于任意冷壁与热壁温度比，消除了先前线性化处理中的线性化和小温差限制。分析确定了近弹道碰撞机制，即即使少量的分子间相互作用也能产生热致环流。具体而言，壁面发射的分子群体携带不等的轨迹加权动量，因此第一阶碰撞修正项驱动了一个净的从冷壁到热壁的滑移，其幅度随Kn减小而增大。该工作补充了Sone等人已知的自由分子解，描述了远离无碰撞极限时产生热致运动的有限Kn修正。结合覆盖所有Kn范围的DSMC和基于粒子的解，该研究提供了全稀薄化率下热致流动的全面描述，包括流动强度极大值位置和大小。在连续介质极限中，气体热力学状态与温度的“纯传导”解进行了比较。HS、BGK和Shakhov描述之间的定量差异也得到了检验。
三、主要技术方法
研究人员采用了多种互补的研究方法。首先，建立了基于气体动理学理论的数学模型，采用圆柱速度空间坐标系描述二维稳态方腔内的气体状态，并应用漫反射Maxwell边界条件模拟气体与壁面的相互作用。其次，在高Kn极限下采用逆Kn幂级数展开方法，推导了自由分子极限解及其一阶修正的渐近解析解。第三，在连续介质极限下获得了“纯传导”解析解。第四，采用直接模拟蒙特卡洛（DSMC）方法在任意稀薄化率下进行数值模拟。第五，采用粒子方案（particle scheme）进行数值计算，与DSMC结果相互验证。研究涵盖了HS、BGK和Shakhov三种碰撞模型，以比较不同动理学模型预测的差异。
四、研究结果
自由分子极限。 在无碰撞极限下，分布函数沿自由飞行轨迹保持恒定，每个粒子携带由其最后一次壁面碰撞确定的分布函数。通过应用壁面不可穿透条件，研究人员确定了各壁面的反射粒子通量，进而得到自由分子极限下的宏观气体状态。结果表明，在无碰撞自由分子极限中，气体处于宏观静止状态，热量以弹道方式输运。
高Kn渐近解。 在自由分子解的基础上，研究人员推导了O(Kn-1)阶的碰撞修正项。分析表明，壁面发射的分子群体携带不等的轨迹加权动量，因此即使极少的分子碰撞也足以产生净的冷壁到热壁滑移，从而驱动方腔内的环流。BGK渐近修正项的分析揭示了该环流的微观起源和旋转方向。渐近结果与数值计算在直至O(1)量级的Kn范围内高度吻合。
连续介质极限。 在Kn趋于零的极限下，气体达到NSF连续介质描述，热量以纯传导方式输运，气体同样处于宏观静止状态。研究人员获得了“纯传导”温度解，并与数值结果进行了比较。
全Kn范围的数值结果。 通过DSMC和粒子方案计算，研究人员获得了任意稀薄化率下的流场。结果表明，气体总动能随稀薄化程度变化，在Kn约为O(10-1)量级时出现极大值，极大值的位置和大小均被量化。
模型比较。 研究人员比较了HS、BGK和Shakhov三种碰撞模型的预测结果。Shakhov模型因能给出单原子气体正确的连续介质极限Prandtl数2/3（而BGK模型给出1），被认为在耦合热与动力学效应的问题中更为可靠。三种模型预测之间的定量差异在文中得到了讨论。
五、讨论与结论
讨论部分指出，本研究填补了热致腔体流动研究中的三个空白：提供了超越纯数值模拟的解析描述；将解析工作从滑移和早期过渡区域扩展至近自由分子极限；将分析从线性化模型和小温差限制推广至完全非线性、任意温度比的情形。研究揭示了近弹道碰撞机制是有限Kn下热致环流产生的微观根源，即壁面发射分子群体间不等的轨迹加权动量在首次碰撞修正中驱动了冷到热的净滑移。
结论部分总结如下：气体在无碰撞自由分子极限和零阶NSF连续介质极限中均处于宏观静止状态，分别对应弹道热输运和纯传导热输运；在高但有限的Kn下，即使少量分子碰撞也足以产生旋涡流动，其强度随接近自由分子极限而趋于零；BGK渐近修正项的分析确定了该环流的微观起源和旋转方向；渐近结果与数值计算在直至O(1)量级的Kn范围内高度吻合；气体总动能随稀薄化程度变化，在Kn约为O(10-1)量级时出现极大值，极大值的位置和大小均被量化；HS、BGK和Shakhov模型预测之间存在定量差异，其中Shakhov模型在耦合热与动力学效应的问题中更为可靠。
