喷水对天然气发动机性能影响的数值模拟

2020-01-01吴静涛吴志军李理光

同济大学学报（自然科学版） 2019年12期

吴静涛, 吴志军, 李理光, 康哲, 邓俊

(1. 同济大学汽车学院, 上海 201804; 2. 重庆大学汽车工程学院, 重庆 400044)

随着汽车工业的蓬勃发展,能源消耗和环境保护问题变得日益严峻.为了应对这些问题,各国政府建立了严格的法规,从而控制汽车的油耗与排放,达到节能减排的目的.严格的法规促进了对发动机节能减排的探索.喷水技术在提高燃油经济性和降低排放等方面有着极大的优势,在汽油机和柴油机上具有较高的应用前景,各国科研人员也对喷水技术开展了大量研究.

目前,喷水技术按照喷射位置的不同可以分为进气总管单点喷水、进气歧管喷水和缸内喷水等方式.进气歧管喷水技术可以使汽油发动机在满负荷下以当量比运行,还可以使发动机的平均有效压力提高5%,热效率提高4%[1].吴志军等[2-3]讨论了喷水发动机的理论循环热效率,随着喷水质量及喷水温度的增加,发动机的热效率增加.于潇等[4-6]在内燃兰金循环发动机上研究了喷水时刻和喷水温度对燃烧特性的影响,实验结果表明,在燃烧过程中喷入的水可控制燃烧速率,喷水时刻对HC和NOx排放特性有较大影响,喷水温度的升高可以提高循环热效率.De Bellis等[7]使用仿真的方法研究了进气道喷水对燃油经济性的影响,结果表明,喷水技术可以有效降低发动机中高负荷的燃油消耗率. Tauzia等[8]在一台柴油机上研究了进气歧管喷水的冷却作用,实验结果发现无论是在大负荷还是小负荷,喷水均可以降低50%的NOx排放.Bozza等[9]对比研究了喷水和废气再循环(EGR)对热力学和燃烧热性的影响,研究表明与EGR类似,喷水对层流火焰速度影响显著,可以有效抑制发动机的爆震.此外,Hoppe等[10-11]在一台配有废气再循环的汽油机上研究了当量比下缸内喷水对发动机油耗和性能的影响,实验结果表明在低油耗区可以将发动机效率提高3.3%~3.8%,在大负荷区域可以将发动机效率提高16%.Iwashiro等[12]研究了喷水对均质压缩燃烧发动机(HCCI)的影响,研究发现HCCI发动机的热效率提高了2%,工作范围从460 KPa拓宽至700 KPa.

天然气因其具有资源丰富、成本低、排放低等特点成为具有前景的替代燃料[13-14].目前,天然气发动机广泛采用稀薄燃烧技术、涡轮增压、高压缩比等技术[15-16].由于排放限值的进一步严格,天然气发动机向当量比燃烧发展,这进一步了天然气发动机的热负荷.面对这些挑战,喷水技术被认为是一种潜在的解决方案.综上所述,目前关于喷水技术的实验与仿真研究主要应用于汽油机和柴油机.因此,本文结合详细的化学反应动力学,运用数值模拟的方法,分别从喷水质量和喷水位置两方面入手,研究分析了喷水对重型天然气发动机燃烧特性和排放的影响.

1 计算模型

1.1 仿真模型

本文研究对象是一台四冲程点燃式涡轮增压天然气发动机,其主要参数见表1.仿真采用三维计算流体动力学(CFD)软件CONVERGE对天然气发动机的工作过程进行模拟.利用有限体积法对质量、动量和能量守恒方程离散求解,湍流模型为RNGk-epsilon双方程模型.使用Kelvin-Helmholtz和Rayleigh-Taylor混合破碎模型以及O' Rourke液膜模型模拟喷雾过程.采用详细化学反应动力学模型(SAGE)计算燃烧反应过程,其中甲烷的反应机理为GRI-mech 3.0,共包含53种组分和325种基元反应.该机理适用于当量比0.1~5下的化学反应.此外,使用扩展的Zeldovich模型预测氮氧化物的生成.此外,仿真所用的天然气燃料的主要成分为:CH4、C2H6和C3H8,其体积分数分别为88%、6%和6%.

