等离子体强化混合燃气燃烧研究

总636期第四期2018年4月河南科技Henan Science and Technology等离子体强化混合燃气燃烧研究纪海龙1冯喜平2（1.中国空空导弹研究院，河南洛阳471009；2.西北工业大学燃烧、热结构和内流场重点实验室，陕西西安710072）摘要：为验证等离子体对固冲发动机中气相组分燃烧的强化特性，通过建立等离子体助燃理论模型，对有无等离子体条件下的混合燃气在固冲补燃室中的燃烧过程进行了数值和实验研究。

研究结果表明：在燃烧室中加入等离子体后，有效促进了混合燃气的链式反应，缩短了燃烧时间；等离子体显著提高了混合燃气的燃烧效率，燃烧效率增加了16.6%。

与实验结果对比可知，数值计算结果误差在10%以内，说明所建立的理论模型较为准确地模拟了混合燃气燃烧过程，具有一定的可靠性。

关键词：等离子体；混合燃气；固冲发动机；燃烧效率中图分类号：V35.12文献标识码：A文章编号：1003-5168（2018）10-0139-04Study on the Plasma Enhancing Combustion of Mixed GasJI Hailong 1FENG Xiping 2（1.China Airborne Missile Academy ，Luoyang Henan 471009；2.Key Laboratory of Combustion,Flow and Thermal-Structure,Northwestern Polytechnical University ，Xi ’an Shaanxi 710072）Abstract:In order to verify the enhancement characteristics of plasma in gas-phase combustion of the solid ramjet rocket,Numerical and experimental studies were carried out by establishing a theoretical model of plasma-assisted combustion.The result indicated that:the addition of plasma in the combustion chamber effectively promoted thechain reaction of the mixed gas and shortens the burning time;the plasma significantly increased the combustion effi⁃ciency of the mixed gas,and the combustion efficiency increased by 16.6%.Compared with the experimental results,the error of numerical calculation was less than 10%which indicated that the theoretical model was relative accurate⁃ly simulate the combustion process of mixed gas and has a certain reliability.Keywords:plasma ；mixed gas ；the solid ramjet rocket ；the combustion efficiency 相比于镁、铝等金属，硼是固冲发动机推进剂的理想高能组分。

但是，硼表面往往存在一层致密的氧化层B 2O 3，且硼粒子具有较高的熔点和沸点，使得硼颗粒点火燃烧困难，抑制着硼所含能量的有效释放，这对固体火箭发动机来说是非常不利的［1，2］。

目前，改善硼点火和燃烧的措施主要有调整燃烧室结构、调节推进剂配方和硼粒子包覆等。

但是，这些措施均会带来新的问题［3-6］。

等离子体助燃技术是一种极具潜力的动力学强化方式［5，7］。

早在20世纪80年代，等离子体助燃技术就被用于解决超音速燃烧中的点火困难和不完全燃烧等问题。

2003年，Lou 等［8］人研究发现，等离子体可以加快乙烷/氧化物的反应速率，提高热量的释放速度。

2011年，SergeyB.Leonov 等［9］通过等离子体助燃技术成功实现了在高马赫气流下CH 4/H 2/空气的点火，并使混合火焰稳定燃烧。

2014年，宋振兴［10］等人研究了等离子体助燃的动力学机理，得到氮系和氧系粒子浓度在不同放电时间和频率条件下的变化规律。

等离子体助燃技术发展迅速，主要用于航空航天所涉及的超音速燃烧领域，在固冲发动机中的应用鲜有报道。

由于固冲发动机含硼富燃燃气不仅含有气相组分，而且含有固体颗粒，研究难度较大。

因此，本文先以固冲发动机一次燃气中的气相组分为研究对象，通过数值计算和点火实验，验证等离子体对补燃室中混合燃气燃烧的强化特性，以此为后续的等离子助燃含硼混合燃气实收稿日期：2018-03-07作者简介：纪海龙（1991—），男，硕士，助理工程师，研究方向：含硼富燃燃气燃烧与流动。

能源与化学·140·第4期验和理论研究提供理论支撑。

1数学物理模型1.1基本假设补燃室中混合燃气燃烧流动非常复杂，为了简化计算，在计算过程中进行如下简化：①混合燃气为理想气体，满足理想气体的状态方程，不考虑辐射换热和体积力；②认为燃气与壁面无热交换，燃气流动为绝热准定常流动；③混合燃气中不可燃烧或质量分数小的成分用氮气代替。

1.2控制方程控制方程为笛卡尔坐标系下的Navier Stokes 方程，其形式为：∂Q ∂t +∂E ∂x +∂F ∂y +∂G ∂z=∂H x ∂x +∂H y ∂y +∂H z ∂z +I （1）其中，H x 、H y 、H z 为粘性通量，I 为源项，Q 、E 、F 、G 等参数具体意义可以参考相关流体力学文献。

