汽车工业的快速发展,在带来工业文明便利与经济蓬勃发展的同时,也带来了能源问题与环境污染.内燃机作为当前车用主要动力源,其包括新型工作循环等在内的节能减排技术研发迫在眉睫.受火力发电站行业成熟运用的清洁能源系统(clean energy system)的启发,著名内燃机学家Robert W Bilger教授于1999年在第五届清洁环境技术及燃烧会议首次提出内燃兰金循环(internal combustion rankine cycle,ICRC)概念[1].该循环的主要特点是采用纯氧代替空气作为助燃剂与碳氢燃料反应,由于燃烧过程中没有氮气的参与从而抑制NOx的形成,通过水作为介质回收尾气废热,控制缸内燃烧温度,水汽化膨胀做功实现高效燃烧和更好的燃油经济性.
同济大学吴志军教授与Robert W Bilger教授合作,使用Stanjan软件计算了不同边界条件下理论循环热效率,证明了喷入缸内的水在循环过程中向系统额外加入工质,提高了缸内做功工质质量[2].理论计算表明,随着喷水质量及喷水温度的增加,ICRC发动机热效率增加.结合零维与三维计算研究,基于一台125 mL单缸汽油发动机,在世界范围内首次完成内燃兰金循环原型机开发,并针对内燃兰金循环开展大量试验研究[3-8].试验研究表明,缸内喷入的高温高压水能够快速蒸发汽化提高缸内工质质量,使得膨胀冲程缸内压力增加,增加示功图面积,提高循环热效率,与未喷水循环工况相比,实际喷水循环工况的指示热效率提高了7.8%.
内燃兰金循环采用纯氧燃烧在规避NOx的同时,也使得缸内燃烧速度提高,缸内压力升高率增加,燃烧过程爆震倾向加剧.与此同时,理论计算证明,提高压缩比可以进一步发掘内燃兰金循环热效率提升潜力,在喷水质量为60 mg时,将压缩比由9提高至14, 可以提高热效率15.7%[9].但高压缩比会使缸内温度和缸内压力升高率增加,同样加剧燃烧过程爆震倾向.因此高压缩比与纯氧燃烧对燃烧过程的控制与优化提出了更高要求.内燃兰金循环在燃烧过程中喷入的高温高压水不仅可以作为额外工质,在一定程度上还抑制燃烧过程的剧烈程度,因此,研究通过缸内喷水技术进行内燃兰金循环发动机爆震控制显得尤为重要.
1 试验台架与方案本文基于一台改造的直列双缸柴油机开展内燃兰金循环爆震工况控制研究.针对使用需求,试验中对柴油机的第一缸进行改造,主要包括进气供给系统、燃料供给系统、高温高压水供给系统,以实现进气氧浓度实时调节、进气道燃料喷射及缸内高温高压水喷射.本文所使用的试验发动机参数如表 1所示,试验系统图见图 1.本文选择丙烷作为燃料开展试验研究,与95号汽油相比,丙烷的研究法辛烷值(research octane number, RON)高达110,其抗爆震能力更强,更有利于内燃兰金循环发动机发挥高热效率的优势.
通过基于单片机自行开发的控制系统,可以对包括喷油量、喷油定时、喷水量、喷水定时、节气门、点火脉宽和点火定时在内的发动机工作参数进行灵活稳定控制.采用瑞士Kislter公司6052C型号的压力传感器和5064B11型号的电荷放大器检测缸内压力变化,并通过美国国家仪器公司的PCI6250信号采集卡结合自行开发的信号采集系统对工作过程中的缸内压力及进排气系统相关参数进行采集.该发动机喷水系统由本文作者自行开发,主要由喷水喷嘴、高压水加热共轨管、气液增压泵3部分组成.高压水加热共轨管可以在不高于45 MPa的压力下将待喷射水加热至最高160 ℃,气液增压泵增压比为50:1,可以在外部输出压力0.7 MPa的情况下将待喷射水由常压加压至35 MPa.缸内高温高压水供给系统示意图如图 2所示.缸内喷水喷嘴布置示意图如图 3所示.
试验中采用的试验工况参数如表 2所示.为有效评估试验过程中的爆震工况,选取爆震特征循环进行单一循环分析,并利用爆震强度指标(knock intensity, KIA)进行评价.该指标的获取是对试验采样的缸内压力转角图使用快速傅里叶变化(fast Fourier transformation, FFT)进行带通滤波,以消除测量过程中混入的压力谐波的影响后,取滤波后压力震荡的最大幅值作为爆震强度的评价指标.该指标能够体现爆震发生时缸内冲击波的强度[10].为保证试验结果是进气氧浓度对爆震倾向的单一影响因素,在试验过程中维持进气温度为25 ℃,以排除因进气温度改变而带来的缸内压缩终了时的温度改变,从而导致的爆震倾向变化.
为描述不同工况时的燃烧过程,基于热力学第一定律,由下式计算缸内燃烧瞬时放热率:
$\frac{{{\rm{d}}Q}}{{{\rm{d}}\theta }} = p \cdot \frac{\kappa }{{\kappa - 1}} \cdot \frac{{{\rm{d}}V}}{{{\rm{d}}\theta }} + V \cdot \frac{1}{{\kappa - 1}} \cdot \frac{{{\rm{d}}p}}{{{\rm{d}}\theta }}$ | (1) |
${Q_{{\rm{cum}}}} = \smallint {\rm{d}}Q = \smallint p \cdot \frac{\kappa }{{\kappa - 1}} \cdot V{\rm{d}}p$ | (2) |
式中:Q为综合放热量;Qcum为累计放热量;p为缸内压力;θ为曲轴转角;κ为比热比;V为气缸容积.同时,采用平均指示压力变动系数γIMEP来度量燃烧的循环变动的程度,其中γIMEP的定义如下式:
${\gamma _{{\rm{IMEP}}}} = \frac{{{\sigma _{{\rm{IMEP}}}}}}{\chi } \times 100\% $ | (3) |
式中:σIMEP为平均指示压力的标准偏差;χ为平均指示压力的平均值.
2 试验结果与分析 2.1 未喷水工况下爆震及燃烧过程分析未喷水工况下,调节发动机压缩比11.6,节气门开度15%,点火时刻对应曲轴转角335°,进气分别为空气和进气氧气体积分数(oxygen fraction, OF)55%的O2和CO2混合气,循环供油量为52.7 mg,转速1 500 r·min-1.
图 4给出了空气和OF55的缸内压力、温度及示功图曲线.由图可知,与助燃剂为空气相比,OF55工况下缸内压力最大值显著提高,点火后混合气滞燃期缩短,缸压升高率较高,且缸内最高温度出现相位显著提前,幅值略有提升,表征燃烧速度得到明显提高.同时,OF55工况缸压曲线在压缩上止点后1°曲轴转角左右起出现了频率较高、幅度较大的锯齿状波动,说明该工况燃烧较为粗暴.对其缸压进行频谱分析和高通滤波,根据爆震信号一般所在频率范围为6~8 kHz[11],滤掉4 kHz以下部分,得到高频压力如图 5所示.从图中可以看到,OF55工况的KIA达到了0.287 MPa,由此可以判定OF55工况发生了爆震.
发生爆震的原因在于末端混合气的快速自燃,在正常火焰到达前,末端混合气内部出现一个或多个火焰中心,从这些中心以100~300 m·s-1直到800~1 000 m·s-1或更高速率传播火焰,迅速将末端混合气燃烧完[12].本文工况出现爆震的主要原因在于进气中氧浓度较高、混合气较浓以及发动机运行转速较低所致,进气中较高的氧浓度使得混合气起始反应所需温度降低,高氧浓度下发动机长时间运行造成机体温度升高,燃烧室壁面温度较高,过热的壁面加快了末端混合气的自燃.此外,燃烧室形状为ω形的非紧凑型燃烧室,造成了火花塞点火后火焰传播距离增加,同样促进了爆震的发生.
图 6给出了空气和OF55工况的放热率曲线.表 3列出了该工况下燃烧和放热参数.表中, pmax为缸内最大爆发压力, θpmax为缸内最大爆发压力对应和曲轴转角, dpmax为最大压力升高率, θdpmax为最大压力升高率对应曲轴转角, dQmax为放热率极值, θdQmax为放热率极值对应曲轴转角, CA10、CA50分别为10%和50%累计放热量对应曲轴转角.由表 3的燃烧和放热参数可知,空气、OF55工况的瞬时放热率峰值分别为32.0 J·(°)-1、35.4 J·(°)-1,对应相位的曲轴转角分别为364°和354°,同时OF55工况的CA10、CA50时刻均早于空气工况,出现这种现象主要是因为进气中氧浓度的提高降低了气缸中非助燃气体量,同时氧浓度的提高促进了丙烷的燃烧化学反应速率,使得点火后CA10时刻较早,缩短了滞燃期,加快了混合气的火焰传播速率,提高了最高燃烧压力,使得缸压升高率峰值增加,燃烧相位前移.由图 4b和表 3可知,相比于空气,OF55爆震工况的最高燃烧压力的提高,增加了示功图的有效面积,OF55工况的指示功增加了75 J,平均指示压力(indicated mean effective pressure, IMEP)增加了0.1 MPa,相同循环供油量下指示功的增加,使得指示热效率提高了3.1%.
综上所述,相比于空气作为助燃剂,高压缩比下高氧浓度工况混合气的燃烧速率加快,燃烧相位前移,同时最高燃烧压力提高,增加了指示功的输出,提高了循环指示热效率.
图 7给出了空气和OF55工况下IMEP的散点图.由图可知,OF55工况的IMEP较高,但循环间波动较大; 而以空气作为助燃剂工况的IMEP整体比较平缓,循环间波动较小.分别计算其变动系数γIMEP,得出空气作为助燃剂时γIMEP为0.009,而OF55工况的γIMEP达到了0.019,氧浓度增加后循环变动系数增加了0.01,说明了高压缩比下氧浓度的提高,增加了IMEP,同时使得燃烧循环变动增加,燃烧稳定性变差.
喷水工况下,保持发动机参数与未喷水工况一致,调节循环喷水量194 mg,喷水压力35 MPa,喷水温度120 ℃,进气OF值55%.
图 8给出了喷水前后的缸压及缸压升高率随曲轴转角的变化曲线.由图可知, 未喷水循环发生了爆震,而喷水后缸压曲线较为圆滑,缸压升高率较为平缓.图 9给出了喷水前后的爆震强度KIA值,分析两种工况下KIA爆震强度可知,未喷水工况的KIA达到了0.287 MPa,而喷水循环的KIA则近似为0 MPa,说明了喷水后爆震消失,由此说明喷水能够对爆震工况进行控制.
表 4列出了爆震工况的燃烧参数,未喷水循环的缸压最大值为6.1 MPa,出现在364°曲轴转角,缸压升高率峰值达到了0.27 MPa·(°)-1,而喷水后缸压最大值降为5.8 MPa,对应相位的曲轴转角后移至369.5°,同时缸压升高率峰值降低为0.19 MPa·(°)-1,说明了喷入缸内的水能够汽化蒸发吸收爆震工况缸内过高的温度,控制缸压升高率及燃烧相位,同时膨胀做功,在此工况下喷水后指示功增加了19 J,提高了2.2%.
图 10为喷水对爆震工况下放热率及累计放热量随曲轴转角变化曲线.计算过程没有排除喷水过程对缸内做功工质增加的影响,故喷水工况的放热计算曲线可看作是缸内“综合放热”,即结合了燃料燃烧放热和高温水汽化膨胀两者对缸内压力变化的影响,这样的计算方法在简化计算过程的同时可以更为方便、直观地分析喷水过程对于燃烧放热的影响.如果喷入高温水汽化膨胀所引起的缸内压力升高量不能补偿由于其汽化吸热及降低混合气反应速率所引起的压力降低,则其综合放热特性将表现为瞬时放热率及累计放热量, 与未喷水工况相比有所降低;反之, 则会使得放热曲线有所提升.
由于爆震的影响,未喷水工况的瞬时放热率在压缩上止点后出现了较大的波动,这是由锯齿状缸压的影响,造成放热不稳定,喷水循环瞬时放热率相位出现后移,波动变得平缓,说明喷水降低了发动机工作的粗暴程度,利用爆震造成的过高的缸内温度迅速汽化了高温高压水,此时水蒸气变成了做功工质,降低了爆震所带来的影响.图 11给出了爆震工况下喷水前后IMEP变化的散点图.可以看出, 喷水开始前的循环IMEP波动严重,图中圆圈标记了较为明显的爆震循环,其明显降低了IMEP,加剧了燃烧的不稳定性,喷水后IMEP有所增加且循环波动较小,喷水使得IMEP的循环变动系数降低了0.006.在喷水结束后的几个循环内IMEP波动变缓,但从第80个循环开始(第二个圆圈中标记)又出现了较大的波动,说明了喷水结束后的若干循环后燃烧开始变差,可见喷水对于爆震工况下的燃烧稳定性具有良好的控制作用.
表 5给出了不同边界条件下喷水对爆震的控制情况.表中工况为氧气体积比例60%,点火时刻对应曲轴转角为330°,发动机转速为1 500 r·min-1,压缩比为10.5,进气压力为0.1 MPa,节气门开度为0%,喷水压力为35 MPa,喷水温度为120 ℃.当喷油量从46.7 mg提高至49.7 mg时,未喷水工况下的爆震概率提高,由6%增加至16%.在各个工况下,喷水都可以有效地抑制爆震,使得爆震概率进一步降低.推迟喷水时刻可以使得排气及缸内残余废气温度降低,从而使得下一个循环缸内温度降低,使得爆震能得到进一步的控制.喷水量的增加也可以有效地降低排气及缸内残余废气, 从而降低缸内温度,更好地控制爆震的产生.
(1) 未喷水工况下,与空气进气相比,OF55进气时缸内燃烧过程产生爆震,KIA值达到0.26 MPa,其生成原因主要是进气氧浓度较高, 混合气较浓, 以及发动机运行转速较低.
(2) 喷水工况下,OF55进气下缸内燃烧爆震消除,KIA值降低至接近0 MPa,说明缸内高温、高压水喷射能够有效控制内燃兰金循环发动机在高氧浓度工况下的爆震倾向.
(3) 与未喷水工况相比,喷水工况下缸内最高爆发压力降低但推迟,缸压压力升高率最大值降低,循环指示功增加,说明喷入缸内的高温高压水能够快速吸热蒸发以提高循环热效率.
[1] |
BILGER R W. Zero release combustion technologies and the oxygen economy[C]//5th International Conference on Technologies and Combustion for a Clean Environment. Lisbon: [s.n.], 1999: 1039-1046.
|
[2] |
BILGER R W, WU Z. Carbon capture for automobiles using internal combustion Rankine cycle engines[J]. Journal of Engineering for Gas Turbines and Power, 2009, 131(3): 034502 DOI:10.1115/1.3077657 |
[3] |
吴志军, 于潇. 基于内燃兰金循环的二氧化碳回收车用动力系统[J]. 吉林大学学报(工学版), 2010, 40(5): 1199 WU Zhijur, YU Xiao. CO2 capture automotive engine system based on internal combustion Rankine cycle[J]. Journal of Jilin University(Engineering and Technology Edition), 2010, 40(5): 1199 |
[4] |
YU X, WU Z J. Simulation on effect of EGR on oxy-fuel IC engine[C]//Applied Mechanics and Materials. [S.l.]: Trans Tech Publications, 2012: 790-795.
|
[5] |
WU Z, YU X, FU L, et al. Experimental study of the effect of water injection on the cycle performance of an internal-combustion Rankine cycle engine[J]. Journal of Automobile Engineering, 2014, 228(5): 580 DOI:10.1177/0954407013511069 |
[6] |
付乐中, 于潇, 邓俊, 等. 内燃兰金循环发动机试验系统开发[J]. 内燃机工程, 2013, 35(6): 87 FU Lezhong, YU Xiao, DENG Jun, et al. Development of internal combustion Rankine cycle engine test system[J]. Chinese Internal Combustion Engine Engineering, 2013, 35(6): 87 |
[7] |
WU Z, FU L, GAO Y, et al. Thermal efficiency boundary analysis of an internal combustion Rankine cycle engine[J]. Energy, 2016, 94: 38 DOI:10.1016/j.energy.2015.10.099 |
[8] |
FU L, WU Z, LI L, et al. Effect of water injection temperature on characteristics of combustion and emissions for internal combustion Rankine cycle engine[C]// SAE International World Congress. Detroit: SAE, 2014: 140126.
|
[9] |
WU Z, YU X, FU L, et al. A high efficiency Oxy-fuel internal combustion engine cycle with water direct injection for waste heat recovery[J]. Energy, 2014, 70: 110 DOI:10.1016/j.energy.2014.03.095 |
[10] |
CHUN K M, HEYWOOD J B. Characteristization of knock in an S.I. engine[C]// SAE International World Congress. Detroit: SAE, 1989: 890156.
|
[11] |
HEYWOOD J B. Internal combustion engine fundamentals[M]. New York: McGraw-Hill, 1988
|
[12] |
PAN J, SHEPPARD C. A Theoretical and experimental study of the modes of end gas auto-ignition leading to knock in S. I. Engines[C]// SAE International Fuels & Lubricants Meeting. Baltimore: SAE, 1994: 942060.
|