沸騰曲線シミュレーター 戻る
伝熱シミュレーター

沸騰曲線シミュレーター — Nukiyama 曲線とプール沸騰

加熱を上げ下げすると動作点が Nukiyama 曲線をリアルタイムに動きます。熱流束制御では CHF 超えで膜沸騰へジャンプ(=バーンアウト)、温度制御ではS字曲線を滑らかに辿る——危険なヒステリシスを目で確かめられます。

パラメータ設定
制御方式

熱流束制御=CHF超えでバーンアウト(膜沸騰へジャンプ)。温度制御=遷移域を含むS字曲線を連続的に辿る。

加熱の目標値(入力レベル)
%

アニメ中は動作点がこの目標へ向かって曲線上を移動します。0%=室温付近、100%=最大入力。

核沸騰係数 C
臨界熱流束 CHF
MW/m²
放射率 ε(膜沸騰時)
プリセット

水・大気圧(T_sat = 100°C)固定。自然対流 h_nc ≈ 1000 W/(m²·K)、膜沸騰 h_fb ≈ 30 W/(m²·K) を仮定。Zuber の CHF は水・大気圧で約 1.1 MW/m²。

計算結果(ライブ)
沸騰域
過熱度 ΔT_e
熱流束 q''
面温度 T_w
臨界熱流束 CHF
Leidenfrost 点 q_min
CHF までの余裕
沸騰曲線の走査と伝熱面の状態

上=沸騰曲線(対数 q'' vs ΔT_e)。黄点=現在の動作点/赤破線=CHF・Leidenfrost/橙矢印=バーンアウトのジャンプ。下=伝熱面(気泡→蒸気膜)。

理論・主要公式

プール沸騰の熱流束 $q''$ は過熱度 $\Delta T_e = T_w - T_\text{sat}$ に対し、Nukiyama(1934)が示した4つの域で異なる挙動を取ります。

自然対流域($\Delta T_e < 5$ K)。$h_{nc}$ は自然対流熱伝達係数:

$$q'' = h_{nc}\,\Delta T_e$$

核沸騰域($5 \leq \Delta T_e < 30$ K)。Rohsenow 式に基づく簡略形:

$$q'' = C\,\Delta T_e^{3}$$

臨界熱流束(Zuber, 水・大気圧で約 1.1 MW/m²):

$$q''_\text{CHF} \approx 0.131\,\rho_v^{1/2}\,h_{fg}\,[\sigma g(\rho_l-\rho_v)]^{1/4}$$

膜沸騰域(Bromley 式 + 放射項):

$$q'' = h_{fb}\,\Delta T_e + \varepsilon\sigma_{SB}(T_w^4 - T_\text{sat}^4)/2$$

熱流束制御で $q''_\text{入力} > q''_\text{CHF}$ になると、核沸騰のまま存在できず膜沸騰枝へ不連続にジャンプし、面温度 $T_w$ が一気に数百〜千度上昇します(バーンアウト=ヒステリシス)。

沸騰曲線シミュレーターとは

