沸騰モデル — 先端技術と研究動向

いくつかの重要な方向性を見ていこう。

RPIモデルはクロージャ相関に依存するマクロスケールのアプローチだが、最近ではDNS（Direct Numerical Simulation）で個々の気泡の成長・離脱を直接計算する研究が進んでいる。

Level Set法やPhase Field法で気泡界面を追跡し、界面での蒸発を直接計算する。Sato & Niceno（2013, PSI）の研究が先駆的で、1つの壁面気泡の成長・離脱・合体を再現している。

ただし計算コストが膨大で、現在は数個〜数十個の気泡に限定される。RPIモデルのクロージャ相関を改良するための基礎データ提供として活用されている。

気泡径の分布を追跡するモデルだ。実際の沸騰流では気泡径は一様ではなく、壁面で生成された小さな気泡が合体（coalescence）して大きくなったり、乱流で分裂（breakup）したりする。

ここで $f(d,t)$ は気泡径 $d$ の数密度関数、$S_{wall}$ は壁面からの核沸騰ソース項だ。MUSIG法やQMOM（Quadrature Method of Moments）で離散化して解く。

RPIモデルのクロージャ相関（気泡離脱径、核サイト密度、離脱頻度）は経験的で汎用性に乏しい。最近ではDNSデータや高解像度実験データを教師データとして、ニューラルネットワークで改良する試みがある。

Bajorek & Walkeworth（NRC, 2020s）はBayesian手法でRPIパラメータの不確かさを定量化し、安全評価に活用する研究を進めている。

マイクロチャネル（水力直径 < 1 mm）での沸騰冷却は、チップ発熱密度の増大に対応する次世代冷却技術として注目されている。ただしマイクロスケールでは閉じ込め効果（Confinement number $Co = \sqrt{\sigma / (g(\rho_l - \rho_v) D^2)}$）が支配的で、通常のRPIモデルの適用範囲外になる。

VOF法やPhase Field法での直接界面追跡が有効で、COMSOLのマイクロ流体モジュールが比較的使いやすいプラットフォームだ。