基于CHEMKIN的H2O2添加剂对正丁醇HCCI发动机着火过程影响的仿真分析

周奥，金慧，张春化，李阳阳，李松峰

（1.兰州交通大学机电工程学院，甘肃 兰州 730070；2.兰州工业学院汽车工程学院，甘肃 兰州 730050；3.长安大学汽车学院，陕西省交通新能源开发、应用与汽车节能重点实验室，陕西西安 710064）

随着日益严苛的排放法规的执行和世界范围内能源危机的加剧，内燃机的高效、清洁燃烧技术受到了学者们的广泛关注。基于此，学者们提出了1种新型燃烧方式—HCCI。HCCI燃烧模式结合了传统柴油机和汽油机的优点，使得内燃机的热效率提高的同时降低了颗粒物和氮氧化物（NOx）排放［1-3］。此外，HCCI发动机的燃料适应性强，近些年备受关注的醇类替代燃料被广泛应用于HCCI燃烧模式的研究，尤其是甲醇、乙醇和正丁醇［4-6］。

但是，HCCI模式的燃烧过程受限于化学反应动力学和燃烧边界条件的控制，导致该种燃烧方式的燃烧相位不可直接控制，同时还存在发动机低负荷易失火，高负荷易爆震以及运行范围窄等技术难题［7-9］。然而，鉴于HCCI模式的显著优点，学者们一直致力于突破该种模式技术壁垒的探索研究，试图通过进气加热、进气增压、可变配气相位、废气再循环、燃用理化特性不同的混合燃料和添加活性添加剂等措施来实现这一目标［10-13］。王帅［14］研究了进气加热对天然气发动机性能的影响，进气温度升高使得着火延迟期缩短，但是热效率下降较为明显。Hasan等［15］研究了乙醚/乙醇混合燃料的燃烧特性，结果表明增加乙醇在燃料中的比例可以延长着火延迟期；但是，对于所有试验燃料，随着混合气当量比的增加，发动机热效率和缸内压力均降低。此外，唐志刚等［16］基于微型内燃机从机理方面研究了过氧化氢（H2O2）添加剂对甲醇燃料燃烧过程的影响，研究表明H2O2添加剂主要是分解生成羟基（OH）来促进甲醇燃料的燃烧。邢一帆等［17］借助于HCCI模型研究了H2O2添加剂对正庚烷HCCI发动机缸内燃烧压力和温度的影响，研究发现H2O2添加剂改善了缸内正庚烷的燃烧，提高了燃油经济性，但是发动机的动力有所下降。本课题组先期研究也发现H2O2添加剂对正丁醇HCCI发动机的燃烧与排放特性均有所改善，并且可以有效拓展正丁醇HCCI发动机的小负荷运行极限［18］。但是，几乎没有研究基于正丁醇HCCI发动机从机理方面详细地阐述H2O2添加剂对正丁醇低温氧化过程的影响。

为了更为深入地了解H2O2添加剂对正丁醇HCCI发动机着火过程的影响，本文基于正丁醇燃烧机理针对H2O2添加剂对于正丁醇低温氧化的促进机制进行研究，主要研究H2O2浓度变化对正丁醇低温氧化过程中关键组分和相关基元反应的影响，揭示其优化正丁醇燃烧的主要机制，为后续相关研究提供理论依据。

1 模型建立与验证

本文选用CHEMKIN PRO软件中的封闭均质反应器模型和封闭内燃机模型进行了H2O2添加剂对正丁醇低温氧化过程影响的研究，文中所用H2O2水溶液中H2O2体积分数为30%。封闭内燃机模型技术参数如表1所示，并且文中使用的试验对比数据为已发表的论文数据［19］。此外，本文使用的正丁醇机理是由Sarathy等［20］提出的详细机理，其包含426种组分，2 335个基元反应。在模拟过程中，正丁醇被选作初始反应物，H2O2作为添加物质输入。同时，模型假定发动机缸内工质浓度、温度和压力等为均匀状态；在模拟计算过程中不考虑工质泄露等因素，发动机缸内处于完全封闭状态；计算过程考虑热损失，瞬时换热系数基于Woschni公式相关参数设定计算而来，其中系数C11、C12和C2分别为2.28、0.308和0.003 24。内燃机模型的模拟工况为进气温度130°C、进 气 压 力0.1 MPa、发 动 机 转 速 为1 200 r/min、当量比为0.333。

表1 原型柴油机主要参数Table 1 Main parameters of the test diesel engine

