试验样机的主要技术参数如表1所示, 试验用天然气燃料的主要成分如下所示:CH₄、C₂H₆、C₃H₈、C₄H₁₀、N₂的体积分数分别为85.45%、4.51%、3.39%、3.71%和2.94%。图1为光学测试系统的示意图, 由进气系统、CNG喷射系统、变频倒拖系统、点火控制系统、光学采集系统以及光学单缸试验机等组成。图2为样机燃烧室的横剖图, 在缸盖上布置两个喷射器和两个火花塞, 喷射压力约为5 MPa, 可实现不同喷射方式和点火方式。通过CCD(Charge coupled device)高速摄像机以563幅/s的速度连续拍摄燃烧室内的火焰传播过程, 观测不同喷射条件下的火焰传播特性。燃烧室的可视率为气缸直径的67%, 试验的同时绘制示功图, 以分析燃烧过程。试验数据是每种工况测量10次并进行平均化处理。

表1

Table 1

表1(Table 1)

表1 试验样机技术参数 Table 1 Technical parameters of experimental prototype 型式 单缸4冲程点燃式 缸径/mm 135 行程/mm 280 排量/L 4.0 压缩比 6.13

表1 试验样机技术参数 Table 1 Technical parameters of experimental prototype

图1

Fig.1

	Figure Option ViewDownloadNew Window
	图1 试验设备示意图Fig.1 Schematic of experimental facilities

图2

Fig.2

	Figure Option ViewDownloadNew Window
	图2 火花塞、喷射器安装位置示意图Fig.2 Schematic of position of spark plug and injector

表2为采用双点点火、点火时刻为(θ _i1, θ _i2)=(356, 357) ° CA、当量比 ϕ为0.7时不同喷射时刻所对应的喷射终了时刻。此时喷射器的喷射持续期间为78 ° CA曲轴转角, 从喷射终了时刻到点火时刻所经历的时间为Δ t, Δ t的大小可直接表示燃料在缸内直喷后混合气形成时间的长短。

表2

Table 2

表2(Table 2)

表2 不同试验条件下的喷射定时 Table 2 Injection timing in different experimental conditions n /(r· min-1) ϕ 喷射开始 时刻/° CA 喷射终了 时刻/° CA 喷射持续 期间/° CA Δ t /ms 200 0.7 180 -102 78 81.7 200 0.7 240 -42 78 31.7 200 0.7 90 -12 78 6.7

表2 不同试验条件下的喷射定时 Table 2 Injection timing in different experimental conditions

在FIRE软件里利用拓扑Topology中的Cylinder生成气缸模型, 并通过Topology自动生成计算网格, 为了减少计算量并保证不同层动网格特性一致, 根据活塞的运动特性把进气下止点到做功行程下止点期间, 划分成如图3所示的4层不同的网格:从进气下止点的180~240 ° CA和对应的膨胀过程后期的480~540 ° CA期间采用60层网格; 压缩过程中期的240 ~270 ° CA和对应的膨胀过程中期的450~480 ° CA期间采用50层网格; 压缩过程后期的270 ~350 ° CA和对应的膨胀过程初期的370 ~450 ° CA期间采用45层网格; 压缩上止点(360 ° CA)附近350 ~370 ° CA期间采用10层网格。由此在保证计算精度的前提下尽可能减小计算工作量。

图3

Fig.3

	Figure Option ViewDownloadNew Window
	图3 计算网格Fig.3 Computational mesh

仿真计算时间步长在喷油和燃烧阶段采用0.5 ° CA曲轴转角, 而膨胀阶段采用1 ° CA曲轴转角; 根据CNG气体燃料的特点, 仿真子模型湍流模型采用k-ζ -f模型; 混合气的燃烧模型采用拟序火焰拓展模型ECFM; 湍流扩散模型采用O'Rourke模型。

作为计算边界条件, 采用样机的实测值并令气缸入口处气流状态为均匀, 其平均压力为9× 10⁴ Pa, 温度为350 K; 活塞顶面作为移动边界, 其温度为593 K, 气缸壁和气缸盖底面为固定边界, 令气缸壁的温度为403 K, 气缸盖底面的温度为593 K。

图4和表3分别为当量比ϕ =0.93、双点点火时刻为(θ _i1, θ _i2)=(356, 357) ° CA、两个喷射器同时喷射、喷射时刻为θ _inj=240 ° CA时的放热规律和火焰传播速度的实测值与仿真计算值的对比。仿真计算的放热规律曲线特性和火焰传播特性与试验值吻合良好, 说明所创建的数值仿真模型符合CNG燃料缸内直喷时实际混合气形成及火焰传播规律。

图4

Fig.4

	Figure Option ViewDownloadNew Window
	图4 放热规律对比Fig.4 Comparation of heat release rate

表3

Table 3

表3(Table 3)

表3 火焰传播过程的仿真和试验结果对比 Table 3 Comparation of simulation and experimental result during flame propagation process 曲轴转角/ ° CA 火焰面积/mm2 试验值 仿真值 357 50.39 54.14 359 180.97 193.29 361 580.12 593.29 363 2457.82 2480.34 366 4872.36 4886.23