1.3燃烧模型混合燃气与空气的反应具备一般同心射流火焰的特征，燃烧室前半段同时存在层流和湍流，后段充分发展为湍流，因此燃烧模型采用有限速率/涡耗散模型描述。

1.4等离子体助燃模型为了方便计算，将等离子体以源项（若干自由基）的形式加入燃烧室中，且保证均匀分布。

将等离子体作为源项引入Navier Stokes 控制方程中，可得出：∂∂t ()ρE +∂∂x i ()u i()ρE +p =∂∂x i æèçöø÷k eff∂T ∂x i -∑j 'h 'j J 'j +u j (τij )eff （2）其中，k eff 为传热系数；J 'j 为扩散流量、h 'j 为显焓；u i ，u j 分别速度分量。

活性基的密度变化采用Boltzmann零阶矩方程表示；动量和能量方程分别采用Boltzmann 一阶矩和二阶矩方程表示。

1.5物理模型补燃室中混合燃气燃烧流动非常复杂。

为了突出补燃室中的湍流扩散燃烧火焰，建立了如图1所示的物理模型，混合燃气由长尾管进入补燃室，空气由燃气入口外部的环形入口进入，等离子体从距离燃气出口15mm 处以源项添加，产生的尾气由喷管出口排出。

空气入口等离子体源图1计算模型与网格1.6边界条件以同等条件下通过固冲燃气发生器点火实验得到的一次燃气主要组分作为计算输入的混合燃气组分，燃气具体成分及比例如表1所示。

混合燃气流量为0.01kg/s ，空气流量为0.3kg/s ，等离子体工质流量为0.001kg/s ，具体边界条件如表2所示。

表1混合燃气主要成分表2边界条件参数空气入口燃气入口喷管出口压力/MPa -0.7010.101温度/K 300-300质量流率（kg/s ）0.300.01-2计算结果图2（a ）和图2（b ）分别为未加等离子体和加活性基后燃烧温度分布云图。

从图2可知，两种工况下的燃烧室火焰形状基本相同，均呈纺锤形，具有纵向受限扩散火焰结构的一般特征；混合火焰最高温度在燃烧室前部，靠近混合燃气出口位置，最低温度为来流空气温度300K 。

对比图2（a ）与图2（b ）两种工况可知，燃烧室中加入等离子体后会对混合燃气产生积极影响，火焰最高温度为2600K ，明显高于图2（a ）中未加等离子工况的2300K ，且火焰锋面沿纵向变宽。

这主要是因为活性粒子诱发了混合燃气与空气之间的链式反应，加快了混合燃气的反应，释放了更多热量。

（a ）未加等离子体工况（b ）加活性基图2燃烧室温度分布云图图3（a ）和图3（b ）分别为未加等离子体和加活性基燃烧室压力分布云图。

从图3可以看出，两种工况下的燃烧室压力分布均匀，只在喷管处压力变化较大。

燃烧室中加入活性粒子后，压强为0.225MPa ，与未加等离子体工况燃烧室压强0.161MPa 相比，提高了39%。

可见，等离子体会对混合燃气产生积极影响，提高了混合燃气的燃烧效率。

（a ）未加等离子体工况等离子体强化混合燃气燃烧研究第4期·141·（b ）加O 活性基图3燃烧室压力分布云图图4为未加等离子体和加活性基两种工况下燃烧室中OH 浓度沿轴线变化曲线。

从图4可知，两种工况下的OH 浓度变化趋势基本相同，均在距离燃气入口处约0.3m 达到最大，然后约在0.44m 处降低为0。

这说明混合燃气在燃烧室前半部分即已燃烧完全。

需要说明的是，图中显示的加等离子体工况下OH 浓度已经除去了由外部加入的那部分OH 自由基。

可以看出，该工况下OH 浓度明显高于未加等离子工况燃烧室中的OH 浓度。

这是因为加入O 等活性粒子后，会与H 2和CH 4等燃气快速反应（O+H 2=OH+H ，H+HO 2=2OH ，O+CH 4=CH 3+OH 等），促进了混合燃气燃烧链式反应，加快了反应进程。

OH 浓度显著增加，混合燃气燃烧完全后又降低为0。

OH 被认为是链式反应的活化中心，OH 浓度增加也从另一个方面说明等离子体对混合燃气燃烧具有强化作用。

Plasma NO-Plasma图4OH 浓度沿轴线变化曲线3实验验证3.1实验系统及实验条件实验系统主要由燃烧室、来流模拟系统、混合燃气供应系统、等离子体系统、数据采集和控制系统以及火焰拍摄装置等组成，燃烧实验装置如图5所示。

实验工况包括混合燃气成分、含量等参数与仿真边界条件相同，其中来流空气流量为0.3kg/s ，总温为300K ，总压为0.27MPa ；混合燃气流量为0.01kg/s ，相应空燃比为30。

混合燃气由储气瓶经管路输送到燃气发生器中，在燃烧室中与来流空气发生扩散燃烧，燃烧产物经尾喷管排出。