此外，为验证模型数据的准确性，本文对所用封闭内燃机模型模拟结果与试验数据进行了对比。图1显示了不同H2O2浓度下缸内压力的计算值和试验值。试验数据表明H2O2溶液的加入促使正丁醇HCCI燃烧相位提前，模拟数据也呈现同样的变化趋势。然而，由试验数据可见，当H2O2水溶液体积浓度达到3%时，缸内峰值压力大约升高了0.5 MPa，而模拟数据的提升效果不明显。这是由于本模型为单区模型，同时忽略了缸内混合气浓度和温度的分层对高温燃烧过程的影响，导致峰值压力变化没有试验数据变化显著。但是，在压缩段截止燃烧开始时刻时间段内，在各个工况下，模拟数据与试验数据的缸内压力曲线基本贴合，能够准确地预测HCCI燃烧模式的着火正时，并且本研究侧重于燃料着火前的低温氧化过程的研究。

图1 不同H2O2浓度下缸内压力的计算与试验数据对比Figure 1 Comparison of calculated and experimental data of cylinder pressure under different H2O2 concentration

图2给出了不同H2O2浓度下，着火时刻的试验与模拟数据的对比［19］。由图可见，随着H2O2浓度的升高，试验所得和模拟计算的着火时刻都不断提前。由试验数据可见，当H2O2水溶液体积浓度达到7%时，相比于纯正丁醇燃料，HCCI发动机的燃烧相位大约提前来了6°CA。当缸内加入H2O2水溶液后，H2O2分解使得高活性OH的浓度得到明显提升，促进了正丁醇的脱氢反应速率，加速了正丁醇的氧化，因此燃烧提前。同时，模拟数据完美复现正丁醇在HCCI模式下的着火时刻变化趋势，并且试验数据和模拟数据基本相同。

图2 H2O2溶液体积分数对着火时刻的影响Figure 2 Effect of volume fraction of H2O2 solution on the ignition time

综上所诉，结合缸内压力数据和着火时刻数据的对比，本模型可以准确计算混合气自燃前的燃烧过程数据，并且本文侧重于不同H2O2浓度对正丁醇低温氧化过程影响的分析。因此，本模型可以用于本文的相关模拟研究。

2 结果与分析

2.1 正丁醇的敏感性

本文分析了不同H2O2溶液体积分数对正丁醇HCCI模式下低温氧化过程的影响。如图3所示，在初始温度800 K和初始压力1.3 MPa条件下，无论H2O2添加剂浓度如何变化，正丁醇低温氧化过程中对温度影响最大的基元反应都是nC4H9OH+OH=C4H8OH-3+H2O。当H2O2溶液加入后，基元反应H2O2（+m）=OH+OH（+m）对低温氧化过程促进作用的影响显著增强，而与C4H8OH-1相关的基元反应对燃烧过程抑制作用的影响被逐步削弱。同时，由图可见，对于正丁醇低温氧化过程影响最为显著的基元反应均有OH和HO2参与，并且OH占据主导地位。

图3 温度800 K时不同H2O2浓度下正丁醇的敏感性分析Figure 3 Sensitivity analysis diagram of n-butanol under different H2O2 concentrations at 800 K

2.2 4种同分异构体的生成速率

图4为不同H2O2浓度下C4H8OH的4种同分异构体的生成速率变化规律。可以看出，C4H8OH-1的生成速率最高，C4H8OH-3和C4H8OH-4次之，C4H8OH-2的生成速率最低约是C4H8OH-1的四分之一。同时，随着H2O2添加剂浓度的增加，C4H8OH-1的生成速率逐渐减小，而C4H8OH-2、C4H8OH-3和C4H8OH-4的生成速率呈现上升趋势，当H2O2溶液体积分数为3%时，变化趋势最为明显，而H2O2溶液体积分数为7%时，变化趋势趋于平缓。因此，结合敏感性分析结果和生成速率分析结果，可以确定添加H2O2溶液后，对正丁醇低温氧化过程中与C4H8OH-1、C4H8OH-3和C4H8OH-4相关的基元反应影响最为明显。

图4 温度800 K时不同H2O2浓度下C4H8OH的四种同分异构体的生成速率分析Figure 4 The formation rate analysis diagram of four isomers of C4H8OH under different H2O2 concentration at 800 K

2.3 OH、HO2和O2的摩尔分数变化

