基于化学反应动力学的DME发动机燃烧与排放特性研究.pdf

华中科技大学博士学位论文基于化学反应动力学的DME发动机燃烧与排放特性研究姓名：莫春兰申请学位级别：博士专业：动力机械及工程指导教师：张煜盛20070512I摘 要 摘 要 能源短缺，排放法规日趋严格，柴油机的可持续发展面临严峻考验近年来，二甲基醚 （dimethyl ether,DME） 由于其良好的着火性能而用作柴油机替代燃料的系列研究引起了人们的加倍关注基于此,本课题结合国家 973 项目的内容,借鉴国内外研究工作的新进展，通过 DME 发动机缸内直喷（DI）和均质充量压燃（HCCI）两种不同的燃烧方式，应用详细化学反应机理进行了模拟计算和实验对比，为实现 DME 发动机的高效、低污染燃烧提供理论依据和实施方向 本文的研究内容和所取得的成果主要有以下几个方面： 对 DME 发动机实验与数值模拟研究的现状进行了综述，全面分析了柴油机燃用DME 的技术措施及其有关问题以 ZS195 型直喷式柴油机为原型机，构建了 DME 发动机专用实验台架，在其上开展了 DME 发动机性能及常规排放物的测试研究针对 DME特殊的理化性能，进行了 DME 经高压燃油泵柱塞套间隙的泄漏量测量实验以及油管压力波实验，分析了 DME 热物性对燃油喷射特性和高压燃油泵泄漏量的影响；探讨了启喷压力、供油提前角、气门开启时刻、喷孔直径和压缩比等对 DME 缸内直喷发动机性能与排放的影响。

台架实验结果表明，与燃用柴油相比，燃用 DME 时缸内直喷发动机在额定转速下，碳烟排放极低，NOX排放随着负荷的增加而增大，但燃用 DME 的 NOX排放比柴油低；在中低转速下，HC 排放比柴油高，且 HC 排放随着有效功率的升高呈下降趋势与 DME 发动机缸内直喷方式相比，HCCI 燃烧方式可有效抑制氮氧化物的排放，甚或降低至零排放水平，但 HC 排放明显增加因此，有必要对 DME 发动机的 HC排放加以重点研究及控制 在国内外首次开展了DME发动机非常规排放物的定量测试研究为此，开展了DME发动机非常规排放物检测方法的研究，包括气相色谱法（GC）和傅立叶变换红外分析法(FTIR)在GC实验中，首先对配有氢火焰离子化检测器的气相色谱仪的分离柱温、载气流速、气化室温度和检测器温度等工作条件进行了优化，采用外标法，经色谱工作站绘制出甲醛和甲酸甲酯的工作标准曲线，继而对发动机尾气进行检测，并将其结果与FTIR检测结果进行对比，确认FTIR可作为检测DME发动机非常规排放物的有效手段进一步的FTIR实验结果表明，DME直喷发动机的排放中存在甲醛、甲酸甲酯和甲酸三种非常规排放物，其中甲醛排放浓度最高，甲酸甲酯排放浓度次之，甲酸排放浓II度最小，甲酸排放浓度最低时低于10ppm。

DME的HCCI燃烧方式与缸内直喷相比，前者甲醛和甲酸甲酯的排放浓度较高 应用详细化学反应动力学计算软件 CHEMKIN，就 DME 发动机着火和燃烧过程开展了性能预测与计算结果的实验验证工作,构建了 DME 发动机的着火数据库，并深入探讨了热力学燃烧模型与着火数据库耦联的方法及其应用于DME发动机循环模拟计算的适用性模拟计算结果表明,滞燃期是缸内温度、压力和燃空当量比的函数，在一定的燃空当量比范围内，着火滞燃期随燃空当量比增大而变小计算结果与验证实验结果对比表明,Wiebe 模型、Watson 模型和 Whitehouse-Way 模型均可应用于 DME 发动机工作过程的模拟计算，其中，以着火数据库与 Watson 燃烧子模型相耦合的模拟计算结果与实验结果最为吻合该项研究对于发展简洁而又实用的 DME 发动机工作过程模拟计算软件具有现实意义 开展了发动机数值计算专用软件KIVA3V的微机化移植、液态DME燃料数据库建库、引入液滴破碎模型以及燃烧模型改进等项基础工作，藉此对DME缸内直喷发动机的性能与非常规污染物排放量进行了三维数值模拟和预测研究其中，液滴破碎过程的模拟分别采用TAB模型、基于表面波不稳定理论的KH-RT模型；燃烧模型分别采用了EDC模型、基于详细化学反应机理的部分搅拌器概念模型（PaSR）。

