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

沸腾曲线模拟器 — 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)确定了四个不同的区域。

自然对流区($\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$ 瞬间上升数百至上千度(烧干=滞后)。

沸腾曲线模拟器简介

🙋
沸腾就是水冒泡吧,打开模拟器才发现有四个区域,太惊讶了。差别这么大吗?
🎓
很好的观察!1934年日本科学家沼上(Nukiyama)通过实验发现了"沸腾有四张脸"这个现象。简单说,逐步升高传热面温度的话,会经历四个阶段:安静的自然对流、冒泡的核沸腾、突然传热变坏的过渡沸腾,最后是蒸汽膜覆盖的膜沸腾。按下"播放"后,黄色工作点会真的沿曲线移动,下面的面板同时显示气泡变成蒸汽膜的过程。图表采用对数坐标是因为热流量的数值变化非常大。
🙋
核沸腾区间里热流量突然上升,然后又下降,这是什么物理原理呢?
🎓
在核沸腾区,传热面的小凹坑中不断产生气泡,剧烈搅动液体,所以传热效果大幅提升。热流量遵循 $q\propto\Delta T^3$ 增长。但过热度超过30度左右时,气泡太多,开始黏在面上形成蒸汽团块,这就是过渡沸腾。蒸汽的传热能力很差,所以尽管面温度更高,传到液体的热量反而减少了。运行"温度控制全S形曲线"预设,你会看到工作点越过CHF峰值后,沿着下降的分支滑下来。
🙋
图上有条红线叫"CHF",超过这条线会怎样?
🎓
那就是可怕的"烧干"!在电热器这类热流量受控的系统中,输入的热流量一旦超过CHF,系统就无法停留在核沸腾状态,会沿着橙色箭头突然跳到膜沸腾分支。这就是滞后现象。运行"升至CHF触发烧干"预设:工作点一到达核沸腾峰值,表面温度(T_w卡片)就瞬间跳升几百度。这样锅炉管或核反应堆燃料棒的包壳会在一瞬间融化。所以设计时要留出CHF的1.5~2倍安全余裕。
🙋
这么说热流量控制和温度控制的行为不一样?温度控制就不会跳吗?
🎓
对,这正是关键。如果能直接控制壁面温度(蒸汽再热、太阳能聚焦等),提高过热度时可以把过渡区当作"下坡"连续地走过去,不会跳变,缓慢进入膜沸腾。而在热流量控制下,过渡区是不稳定的、根本无法存在——所以工作点只能在CHF处跳到膜沸腾分支。同一条曲线,固定什么量决定了会不会烧干。图右侧的Leidenfrost点是膜沸腾的入口,也就是水滴在热锅上滚动的那种状态。蒸汽膜充当隔热层,使传热急剧下降。

常见问题

当热流量超过临界热流量(CHF)时,传热面被蒸汽膜迅速覆盖,导致传热急剧恶化,表面温度跳升数百度,传热面烧毁。这种现象在热流量受控的系统(如电热器)中发生,在锅炉管和电子设备冷却面上会造成灾难性事故。设计中通常对CHF采用1.5~2倍的安全余裕。
在沸腾水反应堆(BWR)和压水反应堆(PWR)的事故冷却(ECCS)中,避免燃料棒表面从核沸腾向膜沸腾过渡的DNB(偏离核沸腾)是安全设计的核心。进入膜沸腾后,传热系数下降几个数量级,燃料棒包壳温度急剧上升,面临破损风险。设计中必须保持DNBR(DNB比)至少为1.3以上。
当水滴落在高温传热面上时,底部瞬间蒸发,形成蒸汽膜使水滴悬浮的现象。在沸腾曲线上对应于Leidenfrost点(q_min),从这点开始提高过热度就进入稳定的膜沸腾区。平底锅中水滴滚动也是同样原理。膜沸腾中蒸汽膜的隔热效应使传热大幅下降。
池沸腾发生在静止液体中,仅传热面被加热,是简单系统,由Nukiyama曲线描述。强制流动沸腾发生在流动液体中,受流速、亚冷度和流路形状影响,复杂得多。实际设备的蒸汽发生器和燃料组件采用强制流动沸腾,但设计基础理解必须先掌握池沸腾的四域模型。

实际应用

火电和核电厂的锅炉和蒸汽发生器:发电厂用高温高压水沸腾产生蒸汽,但若传热管超过CHF就会瞬间破裂。实际运行中用Groeneveld表等相关式评估CHF,设计时至少保留1.3~2.0倍的安全余裕。理解沸腾曲线是电厂安全运行的基础。

电子设备的高密度冷却:CPU、GPU、功率半导体的发热密度已超过100 W/cm²,空冷已无力。直接用水或介电液冷却芯片的"液浸冷却"或"喷流冷却"技术利用核沸腾的高传热系数(可用到CHF前约1 MW/m²级)。数据中心和电动车变流器冷却中应用前景广阔。

金属热处理(淬火):钢铁淬火时放入水或油中,初期表面处于高温膜沸腾,逐渐冷却。温度降至Leidenfrost点下方时,进入过渡→核沸腾的急冷阶段。冷却速率的跃变直接影响金属组织变态,所以淬火油的选择和搅拌条件本质上就是设计沸腾曲线的哪个区域来参与冷却。

极低温技术(液氦、液氮):超导磁体和航天火箭用极低温液体冷却,常温器壁对极低温液体总是膜沸腾态,传热效率很差。设计中通过多孔金属等表面处理让流体尽快进入核沸腾区,提高冷却效率。

常见误解和注意事项

最常见的误解是"面温度越高,传热越好"。实际上核沸腾的峰值(CHF)以上,即使增加过热度,热流量反而会减少——这就是"过渡沸腾"区的特征。运行"温度控制全S形曲线"预设,工作点越过峰值后,热流量卡片的值反而在下降。这种"温度升高,热传递下降"的反直觉行为正是沸腾热传递的最大特征,也是设计事故的根源。

另一个常见误解是把CHF当作"物理绝对上限"。CHF其实是"核沸腾能稳定维持的上限"。在电热器等"热流量受控"系统中超过CHF会立即跳到膜沸腾而烧干;但在"温度受控"系统(蒸汽再热、太阳能聚焦等)中,可能逐步经历过渡区进入膜沸腾,不一定烧干。在本模拟器中切换控制方式,就能看到同一条曲线要么烧干、要么不烧干。把模拟应用到实际设备时,一定要确认控制方式是"热流量模式"还是"温度模式"。

最后要注意,本模拟器是"水、大气压的简化模型"。系数C、CHF、h_fb都是教科书代表值,精确预测实际工况需要单独评估Rohsenow相关式的表面因子$C_{sf}$、Zuber的CHF公式、强制流动修正、流道形状影响等。特别是压力不同于大气压时,CHF会急剧变化(高压增大,极低压降低)。实际工厂设计必须用Groeneveld查表或经过验证的subchannel分析代码(VIPRE、COBRA等)来评估。

使用指南

  1. 选择控制方式(热流量控制 / 温度控制)。热流量控制在越过CHF时烧干,温度控制则连续遍历整条S形曲线
  2. 按"播放",工作点沿沸腾曲线向加热目标值移动,下面的面板同时显示气泡→蒸汽膜的状态
  3. 根据材料和压力条件调整临界热流量CHF(MW/m²)。水、大气压的池沸腾代表值约为1.1 MW/m²(Zuber式)
  4. 改变核沸腾系数和发射率以调整曲线形状,并监视实时数值(ΔT_e、q''、表面温度、CHF余裕)

具体计算示例

水、大气压下,CHF=1.1 MW/m²、核沸腾系数C=100时,核沸腾区达到峰值(CHF)的过热度约为ΔT_e≈22 K。在热流量控制下运行"升至CHF触发烧干",工作点一到达ΔT_e≈22 K就跳到膜沸腾分支,表面温度T_w从约100 ℃骤升到1000 ℃级。而在温度控制下沿同一条曲线行进时,工作点越过ΔT_e=22 K峰值后沿过渡区连续下行:表面更热,但热流量q''却从约1.1 MW/m²开始减少。这种非单调性,以及是否因控制方式而出现滞后(跳变),正是沸腾曲线的核心。

工程实践注意事项