クヌーセン数の概念を理解したうえで、次に重要なのは「実際にどうやって計算するか」という実践的なスキルです。
クヌーセン数の計算は基本公式自体はシンプルですが、平均自由行程の算出・代表長さの設定・各パラメータの取得という各ステップで正確な理解が求められます。
本記事では、クヌーセン数の計算方法を公式の確認から始まり、平均自由行程の計算手順・代表長さの設定方法・具体的な計算例・希薄気体領域の判定まで、実践的に解説していきます。
気体力学・マイクロ流体・真空工学の設計・解析に携わる方にとって役立つ内容です。
目次
クヌーセン数の基本公式と計算手順の全体像
それではまず、クヌーセン数の基本公式と計算手順の全体像を解説していきます。
クヌーセン数の計算は次の基本公式から出発します。
クヌーセン数の基本公式：Kn = λ / L
λ：気体分子の平均自由行程（m）、L：系または流れ場の代表長さ（m）
計算の流れ：①平均自由行程 λ を計算する → ②代表長さ L を設定する → ③Kn = λ / L を計算する → ④流れ領域を判定する
この4ステップを正確に踏むことで、クヌーセン数を精度よく求めることができます。
各ステップの詳細を順を追って確認していきましょう。
平均自由行程の計算方法
続いては、平均自由行程の計算方法を確認していきます。
クヌーセン数計算の最初の重要ステップが平均自由行程（λ）の算出です。
剛体球モデルによる計算式
気体分子を剛体球として扱う最も基本的なモデルでは、平均自由行程は次の式で計算します。
平均自由行程の計算式
λ = k_B × T / (√2 × π × d² × p)
k_B：ボルツマン定数 = 1.381 × 10⁻²³ J/K
T：温度（K）
d：分子の有効直径（m）
p：圧力（Pa）
代表的な気体の分子直径（有効衝突直径）の目安
窒素（N₂）：d ≈ 3.7 × 10⁻¹⁰ m（0.37 nm）
酸素（O₂）：d ≈ 3.5 × 10⁻¹⁰ m（0.35 nm）
空気（平均）：d ≈ 3.7 × 10⁻¹⁰ m（0.37 nm）として近似
ヘリウム（He）：d ≈ 2.6 × 10⁻¹⁰ m（0.26 nm）
この式からわかるように、圧力が低くなるほど・温度が高くなるほど平均自由行程は長くなります。
例えば、大気圧（101325 Pa）・室温（293 K）での空気の平均自由行程を計算します。
空気の平均自由行程の計算例（大気圧・室温）
λ = (1.381 × 10⁻²³ × 293) / (√2 × π × (3.7 × 10⁻¹⁰)² × 101325)
= (4.05 × 10⁻²¹) / (1.414 × 3.1416 × 1.369 × 10⁻¹⁹ × 101325)
≈ 6.1 × 10⁻⁸ m ≈ 61 nm
（実測値はおよそ 65〜70 nm 程度であり、近似として妥当な値です）
圧力・温度の変化と平均自由行程の変化
設計や解析ではさまざまな圧力・温度条件を扱うため、平均自由行程がどのように変化するかを把握しておくことが重要です。
| 圧力（Pa） | 温度（K） | 空気の平均自由行程 λ（目安） | 条件の例 |
|---|---|---|---|
| 101325（1 atm） | 293（20℃） | 約 65 nm | 大気圧・室温 |
| 101325（1 atm） | 1000 | 約 222 nm | 大気圧・高温 |
| 1000 | 293 | 約 6.5 μm | 低圧（約0.01 atm） |
| 1 | 293 | 約 6.5 mm | 高真空（約10⁻⁵ atm） |
| 0.001 | 293 | 約 6.5 m | 超高真空 |
この表から、圧力が1000分の1になると平均自由行程は1000倍になるという比例関係が確認できます。
真空技術の設計では、この圧力依存性の理解が非常に重要です。
粘性係数から平均自由行程を推算する方法
分子直径が不明な場合や、より実用的な計算が求められる場合には、気体の粘性係数（μ）から平均自由行程を推算する方法が使えます。
粘性係数を用いた平均自由行程の推算
λ = μ × √(π / (2 ρ R_specific T))
または：λ ≈ μ / (p) × √(π R T / (2 M))
μ：気体の粘性係数（Pa・s）、ρ：気体密度（kg/m³）、R：一般気体定数（8.314 J/mol・K）、M：モル質量（kg/mol）
空気（20℃、1 atm）：μ ≈ 1.81 × 10⁻⁵ Pa・s → λ ≈ 66 nm（剛体球モデルと一致）
粘性係数は多くの工学データブックに記載されているため、この推算法は実用的に便利です。
代表長さの設定方法
続いては、代表長さの設定方法を確認していきます。
クヌーセン数の計算において代表長さ（L）の選択は、解析の目的と系の幾何学的構造に応じて適切に行う必要があります。
問題の種類別の代表長さの選び方
代表長さは問題の種類によって異なります。
主要な場合の代表長さの選び方を以下に示します。
問題別の代表長さの選び方
円管内の流れ：L = D（内径）または L = r（半径）
平行平板間の流れ：L = h（平板間隔）
マイクロチャネル：L = H（チャネルの短辺方向の長さ）
球・粒子周りの流れ：L = d_p（粒子直径）
翼型・流線形物体：L = c（翼弦長）または L = t（最大厚さ）
多孔質体内の流れ：L = d_pore（細孔直径）または L = d_grain（粒子直径）
一般原則として、気体分子が「感じる」スケールに最も近い幾何学的特徴長さを代表長さとして選ぶことが重要です。
例えばマイクロチャネルでは、チャネルの長さではなく幅・高さなどの断面方向のサイズを使います。
局所クヌーセン数の概念
一様でない流れ場では、場所によってクヌーセン数が異なる局所クヌーセン数の概念が重要になります。
例えば、再突入宇宙機の周囲では先端部（よどみ点付近）と後方（後流域）ではクヌーセン数が異なり、領域によって適用すべき解析手法が変わります。
また、圧縮性流れでは衝撃波後方で密度が急増し、局所的な平均自由行程が変化するため局所クヌーセン数の評価が必要です。
代表長さとクヌーセン数の感度
代表長さが10倍変わればクヌーセン数も10倍変わるため、代表長さの選択は流れ領域の判定に直接影響します。
代表長さの選択によるクヌーセン数の変化例
条件：空気（λ = 65 nm）、マイクロチャネル（幅 10 μm、高さ 50 μm、長さ 1 mm）
幅を代表長さとした場合：Kn = 65 × 10⁻⁹ / 10 × 10⁻⁶ = 6.5 × 10⁻³ → すべり流領域
高さを代表長さとした場合：Kn = 65 × 10⁻⁹ / 50 × 10⁻⁶ = 1.3 × 10⁻³ → すべり流領域（連続体に近い）
長さを代表長さとした場合：Kn = 65 × 10⁻⁹ / 1 × 10⁻³ = 6.5 × 10⁻⁵ → 連続体領域
→ 断面方向の最短寸法（ここでは幅 10 μm）を代表長さとすることが適切です。
この例からわかるように、代表長さの選択が流れ領域判定を左右するため、問題の物理を理解したうえで適切に選ぶことが求められます。
クヌーセン数の具体的な計算例
続いては、クヌーセン数の具体的な計算例を確認していきます。
実際の問題に対してクヌーセン数を計算する手順を複数の例で確認します。
計算例1：MEMS デバイスの流路
【計算例1：MEMSマイクロポンプの流路】
条件：窒素ガス（N₂）、温度 300 K、圧力 50000 Pa（0.5 atm）、流路幅 5 μm
Step1：窒素の平均自由行程を計算
d(N₂) = 3.7 × 10⁻¹⁰ m
λ = (1.381 × 10⁻²³ × 300) / (√2 × π × (3.7 × 10⁻¹⁰)² × 50000)
≈ 4.14 × 10⁻²¹ / (6.11 × 10⁻²⁰) ≈ 0.068 μm ≈ 130 nm（0.5 atm での値）
Step2：代表長さ L = 5 × 10⁻⁶ m（流路幅）
Step3：Kn = 130 × 10⁻⁹ / 5 × 10⁻⁶ = 0.026
Step4：0.001 ≤ Kn = 0.026 < 0.1 → すべり流領域
→ ナビエ・ストークス方程式にすべり境界条件を適用する必要があります。
計算例2：低軌道衛星の周囲
【計算例2：高度300 km での衛星周りの気体流れ】
条件：高度 300 km の大気（数密度 n ≈ 3 × 10¹⁵ /m³、温度 ≈ 1000 K）、衛星の代表長さ L = 1 m
高高度大気の平均自由行程（近似）：λ ≈ 1 / (√2 × π × d² × n)
d(N₂) = 3.7 × 10⁻¹⁰ m として：
λ ≈ 1 / (1.414 × 3.1416 × (3.7 × 10⁻¹⁰)² × 3 × 10¹⁵) ≈ 550 m
Kn = 550 / 1 = 550 ≫ 10 → 自由分子流領域
→ 連続体流体力学は完全に不適用。分子論的手法が必要。
計算例3：CVD プロセスチャンバー
【計算例3：半導体 CVD プロセスチャンバー】
条件：シラン（SiH₄）ガス、温度 900 K、圧力 100 Pa、チャンバー寸法 L = 0.3 m
SiH₄ の分子直径：d ≈ 4.1 × 10⁻¹⁰ m（近似値）
λ = (1.381 × 10⁻²³ × 900) / (1.414 × 3.1416 × (4.1 × 10⁻¹⁰)² × 100)
≈ 1.24 × 10⁻²⁰ / (7.48 × 10⁻¹⁹) ≈ 0.017 m ≈ 17 mm
Kn = 0.017 / 0.3 = 0.056
0.001 ≤ Kn = 0.056 < 0.1 → すべり流領域
→ 連続体方程式が使えますが、すべり境界条件と温度ジャンプの考慮が推奨されます。
クヌーセン数計算の注意点と実用的なヒント
続いては、クヌーセン数計算の注意点と実用的なヒントを確認していきます。
正確なクヌーセン数の計算と適切な流れ領域判定のために、いくつかの実用的な注意点があります。
混合気体での平均自由行程の扱い
実際のガスは純粋な単一成分ではなく、複数の成分からなる混合気体であることが多いです。
混合気体では各成分分子の直径・質量・数密度が異なるため、平均自由行程の計算がより複雑になります。
実用的な近似として、混合気体の平均粘性係数から推算した平均自由行程を使う方法があります。
空気は主に窒素（78%）・酸素（21%）の混合気体ですが、近似的に単一成分として扱い d ≈ 3.7 × 10⁻¹⁰ m を使うことが工学的に広く行われています。
温度・圧力が変化する流れでの扱い
圧縮性流れや加熱・冷却を伴う流れでは、温度・圧力が場所によって変化するため局所的な平均自由行程も変化します。
この場合、局所条件（局所温度・局所圧力）を用いて各点での λ を計算し、局所クヌーセン数を評価することが適切です。
特に衝撃波前後では密度・温度が急変するため、クヌーセン数の急激な変化に注意が必要です。
判定がボーダーラインの場合の対応
Kn が各領域の境界付近（例：Kn ≈ 0.001 または Kn ≈ 0.1）にある場合は、両側の解析手法を両方適用して比較することが推奨されます。
より高精度な解析が必要な場合は、遷移流領域（0.1 ≤ Kn < 10）を直接扱えるDSMC法や格子ボルツマン法を使うことで、連続体近似の誤差を排除した評価が可能です。
まとめ
本記事では、クヌーセン数の計算方法について、基本公式・平均自由行程の計算手順・代表長さの設定方法・具体的な計算例・注意点と実用的ヒントまで詳しく解説しました。
Kn = λ / L という公式を正確に計算するためには、平均自由行程を気体の温度・圧力・分子直径から算出し、問題の物理に即した代表長さを設定することが不可欠です。
大気圧・室温の空気では λ ≈ 65 nm という基準値を念頭に置きつつ、設計対象の圧力・温度条件と代表スケールに応じてクヌーセン数を評価することが実務の第一歩です。
本記事がクヌーセン数の計算スキル習得に役立てば幸いです。
