0 引言

航空煤油和高压氧气按化学反应式的恰好完全反应的摩尔比在燃烧室内混合,航空煤油被氧气雾化后形成混合液雾,将液雾点火并通过拉瓦尔喷管加速后喷出,可以形成高温、高速的射流。改变航空煤油和高压氧气输入的混合比,即改变混合液雾的富氧度,可改变射流的速度和温度的分布。该射流具有热值高、比冲大等特点,在喷涂、烧蚀、航天发动机等工业和航空航天等领域中具有广阔的应用前景[1]。

国内外在煤油和氧气混合火焰的研究等方面做了大量的工作:Read等人对低温和低压下航空发动机的点火过程和内流场的结构进行了研究[2]。Vasu等人研究了激波管中航空煤油燃烧的延迟时间,并对混合气的当量比的影响进行了归一化分析[3]。国内方面,雷向东等研究了高空氧气/煤油点火装置[4]。梁金虎等人对RP-3航空煤油的点火特性开展了研究[5]。然而,由于发生器燃烧室内的温度高、出口射流的速度快,故对射流的温度和速度等参数都很难直接测量,而超音速射流的速度场和温度场分布的不确定是制约其更广泛应用的关键问题。

文中研究通过理论仿真计算和实验验证相结合的方法[6-7],对航空煤油/氧气超音速射流场的特征参数开展研究。通过理论计算获得射流轴心线上速度和温度的分布规律,数值仿真计算获得整个射流场的温度和速度等参数的分布。实验中,首先采用金属熔点法,定点测量射流场的温度分布;其次为了能够实时测量射流场参数的变化情况,在出口射流中加入三氧化二铝粉末,采用Spray watch粒子测试仪对射流中粉末的温度和速度等参数进行实时测量;再次改变煤油和氧气的混合比,测量在不同的富氧度下,距离发生器出口不同距离处射流的温度和速度的分布。比较仿真计算结果和测试结果,可以获得航空煤油/氧气超声速射流场的参数分布。

1 航空煤油/氧气超声速射流系统

通过高压油泵将作为燃料的航空煤油输送至系统的燃烧室内,同时通过减压阀将作为助燃剂的氧气输送到燃烧室中,输入的航空煤油在氧气的吹射下雾化成可燃液雾,液雾经火花塞点火后燃烧,形成高温高压的燃气,燃气通过拉瓦尔喷管加速至超音速从出口喷出[8-9],实验系统如图1所示。

图1 航空煤油/氧气超音速射流实验系统组成

为了实时地测量射流场的温度和速度的变化情况,将三氧化二铝粉末从拉瓦尔喷管喉部的负压区送入超音速射流中,粉末在加长喷管中经射流加温加速后喷出,形成带颗粒的超音速气固两相流,采用Spray watch粒子测试仪可以实时地测量到高温射流中粉末颗粒的温度和速度等参数。气固两相射流的图片如图2所示。

图2 航空煤油/氧气超音速火焰射流图

2 超音速射流场的计算

2.1 理论计算

式中:Te为喷管出口处射流的温度;Ta为周围大气的温度。

由于钼丝的熔点相对最高,测试得到钼丝的熔点的距离也最近,而铁丝的熔点相对最低,其测得的熔点的距离也最远,因而可以证明随着距离的增加,射流的温度也逐渐递减。另外根据测试结果可以发现:随着富氧度增加,金属丝在射流中熔点的距离也越来越近,说明随着富氧度的增加射流的温度也逐渐降低。

何谓“新工科”，通过对《教育部高等教育司关于开展新工科研究与实践的通知》（教高司函〔2017〕6号）文件的学习，主要内容被归纳为“五个新”，即工程教育的新理念、学科专业的新结构、人才培养的新模式、教育教学的新质量、分类发展的新体系[3]。

(1)

根据仿真计算的结果可以得出:当射流在流过拉瓦尔喷管的喉部区域时,速度达到了当地的音速;当射流继续流过喷管的扩张段时,由于此时过流断面不断增大,使射流的速度在此处得以继续增加,达到了超音速;当射流在加长的喷管中流动时,由于射流的外沿与管壁之间存在速度差和能量交换的相互作用,使射流外沿的速度不断减小,而温度则有所上升;为了防止系统在工作时的高温烧穿燃烧室和喷管,需要对燃烧室和喷管采用冷却水进行强制冷却,因而射流在燃烧室和喷管中流动时温度将不断地下降,这种冷却过程持续到射流从喷管的出口喷出为止。通过仿真计算的结果可以得出:射流在喷管的出口处的速度约为V=2 150 m/s,温度约为T=2 780 K。当射流进入大气后,与周围的空气发生剧烈的膨胀、压缩和阻滞作用,射流因而发生连续的扩张和收缩等现象,沿轴线形成一系列的马赫锥;随着距离出口距离的增加,射流的速度不断地降低。射流在流动的过程中,与周围的大气发生强烈的对流和热辐射等热交换现象,因而射流的温度也急剧下降。

