混合预混铁粉尘-甲烷火焰：数值分析与表格化火焰面流形方法的验证

《Advanced Powder Technology》：Hybrid premixed iron dust–methane flames: numerical analysis and validation of a tabulated flamelet manifold method

　　

随着清洁能源转型带来的挑战，粉状金属燃料因其作为循环能源载体的潜力近年来受到关注。为确保稳定和安全的燃烧，目前运行的实验性铁粉燃烧器通常依赖气体辅助。因此，在设计此类系统时，理解不同颗粒尺寸、铁与甲烷比例以及组合当量比下气相与铁燃烧的耦合至关重要。此外，降阶模型（reduced-order models）的开发对于实现高效数值模拟和促进快速设计过程至关重要。在本研究中，研究人员使用拉格朗日粒子模型（Lagrangian particle model）在通道中研究了层流预混铁粉-甲烷火焰。一个目标是通过详细化学模拟（detailed-chemistry simulations）阐明铁与甲烷燃烧之间的耦合。因此，研究人员改变了铁和甲烷的当量比以及颗粒直径。研究揭示了以甲烷或铁开始消耗氧气为特征的耦合机制。这些机制影响火焰速度、颗粒温度和转化行为。第二个目标是全面验证基于一维火焰开发的表格化火焰面流形方法（tabulated flamelet manifold method），该方法在先前的工作中开发。在所有研究参数下，详细化学与表格化化学的比较显示，在耦合行为、铁转化、氧气消耗、温度和火焰速度方面具有良好的一致性。

论文解读文章

**研究背景与问题**

清洁能源转型面临太阳能和风能间歇性及分布不均的挑战，需要可储存和可运输的能源载体。粉状金属燃料，特别是铁，因其丰富、低成本、无毒和高体积能量密度，作为循环能源载体近年受到关注。然而，铁燃烧在实际规模上的理解仍处于初期阶段。目前运行的实验性铁粉燃烧器通常依赖甲烷等气体辅助来稳定火焰，这引入两个关键问题：首先，必须理解气体燃料与固体铁燃料之间的耦合机制；其次，需要开发降阶模型以实现高效数值模拟。因此，本研究旨在通过系统数值分析，探究层流预混铁粉-甲烷混合火焰中颗粒尺寸、铁-甲烷当量比和组合当量比对火焰结构与速度的影响，并验证一种基于一维火焰的表格化火焰面流形方法（tabulated flamelet manifold method）的准确性。论文发表在《Advanced Powder Technology》。

**主要技术方法**

研究人员采用欧拉-拉格朗日框架模拟颗粒负载流动。颗粒模型包括一阶表面动力学（FOSK）模型（描述铁氧化为FeO的动力学与扩散限制）和颗粒焓方程（包括对流和辐射传热）。气相模型求解连续性方程和纳维-斯托克斯方程，详细化学采用简化GRI 3.0机制（28种物种、112个反应）。表格化火焰面流形方法构建四维表格（燃料混合分数Z_f、氧消耗参数Y_O、绝对焓h、进程变量Y_c），通过一维自由传播火焰预计算。所有模拟在OpenFOAM^?中进行，使用二阶精度格式。

**研究结果**

**4. Comparison with Experiments and Analysis of Fundamental Flame Features**
通过与实验层流火焰速度数据对比，验证了数值模型。模拟结果再现了实验趋势：火焰速度随甲烷当量比增加而增加，铁颗粒添加通常降低火焰速度。对于高甲烷当量比（?_CH4=0.9和1.0），模拟与实验吻合良好；低当量比下存在偏差，可能源于FOSK模型的高氧条件下动力学率过早触发铁氧化，以及忽略颗粒-自由基相互作用。

**4.1 Comparison with experimentally measured laminar flame speed**
通过准一维模拟与实验热流燃烧器数据对比，发现纯甲烷火焰速度一致性良好，铁颗粒火焰速度趋势一致，但低甲烷当量比时模拟值偏高。

**4.2 Features of typical premixed hybrid flames**
分析了典型混合火焰的温度分布：小颗粒（2μm）温度梯度陡峭，颗粒与气相温度几乎重合；中等颗粒（5μm）气相先升温，颗粒后反应并产生温度过冲；大颗粒（10μm）延迟更明显，过冲温度更高，火焰结构更厚。颗粒尺寸显著影响加热动力学和火焰结构。

**5. Detailed analysis of coupling regimes and its prediction by a tabulated flamelet manifold method**
基于氧气消耗的空间范围，识别出四种耦合机制。Regime 1a和1b中甲烷先消耗氧气，铁随后反应；Regime 2a和2b中铁先消耗氧气。颗粒尺寸和当量比决定耦合机制，小颗粒始终处于Regime 2，中等和大颗粒在低当量比时多为Regime 1。表格化化学方法准确预测了这些机制。

**5.1 Hybrid flame regimes**
通过氧气消耗和焓变分析，分类了四种火焰机制。耦合机制由热惯性和点火温度平衡决定，表格化化学与详细化学结果高度一致。

**5.2 Ignition distance based on oxygen consumption**
大颗粒和高当量比增加点火距离。小颗粒点火距离短，且铁先反应。表格化化学方法在点火距离和机制预测上与详细化学吻合良好。

**5.3 Ignition conditions**
铁在相同氧质量分数下点火温度低于甲烷，但大颗粒热惯性大导致延迟点火。表格化化学方法在颗粒点火条件上一致，气相点火条件略有差异。

**5.4 Global features**
最大颗粒温度在大颗粒、贫燃条件下最高；小颗粒在富燃条件下因机制2b产生高温。火焰速度通常随颗粒尺寸增大而降低，小颗粒在高当量比时出现速度上升（机制2b）。表格化化学方法在温度与速度上误差小于50K和0.025 m s^-1。

**5.5 Comparison of iron conversion and major species**
表格化化学方法在铁转化和主要物种质量分数上平均误差小，低当量比时最佳，高当量比时因氧气竞争出现较大偏差（铁质量分数误差0.1，H₂O和CO₂误差0.01）。

**讨论与结论**

讨论部分指出，FOSK模型在高温富氧条件下可能过早触发铁氧化，且其自燃温度低于最新实验数据，但这不影响表格化方法本身的验证。未来需改进铁氧化模型，并考虑颗粒与自由基相互作用。结论部分翻译如下：
在本研究中，通过详细化学模拟和表格化火焰面流形方法研究了预混铁-甲烷混合火焰。进行了系统参数研究，变化甲烷和铁当量比以及单分散颗粒直径（2μm至10μm），计算了48个案例。结论如下：RQ1：全局当量比主要影响平衡温度，最大火焰速度出现在全局当量比为1时；小颗粒增加火焰速度，中等和大颗粒降低火焰速度。RQ2：基于氧气消耗空间范围识别出两种主要耦合机制（Regime 1和2），热耦合和氧耦合是主要因素，取决于点火温度和热惯性。RQ3：表格化火焰面流形方法在大多数条件下准确表示火焰结构，点火距离、温度、速度、物种质量分数和铁转化均与详细化学结果良好一致，仅富燃条件下有轻微误差。未来工作应改进铁氧化模型，并转向应用导向研究，如湍流中多分散金属粉末行为。