模拟计算考虑了气相可压缩性、流场不均匀性及液滴破碎对发动机动力及排放性能的影响三维湍流燃烧模拟关于直喷发动机缸内压力随曲轴转角的变化，以及NO排放浓度的预测值与实验结果吻合较好同时，计算分析了缸内流场速度、温度和组分浓度分布随曲轴转角的变化历程；就进气加热、缸内温度、气体流动、壁面传热、缝隙效应，以及喷雾模型和化学反应机理的选择对甲醛、甲酸甲酯生成及其氧化过程的影响进行了参数分析；进一步探讨了预测值小于实验值的原因以及甲酸甲酯的生成机理 通过构建可应用于HCCI发动机燃烧排放研究的DME详细化学反应动力学燃烧模型，对DME均质充量压燃发动机的非常规排放物生成机理与控制策略进行了深入研究详细模型包含97种物种和457个基元反应通过模拟及其计算结果的实验验证，证实了详细模型的实用性和有效性，判明了DME氧化反应的主要历程与甲醛（CH2O）、甲酸（HCO2H） 和NOX主要生成机理 计算结果表明， NOX排放中NO生成量达到最大值后有 “冻结”现象出现，NO与N2O最终排放浓度极少且受缸内温度影响不大；随着缸内温度的增加,NOX排放中NO所占比重逐渐增加基于化学反应速率及敏感度分析得知，DME发动机HCCI燃烧的NO生成主要有扩充的切尔多维奇机理和N2O生成机理这两个途径。

III为了节省计算机时，并为化学反应动力学模型与CFD多维模型耦合的燃烧计算提供一个行之有效的途径， 在详细化学反应动力学模型的基础上， 构建了包括36种组分、73个基元反应组成的简化机理模型简化模型包括低温反应和负温度系数区子模型、高温裂解和高温氧化子模型以及NO排放子模型利用该模型，可以计算着火点、燃烧速率、甲醛、甲酸和NO排放将简化模型与FLUENT相耦合，进行了DME发动机HCCI燃烧数值二维模拟计算，结果表明，DME发动机HCCI燃烧的甲醛和甲酸排放物浓度较高，其排放浓度受燃空比、活塞环与气缸壁处缝隙和压缩比的影响较为明显 关键词： 二甲基醚（DME），化学反应机理，湍流燃烧，非常规排放物， 傅立叶变换红外分析法（FTIR），气相色谱分析法（GC） IVAbstract Environmental and human health concerns over emissions from internal combustion engines continue to bring about increasingly stringent emission standards and drive research into the development of cleaner-burning fuels. Aiming at this current research situation together with the contents of National Basic Research Priorities Program (2001CB209207) administrated by the State Ministry of Science kVc,为常压下第k种组分的定容比热。

