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

克努森数 模拟器 — 稀薄气体流动分类

通过温度 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）；4种色带 = 连续体、滑流、过渡、自由分子；黄圆 = 当前 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} \lt 0.01 \,\,(\text{连续体}),\quad 0.01 \le \mathrm{Kn} \lt 0.1 \,\,(\text{滑流})$$ $$0.1 \le \mathrm{Kn} \lt 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 吗？
🎓
没错。在溅射和分子束外延（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 气压空气中的纳米通道（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.6~0.7。这个区间需要 DSMC 或格子 Boltzmann 方程，"气体是连续体"的朴素直觉彻底失效。

常见问题

这些边界值不是严格的物理常数，而是由实验和数值计算经验确立的经验指标。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）自动导致大 Kn。如航天器再入大气层（M ≈ 25、Re ≈ 数百）时，高度 110 km 处 Kn ≈ 10（自由分子），85 km 处 Kn ≈ 0.1（过渡），60 km 以下 Kn < 0.01（连续体）。各高度需要不同的空气动力学模型（自由分子→桥接→Navier–Stokes），CFD 与 DSMC 混合求解，这就是"为何再入热保护计算如此困难"的根本原因。
DSMC 由 G. A. Bird 在 1960 年代开发，直接数值模拟 Boltzmann 方程。一个模拟粒子代表约 10⁶~10²⁰ 个真实分子。每步迭代：(1) 粒子自由飞行 Δt；(2) 在网格单元内随机选择碰撞对，蒙特卡罗处理分子碰撞；(3) 统计处理壁面反射（镜面 vs 漫反射）。对任何 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 后转为压力相关性降低，这是魔瓶（高真空）的隔热原理。本工具中固定 T、d、L，把 P 从 100 kPa（连续体）变到 1 Pa（自由分子），清晰观察分类转变。装置设计要确保工作压力全范围内流动分类不变。

实际应用

MEMS 器件设计：微加速度传感器、微陀螺、微泵等 MEMS 中可动部件的间隙通常 1~10 μm，常压空气下 Kn ≈ 0.007~0.07，处于滑流~过渡流边界。粘性阻尼系数偏离 Navier–Stokes 预测，影响共振频率和品质因数。业界采用"修正 Reynolds 方程"（Burgdorfer 1959、Fukui-Kaneko 1988），在 HDD 浮动滑块和 MEMS 麦克风设计中标配，本工具设定 L = 1~10 μm、P = 10⁵ Pa 能直观说明为何 MEMS 必须修正。

真空泵和工艺腔体设计：半导体溅射、CVD、外延等工艺需要 Kn > 1 的自由分子流，保证靶离子向基板的直线飞行和表面反应平衡。PVD 溅射腔（直径 30 cm）中，要满足 L = 0.3 m 对应 λ > 0.3 m，需 P < 0.02 Pa（10⁻⁴ Torr）。本工具设定 L = 3×10⁵ μm、扫动 P，直观看到进入完全分子流的压力门槛。电子束刻蚀、MBE 等更低压应用达 P < 10⁻⁵ Pa 的超高真空。

航天器再入热分析：地球大气再入中，高度 80 km 以上进入稀薄区（Kn > 0.01）。航天飞机再入轨迹中，110 km 高度 Kn ≈ 10（自由分子）、85 km 处 Kn ≈ 0.1（过渡）、60 km 以下 Kn < 0.01（连续体）。各区段需用不同空气动力学模型（自由分子→桥接→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、纳米流体、催化孔、微细加工光学器件中，即使是常压空气也要加 Knudsen 补正。本工具中固定 P = 10⁵ Pa（常压），把 L 拉到 0.1 μm，会看到 Kn 瞬间超过 0.1，进入过渡流。装置尺度评估不周、武断认为"气体选择 NS 方程就 OK"，对纳米尺度设计是致命错误。

次常见误区是"Kn 大表示气体稀薄"的混淆。Kn = λ/L 是相对量，同样气体密度，改变 L 后 Kn 也变。反之亦然。常压微通道中 Kn 可能很大，低密度大型天体现象（地球大气风场）Kn 可能很小。本工具独立调控 L 和 P，验证两者对 Kn 的平等地位。

