気体の流れを正確に解析するためには、その流れがどの領域(連続流・すべり流・遷移流・自由分子流)に属するかを正確に判定することが出発点となります。
この判定を可能にするのがクヌーセン数(Kn)です。
特に Kn = 0.001・0.1・10 という境界値の意味と、各領域での気体の性質・適用できる解析手法・工学的含意を理解することは、マイクロ流体・宇宙工学・真空技術の実務において非常に重要です。
本記事では、クヌーセン数による流れ領域の判定基準を中心に、連続流から自由分子流までの各領域の特徴・気体の性質の変化・適用手法の選択まで丁寧に解説していきます。
目次
クヌーセン数による流れ領域の判定:全体像と結論
それではまず、クヌーセン数による流れ領域の判定の全体像について解説していきます。
クヌーセン数の値によって気体の流れは4つの領域に分類されます。
クヌーセン数による流れ領域の分類
Kn < 0.001:連続流(連続体近似が十分成立)
0.001 ≤ Kn < 0.1:すべり流(連続体方程式+すべり境界条件)
0.1 ≤ Kn < 10:遷移流(ボルツマン方程式・DSMC 法が必要)
Kn ≥ 10:自由分子流(分子論的手法が必要)
この分類は気体の圧力・温度・系のスケールのさまざまな組み合わせに対して適用できる普遍的な枠組みです。
日常の大気圧環境では代表長さが 1 mm 以上あれば Kn ≪ 0.001 であり、連続流として扱えます。
しかしマイクロスケール・低圧・高温という条件になるにつれて Kn が増大し、希薄気体効果が無視できなくなります。
連続流領域(Kn < 0.001)の特徴
続いては、連続流領域の特徴について確認していきます。
Kn < 0.001 の連続流領域は、日常的な流体工学の問題のほとんどが該当する最も一般的な領域です。
連続体近似と支配方程式
連続流領域では、分子間の衝突が非常に頻繁に起こるため気体は集団として連続的な媒体として振る舞います。
この状態では、ナビエ・ストークス方程式が気体の運動量保存を正確に記述します。
また、境界条件として壁面での「すべりなし条件(no-slip condition)」が成立します。
すべりなし条件とは、壁面に接する流体の速度が壁面の速度と等しい(相対速度ゼロ)という条件であり、連続体流では非常によい近似です。
エネルギー方程式についても、壁面での温度は壁面の温度と等しい(熱平衡)という境界条件が成立します。
連続流における気体の粘性・熱伝導の特性
連続流領域では、気体の粘性係数・熱伝導率は圧力に依存せず、温度のみの関数です。
これは分子運動論から導かれる重要な性質であり、粘性係数と熱伝導率は温度の平方根に比例(サザーランドの式など)します。
この圧力独立性は、広い圧力範囲にわたって連続流として扱える範囲内では成立しますが、希薄気体領域(Kn が大きい領域)では崩れてしまいます。
連続流での代表的な事例
連続流が成立する代表的な工学事例を確認します。
| 事例 | 代表長さ | 条件 | Kn の目安 |
|---|---|---|---|
| 自動車周りの流れ | 約 2 m | 大気圧・室温 | 約 3 × 10⁻¹¹ |
| 航空機翼周りの流れ | 約 2 m(翼弦長) | 高度 10 km | 約 10⁻⁸ |
| 配管内の流れ(内径 10 mm) | 10 mm | 大気圧・室温 | 約 6.5 × 10⁻⁶ |
| 小型ファンの流れ | 数 cm | 大気圧・室温 | 約 10⁻⁶ 以下 |
これらはすべて Kn ≪ 0.001 であり、通常の流体力学の教科書で学ぶ連続体理論が完全に適用できます。
すべり流領域(0.001 ≤ Kn < 0.1)の特徴
続いては、すべり流領域の特徴を確認していきます。
すべり流領域は、連続体方程式が基本的に成立するが壁面での境界条件に修正が必要な領域です。
速度すべりと温度ジャンプ
すべり流領域では、壁面近傍で気体分子が壁と完全な熱平衡に達しないため、速度すべり(velocity slip)と温度ジャンプ(temperature jump)が生じます。
速度すべりとは、壁面直上での流体の速度が壁速度と一致しなくなる現象です。
速度すべりの境界条件(Maxwell のすべりモデル)
u_slip = (2 − σ_v) / σ_v × λ × (∂u/∂n)|_wall
σ_v:速度調整係数(accommodation coefficient、0〜1、通常 0.85〜1.0 程度)
n:壁面法線方向
温度ジャンプの境界条件:
T_slip − T_wall = (2 − σ_T) / σ_T × (2γ / (γ+1)) × λ / Pr × (∂T/∂n)|_wall
σ_T:熱調整係数、γ:比熱比、Pr:プラントル数
これらのすべり境界条件をナビエ・ストークス方程式と組み合わせることで、すべり流領域の流れを比較的精度よく解析できます。
すべり流領域での工学的事例
すべり流領域に該当する代表的な工学問題を確認します。
マイクロチャネル(幅数 μm〜数十 μm)内の大気圧での気体流れ、低圧環境(10²〜10³ Pa 程度)での通常サイズの配管内の流れ、高高度(50〜80 km 程度)での小型飛行体の流れなどが代表例です。
MEMSデバイスの設計ではすべり流効果を考慮した流量補正が重要であり、通常の連続体理論で計算した流量よりも実際の流量が大きくなることが知られています。
これをクヌーセン効果による流量増加(Knudsen minimum)と呼び、マイクロポンプやマイクロバルブの設計に直接影響します。
遷移流領域(0.1 ≤ Kn < 10)と自由分子流領域(Kn ≥ 10)
続いては、遷移流領域と自由分子流領域について確認していきます。
これらの領域では連続体近似が成立しないため、より高度な解析手法が必要です。
遷移流領域の特性と解析手法
遷移流領域(0.1 ≤ Kn < 10)は、連続体流の近似も自由分子流の近似も適用できない最も解析が難しい領域です。
この領域では、分子間の衝突と壁面との衝突がほぼ同程度の頻度で起こり、気体の振る舞いが複雑になります。
適用できる解析手法として、ボルツマン方程式の数値解、DSMC(Direct Simulation Monte Carlo)法、格子ボルツマン法(LBM)があります。
DSMC 法は各気体分子の運動・衝突をモンテカルロ法的に統計的に模擬する手法であり、遷移流〜自由分子流領域で最も広く使われる強力な方法です。
ただし計算コストが高く、大規模な系への適用には高性能計算資源が必要です。
自由分子流領域の特性
自由分子流領域(Kn ≥ 10)では、分子同士の衝突がほとんど起こらず、各分子は独立して運動します。
この領域では気体の巨視的物性(粘性・熱伝導率)の概念自体が変化します。
特徴的な現象として、熱泳動力(thermophoresis)やクロークス力(Crookes force)などの自由分子流特有の力が重要になります。
これらの力は真空計の設計・粒子分離装置・宇宙機の姿勢制御などで考慮が必要です。
また、自由分子流領域では熱伝達が圧力に依存するようになり、真空断熱パネル(VIP)の断熱性能がこの効果を利用しています。
クヌーセン最小値現象
管内を流れる気体の体積流量をクヌーセン数の関数として描くと、Kn ≈ 1 付近で流量が最小値(クヌーセン最小値)を示すという興味深い現象があります。
連続流領域(Kn → 0)でも自由分子流領域(Kn → ∞)でも流量は大きく、遷移流付近で最小となります。
この非単調な挙動は、マイクロ流体デバイスの設計において意外な落とし穴になることがあり、注意が必要です。
まとめ
本記事では、クヌーセン数による流れ領域の判定(連続流・すべり流・遷移流・自由分子流)について、各領域の特徴・気体の性質・適用できる解析手法・工学的事例まで幅広く解説しました。
Kn = 0.001・0.1・10 という3つの境界値を頭に入れておくことで、気体流れの問題に対して適切な解析手法を迅速に選択できるようになります。
マイクロ流体・真空技術・宇宙工学などの先端分野では、連続体近似の適用限界を正確に把握し、必要に応じて希薄気体効果を考慮した解析を行うことが設計品質の鍵となります。
本記事がクヌーセン数と流れ領域の判定理解に役立てば幸いです。