混合物的平均定容比热Vc为: kVKkkVcYc, 1∑ == (4-9)将(4-7)式代入(4-5)式得: 011=++∑∑ ==pdvdYedeYkKkkkKkk(4-10)将(4-8)式代入(4-10)式可得: 01, 1=++∑∑ ==pdvdYedTcYkKkkkVKkk(4-11)将(4-9)式代入(4-11)式可得: 01=++∑=pdvdYedTckKkkV(4-12)(4-12)式对时间t求导得: 01=++∑=dtdvpdtdYedtdTckKkkV(4-13)将(4-4)式代入(4-13)式得: 01=++•=∑dtdvpWevdtdTckkKkkVω (4-14)由于系统容积随时间变化,即)(tVV =,且比容mtVtv)()(= ,则: dtdV mdtdv1= (4-15)将(4-15)式代入(4-14)式得: 5301=++•=∑dtdV mpWevdtdTckkKkkVω (4-16)根据理想气体状态方程: WRTpρ= (4-17)式中，R为普适恒量；W为混合气的平均分子量，∑ ==Kkkk WYW11；ρ为混合气密度,∑ ==KkkkWX1][ρ ，其中][kX是第k种组分的摩尔浓度。

将(4-17)式和(4-9)式代入(4-16)式得到能量方程: 01, 1=++•==∑∑dtdV WmRTWevcYdtdTkkKkkkVKkkρω (4-18)方程（4-18）和（4-4）组成了为由kY和T组成的1+k个未知数的一阶非线性方程组 4.1.1.2 发动机几何模型 方程（4-18）中的dtdV可由发动机气缸容积随曲轴转角的变化函数求得气缸容积随发动机参数的变化函数为： ]sincos1[211)(22θθ−+−+−+=RRC VtVc(4-19)式中，cV为气缸的余隙容积，C为压缩比; 连杆长度/曲轴半径比Ac LLR= ，cL为连杆长度，AL为曲柄半径 对（4-19）式求导，得到气缸容积的变化率： 54] ]sincos1 [sin)21())((22θθθθ−+−= RC dtd dtVtVdc(4-20)联立方程组（4-18）和（4-20）式就能解决发动机中的燃烧问题，并得出整个燃烧过程加的温度、压力和各组分浓度的变化过程 4.1.1.3 化学反应子模型 方程（4-19）中的基元反应中第k种组分的摩尔生成速率•kω由所有和它相关的化学反应的综合作用决定 kKkr kikkkKkf kixxfiri∑∑ ==⇔11νν ),, 1(sNi⋅⋅⋅ ⋅ ⋅ ⋅= (4-21)式中，f kiν和r kiν分别为第i个反应第k种物质的正向和逆向反应计量系数；fik和rik分别为第i个反应的正向和逆向反应速率常数ik。

1889年Arrhenius通过大量实验与理论的认证，提示了在恒定浓度的条件下，反应速率常数对温度的信赖关系，建立了著名的Arrhenius定理常用的是它的指数表达式： )exp(RTETAki iii−=β(4-22)式中,iA和iE是由反应本性决定而与温度无关的常数iA为第i个反应的指前因子，又叫频率因子， 代表分子相互碰撞次数iE为阿累尼乌斯活化能(单位为J/(mol·K)),分子要进行反应必须要有iE以上的能量才行R为通用气体常数其值为8.3144kJ/(mol·K)T为热力学温度RTEi−exp叫波兹曼因子，它相当于能量 在iE以 上的分子数占总分子数的比例 控制方程组中出现的基元反应中第k种组分的摩尔生成速率如下： ⎪⎩⎪⎨⎧⎭⎬⎫−−=∏∏∑ ===•ssr kif kiNiNiiriifiKkf kir kikXkXk111][][)(ννννω(4-23)在某些条件下，一些基元反应受到高压和低压条件的限制，即高压限和低压限时55有着截然不同的反应速率，从而形成所谓的依赖压力的“衰退（fall-off） ”反应，其反应速率常数有如下形式： FPPkkrr i+=∞1(4-24)式中，∞k和0k分别为高压限的速率常数和低压限的速率常数，表达式符合（4-22）式的形式；][X为混合物浓度。

F为速率常数修正函数 ∞=kXkPr][0(4-25)根据F不同的取值， （4-23）式可分为Lindemann形式、Troe形式和SRI形式，可对衰退反应的处理可更真实地反映一些化合反应和含有增强第三体（enhanced third-body）反应的进程 4.1.1.4 热力学参数 比热RCp,焓RTH0,熵RS0的多项式计算如下因为下面三个式子适用于两种温度范围,除了个别比较特殊外，一般断点为1000度，因此输入文件中根据7个系数乘以2种温度范围需要14个系数进行输入 4 53 42 。