甲烷MILD燃烧中H2/NH3/蒸汽添加的空间热负荷管理与N2O感知多目标设计

《Energy Conversion and Management-X》：Spatial thermal load management and N 2O-aware multiobjective design of methane MILD combustion with hydrogen, ammonia and steam

　　

MILD燃烧依赖内部回流在显著反应发生前对反应物进行预热和稀释。由此形成的分布式放热可限制温度峰值和热力型NO生成，但其运行区间仍取决于混合、停留时间、壁面取热和燃料化学特性。氢气与氨以不同方式改变上述过程：氢气加速链式分支反应并减少替代燃料组分中的碳含量，但可能强化局部反应并增加壁面负荷；氨同样不含碳，但其氧化速率较慢且燃料氮引入了通往NO、NO₂和N₂O的竞争路径。因此，针对单一组分获得的结果不能直接移植至甲烷-氢气-氨混合体系。研究采用实验室尺度MILD炉膛的二维轴对称模型，以GRI-Mech 2.11为甲烷/氢气氧化骨架，按路径不重叠原则引入Stagni等（Stagni et al.）的N-H-O化学，形成含61组分、418个气相反应的合并机理，并以Pan等（Pan et al.）开发的简化机理进行边界匹配的化学与物性数据敏感性对比。以CO、CO₂强度、NOx和N₂O为四个环境目标，采用响应特异性径向基函数（RBF）代理模型、多目标粒子群优化（MOPSO）和逼近理想解排序法（TOPSIS）开展优化，并结合壁面热流、放热率（HRR）轴向质心与中心跨度、反应-壁面输运距离等空间指标进行热负荷评价。研究识别出碳-氮帕累托结构，量化了氢气、氨和蒸汽对壁面热输运与空间负荷集中的影响，并通过网格细化、机理替换、EDC时间尺度扰动和壁温变化评估了结果的稳健性。研究结果为含氨MILD燃烧的环境-热负荷协同筛选提供了框架，同时表明绝对氮排放与局部极值仍强烈依赖化学机理与热边界设定。

研究背景：MILD（Moderate or Intense Low-oxygen Dilution）燃烧通过内部回流在大量反应发生前对反应物进行预热与稀释，形成分布式放热，从而抑制温度峰值与热力型NO生成；但其运行区间受混合、停留时间、壁面取热与燃料化学共同制约。氢气加速链式分支并去除替代燃料中的碳，却可能强化局部反应、提高壁面负荷；氨同样无碳，但氧化缓慢且燃料氮引入通往NO、NO₂与N₂O的竞争路径，已有结果表明含氨计算对稀释、湍流模型与动力学模型高度敏感。蒸汽增加热容并改变自由基池，可降低温度与NO，同时使氮向其他产物迁移，因此N₂O响应必须显式解析。单一组分的结果无法直接外推至CH₄/H₂/NH₃混合体系；且仅凭出口排放不能确立热适应性，同一混合气可在降低全局污染物的同时将放热集中于小体积或提高峰值壁面热流。基于此，研究人员在同一实验室炉膛构型上重新表述环境响应，并引入壁面传热、反应局部化与轴向反应-壁面输运分析，开展空间热负荷管理与N₂O感知的多目标设计研究。研究目标包括：识别氢、氨与蒸汽添加的碳-氮帕累托结构；比较所选组成在壁面热输运与空间负荷集中方面的差异；通过单元与面分辨指标量化反应-壁面输运；并在网格细化、机理替换、EDC时间尺度扰动、壁温变化与决策归一化变化下评估所选结果。

关键技术方法：计算域为Rebola等（Rebola et al.）实验的实验室尺度MILD炉膛（腔体直径150 mm、长300 mm，燃料经4 mm中心喷嘴进入，氧化剂经同心环缝进入，下游设85 mm直径、150 mm长烟囱），采用二维轴对称、结构化网格（生产网格9 600个控制元，另设19 760与39 726网格用于三网格检验）；湍流采用标准k-ε模型配增强壁面处理，湍流-化学相互作用采用涡耗散概念（EDC）；辐射采用离散坐标（DO）法与灰气体加权和（WSGG）模型；化学为GRI-Mech 2.11与Stagni N-H-O化学按路径不重叠规则合并的61组分、418反应机理，并以Pan简化机理作独立敏感性对比；代理模型为响应特异性正则化高斯RBF，优化采用MOPSO，排序采用TOPSIS并进行权重与归一化敏感性检验；数据库含50个CFD案例（24个复用因子设计案例与26个新组成），其中45例用于建模、5例预先指定为留出验证，另以3个全新初始化CFD案例确认优化选择。

