連続体近似は、関連する最小の流動長さが分子スケールよりも十分に大きく、かつ局所的な分子平衡が維持されている場合に妥当となります。 実際には、これは主にクヌーセン数（(\mathrm{Kn}=\lambda/L)）を用いて評価されます。ここで、(\lambda) は分子の平均自由行程、(L) は速度、温度、濃度、または圧力が著しく変化する最小の特性長さです。一般的な化学工学やバイオプロセスのパイロット装置では、連続体輸送は通常信頼できますが、主な例外として高真空、マイクロスケールのチャネル、多孔質構造、薄膜、および急峻な界面領域が挙げられます。
連続体モデルは、(\mathrm{Kn}\ll 1) であり、局所平均化体積内に揺らぎを抑制するのに十分な分子が含まれ、かつ局所熱力学平衡が成立している場合に信頼できます。 関連する長さとは、必ずしも容器や配管の全体直径ではなく、最小の輸送または勾配スケールです。
連続体近似に必要な条件
分子平均化が意味を持つこと
密度、速度、濃度などの連続体場は、個々の分子の特性ではなく、小さな体積内での平均値を表します。
その平均化体積には、ランダムな揺らぎが小さくなるのに十分な分子が含まれている必要があります。相対的な平衡揺らぎは、通常 (N^{-1/2}) に比例して減少します（(N) は平均化体積内の分子数）。
(N\gtrsim 10^4) といった基準は実用的な工学的推定値として役立ちますが、普遍的な物理的閾値ではありません。許容される値は、要求される精度、測定対象の特性、およびプロセスの統計的感度に依存します。
平均化スケールが局所的であること
平均化体積は、分子の寸法と比較して大きく、かつマクロな場が変化するスケールと比較して小さくなければなりません。
これにより、スケール分離の要件が生じます：
[ \text{分子スケール} \ll \ell_{\text{average}} \ll L_{\text{gradient}} ]
中間スケールが存在しない場合、滑らかな局所密度、速度、または濃度場を一貫して定義することはできません。
局所平衡が成立していること
ニュートンの粘性法則、フーリエの熱伝導法則、フィックの拡散法則などの構成式は、分子が局所的な熱力学的状態を確立するのに十分な衝突を繰り返していることを前提としています。
この要件は、装置が幾何学的に大きくても、気体が十分に希薄である場合や、プロセスが極めて急速な時間的・空間的変化を生じさせる場合には満たされないことがあります。
クヌーセン数が主要な制限となる
定義と物理的意味
クヌーセン数は以下の通りです：
[ \mathrm{Kn}=\frac{\lambda}{L} ]
ここで、(\lambda) は分子の平均自由行程、(L) は関連する最小のマクロな長さです。
気体の場合、(L) はチャネル幅、境界層厚さ、混合ゾーン厚さ、粒子径、または温度や濃度勾配の長さスケールとなります。したがって、容器の直径のみを使用すると、クヌーセン数を過小評価し、誤った判断を招く可能性があります。
連続流領域
一般的な工学的分類は以下の通りです：
- (\mathrm{Kn}<0.01)： 通常の連続流。滑りなし条件のナビエ・ストークス方程式および標準的な輸送法則が一般的に適切です。
- (0.01\lesssim\mathrm{Kn}<0.1)： すべり流領域。連続体方程式が依然として有用な場合がありますが、壁面すべりや修正された熱・物質移動境界条件が必要になることがあります。
- (0.1\lesssim\mathrm{Kn}<10)： 遷移領域。運動論的補正なしでは連続体方程式は信頼できなくなります。
- (\mathrm{Kn}\gtrsim10)： 自由分子流領域。分子論的または運動論的な記述が必要です。
これらの境界は鋭い物理的断絶ではなく、実用的な範囲です。
気圧が制限を制御する
気体の場合、平均自由行程は圧力が低下するにつれて、概ね以下の関係で増加します：
[ \lambda \propto \frac{T}{P} ]
（気体組成が一定で、異常な非理想状態から離れている場合）
その結果、大気圧下で安全に連続体領域にある反応器や配管でも、真空下ではすべり流や遷移流に入る可能性があります。したがって、装置サイズが一定であっても、圧力を下げることで、バルク方程式自体が無効になる前に、通常の連続体境界条件が無効になることがあります。
液体は分子スケールがはるかに小さい
液体では、分子間隔はおよそ (10^{-10}) 〜 (10^{-9}\,\mathrm{m}) のオーダーであり、有効な分子平均自由行程は一般的な装置寸法と比較して極めて短いです。
これが、パイロットスケールの容器、配管、熱交換器、膜モジュールにおける液体の流れが通常、連続体方程式で十分に記述される理由です。それでも、ナノメートルスケールの閉じ込め、非常に薄い液膜、高度に構造化された界面、または非局所的な分子効果が支配的なシステムでは、連続体の限界が現れることがあります。
パイロット装置における実用的な長さスケール
従来のパイロットプラントの寸法
配管、攪拌槽、塔、熱交換器、バイオリアクターは、一般的にミリメートルからメートルの特性寸法を持っています。
これらの寸法は、通常の圧力下での気体の平均自由行程や液体の分子スケールよりも何桁も大きいです。したがって、従来のパイロットプラントの輸送設計において、連続体近似が制限要因となることは通常ありません。
小さなチャネルと多孔質媒体
マイクロチャネルにおける関連する (L) は、接続されたプラントのサイズではなく、その流体力学的寸法です。
充填層、触媒細孔、膜細孔、および多孔質支持体では、細孔径や局所的な狭窄部が妥当性を左右する可能性があります。したがって、マクロな装置であっても、連続体輸送が局所的に失敗する領域が含まれることがあります。
界面と境界層
気固、気液、または液固界面の近くでは、特性勾配長さが装置寸法よりもはるかに小さくなる場合があります。
気体の場合、厚さが平均自由行程に近づく境界層や濃度遷移領域では、すべり、運動論的、または界面補正が必要です。液体の場合、バルク相が連続体であっても、分子スケールの界面領域では分子シミュレーションや特殊な構成モデルが必要になることがあります。
その他の物理的条件の重要性
強い勾配は局所的な構成法則を無効にする
連続体仮説だけでは、ニュートン粘性、フーリエ伝導、フィック拡散が正確であることは保証されません。
特性が分子平均自由行程と同程度の距離で大幅に変化する場合、これらの法則の基礎となる局所勾配近似は疑わしいものとなります。気体では、これはしばしば (\lambda |\nabla T|/T) や (\lambda |\nabla \rho|/\rho) といった勾配長さクヌーセン数を通じて表現されます。
時間スケールの分離が必要
分子は、マクロなプロセスが変化するよりも速く、局所平衡に向かって緩和しなければなりません。
急速な圧縮、膨張、加熱、反応、または混合は、この仮定を破る可能性があります。そのような場合、流体は依然として有用な平均的記述を持つかもしれませんが、平衡熱力学的特性や通常の構成係数は局所的に適用できなくなる可能性があります。
連続体はニュートン流体を意味しない
流体は、非ニュートン性、粘弾性、異方性、反応性、または多相性であっても、連続体であり得ます。
例えば、ポリマー溶液や高濃度細胞懸濁液は、速度場や濃度場が連続体として扱える場合でも、ニュートン粘性よりも複雑な構成モデルが必要になることがあります。
多相・生物学的システムには追加の注意が必要
気液分散系、エマルション、結晶化流体、細胞培養、バイオフィルムは、相または混合物スケールで連続体としてモデル化される場合がありますが、離散相は追加の長さと時間スケールを導入します。
各バルク相に対して連続体近似が妥当であっても、単純な単相モデルでは不十分な場合があります。界面抗力、粒径分布、すべり速度、凝集、細胞スケールの影響、および界面を通じた物質移動には、明示的な扱いが必要になることがあります。
トレードオフの理解
分子数ルールはあくまで推定
(N\ge 10^4) といった規定の最小値は統計的ノイズを判断するのに役立ちますが、連続体妥当性の普遍的な境界として提示されるべきではありません。
分子数、クヌーセン数、局所平衡時間スケールは、それぞれ異なる失敗モードに対処します。ある領域に多くの分子が含まれていても、平均自由行程が勾配スケールに対して小さくない場合は、運動論的モデリングが必要になることがあります。
分子径は一般的な気体の基準ではない
気体の場合、関連する相互作用の長さは通常、圧力と温度に強く依存する平均自由行程です。
分子径は分子スケールの基準として有用ですが、装置サイズを分子径と直接比較しても連続体の妥当性は確立されません。低圧条件下では、気体分子は衝突の間に多くの分子径分を移動できます。
標準方程式は境界条件の修正が先に必要になる場合がある
(\mathrm{Kn}) が非常に小さい連続体限界から増加する際、最初に失敗するのは多くの場合、バルク全体のバランス法則の枠組みではなく、滑りなし条件や局所平衡境界条件です。
すべり流モデルは、連続体計算を中程度のクヌーセン数領域まで拡張できる場合があります。遷移領域や自由分子流領域では、ボルツマン方程式、直接シミュレーションモンテカルロ法（DSMC）、その他の分子スケール手法といった運動論的アプローチが必要です。
輸送係数は自動的に一定ではない
連続体領域内であっても、粘度、熱伝導率、拡散係数、密度、比熱は、温度、圧力、組成、または濃度によって大きく変化する可能性があります。
不適切な条件で評価されたバルク係数を使用すると、連続体仮定自体が妥当であっても大きな誤差が生じる可能性があります。したがって、連続体の妥当性と特性モデルの妥当性は個別にチェックする必要があります。
目標に向けた正しい選択
設計や解析において最小の物理スケールを使用し、パイロットプラント全体のサイズだけに頼らず、局所的な領域を評価してください。
- 主な焦点が従来の液体パイロット装置である場合： 流体を連続体として扱い、選択したレオロジー、混合モデル、および特性相関が実際の液体や生物学的懸濁液を表現していることを確認してください。
- 主な焦点が気圧付近の気体輸送である場合： 最小のチャネルや勾配スケールが平均自由行程よりもはるかに大きい場合、標準的な連続体ナビエ・ストークス方程式、エネルギー方程式、および種保存式が一般的に適切です。
- 主な焦点が真空または低圧気体操作である場合： 動作圧力と温度で (\lambda) を再計算し、(\mathrm{Kn}) を評価し、結果が (0.01) をはるかに下回らなくなった場合は、すべりや運動論的補正を含めてください。
- 主な焦点がマイクロチャネル、ナノ細孔、膜、または薄膜である場合： 局所的な細孔、隙間、または膜厚を (L) として使用し、クヌーセン数と分子平均化の要件の両方をテストしてください。
- 主な焦点が多相またはバイオプロセス輸送である場合： 各相の連続体の妥当性を個別に確認し、必要に応じて追加の界面、粒径分布、レオロジー、または生物学的モデルを選択してください。
- 主な焦点が急速な反応や極端な勾配である場合： 装置が幾何学的に大きくても局所平衡が失敗する可能性があるため、(\mathrm{Kn}) に加えて分子緩和時間や勾配の時間スケールをチェックしてください。
連続体近似は本質的にスケール分離の仮定です。分子運動がモデル化対象のマクロな輸送現象よりもはるかに小さく、かつ速い場合に信頼できます。
要約表：
| 条件または領域 | 典型的な基準 | 推奨されるアプローチ |
|---|---|---|
| 通常の連続流 | Kn < 0.01 | ナビエ・ストークス方程式および標準的な輸送法則が一般的に適切 |
| すべり流領域 | 0.01 ≤ Kn < 0.1 | 壁面すべりや修正された境界条件を用いた連続体方程式を使用 |
| 遷移流 | 0.1 ≤ Kn < 10 | 運動論的補正または分子スケールモデリングを適用 |
| 自由分子流 | Kn ≥ 10 | ボルツマン方程式やDSMCなどの運動論的手法を使用 |
| 分子平均化 | 平均化体積内に十分な分子数 | 意味のある低ノイズな連続体特性を確保 |
| 局所平衡 | プロセス変化よりも速い分子緩和 | 急速な反応、圧縮、加熱、混合をチェック |
| 臨界長さスケール | 最小の細孔、隙間、膜、勾配厚さを使用 | 容器や配管の全体直径のみに頼らない |
化学工学やバイオプロセスのパイロットプラントを計画中ですか？適切な輸送モデルと装置の選択に関する専門的なサポートについては、LABPARKにお問い合わせください。LABPARKは、化学工学、バイオプロセス・バイオテクノロジー、環境・水処理用途向けの教育・職業訓練用単位操作パイロットプラントを提供しています。当社のシステムは、大学、研究機関、企業が流体の流れ、熱・物質移動、混合、ろ過、膜プロセス、その他の単位操作を、実用的で応用重視の装置を用いて評価するのを支援します。パイロットプラントのスケール、動作圧力、流体特性、およびプロセス目標を、教育、研究、プロセス開発、人材育成のための信頼できる実験プラットフォームと一致させるお手伝いをいたします。