表3 火焰传播过程的仿真和试验结果对比 Table 3 Comparation of simulation and experimental result during flame propagation process

图5为不同喷射时刻对缸内平均湍流动能动态变化特性的影响。当早期喷射时, 因缸内背压低, 因此湍流动能随喷射过程快速增加并达到最高峰值, 喷射后期湍流动能迅速衰减; 随喷射时刻的退后, 缸内压力增大, 随喷射过程的进行湍流动能上升速率和峰值逐渐减小; 在上止点附近点火后, 由于火焰传播, 使湍流动能有所增大, 但火焰传播引起的湍流动能变化量与喷射过程对湍流动能的影响相比可忽略。

图5

Fig.5

	Figure Option ViewDownloadNew Window
	图5 不同喷射时刻缸内湍流动能对比Fig.5 Comparation of TKE in cylinder in different injection timing

图6为在压缩上止点附近缸内流场随不同喷射时刻的变化特性。在横截面A上由于进气涡流的导向作用不同喷射时刻缸内气流整体绕气缸中心线顺时针方向旋转, 形成外强内弱的气流分布特性, 且随喷射时刻的不同, 各瞬态气流强度分布特性不同, 峰值出现点也不一样, 根据火花塞点火位置和火焰传播方向, 这种气流的分布特性直接影响火焰的传播速度; 而在纵截面B上, 进一步证实气流强度除了具有外强内弱的分布特性外, 随着喷射时刻的滞后, 燃烧室顶部的气流强度增强。特别是当喷射时刻为270 ° CA时, 由于喷射期间活塞已快速向上止点移动, 因此喷射的气柱在活塞顶部产生碰壁现象, 导致纵向气流比较活跃, 在点火之前缸内形成上、下及边缘气流强度大而中心气流强度弱的分布特性。这种不同的速度场动态分布特性决定缸内混合气浓度场的分布特性, 因此点火后造成火焰传播速率不同, 导致膨胀过程中形成如图5所示的湍流特性。

图6

Fig.6

	Figure Option ViewDownloadNew Window
	图6 不同喷射时刻对速度场的影响Fig.6 Effect of different injection timing on velocity field

天然气的主要成分为甲烷(CH₄), 其相对空气的密度约为0.6, 因此其混合气形成机理和混合气浓度分布特性取决于天然气和空气的相对运动特性。图7为不同喷射时刻对不同密度的两种气体的相对运动规律以及由此决定的混合气浓度分布特性的影响。由此可知, 在压缩行程为250 ° CA、喷射时刻为180 ° CA BTDC时, 喷射过程快要结束, 喷气动能减弱, 导致喷气的贯穿能力也弱, 同时天然气的密度远比空气小, 所以随活塞的上移大多数燃气与空气分离并浮上燃烧室顶部; 当喷射时刻推迟到240 ° CA BTDC时, 借助喷气的动能使喷气贯穿到空气层之中, 减缓了压缩过程中天然气与空气的分离现象, 使缸内纵向混合气流分布趋于均匀; 当喷射时刻推迟到270 ° CA BTDC时, 在点火之前喷气过程尚未结束, 使喷注贯穿到活塞顶, 从而在缸内纵向气流运动区域拓宽, 湍流强度增大, 同时有效地控制天然气与空气的分离速度。

图7

Fig.7

	Figure Option ViewDownloadNew Window
	图7 不同喷射时刻对缸内流场及浓度分布规律的影响Fig.7 Effect of different injection timing on flow field and concentration distribution law

图8为喷射时刻对上止点附近混合气浓度场分布特性的影响。

图8

Fig.8

	Figure Option ViewDownloadNew Window
	图8 不同喷射时刻对上止点附近浓度场演变规律的影响Fig.8 Effect of different injection timing on variation law of concentration field near TDC

在火花塞布置的横截面A上(见图1), 喷射时刻越早, 在点火时刻附近混合气浓度横向分布越均匀, 推迟喷射时刻, 火花塞附近混合气浓度梯度增大; 而在纵截面B上混合气浓度分布特性则相反, 喷射时刻过早, 由于天然气密度比空气小, 所以混合气浓度呈上浓下稀的云彩状分布, 在火焰传播期间燃料与空气分层明显; 当推迟喷射时刻到240 ° CA BTDC时, 喷射时期喷气动能与缸内空气流匹配良好, 使缸内形成较均匀的混合气; 当喷射时刻过晚(为270 ° CA BTDC)时, 因喷注与活塞顶相碰, 造成混合气上稀下浓的分布特性。