🙋
沸騰って、お湯がポコポコ沸くだけのことだと思ってたんですけど、シミュレーターを開いたら4つも領域があってびっくりしました。何がそんなに違うんですか?
🎓
いいところに気づいたね。1934年に日本の沼上(Nukiyama)博士が「沸騰には4つの顔がある」ことを実験で示したんだ。ざっくり言うと、伝熱面の温度をちょっとずつ上げていくと、まずは静かな自然対流、次にポコポコの核沸騰、急に伝熱が悪くなる遷移、最後は蒸気膜に覆われる膜沸騰、という4段階を経るんだよ。シミュレーターで「再生」を押すと、黄色い動作点が曲線の上を実際に動いていく。下のパネルでは気泡が膜に変わる様子も見えるよ。グラフが対数軸なのは、熱流束が桁違いに変わるからだ。
🙋
えっ、核沸騰のところで急に上がって、そのあと急に下がってますね。これってどういう物理現象なんですか?
🎓
核沸騰域では、面の小さな凹み(キャビティ)から気泡が次々生まれて、激しく液体をかき混ぜるから伝熱が一気に良くなる。$q\propto\Delta T^3$ で増えるんだ。でも過熱度が30度くらいを超えると、気泡が多すぎて面に貼りついた蒸気の塊になり始める。これが遷移沸騰だ。蒸気は熱を通しにくいから、面温度が上がっているのに伝わる熱量はむしろ減る。「温度制御で全S字曲線」プリセットを再生すると、動作点が頂点(CHF)を越えていったん下り坂を辿るのがはっきり見えるよ。
🙋
「CHF」って赤い線が引いてありますけど、これを超えたら何が起こるんですか?
🎓
それが恐ろしい「バーンアウト」だ。電気ヒーターのように熱流束(=入力電力)を一定に制御している場合、入力が CHF を超えた瞬間、システムは核沸騰のままでいられず、橙色の矢印のように一気に膜沸騰側へ飛ぶ。これがヒステリシスだ。「CHFまで上げてバーンアウト」プリセットを再生してごらん。動作点が核沸騰の頂点まで来た途端、面温度(T_w カード)が一気に数百度跳ね上がる。これでボイラチューブや原子炉燃料棒の被覆管が一瞬で溶ける。実務では CHF の半分くらいまでしか運転しないように設計するよ。
🙋
熱流束制御と温度制御で挙動が違うんですね。温度制御だと飛ばないんですか?
🎓
そう、そこが核心だ。壁の温度そのものを制御できる系(蒸気再加熱や太陽炉など)では、過熱度を上げると遷移域を「下り坂」として連続的に辿れる。だからジャンプは起きず、ゆっくり膜沸騰へ移る。一方、熱流束を制御する系では遷移域は不安定で存在できない——だからCHFで膜沸騰枝へ飛ぶしかない。同じ曲線でも、何を固定するかでバーンアウトするかどうかが変わるんだ。グラフ右側の Leidenfrost 点は膜沸騰の入口で、フライパンの上で水滴がコロコロ転がるあの状態。蒸気膜が断熱材になって伝熱がガクッと落ちるんだよ。

よくある質問

熱流束が臨界熱流束(CHF)を超えると、伝熱面が一気に蒸気膜で覆われて熱伝達が急激に悪化し、面温度が数百度跳ね上がって伝熱面が焼損する現象です。電気ヒーターのように熱流束を制御している場合に発生し、ボイラチューブや電子機器の冷却面で致命的な事故につながります。設計では CHF に対して 1.5〜2 倍の安全余裕を取るのが一般的です。
BWR(沸騰水型原子炉)や PWR の事故時冷却(ECCS)では、燃料棒表面が核沸騰域から膜沸騰域へ遷移する DNB(Departure from Nucleate Boiling)を避けることが安全設計の中心課題です。膜沸騰に入ると熱伝達係数が桁違いに落ち、燃料棒の被覆管温度が急上昇して破損リスクが生じます。設計では DNBR(DNB 比)という安全指標を必ず 1.3 以上に保ちます。
高温の伝熱面に水滴が落ちると、底面が瞬時に蒸発して水滴を蒸気の層で浮かせる現象です。沸騰曲線上ではちょうど Leidenfrost 点(q_min)に対応し、ここから過熱度を上げると安定した膜沸騰域に入ります。フライパンの熱さチェックで水滴がコロコロ転がるのも同じ原理です。膜沸騰では蒸気膜の断熱効果で伝熱が著しく低下します。
プール沸騰は静止液体中で伝熱面のみが熱せられる単純な系で、Nukiyama 曲線で記述されます。強制対流沸騰は流動する液体内で起こり、流速・サブクール度・流路形状の影響を受けるため、はるかに複雑です。実機の蒸気発生器や燃料集合体は強制対流沸騰ですが、設計の基礎理解として、まずプール沸騰の4域モデルを把握することが重要です。

実世界での応用

火力・原子力発電のボイラ・蒸気発生器:発電プラントでは水を高温高圧で沸騰させて蒸気を作りますが、その伝熱管が CHF を超えれば一瞬で破断します。実機では CHF を直接相関式(Groeneveld の表など)で評価し、最低でも 1.3〜2.0 倍の安全余裕を持たせて運用します。沸騰曲線の理解は、まさに発電所の安全運転の出発点です。

電子機器の高密度冷却:CPU・GPU・パワー半導体の発熱密度はすでに 100 W/cm² を超える領域に達しており、空冷では限界です。水や誘電液を直接素子に当てる「液浸冷却」や「噴流冷却」では、核沸騰の高い熱伝達能力(CHF まで使えれば 1 MW/m² 級)を活用します。データセンターや EV のインバータ冷却で実用化が進んでいます。