(2)

由于射流在空气中要与周围的环境大气产生传热、对流、热辐射等气动热效应,射流外围的参数变化剧烈,通常以射流轴心线上的气动参数来表征射流,射流轴心线上的速度可由式(1)表达[10]。

要做好金融扶贫工作，对相关风险的监测与管控力度要加强。一是要加强相关客户风险分析，要紧密关注农产品价格变化，做好风险管控预案。二是要防止农产品过剩风险。由于农产品不便保存的特点，要高度关注其市场趋势，加强监测，对发现的问题及时汇报、分析并解。

(3)

式中:ΔTg=Tg-Tα;ΔTe=Te-Tα;Tg为射流的温度。

实验测得超音速射流的出口速度和温度值分别为2 150 m/s和2 800 K,射流的扩展角约为6°。假定实验时周围环境的温度为恒温,取为303 K(30 ℃),通过计算可得焰流轴心的速度和温度的分布如图3所示。

采用TRIAX550光谱仪对射流的成分进行测量,其成分主要为CO、CO2和H2O,故在燃烧室中简化的化学反应方程式如下:

图3 理想状态下射流轴心线速度和温度分布

2.2 FLUENT仿真计算

高洪臣，出生年月：1980年2月25日，性别：男，籍贯：黑龙江大庆，民族：汉，最高学历：本科，职称：中级。

2.4 患者注射胰岛素行为得分情况 300例患者中，行为得分最高96.0分，最低52.0分，平均得分(69.64±10.02)分，处于中等水平。

C12H24+6O212CO+12H2

(4)

2CO+O22CO2

(5)

2H2+O22H2O

(6)

采用压力传感器测定燃烧室工作时的压力约为1.5 MPa。应用Fluent数值仿真计算软件对整个射流场进行仿真计算,采用有限体积法(FVM)、稳态、隐式解法步进求解[11-13],根据N-S方程,采用标准的k-ε湍流模型,结合燃烧室和喷管的边界条件,计算出航空煤油和高压氧气在充分燃烧时射流场参数的分布,计算结果如图4所示。

图4 仿真计算射流场温度和速度的分布图

式中:Vg为射流的速度;Ve为喷管出口处射流的速度;re为喷管出口的半径;x为离喷管出口的轴向距离;α为射流的扩展角,即射流外边界夹角称射流的扩展角;θ为无因次温度系数,为喷管出口处的射流温度与周围大气的温度之比,即:

射流轴心线处的温度可由式(3)表达[10]:

3 实验测试

3.1 熔点法测试

采用熔点不同的金属丝对射流的温度进行实验标定。分别将直径为1 mm的钼丝(熔点为2 617 ℃)、锆丝(熔点为1 852 ℃)和铁丝(熔点为1 534 ℃)置入射流的中心线上,并沿中心线由射流的低温区向高温区方向缓慢移动,直至金属丝的尖端熔化为止,测量此时不同的金属丝距射流出口端面的距离。射流在此处的温度等于对应着金属丝的熔点;通过流量控制阀调节输入的航空煤油和氧气的流量可以改变燃烧室内航空煤油和氧气的混合比,得到富氧度约为0%～15%的混合比,采用熔点法重复测试不同富氧度时射流的温度,由此标出了不同富氧度时射流的温度变化如图5所示。

图5 富氧度与熔点距离的关系示意图

系统功能语言学研究认为口头语言具有“信息索取”与“提供产品或服务”的功能。近年来，一些国外学者认为图片本质上与口头语言类似，也带有语言的这两种功能，画中人面部和身体姿势可用来实现潜势意义（Kress&Van Leeuwen 2006；Painter，Martin&Unsworth 2013）。画中人注视读者意味着画中人向读者“索取”（demand）；画中人没有注视读者，意为“提供信息”（offer of information）。面部和身体姿势可以展示人物的真实情感。该公益广告视频中的人物的目光多数没有注视读者，表明该广告的目的是提供信息，而非索取信息。

3.2 粒子测试仪测试

为了获得射流的瞬时速度[14-15],通过粒子输送装置在拉法尔喷管喉部的负压区输入粒径约为50～70 μm的三氧化二铝粉末,此时氧气的流量约为4.5 m3/h,实验中调节航空煤油的流量控制阀,使航空煤油的流量由5 L/h降至4 L/h,从而获得富氧度为0%～15%的不同的实验工况,此时燃烧室的压力稳定在1.6 MPa左右。采用Spray watch粒子在线测试仪测量射流中粒子的温度和速度,获得实验参数如图6所示。