表1 发动机参数

1.2 计算网格

为了兼顾计算效率和计算结果的精确性,体网格基本尺寸为4 mm.仿真过程中采用了速度及温度自适应加密体网格(AMR),加密等级是3,即最小网格尺寸为0.5 mm.同时,在燃烧室、火花塞、活塞顶面、气门封闭面等关键位置采用了不同等级的嵌入式固定加密.计算过程中,每个时间步长都会重建网格,因此在不同时刻,缸内网格数目与分布也不相同,最大网格数为300万.图1为上止点前320°曲轴转角时生成的计算网格.

1.3 边界条件

为了研究天然气发动机的进气、压缩及做功过程,本文的数值模拟计算区域选定为进、排气道和燃烧室三部分.进、排气道部分分为气门座、气门体和气道体区域,燃烧室部分分为气缸盖、缸壁和活塞顶.借助GT-Power软件,完成天然气发动机的热力循环计算,获取三维数值模拟的瞬态边界条件.如进口、出口的温度和压力等,其他边界条件的设置见表2.

图1 计算网格

1.4 模型标定

为了验证本文模型的准确性,保证计算结果可靠性,选取当量比为1.21,1 100 r·min-1全负荷工况下的台架试验结果与仿真结果进行对比.图2为实验测得的缸压曲线与仿真结果的对比.从图中可以看出,两条曲线拟合程度良好,误差范围在5%以内,说明模型的选择、边界条件和初始条件的设置能够满足计算的要求.

表2 边界条件

图2 计算模型标定

为了研究不同喷水位置对天然气发动机燃烧和排放特性的影响,选取当量比为1.21,喷水和燃气质量比为0.5,1 100 r·min-1全负荷工况下对进气总管及进气歧管两个喷水位置的仿真结果进行对比研究.对于进气总管喷水,在低速大负荷条件下,喷入进气系统的水以水蒸气的形式从进气边界进入.对于进气歧管喷水,喷嘴位置如图 1所示.为了排除水喷雾蒸发过程对发动机性能的影响,选取在总管喷水条件下研究喷水质量对天然气发动机燃烧过程的影响.由于当喷水与燃气质量比大于1:1时易发生失火现象,因而在当量比为1.21, 1 100 r·min-1全负荷工况下对质量比为0:1、0.5:1和1:1三种工况下的仿真结果进行对比研究,分析喷水质量对天然气发动机缸内流动和燃烧过程的影响.

2 计算结果与讨论

2.1 喷水位置对天然气发动机的影响

图3为不同喷水位置下燃烧室内的平均压力及平均放热率曲线的变化曲线.从图中可以看出,当采用进气总管喷水时,燃烧室内的峰值压力为8.964 MPa,而采用进气歧管喷水时,燃烧室内的峰值压力为9.196 MPa.当采用进气总管喷水时,燃烧室内的峰值放热率为546 J·(°)-1,而采用进气歧管喷水时,燃烧室内的峰值放热率为524 J·(°)-1.从两图可以得到,与不喷水的工况相比,进气歧管喷水方案的缸内压力和放热率的下降较小,表明在相同水油比下,进气歧管喷水形成的水雾对发动机输出性能及缸内燃烧过程的负面影响更小.分析认为,造成上述现象的原因主要在于进气歧管喷水对天然气发动机进气过程充量系数的影响.

图3 不同喷水位置的缸压及放热率对比图

Fig.3 Comparison of cylinder pressure and heat released rate in different water injection positions

