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

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

水・大気圧でのプール沸騰における熱流束を、過熱度の関数として4域モデルで可視化。CHF とバーンアウト現象を学べます。

パラメータ設定
過熱度 ΔT_e
K
核沸騰係数 C
—
臨界熱流束 CHF
MW/m²
放射放射率 ε（膜沸騰時）
—

水・大気圧（T_sat = 100°C）固定。自然対流域 h_nc ≈ 1000 W/(m²·K)、膜沸騰 h_fb ≈ 30 W/(m²·K) を仮定。

計算結果
—
沸騰域
—
熱流束 q
—
沸騰熱伝達係数 h_boil
—
CHF までの余裕
沸騰曲線（対数 q vs ΔT_e）

横軸＝過熱度 ΔT_e (K, 対数)／縦軸＝熱流束 q (W/m², 対数)／青実線＝4域モデル／赤水平線＝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 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$$

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

🙋
沸騰って、お湯がポコポコ沸くだけのことだと思ってたんですけど、シミュレーターを開いたら4つも領域があってびっくりしました。何がそんなに違うんですか？
🎓
いいところに気づいたね。1934年に日本の沼上（Nukiyama）博士が「沸騰には4つの顔がある」ことを実験で示したんだ。ざっくり言うと、伝熱面の温度をちょっとずつ上げていくと、まずは静かな自然対流、次にポコポコの核沸騰、急に伝熱が悪くなる遷移、最後は蒸気膜に覆われる膜沸騰、という4段階を経るんだよ。グラフで対数軸になっているのは、熱流束の値が桁違いに変わるからだ。
🙋
えっ、核沸騰のところで急に上がって、そのあと急に下がってますね。これってどういう物理現象なんですか？
🎓
核沸騰域では、面の小さな凹み（キャビティ）から気泡が次々生まれて、激しく液体をかき混ぜるから伝熱が一気に良くなる。$q\propto\Delta T^3$ で増えるんだ。でも過熱度が30度くらいを超えると、気泡が多すぎて面に貼りついた蒸気の塊になり始める。これが遷移沸騰だ。蒸気は熱を通しにくいから、面温度が上がっているのに伝わる熱量はむしろ減る。シミュレーターで ΔT_e を20→50→100と動かすと、熱流束カードの値がいったん上がって減るのが見えるよ。
🙋
「CHF」って赤い線が引いてありますけど、これを超えたら何が起こるんですか？
🎓
それが恐ろしい「バーンアウト」だ。電気ヒーターのように熱流束（=入力電力）を一定に制御している場合、入力が CHF を超えた瞬間、システムは核沸騰のままでいられず、いきなり膜沸騰の対応点まで飛ぶ。グラフを縦に見て——CHF（約 1 MW/m²）の高さで、核沸騰側は ΔT_e≈30K、膜沸騰側は ΔT_e≈600K以上。つまり面温度が一気に数百度跳ね上がる。これでボイラチューブや原子炉燃料棒の被覆管が一瞬で溶けるんだ。実務では CHF の半分くらいまでしか運転しないように設計するよ。
🙋
フライパンで水滴がコロコロ転がるあれも、関係あるんですか？
🎓
まさに Leidenfrost 効果だ。グラフの右側の Leidenfrost 点が、その入口にあたる。水滴の底面が瞬時に気化して、水滴を蒸気の薄い膜が浮かべる状態だ。膜沸騰域では蒸気膜が断熱材として働くから、伝熱がガクッと落ちる。シミュレーターで ΔT_e=200K に設定すると、核沸騰域の同じ過熱度（実際は届かないけど）と比べて、熱流束が桁違いに小さくなることがわかるよ。極低温の液体窒素を素手に少しかけても凍傷にならないのも、皮膚との温度差が極端で 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）を境に、過熱度を上げてもむしろ熱流束が減る「遷移沸騰」の領域が存在します。シミュレーターで ΔT_e を 30→60→90 と動かしてみてください。熱流束カードの値が増えるどころか減っていきます。この「面温度を上げたのに伝わる熱が減る」という直感に反する挙動こそ、沸騰熱伝達の最大の特徴であり、設計事故の温床でもあります。

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

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

関連ツール

沸騰熱伝達計算機 Nukiyama曲線（沸騰曲線）をリアルタイム描画。表面過熱度ΔTeを入力して核沸騰・遷移沸騰・膜沸騰レジームを自動判定。Rohsenow相関による熱流束・熱伝達率・CHF・ライデン... 黒体放射スペクトルシミュレーター プランクの法則による黒体放射スペクトルをリアルタイム可視化。温度を変えてピーク波長移動・ウィーンの変位則・ステファン-ボルツマン則・古典的紫外破綻との比較ができます。 ボーア水素原子モデル ボーア水素原子モデルをインタラクティブに可視化。主量子数nを変えてエネルギー準位・軌道半径・スペクトル線（バルマー系列等）をリアルタイム計算。 荷電粒子トラップシミュレーター（ポールトラップ） NovaSolverの荷電粒子トラップシミュレーターで、ポールトラップの原理をリアルタイムに可視化。電圧や周波数を調整し、荷電粒子の捕捉・脱出の動きを体験しながら、イオントラップや質... コンプトン散乱 シミュレーター — 光子-電子散乱の波長シフト コンプトン散乱シミュレーターは、入射光子エネルギー・散乱角・標的質量から散乱光子のエネルギー、波長シフト Δλ、反跳電子の運動エネルギーと散乱角をリアルタイム計算します。Δλ(θ) ... ハイゼンベルクの不確定性原理シミュレーター 位置の不確定性Δxと運動量の不確定性Δpの積がℏ/2以上であることを波動関数（ガウス波束）で視覚化。量子力学の根本原理を体験。