研究结果：（1）验证层级与数值一致性——Rebola甲烷数据用于评估甲烷参考场（210 mm处温度MAE/RMSE为74.56/93.41 K，CO为0.00114/0.00137干基vol.%，CO₂为0.571/0.584干基vol.%，NO为1.217/1.467干基ppm），三网格试验度量含氨V2点的离散敏感性，守恒报告检验数值闭合，留出与全新启动案例检验代理与优化流程，机理、EDC与壁温变化属模型敏感性检验；守恒与插值精度不等于实验精度。（2）修正CFD响应趋势——在零氨、零蒸汽下提高H₂通过促进CO氧化并置换甲烷降低CO与CO₂强度；20% H₂、零氨时蒸汽使CO出现浅极小值、显著降低NOx并略增N₂O；20% H₂、20%蒸汽下增加NH₃降低CO₂强度但使NOx与N₂O陡增，CO在较大氨添加时亦上升，修正NH₃逃逸相对较小，数据库因此存在真实冲突：碳强度最低的组成未必氮负荷最低。（3）代理模型性能——拟合对数域交叉验证R²为0.9807–0.99996，5个留出案例物理域MAPE为CO 1.145%、CO₂强度0.053%、NOx 7.299%、N₂O 8.340%，最大单点偏差为低氨组成下NOx的16.115%，与零氨边界附近陡梯度一致。（4）四目标帕累托行为——无约束档案含极低NOx与N₂O的宽无氨分支，以及以更高氮排放换取更低CO₂强度的含氨分支，H₂在低碳区大部分逼近上界，蒸汽调控CO氧化与两种氮产物间的局部平衡；施加NH₃≥5%约束后近零氮分支被移除，暴露含氨混合气内部权衡，较大氨添加降低碳强度但提高NOx与N₂O。（5）TOPSIS选择与权重敏感性——等权TOPSIS在无约束档案选择0% NH₃、20% H₂、约9.8% H₂O，等权、碳优先与平衡排放权重返回同一档案成员，氮优先权重转向0% NH₃、约19.6% H₂、20.0% H₂O，避氨与高氢偏好跨四种预设情景稳定而蒸汽份额依赖偏好；对数归一化下四种权重在含氨约束档案均选择5% NH₃、20% H₂、零蒸汽，其确认CO₂强度低于无约束折中但预测NOx超过1200 ppm；线性极差归一化将无约束等权选择改为约0.03% NH₃、20% H₂、零蒸汽，而约束等权选择不变，四种权重下对数与线性排序的Spearman相关系数在无约束档案为0.640–0.883、约束档案为0.977–0.998；甲烷参考缩放使无约束档案四种权重均选择约0% NH₃、19.55% H₂、19.98%蒸汽，约束档案等权选择约5% NH₃、1.88% H₂、零蒸汽，表明中间蒸汽选择依赖归一化方式。（6）优化后CFD确认——12个代理预测经全新启动CFD在15%以内复现，CO₂强度偏差不大于0.11%，CO、NOx与N₂O最大偏差分别为11.47%、13.28%与13.05%，V2的NOx误差位于NH₃≥5%约束边界并与留出试验中的低氨陡响应一致；V2最大温度为1722.7 K，V3引入5%蒸汽后降至1648.0 K，含氨案例N₂O集中于上游反应区而无氨V1场中几乎不存在。（7）氢富集：更高壁面热输运伴随更强局部化——零氨零蒸汽下H₂由0%增至20%使总壁面传热由4.603 kW升至5.256 kW（+14.2%），平均壁面热流同比例增加，最大壁面热流由60.47升至86.15 kW m^?2（+42.5%），平均温度升高91.8 K、标准差升高29.3 K；q99由1.548×10⁷升至2.624×10⁷ W/m³，CVHRR由1.97升至10.95，C007仅10个单元（占炉膛体积0.0066%）超过其最大值90%，单单元极值仅作警示指标，轴向质心与跨度为主要证据。（8）蒸汽缓和峰值负荷且壁面热代价有限——20% H₂、零氨下添加10%蒸汽使平均温度降低62.3 K、最大壁面热流降低23.0%、CVHRR降低31.5%，总壁面传热降低3.9%；20%蒸汽时CVHRR与最大壁面热流降幅达52.9%与36.2%，壁面传热代价增至7.7%；交叉细网格结果显示固定U007时C008使峰值壁面热流降低19.3%、CVHRR降低26.2%、平均温度降低50.3 K、HRR质心下移15.1 mm，固定U008时相应为19.0%、26.8%、48.6 K与15.9 mm，而单纯提高速度仅使峰值热流降低4.9%（C007）与4.5%（C008）；固定速度下总壁面负荷变化仅0.003%（U007）与0.249%（U008），取热发生再分布，蒸汽替代案例腔体负荷增加7.2–7.4%、过渡段增加3.4–3.5%、烟囱降低11.8–12.2%，四种案例最大壁面热流均位于x=489.3 mm烟囱处；q99由无蒸汽的2.624×10⁷升至10%蒸汽的3.610×10⁷ W/m³，而固定阈值体积分数f100几乎不变（0.324%对0.329%），与再分布一致；环境响应亦为混合，10%蒸汽大幅降低NOx却在所采用机理内增加预测N₂O。（9）环境选择的热积分筛选——无约束环境选择V1向壁面传递5.057 kW、CVHRR=7.57，与相邻10%蒸汽数据库案例C008（5.049 kW；CVHRR=7.51）在壁面传热、最大壁面热流、CVHRR与平均温度上分别相差0.16%、0.55%、0.82%与0.08%；含氨约束选择V2保持高局部化（CVHRR=10.55）并传递5.113 kW，V3添加5%蒸汽使CVHRR降低18.2%、最大壁面热流降低13.3%、平均温度降低34.1 K、壁面传热降低1.9%；氨本身热影响相对温和但使预测氮排放增加数量级，氨含量主要由环境权衡决定，蒸汽同时调控氮分配与空间热缓和。（10）单元分辨证据与局部化解释的稳健性——表7、表8给出代表性案例的壁面传热、温度、积分放热与局部化统计；最大参比指标依赖未解析极值，不能独立证明反应体积扩展，CVHRR亦依赖稀疏上尾，轴向HRR质心与中心跨度为补充空间证据；V2三网格试验不证明所有场收敛，交叉19 760网格C007-C008计算提供配对细网格趋势检验而非完整渐近收敛研究。（11）轴向反应布置与壁面取热——由C001至C007，20% H₂使HRR加权质心由231 mm前移至117 mm，10–90%轴向跨度由164 mm缩至81 mm，回流体积相近（56.7%与54.3%）但其正HRR份额由16.5%降至1.6%，模型预测上游反应集中同时保留回流区，该体积分数不度量回流量或确立MILD运行；蒸汽部分逆转收缩，C008相对C007使质心下移20 mm、轴向跨度扩大16%，同时最大壁面热流降低23.0%、总壁面传热仅降低3.9%；5% NH₃下5%蒸汽使质心由119移至129 mm，轴向跨度维持约85 mm，CVHRR降低18.2%，壁面负荷集中比由3.60降至3.19；HRR与壁面热流质心间距在C001为46 mm、C007为172 mm、C008为141 mm、V2为171 mm、V3为154 mm，氢富集案例正HRR质心显著位于面积加权壁面热流质心上游，蒸汽缩短该距离并降低峰均比，指示反应与壁面负荷耦合更为分布；轴向壁面剖面显示积分变化并非均匀偏移，氢抬高下游热流峰值并加陡空间变化，蒸汽降低该峰值并将取热展布于更宽轴向区间，x≈490 mm处最大值位于压力出口旁烟囱内，在缺乏出口延伸检验与分区负荷的情况下不能定义腔体设计负荷或有效工艺热回收，42–49%的化学机理差异超过所报道23%的蒸汽峰值降幅。（12）受控化学机理敏感性——PAN计算保持入口状态足够接近以限制总流不匹配，V2燃料、氧化剂与出口质量流差异均低于0.001%，V3最大差异0.17%，能量不平衡为V2-PAN 0.0199%、V3-PAN 0.0299%；积分量相对不敏感，PAN使CO₂变化不超过0.33%、积分HRR不超过0.16%，但N₂O在V2降低72.5%、V3降低87.8%，NO变化?5.5%与+3.5%，CO增加70.7%与56.7%；PAN使V2与V3的CVHRR分别增加10.1%与13.7%、最大HRR增加21.6%与29.8%、最大绝对壁面热流增加42.0%与48.6%，体积平均温度分别升高92.8 K与100.5 K，壁面热流质心下移并使反应-壁面质心距离由171增至211 mm（V2）、由154增至200 mm（V3），PAN还在平均温度更高情况下使总壁面负荷降低5.79–8.02%，表明平均气温单独不能决定壁面热流；两种机理下蒸汽均产生相同方向的热响应（V2至V3平均与最大壁面热流、最大HRR与CVHRR均降低），但氮响应未被定量保持，合并机理预测5%蒸汽使N₂O增加254%而PAN在更低基线下预测增加57.9%并逆转小幅NO下降，CVHRR降幅分别为18.2%与15.5%，最大壁面热流降幅为13.3%与9.2%，局部化方向比绝对或路径特定的氮论断更稳健，该方向与交叉细网格C007-C008比较一致但未确立对所有组成的稳健性。（13）EDC参数与壁温敏感性——EDC时间尺度常数扰动对主要积分响应影响较小，V2与V3间修正排放变化至多3.22%、总壁面负荷0.04%、积分HRR 0.002%，q99至多变化4.00%、CVHRR 2.53%、单单元最大值8.22%，两个扰动Cτ值下蒸汽诱导的CVHRR降低均保留；Cγ扫描中积分HRR偏差至多0.022%、总壁面负荷0.056%、体积平均温度0.097%、峰值壁面热流0.217%、修正NOx 0.583%，而最大温度（至多4.96%）、CVHRR（37.1%）、单单元HRR最大值（98.4%）、修正CO（198.0%）与修正N₂O（119.3%）的基线相对偏差更大，故积分热响应与NOx在测试Cγ范围内局部不敏感而微量组分与上尾反应统计敏感，两种扫描均未覆盖扩展EDC、PaSR或输运PDF的模型形式不确定性；壁温由1173.15 K改为1073.15与1273.15 K使积分HRR仅变化+0.014%与?0.015%，壁面负荷变化+6.81%与?7.61%，修正NOx变化小于2.4%，但N₂O在1073.15 K增加154%、在1273.15 K降低54%，壁温因此是氮分配的 influential 边界变量而非已证明的控制策略，更热壁面还降低总壁面负荷、改变局部负荷且不代表共轭传热；近常数积分HRR不意味着反应结构、壁面负荷或氮分配不变，回流与分布式HRR本身不证明MILD区间保持。

