Mechanism of in-cylinder turbulence on the distribution of fuel activity in hybrid combustion
-
摘要: 缸内受限条件下燃料与湍流的相互作用是燃料分层控制复合燃烧的关键问题。针对该问题,通过向缸内直喷高活性燃料二甲醚(Dimethyl ether,DME),形成高活性燃料浓度分层。基于光学可视化发动机实验平台,利用粒子图像测速(Particle Image Velocimetry,PIV)、Rayleigh散射、Mie散射以及高速摄影结合放热分析等手段对复合燃烧这一缸内受限空间下的流动及燃烧过程进行了观测,并通过三维计算流体力学(Computational Fluid Dynamics,CFD)仿真手段对观测到的现象进行解释。结果表明:缸内存在大范围逆时针涡流场,DME的蒸发和扩散过程受到流场的作用;在流场的作用下,缸内燃烧过程呈现DME集聚区域自燃-火焰传播-多点自燃放热特征。Abstract: The interaction between fuel and turbulence under the in-cylinder limited conditions is the key issue for hybrid combustion controlled by fuel activity stratification. Dimethyl ether (DME) is injected to the cylinder to produce high activity fuel stratification. Particle image velocimetry, laser Rayleigh scattering, Mie scattering and high speed imaging combined with heat release analysis on the optical engine experiment platform are used to observe the flow field and combustion process of hybrid combustion in the limited space of cylinder. 3D Computational Fluid Dynamics (CFD) simulation are used to explain the experimental phenomena. Result shows that there is a large range of counter-clockwise vortex field in the cylinder, and the diffusion and evaporation process of DME is influenced by the flow. Under the flow field, the combustion process in the cylinder shows characteristics of DME auto-ignition in the distribution area, flame propagation, multi-point auto-ignition.
-
Keywords:
- dimethyl ether /
- micro flame /
- hybrid combustion /
- visualization /
- 3D-simulation
-
0 引言
作为汽车的主要动力源之一, 汽油发动机具有动力响应性好、升功率高、成本低等优点, 但是其燃油经济性受到汽油特性的限制, 相比柴油机要低20%~ 30%.针对该问题, 学界提出火花点火-可控自燃复合燃烧(Spark ignition-Controlled auto-ignition Hybrid Combution, SCHC)这一概念, 将汽油和柴油发动机可控自燃(Controlled Auto-Ignition, CAI)的燃烧方式加以组合, 不仅可以实现汽油机均质压燃, 还可以通过火花点火的火焰传播对CAI燃烧进行有效调控, 减轻中高负荷下CAI燃烧的爆震问题, 拓宽CAI燃烧的运行极限.但是, 随着研究的深入, SCHC燃烧暴露出高废气率条件下的平均指示有效压力(Indicated Mean Effective Pressure, IMEP)循环变动以及可控性差等问题[1-2].研究表明, 缸内湍流对初期的火核形成及火焰传播起到重要作用, 被认为是引起循环变动的重要原因之一[3-4]。为解决该问题, 本课题组提出二甲醚(Dimethyl ether, DME)微火源引燃复合燃烧概念[5], 利用高活性燃料二甲醚自燃形成高活性燃料分层, 形成微火源并引燃油气混合气进而引发后期自燃。微火源可以增强高废气稀释条件下火核形成及前期火焰传播过程, 提高燃烧稳定性, 从而实现发动机热效率的提高。
韩国汉阳大学Junepyo Cha等[6]对比了DME-汽油双燃料与DME单燃料下可控自燃(CAI)燃烧过程。结果表明, DME喷射时刻可以直接控制着火时刻, 当着火时刻在上止点之前时, 采用DME-汽油双燃料具有更好的经济性。Brunel大学的Zhang等[7]在一台单缸汽油机上同样利用DME-汽油双燃料喷射技术实现了CAI燃烧。研究结果表明, 不同的DME喷射策略下的燃烧特征呈现较大差异。当DME单次早喷时, DME在气缸内与汽油混合气混合较为均匀, 因此整缸混合气活性分布较为一致, 燃烧过程呈现CAI燃烧特征; 而当DME单次晚喷时, 燃烧过程呈现火焰传播的特征。徐康[8]在单缸光学发动机上开展了DME微火源引燃复合燃烧的可视化实验。结果表明, 微火源的分布受DME分布影响。微火源引燃复合燃烧的初始反应区位置(微火源)呈现集中分布的特征。但是, 已有的实验研究对高活性燃料DME在发动机缸内这一受限空间受流动影响的扩散和分布机制尚不明确, 仅通过实验研究无法对缸内微火源产生及引燃汽油的作用机制进行解释。因此, 需要通过实验结合仿真手段对实验现象中的燃烧过程进行分析。
本文基于光学发动机平台, 通过粒子图像测速(Particle Image Velocimetry, PIV)手段探究微火源引燃特征最为明显的单次晚喷工况下缸内流场分布, 并对其成因进行分析。通过Rayleigh散射测量缸内DME分布, 研究缸内复杂的流动条件对DME活性分布的影响。同时, 通过高速摄影结合三维计算流体力学(Computational Fluid Dynamics, CFD)仿真手段对DME分布及流动特征对燃烧放热过程影响的机理进行解释。
1 实验及仿真设置
1.1 光学发动机实验平台
实验均在单缸可视化光学发动机实验平台上完成。光学发动机实验平台系统如图 1所示, 其中可视化窗口的直径为71mm。在光学测量系统方面, 配备4倍频高能钇铝石榴石Nd:YAG (Neodymium-Yttrium-Aluminium Garnet)激光器、双脉冲532nm波长的Nd:YAG激光器、增强型光电耦合探测器(Intensified Charge-Coupled Device, ICCD)和高速摄影机。光学发动机主要参数如表 1所示。发动机采用的进气道及燃烧室设计可以使进气在缸内形成涡核, 从而增强缸内涡流运动。发动机配备了不同升程和持续期的凸轮轴系, 可以根据实验需求灵活调整进排气升程, 并通过机械错齿改变气门相位。由于光学发动机较低的压缩比, 采用辛烷值较高的汽油燃料难以在光学发动机上实现混合气自燃, 因此, 实验中采用的气道喷射燃料为辛烷值较低的基础燃料(Primary Reference Fuel, PRF)PRF40(PRF40表示在基础燃料中, 异辛烷的占比为40%, 余数为正庚烷的比重)。
表 1 光学发动机主要参数Table 1 Specifications of the optical engine发动机形式 四冲程单缸汽油机 缸径 95mm 冲程 95mm 气门升程及相位 可变 排量 0.67L 压缩比 9.24 燃烧室结构 棚顶室 燃油喷射方式 气道喷射+缸内直喷 气道喷射燃料 PRF40 气道喷油压力 300kPa 直喷燃料 DME 直喷喷油压力 4MPa 进气方式 自然吸气 发动机转速 600r/min 节气门 有节气门 由于DME微火源引燃要求较低的DME循环喷射量, 因此DME油轨内压力不宜过高。综合考虑涉及工况点的DME喷射背压, 在实验中, DME共轨压力保持为4MPa。DME直喷喷油器采用的是6孔直喷汽油机喷油器, 单个喷孔直径为0.2mm, 该喷油器采用侧置式布置方案安装于进气道的正下方, 喷油器轴线与气缸横截面呈15°夹角。
由于光学发动机较短的运行时间, 常规的外部废气再循环(Exhaust Gas Recirculation, EGR)系统难以从排气中将废气引入进气系统中。因此, 实验中使用氮气(N2)和二氧化碳(CO2)按照一定比例混合来模拟外部EGR。其中, N2和CO2的比例根据实验中所采用的燃料化学式, 在当量燃空比和完全燃烧的假设下计算得出。虽然在模拟的外部EGR中不包含燃烧产物H2O, 且完全燃烧的假设也忽略了外部EGR中微量活性小分子物质的存在(这会对结果产生一定的影响), 但是, 本课题组之前的燃烧采样分析结果显示[9], 外部EGR中小分子组分的变化对燃烧相位并无显著影响。因此, 实验中对模拟外部EGR组分的简化假设并不会对结果产生明显影响。
1.2 CFD仿真平台介绍
使用由美国Convergent Science公司开发的计算流体力学软件Converge构建了发动机三维模型, 如图 2所示(图中A, B为截面位置)。其中, 发动机内流场包含燃烧室和进、排气道等区域。计算起始时间为进气门开启前10°曲轴转角(Crank Angle, 后文均以℃A表示), 计算结束时间为放热终止时刻(压缩上止点后40℃A)。在整个计算过程中, 最大网格尺寸为基础网格尺寸4mm, 最小网格尺寸为自适应网格加密极值0.125mm, 网格数量为120~440万。
在模拟缸内直喷DME过程时, 为了更好地模拟液滴撞壁现象, 采用Wall Film[10]模型, 该模型可模拟燃油撞壁后的黏附、反弹、摊布和飞溅等现象; 液滴之间的相互作用采用NTC Collision模型[11]进行模拟; 使用KH-RT Breakup Length模型[12]可以较好地模拟液滴的雾化破碎过程, 其优势在于可以较准确地预测油束的贯穿度; 而Frossling Drop Evaporation模型耦合Boiling模型[13]的使用可以准确模拟DME的蒸发过程。为验证所选用模型的准确性, 基于定容弹实验平台对喷雾模型进行了标定。实验与仿真贯穿距测量验证结果如图 3所示。结果表明:采用的DME喷雾和蒸发模型喷雾贯穿距、喷雾形态及锥角能够较好地模拟DME的发展情况, 满足本文研究的需求。
在DME-PRF复合燃烧的化学反应动力学机理方面, 基于Cai和Pitsch等[14]构建的PRF机理和Burke等[15]构建的DME化学反应动力学机理, 构建了包含348种组分和1125步反应的DME-PRF详细化学反应动力学机理, 并基于雅可比辅助的直接关系图(Direct Relation Graph with Error Propagation, DRGEP)方法[16-17], 结合目标检索算法(Targeted Search Algorithm, TSA)、准稳态假设(Quasi-Steady State Assumption, QSSA)和氢离子敏感性分析对反应机理进行了简化, 构建包括了143种组分、729步反应的DME-PRF骨架反应机理。详细机理的点火延迟和火焰传播速度如图 4所示(相关信息请参考文献[21-26]), 详细机理与骨架机理的对比如图 5所示。结果表明, 构建的DME-PRF机理与实验验证结果吻合较好, 可以用于本文所进行的CFD仿真研究。
2 实验及仿真结果分析
2.1 缸内流动特征分析
针对发动机燃烧室内的复杂流动条件, 采用LaVison开发的FlowMaster PIV系统测量缸内流动特征, 如图 6所示。该系统以2台SOLO120 Nd:YAG激光器作为光源, 以Laskin喷雾器将癸二酸二异辛酯溶液以300~500kPa的喷射压力喷入气流中。采用CCD照相机作为传感器进行图像采集, 并使用DaVis 7.2软件进行图像处理。根据之前的研究结果[18], 选取上止点前20°曲轴转角(Crank Angle Before Top Dead Center, ℃A BTDC)作为测量时刻, 此时缸内流场分布会对DME喷雾扩散及后续的燃烧过程产生显著影响, 且在该时刻活塞上行对测量没有影响。
除DME喷射外, 本研究所采用的发动机运行工况与PIV实验工况条件均一致, 如表 2所示。众所周知, 火花点火发动机的缸内流动主要受进气和压缩过程中缸内滚流的影响。为了更好地形成缸内涡流, 对采用的光学发动机进气道及燃烧室进行了增强缸内涡流的设计。随着活塞上行, 滚流在燃烧室及进气的影响下逐渐转化为涡流。实验测量了100个循环, 图 7展示了其中的8个随机循环的缸内流场及其涡心分布, 而图 8显示了该工况下PIV测量获得的100个循环的平均流场。缸内流场从进气侧开始, 向光学视窗的右上流动, 并沿缸壁从窗口左侧流出, 以窗口中心偏左侧为中心形成了涡流结构。
表 2 发动机运行工况Table 2 Engine operation setup发动机转速/(r·min-1) 600 排气门开启时刻[CA]/(°) 180 排气门关闭时刻[CA]/(°) 328 排气门升程/mm 2 进气门开启时刻[CA]/(°) 396 进气门开启时刻[CA]/(°) 508 进气门升程/mm 1 循环气道喷油量/(mg·cycle-1) 14.8 气道喷油时刻[CA]/(°) 330 循环DME直喷油量/(mg·cycle-1) 1.8 DME直喷时刻(Start of Injection, SOI)[CA BTDC]/(°) 25 过量空气系数(lambda) 1 外部EGR率 27% 进气温度/℃ 40±1 冷却水温/℃ 85±2 三维CFD仿真所得到的缸内流场形成过程如图 9所示。进气结束时刻(200℃A BTDC)缸内主要流动形式为逆时针滚流, 其滚流比为1.1。同时, 在缸内进排气中轴线一侧出现约4mm的涡流核心, 该核心位置由进气道及燃烧室形状共同决定。随着活塞运行, 在燃烧室形状的影响下, 缸内滚流逐渐减弱, 滚流涡核逐渐消失, 滚流在壁面和初始涡核引导下逐渐转化为涡流, 缸内涡流比逐渐提升至0.8, 形成以涡流为主的流动特征。
针对所测得的PIV流动图像及100个循环的平均流场, 对缸内流动的循环变动进行了分析。其中, 通过计算每循环缸内湍动能与平均流场的标准差来计算缸内流场的循环变动(COV)。针对缸内湍动能, 首先通过二维快速傅里叶变换(Fast Fourier Transform, FTT)将二维瞬时速度变换为空间频率域, 然后再将高于给定空间截止频率的所有傅里叶系数设置为0后执行逆FTT。截止频率通过每个测量周期的速度波动分量的功率谱密度(Power Spectral Density, PSD)确定[19-20]。缸内湍动能及其循环变动如表 3所示。结果表明, 低频湍动能约占总湍动能的90%, 此时缸内涡流以大涡为主, 同时伴随着部分小尺度涡流在废气中耗散。结合PIV分析结果, 高频湍动能的循环变动表明缸内湍流循环变动主要是由小尺度涡核位置的循环变动产生的。
表 3 缸内湍动能及其循环变动Table 3 Turbulent kinetic energy and COV in cylinder湍动能/(m2·s-2) 循环变动(COV)/% 总湍动能 2.12275 23.20492 高频湍动能 0.21347 26.52597 低频湍动能 1.86785 25.39271 2.2 缸内流动对DME分布的影响
缸内DME分布及扩散过程受到流场的调控作用。由于实验条件限制, 无法在DME喷射时利用PIV手段测量缸内流场。故通过Mie散射及Rayleigh散射对DME喷雾进行测量, 并使用仿真手段解释流动对DME分布的影响。其中, Mie散射是入射光波长远小于测量介质粒子直径时出现的一种散射现象, 主要用于测量液相喷雾形态; 而Rayleigh散射在前期液相Mie散射较强时被Mie散射所覆盖, 因此主要用于测量小液滴及气相信号。图 10展示了喷雾发展历程中50个循环的平均值。在喷射初期DME未完全雾化时, 测量结果中会同时存在Rayleigh散射和Mie散射, 但由于Mie散射信号较强, 在直喷开始时刻后(After Start of Injection, ASOI)1.44℃A和2.16℃A的拍摄角度下, 可能存在Rayleigh散射信号被Mie散射所掩盖的情况。从图中只能观测到液相喷雾的形态, 而在2.88℃A ASOI时, 由于DME迅速雾化, 此时Mie散射信号大幅削弱, 已经能从图中观测到部分散射光信号由强至弱的过渡现象。在4.32℃A ASOI后(对应图 7中的结果), 散射光信号图像中已经观测不到明显的局部高亮度区, 而呈现连续的均匀过渡, 表明此时DME已经基本雾化完毕, Mie散射信号基本消失, 此时测量平面内采集到的是入射光在气相的DME中产生的Rayleigh散射信号。
三维CFD仿真得到的5℃A ASOI缸内DME分布及流场分布如图 11所示。其中缸内整体流场形态仍为逆时针涡流, 同时, 在DME喷雾的影响下, 大涡核心逐渐向气缸的左侧移动, 使缸内左侧部分的流动增强。此时缸内DME喷雾受到缸内逆时针涡流影响而形成逆时针扩散, 而缸内涡核位置分布在DME喷雾外围, 影响DME边缘的扩散及其浓度梯度, 但不会对DME喷雾的整体形态产生影响。二者共同作用使DME移动到进排气中轴线位置形成集聚型分布, 此时集聚区域DME浓度与其他区域的比值为5:1。
2.3 缸内流动及DME分布对燃烧过程的影响
在所选工况下, DME微火源引燃复合燃烧呈现多阶段放热的特征。图 12中展示了光学发动机实验的100个循环平均缸压、放热率以及三维仿真获得的缸压、放热率及缸内温度历程数据。主燃烧放热过程可分为3个主要阶段。第一阶段为剧烈放热阶段, 燃烧放热率在不到5℃A的时间内超过30J/℃A; 接下来的第二阶段放热率小幅下降并出现平台期, 持续时间略短于第一阶段的快速放热; 而在主燃烧放热的第三阶段, 放热率曲线斜率猛增, 缸内出现迅猛的燃烧放热反应, 此时放热达到峰值。通过多阶段放热, 可以对缸内燃烧的放热率及最大压升率进行调控, 从而抑制燃烧过程中的爆震现象, 提高缸内燃烧的稳定性。
为进一步了解DME、MFI复合燃烧的燃烧放热过程, 使用高速摄影采集的数据对放热过程进行分析, 如图 13所示。由于其循环变动特征十分微弱, 因此仅选择其中一个循环的放热过程进行分析。缸内首次观测到着火点的时刻为1.4℃A BTDC, 初始自燃点首先出现在可视区域的中部偏左区域, 但是该阶段的自燃反应强度较弱, 并伴随着较弱的自燃化学发光。随后, 自燃着火面积迅速扩大, 逐渐向视窗的右下方发展。从0.7℃A ATDC之后, 自燃着火的放热逐渐减缓, 此时从高速摄影的结果可以看出火焰锋面正在逐渐成型, 在2.2℃A ATDC时, 已经可以看到一个带着明亮反应锋面的火焰结构形成, 并在接下来的一段时间内, 火焰传播覆盖区域缓慢扩大, 火焰锋面的亮度也有一定程度的增加, 但是由于较慢的火焰传播速度和活塞的下行, 导致了此阶段燃烧放热率的小幅度下降, 放热率出现平台区域。在5.8℃A ATDC附近, 未燃区的PRF40燃料达到自燃条件, 开始自燃。此时, 图中火焰锋面外侧(视窗右上角)出现了十分明亮的着火点, 其亮度远高于火焰传播覆盖的区域。随后外围的混合气迅速自燃, 并覆盖几乎整个火焰锋面以外的区域。通过上述分析可知, DME微火源引燃复合燃烧的主放热阶段呈现自燃-火焰传播-自燃的3阶段顺序放热特征。
三维CFD仿真获得的缸内DME分布及火焰锋面的结果如图 14所示。其中以1000K等温面来体现火焰锋面发展过程。结果表明, 在2℃A BTDC时, 初始火核在缸内高活性燃料DME集聚分布区域的中心产生, 随后向未燃区进行火焰传播, 其初始自燃区域与DME集聚区域吻合。可以看到, 在2℃A ATDC时刻放热量为总放热量的37.5%, 原高活性区域内的DME几乎已经被完全消耗, 高活性燃料DME自燃起到高能点火源的作用。随后火焰向燃烧室外围传播, 但传播速度减缓, 同时, 火焰传播使缸内的温度从850K提高到1250K。5℃A ATDC时刻, 较高的缸内温度诱发了已燃反应区外围的自燃, 形成如图 11所示的多阶段放热特征。
根据前文分析可知, 缸内涡核位置会对DME集聚区域边缘的浓度梯度产生影响, 而随着涡核位置接近DME集聚区域位置, DME分布区域边缘流动增强, 扩散更为充分, DME分布区域边缘浓度梯度降低, 使DME自燃-火焰传播这一转换过程更加平顺, 从而稳定复合燃烧过程。
3 结论
基于光学可视化和计算流体力学仿真平台, 研究了缸内湍流运动对DME微火源引燃高活性燃料分层及燃烧过程的影响, 主要结论如下:
(1) 进气结束时刻缸内初始流场为滚流, 滚流比为1.1, 同时燃烧室和进气道结构的改变在缸内产生了4mm的涡流涡核。随着活塞移动, 滚流在燃烧室形状及涡流涡核的引导下逐渐转化为涡流, 涡流比为0.8, 此时缸内主要循环变动形式为涡核位置的循环变动。
(2) 在缸内涡流作用下, DME喷雾随流动向逆时针方向扩散, 在进排气轴线靠近进气侧分布, 该区域DME浓度与其他区域比约为5:1。此时缸内涡核位置对DME分布区域边缘的浓度梯度产生影响。
(3) 随着活塞上行, DME活性浓区由于自燃需求较低而开始自燃, 其放热率在不到5℃A的时间内超过30J/℃A, 该自燃过程放热占总热值的37.5%, DME几乎被完全消耗。活性燃料DME的自燃及火焰传播过程使缸内温度提高了400K, 引发了缸内未燃混合气自燃。
-
表 1 光学发动机主要参数
Table 1 Specifications of the optical engine
发动机形式 四冲程单缸汽油机 缸径 95mm 冲程 95mm 气门升程及相位 可变 排量 0.67L 压缩比 9.24 燃烧室结构 棚顶室 燃油喷射方式 气道喷射+缸内直喷 气道喷射燃料 PRF40 气道喷油压力 300kPa 直喷燃料 DME 直喷喷油压力 4MPa 进气方式 自然吸气 发动机转速 600r/min 节气门 有节气门 表 2 发动机运行工况
Table 2 Engine operation setup
发动机转速/(r·min-1) 600 排气门开启时刻[CA]/(°) 180 排气门关闭时刻[CA]/(°) 328 排气门升程/mm 2 进气门开启时刻[CA]/(°) 396 进气门开启时刻[CA]/(°) 508 进气门升程/mm 1 循环气道喷油量/(mg·cycle-1) 14.8 气道喷油时刻[CA]/(°) 330 循环DME直喷油量/(mg·cycle-1) 1.8 DME直喷时刻(Start of Injection, SOI)[CA BTDC]/(°) 25 过量空气系数(lambda) 1 外部EGR率 27% 进气温度/℃ 40±1 冷却水温/℃ 85±2 表 3 缸内湍动能及其循环变动
Table 3 Turbulent kinetic energy and COV in cylinder
湍动能/(m2·s-2) 循环变动(COV)/% 总湍动能 2.12275 23.20492 高频湍动能 0.21347 26.52597 低频湍动能 1.86785 25.39271 -
[1] Galloni E. Analyses about parameters that affect cyclic variation in a spark ignition engine[J]. Applied Thermal Engineering, 2009, 29(5-6):1131-1137. DOI: 10.1016/j.applthermaleng.2008.06.001
[2] Xu K, Xie H, Chen T, et al. Effect of flame propagation on the auto-ignition timing in SI-CAI hybrid combustion (SCHC)[R]. SAE Technical Papers, 2014-01-2672, 2014.
[3] Joelsson T, Yu R S, Bai X S. Large eddy simulation of turbulent combustion in a spark-assisted homogenous charge compression ignition engine[J]. Combustion Science and Technology, 2012, 184(7-8):1051-1065. DOI: 10.1080/00102202.2012.663997
[4] Yoo C S, Luo Z Y, Lu T F, et al. A DNS study of ignition characteristics of a lean iso-octane/air mixture under HCCI and SACI conditions[J]. Proceedings of the Combustion Institute, 2013, 34(2):2985-2993. http://cn.bing.com/academic/profile?id=3d13346ba0c59b3eafe6d8ff093d61a2&encoded=0&v=paper_preview&mkt=zh-cn
[5] Xie H, Xu K, Chen T, et al. Investigations into the influence of Dimethyl Ether (DME) micro flame ignition on the combustion and cyclic variation characteristics of flame propagation-auto-ignition hybrid combustion in an optical engine[J]. Combustion Science and Technology, 2016, 189(3):453-477. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=10.1080/00102202.2016.1223060
[6] Cha J, Kwon S, Kwon S, et al. Combustion and emission characteristics of a gasoline-dimethyl ether dual-fuel engine[J]. Proceedings of the Institution of Mechanical Engineers, Part D:Journal of Automobile Engineering, 2012, 226(12):1667-1677. DOI: 10.1177/0954407012450122
[7] Zhang H F. Experimental investigation of gasoline-dimethyl ether dual fuel CAI combustion with internal EGR[D]. London: Brunel University, 2011.
[8] 徐康. DME微火源对汽油机高稀释混合燃烧调控机理研究[D].天津: 天津大学, 2017. Xu K. Study of DME micro flame ignition hybrid combustion under high dilution conditions in a gasoline engine[D]. Tianjin: Tianjin University, 2017.
[9] Xie H, Lu J, Chen T, et al. Chemical effects of the incomplete-oxidation products in residual gas on the gasoline HCCI auto-ignition[J]. Combustion Science and Technology, 2014, 186(3):273-296. DOI: 10.1080/00102202.2013.858714
[10] O'Rourke P J, Amsden A A. A spray/wall interaction submodel for the KIVA-3 wall film model[R]. SAE Technical Paper, 2000-01-0271, 2000.
[11] Schmidt D P, Rutland C J. A new droplet collision algorithm[J]. Journal of Computational Physics, 2000, 164(1):62-80. DOI: 10.1006-jcph.2000.6568/
[12] Richards K J, Senecal P K, et al. Converge (Version 2.3) Manual[M]. Middleton, WI:Convergent Science Inc, 2016.
[13] Turns S R. An introduction to combustion:concepts and applications[M]. New York:McGraw-Hill Education, 2012.
[14] Cai L M, Pitsch H. Optimized chemical mechanism for combustion of gasoline surrogate fuels[J]. Combustion and Flame, 2015, 162(5):1623-1637. DOI: 10.1016/j.combustflame.2014.11.018
[15] Burke U, Somers K P, O'Toole P, et al. An ignition delay and kinetic modeling study of methane, dimethyl ether, and their mixtures at high pressures[J]. Combustion and Flame, 2015, 162(2):315-330. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=45f4569671a0afe3e958124c08d9b46a
[16] Chen Y L, Chen J Y. Application of Jacobian defined direct interaction coefficient in DRGEP-based chemical mechanism reduction methods using different graph search algorithms[J]. Combustion and Flame, 2016, 174:77-84. DOI: 10.1016/j.combustflame.2016.09.006
[17] Chen Y L, Chen J Y. Towards improved automatic chemical kinetic model reduction regarding ignition delays and flame speeds[J]. Combustion and Flame, 2018, 190:293-301. DOI: 10.1016/j.combustflame.2017.11.024
[18] Xie H, Li L, Chen T, et al. Study on spark assisted compression ignition (SACI) combustion with positive valve overlap at medium-high load[J]. Applied Energy, 2013, 101:622-633. DOI: 10.1016/j.apenergy.2012.07.015
[19] Li Y, Zhao H, Ladommatos N. Analysis of large-scale flow characteristics in a four-valve spark ignition engine[J]. Journal of Mechanical Engineering Science, 2002, 216(9):923-938. DOI: 10.1177/095440620221600906
[20] Liu D M, Wang T Y, Jia M, et al. Cycle-to-cycle variation analysis of in-cylinder flow in a gasoline engine with variable valve lift[J]. Experiments in Fluids, 2012, 53(3):585-602. DOI: 10.1007/s00348-012-1314-4
[21] Huang Y, Sung C J, Eng J A. Laminar flame speeds of primary reference fuels and reformer gas mixtures[J]. Combustion and Flame, 2004, 139(3):239-251. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=87b3f47808eb35330dfb30f8dd6d2633
[22] Jerzembeck S, Peters N, Pepiot-Desjardins P, et al. Laminar burning velocities at high pressure for primary reference fuels and gasoline:experimental and numerical investigation[J]. Combustion and Flame, 2009, 156(2):292-301.
[23] Daly C A, Simmie J M, Würmel J, et al. Burning velocities of dimethyl ether and air[J]. Combustion and Flame, 2001, 125(4):1329-1340. DOI: 10.1016-j.proci.2004.08.251/
[24] Zhao Z, Kazakov A, Dryer F L. Measurements of dimethyl ether/air mixture burning velocities by using particle image velocimetry[J]. Combustion and Flame, 2004, 139(1-2):52-60. DOI: 10.1016/j.combustflame.2004.06.009
[25] Qin X, Ju Y G. Measurements of burning velocities of dimethyl ether and air premixed flames at elevated pressures[J]. Proceedings of the Combustion Institute, 2005, 30(1):233-240. DOI: 10.1016-j.proci.2004.08.251/
[26] Wang Y L, Holley A T, Ji C, et al. Propagation and extinction of premixed dimethyl-ether/air flames[J]. Proceedings of the Combustion Institute, 2009, 32(1):1035-1042.