› 沸騰熱伝達計算機 戻る
熱解析

沸騰熱伝達計算機

Nukiyama曲線をリアルタイム描画。核沸騰・遷移沸騰・膜沸騰レジームを自動判定し、臨界熱流束（CHF）・ライデンフロスト点・Rohsenow相関による熱流束を計算。

パラメータ設定
流体選択
加熱面・流体組合せ (Csf)
小さいCsf＝核生成サイトが多く活発＝同じΔTでhが高い
表面過熱度 ΔTe
K
T_s − T_sat。過熱度が上がると核生成サイト密度・離脱頻度が増えhが急上昇
プリセット（熱流束）
オーバーレイ
沸騰レジーム： 核沸騰
計算結果
—
熱流束 q'' [kW/m²]
—
熱伝達率 h [kW/m²K]
—
壁面過熱度 ΔTe [K]
—
気泡離脱径 Db [mm]
—
核生成サイト密度 [1/cm²]
—
離脱頻度 f [Hz]
—
ヤコブ数 Ja
—
CHF [MW/m²]
気泡核沸騰アニメーション（リアルタイム）
核生成中 成長中の気泡 離脱・上昇（潜熱輸送） 熱流 q''
h vs ΔTe（動作点）
理論・主要公式

核沸騰熱流束（Rohsenow, 1952）：

$$q''=\mu_l h_{fg}\left[\frac{g(\rho_l-\rho_v)}{\sigma}\right]^{1/2}\left[\frac{c_{pl}\Delta T_e}{C_{sf}h_{fg}\mathrm{Pr}^n}\right]^3$$

熱伝達率：$h=q''/\Delta T_e$（核沸騰では $q''\propto\Delta T_e^{3}$ なので $h\propto\Delta T_e^{2}$ で急増）

気泡離脱径（Fritz）：$D_b=0.0208\,\beta\sqrt{\dfrac{\sigma}{g(\rho_l-\rho_v)}}$、離脱頻度 $f$ は $f\,D_b\approx0.59\left[\dfrac{\sigma g(\rho_l-\rho_v)}{\rho_l^2}\right]^{1/4}$

ヤコブ数（過熱の指標）：$\mathrm{Ja}=\dfrac{\rho_l c_{pl}\Delta T_e}{\rho_v h_{fg}}$。臨界熱流束（Zuber）：$q''_{max}=0.131\,h_{fg}\rho_v\left[\dfrac{\sigma g(\rho_l-\rho_v)}{\rho_v^2}\right]^{1/4}$

サイト密度 N は離脱フラックスと潜熱輸送のバランス $q''\sim N\,f\,(\tfrac{\pi}{6}D_b^3\rho_v h_{fg})$ から逆算（潜熱輸送＋微小液膜蒸発の代表値）。

沸騰熱伝達計算機とは

