› 沸腾曲线模拟器 返回
传热模拟器

沸腾曲线模拟器 — Nukiyama 曲线与池沸腾

将水在大气压下池沸腾的热流密度按过热度的四区域模型可视化，学习 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², 对数)／蓝实线＝四区域模型／红虚线＝CHF 与莱顿弗罗斯特点／黄点＝当前状态

理论与主要公式

池沸腾热流密度 $q''$ 是过热度 $\Delta T_e = T_w - T_\text{sat}$ 的函数，呈现 Nukiyama（1934）所揭示的四个不同区域的行为。

自然对流区（$\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$$

什么是沸腾曲线模拟器

🙋
我一直以为沸腾就是水冒泡而已，可是打开模拟器看到居然有四个区域，区别这么大吗？
🎓
观察得很细。1934 年日本的努木山（Nukiyama）教授通过实验首次揭示「沸腾其实有四副面孔」。简单说，慢慢升高表面温度时，液体先做安静的自然对流，然后是剧烈冒泡的核沸腾，紧接着是传热反而变差的过渡区，最后是被蒸汽膜笼罩的膜沸腾。坐标用对数是因为热流密度变化跨好几个数量级。
🙋
核沸腾区曲线急剧上升，然后过渡区又突然下降，这到底是什么物理机制？
🎓
核沸腾区里气泡不断从表面微小凹坑产生并剧烈搅动液体，传热以 $q\propto\Delta T^3$ 的速度迅猛上升。但是过热度超过 30 K 左右后，气泡过多并聚合成贴在壁面的蒸汽块，这就是过渡沸腾。蒸汽几乎不导热，所以即使壁温在涨，传到液体的净热量反而下降。把 ΔT_e 从 20 调到 50、100，看热流密度卡片先升后降就是这个意思。
🙋
图里有条标着 CHF 的红线，越过它会发生什么？
🎓
那就是可怕的「烧毁」（burnout）。如果像电加热器那样控制的是热流密度，一旦输入超过 CHF，系统就再也维持不住核沸腾，瞬间跳到对应热流密度的膜沸腾点。竖着看图：CHF≈1 MW/m² 对应核沸腾侧 ΔT_e≈30 K，但膜沸腾侧要 ΔT_e≈600 K 以上，于是表面温度一瞬间上升几百度，锅炉管或燃料棒包壳被烧穿。所以工程上对 CHF 留 1.5～2 倍的安全余量。
🙋
煎锅上水珠会乱跑那个现象，也跟这条曲线有关吗？
🎓
正是莱顿弗罗斯特（Leidenfrost）效应。图右侧的莱顿弗罗斯特点就是稳定膜沸腾的入口。水滴底部瞬间蒸发，托起水滴的就是一层薄薄的蒸汽。蒸汽膜起隔热作用，传热急剧下降。在模拟器中设 ΔT_e=200 K，热流密度会比同样过热度的核沸腾分支低好几个数量级。把液氮甩到手上不会冻伤，也是因为温差极大下出现莱顿弗罗斯特状态。

常见问题

当热流密度超过临界热流密度（CHF）时，传热面会被蒸汽膜覆盖，传热急剧恶化，表面温度瞬时上升数百度，造成传热面烧毁。这是电加热器、锅炉管和电子器件冷却面等热流控制系统的致命失效模式，设计上通常对 CHF 留出 1.5～2 倍的安全余量。
在 BWR（沸水堆）和 PWR 事故冷却（ECCS）中，避免燃料棒表面从核沸腾向膜沸腾过渡的偏离核沸腾（DNB）现象，是安全设计的核心问题。一旦进入膜沸腾，传热系数下降几个数量级，包壳温度急升，存在燃料破损风险。安全分析要求 DNBR（DNB 比）在所有工况下保持 1.3 以上。
当水滴落到高温传热面上时，底部瞬时蒸发形成蒸汽层将水滴抬起的现象。这正对应沸腾曲线上的莱顿弗罗斯特点（q_min），越过该点进入稳定膜沸腾区。煎锅上水珠四处滚动也是同一原理。蒸汽膜的隔热作用使膜沸腾区的传热显著下降。
池沸腾是仅由传热面加热静止液体的简单体系，用 Nukiyama 曲线描述。强迫对流沸腾发生在流动液体中，还受流速、过冷度与流道几何影响，复杂得多。实际工程中的蒸汽发生器和燃料组件属于强迫对流沸腾，但作为设计基础理解，首先掌握池沸腾的四区域模型至关重要。

