上个月帮一家环保工程公司做焚烧炉炉膛改造的仿真评估调了三轮边界条件最后模拟得到的炉膛出口烟气温度和现场实测只差了不到30℃。这种事放在十年前我是不太敢想的但现在CFD燃烧仿真在固废焚烧领域的介入程度确实越来越深。焚烧炉燃烧仿真不是新课题行业里叫了十几年可真正能把一个炉子从几何建模到结果评估完整做下来、并且让结果经得起现场检验的工程师数量还是不多。这篇文章就把这类项目从头到尾拆一遍聊聊焚烧炉燃烧仿真的模型怎么选、边界条件怎么给、二次风配比怎么调、最容易踩的坑有哪几个。适合已经在用CFD做流体仿真、想转做带化学反应的燃烧仿真的人也适合给垃圾焚烧发电厂、危废焚烧线做技术支持的设计工程师。1. 焚烧炉燃烧仿真到底解决什么问题1.1 设计选型阶段的不可视风险焚烧炉是固废处理链条里最吃工况的核心设备它不像普通锅炉烧的是单一燃料垃圾的热值、水分、灰分随时在变。设计一个焚烧炉光靠经验公式和前人的图纸风险集中在三个地方炉膛温度分布是否均匀、烟气在850℃以上区域停留时间是否满足环保要求、炉排上方有无局部高温导致结焦。这三个问题几乎没有哪个仪表能在设计阶段直接测出来但燃烧仿真可以提前给出可视化的温度场和流场分布。举个例子一个日处理量750吨的炉排炉设计时的容积热负荷、机械热负荷都通过了校核但实际投运后左右两侧的温差总在80℃以上。问题往往不在耐材和炉排选型而在后拱的形状与二次风射流的角度一个是结构件一个是流体动力学问题。仿真在这个阶段的介入就是把不可见的流场、温度场变成可量化的云图让工程师在制造之前就能评估炉膛型线的合理性。1.2 运行调试阶段的操作边界确定焚烧厂投运后运行工况也常需要调整。比如夏天垃圾热值偏低、水分偏高一次风量要不要加大冬季掺烧污泥后炉膛温度下降明显二次风比例要不要同步调整。这些参数如果都在实体炉上试代价极高试一次就要烧掉几百吨垃圾还有环保数据超标的压力。仿真可以把整条调试曲线在计算机上先跑一遍圈定一次风、二次风、辅助燃烧器的安全操作区间。所以一个合格的焚烧炉燃烧仿真项目交付的绝不止几张漂亮的温度云图。它要回答的问题是“这个炉子在各种边界条件下怎么运行最稳”具体到炉膛出口温度、氧含量、CO和NOx排放的趋势判断这些都是运行人员的操作依据。2. 仿真模型选型烧得准不准模型组合说了算2.1 湍流模型工业炉膛里我会优先选Realizable k-epsilon焚烧炉炉膛内基本都是高雷诺数、强旋流或强射流主导的湍流流动。常见的湍流模型里标准k-epsilon在旋流修正上偏弱k-omega SST在近壁区表现好但整个炉膛全用SST算网格和计算成本会明显上升收敛稳定性也是一个问题。Realizable k-epsilon在旋转流、射流、混合层这几类典型炉膛流动里表现均衡对网格质量的要求也不像LES那样苛刻是工程仿真里的稳妥选择。如果你做的是带旋流燃烧器的小型危废焚烧炉旋流强度大雷诺应力各向异性明显那可以再考虑加一个雷诺应力模型RSM来复核但RSM对边界条件很敏感计算余量不足的话容易振荡。我的做法是先用Realizable k-epsilon算出全局流场再针对局部重点关注区域切出来做二次验证性价比最高。2.2 燃烧模型涡耗散模型是焚烧炉仿真的主力但不是万能焚烧炉内的垃圾燃烧严格来说包含水分蒸发、挥发分析出、挥发分气相燃烧、焦炭表面燃烧等多个子过程。工程级仿真不可能把每个垃圾颗粒的详细反应都建模出来成本和时间都不现实。主流方案是把固体废物的热解和挥发分析出简化成源项把气相空间里的燃烧用涡耗散模型EDM来处理。涡耗散模型的核心假设是“化学反应速率由湍流混合速率控制”反应物一到就烧所以计算稳定、收敛性好非常适合大尺度工业燃烧室。焚烧炉火焰正好是典型的湍流扩散火焰这比使用有限速率化学动力学模型更贴合实际。但要注意EDM模型不能精确预测着火延迟和局部熄火所以在点火阶段和低负荷工况下结果只能做趋势参考。如果你需要更精确的CO、NOx分布可以在EDM基础上叠加有限速率/涡耗散FR-EDM两种机制取小值。2.3 辐射模型DO模型是焚烧炉温度场精度的保障燃烧室烟气温度动辄上千摄氏度辐射传热占总换热量的比例可以超过八成。选错辐射模型温度场误差一两百度都很正常。焚烧炉仿真首选DO模型也就是离散坐标法。它能处理各种光学厚度、能考虑气灰介质的散射和吸收对炉内温度分布和各壁面热流的模拟能力最强。P1模型计算开销小但它在光学厚度很薄的区域误差大而且对局部辐射源的各向异性描述不够S2S模型适合封闭腔内纯表面辐射不适用于含有吸收性介质的燃烧室。有些同事为了省迭代时间用P1结果炉膛出口温度偏差偏大后面对标现场数据时非常被动。DO模型虽然内存开销高一点以现在的计算能力完全可以接受。2.4 污染物生成模型NOx分析要看是定性还是定量焚烧炉氮氧化物来源包括热力型、快速型和燃料型三类。热力型NOx取决于高温区温度和氧气浓度在焚烧炉的富氧高温区占主导快速型在低温富燃料区重要燃料型则来自垃圾中的含氮组分机理复杂。Fluent自带的NOx模型可以耦合或者后处理方式计算热力型和快速型。工程上我常用后处理方式先收敛温度场和组分场再关掉其他方程单独算NOx生成速率速度快对趋势判断足够。但垃圾中的燃料氮含量波动大如果客户要求精确的NOx绝对值就得考虑燃料型NOx的用户自定义方程或专门的耦合模拟否则结果偏差会很明显。3. 实操过程一个完整焚烧炉仿真项目的核心环节3.1 几何处理把一台炉子变成能画网格的流体域焚毁炉几何建模的第一步是拿到炉膛结构图。不管原始模型是二维图纸还是三维装配体仿真需要的是烟气流动的流体域所以要对结构做简化处理。炉排以炉排缝的形式简化成多孔进气面挡火拱按实际形状保留过热器区域如果只模拟炉膛段就用一个交界面截断。我第一次做这类项目时踩过一个坑炉排下的风箱被我从计算域里删掉了但风箱出口风速对炉膛下部流场影响很大结果每次算出来炉排上方流速都偏低。后来我把风箱和炉排作为一个整体简化成带阻力特性的多孔介质入口问题就解决了。简化几何的原则是“保留对流动和传热有影响的特征去掉小圆角小倒角”比如看火孔、检修门尺寸远小于网格尺度留下只会增加网格量和收敛难度。3.2 网格划分多面体网格和边界层的配合焚烧炉几何复杂全六面体网格几乎不可能我优先使用多面体网格。多面体网格相比四面体网格在相同计算精度下数量能少一半以上收敛性也好。炉排上方近壁区域和二次风喷口处需要局部加密喷口射流核心区的网格尺寸要足够细才能捕捉到射流穿透和卷吸。壁面区域要不要画边界层取决于你关心不关心壁面换热细节。如果主要分析炉膛内部流场和温度场可以用标准壁面函数第一层网格的高度控制在y在30到100之间即可不需要刻意粘性底层。但如果你要计算水冷壁的热流密度那就得在壁面生成棱柱层保证y在1附近。网格无关性验证至少做三套从粗到细对比炉膛出口平均温度和壁面平均热流连续两套的偏差小于5%才算合格。偷懒不做网格无关性验证的仿真结果后面审查时容易被一票否决。3.3 边界条件与参数设定量级算错全盘皆输边界条件不是随意填数字最核心的是入口风量和温度这些要根据垃圾处理量和元素分析来算。以一个处理量1000吨/天的炉排炉为例垃圾收到基热值按10 MJ/kg估算总输入热量就大约是11574 kW。按过量空气系数1.8左右理论空气量可以按经验值每千克垃圾3.5标准立方米估算实际烟气量总质量流量大约在86 kg/s的量级。一次风和二次风的比例国内炉排炉常见设置为一次风占总风量的6到7成二次风占3到4成。一次风从炉排下方穿过垃圾层风温通常预热到120到180℃二次风在炉膛喉部喷入速度要达到20到40米每秒保证穿透卷吸效果。风温、风量、二次风速度这三个参数直接决定了炉膛内的主燃区位置和烟气扰动程度。边界条件设置还有一个容易忽略的细节入口的湍流强度和水利直径。风管入口的湍流强度按经验给5%到10%即可水力直径按风管截面折算。如果入口边界给得太粗糙射流角度和扩展角都会失真对炉膛中部的温度分布影响很大。3.4 求解器设置与收敛控制让残差掉下去是门手艺焚烧炉燃烧仿真运行时压力基求解器搭配SIMPLE算法是主流配置适合低速不可压缩到可压缩的工程问题。动量方程和湍流方程的离散格式尽量用二阶迎风能量和组分方程同样用高阶格式可以降低数值耗散。亚松弛因子要在计算过程中动态调整压力取0.3动量取0.7湍动能和耗散率取0.4到0.5能量取0.9组分取0.6到0.8。如果残差振荡先把能量和组分的亚松弛调到0.7附近等流场稳定后再恢复。点火策略我强烈建议分段操作。第一次启动直接全开燃烧很容易因为初场温度太低导致反应不启动、温度场发散。我是先给炉膛赋一个1300K左右的均匀初场并把能耗方程的亚松弛调低先算几百步建立基本流场再恢复边界条件和标准松弛。这样做下来炉膛温度起得快残差曲线也平顺很多。收敛判据除了残差更要看整炉温度是否稳定以及炉膛出口的氧含量、CO浓度是否不再单调变化。我一般迭代到3000步时看一眼出口温度历史曲线稳定了才算收敛。3.5 后处理与结果验证用数据而不是颜色判断方案后处理用温度云图、流线图和氧浓度分布切片但在项目交付时我会把仿真数据和现场设计值放到一张表里对比。重点看炉膛出口烟气温度是否在850到1100℃范围内炉膛中部的氧含量是否在6%到10%之间烟气在高温区的停留时间是否满足2秒以上。这些是焚烧炉设计和环保达标的核心考核指标仿真结果要能体现这些关键数据。如果一个方案在仿真中炉膛出口温度只有780℃那无论流场云图多漂亮环保排放已经挂在悬崖边了。这时需要用流线工具检查是否存在短路流烟气沿左右侧墙直接涌向后墙出口绕过中上部的主燃烧区于是停留时间严重缩水。这类问题在云图上看不直观但拉出流线后一目了然。所以后处理不要只截图而要针对关键指标做量化分析。4. 常见问题与排查技巧实录4.1 残差不降、温度发散的火爆问题燃烧仿真最头疼的就是算着算着温度飞了。温度发散的最常见原因有三个初场温度给得太低导致点火失败、局部网格质量差导致能量方程数值振荡、亚松弛因子给定不合理。点火失败的典型表现是氧含量呆在高位、温度从初场一路回落接近300K。这时要果断检查能量方程的亚松弛把它调到0.5以下把初场温度区间拉高并在炉膛中部做一个临时的高温点火楔形区算到火着起来之后再移除。另一种麻烦是出口回流导致负温。烟气出口如果位置紧邻结构突变区回流会带进冷空气再混合热烟气时就容易出现温度反插值。解决方法是加一段虚拟出口延长段让回流发生在计算域外很多情况加个一米长的扩展段就稳定了。4.2 二次风射流穿透度不足、混合质量差二次风的作用是强化炉膛内烟气扰动、延长烟气流程、补充氧量。射流穿透力不够的表现是烟气流线在上部炉膛几乎不存在横向卷吸温度场分布呈现出明显的高低温分层。遇到这种情况首先检查二次风边界条件中速度是否按实际设计风速输入而不是只给了流量。二次风喷口数量多、直径小单独每个口的速度要按实际几何面积折算只给总流量容易把速度算小射流自然穿不透。然后看喷口位置与后拱的配合。后拱的存在是为了把前端烟气导向中部燃烧区二次风若安装在后拱下方射流方向要向上倾斜角度常在30度到45度之间。角度偏平射流会直接撞向前墙卷吸区域变小。这一块往往是炉膛改造时最容易通过仿真验证后修正的地方。4.3 仿真结果与现场实测数据偏差大的排查思路仿真值和现场值的偏差通常集中在两类。一类是测点本身代表性不足现场的高温热电偶经常因为辐射误差比真实烟气温度高或者套管积灰导致响应迟钝拿这种数据对标仿真越对越慌。建议取现场数据时优先用网格测点或在线分析仪的氧含量数据来对比流场趋势温度数据作为参考带而不是硬点。另一类偏差出现在边界条件输入上。焚烧炉实际运行中垃圾热值随时漂移如果仿真全工况使用固定热值偏差是必然的。实操中我会做三个典型工况高温高负荷工况、平均设计工况、低热值最大连续工况。用这三个Case覆盖运行区间和现场数据取交集才能说明模型趋势可靠。4.4 常见异常的综合速查表问题现象可能原因优先排查方向温度场发散残差飞升初场温度过低点火失败提高初场温度、降低能量亚松弛出口烟气温度偏低但仍在下降燃烧不充分一次风过大核查空气总量与过量空气系数CO浓度异常升高二次风混合不足尾部短路检查二次风射流速度与喷口角度氧浓度分布严重偏置入口风量分配不均核查各风箱入口流量设定NOx预测偏高温度场偏热辐射偏弱检查辐射模型是否选用DO多工况收敛困难网格存在局部负体积检查翻面网格并重建局部区域炉内负压边界回流出口截面短产生反向流加设出口延长段4.5 搭建仿真模板把重复劳动变成标准化焚烧炉仿真项目做多了以后重要的事情是把边界条件参数做成模板而不是每次都从头点一遍鼠标。Fluent的Case文件可以保存液相、气相、壁面所有边界参数的命名规则我通常在第一个收敛的Case基础上把垃圾热值、空气量、二次风速度做成工作簿参数后续新工况只需改参数再重新迭代就行。一个炉型的Case模板搭配Excel参数表能把仿真周期从两周压到三天而且参数变更记录清清楚楚评审时能解释每个边界值怎么来的。用参数化脚本驱动的话还可以用Fluent的TUI命令把收敛判据和监控点都写进脚本批量跑工况时省得盯屏幕。这类工作的大量重复点在边界条件管理和数据提取花点时间做一套自己的工具链后面项目会越做越顺。5. 燃烧仿真在焚烧炉项目里的经验沉淀5.1 先搞懂工艺再动软件燃烧仿真最看不惯“软件操作很强、锅炉工艺不熟”的做法。模型选型、边界条件、收敛判断本质上全是工艺问题。如果不知道垃圾焚烧“850℃、2秒”的要求不知道一次风从炉排下方穿过的意义不知道怎么通过后拱和二次风组织烟气流程残差收敛得再漂亮也判断不了方案好坏。做焚烧炉仿真前花半天时间把炉排炉和流化床炉的工艺书本翻一遍比调试三天软件都有用。5.2 简化模型要有底线把复杂焚烧过程简化成EDM燃烧加DO辐射这是工程权衡是允许的。但简化要有底线垃圾热值、水分不能随意拍脑袋至少需要拿到设计院提资的元素分析数据或参考典型生活垃圾元素分析值。边界条件里每一个数据都要记录来源算完之后可以和现场运行报表对照。一次仿真项目如果连垃圾元素分析都拿不到只靠默认烟气组分估算出来的温度场只能算演示级结果不能支撑改造决策。5.3 工况覆盖比单点精度更重要一个焚烧炉项目的成败不在于某个工况算得准而在所有可能工况下都不会错得离谱。燃烧仿真不能只算一个设计工况就算完高负荷、低负荷、低热值垃圾三个典型工况至少都要覆盖。过程优化方案更是要避免把温度场调平以后轻微变工况立刻出现局部超温或温度不足。好的方案要有余量仿真在这个环节扮演的角色就是提前把余量的边界试出来。最后再说一个个人的习惯每个焚烧炉仿真项目结束后我会把“边界条件表、网格信息、收敛曲线、结果对比、结论和建议”打包成一份完整的计算报告存起来。下次再碰到同类型炉子直接在这个模板上改参数效率和稳定性都会好很多。燃烧仿真这个行当经验就是靠一个炉子一个炉子算出来的没有捷径。