图6 不同富氧度时粒子的温度和速度测试结果

由图6(a)可以发现:随着富氧度的增加,射流中温度相同的点越来越靠近出口,即射流场的温度随着富氧度的增加而降低,原因是过量的氧气没有参与到化学反应中去,因而随着射流喷射出来时,和射流中的高温燃气不断的进行对流换热,从而使射流的温度下降。

由图6(b)可以发现,在测试范围内,随着距离的增加,射流中粒子的速度也增加。原因是粒子的速度是射流与粒子之间拖动力互相作用的结果,而拖动力的大小取决于射流与粒子之间的相对速度,其大小可表示为：

(7)

式中:Ad为粒子的表面积;Cdrg为拖动系数;ρg为气流的密度;Vg为气流的速度;Vd为粒子的速度。由式(7)可以看出,射流与粒子的相对速度越高,其拖动力也就越大,射流对粒子的加速作用也就越强。在拉瓦尔喷管的喉部,射流的速度约为音速,而粒子的速度约为零,故此时的拖动力最大;当射流在加长喷管中流动时,随着距离的增加,粒子在射流的拖动力作用下,速度逐渐增加;当射流由喷管喷出进入大气后,射流在周围大气的作用下,由层流变为湍流,并形成一系列的马赫波,射流的速度逐渐降低,但此时射流的速度依然大于粒子的速度,射流依然对粒子施加拖动力;当射流的速度最后和粒子的速度趋于一致时,射流不再对粒子施加拖动力,此时粒子的速度达到最高。另外,由于随着富氧度的增加,燃气的初速度降低,因而射流中粒子的速度也随着富氧度的增加呈下降趋势。

4 结论

通过研究,得出了航空煤油/氧气超音速射流的出口速度和温度约为2 150 m/s和2 800 K;随着距出口距离的增加,射流的温度和速度都呈抛物线型下降;随着富氧度的增加,射流中同一点的温度和速度都显著降低。

子贡问曰：“孔文子何以谓之‘文’也？”子曰：“敏而好学，不耻下问，是以谓之‘文’也。”(子贡问道：“孔文子凭什么能得到‘文’的谥号呢？”孔子说：“他这个人头脑敏捷而爱好学习，并且不以向学问和地位不如自己的人求教为耻，所以被封为‘文’。”)

参考文献:

[1] 李斌, 谭永华. 液氧/煤油富氧补燃发动机 [J]. 火箭推进, 2003, 29(2): 1-6.

[2] READ R W, ROGERSON J W，HOCGREB S. Relight imaging at low temperature, low pressure conditions：AIAA 2008-0957 [R].[S.l.]:AIAA, 2008.

[3] Vasu S S, Davidon D F. Jet fuel ignition delay times: shock tube experiments over wide conditions and surrogate model predictions [J]. Combustion and Flame, 2008, 152(1/2): 125-143.

[4] 雷向东, 余协正, 杨京军, 等. 氧气/煤油点火装置高空点火试验研究 [J]. 火箭推进, 2014, 40(6): 37-43.

[5] 梁金虎, 王苏, 张灿, 等. RP-3航空煤油点火特性研究 [J]. 力学学报, 2014, 46(3): 352-360.

[6] 吴宝元, 张民庆. 液氧/煤油富氧预燃室的理论分析和试验研究 [J]. 推进技术, 1999, 20(5): 11-16.

[7] 雷震, 代玉东. 液氧/煤油发动机试车热敏电阻测温技术 [J]. 火箭推进, 2003, 29(4): 35-39.

[8] 王爱玲, 吴建军. 气氧/煤油发动机电点火技术研究 [J]. 上海航天, 2006(6): 14-17.

[9] 赵正社. 液氧/煤油发动机试验低温测量系统的设计与实现 [J]. 火箭推进, 2001, 27(2): 36-41.

[10] 张福祥. 火箭燃气射流动力学 [M]. 哈尔滨: 哈尔滨工程大学出版社, 2004.

[11] 曹再勇, 蔡体敏, 谭永华. 富氧燃气发生器三维燃烧流场的数值模拟 [J]. 西北工业大学学报, 2005, 23(5): 557-561.

[12] 翟鹏，杨月诚.基于FLUENT的烧蚀喷枪两相流流场数值模拟 [J].弹箭与制导学报， 2012，32(3)：133-136.

[13] 程钰锋, 聂万胜, 丰松江. 氢氧火箭发动机不稳定燃烧数值研究 [J]. 装备指挥技术学院学报, 2009, 20(4): 69-73.

[14] 饶琼，周香林，张济山，等.超音速喷涂技术及其应用[J].热加工工艺， 2004(10)：49-52.

[15] 于永志. 固体火箭发动机尾焰温度场特性研究 [D]. 沈阳: 沈阳理工大学, 2016.