第三是过度信任 Kn 边界值的"严格性"。Kn = 0.01, 0.1, 10 是经验指标，实际可因气体种类、温度、几何、壁面粗糙度等偏差 ±1 个数量级。滑光壁即使 Kn < 0.1 也能见显著滑移，粗糙壁即使 Kn > 0.1 有时连续性好。重要设计应在本工具确认 Kn 后，用 DSMC 或实验验证，Kn 只是模型选择的一级指标，物理量精度最终由模型本身和边界条件决定。

相关工具

平均自由程 模拟器 — 气体分子运动论 平均自由程模拟器通过温度、压力、分子直径和分子量，实时计算刚球气体的平均自由程 λ、平均速率 ⟨v⟩、碰撞频率 ν 和数密度 n，可视化分子... 格拉晓夫数 模拟器 — 自然对流特性数 实时计算浮力与粘性之比 Gr = gβΔTL³/ν²。瑞利数 Ra=Gr·Pr、Churchill-Chu Nu 相关式、层流/湍流判定，用于竖平板自然对流... 范德瓦尔斯状态方程 模拟器 — 实际气体压缩因子 通过温度、压力、分子间吸引力 a 和排斥体积 b，计算实际气体摩尔体积、压缩因子、临界温压... 流量堵塞 模拟器 — 收缩喷管临界条件和质量流量 通过上游全压、全温、喉部面积、背压，计算收缩喷管临界压力比、最大质量流量、喉部音速和堵塞状态... 范诺流 模拟器 — 摩擦压缩管道流 通过入口马赫数、比热比、Darcy 摩擦因子、长径比，计算出口马赫、温度比、压力比... 正激波关系式 模拟器 — 压缩流 通过上游马赫数、比热比、上游静温、上游静压，用 Rankine–Hugoniot 关系式计算下游 M_2、P_2/P_1、T_2/T_1...

使用指南

  1. 温度滑块（slT）在 100~1000 K 范围调节。以氮气 N₂ 为例，从室温 300 K 到加热条件下观察分子热速变化。
  2. 压力滑块（slLogP）在 10⁻⁶~10⁵ Pa 对数刻度调节，追踪从大气压 101325 Pa 到真空区的过渡。
  3. 分子直径（slD_pm）以皮米单位输入。N₂ 约 370 pm，He 约 220 pm。
  4. 特征长度 L（slL_um）以微米设定，实时计算平均自由程 λ 和克努森数 Kn = λ/L。
  5. 输出流动分类判定：Kn < 0.01 为连续体流；0.01 ≤ Kn ≤ 10 为过渡流；Kn > 10 为自由分子流。

计算示例

微流体器件设计，氮气流道（d = 370 pm、T = 300 K、P = 101325 Pa、L = 100 μm）：平均自由程 λ ≈ 68 nm，Kn ≈ 0.00068，为连续体流。同一流道减压至 P = 1 Pa 时，λ ≈ 6.8 μm，Kn ≈ 0.068，进入过渡流。再减压至 P = 0.01 Pa 时，λ ≈ 680 μm，Kn ≈ 6.8，接近自由分子流。本工具还显示连续体极限长度 Lc = 0.01 × λ，即连续体假设失效的临界尺度。

工程应用注意

  1. MEMS 真空泵和微热交换器设计中，过渡流域（Kn = 0.01~10）使传统 Navier–Stokes 方程失效，出现速度滑移。Kn ≥ 0.1 时强制加入滑移系数补正。
  2. 航天器再入热防护评估，高空（约 100 km）进入自由分子流，冲击波厚度消失，温度分布非平衡，必须切换到分子动力学模拟。
  3. 压力低于 10⁻⁶ Pa 的超高真空区，分子碰撞频率极低，水银扩散泵失效，改用离子泵。
  4. 温度变化时平均自由程修正：λ ∝ √T，T = 600 K 时 λ 增加 1.41 倍，流动分类可能从连续体流跃入过渡流。