› 克努森数模拟器 返回
稀薄气体动力学模拟器

克努森数模拟器 — 稀薄气体的流动区域

由温度 T、压力 P、分子直径 d 和特征长度 L 实时计算平均自由程 λ 与克努森数 Kn = λ/L，自动判定连续流、滑移流、过渡流、自由分子流四个区域。尺度比较与 Kn 区域图双画面，让你直观掌握流体模型选择的物理。

参数设置
温度 T
K
log₁₀(压力 P)
log Pa
P = 1000 Pa
分子直径 d
pm
特征长度 L
μm

默认值：T = 300 K、P = 1000 Pa（中真空）、d = 370 pm（相当于 N₂）、L = 10 μm（MEMS 尺度）。使用 Boltzmann 常数 k = 1.380649×10⁻²³ J/K。

计算结果
—
平均自由程 λ
—
克努森数 Kn
—
流动区域
—
连续介质极限 L (Kn=0.01)
尺度比较（特征长度 L 与 平均自由程 λ）

左侧蓝色条＝特征长度 L／右侧橙色条＝平均自由程 λ（同一对数尺度）／背景点＝分子（数密度 ∝ P/T）／背景色＝流动区域（绿=连续流，黄=滑移流，橙=过渡流，红=自由分子流）。

Kn 区域图（流动区域判定）

横轴＝log₁₀(Kn)（−3〜2）／四条彩带＝连续流、滑移流、过渡流、自由分子流／黄圆＝当前 Kn 位置／边界线为 Kn = 0.01、0.1、10。

理论与主要公式

克努森数是平均自由程与装置尺寸之比，定量化稀薄气体流动的"稀薄程度"。通过硬球气体的平均自由程由下式计算。

克努森数的定义：

$$\mathrm{Kn} = \frac{\lambda}{L}$$

平均自由程（硬球模型，含相对速度修正 √2）：

$$\lambda = \frac{k\,T}{\sqrt{2}\,\pi\,d^2\,P}$$

流动区域边界：

$$\mathrm{Kn} < 0.01 \,\,(\text{连续流}),\quad 0.01 \le \mathrm{Kn} < 0.1 \,\,(\text{滑移流})$$ $$0.1 \le \mathrm{Kn} < 10 \,\,(\text{过渡流}),\quad \mathrm{Kn} \ge 10 \,\,(\text{自由分子流})$$

$k = 1.380649\times10^{-23}$ J/K 为 Boltzmann 常数，$T$ 为温度 [K]，$d$ 为分子直径 [m]，$P$ 为压力 [Pa]，$L$ 为装置特征长度 [m]。$\lambda$ 是分子到下次碰撞前走过的平均距离，$\mathrm{Kn}$ 是从装置看气体稀薄程度的指标。

什么是克努森数模拟器

🙋
流体力学教科书说"Navier–Stokes 方程要在连续介质假设成立时才能用"，但这个界线到底在哪？感觉对空气来说什么情况都能用。
🎓
决定它的就是克努森数 Kn = λ/L。λ 是气体分子的平均自由程，L 是装置的特征长度。Kn < 0.01 可以放心用教科书上的连续介质 + 无滑移边界条件；0.01〜0.1 要加滑移边界延寿；0.1〜10 是过渡流，连续介质假设崩坏；10 以上是完全分子流，必须用 Boltzmann 方程或 DSMC。常压空气在 1 m 管道里 Kn ≈ 10⁻¹⁰，完全连续介质；但 MEMS 的 10 μm 通道、中真空 1000 Pa 下 Kn ≈ 0.7，就跑到过渡流去了。
🙋
本工具的默认值（T = 300 K、P = 1000 Pa、d = 370 pm、L = 10 μm）正好是"过渡流"，这是 MEMS 中真空的典型例子吗？
🎓
是的。λ = kT/(√2 π d² P) ≈ 6.81 μm，相对 L = 10 μm 得到 Kn = 0.68。MEMS 微加速度传感器或半导体真空腔内小缝隙的典型情形。这一区域需要在 Navier–Stokes 上加"滑移边界 + 温度跳跃"的扩展模型（DSMC 或 Burnett 方程），普通 CFD 会有较大误差。本工具里滑动 L 滑块，能在 Kn 区域图上看到流动区域由连续流（绿）→ 滑移流（黄）→ 过渡流（橙）→ 自由分子流（红）的变化。
🙋
半导体溅射腔会有意降到低压，是为了提高 Kn 吗？
🎓
正是。溅射或 MBE 追求 Kn ≫ 1 的完全自由分子流。在本工具把 P 设为 0.1 Pa（log P = −1）、d = 250 pm（Ar）、L = 30 cm = 3×10⁵ μm 试试，λ 仍只有几十 mm、Kn ≈ 0.1，没那么高。所以实际设备会降到 0.001〜0.01 Pa（UHV）让 λ 达到米级，保证靶到基板的直线飞行。但为了维持等离子体密度又需要一定压力，最终常在 0.5〜2 Pa 的过渡-分子流交界附近运行。
🙋
那 1 atm 的空气在 MEMS 的纳米通道（如 L = 100 nm）里，连续介质假设就不能用了吗？
🎓
观察得很好。常压下 λ ≈ 67 nm（N₂），L = 100 nm 的纳米通道 Kn ≈ 0.67，刚好进入过渡流区。这正是纳米流体研究中"滑移流""Knudsen 修正""Cunningham 修正"频繁出现的原因。在本工具中把 P = 10⁵ Pa（log P = 5）、L = 0.1 μm 设一下，Kn 立刻越过 0.1。Bird 的 DSMC 或格子玻尔兹曼是这类典型场景的常用工具——"气体是连续介质"的朴素直觉在这种尺度下崩坏，非常有意思。