🙋
アニメーションで加熱面からプツプツ気泡が出てますね。これが「核沸騰」なんですか？ 何が起きてるんでしょう？
🎓
そう、それが核沸騰のしくみそのものだよ。加熱面の小さな傷やくぼみ（核生成サイト）で気泡が生まれて、ふくらんで、ある大きさ（離脱径 Db）になると浮力で離脱して上昇していく。これを①核生成→②成長→③離脱の「気泡サイクル」と呼ぶんだ。離脱した気泡は内部に蒸気（潜熱）を抱えて運び去る。さらに気泡が液体をかき混ぜるので、加熱面が冷たい液で次々と洗われる。だから核沸騰は熱伝達率 h が桁違いに高い。右上に出てる h と q'' の数字に注目してね。
🙋
過熱度ΔTのスライダーを上げると、気泡の数も h の数字も一気に増えますね！ なぜそんなに効くんですか？
🎓
いい観察だ。過熱度が上がると、活性化する核生成サイトの密度が増え、しかも各サイトの離脱頻度（1秒あたり何個気泡を出すか）も上がる。サイト密度×頻度×1個あたりの潜熱、で運ばれる熱が決まるから、効果が掛け算で効く。Rohsenow相関だと核沸騰は $q''\propto\Delta T_e^3$、つまり $h=q''/\Delta T_e\propto\Delta T_e^2$ で急上昇する。だから「中」プリセットから「高」に上げるだけで h が何倍にも跳ね上がる。ただし上げ続けると気泡が合体して蒸気膜になり、CHF（臨界熱流束）を超えると一気に熱が伝わらなくなる——バッジが赤くなったら危険信号だよ。
🙋
なるほど！「加熱面・流体組合せ」でCの値を選ぶのはなぜですか？ あと、計算に使われてる「Rohsenow相関」って何を計算してるんですか？
🎓
良いところに気が付いたね。Csfは「加熱面の材質と表面状態」と「流体の組み合わせ」で決まる実験定数で、気泡の生まれやすさ（核生成サイトの活発さ）を表している。Csfが小さいほどサイトが多く活発で、同じΔTでも気泡がたくさん出て h が高くなる。アニメで「銅-水(0.013)」から「ニッケル-水(0.020)」に変えると、サイト密度と h の数字が下がるのが見えるはずだ。計算の中身はRohsenow相関で、君が設定したΔTから核沸騰熱流束 $q''$ を出し、$h=q''/\Delta T_e$ を求めている。右の「h vs ΔTe」グラフの赤い動作点も一緒に動くから、効果を見比べてみて。

よくある質問

気泡が湧き出す点の数は核生成サイト密度（過熱度Csfから推定）を、気泡の最大サイズはFritz式による離脱径Dbを、湧き出すテンポは離脱頻度fを反映しています。これらの積が潜熱輸送量＝熱流束q''に対応します。過熱度ΔTeを上げるとサイト密度と頻度がともに増え、熱伝達率hが急上昇する様子が視覚的に確認できます。バッジが赤（CHF接近）になると気泡が密集し、膜沸騰への移行が近いことを示します。
C_sfとnは加熱面の材質と流体の組み合わせに依存する実験定数です。一般的な値（例：水-銅でC_sf=0.013、n=1.0）がデフォルトで設定されていますが、文献値を参考に手動で変更可能です。精度向上には実測データとのフィッティングをお勧めします。
CHFはZuberの式（q''_CHF=0.131ρ_g^{1/2}h_fg[gσ(ρ_l-ρ_g)]^{1/4}）を、ライデンフロスト点は飽和温度に基づく経験式を使用しています。計算結果はグラフ上にマーカーで表示され、過熱度の閾値として参照できます。
本ツールでは飽和温度における代表値を用いています。温度依存性が強い系（高過熱度など）では誤差が生じる可能性があるため、必要に応じて飽和温度近傍の値を入力してください。より高精度な解析には、温度依存モデルを別途組み込むことを検討ください。

実世界での応用

原子炉冷却系の安全解析：核燃料棒の冷却性能を評価する上で、臨界熱流束（CHF）を超えて膜沸騰に陥る「バーンアウト」は絶対に避けなければなりません。CAE解析では、炉心内の局所的な熱流束がCHFを下回っていることを詳細に検証します。

高性能コンピューティングの液浸冷却：AIサーバーやスーパーコンピューターのCPU/GPUは、核沸騰領域を積極利用した液浸冷却で発熱を処理します。沸騰熱伝達計算は、冷却液の選定や沸騰が起こる最適温度の設計に不可欠です。

ボイラーや蒸発器の熱設計：効率的に蒸気を発生させるため、伝熱管の表面温度を核沸騰が活発に起こる範囲に保つ設計がなされます。過熱度が低すぎると自然対流のみ、高すぎるとCHF超過のリスクがあります。

金属の焼入れ（急冷）プロセス：高温の金属部品を液体中で急冷する際、部品表面では膜沸騰（蒸気膜形成）から核沸騰へと移行します。この冷却曲線を制御することで、金属の硬さや残留応力を最適化できます。

よくある誤解と注意点

