気体の流れを扱う工学では、スケールの違いによって適用できる理論が大きく変わります。
その境界を判断するうえで非常に重要な役割を果たすのがクヌーセン数（Knudsen number）です。
クヌーセン数は、気体分子の平均自由行程と流れの代表長さの比を表す無次元数であり、連続体近似が成立するかどうかを判定する指標として広く用いられています。
しかし「具体的に何を意味するのか」「どのような場面で重要になるのか」といった疑問を持つ方は少なくないでしょう。
本記事では、クヌーセン数の意味・定義・物理的背景から、分子運動論との関係・気体流れの領域判定・工学的応用まで、わかりやすく丁寧に解説していきます。
気体力学・マイクロ流体・真空技術などに携わる方や、これから学ぼうとしている学生・研究者の方にとって役立つ内容です。
目次
クヌーセン数とは何か：定義と基本概念の結論
それではまず、クヌーセン数の定義と基本概念について解説していきます。
クヌーセン数（Kn：Knudsen number）は、気体分子の平均自由行程（λ）と流れ場の代表長さ（L）の比として定義される無次元数です。
デンマークの物理学者マーティン・クヌーセン（Martin Knudsen、1871〜1949）にちなんで命名されました。
クヌーセン数の定義式：Kn = λ / L
λ：気体分子の平均自由行程（m）、L：流れ場または系の代表長さ（m）
クヌーセン数は、気体を「連続体」として扱えるかどうかを判定する根本的な基準です。
クヌーセン数が小さい（Kn ≪ 1）場合、分子間の衝突が非常に頻繁に起こり、気体は連続体として振る舞います。
この状態では、ナビエ・ストークス方程式をはじめとする連続体力学の方程式が適用できます。
一方、クヌーセン数が大きい（Kn ≫ 1）場合、分子間衝突よりも壁面との衝突が支配的になり、気体は希薄状態となります。
このとき連続体近似は成立せず、分子動力学（Molecular Dynamics）やボルツマン方程式などの分子論的アプローチが必要になります。
| クヌーセン数 Kn の範囲 | 流れの領域 | 適用できる理論・手法 |
|---|---|---|
| Kn < 0.001 | 連続流（continuum flow） | ナビエ・ストークス方程式（すべり無し境界条件） |
| 0.001 ≤ Kn < 0.1 | すべり流（slip flow） | ナビエ・ストークス方程式（すべり境界条件） |
| 0.1 ≤ Kn < 10 | 遷移流（transition flow） | ボルツマン方程式・DSMC 法 |
| Kn ≥ 10 | 自由分子流（free molecular flow） | 分子論的手法・運動論的手法 |
この分類はマイクロ・ナノスケールの流体工学や宇宙空間での気体挙動を理解するうえで基本的な枠組みとなっています。
平均自由行程と分子運動論の基礎
続いては、平均自由行程と分子運動論の基礎について確認していきます。
クヌーセン数を理解するためには、平均自由行程（mean free path）という概念を正確に把握することが不可欠です。
平均自由行程の定義と計算式
平均自由行程とは、気体分子が他の分子と衝突するまでに移動する平均的な距離のことです。
気体分子は常に高速で運動しており、周囲の分子と頻繁に衝突を繰り返しています。
この衝突間の平均移動距離が平均自由行程です。
平均自由行程の計算式（剛体球モデル）
λ = 1 / (√2 × π × d² × n)
d：分子の直径（m）、n：数密度（分子数/m³）
理想気体の状態方程式（p = n × k_B × T）を用いると：
λ = k_B × T / (√2 × π × d² × p)
k_B：ボルツマン定数（1.38 × 10⁻²³ J/K）、T：温度（K）、p：圧力（Pa）
この式から、圧力が低いほど・温度が高いほど平均自由行程は大きくなることがわかります。
大気圧・室温（25℃）での空気の平均自由行程は約68 nm（0.068 μm）程度です。
一方、高真空状態（10⁻³ Pa 程度）では平均自由行程は数メートル以上に達することがあります。
分子運動の速度と衝突頻度
分子運動論によれば、気体分子はマクスウェル・ボルツマン分布に従う速度で運動しています。
平均分子速度（v̄）は次の式で与えられます。
平均分子速度：v̄ = √(8 k_B T / (π m))
m：分子の質量（kg）
室温（25℃）での空気（平均分子量≈29 g/mol）の場合：v̄ ≈ 468 m/s
単位時間あたりの衝突頻度：ν = v̄ / λ（Hz）
大気圧・室温では空気分子の衝突頻度は約10⁹ 回/s（10億回/秒）のオーダーです。
このように、常温・常圧では気体分子は非常に頻繁に衝突しており、巨視的な連続体的挙動が自然に現れます。
しかし系のスケールが分子の平均自由行程に近づくと、この集団的な挙動が崩れて分子論的な効果が顕在化します。
代表長さの選び方と影響
クヌーセン数の計算において代表長さ（L）の選び方は非常に重要です。
代表長さとして何を選ぶかは問題の種類によって異なります。
管内流れでは管の直径や半径、翼型周りの流れでは翼弦長や翼の厚さ、マイクロチャネル内の流れではチャネルの幅または高さが代表長さとして用いられます。
同じ気体条件でも代表長さが変わればクヌーセン数が変わるため、問題設定に応じた適切な代表長さの選択が重要です。
また、局所的なクヌーセン数（流れの中の特定の場所での λ/L）が場所によって異なる場合もあり、このとき流れの解析が特に複雑になります。
クヌーセン数の物理的意味と気体の挙動
続いては、クヌーセン数の物理的意味と気体の挙動を確認していきます。
クヌーセン数は単なる計算指標ではなく、気体が「どのように振る舞うか」という本質的な性質を表しています。
連続体近似の成立条件
流体力学では通常、気体を分子の集まりではなく連続的な媒体として扱います。
この連続体近似が成立するためには、流れ場の代表スケールが平均自由行程よりも十分に大きい（Kn ≪ 1）という条件が必要です。
連続体近似が成立するとき、気体の状態（密度・圧力・温度・速度）は空間的に連続的な場として定義できます。
この場合、ナビエ・ストークス方程式・フィックの法則・フーリエの熱伝導則などの連続体方程式が適用可能です。
日常的に遭遇するほとんどの流体現象（室内の空気の流れ・配管内の流れ・自動車周りの流れなど）では Kn ≪ 0.001 であり、連続体近似が十分に成立します。
希薄気体効果と分子論的挙動
Kn が増大するにつれて、希薄気体効果（rarefied gas effects）が顕在化します。
すべり流領域（0.001 ≤ Kn < 0.1）では、固体壁面での「すべり（slip）」現象が生じます。
通常の連続体流では壁面での流速はゼロ（すべりなし条件）ですが、希薄効果が現れると壁面近傍の分子が壁と完全に熱平衡に達せず、壁面での流速が有限の値を持つようになります。
このすべり速度（velocity slip）および温度ジャンプ（temperature jump）が希薄気体特有の境界条件として重要です。
遷移流領域（0.1 ≤ Kn < 10）では連続体方程式も分子流の近似も成立せず、ボルツマン方程式の数値解法やDSMC法（Direct Simulation Monte Carlo）が必要となります。
自由分子流の特性
Kn ≥ 10 の自由分子流領域では、分子間の衝突よりも壁面との衝突が支配的となります。
このとき分子は互いに独立して運動し、気体の粘性・熱伝導などの巨視的物性は通常とは大きく異なる挙動を示します。
例えば、自由分子流領域では気体の熱伝導率が圧力に依存するようになります（通常の連続体では圧力依存性がない）。
この特性は真空断熱パネル（VIP）や低圧センサーの設計において重要な役割を果たしています。
クヌーセン数の工学的応用分野
続いては、クヌーセン数の工学的応用分野について確認していきます。
クヌーセン数の概念は、現代の先端技術のさまざまな分野で実際に活用されています。
マイクロ・ナノ流体工学への応用
半導体製造・バイオチップ・マイクロ熱交換器などのマイクロ電気機械システム（MEMS）では、流路の代表スケールがマイクロメートル（10⁻⁶ m）からナノメートル（10⁻⁹ m）オーダーになります。
大気圧下でもこのスケールでは Kn が 0.01〜1 程度になることがあり、希薄気体効果（すべり流・遷移流）が無視できなくなります。
マイクロチャネルでの流量・熱伝達・ガスの分離特性は、クヌーセン数を考慮した設計なしには正確に予測できません。
また、シェールガス採掘では岩盤のナノポア内での気体流れがクヌーセン数の大きい希薄流領域となることが多く、クヌーセン数を考慮した浸透率モデルが重要視されています。
宇宙・高高度飛行への応用
高高度飛行（80 km 以上の高度）や宇宙空間では、大気が非常に希薄であるためクヌーセン数が大きくなります。
低軌道衛星（高度200〜400 km）の周囲では Kn が非常に大きく、自由分子流領域に近い状態です。
この条件下での空気抵抗（エアロダイナミクス）は通常の連続体流体力学とは全く異なる計算手法が必要です。
再突入宇宙機の設計では、高高度（希薄気体）から低高度（連続体流）へと変化するさまざまなKn領域を通過するため、各領域に適した解析手法を組み合わせた設計が行われています。
真空技術・薄膜プロセスへの応用
半導体製造における化学気相成長（CVD）・物理気相成長（PVD）・エッチングプロセスなどの真空プロセスでは、低圧環境でのクヌーセン数が設計パラメータとして重要です。
プロセスチャンバー内の圧力・温度・ガス組成を調整することでクヌーセン数を制御し、薄膜の均一性・成膜レート・エッチングの選択比などを最適化します。
真空ポンプの設計においても、内部流路でのクヌーセン数を考慮した性能予測が必要です。
クヌーセン数と他の無次元数との関係
続いては、クヌーセン数と他の無次元数との関係について確認していきます。
クヌーセン数は気体力学の他の重要な無次元数と密接な関係を持っており、流れ場の総合的な評価にはこれらの関係を理解することが重要です。
マッハ数・レイノルズ数との関係
クヌーセン数はマッハ数（Ma）とレイノルズ数（Re）と次の関係があります。
クヌーセン数・マッハ数・レイノルズ数の関係
Kn ∝ Ma / Re（比例関係）
より正確には：Kn = C × Ma / Re
C：比熱比・分子モデルに依存する定数（理想気体では √(π γ / 2) 程度、γは比熱比）
この関係から、レイノルズ数が小さく（粘性支配）、かつマッハ数が大きい（高速流）場合に Kn が大きくなることがわかります。
この関係は、希薄気体効果が高マッハ数・低レイノルズ数の流れで特に重要になることを示しており、極超音速飛行体の設計に直接関係します。
プラントル数・シュミット数との関連
気体の輸送係数（粘性・熱伝導・拡散）は分子運動論から導かれるため、クヌーセン数が大きい領域ではこれらの値が変化します。
希薄気体ではプラントル数（Pr = μ cp / λ）も通常の連続体近似での値から外れることがあり、熱伝達解析に注意が必要です。
このような高Kn領域での輸送係数の修正は、MEMS デバイスの熱設計や宇宙機の熱防護システム設計で重要な課題となっています。
クヌーセン数と気体の圧縮性の関係
気体の圧縮性効果と希薄効果（クヌーセン数効果）は異なる現象ですが、同時に現れることがあります。
超音速・極超音速流れでは圧縮性（マッハ数効果）と希薄効果（クヌーセン数効果）の両方が重要となります。
これらを同時に扱う必要がある場合は、DSMC（直接シミュレーションモンテカルロ）法や格子ボルツマン法などの高度な数値解法が必要となります。
まとめ
本記事では、クヌーセン数の意味・定義・物理的背景から、平均自由行程の計算・気体流れの領域分類・工学的応用・他の無次元数との関係まで幅広く解説しました。
Kn = λ / Lという基本式が示すように、気体分子の平均自由行程と系のスケールの比が、連続体流から自由分子流まで気体の挙動を決定する根本的な指標です。
マイクロ流体・宇宙工学・真空技術・半導体製造など多くの先端分野でクヌーセン数の理解が不可欠であり、その応用範囲は今後さらに広がっていくでしょう。
本記事がクヌーセン数への理解深化に役立てば幸いです。