常见问题

这些边界并非严格的理论值，而是由实验和数值计算总结出的经验阈值。Kn = 0.01 是无滑移边界与 Navier–Stokes 解的误差超过约 1% 的尺度，也是 Maxwell 一阶滑移修正能延寿的下限。Kn = 0.1 是滑移修正也开始出现 10% 以上误差的尺度，需要 Burnett 方程或 DSMC。Kn = 10 是壁碰撞完全压过分子间碰撞的区域，可以忽略碰撞项，用自由分子近似（Knudsen 解）即可。这些以 Schaaf-Chambré 图的形式从 1960 年代起就出现在教科书中，被广泛用作半导体工艺和真空工程的设计指导。
三个无量纲数之间有著名的关系：Kn = M/Re × √(πγ/2)（γ 是比热比，M 是马赫数，Re 是雷诺数）。由此可知，高速低密度（高 M、低 Re）的流动 Kn 自然变大。例如超音速高高度再入舱（M ≈ 25、Re ≈ 数百）下 Kn ≈ 0.1〜1，进入过渡-分子流区域。这就是"再入的气动加热不能用普通 CFD 计算"的原因，需要 DSMC 或 Direct Boltzmann Solver。
DSMC 由 G. A. Bird 在 1960 年代开发，是一种粒子追踪型方法，直接模拟 Boltzmann 方程的解。一个仿真粒子代表约 10⁶〜10²⁰ 个真实分子，反复执行三步：(1) 在 Δt 内让粒子自由飞行；(2) 在单元内随机选碰撞伙伴、按蒙特卡罗处理碰撞；(3) 统计性地处理与壁面的反射（镜面反射或漫反射）。它适用于所有 Kn > 0.01 的稀薄气体流动，特别在 Kn > 1 时计算量比连续介质 CFD 更小。商用软件代表有 SPARTA、dsmcFoam（OpenFOAM 扩展）、Bird 的 DS2V/DS3V 等。
最显著的是粘性系数与导热系数对压力的依赖性出现。连续流（Kn < 0.01）下 Maxwell 证明粘性系数与压力无关（η = (1/3) ρ ⟨v⟩ λ 中 ρλ 不变）。但 Kn > 0.1 时壁碰撞占主导，表观粘性系数与压力成正比下降。导热系数同样在 Kn > 0.1 下随低压而降低，这正是真空保温瓶（高真空）隔热的原理。在本工具把 P 从 100000 Pa（连续流）变到 1 Pa（自由分子流），同一 T、d、L 下区域会大幅变化。装置设计中需要保证整个运行压力范围内都在同一区域，这往往决定了特征长度 L 的选取。

现实应用

MEMS 器件设计：微加速度传感器、微陀螺仪、微泵等 MEMS 器件中，活动结构间隙为 1〜10 μm，常压空气下 Kn ≈ 0.007〜0.07，正好落在滑移-过渡边界。这导致粘性阻尼系数与 Navier–Stokes 的预测产生较大偏差，影响共振频率与 Q 值设计。给 Reynolds 方程加入 Knudsen 修正的"修正 Reynolds 方程"（Burgdorfer 1959、Fukui-Kaneko 1988）已成业界标准，用于 HDD 滑块悬浮计算和 MEMS 麦克风设计。在本工具尝试 L = 1〜10 μm、P = 10⁵ Pa，就能体会 MEMS 为什么离不开 Knudsen 修正。

真空泵与工艺腔设计：半导体的溅射、CVD、外延生长中，靶到基板的分子飞行与表面反应需要平衡，所以有意在 Kn > 1 的自由分子流区运行。例如直径 30 cm 的 PVD 溅射腔，L = 0.3 m 要满足 λ > 0.3 m，需要 P < 0.02 Pa（10⁻⁴ Torr）。在本工具设 L = 3×10⁵ μm（30 cm）、改变 P，就能直观看到何种压力才进入完全分子流。更低压的电子束光刻和 MBE 在 P < 10⁻⁵ Pa 的超高真空下实现完全分子流。

