能量方程的账本里,辐射是最容易被跳过、也最不该被跳过的路径:它以温度四次方的速度增长,高温场景一夜之间从配角变主角。本节先给一个"要不要算辐射"的快速判断法,再过一遍 Fluent 的五个辐射模型,最后用炉膛案例把选择过程走一遍。
第一步算黑体辐射量级:单位面积辐射热流等于发射率乘 Stefan-Boltzmann 常数乘绝对温度的四次方。八百摄氏度的表面,黑体辐射热流约每平方米七十五千瓦——同条件下对流换热往往只有几十千瓦,辐射已与对流平起平坐。第二步看介质:干燥空气基本不吸收辐射(透空腔体间辐射),含二氧化碳、水蒸气的燃烧产物吸收显著(介质参与辐射),颗粒群散射加剧复杂性。第三步看温差与几何视角因子:低温差大视角的场景(室温设备间)辐射可并入综合换热系数;高温差强视角(炉墙对面是工件)必须单独建模。三步走完,"算不算"有了量化答案,写进报告就是一条证据。
| 模型 | 原理 | 适用 | 忌讳 |
|---|---|---|---|
| Rosseland | 光学厚扩散近似 | 光学厚度远大于 1 的浓介质 | 光学薄时严重失真 |
| P1 | 球谐一阶展开 | 光学厚度中到大的炉膛、燃烧室 | 强方向性辐射(束流) |
| DO 离散坐标 | 按立体角逐方向解辐射输运 | 全景:透空、半透明、参与介质、散射 | 立体角加密带来计算量 |
| DTRM | 特征线追踪 | 简单几何透空辐射 | 复杂角系数与散射 |
| S2S | 面到面视角因子 | 干燥空气透空腔体,介质不参与 | 不能用于参与性介质 |
选择逻辑链条其实很短:介质不参与且几何不复杂,S2S 或 DTRM;介质参与或半透明界面在案,DO 是通用答案;光学特别厚且只需要能量源项,P1 或 Rosseland 省算力。DO 模型以通用性成为现代默认候选,代价是角分辨率参数(每像素立体角划分数)要按精度与算力平衡。
背景。某玻璃退火炉改造,需要预测玻璃板内温度梯度以控制残余应力。玻璃在红外波段半透明,辐射既发生在炉内气体与板面之间,也穿透玻璃板体内部——这是辐射仿真里最刁钻的一类。
操作。能量方程开启,辐射模型选 DO:它支持半透明介质(半透明界面与介质内吸收系数按波段设定)。把玻璃的吸收系数按红外波段分档输入,炉墙按发射率零点八设定,加热元件用辐射源面表达。与全 opaque 壁面、P1 模型两套简化方案做横向对比。
结果。DO 半透明方案预测的板内厚度方向温差与在线红外实测吻合良好;opaque 简化把板体当不透明体,板心温度被高估近三十摄氏度,据此设计的退火曲线会导致应力超标;P1 方案趋势对但细节糊,能量源项尚可用于初算。
解读。这个案例展示了辐射模型选择的另一维度——不光选模型,还要选物性表达:同一 DO 模型,把玻璃当不透明还是半透明,结论天差地别。辐射参与介质的物性(吸收系数、散射系数、波段划分)与模型本身同等重要,物性无出处,模型再对也白算。
变式。燃烧炉膛里介质物性来自烟气组分(二氧化碳与水蒸气的加权灰气体模型 WSGGM),Fluent 支持在组分输运求解时自动估算;含颗粒的煤粉炉还要叠加颗粒散射,此时 DO 几乎是唯一现实选项。
DO 模型的角分辨率参数设多少。透空腔体与强方向性问题从较低分辨率起步做收敛性检查,参与介质与半透明界面的精度随分辨率单调改善——它的正确设置方法只有一条:敏感性扫描,而不是抄参数。
发射率取值哪里最容易错。两处:氧化与抛光表面的发射率差数倍(金属表面状态敏感);灰体假设在宽波段强选择性介质(玻璃、塑料膜)失真——后者回到波段物性而不是换个发射率数字。
辐射算出来能量不平衡先查哪。先查参与介质的吸收系数量纲与波段定义,再查半透明边界的透射反射设置,最后查辐射与能量方程的耦合频率——三类错误合计覆盖九成不平衡。
辐射算例高频异常对照:
| 异常现象 | 多半原因 | 处置 |
|---|---|---|
| 温度系统性偏低 | 辐射散热被低估 | 核对发射率与视角配置 |
| 局部热点过锐 | 角分辨率不足 | 加密立体角并复算 |
| 迭代早期发散 | 辐射与能量耦合过强 | 降耦合频率或欠松弛 |
| 视角因子明显不对 | 遮挡计算配置错 | 核对遮挡面与几何精度 |
把辐射与前面章节的接口再理一遍,因为它的耦合方式决定收敛策略:辐射通过源项挂进能量方程,主温度场变化反过来改变辐射源,两者迭代耦合。松耦合(辐射每若干步更新一次)省算力但拖慢最终收敛,紧耦合每步同步代价高。工程节奏是:初算用松耦合快速逼近,收尾阶段切紧耦合把能量账轧平。另外两条纪律:辐射算例的能量残差判据应比纯流动更严(辐射源项对温度四次方敏感,残差标准要相应收紧);含辐射的确认对表优先选壁面热流密度而不是气温——热流对辐射配置的辨识度远高于温度,配错一眼见底。
关于多表面封闭腔的工程速算再补一笔:两个大平面之间的净辐射换热有经典的封闭腔公式,系数只取决于两表面的发射率——这个公式是辐射仿真的"手算锚点"。建模完成后,取两块代表性表面把仿真热流与手算公式对比,量级吻合说明基础配置(发射率、面积、视角)没有低级错误,再谈复杂几何与介质效应。这个小动作一分钟完成,却能拦住"介质模型很高级、基础发射率却填错"的尴尬组合,是辐射任务验收单上性价比最高的一行。
最后把辐射任务的"何时可以降级"说透,因为降级决策同样是专业能力:介质不参与、几何不复杂、温度低于五百开尔文三个条件同时满足时,辐射对结论的贡献可以并入综合换热系数,省下的算力与工时相当可观——但降级声明要进报告,且降级的适用范围(温度、几何)要写明。反过来,燃烧室、玻璃窑、太阳能接收器这类场景,任何降级都是拿结论冒险。判断的本质还是那三步快判加一条纪律:降级是有意识的工程决定,而不是默认的疏忽。把降级与升级都写进验收单的模型声明栏,这份诚实在评审时的价值会超出你的预期。