相恒升, 陈 雄, 周长省, 赖华锦
( 南京理工大学机械工程学院, 江苏 南京 210094 )
固体燃料冲压发动机以贫氧推进剂为燃料,工作时利用大气中的氧气助燃[1],由于随着海拔高度增加大气中氧含量减小,从大气中获取氧气使固体燃料冲压发动机性能受大气条件限制。贫氧推进剂自身不含氧化剂,或者氧化剂含量大大低于常规固体火箭推进剂中的氧化剂含量[2],如铝镁贫氧推进剂。大气中氧含量变化势必影响贫氧推进剂燃烧性能,因此研究环境氧含量变化对贫氧推进剂燃烧性能的影响有重要意义。
铝镁贫氧推进剂可以大幅度提高冲压发动机的性能,已被广泛应用于固体燃料冲压发动机[3]。陈超等[4]研究了铝粉粒径对高铝含量贫氧推进剂一次燃烧性能的影响,揭示了铝粉在推进剂中的燃烧机理。郑磊等[5]研究了镁粉粒度对镁/聚四氟乙烯(Mg/PTFE)贫氧推进剂燃烧性能的影响,发现铝粉粒度影响推进剂燃速和燃烧火焰温度。聂芝侠等[6]实验研究了氧化剂级配、催化剂及金属总含量等对推进剂燃速与压力指数的影响。冯喜平等[7]采用数值模拟的方法分析了旋流对贫氧推进剂燃烧效率的影响。然而关于环境氧含量与压力对铝镁贫氧推进剂燃烧性能的影响还少有报道。
本研究采用激光辐射点火,使用高速摄影仪记录推进剂点火与燃烧过程,利用红外测温仪测量推进剂表面温度以及火焰温度,讨论了不同压力和环境氧含量下推进剂点火过程,利用B数理论分析了环境氧含量和压力对燃速的影响,研究结果可为大海拔跨度使用的固体燃料冲压发动机设计提供参考。
铝镁贫氧推进剂由西安近代化学研究所提供,其主要配方(质量分数)为: 铝粉20%,球形铝粉FLQT2,平均粒径24 μm; 镁粉20%,平均粒径29.97 μm; 端羟基聚丁二烯(HTPB)体系20%; 高氯酸铵(AP)36%,粒径 100~120 μm; 小组分物质4%; 试样中不含催化剂。样品尺寸为Φ4 mm×7 mm,为防止样品侧面燃烧,样品圆柱面用硅橡胶包覆。
实验装置主要由 CO2激光器、燃烧室和数据采集设备组成,其组成示意图如图1所示。燃烧室两侧为钢化玻璃透明窗(燃烧室详细介绍见文献[8]),燃烧室激光透镜下有火焰阻焰装置[9],燃烧室与高压气瓶相连,通过更换不同氧含量的高压气瓶调节燃烧室内氧含量。CO2激光器,最高输出功率为200 W,激光波长10.6 μm,光斑直径3.7 mm,上海御鸿激光制造; 红外测温仪,MODLINE 5R-3015型; 高速摄影仪,三洋HD2000,拍摄速度为240 fps。
图1 实验装置示意图
Fig.1 Schematic diagram of the experimental setup
CO2激光器采用点射模式将推进剂样品点燃,激光器工作时间为(100±5)ms,由于激光工作时间远小于推进剂燃燃烧时间,所以可忽略激光对推进剂燃速的影响。高速摄影仪记录推进剂点火与燃烧过程,红外测温仪记录推进剂燃烧过程中红外测温仪镜头焦点位置处温度变化过程,点火前红外测温仪镜头焦点在推进剂表面。由于红外光线无法透过钢化玻璃,红外测温仪测温时采用比色测温模式,具体原理见文献[10]。通过更换与燃烧室相连的高压气瓶使燃烧室内氧含量依次为10%、14%、18%和21%。试验过程中恒压阀可保证在推进剂燃烧过程中不断有新鲜环境气体流入燃烧室的同时燃烧室压力不变。
铝镁贫氧推进剂被激光点燃后火焰传播至整个样品表面计为推进剂开始燃烧时刻,推进剂火焰高度降至稳定燃烧时火焰高度的一半时认为推进剂燃烧结束,从高速摄影仪拍摄到的推进剂燃烧过程中计算出推进剂燃速[5,11-13]。每种氧含量下分别调节燃烧室压力为0.1(常压),0.35,0.6,0.85,1.1,1.3,1.5 MPa。燃烧室内氧含量记为φO2,燃烧室压力记为p。每种工况进行6次试验,试验温度为室温,25 ℃。
不同压力与不同环境氧含量下铝镁贫氧推进剂点火与燃烧过程如图2所示,图中0 ms是视频中推进剂初次出现火焰的时刻,视频中初次出现的火焰简称初现焰,假设初现焰即为推进剂点火时的第一个火焰。
由图2可见,当p=0.1 MPa、φO2=10%时,初现焰紧靠推进剂表面出现,火焰从推进剂表面开始向外传播; 当p=0.1 MPa、φO2=21%时,球形初现焰出现在推进剂上方,距推进剂表面约15 mm,火焰向推进剂表面传播; 当p=0.1 MPa、φO2=18%时,初现焰同时出现在推进剂表面与推进剂上方,之后这两个位置处的火焰相互靠近,在6.25 ms时两处火焰会合。然而,随着压力的增加,如图2中p=0.35 MPa与1.1 MPa时,即使在高环境氧含量条件下点火铝镁贫氧推进剂的初现焰距离推进剂表面也较近。
利用图3所示的推进剂热解产物与环境气体的扩散模型可以解释上述现象。图3中A区是在激光辐射作用下铝镁贫氧推进剂热解产物的分散区,C区是环境气体,B区是推进剂的热解产物与环境气体的扩散区。基于图3的扩散模型分析出现上述现象的原因为: 由于铝镁贫氧推进剂含有一定量的AP,其热解产物中既含有氧化性气体也有还原性气体,根据该铝镁贫氧推进剂的配方计算得热解产物中氧含量为18.7%。当环境气体氧含量低于推进剂热解产物中氧含量时,B区中的氧化性气体被环境气体稀释,使B区中的氧含量比A区低,初现焰出现在氧含量高的A区,所以p=0.1 MPa、φO2=10%时初现焰出现在推进剂表面。随着铝镁贫氧推进剂不断热解,火焰从A区向上传播; 当环境气体氧含量高于推进剂热解产物中氧含量时,B区氧含量比A区高,所以铝镁贫氧推进剂在p=0.1 MPa、φO2=21%环境中点火时初现焰出现在远离推进剂表面的B区,火焰从B区向A区传播; 当环境气体氧含量为18%时,与热解产物氧含量18.7%相差不大,所以初现焰同时出现在A区与B区。但随着压力增加,如p=0.35 MPa与1.1 MPa时,推进剂热解产物扩散过程减弱,B区与推进剂表面的距离减小,使初现焰靠近推进剂表面。
图2中铝镁贫氧推进剂稳定燃烧时的火焰是在高速摄影仪镜头前加滤光片拍摄到的。铝镁贫氧推进剂稳定燃烧时火焰基本为柱状。相同压力下随着环境氧含量的增加,推进剂火焰尺寸逐渐变大,燃烧更剧烈。
p=0.1 MPa、φO2为10%~21%时铝镁贫氧推进剂无法稳定燃烧,点火之后慢慢熄灭。分析原因为: 当压力小时,推进剂燃烧时火焰“吹离效应”[14]明显,燃烧产生的热量反馈到推进剂表面少,推进剂表面热损失大于火焰反馈热,推进剂表面温度逐渐降低,最终导致熄火。
图2 不同压力与环境氧含量下铝镁贫氧推进剂点火与燃烧过程
Fig.2 Ignition and combustion process of aluminum-magnesium fuel-rich propellant at different pressure and different environment oxygen content
图3 推进剂热解产物与环境气体的扩散模型
Fig.3 Model of the diffusion process between propellant thermal decomposition products and environment gas
铝镁贫氧推进剂在不同压力与不同环境氧含量条件下点火过程中红外测温仪焦点处温度变化如图4所示,由图4可见,在激光辐射25 ms后推进剂表面达到1000 ℃(因红外测温仪测量范围为1000~3000 ℃,当被测物体温度小于1000 ℃时红外测温仪电信号一直输出1000 ℃的值),之后约6 ms推进剂表面温度上升到1150~1200 ℃,此时出现一个温度平台,平台温度维持约12 ms。形成这个温度平台的原因为铝镁贫氧推进剂在1150~1200 ℃时开始热解,推进剂热解吸收的能量约等于推进剂吸收的激光能量,所以此时推进剂表面温度基本不上升。这说明该铝镁贫氧推进剂的燃面温度约1150~1200 ℃,此数值与热电偶测到的温度相近(热电偶测到的燃面温度约1135~1282 ℃)。经过70 ms后火焰温度逐渐到达点火过程中的最高温度。100 ms时激光辐射消失,由于火焰温度梯度不能立即改变,原有的热平衡被打破,散热比产热多,火焰温度降低。经过约100 ms后热平衡重新建立,火焰温度开始上升,直至达到稳定。
图4 不同压力与环境氧含量下铝镁贫氧推进剂点火过程中温度变化
Fig.4 Temperature variation in the ignition process of aluminum-magnesium fuel-rich propellant at different pressure and different environment oxygen content
不同压力与不同环境氧含量下铝镁贫氧推进剂稳定燃烧时火焰温度Tf如图5所示,由图5可见,相同压力下,稳定燃烧时的火焰温度随环境氧含量增加线性增加,这个变化趋势与Quintiere等[15]的研究结论相符。
(1)
式中,T∞为环境温度,K;Xr为火焰辐射系数;Kp为系数;hc为燃烧热,J·kg-1;r为化学计量燃料氧气质量比;cp为定压比热, J·kg-1·K-1。
图5 不同压力与环境氧含量下铝镁贫氧推进剂火焰温度
Fig.5 Flame temperature of aluminum-magnesium fuel-rich propellant at different pressure and different environment oxygen content
目前关于压力对推进剂火焰温度影响的相关理论研究较少,观察图5发现环境氧含量相同时火焰温度与压力成正相关。火焰温度随压力增加而增加主要是因为环境氧含量相同时增加压力实际增加了氧浓度,燃烧时反应更剧烈,致火焰温度升高,所以推进剂火焰温度与压力之间的关系应与火焰温度与环境氧含量的之间的关系相似,即火焰温度与压力正相关,结合式(1)可知推进剂火焰温度与压力和环境氧含量线性正相关,式(2)是推进剂火焰温度与压力和环境氧含量采用线性关系拟合的结果
Tf=0.87637φO2+228.21p+1366.23
(2)
式中,Tf为火焰温度,℃;p为绝对压强,MPa。
式(2)与实验结果拟合得较好,图5中黄色半透明面是式(2)的拟合平面。式(2)中含有常数项1366.23 ℃,分析原因可能是,铝镁贫氧推进剂含有一定量的氧化剂,在环境氧含量为0的氛围中也可以短暂燃烧。本次实验测到的铝镁贫氧推进剂在纯氮气氛围、0.1 MPa下燃烧时火焰最高温度为1397.3 ℃,此数值与按式(2)计算得到的数值非常接近,说明可用式(2)计算不同压力和环境氧含量下的铝镁贫氧推进剂火焰温度。
3.3.1 燃速B数理论分析
常用的推进剂燃速公式是Vieille燃速公式和Summerfield燃速公式[2,11],但是这两种公式皆仅考虑了压力对推进剂燃速的影响。Esfahani[16]提出了一种考虑环境氧含量的改进型的阿累尼乌斯燃速公式,但并未考虑压力对推进剂燃速的影响。Quintiere[15]提出了一种综合考虑压力和环境氧含量对含能材料燃速影响的理论——B数理论。B数理论定义了单位质量空气参与燃烧时火焰生成热量与散热量之比,考虑了空气中氧含量及压力变化对推进剂燃烧时产热与散热的影响,得到了不同压力与不同环境氧含量下燃面处能量平衡方程,进而得到燃速与压力和环境氧含量的关系。
根据B数理论,无外加热源加热,推进剂稳定燃烧时推进剂燃面处能量方程为[15]:
(3)
假设推进剂稳定燃烧时火焰形状为圆柱体(此假设与图2稳定燃烧时的火焰形状基本相符),火焰高度按照湍流火焰计算,火焰对推进剂燃面瞬时辐射热流密度qf,r为[17]:
(4)
式中,lm,f为火焰平均辐射长度,m;κf为火焰吸收系数。
压力影响火焰吸收系数κf,κf与压力的平方有关[18]:
(5)
式中,p0为0.1 MPa;κ0为压力等于0.1 MPa时火焰吸收系数。
自然对流换热水平热面向上(冷面向下)时,根据努塞尔数、格拉晓夫数与普朗特数之间的关系[19]得
(6)
式中,k为导热系数,W·m-1·K-1;l为特征长度,m;ρg为火焰密度,kg·m-3;μ为动力粘度,Pa·s。
式(1)、式(5)与式(6)说明qf,r为p1/2φO2的函数[20]; 式(3)右边第一项为火焰对燃面的对流换热热流密度qf,c,根据式(6)与式(3)得qf,c为p1/2与p1/2φO2的函数。假设汽化能不受压力与环境氧含量影响,为常数,图2表明Tv受压力与氧浓度的影响较小,按常数计算。根据以上分析结合式(3)得到不同压力与环境氧含量下铝镁贫氧推进剂燃速经验公式形式为
(7)
式中,C为常数。
不同压力与环境氧含量下镁贫氧推进剂的燃速如图6所示,采用式(7)的形式对实验数据进行回归分析,得式(8)
(8)
图6 不同压力与不同环境氧含量下铝镁贫氧推进剂燃速
Fig.6 Burning rate of aluminum-magnesium fuel-rich propellant at different pressure and different environment oxygen content
3.3.2 压力与氧含量对铝镁贫氧推进剂燃速的影响
式(8)表明压力和氧含量对铝镁贫氧推进剂燃速皆有影响,为分析压力和环境氧含量对铝镁贫氧推进剂燃速影响程度,定义燃速的压力敏感系数αp、燃速的环境氧含量敏感系数αφO2及燃速的压力敏感系数与燃速环境氧含量敏感系数比γ,称γ为燃速压氧敏感系数比,
(9)
(10)
(11)
根据铝镁贫氧推进剂的燃速计算公式(8)得:
(12)
不同压力与环境氧含量下燃速压氧敏感系数比γ计算结果如图7所示,由图7可见,p=0.1 MPa时,γ最大,约为200,随着压力增大,γ不断减小,当压力增加到1.5 MPa时γ下降到约40。但是,任意一实验条件下γ≫1,说明任意一实验条件下压力对燃速的影响大于环境氧含量对燃速的影响。随着压力增大,γ减小,这意味着随着压力的增加,环境氧含量对燃速影响加强,压力对燃速的影响开始相对减弱,分析原因为: 压力对燃速的影响主要是影响推进剂燃烧时火焰“吹离效应”,随着压力增大,吹离效应减弱,火焰对推进剂燃面热反馈增加,从而燃速增加,但当压力增加到一定程度时,压力的增加已经不能明显减弱火焰的“吹离效应”,然而随着压力增大,氧浓度增加,火焰温度上升,火焰辐射对推进剂燃速影响相对加强,从而使的高压下环境氧含量对燃速的影响相对加强。
图7 不同环境氧含量时γ与p的关系
Fig.7 Relationship betweenγandpwith different environment oxygen content
压力和环境氧含量是通过改变火焰对推进剂固相区的热反馈来影响推进剂燃速,设火焰对推进剂固相区反馈热流密度为Λs,假设火焰热反馈过程为一维传热,以推进剂燃面为坐标原点建立动x坐标系,设推进剂初温与环境温度相同,得Λs的表达式为:
(13)
式中,T0为推进剂初温,℃;Ts为推进剂燃面温度,℃; 实验所用铝镁贫氧推进剂ρs=1.64×103kg·m-3,cp=1300 J·kg-1·K-1,Tv=1200 ℃,T∞=25 ℃,计算得:
Λs= 2.505×106(0.01436p1/2φO2+
0.58739p+1.88625)
(14)
式中,Λs的单位为W·m-2。
从式(14)中可以发现压力是影响热反馈的主要因素,这从侧面说明压力是影响铝镁贫氧推进剂燃速的主要因素。
(1) 环境气体氧含量高于推进剂热解产物中氧含量时,铝镁贫氧推进剂点火时气相反应发生在推进剂热解产物与环境气体的扩散区,初现焰远离推进剂表面。但是随着压力增加,扩散区与推进剂表面距离减小。
(2) 激光辐射结束后铝镁贫氧推进剂火焰温度降低,但不会熄火(压力高于0.25 MPa时),70 ms后火焰温度逐渐升高; 推进剂火焰温度与环境氧含量和压力线性正相关。
(3) 环境氧含量和压力对铝镁贫氧推进剂燃速的影响符合B数理论,环境氧含量与压力增加时铝镁贫氧推进剂燃速增加; 通过引入燃速压氧敏感系数比发现压力是影响推进剂燃速的主要因素,但是随着压力增加,压力对火焰“吹离效应”的减弱不明显,致使压力对燃速的影响相对减小; 压力从0.1 MPa增加到1.5 MPa时,燃速压氧敏感系数比从200下降到40。
参考文献:
[1] 彭新. 富燃料固体推进剂瞬态燃烧试验研究[D]. 长沙: 国防科学技术大学,2011.
PENG Xin. Experimental research on transient burning for rich fuel solid propellant[D]. Changsha: National University of Defense Technology, 2011.
[2] 张炜, 朱慧. 铝镁贫氧推进剂低压燃烧性能表征方法研究[J].含能材料,2002,10(3): 125-127.
ZHANG Wei, ZHU Hui. Characterization methods of combustion properties of aluminum-magnesiumfuel-rich propellant at low pressure[J].ChineseJournalofEnergeticMaterials(HannengCailiao), 2002,10(3): 125-127.
[3] 朱国强, 薛谈顺, 周长省, 等. 铝镁贫氧推进剂激光点火特性[J]. 弹道学报, 2013, 25(1): 85-88.
ZHU Guo-jiang, XUE Tan-shun, ZHOU Chang-sheng, et al. Laser ignition characteristics of aluminium-magnesium fuel-rich propellant[J].JournalofBallistics, 2013, 25(1): 85-88.
[4] 陈超, 王英红, 张放利. 铝粉粒径对高铝含量富燃料推进剂一次燃烧性能的影响[J]. 固体火箭技术, 2010, 33(6): 670-674.
CHEN Chao, WANG Ying-hong, ZHANG Fang-li. Effect on primary combustion performance of aluminum particle size in rich-fuel propellant with high alumina content[J].JournalofSolidRocketTechnology, 2010, 33(6): 670-674.
[5] 郑磊, 潘功配, 陈昕, 等. 镁粉粒度对Mg/PTFE贫氧推进剂燃烧性能的影响[J]. 含能材料, 2010,18(2): 180-183.
ZHENG Lei, PAN Gong-pei, CHEN Xin, et al. Effect of magnesium powder particle size on combustion properties of Mg/PTFE fuel-rich propellant[J].ChineseJournalofEnergeticMaterials(HannengCailiao), 2010, 18(2): 180-183.
[6] 聂芝侠, 李葆萱, 胡松启, 等. 镁铝富燃料推进剂燃烧性能研究[J]. 固体火箭技术, 2008, 31(3): 243-246.
NIE Zhi-xia, LI Bao-xuan, HU Song-qi, et al. Study on combustion performance of magnesium-aluminum fuel-rich propellants[J].JournalofSolidRocketTechnology, 2008, 31(3): 243-246.
[7] 冯喜平, 周晓斯, 郑亚, 等. 铝镁推进剂固冲发动机两相燃烧数值模拟[J]. 固体火箭技术, 2011, 34(4): 431-435.
FENG Xi-ping, ZHOU Xiao-si, ZHENG Ya, et al. Numerical simulation of two-phase combustion for Al/Mg-contained propellant solid rocket ramjet[J].JournalofSolidRocketTechnology, 2011, 34(4): 431-435.
[8] 王鸿美, 陈雄, 赵超, 等. NEPE推进剂激光辐照下点火燃烧性能研究[J]. 推进技术, 2015, 36(8): 1262-1267.
WANG Hong-mei, CHEN Xiong, ZHAO Chao, et al. Study on ignition and combustion characteristics of NEPE propellant under laser[J].JournalofPropulsionTechnology, 2015, 36(8): 1262-1267.
[9] 相恒升, 陈雄, 李映坤, 等.一种固体推进剂激光点火实验时的阻焰装置: 中华人民共和国, CN105021765A[P].2015-07-31.
XIANG Heng-sheng, CHEN Xiong, LI Ying-kun, et al. A kind of flameproof device of solid propellant laser ignition: People′s Republic of China, CN105021765A[P]. 2015-07-31.
[10] 徐朝启, 何国强, 刘佩进, 等. 固体推进剂燃烧温度的双波长测试方法[J]. 固体火箭技术, 2010, 33(5): 594-598.
XU Chao-qi, HE Guo-jiang, LIU Pei-jin, et al. A two-wavelength optical method for the high combustion temperature measurement of solid propellant[J].JournalofSolidRocketTechnology, 2010, 33(5): 594-598.
[11] 张炜, 朱慧, 方丁酉, 等. 低压下贫氧推进剂燃烧性能测试方法研究[J]. 含能材料, 1999, 7(3): 118-121.
ZHANG Hui, ZHU Hui, FANG Ding-you, et al. A study on measurement of fuel-rich propellant combustion property at low pressure[J].EnergeticMaterials(HannengCailiao), 1999, 7(3): 118-121.
[12] Kakami A, Hiyamizu R, Shuzenji K, et al. Laser-assisted combustion of solid propellant at low pressures[J].JournalofPropulsionandPower, 2008, 24(6): 1355-1360.
[13] Li J, Litzinger T A. Laser-driven decomposition and combustion of BTTN/GAP[J].JournalofPropulsionandPower, 2007, 23(1): 166-174.
[14] 郭K, 萨默菲尔德, M. 固体推进剂燃烧基础[M]. 北京: 宇航出版社,1998.
Kuo K K, Summerfield M. Fundamentals of soild-propellant combustion[M]. Beijing: CAPH, 1998.
[15] Quintiere J G. Fundamentals of Fire Phenomena[M]. England: J. Wiley, 2006.
[16] Kashani A, Esfahani J A. Interactive effect of oxygen diffusion and volatiles advection on transient thermal degradation of poly methyl methacrylate (PMMA)[J].HeatandMassTransfer, 2008, 44(6): 641-650.
[17] Heskestad G. Luminous heights of turbulent diffusion flames[J].FireSafetyJournal, 1983, 5(2): 103-108.
[18] John L R, Peter K W, Heskestad G. Radiation fire modeling[J].ProceedingsoftheCombustionInstitute, 2000, 28(2): 2751-2759.
[19] 杨世铭, 陶文铨. 传热学[M]. 第4版. 北京: 高等教育出版社,2006.
YANG Shi-ming, TAO Wen-quan. Heat Transfer[M]. 4th ed. Beijing: Higher Education Press, 2006.
[20] Mcallister S, Fernandez-pello C, Urban D, et al. The combined effect of pressure and oxygen concentration on piloted ignition of a solid combustible[J].CombustionandFlame, 2010, 157(9): 1753-1759.