航天器再入热分析：地球大气再入过程中，80 km 以上进入稀薄区（Kn > 0.01）。航天飞机再入轨迹上，110 km 高度 Kn ≈ 10（自由分子流），85 km Kn ≈ 0.1（过渡流），60 km 以下 Kn < 0.01（连续流），呈阶段性变化。各区域需要切换不同的气动模型（自由分子→Bridging→Navier–Stokes），CFD 与 DSMC 结合的混合计算是标准做法。在本工具改变气体参数追踪 Kn 的变化，能体会再入舱气动加热为什么这么难算。

纳米技术与碳纳米管流体：气体在直径 1〜10 nm 的碳纳米管内部流动时，常压下 Kn ≫ 10，进入完全自由分子流区。此时分子-壁面相互作用（漫反射 vs 镜面反射）决定流量，称为 Knudsen 扩散现象。这在膜分离、气体传感、催化剂载体设计中很重要。在本工具尝试 L = 0.001〜0.1 μm，能确认常压下也能实现完全分子流。Maxwell 漫反射系数 σ 通常 0.8〜1，但纳米管的光滑壁面下 σ < 0.5，需要在 Knudsen 修正之上再修正。

常见误解与注意事项

最常见的误解是 认为"常压空气在任何装置中都可以用连续介质假设"。实际上 L 到 100 nm 量级时常压下 Kn 也能达到约 0.7，进入过渡流区。MEMS、纳米流体、催化孔、微加工光学元件即使在常压也需要连续介质假设的修正。在本工具中固定 P = 10⁵ Pa（log P = 5），把 L 滑到 0.1 μm，Kn 会一下子超过 0.1。在纳米尺度设计中忽视装置尺寸评估、断言"空气用 NS 没问题"，会造成致命错误。

其次常见的是 忘记"Knudsen 数大≠气体稀薄"。Kn 是 λ/L 的相对量，同一气体密度下改变装置尺寸 Kn 也会变化。反过来高密度下装置极小 Kn 也会很大，低密度下装置巨大 Kn 反而很小。例如常压空气在微通道中连续介质假设破坏，深空稀薄气体在地球规模的天文现象下却能视作连续介质。在本工具独立调节 L 和 P，就能确认两者的独立性。

最后是 过度信任"Kn 边界（0.01、0.1、10）是严格的物理边界"。这些只是经验阈值，实际会随气体种类、温度、几何、表面状态偏移大约 ±1 个数量级。例如光滑壁面下即使 Kn < 0.1 也可能出现明显滑移；粗糙壁面下即使 Kn > 0.1 仍可能呈现连续介质行为。重要设计中先用本工具确认 Kn，然后用 DSMC 或实验验证较稳妥。Kn 是"该选哪种模型"的一次性指标，最终物理量的精度由模型本身的选择和边界条件处理决定。

相关工具

平均自由程模拟器 — 气体分子动理论 平均自由程模拟器：基于温度 T、压力 P、分子直径 d、分子量 M_g 实时计算硬球气体的平均自由程 λ、平均速率 ⟨v⟩、碰撞频率 ν 与数密度 n，并以分子箱与 λ-P 对数图直... 格拉斯霍夫数模拟器 — 自然对流的特征数 实时计算自然对流的浮力与粘性之比 Gr = gβΔTL³/ν²。同时给出瑞利数 Ra=Gr·Pr、Churchill-Chu 关联的 Nu、层流/湍流判定，并可视化竖直平板自然对流，... 范德瓦尔斯气体模拟器 — 实际气体压缩因子 范德瓦尔斯气体模拟器基于状态方程 (P + a/V_m²)(V_m - b) = RT，从温度 T、压力 P、分子吸引常数 a 与排斥体积常数 b 实时计算实际气体摩尔体积 V_m、... 壅塞流模拟器 — 收缩喷管的临界条件与质量流率 壅塞流模拟器根据上游总压 P_0、总温 T_0、喉部面积 A* 和背压 P_b，使用等熵流关系实时计算收缩喷管的临界压力比、最大质量流率与喉部音速，并可视化壅塞状态，是可压缩流动的教... 法诺流动模拟器 — 带摩擦的可压缩管道流 法诺流动模拟器基于 Fanno 关系，从入口马赫数 M_1、比热比 γ、达西摩擦系数 f 与管道长径比 L/D 实时计算出口马赫数 M_2、温度比 T_2/T_1、压力比 P_2/P... 正激波关系模拟器 — 可压缩流动 正激波关系模拟器使用兰金-雨贡尼奥关系式，根据上游马赫数 M_1、比热比 γ、上游静温 T_1 和上游静压 P_1，实时计算下游 M_2、P_2/P_1、T_2/T_1 和密度比 ρ...