图5显示了不同H2O2浓度下OH、HO2和O2的摩尔分数变化规律。如图所示，随着H2O2浓度的升高，OH和HO2的浓度明显升高。在相同的曲轴转角下，当H2O2添加剂体积分数增加到3%时，OH和HO2的浓度呈现大约2个量级的增加，而H2O2添加剂体积分数继续增加时，OH和HO2的浓度提升幅度减弱，O2的浓度整体上呈现下降趋势。这是由于H2O2溶液的加入，H2O2的分解生成OH，OH再参与正丁醇脱氢反应消耗，而在后续反应过程中又继续生成，最终促使低温氧化过程中OH浓度的提升。对于HO2而言，正丁醇由脱氢反应生成C4H8OH-1，其再与O2发生反应绝大部分生成HO2和丁醛，使得HO2浓度增加。此外，在正丁醇着火之前，O2浓度的下降则说明在正丁醇低温氧化过程中O2的消耗程度加剧，抑制了正丁醇后续反应过程中相关的裂解反应的发生，使得有O2参与的氧化反应比例增加。

图5 不同H2O2溶液体积分数下OH、HO2和O2的摩尔分数Figure 5 Mole fractions of OH，HO2 and O2 under different volume fractions of H2O2 solution

2.4 对正丁醇低温反应过程的影响

为了进一步揭示不同H2O2溶液体积分数对正丁醇低温反应过程的影响，使用反应路径分析法对不同H2O2浓度影响下正丁醇的反应路径进行了分析。图中各个反应路径的占比是从压缩起点到燃料消耗20%时间段内对应相位反应速率曲线的积分。如图6所示，随着H2O2溶液体积分数的上升，正丁醇由OH和HO2主导的脱氢反应总消耗比例逐渐增加，其中，OH参与的脱氢反应整体比例上升，而由HO2主导的脱氢反应整体比例下降。在同一工况条件下，正丁醇脱氢产物的整体生成比例增加，但是C4H8OH-1的比例略有下降。同时，C4H8OH-1的裂解反应比例下降，而C4H8OH-1的加氧反应比例明显增加，继续反应生成物中ROO和HO2比例增加。ROO继续反应，一部分异构化生成QOOH，一部分通过分解反应生成丁烯醇和HO2，QOOH的生成比例先增大后减小，而丁烯醇和HO2的生成比例呈现单调递减的趋势。相比于C4H8OH-1，C4H8OH-2、C4H8OH-3和C4H8OH-4这3种自由基加氧反应生成ROO的比例增加显著，ROO发生分解反应生成的丁烯醇和HO2比例均有所下降，其中，由C4H8OH-2继续反应生成QOOH的比例有所下降，C4H8OH-3和C4H8OH-4继续反应生成QOOH的比例增加。但是，C4H8OH-2继续反应产物中丙醛和OH的比例有所上升。所以，正丁醇在低温氧化过程中最终生成的HO2和OH量有所增加，同时醛类产物增多，而烃类产物减少。此外，当H2O2溶液体积分数进一步提高时，C4H8OH-1继续反应产物中QOOH比例有明显降低，这也使得正丁醇着火延迟期缩短的趋势减缓，和实验现象相吻合。

图6 H2O2溶液体积分数对正丁醇HCCI发动机的低温反应路径的影响［19］Figure 6 Effect of volume fraction of H2O2 solution on the reaction path of n-butanol HCCIengine at low temperature

综上所述，正丁醇低温氧化主要始于由OH和HO2参与的脱氢反应。当H2O2溶液加入后，添加剂中的H2O2分解生成OH，在高活性自由基OH的作用下，由OH主导的脱氢反应消耗了更多的正丁醇燃料，再经历加氧反应和异构化反应生成QOOH以及后续一系列反应来完成正丁醇的消耗，并且主要集中于有C4H8OH-1、C4H8OH-3和C4H8OH-4参与的3条反应路径上。

3 结论

1）随着H2O2溶液体积分数的增加，与C4H8OH-1相关的基元反应对燃烧过程抑制作用的影响被逐步削弱，基元反应nC4H9OH+OH=C4H8OH-3+H2O对温度的影响程度始终最大。

2）H2O2添加剂极大地提高了正丁醇后续反应中有O2参与的氧化反应比例，减少了裂解反应的比例，使得生成物中的醇类和烃类产物减少。

3）H2O2添加剂促使缸内OH和HO2的浓度升高，同时使得正丁醇低温氧化过程中由OH主导的脱氢反应比例增加，而由HO2主导的脱氢反应整体比例下降。