讨论部分总结：无约束等权解V1避开了燃料氮、采用H₂上界并引入中间蒸汽稀释，保留大部分氢相关壁面热增益同时降低峰值壁面热流与CVHRR；含氨约束V2回答不同设计问题，其确认CO₂强度比V1低4.97 g kWh^?1(thermal)（4.0%），但确认NOx高约1063 ppm，表明碳收益不能脱离氮代价呈现，且5 mass% NH₃按示例LHV仅贡献约1.49%燃料热输入，该结果为约束边界内模型预测权衡而非实际脱碳最优。蒸汽并非均匀污染物抑制劑：5% NH₃下V3添加5%蒸汽相对V2降低CO、Tmax、壁面传热与CVHRR，却使修正N₂O由0.326升至1.315 ppm；零氨切片中10%蒸汽相对C007使HRR质心下移20 mm、中心80%轴向HRR跨度扩大13 mm，支持再分布而非简单抑制。入口审计表明质量流与动量须区分：交叉案例空气动量流率为0.4097 N（C007-U007）、0.4579 N（C007-U008）、0.3876 N（C008-U007）、0.4331 N（C008-U008），固定速度下蒸汽替代因密度与质量流变化使动量降低5.41%，速度由U007升至U008使任一组成动量增加11.76%，但匹配速度的蒸汽替代对比仍保留19.0–19.3%峰值壁面热流降幅而单纯速度仅4.5–4.9%，表明峰值负荷趋势主要不是射流速度假象，但因蒸汽替代同时改变氧化剂组成与氧输送，未能分离纯蒸汽热或化学效应。NH₃氧化主要经H、O、OH夺氢生成NH₂，其后氮通量可经NH、HNO与成NO路径或经还原反应返回N₂，N₂O可由含NH中间体与NO耦合反应生成而其破坏强烈依赖温度与H/O/OH自由基池，蒸汽可同时降低温度并改变自由基可用性，故NO降低未必伴随N₂O降低；现有场支持该解释但不构成反应路径计算，50案例数据库未归档物种生成速率与基元反应通量，主导路径不能定量指派，机理用于解释可能的氮分配而非宣称反应特定因果，未来路径分析应对V2与V3在主反应区积分基元速率。Minamoto等（Minamoto et al.）描述相互作用烃类MILD反应区，Sorrentino等（Sorrentino et al.）识别氨扩散点火中燃料依赖结构，本研究的氢诱导上游偏移与低氨替代较小空间效应与反应性变化一致，但混合分数条件HRR、自由基场与基元速率不可得，本CFD场不重现也不验证这些详细结构分析，PAN同样表明近等积分HRR可伴随不同局部化且速率与物性效应仍合并。轴向质心区分反应布置与取热：氢使反应质心由231前移至117 mm而壁面热流质心保持约290 mm，反应-壁面距离由46增至172 mm；蒸汽使反应质心下游移并将该间距降至141 mm同时降低壁面负荷集中；CFD结果因此将正HRR质心显著置于面积加权壁面热流质心上游，总壁面传热单独不能揭示该模拟输运路径；因施加等温壁，本研究不能预测壁温非均匀性或其与质心间距的相关性，多轴向分辨壁面热流或外壁温度测量可提供适当实验检验。氨在所考察切片中空间效应小于氢或蒸汽：零蒸汽下NH₃由0增至10%仅使HRR质心由117移至121 mm、中心跨度由81变为82 mm，10%蒸汽下相应质心由137变为142 mm，该数据集中氨的主要代价是氮物种生成而氢与蒸汽控制大部分反应布置。所报道比较各有范围：甲烷场在一个实验站评估，主要响应与积分热量的网格依赖性适中，质能闭合小，留出与全新启动误差量化代理插值，边界匹配PAN计算暴露化学与物性数据合并敏感性；这些确定性检验非统计重复，不能合并为置信区间，百分比变化是数值模型内效应量而非统计显著性检验。PAN比较缩小一种不确定性而未消除：积分HRR与CO₂稳定而N₂O、CO、温度、CVHRR与峰值壁面热流对机理敏感，蒸汽响应在局部化与壁面负荷上的共同方向因此比绝对氮水平更可辩护；±20% Cτ试验变化远小于机理替换并保留蒸汽诱导的局部化降低，Cγ扫描同样使积分HRR、总壁面负荷、平均温度、峰值壁面热流与NOx近乎不变却使最大温度、CVHRR、单单元HRR最大值、CO与N₂O大幅变化，故所测EDC常数支持积分热趋势的局部稳健性而非局部极值或微量组分，且未涉及替代EDC形式或确立氨特定校准；壁温试验同样显示积分HRR稳定而壁面负荷、反应局部化与N₂O显著变化，这些发现将化学形式与热边界设定识别为局部与微量组分预测的主导限定；因缺乏可比MILD条件下含氨测量，两种机理均不能排序，PAN差异仍为敏感性度量而非精度改进。氢在基线模型内收缩并前移反应区，蒸汽对CVHRR与最大壁面热流的缓和在V2-V3两种化学包下以及交叉细网格C007-C008两种固定空气速度下均复现，后者为配对细网格趋势检验而非完整网格收敛证明，出口独立性仍未解决；经典MILD判据需局部混合预反应温度、混合物依赖点火温度与适当定义的反应温升，入口温度不描述回流后局部混合状态，这些量不可得，故MILD指参考炉膛概念，所有混合气下区间保持仍未验证，绝对氮水平与局部极值仍以化学、湍流/EDC、辐射、等温壁与压力出口为条件。结果识别了检验设计结论最有效的测量：含氨 campaig n应同时报告干基O₂修正NO、NO₂、N₂O、NH₃与CO及轴向或径向温度剖面，壁面量热或分辨热流测量应在多个轴向位置获取而非仅积分负荷，以检验氢驱动上游局部化与壁面取热下游位移是否在共同CFD模型之外持续；对燃烧器评估，环境档案应首先识别候选碳-氮折中，随后以轴向HRR跨度、回流HRR参与、壁面热流质心与峰均壁面负荷筛选候选组成，该序列不规定几何特定运行点但指示在组成进入新燃烧器或壁面热条件前应监测的量，入口核算、载入化学、区间分类与下游边界影响需在用于燃烧器设计前验证。

