光风车的推力，一直算到边缘

克鲁克斯在1874年描述的克鲁克斯辐射计，是物理学家所说的辐射计效应的经典例子。一块两面温度不同的薄叶片放在稀薄气体中，会受到一个力，并以较热的一面落在后方的方式移动。作者给出的解释是：热面弹回的气体分子平均速度比冷面更快，因此气体对热面的反推更强。

这一效应在微型流体器件、气溶胶以及高层大气极稀薄空气中的飞行等领域重新受到关注。它的表现在很大程度上取决于气体的稀薄程度，用克努森数来衡量：即平均自由程（分子在两次碰撞之间平均行进的距离）与叶片尺寸之比。在非常稀薄的气体中，力分布在整个表面上；在较稠密的气体中，力集中在边缘附近。

最简单的叶片

京都大学信息学研究科的富田拓真（Takuma Tomita）、田口智（Satoshi Taguchi）和辻徹郎（Tetsuro Tsuji）着手描述这种力在整个稀薄程度范围内的表现，并特别关注边缘。他们的模型叶片是保留了边缘的最简单形状：一块无限薄的圆盘，可在无限大的气体中沿其轴线自由移动，一面比另一面稍热。

他们没有把气体当作流体，而是当作分子群体来处理，采用了玻尔兹曼方程的一个简化版本（BGK模型），并针对小温差和慢速运动作了线性化。难点在于锋利的边缘：它会在分子速度分布中造成突变，这些突变扩散到气体中，在三维情况下尤其棘手。团队采用了一种数值方法，沿直线路径追踪分子并显式处理这些突变，计算在京都大学的超级计算机上完成。

得益于线性，问题可以拆成几部分：一面受热的圆盘所受的辐射计力，以及运动圆盘所受的阻力。当二者平衡时，圆盘达到终端速度。

计算揭示了什么

（a）辐射计力随稀薄程度的变化（圆点），以及Selden等人在氩气中的测量值（三角形）；（b）圆盘的终端速度，有一个轻微的极大值（插图）。——图3，Tomita、Taguchi & Tsuji（2026），arXiv:2610.02896。

• 以无量纲形式表示时，力随稀薄程度稳定增长。在较稠密的气体中，它按**Kn^(3/2)**变化；在非常稀薄的气体中，它趋于一个常数值−π/2。
• 这解释了为什么通常改变压强的实验会发现实际力在中等压强下最大。根据Selden及其同事在氩气中的测量推得的数值，与计算结果画在一起。
• 终端速度并非稳定上升：它在中等稀薄程度时达到峰值，略高于在最稀薄气体中的值。让气体更稀薄并不总能使圆盘更快。
• 在圆盘周围，靠近热面的气体更热、密度更低，热面与相邻气体之间存在温度跳跃，沿热表面有一股流动，而在边缘周围的狭窄区域内集中着一股强烈的流动。

真正出力的是边缘

最引人注目的结果与两面之间的压力差有关。对于固定不动的圆盘，整个盘面都在推动，气体越稠密，边缘的推动就越强。而对于在稠密气体中以终端速度自由运动的圆盘，情形一分为二：边缘处的一条窄带推动圆盘前进，而盘面其余部分则在刹车，二者相互平衡。

这条边缘带的宽度与平均自由程成正比增长，因此在稠密气体中，它约为一个分子步长。作者指出，这与爱因斯坦1924年一篇关于辐射计力的论文中所作的假设一致，近期的估算仍在使用这一假设。

换一种方法验证

为了检验数值结果，作者利用了玻尔兹曼方程的一种与昂萨格倒易关系相关的对称性：辐射计力必须等于另一个问题——温度均匀的运动圆盘——中算出的热流。在每一种情况下，两者都至少吻合到三位有效数字——这也证明了这种对称性在存在锋利边缘时依然成立。

该模型没有考虑真实叶片的厚度，并采用了简化的碰撞定律。作者把他们的数据表视为带锋利边缘的动理学问题以及非球形气溶胶颗粒的参考数据，并计划下一步研究完整的玻尔兹曼方程和有限厚度的圆盘。