图4为不同喷水位置下进气质量对比图.当采用进气总管喷水方案时,进气质量为4 436 mg,当采用进气歧管喷水方案时,进气质量增加至4 513 mg.在相同水油比下,喷水位置对进气质量的影响十分显著,进气歧管喷水的进气质量比进气总管喷水进气质量大5%.此外,从图5进气至压缩过程的温度分布图可以看出,在进气及压缩过程中,进气歧管喷水的燃烧室内的温度低于进气总管喷水的燃烧室内温度,在上止点前20°曲轴转角的时候,采用进气歧管喷水方案的燃烧室内平均温度为700 K,而采用进气总管喷水方案的燃烧室内的平均温度为750 K.这表明,在低转速大负荷工况下,发动机的热负荷较高,对与进气道喷水而言,喷水位置距离燃烧室距离越近,水雾蒸发吸热降低燃烧室温度的作用越明显,从而增加进气过程的缸内充量密度.

图4 不同喷水位置进气质量对比图

2.2 喷水质量对天然气发动机的影响

图6为进气行程在不同喷水质量时湍动能及流场分布.进气行程初期,混合气被挤入燃烧室内,在上止点前344°曲轴转角时,受到气门节流作用的影响,在气门间隙处的最大湍动能仅为50 m2·s-2;活塞下行到上止点前318°曲轴转角时,活塞运动加快,缸内产生负压,同时气门处的节流损失减小,空气经螺旋气道后沿进气门面流入缸内.从图中可以看出,此时随着喷水质量的增加,进气质量增加,使得气流在进气门下方的两个滚流的湍动能增加,这有利于组织更加均匀的混合气;在上止点前263°曲轴转角时,不同滚流区之间相互影响,缸内滚流的结构发生改变.随着喷水质量的增加,右侧的滚流逐渐变大,但由于进气过程的不对称性,进气门下方的滚流区逐渐变小;随着活塞继续运动,喷水质量对湍动能的影响逐渐减弱,但是缸内仍然存在多个较小尺度的滚流区域,各个滚流区之间互相交换能量,湍动能逐渐减弱至25 m2·s-2以下,形成均匀的混合气.

a -344°

b -318°

c -263°

d -182°

e -120°

Fig.6 Distribution of turbulence kinetic energy in the intake stroke

图7为在不同喷水质量下温度场分布图,从图中可以看出,在火花塞点火之前,缸内温度分布均匀,为800 K;由于天然气的着火延迟期较长,火花塞点火后,在上止点前8°曲轴转角时,火焰在火花塞附近形成明显的火核;在燃烧上止点时,火焰处于明显发展期,以球形向外扩展,从图中可以看出,在压缩过程中缸内滚流的作用下,使得缸内已燃混合气团向进气门侧(左侧)偏移.随着燃烧过程的进行,火焰从燃烧室中心的火花塞迅速传播至燃烧室边缘的挤流区,图中温度梯度较大的分界面可以看作火焰前封面,相同时刻火焰前封面距离燃烧室中心火花塞的距离可以用来在衡量火焰传播速度.从上止点后8°曲轴转角时的对比图可以看出,随着喷水质量的增加,火焰传播速度明显减慢,从而推迟了燃烧相位.在上止点后24°曲轴转角时,火焰前封面继续传播至燃烧室的凹坑及挤流区域,随着喷水质量的增加,火焰传播速度减慢,这容易使气缸壁缝隙处残余的混合气增加,燃烧不充分,HC和CO的排放量增加.

图8为燃烧室内的平均温度及放热率曲线,从图中可以看出,随着喷水质量的增加,燃料的放热率减缓,燃烧持续期延长,燃烧室内的最高温度降低,最高温度的相位也被推迟.当没有水喷入的时候,燃烧室内的最高温度为2 490 K(上止点后15°曲轴转角);当喷水质量与天然气质量相同时,燃烧室内的最高温度为2 200 K(上止点后30°曲轴转角).

图9为不同喷水质量下燃烧室内的平均压力的变化曲线,从图中可以看出,随着喷水质量的增加,燃烧室内的峰值压力从10 MPa下降到8 MPa,峰值压力出现的相位从上止点后19°曲轴转角推迟到上止点后23.5°曲轴转角.这意味着天然气发动机采用喷水技术后,在相同的点火提前角下,发动机的做功能力会略微下降.但是,随着喷水质量的增加,火焰传播速度减慢,燃烧持续期延长,爆震倾向也会减弱,这为天然气发动机的当量比燃烧提供了可能性.

a -20°

b -8°