このツールを使い始める際、特に現場からCAEを学ぶ方によく見られる誤解がいくつかあります。まず第一に、「Rohsenow相関は万能ではない」という点。この式は核沸騰領域の目安を出すためのもので、実際の機器設計ではこれだけに頼るのは危険です。例えば、流路が狭い場合や加熱面が傾いている場合、計算結果と実測値が大きく乖離することがあります。第二に、「実験定数Csfの選択は慎重に」。ツールでは代表的な組み合わせを選べますが、実際の表面粗さや汚れ（スケーリング）は千差万別です。磨かれた銅（Csf~0.013）の値で設計した熱交換器が、製造段階で微妙に表面状態が変わるだけで性能が予測より落ちる、ということは珍しくありません。第三の落とし穴は「CHFを安全マージンなしで目標値に設定しない」こと。沸騰曲線のピークぎりぎりで運転すると、わずかな条件変動（例えば、システム圧力の変動）で遷移沸騰に突入し、急激な温度上昇（焼損）を引き起こします。実務では、CHFに対して少なくとも1.5〜2.0の安全率（実際の熱流束 ≦ CHF / 安全率）を見込むのが常識です。

関連ツール

沸騰曲線シミュレーター — Nukiyama 曲線とプール沸騰 沸騰曲線（Nukiyama 曲線）シミュレーターで、過熱度に対するプール沸騰の熱流束を自然対流・核沸騰・遷移・膜沸騰の4域でリアルタイム計算。CHF とバーンアウト現象を直感的に学べ... フィン熱伝達シミュレーター NovaSolverのフィン熱伝達シミュレーターは、矩形、三角形、環状フィンからフィンアレイまで、様々な形状の効率と温度分布をリアルタイムで可視化。最適なフィン設計と間隔をシミュレー... 輻射熱伝達計算 輻射熱伝達計算ツールでStefan-Boltzmann則に基づく放射熱流束を簡単シミュレーション。表面温度や放射率、面積を入力するだけで熱流束をリアルタイム計算。熱流束曲線と輻射抵抗... ヒートシンク設計計算機（電子冷却・ジャンクション温度） NovaSolverの無料オンライン計算機で、ヒートシンクの熱設計を効率化。電子部品のジャンクション温度を、自然対流/強制空冷に対応したフィンアレイ熱抵抗モデルでリアルタイム計算。B... 有限差分熱伝導シミュレーター 陽解法有限差分法 T^(n+1)_i = T^n_i + r(T^n_i+1-2T^n_i+T^n_i-1) をリアルタイムアニメーション。フーリエ数・安定性条件r≤0.5・ディリク...

使い方ガイド

  1. 過熱度ΔTe(K)スライダーまたはプリセット（低/中/高）を動かすと、アニメーションの核生成サイト数・気泡の大きさ・離脱テンポが即座に変化します。核沸騰域は概ね4～30Kが対象です
  2. 流体と「加熱面・流体組合せ(Csf)」を変えると、同じ過熱度でもサイト密度・離脱径Db・熱伝達率hがどう変わるか比較できます
  3. 右の「h vs ΔTe」グラフの赤い動作点と、上部のq''・h・ΔTe・Db・サイト密度・離脱頻度・ヤコブ数Ja・CHFのライブ数値で物理量を確認できます。再生/一時停止/最初からで観察を制御できます

具体的な計算例

飽和水(100℃、1気圧)が銅面上で核沸騰する場合、ΔTe=15K、Rohsenow相関式(Csf=0.013、n=1.0)によりq''≈466.8kW/m²、h≈31.1kW/m²Kと算出されます（CHF≈1.111MW/m²）。このときFritz式の気泡離脱径Db≈2.34mm、離脱頻度f≈39Hz、ヤコブ数Ja≈45で、潜熱輸送バランスから核生成サイト密度はおよそ130/cm²（≈1.3×10⁶/m²）と逆算されます。ΔTe=8Kに下げるとq''≈71kW/m²・h≈8.9kW/m²Kまで落ち、サイト密度も約20/cm²に減ってアニメの気泡がまばらになります。

実務での注意点

🎬 動画で見る

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