当平均当量比ϕ =0.7、两个喷油器同时喷射、双点点火时刻为(356, 357) ° CA时, 喷射时刻对火焰传播特性的影响如图9所示。当天然气向缸内直喷时, 火焰传播速率取决于混合气的浓度分布特性及天然气向空气层的扩散速度:①当喷射时刻为180 ° CA BTDC时, 根据此时的混合气浓度分布特性(见图8), 在横截面A上因混合气形成时间较长, 混合气浓度分布相对比较均匀, 而在纵截面B上因天然气密度比空气小, 难以向底层空气层扩散, 因而形成上浓下稀的梯度分布(天然气与空气分层), 所以火焰传播速度缓慢, 使得放热速率降低, 放热率峰值减小, 燃烧期间相对迟后(见图10); ②当喷射时刻为240 ° CA BTDC时, 在横截面上形成有一定梯度的混合气浓度分布的同时, 在纵截面上形成较均匀的混合气, 因此火焰传播速度得到改善, 提高了燃烧放热速率, 放热率峰值得到提高, 放热率曲线中心有所提前; ③当喷射时刻为270 ° CA BTDC时, 燃烧室横截面上形成的混合气浓度梯度最大, 而在纵截面上形成下浓上稀的混合气, 密度小的天然气易向上扩散, 有利于提高火焰传播速度, 因此燃烧放热速率及放热率峰值明显提高, 燃烧持续期明显缩短(见图10), 同时用平均指示压力的平方差表示的循环变动也明显减小(见图11), 表明稀薄燃烧过程更加趋于稳定。

图9

Fig.9

	Figure Option ViewDownloadNew Window
	图9 不同喷气正时的火焰传播情况Fig.9 Flame propagation condition in different injection timing

图10

Fig.10

	Figure Option ViewDownloadNew Window
	图10 喷射时刻对燃烧压力及放热规律的影响Fig.10 Effect of injection timing on combustion pressure and heat release rate

图11

Fig.11

	Figure Option ViewDownloadNew Window
	图11 喷射方式对平均指示压力及其偏差的影响Fig.11 Effect of injection timing on IMEP and its deviation

图12

Fig.12

	Figure Option ViewDownloadNew Window
	图12 不同喷射时刻的温度场Fig.12 Temperature field in different injection timing

NO的生成条件是高温、富氧, 因此在燃烧过程中缸内混合气浓度场和温度场及其动态变化特性直接影响NO的生成规律。图12为喷射时刻对NO开始生成期间缸内温度场分布特性的影响。由于燃烧后的高温区生成于火焰传播后的已燃区, 因此NO易生成于火焰带之后的高温区, 而且其生成速率取决于当时当地的温度和浓度条件。

当喷射时刻为180 ° CA BTDC时, 虽然混合气浓度稀薄区较宽, 但由于燃烧放热速率缓慢, 燃烧初期缸内温度普遍较低; 当喷射时刻为240 ° CA BTDC时, 缸内混合气浓度趋于均匀, 燃烧放热速率提高, 高温区迅速增大, 不仅NO生成速率增加, 而且NO快速生成的反应区域面积增大, 持续时间延长(见图13), 因此NO最终生成量最多(见图14); 当喷射时刻为270 ° CA BTDC时, 如前所述燃烧放热速率最大, 因此缸内温度迅速升高, 同时高温区也扩大, 但由于混合气形成下浓上稀的分布特性, 因此虽然某一时刻(图13中372 ° CA时), 局部地区的NO生成速率最大, 但NO生成的反应区域面积较小(带状), 持续时间也较短, 所以最终NO的生成量较少。

图13

Fig.13

	Figure Option ViewDownloadNew Window
	图13 不同喷射时刻的NO生成速率Fig.13 NO formation rate in different injection timing

图14

Fig.14

	Figure Option ViewDownloadNew Window
	图14 不同喷射时刻NO的质量分数Fig.14 NO mass fraction of different injection timing

图15表示在喷射时刻为240 ° CA BTDC、曲轴转角为376 ° CA时, 缸内混合气浓度场、温度场和NO反应速率之间的对应关系。图中用A、B、C₁~C₃分别表示不同NO反应速率及其对应的混合气浓度范围和温度范围, 其中C₁~C₃表示NO反应速率为较小的区域。

图15

Fig.15

	Figure Option ViewDownloadNew Window
	图15 NO的生成条件Fig.15 NO formation condition

表4为不同NO反应速率及其对应的混合气浓度和温度情况。结果表明:NO的生成速率是当量比和温度的共同函数, 并非温度(或氧浓度)越高NO反应速率就越快, 而是在一定的氧浓度和温度下(如当量比为0.81~0.86, 温度为2468~2536 K时)NO生成速率最高, 当 ϕ为0.86~0.9时, 即使温度提高到2655~2700 K, NO的生成速率明显降低; 当 ϕ为0.65~0.75时, 温度超过2000 K时, 随温度的升高NO生成速率有所增大; 但当温度低于2000 K时不管是否为富氧条件, NO生成速率几乎为零。

表4

Table 4

表4(Table 4)

表4 A、B、C1~C3处的当量比和温度范围 Table 4 Equivalence ratio and temperature scope in A、B、C1~C3 NO生成速率 /(10-3 s-1) 当量比 温度/K A 2.15~2.64 0.80~0.86 2468~2536 B 0.49~0.63 0.86~0.9 2655~2700 C1 0.1~0.2 0.66~0.69 2221~2243 C2 0.02~0.07 0.65~0.74 1929~2220 C3 0~0.01 0.53~0.67 1559~2000

表4 A、B、C1~C3处的当量比和温度范围 Table 4 Equivalence ratio and temperature scope in A、B、C1~C3