第三节:对流与辐射换热——牛顿冷却定律、斯蒂芬-玻尔兹曼定律、辐射换热系数

各位同学,咱们今天聊点实在的。热防护系统设计,说白了就是跟热量“斗智斗勇”。你想想看,航天器再入大气层时,表面温度能飙到几千度,热量怎么进来、怎么出去,必须算得明明白白。这一节,我重点讲两个最基础的换热机制:对流和辐射。

一、牛顿冷却定律——对流换热的基本公式

对流换热,说白了就是流体流过固体表面时带走(或带来)热量的过程。我在项目中遇到过不少新手,一上来就套公式,结果算出来的热流密度跟实测差了好几倍。为什么?因为对流换热系数 h 这个参数,太容易出错了。

牛顿冷却定律的表达式很简单:

q = h · (T_w - T_f)

其中:

  • q —— 热流密度,单位 W/m²
  • h —— 对流换热系数,单位 W/(m²·K)
  • T_w —— 壁面温度,单位 K
  • T_f —— 流体主流温度,单位 K

嗯,这里要注意:公式看着简单,但 h 可不是常数。它跟流体性质(黏度、导热系数、比热容)、流动状态(层流还是湍流)、几何形状都有关系。我个人的习惯是,先判断流态,再选经验关联式,最后才代入计算。

核心要点:对流换热系数 h 是“算”出来的,不是“查”出来的。工程上常用努塞尔数 Nu 来关联:

Nu = h · L / k

其中 L 是特征长度,k 是流体导热系数。然后根据流态选择对应的 Nu 经验公式,比如 Dittus-Boelter 公式(湍流管内流动):

Nu = 0.023 · Re^0.8 · Pr^n

n = 0.4(加热流体)或 0.3(冷却流体)。

我的经验:在航天器热防护设计中,对流换热往往发生在高超声速条件下。这时候气体已经部分电离,物性参数变化剧烈。我曾经用常温下的空气物性去算,结果热流密度低估了30%以上。后来我改用参考温度法(Reference Temperature Method),才把误差控制在10%以内。

二、斯蒂芬-玻尔兹曼定律——辐射换热的基础

辐射换热跟对流完全不同。它不需要介质,真空中也能传热。说白了,任何温度高于绝对零度的物体都在向外辐射能量。斯蒂芬-玻尔兹曼定律给出了黑体的辐射力:

E_b = σ · T^4

其中:

  • E_b —— 黑体辐射力,单位 W/m²
  • σ —— 斯蒂芬-玻尔兹曼常数,5.67×10⁻⁸ W/(m²·K⁴)
  • T —— 绝对温度,单位 K

你想想看,温度从1000K升到2000K,辐射力不是翻倍,而是变成16倍!这就是为什么再入段高温下,辐射换热占绝对主导地位。

但实际物体不是黑体,所以需要引入发射率 ε:

E = ε · σ · T^4

发射率 ε 在0到1之间。我见过有人把碳化硅涂层的发射率直接取0.9,结果算出来的辐射冷却效果偏乐观。实际上,高温下发射率会随温度和表面状态变化。我记得有一次做试验,涂层表面被烧蚀后,发射率从0.85降到了0.6,辐射散热能力直接掉了近30%。

避坑指南:我曾经在计算两个表面之间的辐射换热时,直接用了两个表面的发射率相乘。这是错的!正确的做法是用系统发射率(或有效发射率),对于两个无限大平行平板:

1/ε_eff = 1/ε_1 + 1/ε_2 - 1

否则你会高估辐射换热量,导致热防护设计偏薄。

三、辐射换热系数——一个“伪”系数

为了工程计算方便,我们经常把辐射换热也写成类似牛顿冷却定律的形式:

q_rad = h_rad · (T_1 - T_2)

其中 h_rad 就是辐射换热系数。但你要明白,这只是一个线性化的处理方式。真正的辐射换热跟温度的四次方成正比,不是线性关系。

h_rad 的表达式可以从斯蒂芬-玻尔兹曼定律推导出来。假设两个表面温度分别为 T_1 和 T_2,系统发射率为 ε_eff,则:

q_rad = ε_eff · σ · (T_1^4 - T_2^4)

利用因式分解:

T_1^4 - T_2^4 = (T_1^2 + T_2^2)(T_1 + T_2)(T_1 - T_2)

所以:

h_rad = ε_eff · σ · (T_1^2 + T_2^2)(T_1 + T_2)

你看,h_rad 本身是温度的函数。温度越高,h_rad 越大。在高温下,辐射换热系数可以轻松超过对流换热系数。

实用技巧:在热防护系统初步设计阶段,我习惯先估算一个平均的辐射换热系数。比如,假设壁面温度1500K,环境温度300K,发射率0.8,那么:

h_rad ≈ 0.8 × 5.67e-8 × (1500^2 + 300^2) × (1500 + 300)
     ≈ 0.8 × 5.67e-8 × (2.25e6 + 9e4) × 1800
     ≈ 0.8 × 5.67e-8 × 2.34e6 × 1800
     ≈ 191 W/(m²·K)

这个值已经相当可观了。如果对流换热系数只有50 W/(m²·K),那辐射就是散热的主力。

四、知识体系总览

为了让你更直观地理解这一节的内容,我画了一张图,把对流和辐射换热的核心逻辑串起来:

对流与辐射换热知识体系 对流换热 牛顿冷却定律:q = h·(T_w - T_f) 核心参数:对流换热系数 h h 的确定:Nu = h·L/k 经验关联式:Dittus-Boelter 等 辐射换热 斯蒂芬-玻尔兹曼定律:E = ε·σ·T⁴ 核心参数:发射率 ε 系统发射率:1/ε_eff = 1/ε₁ + 1/ε₂ - 1 辐射换热系数:h_rad = ε_eff·σ·(T₁²+T₂²)(T₁+T₂) 工程应用:热防护系统设计 总热流密度:q_total = q_conv + q_rad 高温区:辐射主导(T⁴效应) 低温区:对流主导(需精确计算 h) 注意:辐射换热系数 h_rad 是温度的函数,高温下不可视为常数

我的建议:在实际工程中,对流和辐射往往是同时存在的。我习惯先分别计算,再叠加。但要注意,如果壁面温度很高,辐射冷却会显著降低壁面温度,从而影响对流换热的温差。所以,有时候需要迭代计算。别嫌麻烦,这一步省了,后面试验数据对不上才真头疼。

好了,这一节的内容就到这里。记住三个核心:牛顿冷却定律算对流,斯蒂芬-玻尔兹曼定律算辐射,辐射换热系数是线性化的工具。下次你设计热防护系统时,希望你能想起我今天说的这些“坑”和“经验”。


公众号:蓝海资料掘金营,微信deep3321