辐射换热 从绝对温度到等效系数与线性化边界
把概念、模型计算与实测证据分开，阅读时留意条件、单位和比较口径。
章节速览9 节 · 展开选择
集总热容把物体近似成一个温度，用对流换热得到指数冷却。但同一表面还会发出辐射，也会吸收周围的辐射。当辐射与对流处于相近量级时，仅把旧公式中的h换成一个猜测常数，可能同时漏掉温度依赖和环境温度的区别。
本篇研究一个边界清楚的模型：温度均匀、不透明的灰表面，完全看到温度均匀的黑体环境；介质不参与辐射，发射率作为给定常数。 我们先分清发出、吸收和净换热，再推导等效系数及其近似误差。全文数值均为合成算例，不对应某个真实材料牌号、温度传感器或寿命试验。
同一表面一边发出、一边吸收；净功率的符号应由两边温度决定。
区分精确因式分解、平均温度近似与把系数冻结在环境温度的近似。
先固定温度、面积与表面条件
表面温度记作Tₛ，辐射环境温度记作Tᵣ，辐射面积为A。本文环境是等温黑体包围面，表面全部辐射视域都落在这个环境上，不存在一部分朝向别的温度的物体。灰表面意味着在所采用的模型中，用一个不随波长变化的发射率ε表征表面，并使用相应的吸收率关系。
黑体单位面积发射功率为σT⁴；灰表面发射为εσTₛ⁴。在这里，环境给表面的辐照度为σTᵣ⁴，表面吸收其中的ε部分。因此净向外功率是发出减去吸收。基础关系与黑体环境的说明见文献1。
Q̇_rad = εσA(Tₛ⁴ − Tᵣ⁴)
以离开物体为正，Q̇单位W，A单位m²，Tₛ与Tᵣ必须用K。ε无单位，取0至1；本文计算用σ≈5.670374419×10⁻⁸W/(m²·K⁴)，是常数十进制值的截取，原常数信息见文献2。
- Tₛ、Tᵣ｜绝对温度
- 分别属于表面与辐射包围环境，均用K。空气温度Tₐ是另一项边界条件，不能自动代替Tᵣ。
- ε、σ｜表面参数与常数
- ε在本例固定为0.8，仅为合成参数。σ带有K⁻⁴单位；改变温度零点后不能沿用四次方数值。
- A、Q̇｜面积与功率
- 面积用于辐射能量账本。Q̇为W，累计传递能量还需对时间积分，不能把瞬时功率直接当作J。
这不是任意两块灰表面都能直接使用的表达式。若面对的是另一块有限、非黑表面，几何视角与多次反射会进入能量账；MIT的两无限灰平板结果就含两个发射率组成的分母，见文献3。把那一分母省去，再把ε随意换成两个发射率的乘积，并不一般成立。
为什么摄氏温差可用，摄氏温度的四次方不可用
400K与300K分别是126.85°C与26.85°C，温差都是100K或100°C。温差数值相同，是两种标度的单位间隔相同；绝对温度的零点却相差273.15K，换算依据见文献4。
所以在线性的hA(Tₛ−Tₐ)中，温差可用摄氏度数值计算；而Tₛ⁴−Tᵣ⁴必须先把两边转换成K。不能把126.85⁴−26.85⁴放进R1后仍使用同一个σ，也不能把四次方之差写成(400−300)⁴。
检查量纲也能发现问题：σ的单位是W/(m²·K⁴)，乘上m²和K⁴才能得到W。量纲检查不能识别所有错误，但它能提醒我们，“温度差相同”不是任意非线性表达都不变的理由。
一个400K表面的完整手算
取ε=0.8，A=0.1m²，Tₛ=400K，Tᵣ=300K。先分别计算两笔正的功率，再按约定相减：
| 能量项 | 代入关系 | 结果，W |
|---|---|---|
| 表面发出 | 0.8×σ×0.1×400⁴ | 116.129268 |
| 吸收环境辐射 | 0.8×σ×0.1×300⁴ | 36.744026 |
| 净向外 | 发出−吸收 | 79.385242 |
若只写116.129268W，就把环境吸收这一笔漏掉了。若交换两边温度，净功率变为−79.385242W，说明物体净吸热；它仍然在发射辐射，只是吸收得更多。两边同温时净功率为0，同样不代表两笔功率分别为0。
自检1：300K表面对400K环境，是否因为“比较冷”就停止发射？
没有停止。发出与吸收各自仍可为正；按离开物体为正的约定，净辐射为负，表示吸收大于发出。要报告净换热，应保留两笔能量，而不是只判断表面有没有发射。
等效辐射系数：先做精确分解
四次方之差可以分解为(Tₛ−Tᵣ)(Tₛ+Tᵣ)(Tₛ²+Tᵣ²)。因此在同一组表面与环境温度下，可以把R1精确地改写成类似对流的形式：
Q̇_rad = A h_rad (Tₛ − Tᵣ)
h_rad单位W/(m²·K)。这是代数恒等式，不是把四次方近似成一次方；当Tₛ=Tᵣ=T时，上式连续给出h_rad=4εσT³，但净功率仍为0。
本例h_rad=7.938524W/(m²·K)，于是0.1×7.938524×100≈79.385242W，与四次方直接计算一致。在这个温度点，还可以写辐射热阻R_rad=1/(Ah_rad)，但它随温度改变；不能把指定温度点的等效热阻当作整个冷却过程的常数。
当ε=0时，模型中的辐射通道关闭，h_rad=0。此时热阻形式应理解为无穷大，而不是继续做除以0的数值计算。相等温度时，直接从因式表达取得连续系数，也比用Q̇/[A(Tₛ−Tᵣ)]得到0/0更清楚。
平均温度近似的误差能算出来
令Tₘ=(Tₛ+Tᵣ)/2，ΔT=Tₛ−Tᵣ。代入R2，可得h_rad=4εσTₘ³+εσTₘ(ΔT)²。只有在第二项足够小时，才近似采用hₘ=4εσTₘ³。
|Q̇ₘ − Q̇_rad| / |Q̇_rad| = (ΔT)² / [4Tₘ² + (ΔT)²]
要求A>0、ε>0、两边温度为正且不相等；Q̇ₘ=AhₘΔT。右侧无单位。若净功率为0，相对误差本身未定义，不能把0/0直接填成0；其小温差极限为0。该式只量化这一项近似，不能代表灰体假设的总误差。
400K与300K时，Tₘ=350K，hₘ≈7.779754W/(m²·K)，净辐射近似值约77.797537W，比精确值低2%。这2%来自100²/[4×350²+100²]，不是观察图形后估计的百分比。
若改为把系数冻结在环境300K，得到h₀=4εσ300³≈4.899203W/(m²·K)，对同一100K温差只预测48.992035W，低约38.29%。两种计算都常被叫作“线性化”，参照温度不同，误差也明显不同。应写明参照温度与允许误差，而不是只写“温差不太大”。
自检2：写成AhΔT以后，是否就可以沿用一个固定时间常数？
不能直接沿用。精确h_rad同时依赖表面与环境温度，冷却时通常在变化。只有限定温度附近、固定边界并明确采用局部常系数近似后，才有相应的近似时间常数。
与对流并存：两个通道相加，而非串联
若同一面积还与空气对流换热，空气温度为Tₐ，给定对流系数h，则瞬时对流功率为hA(Tₛ−Tₐ)。两个通道同时从表面交换能量，功率应相加。本文为便于核算，假定对流面积与辐射面积相同；真实几何不同时，要分别保留各自面积。
C为恒定热容，单位J/K，时间用s；没有内部热源、相变和其他传热通道，且物体内部近等温。h单位W/(m²·K)。空气与辐射环境分别由Tₐ、Tᵣ给出；内部近等温条件仍须核对，不能直接照搬仅含对流的模型适用结论。
继续取400K表面，令Tₐ=Tᵣ=300K，h=10W/(m²·K)。对流向外100W，辐射向外79.385242W，总向外179.385242W。若C=1000J/K，初始降温速率约−0.179385K/s。这里计算的是这一时刻的斜率，没有假定整个过程保持此速率。
在300K附近，辐射的局部系数约4.899203W/(m²·K)。相同环境温度下的小扰动可用h与该局部系数之和，得到近似时间常数约671.176818s；若完全忽略辐射，只考虑对流，时间常数为1000s。二者的差异属于这组模型参数，不是实际器件响应的实测结论。
flowchart LR A[同一个表面温度] --> B[向空气对流] A --> C[与包围环境交换辐射] B --> D[按同一正方向相加功率] C --> D D --> E[由热容决定瞬时温度变化]
空气温度与辐射环境温度不同会怎样
保持h=10、ε=0.8，但令空气300K，辐射环境400K。一个温度在两者之间的表面，一边从热环境吸收净辐射，一边向冷空气对流放热。稳态时净功率总和为0，物体不必与空气同温。
消去共同面积后，稳态方程为10(Tₛ−300)+0.8σ(Tₛ⁴−400⁴)=0。左端在300K为负，在400K为正；其对Tₛ的导数为10+4×0.8σTₛ³，在正温度下严格为正，因此区间内恰有一个根。数值解约348.904614K，即75.754614°C，高于空气温度。
这个例子说明，表面或探头的温度可能同时受空气和周围壁面的影响。它并未实现真实温度计模型：支撑导热、真实表面光谱特性、流动与内部温度分布均未纳入。MIT的热电偶示例也通过两条传热通道平衡解释这种偏离，见文献3，但其参数与本文算例不同。
flowchart TD
A[准备使用辐射公式] --> B[先转换为绝对温度]
B --> C[核对灰表面与包围环境条件]
C --> D[先算精确净功率和方向]
D --> E{是否需要常系数近似}
E -->|需要| F[写参照温度与误差范围]
E -->|不需要| G[保留温度依赖关系]
自检3：空气与壁面温度不同，能否只保留一个“环境温度”？
不能一般地合并。对流使用空气温度Tₐ，辐射使用包围环境温度Tᵣ，两者进入不同函数。只有附加条件允许时，才能定义另一个明确的等效环境，不能直接把其中一个代入所有项。
自检4：合成参数下近似只差2%，是否说明真实表面测温也只差2%？
不说明。2%只比较同一灰体模型下的精确四次方差与平均温度近似，未包含发射率不确定、几何、参与介质、传感器导热或内部温差。模型内近似误差与模型对现实的偏差，应分别评估。
使用前的一张核算清单
先记录每个温度属于哪里、是否已转换成K；再写明辐射面积、包围条件与ε的来源。真实ε可能依赖波段、温度、方向和表面状态，本文固定ε是建模选择，不是从材料名字直接读取的常数。
随后分别列出发出与吸收，确认净功率符号；与对流并存时，把各通道按同一正方向相加。最后才讨论是否冻结等效系数，并把参照温度、温度区间及误差口径一起保留在笔记里。若转向稳态导热或热应力问题，还需补上内部温度场和边界耦合，不应直接复用单温度假设。
本篇内部验证用独立高精度四次方计算核对因式分解、平均温度误差、温度交换、面积缩放和边界；稳态根另检查区间与残差。它核对的是合成算术与选定方程，不测量真实辐射、材料参数或设备性能。
参考文献与适用位置
- MIT Unified Engineering：19.1 Ideal Radiators；19.2 Kirchhoff's Law and Real Bodies。 黑体与腔内辐照、灰体表面关系。用于本篇能量模型的基础条件；R1据这些条件直接列账。返回公式R1。
- NIST：2022 CODATA Fundamental Physical Constants Complete Listing。 完整常数表，Stefan-Boltzmann条目。σ的完整值由定义常数关系确定，表列5.670374419…×10⁻⁸W/(m²·K⁴)；本文使用截取数值，不将9位小数当作完整常数。
- MIT Unified Engineering：19.3 Radiation Heat Transfer Between Planar Surfaces。 灰平板与热电偶能量平衡。用于说明任意两灰表面不能直接套单表面黑体环境表达式，以及空气/壁面温度不同的影响；本文未复制其数值算例。
- NIST SP330，第2节。 开尔文与摄氏温度关系。核对273.15K偏移及温差单位间隔相同的条件。