参数设置
假设 A_1 = A_2 = 1 m²、角系数 F_12 = 1(无限大平行平板)。
辐射热阻网络
左节点=σT_1^4 / 右节点=σT_2^4 / 中间为表面热阻 R_1, R_2 与空间热阻 R_12 串联
理论与主要公式
把 σT^4 视为驱动势,将表面热阻和空间热阻串联,与电路完全相同地计算辐射热流。
表面热阻(ε < 1 造成的反射损失):
$$R_s = \frac{1-\varepsilon}{A\,\varepsilon}$$
空间热阻(由角系数 F_12 决定的几何损失):
$$R_{12} = \frac{1}{A_1\,F_{12}}$$
两表面之间的辐射热流。σ = 5.670×10⁻⁸ W/(m²·K⁴):
$$Q_{12} = \frac{\sigma\,(T_1^4 - T_2^4)}{R_1 + R_{12} + R_2}$$
黑体极限(ε_1 = ε_2 = 1)和灰体效率:
$$Q_{BB} = \sigma A (T_1^4 - T_2^4),\quad \eta = \frac{Q_{12}}{Q_{BB}}$$
无限大平行平板(F_12 = 1, A_1 = A_2 = A)时简化为 η = 1/(1/ε_1 + 1/ε_2 - 1)。
什么是辐射网络模拟器
🙋
听说真空中也能传热,那就是辐射对吧。但是为什么要特意叫「网络」呢?
🎓
对,辐射是电磁波,所以真空中也能传播。简单说,把辐射换热的分析改写成与电路完全相同的形式,就能解得非常清晰。「电压」是 σT^4,「电阻」是表面热阻和空间热阻,「电流」就是热流 Q——这种对应关系本身就叫「辐射网络」。
🎓
表面热阻 R_s = (1−ε)/(Aε) 是因为表面不是黑体造成的「反射损失」。发射率 ε 越小、越像镜面的表面,辐射就越容易被反射回自己。空间热阻 R_12 = 1/(A_1·F_12) 是由角系数 F_12 决定的「几何上的不易到达」。在模拟器里把发射率滑块往下调,会看到表面热阻急剧上升、热流一下子下降。
🙋
温度的四次方影响这事经常听说。一动滑块 Q_12 就增长得特别快。
🎓
这就是斯特藩-玻尔兹曼定律 E_b = σT^4 的威力。把 T_1 从800K加到1600K,T^4 就变成16倍。同样是「800K温差」,800K vs 0K 与 1600K vs 800K 的辐射强度差好几个数量级。所以炉子和燃烧器这类高温装置里,辐射会比对流还要主导。
🙋
那灰体效率 η 只有66.7%,剩下的能量「丢」到哪里去了?
🎓
没有真的丢。在表面热阻处「被反射回自己表面」的那部分,结果上不贡献到净热流。所以保温瓶和MLI(多层隔热)反过来利用这一点。把 ε 极小的镀铝箔叠很多层,故意把表面热阻做得很大来阻挡辐射热。在模拟器里把 ε_1 = ε_2 = 0.05 试试,Q 会一下子掉到原来的1/10以下。
常见问题
由材质和表面状态决定。抛光铝为0.04~0.07,氧化铝为0.20~0.30,氧化钢为0.60~0.85,白漆为0.80~0.95,接近黑体的黑色涂料在0.95以上。同种材料随温度和氧化状态变化很大,所以设计时要使用实测值或可靠手册数据。波长依赖性强时,需要考虑光谱发射率。
是「发射率与波长无关」的近似,能大幅简化分析。对大多数工业材料,在中红外波段(典型工业温度300~2000K的辐射峰值附近)基本成立。但像玻璃、陶瓷这类光谱特性强的材料,或选择吸收涂层(太阳能集热器的选择面),需要按波长带分别处理。
随温度区间差异很大。室温附近自然对流的传热系数约5~10 W/(m²·K),辐射的等效系数也大致相当,二者势均力敌。500K以上辐射开始迅速主导,超过1000K后达到对流的几倍至十几倍。反之低温下导热和对流主导,辐射有时被忽略。
主要有三种方法。S2S(surface-to-surface)法预先计算角系数矩阵,用与本工具相同的网络求解,适合介质透明的场合,速度最快。DO(discrete ordinates)法是能处理介质吸收散射的有限体积方法,用于燃烧室等参与性介质。蒙特卡洛法用光线追踪,对复杂形状和镜面反射很强,常用于校核和高精度计算。
实际应用
炉子和燃烧器设计:炼钢炉、玻璃熔炉、锅炉炉膛中辐射占传热的60~90%。精确求解从燃烧器到炉壁、再从炉壁到被加热物的辐射网络,直接决定燃料消耗、温度均匀性和耐火材料寿命。CFD程序中的S2S和蒙特卡洛辐射模型,正是把这种辐射网络放大到大规模情形。
航天器热设计:真空中的人造卫星不存在对流,几乎全部排热都靠辐射。用MLI(多层隔热,发射率0.05的镀铝薄膜十几层)把设备包住,以放大表面热阻;通过散热面(高发射率)向太空背景排热——这种设计本身就是辐射网络分析。通过表面处理(白色涂装、镀银、OSR)的选择来控制热平衡。
建筑窗户和保温:Low-E(低发射率)玻璃在玻璃内表面镀氧化锡等薄膜,把发射率从0.85降到0.10左右,大幅减少室内外的辐射换热。中空玻璃中夹层的氩气抑制对流,Low-E镀膜抑制辐射,使热透过率(U值)降到单层玻璃的1/3以下。
红外热成像的校准:热成像把表面发射的辐射能量换算为温度,但严重依赖 ε。同样温度下,黑漆面(ε≈0.95)和抛光不锈钢面(ε≈0.10)的检测强度差近10倍,所以在测量对象上贴发射率胶带或喷已知发射率涂料来校准。本工具改变 ε 滑块时的变化,正对应这种误差源。
常见误解与注意事项
最常见的误解是把「越黑的物体越吸热所以越烫」简化得过头。可见光波段「黑色=高吸收率」确实成立,但传热中重要的是中红外波段(典型工业温度的辐射峰值)的发射率,并不与可见颜色一致。例如白色涂料反射可见光,但在中红外波段 ε ≈ 0.9,行为接近黑体。请记住模拟器的 ε 滑块要输入的是「目标温度区间的发射率」,不是可见颜色。
其次常见的是把角系数 F_12 误解为「只要看得到对方就是1」。F_12 是「面1发射的辐射中到达面2的比例」,只有像无限大平行平板或同心圆筒内→外这种完全包围的情况才等于1。一般凸形状下 F_12 < 1,剩余的 (1−F_12) 流向其他面或环境。本工具简化为无限大平行平板(F_12 = 1)所以网络解很简单,但实际机器中角系数的精确评估直接决定分析精度。
最后请注意,不能用温差或绝对温度来除表面热阻。辐射网络是「σT^4 为电位、Q 为电流」的欧姆定律类比,电位差是 σ(T_1^4 − T_2^4),不是 T_1 − T_2。如果照对流换热中「传热系数 × 温差」的感觉把 R 改成基于温差的形式,低温下看似对,高温下偏差会很大。辐射的本质是四次方定律,网络是严格地以这个四次方差为驱动力的线性电路——这一点必须牢记。