c 0°

d 8°

e 16°

f 24°

g 32°

h 40°

图8 不同喷水质量下缸内温度及放热率对比图

图9 不同喷水质量下缸内压力对比图

高温、富氧以及在高温富氧区的停留时间是氮氧化物生成的3个条件.图10为仿真不同喷水质量的氮氧化物的分布.从温度场分布可以看出,高温区域主要存在于已燃区和火焰前峰面处,而氮氧化物的形成区与温度场的高温区分布一致.在上止点附近的快速燃烧期,火焰传播速度急剧加快,燃烧室内火花塞附近温度达到 2 400 K 以上,有明显的氮氧化物生成迹象;在上止点后24°曲轴转角时,火焰前封面继续传播至燃烧室的凹坑及挤流区域,在温度很高的已燃区域生成了大量的氮氧化物,在没有喷水的工况下,局部最高的氮氧化物质量分数高达0.012;随着活塞的继续下行,燃烧进入后燃期,氮氧化物的生成速率降低,并总生成量维持稳定.与没有喷水的工况相比,随着喷水质量的增加,燃烧室内的高温区域明显减少,使得缸内氮氧化物的生成面积和峰值显著减小.

a 0°

b 8°

c 16°

d 24°

e 32°

f 40°

g 48°

图10 氮氧化物生成分布

图11为不同喷水质量下氮氧化物的生成量.在不改变点火时刻和空燃比的情况下,不喷水的原机工况下,单缸氮氧化物的生成量为805×10-6;当喷水质量为天然气质量的一半时,氮氧化物的生成质量下降至原机的59.94%;而喷水质量等于天然气质量时,氮氧化物的生成质量下降至原机的23.08%.通过上述仿真结果可以发现,随着喷水质量的增加,氮氧化物的成量明显降低.因此,进气总管喷水可以显著抑制燃烧室内氮氧化物的生成,从而大幅缓解稀燃工况下天然气发动机对后处理装置的依赖.

图11 不同喷水质量下氮氧化物的生成量

2.3 喷水对反应过程的影响

为了进一步研究喷水对天然气发动机化学反应过程的影响,在当量比为1.21,喷水和燃气质量比为0.5,1 100 r·min-1全负荷工况下,对比了喷水和未喷水情况下几种物质浓度和温度曲线的关系.在火焰前锋之前的末端气体中产生CH2O和OH自由基可以作为燃烧反应和已燃区的标记[17].如图12所示,在温度低于1 500 K的低温反应阶段,随着燃料的消耗,在两种情况下的温度均缓慢升高,CH2O和H2O2中间自由基在此阶段都显著上升,且无明显差别.在低温反应阶段主要是脱氢反应,即CH4被自由基(H,O和OH等)氧化成CH3,CH3与氧原子结合并被氧化成CH2O.仿真结果表明,喷水对低温反应阶段CH2O和H2O2的生成没有显著影响.

在温度高于1 500 K的高温反应阶段,喷水工况的OH和HO2自由基浓度明显低于未喷水工况的浓度,这些变化主要涉及以下反应:

H2O2+ (M) = OH + OH + (M)

(R1)

CH2O + 2OH = CO + 2H2O

(R2)

HCO + OH = CO + H2O

(R3)

CO + OH = CO2+ H

(R4)

如图12所示,当燃烧室内的温度高达1 500 K时,CH2O和H2O2的浓度上升到相对较高的水平.此时H2O2开始分解产生大量OH自由基,这也是高温反应的开始指标.同时,R2反应消耗OH和CH2O形成CO,而CO进一步被氧化成CO2放出大量的热量.当喷水以后,R1反应被H2O抑制,导致生成的OH浓度降低,因此燃烧持续期被延长.而R3反应速率(指前因子A= 50×1013)比R4反应速率快得多(指前因子A= 4.76×107),这导致没有足够的OH来进行R4反应,因此喷水以后的最高温度明显降低.因此,在高温反应阶段,喷水以后,OH自由基的生成速率被抑制,致使没有足够的氧化剂氧化CO,放热持续期延长,反应温度降低.