实际应用

火电与核电的锅炉、蒸汽发生器：发电厂在高温高压下将水汽化，传热管一旦越过 CHF 就会瞬间破裂。实际工程中用 Groeneveld 查表等关联式评估 CHF，留出至少 1.3～2.0 倍的安全余量。理解沸腾曲线正是电站安全运行的起点。

电子器件的高密度冷却：CPU、GPU 与功率半导体的发热密度已超过 100 W/cm²，空气冷却已经力不从心。直接将水或介电液喷射到芯片上的「浸没冷却」「喷射冷却」充分利用核沸腾的高传热能力（CHF 可达 MW/m² 量级），目前在数据中心和电动车逆变器中开始普及。

金属热处理（淬火）：把炽热的钢件浸入水或淬火油时，表面先进入膜沸腾，冷却很慢；温度降到莱顿弗罗斯特点以下后转入过渡沸腾再到核沸腾，进入急冷阶段。冷却速度的急变直接影响相变与组织，因此淬火油的选择与搅拌条件本质上是「让工件经历哪一段沸腾曲线」的设计问题。

低温工程（cryogenics）：使用液氦、液氮的超导磁体和航天发动机，低温侧的沸腾特性至关重要。常温壁面对低温液体往往一开始就处于膜沸腾，冷却效率较低。工程上通过多孔金属等表面处理，让系统尽快过渡到核沸腾区。

常见误解与注意事项

最常见的误解是认为「表面温度越高，传热越好」。事实上越过核沸腾的 CHF 峰值后，过渡沸腾区里温度升反而热流密度降。把模拟器的 ΔT_e 从 30 调到 60、90，热流密度卡片会发现数值反而减小。这种「明明温度升高，传热却减少」的反直觉行为，正是沸腾传热最大的特点，也是众多设计事故的温床。

其次常见的错误是把 CHF 当成「绝对的物理上限」。CHF 是「可以稳定维持核沸腾的最大热流密度」，越过它不一定立刻烧毁。在「热流控制系统（电加热器等）」中确实会瞬间跳到膜沸腾烧毁；但在「壁温控制系统（再热器、太阳能集热器等）」中会经过过渡区平稳过渡，不一定造成毁灭性故障。把模拟结果用于实机时，请务必确认控制方式是热流控制还是温度控制。

最后请注意这个模拟器使用的是「水在大气压下的简化模型」。系数 C、CHF、h_fb 都是基于教科书的代表值，工程精确预测需要分别评估 Rohsenow 关联式中的表面系数 $C_{sf}$、Zuber 的 CHF 公式、强迫对流修正、流道几何影响等。尤其压力偏离大气压时 CHF 变化很大（高压下增大，极低压下减小）。实际电厂设计始终依赖 Groeneveld 查表或已验证的子通道程序（VIPRE、COBRA 等）。

相关工具

沸腾传热计算器 沸腾传热计算器：免费在线绘制池沸腾曲线，自动识别核态沸腾、过渡沸腾等区域，计算热流密度、传热系数及临界热流密度（CHF），采用Rohsenow关联式，助力工程热分析。 黑体辐射模拟器 基于普朗克定律的黑体辐射光谱实时可视化。可调节温度，观察峰值波长移动、维恩位移定律、斯特藩-玻尔兹曼定律以及与经典紫外灾变的对比。 玻尔氢原子模型模拟器 交互式可视化玻尔氢原子模型。改变主量子数n，实时计算能级、轨道半径和谱线（巴尔末系等）。 带电粒子阱模拟器（保罗阱） 探索带电粒子在保罗阱中的捕获与逃逸过程。本模拟器通过调节交流电压、频率和直流偏压，实时演示马修方程稳定区。直观的可视化界面以红点标记当前a-q参数位置，蓝色区域清晰展示稳定捕获区，助... 康普顿散射模拟器 — 光子-电子散射的波长偏移 康普顿散射模拟器，根据入射光子能量、散射角与靶质量实时计算散射光子能量、波长偏移 Δλ、反冲电子的动能与散射角。配合散射几何示意图与 Δλ(θ) 曲线，直观理解光子-电子的量子碰撞。 海森堡不确定性原理模拟器 — Δx·Δp可视化 演示海森堡不确定性原理（Δx·Δp≥ℏ/2），展示位置与动量不确定性的此消彼长。调节波包宽度，观察空间局域化对动量分布的影响。