研究结论翻译：研究人员将四目标CFD-代理框架与壁面传热及空间热负荷分析、三网格试验、守恒检验、边界匹配的化学与物性数据比较以及受控EDC参数与壁温敏感性研究相结合，用于含氢、氨与蒸汽的甲烷MILD燃烧。响应特异性代理以0.053–8.34%的MAPE预测5个未触碰案例，3个全新启动计算在15%以内复现全部12个优化后响应。等权TOPSIS在无约束档案选择0% NH₃、20% H₂与约9.8%蒸汽，在NH₃约束至至少5%时选择5% NH₃、20% H₂与零蒸汽；线性极差重排序保留约束等权选择但将无约束选择改为约0.03% NH₃、20% H₂与零蒸汽，中间蒸汽选择因此依赖归一化，基于影响的排序仍未评估；甲烷参考基准还将约束等权选择改为约5% NH₃、1.88% H₂与零蒸汽，确认对归一化参考的依赖。生产网格在V2处对四个优化响应的预测与最细网格值相差5.55%以内而峰值壁面热流相差10.82%；19 760网格上的交叉C007-C008检验在两种固定空气速度下复现蒸汽缓和方向，但两项计算均不构成完整渐近收敛研究。添加20% H₂使壁面传热增加14.2%，但使峰值壁面热流增加42.5%、CVHRR由1.97升至10.95，并使HRR质心上移49%、中心80%轴向HRR跨度减半、回流HRR参与由16.5%降至1.6%。在所记录生产网格条件下，添加10%蒸汽使峰值壁面热流降低23.0%、CVHRR降低31.5%、壁面热代价3.9%，同时使反应质心下移并展宽轴向HRR分布；在19 760网格上匹配空气速度保留19.0–19.3%峰值热流降幅与26.2–26.8% CVHRR降幅，而单纯速度产生4.5–4.9%峰值降幅，固定速度对比仍含耦合氧化剂组成与氧输送效应，故非孤立纯蒸汽效应。壁面热流质心保持于反应质心下游，其间距由甲烷46 mm增至加氢后172 mm、加10%蒸汽后降至141 mm，该反应-壁面输运度量将反应局部化与下游壁面取热联系起来；出口邻近烟囱峰值与规定等温壁阻止其解释为腔体设计极限或有效工艺负荷。PAN使积分HRR变化小于0.16%但使最大壁面热流变化42.0–48.6%、N₂O变化72.5–87.8%；相对基线Cγ=2.1377，±20% Cγ扫描使积分HRR至多变化0.022%、总壁面负荷0.056%、峰值壁面热流0.217%、NOx 0.583%，但使CVHRR至多变化37.1%、单单元HRR最大值98.4%、修正CO 198.0%、修正N₂O 119.3%。两种化学包下蒸汽均使V2至V3的CVHRR与最大壁面热流降低，交叉细网格案例在匹配空气速度后保持相同方向；绝对氮预测、微量组分与局部极值仍依赖模型。含氨炉膛验证、替代湍流-化学闭合、分离速率/物性效应、完整配对网格收敛、出口延伸敏感性与逐案例MILD评估仍需进一步工作，所报道候选并非经验证的燃料灵活运行点。