金属の熱処理(焼入れ):鋼を水や油に焼入れすると、最初は表面が高温で膜沸騰状態になり、ゆっくり冷却されます。温度が下がって Leidenfrost 点を割ると遷移→核沸騰へ移り、急冷段階に入ります。冷却速度の急変は組織変態に直結するため、焼入れ油の選定や攪拌条件は、まさに沸騰曲線のどの領域を通すかの設計といえます。

クライオジェニクス(極低温機器):液体ヘリウム・液体窒素を使う超伝導磁石や宇宙ロケットでは、極低温液体側の沸騰特性が重要です。常温の機器壁に対しては膜沸騰になりやすく、冷却効率が悪い領域から始まります。設計では伝熱面処理(多孔質金属など)で核沸騰域へ早く遷移させる工夫が行われます。

よくある誤解と注意点

最も多い誤解は、「面温度を上げれば上げるほど熱がよく伝わる」と思ってしまうことです。実際には核沸騰のピーク(CHF)を境に、過熱度を上げてもむしろ熱流束が減る「遷移沸騰」の領域が存在します。「温度制御で全S字曲線」プリセットを再生してみてください。動作点が頂点を越えると、熱流束カードの値が増えるどころか減っていきます。この「面温度を上げたのに伝わる熱が減る」という直感に反する挙動こそ、沸騰熱伝達の最大の特徴であり、設計事故の温床でもあります。

次に多いのが、CHF を「物理的な絶対上限」と誤解することです。CHF は「核沸騰のまま安定して維持できる上限」であって、それを超えると突然焼損するわけではありません。「熱流束制御系(電気ヒーター等)」では超えた瞬間に膜沸騰へジャンプして焼損しますが、「壁温制御系(蒸気再加熱、太陽炉等)」では遷移域を経由してゆっくり膜沸騰に移行し、必ずしも焼損しません。このシミュレーターで制御方式を切り替えると、同じ曲線でもバーンアウトするかどうかが変わることが確かめられます。実機に当てはめる際は、制御方式が「熱流束系」か「温度系」かを必ず確認してください。

最後に、このシミュレーターの式は「水・大気圧の単純化モデル」であることに注意してください。係数 C・CHF・h_fb は教科書値に基づく代表値であり、実機の正確な予測には Rohsenow 相関の表面係数 $C_{sf}$、Zuber の CHF 式、強制対流補正、流路形状の影響などを個別に評価する必要があります。特に圧力が大気圧と異なる場合、CHF は急激に変化します(高圧で増え、極低圧で減る)。実プラントの設計は必ず Groeneveld の look-up table や、検証済みの subchannel コード(VIPRE、COBRA など)で評価します。

使い方ガイド

  1. 制御方式(熱流束制御 / 温度制御)を選びます。熱流束制御はCHF超えでバーンアウト、温度制御は全S字曲線を連続的に辿ります
  2. 「再生」を押すと、動作点が加熱の目標値へ向かって沸騰曲線上を移動します。下のパネルで気泡→蒸気膜の状態が同時に見えます
  3. 臨界熱流束CHF(MW/m²)を材料・圧力条件に応じて調整。水・大気圧のプール沸騰では約1.1MW/m²が代表値です(Zuber式)
  4. 核沸騰係数と放射率を変えて曲線形状を調整し、ライブ数値(ΔTe・q''・面温度・CHF余裕)を監視します

具体的な計算例

水・大気圧、CHF=1.1MW/m²、核沸騰係数C=100の場合、核沸騰域が頂点(CHF)に達する過熱度はΔTe≈22Kです。熱流束制御で「CHFまで上げてバーンアウト」を再生すると、動作点がΔTe≈22Kに達した瞬間に膜沸騰枝へジャンプし、面温度T_wが100℃台から1000℃級へと一気に上昇します。一方、温度制御で同じ曲線を辿ると、ΔTe=22Kの頂点を越えても遷移沸騰域を連続的に下り、面温度を上げているのに熱流束q''が約1.1MW/m²から減少に転じます。この非単調性と、制御方式によるヒステリシスの有無が沸騰曲線の核心です

実務での注意点

🎬 動画で見る

物質の状態変化(相転移)とは|蒸発・凝固・昇華を物理シミュレーションで可視化
物質の状態変化(相転移)とは|蒸発・凝固・昇華を物理シミュレーションで可視化