1. 辐射传热燃烧仿真中最容易忽略的主角做燃烧仿真的人都知道一个矛盾大多数人入门的第一个燃烧算例火焰温度总是偏高。不管湍流模型怎么调、化学反应机理怎么换炉膛出口温度就是压不下来。直到某天某位前辈淡淡说了一句“你开辐射了吗”你才猛然发现——你的燃烧器在纯真黑暗中燃烧热量只靠流动和对流带走壁面就像一块完美的隔热海绵什么东西都吸收不了。这个场景我见过太多次了。尤其是在导师项目组里做模拟项目X的燃烧室换热分析时不开启辐射传热模型和开启之后壁面热流分布差异可以达到30%到50%。换句话说燃烧仿真里如果漏掉了辐射传热你做出来的不是燃烧室只是一个“绝热搅拌反应器”。辐射传热模型在燃烧仿真中承担的任务简单说就是计算高温烟气、火焰碳烟、颗粒相等通过电磁波辐射方式向壁面或低温区域传递的热量。燃烧反应温度动辄1500K到2000K在这个温度区间辐射换热占总换热量的比例通常超过70%。它不是锦上添花的细节而是主导物理过程。这篇内容我想系统梳理一下燃烧仿真中辐射传热模型的原理、选型、实操设置和排查经验所有思路都基于我在几个燃烧室、锅炉和烧嘴仿真项目里的实际迭代过程希望能给正在做燃烧仿真、尤其是卡在温度偏高这个问题上的朋友一点启发。2. 燃烧仿真里为什么必须处理辐射先搞懂物理本质再谈模型2.1 燃烧产物为什么能“发光发热不靠摸”先聊一个基础问题气体凭什么能辐射热量我们理解对流换热靠分子碰撞理解导热靠分子微观扩散辐射完全不是这个逻辑——它靠的是电磁波。分子内部能级跃迁时释放或吸收光子高温烟气里的二氧化碳、水蒸气、碳烟颗粒就是靠这个过程和壁面交换能量。这里有个和生活经验强相关的点你以为火焰的“亮”只是化学发光其实火焰发光主要来源于碳烟颗粒的连续辐射接近黑体辐射谱。而二氧化碳和水蒸气的辐射则是有选择性的——只在特定波段吸收和发射。这就是为什么燃烧仿真里不能简单地把烟气当作一个“灰体”处理而必须考虑参与介质的光谱选择性。辐射传输方程Radiative Transfer Equation简称RTE是所有辐射模型的物理起点。它描述的是辐射强度沿着一条视线方向传播时被介质吸收而衰减、被介质发射而增强、被颗粒散射而改变方向的过程。写出来大概是dI(r,s) / ds -(ka ks) * I(r,s) ka * Ib(r) (ks / 4π) * ∫ I(r,s) * Φ(s,s) dΩ这里每一项的物理意义可以用一个类比来理解一束光在烟气中往前走就像一个人在黑暗的走廊里举着手电筒——走廊墙体在吸收光吸收项走廊里有人在不断点亮新的灯发射项还会有人拿着镜子改变光线方向散射项。燃烧仿真要做的事情就是把这个辐射强度场在整个计算域内求解出来然后换算成辐射热流、辐射源项耦合进能量方程。2.2 七成热量靠辐射这个数字是怎么来的很多刚接触燃烧仿真的朋友会问辐射有那么重要吗我开不开辐射温度也没差多少啊——这种情况多半发生在小火焰、低温度或壁面温度接近火焰温度的工况。一旦温度超过1200K辐射换热的占比会急剧上升。道理可以从斯特藩-玻尔兹曼定律看出来辐射热流正比于温度的四次方。温度从1000K升到2000K辐射换热能力提升16倍。对流换热大致正比于温度的一次方或二次方还要看流动状态此消彼长之下高温区间必然是辐射主导。我曾经做过一个简化验证同样的500 kW燃气燃烧器分别开辐射和不开放辐射。不开放辐射时出口烟气温度预报值比实测高了大约180K壁面热流几乎为零。而开启DO辐射模型配合加权灰气体模型之后出口烟温预报值落在实测误差带内壁面热流分布也与红外热像仪测到的趋势吻合。温差180K是什么概念它可以改变整个换热器的设计裕度可以直接决定后墙耐火材料选型对不对。所以不是“开不开都行”而是“不开就别谈准确预测”。3. 燃烧仿真主流辐射模型选型DO、P1和那些“看起来简单但藏坑”的方案3.1 模型家族梳理从六通量法到蒙特卡洛工程上做燃烧仿真辐射模型的选择范围其实很窄。早期火电锅炉设计里流行过区域法Zone Method和六通量法但这些方法在复杂几何和三维问题里局限性很大现在商业软件里基本退居二线。目前主流CFD软件里能选的主要是P1法、离散坐标法DODiscrete Ordinates和蒙特卡洛法MC。P1近似属于球谐函数展开的一阶近似求解的是一个关于入射辐射强度的扩散方程。计算量小对光学厚介质烟气层很厚、吸收很强预测不错。但对光学薄区域、各向异性散射、局部热源附近的表现比较差。DO模型把辐射传输方程按照有限个方向离散对每一个方向独立求解。计算量大但适用范围广从光学薄到光学厚都行是燃烧仿真的主力选项。蒙特卡洛法用大量光子追踪统计精度高但计算成本极高目前在工程燃烧仿真中用得比较少更多用于基准验证。从我的实际使用体会来说大部分燃烧室、锅炉、燃气轮机燃烧室仿真直接上DO模型即可。只有在计算资源非常紧张、且燃烧室空间尺度大、烟气流场充分发展、光学厚度较大的锅炉炉膛这类场景P1才是性价比合理的选择。3.2 为什么我推荐DO模型加WSGGM一个关键组合选DO模型只是第一步真正决定辐射计算准不准的是介质辐射特性的处理方式。燃烧烟气里的二氧化碳和水蒸气都是“非灰气体”——它们在特定波段的吸收能力很强在其他波段又几乎是透明的。如果用灰体假设统一处理误差会相当可观。工程上解决这个问题最成熟的方法是加权灰气体模型WSGGMWeighted Sum of Gray Gases。它把非灰气体的辐射特性用若干个“假想灰气体”加权叠加来近似每个灰气体有自己的吸收系数和权重权重随温度变化。商业CFD软件里内置了常用的WSGGM系数一般适用于燃烧产物中二氧化碳与水蒸气的分压比在一定范围内的情况。我踩过一个很典型的坑在某个富氧燃烧仿真里烟气中水蒸气分压特别高远超默认WSGGM系数的适用范围。用默认参数算出来的壁面辐射热流明显不合理——偏低。后来改用预设的富氧燃烧WSGGM参数组合某些软件里有针对高水蒸气分压的系数集结果才回到合理区间。所以当你做富氧燃烧、高湿度烟气或特定燃料燃烧时一定要去核查WSGGM系数的适用分压范围别做默认党。3.3 DO模型离散参数怎么定方向数不是越多越好DO模型里有Theta Divisions和Phi Divisions两个参数分别控制极角和方位角的离散数量。默认值通常是几乘几但好多人不知道这个值会显著影响计算量和精度。以我常用的三维算例来说从3×3提升到5×5计算量提升非常明显辐射热流结果却只变化了百分之几。而某些细长体几何中如果角度离散太粗会出现“射线效应”——辐射能量沿着离散方向“漏出去”壁面热流出现非物理的条带状分布。解决办法不是盲目增加方向数而是先做一次角度离散无关性验证。我一般会选择3×3、5×5、8×8各算一遍比较关键壁面热流的变化如果5×5和8×8差异小于2%就选用5×5。这个习惯能帮你在精度和成本之间找到平衡点。另外DO模型与壁面发射率配合时还有一个容易被忽略的设置壁面角离散数。它影响壁面附近辐射强度的计算精度默认值是多少就用多少但如果是半透明介质边界或反射较强的壁面需要考虑增加。这部分我在下面的实操章节会展开说。4. 实操案例某小型燃气燃烧室的辐射传热仿真全流程4.1 几何简化与网格策略辐射计算对网格的隐性要求我用一个真实做过的模拟项目X来走一遍流程。这是一个中试规模的圆柱形燃烧室直径0.6米长度1.5米甲烷-空气旋流扩散燃烧设计热功率500千瓦。仿真目标很明确获取壁面辐射热流分布校验燃烧器喷口结构改进方案。几何建模上考虑到旋流器和燃烧室是轴对称结构我选择了二维轴对称模型。这里要先提醒一句辐射模型在二维轴对称下是完全支持的但后处理时壁面热流是按单位面积计算的你需要手动把轴向热流分布转换成总面积热流。如果换热结构有明显的三维特征比如多个燃烧器交错布置、水冷壁鳍片那就老实做三维不要贪二维的便宜。网格划分上很多人觉得辐射计算对网格要求不高其实不完全对。辐射求解本身对网格没有那么敏感但如果辐射模型和燃烧模型耦合求解火焰面附近温度梯度极大网格太粗会导致火焰面分辨率不足温度场失真然后辐射热流自然也不准确。我的经验是火焰面附近最小网格尺寸控制在2毫米左右燃烧室主体网格做渐变拉伸。总网格量在二维下大约5万左右三维下可能需要300万到500万取决于几何复杂度。最重要的是要保证壁面第一层网格的y满足湍流壁面函数的要求否则近壁温度不对辐射热流的边界条件就错了。4.2 燃烧与湍流模型设置辐射不是孤岛它会和燃烧互相反馈在解算之前先把燃烧仿真本身的模型选好。我选的是Realizable k-ε湍流模型配合EDCEddy Dissipation Concept燃烧模型。EDC能考虑有限速率化学反应效应对甲烷空气扩散燃烧的预报比简单的涡耗散模型EDM更可靠尤其是对熄火极限和污染物排放预测有优势。如果你只是看热流分布和温度场EDM也够用但如果项目中涉及CO等中间产物建议用EDC并配合一套简化的甲烷反应机理。这里有一个重要的耦合逻辑燃烧反应放热决定温度场温度场通过辐射模型决定辐射源项辐射源项反过来修正能量方程最终改变温度场。所以辐射计算不是独立的一次性计算而是在每次迭代中与流动方程、能量方程同步推进。开启辐射模型后收敛速度和稳定性都会受到明显影响。这也是后面讲收敛策略时要重点说的事。材料物性设置上除了常规的密度、比热、导热系数关键是参与辐射介质的属性。Fluent中DO模型配合WSGGM之后气相的吸收系数会自动根据CO2和H2O的分压与温度计算。这里注意一下“吸收系数”的呈现形式WSGGM是非灰气体模型它计算得到的是基于灰度加权后的等效吸收系数你不需要手动给烟气填一个常数。相反如果你用的是灰气体假设就需要手动估算一个平均吸收系数这个过程的主观性很大我试过给同一个炉膛填0.1和0.5的常数吸收系数壁面热流差了将近一倍。4.3 边界条件与材料参数壁面发射率比你想的更重要辐射传热的边界条件主要涉及壁面发射率Emissivity和壁面温度或热流条件。很多仿真案例里大家习惯把壁面当作绝热或定温边界但燃烧室壁面实际是通过耐火材料和冷却水套带走热量的因此通常应设置对流换热边界条件或者给定壁面热流。壁面发射率这个数看似简单取值影响极大。大多数耐火材料在高温下的发射率在0.8到0.95之间粗糙的耐火砖偏上限光滑的金属壁面或水冷壁偏低可能在0.4到0.7。我在某个案例里把燃烧室内壁发射率从默认的1黑体改为实测的0.85辐射热流下降约15%。这个量级足以影响水冷壁设计。入口边界条件甲烷入口设为质量流量进口质量流量0.0093 kg/s对应500 kW甲烷低热值约50 MJ/kg温度300K。旋流空气入口质量流量约0.10 kg/s过量空气系数取1.1左右温度600K模拟实际预热空气工况。出口用压力出口表压为0回流温度设为环境温度。这些参数直接影响燃烧温度进而影响辐射——入口条件不准后面全都是白算。4.4 求解控制与收敛策略开着辐射模型跑不起来怎么办把辐射模型和燃烧模型耦合在同一个算例里最常见的现象就是发散或收敛速度极慢。这里我分享一套屡试不爽的操作流程第一步先关闭辐射模型只算冷态流场或者先用一个简单的点燃温度分布把流动和燃烧迭代到初步收敛。这个阶段用一阶迎风格式亚松弛系数保守一点让残差下降两个量级以上。第二步打开辐射模型把能量方程的亚松弛因子从默认值略微调低比如从1.0降到0.95或0.9。同时DO辐射方程本身的亚松弛因子保持默认或稍降低。继续迭代这时残差通常会先上升尤其是能量残差别慌让它跑几百步再看看趋势。第三步当残差重新掉头下降并稳定后把动量、能量、湍流方程的离散格式从一阶迎风切换到二阶迎风或QUICK辐射方程维持一阶即可。再次迭代到残差收敛同时监测壁面辐射热流和出口温度是否基本不变。关于残差收敛标准能量方程残差我一般要求低于1e-6辐射方程低于1e-4至1e-5连续性低于1e-4。但仅看残差不够还必须看关键物理量——出口烟气温度、壁面总热流、火焰峰值温度——这些量在最后几百步迭代里不再发生趋势性变化才叫真收敛。还有一个小技巧开辐射后如果发散多半是能量方程亚松弛太大或者初始温度场不合理。我的做法是先给定燃烧室内部一个合理的初始温度场比如1200K而不是默认的300K这相当于给点火过程一个热身能显著降低前期的不稳定性。5. 从结果里读出问题辐射热流数据的分析与验证5.1 壁面辐射热流怎么看曲线、云图与三个判断基准计算收敛之后第一步是提取壁面辐射热流分布。在轴对称模型里我会沿着燃烧室长度方向画一条轴向位置-热流密度曲线在三维模型里可以生成热流云图。判断结果是否合理我习惯用三个基准第一分布趋势。正常燃烧室的壁面热流峰值应当出现在火焰高温区附近也就是燃烧器喷口下游一段距离的位置而不是紧贴着喷口壁面或出现在出口附近。如果热流峰值跑到出口附近多半是燃烧区域位置预测不准要回头看燃烧模型和入口旋流强度。第二量级合理性。500千瓦热功率的燃烧室壁面平均热流一般在小几十千瓦每平方米量级峰值能达到平均值的两到三倍。如果你算出来的壁面热流达到几百千瓦每平方米只能说明辐射模型或边界条件有重大问题。此时先检查壁面发射率是不是被设成了1烟气吸收系数是不是被手动填了一个离谱的值。第三与实验或经验公式对比。这种对比是提升仿真可信度的关键一步。我们项目中用红外热像仪测了石英视窗外壁温度再通过一维导热反推内壁热流和仿真结果的偏差控制在15%以内这个误差水平对于燃烧工程是完全可以接受的。5.2 有没有辐射模型结果差异到底有多大过程记录一个直观的对比结果——同一个模拟项目X同一套网格和燃烧设置分别开启和关闭DO辐射模型。不开放辐射火焰峰值温度约1950K出口烟气温度约1250K壁面热流几乎为0。开启DOWSGGM辐射火焰峰值温度降到1800K左右出口烟气温度约1120K壁面辐射热流峰值约85 kW/m²。差异不只是温度数值。没有辐射时能量只能通过对流传递到壁面壁面附近温度边界层被拉得很高而实际物理过程中辐射可以直接把能量“跨过边界层”送到壁面导致壁面附近的温度梯度明显变化。也就是说辐射模型还影响了近壁对流换热的计算。如果你在某天需要做流固耦合换热分析不在意这个耦合逻辑后面传热计算会差得很远。6. 常见问题与排查技巧辐射传热仿真“翻车”集中营6.1 典型问题速查从发散到热流失真我做燃烧仿真这几年辐射相关的问题基本可以归成几类。下面用表格整理一下方便各位对照排查。症状可能原因排查与解决方案开辐射后能量残差骤升、发散亚松弛过大、初始温度场不合理降低能量亚松弛到0.8左右初始化一个合理的高温场先跑一阶格式壁面热流明显偏高壁面发射率被默认成了1吸收系数被填得过大核对壁面发射率耐火材料取0.8到0.9确认气态吸收系数是否由WSGGM自动计算壁面热流出现条带状非物理分布DO模型角度离散数不足射线效应增加Theta/Phi Divisions做方向无关性验证燃烧室出口温度比实测高约200K没有开启辐射模型开启DO模型并勾选WSGGM检查能量方程收敛富氧燃烧或高水蒸气工况辐射热流偏低默认WSGGM系数不适用查找更适用于高水蒸气分压的WSGGM系数集流场计算稳定但温度场飘移不定辐射与燃烧耦合后收敛判据不足监测出口温度和壁面热流的动态变化不要只看残差火焰熄灭或温度过低燃烧模型点火不稳定EDC在高拉伸率区chemistry偏慢检查入口温度、点火源设置或改用稳态火焰面模型作对比壁面温度分布与热流分布矛盾壁面边界条件设置不正确区分固定温度、固定热流和外部对流三类边界条件结合物理实际选择6.2 DO模型和P1模型选错导致的经典失败案例还有一类问题不是参数设错而是模型本身就选错了。我曾做过一个高温贫燃燃烧室案例介质光学厚度比较薄散射也很弱。当时图省事选了P1模型结果火焰下游区域的壁面局部热流出现了“凹陷”——在物理上完全说不通。P1模型在光学薄区域和辐射强度定向性强的位置会产生较大误差它本质上是一个扩散近似更适合那些辐射可以近似各向同性传播的场景。换成DO模型之后同样条件算出来的热流分布非常平滑火焰附近的峰值和下游的衰减都自然合理。所以你在燃烧仿真里如果拿不准选哪个模型我的建议永远是默认选DO除非你明确知道这个算例的光学厚度很大、计算资源又紧张到必须压缩计算量。这个选择逻辑能避免90%以上的模型适用性问题。6.3 壁面发射率取值现场实测与软件默认值的博弈最后再聊聊壁面发射率。很多软件的默认边界条件里发射率默认是1但工程实体几乎不可能是黑体。耐火浇注料、莫来石砖、碳化硅材料、金属水冷壁发射率各不相同。我从现场经验里学到的做法是有条件就用高温红外测温仪测量实际壁面温度再结合热电偶数据反算出发射率的有效范围没条件时至少查阅材料手册得到合理的取值范围然后取中间值并做敏感性分析——把发射率从0.75调到0.9看热流和温度场的变化幅度。如果变化在10%以内说明你的计算对该参数不敏感结论比较可靠如果变化超过30%那这个参数就是你项目里的重大不确定性来源需要在报告里如实说明。我个人在实际做燃烧仿真时的体会是辐射传热模型是整个燃烧仿真链条里“投入产出比”最高的一个环节。它不像详细化学反应机理那样需要消耗大量计算资源也不像湍流模型那样有几十个变种让你纠结它只需要你理解物理、选对模型、认真设边界条件。但就是这么简单的一件事很多仿真报告里就是不愿意把那个勾选打开。如果你现在手头的燃烧算例温度偏高、热流离谱先别急着换湍流模型好好审视一下你的辐射设置——那行小小的“Radiation: Discrete Ordinates”勾选框可能就是你和真实物理世界之间的桥梁。这个内容后续做流固耦合、热应力分析时辐射热流数据还能继续延续使用所以这一关值得多花心思。