沸腾热传递计算机 返回
热分析

沸腾热传递计算机

核基亚玛曲线实时绘制。核沸腾·过渡沸腾·膜沸腾区域自动判定,计算临界热流密度(CHF)、莱顿弗罗斯特点、Rohsenow相关热流密度。

参数设置
流体选择
加热面·流体组合 (Csf)
Csf越小=核生成位点越多越活跃=相同ΔTe下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 的数值一下子都猛增!为什么这么敏感?
🎓
观察得好。过热度升高会激活更多核生成位点(位点密度增大),而且每个位点的脱离频率(每秒冒出几个气泡)也提高。位点密度×频率×每个气泡的潜热,决定了输运的热量,所以效果是相乘的。Rohsenow相关中核沸腾 $q''\propto\Delta T_e^3$,因此 $h=q''/\Delta T_e\propto\Delta T_e^2$,急剧上升。但继续升高,气泡会合并成蒸气膜——一旦超过CHF(临界热流密度),热传递就崩溃。徽章变红就是烧损的危险信号。
🙋
明白了!"加热面·流体组合"为什么要选Csf值?另外计算用的"Rohsenow相关"是算什么的?
🎓
好眼光!Csf是由"加热面材质·表面状态"和"流体"组合决定的实验常数,表示气泡产生的容易程度(核生成位点的活跃度)。Csf越小,位点越多越活跃,相同ΔT下冒出的气泡更多、h更高。在动画里把"铜-水(0.013)"换成"镍-水(0.020)",你会看到位点密度和 h 数值都下降。计算内核是Rohsenow相关:先求核沸腾热流密度 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超过风险。

金属淬火(急冷)工艺:高温金属部件在液体中急冷时,表面经历膜沸腾(蒸汽膜形成)→核沸腾的转变。控制这条冷却曲线能优化金属硬度和残余应力。

常见误解与注意事项

初学者使用这个工具时有几个常见误解。第一,"Rohsenow相关不是万能的"。这个公式只是核沸腾区的目测工具,实际机器设计单靠它很危险。比如流路狭窄或加热面倾斜,计算值和实测值差很大。第二,"选择实验常数Csf要谨慎"。工具里能选典型组合,但实际表面粗糙度、污垢(结垢)千差万别。磨光铜(Csf~0.013)设计的换热器,制造过程表面状态微妙变化,性能就可能低于预测。第三个陷阱是"不能不留余量地按CHF峰值设计运转"。在沸腾曲线峰值附近运转,条件稍变(如系统压力波动)就会突进过渡沸腾,导致急速升温(焼损)。实务上对CHF必须有1.5~2.0的安全系数(实际热流密度 ≤ CHF / 安全系数)。

使用指南

  1. 拖动表面过热度ΔTe(K)滑块或使用低/中/高预设,动画的核生成位点数、气泡大小、脱离节奏会即时变化(核沸腾区域约4~30K)
  2. 改变流体和"加热面·流体组合(Csf)",可对比相同过热度下位点密度、脱离径Db、热传递系数h如何变化
  3. 查看实时数值q''、h、ΔTe、Db、位点密度、脱离频率、雅各布数Ja、CHF,以及"h vs ΔTe"图上的红色工作点。用 播放/暂停/重置 控制观察

具体计算示例

饱和水(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²,动画中气泡变稀疏。

实务中的注意事项