固体燃料冲压发动机自持燃烧的影响因素

巩伦昆,陈雄*,李唯暄,杨海涛,周长省

南京理工大学 机械工程学院,南京 210094

固体燃料冲压发动机自持燃烧的影响因素

巩伦昆,陈雄*,李唯暄,杨海涛,周长省

南京理工大学 机械工程学院,南京 210094

采用试验和仿真方法,研究了固体燃料冲压发动机(SFRJ)自持燃烧的主要影响因素。研究结果表明:SFRJ内为典型的扩散火焰,化学反应特征时间与反应物停留特征时间的比值是影响发动机自持燃烧的关键参数,该值越小,扩散火焰越稳定;本文研究工况下,较小的后台阶高度即可保证发动机实现自持燃烧,后台阶稳定火焰的物理原因是在固体燃料表面附近形成一个低速区,保证了固体燃料的分解燃烧,当台阶高度为0时,由于化学反应特征时间与反应物停留特征时间的比值迅速增加,SFRJ无法自持燃烧;台阶高度较小时,点火过程会对SFRJ的稳定燃烧产生影响,SFRJ可能会出现熄火复燃现象。

冲压发动机;固体燃料;自持燃烧;后台阶;直连式试验;数值仿真

固体燃料冲压发动机(Solid Fuel Ramjet,SFRJ)结构简单、比冲高,在超声速导弹或者增程炮弹中有着较大的应用潜力。但是,SFRJ存在一些不足,阻碍了其在实际工程中的应用。SFRJ对飞行速度、高度等飞行状态比较敏感[1],工作性能在不同飞行状态下的变化较大。在20世纪90年代,国外针对该问题进行了大量的探索研究。Elands等[2]对聚乙烯(Polyethylene,PE)的燃烧性能进行了理论与试验研究;Ferreira等[3]研究了PE在不同结构发动机中的燃烧性能;Zvuloni等[4]研究了小尺寸SFRJ的燃烧性能。而国内在20世纪末开始进行探索研究。其中,郭健[5]对丁羟(Hydroxyl-Terminated Polybutadiene,HTPB)及贫氧推进剂的燃烧性能进行了研究;夏强[6]和成红刚[7]分别对有机玻璃(Polymethyl Methacrylate,PMMA)和PE的燃烧性能进行了仿真与试验研究;谢爱元和武晓松[8]对外侧面结构固体燃料的燃烧性能进行了试验与仿真研究;巩伦昆等[9-10]对固体燃料PE燃速的影响因素进行了分析。而SFRJ较差的火焰稳定性是阻碍其工程应用的另一个重要因素,为了保证SFRJ能够持续稳定工作,往往需要在发动机入口处放置一后台阶结构,如图1所示。SFRJ能否实现自持燃烧是SFRJ设计中需要考虑的关键问题。所谓自持燃烧,是指在点火器工作结束之后,SFRJ中的固体燃料能够持续稳定地燃烧,发动机燃烧室能够建立稳定的压力场。

图1 SFRJ示意图Fig.1 Diagram of SFRJ

针对后台阶结构对SFRJ自持燃烧的作用问题,国外进行了一些试验研究,得到了SFRJ自持燃烧需要的最小台阶高度。Schulte[11]研究了不同燃料内径以及不同来流总温下的SFRJ自持燃烧性能,研究结果表明:当来流总温和燃料内径固定时,自持燃烧需要一个最小突扩比(燃料内径与入口内径的比值),来流总温和燃料内径的增加会使该值减小,当来流总温大于600℃时,无台阶情况下,SFRJ依然可以实现自持稳定燃烧。Netzer和Gany[12]研究了燃料内径dp=10 mm工况下的SFRJ自持燃烧问题,并对比了Schulte的研究结果,结果显示:燃料内径的减小会提高自持燃烧需要的最小突扩比。在仿真研究方面,Elands[13]和Liou[14]等对PE和PMMA固体燃料的自持燃烧条件进行了计算。仿真中,固体燃料燃速均采用固定值,仿真结果显示:随着来流空气质量流率的增加,特征速度下降,当特征速度达到某一固定值时,开始迅速下降,该值被认为是SFRJ自持燃烧的特征速度下限。两位研究者得到的特征速度下限分别为800～900 m/s和775 m/s。国内在这方面的研究较少,朱国强[15]采用数值仿真以及试验方法对SFRJ的点火过程进行了详细的研究,但是,并没有对SFRJ自持燃烧问题进行分析。而在目前公开文献中,并没有其他的相关研究成果。

综上所述,国内外针对SFRJ的自持燃烧问题的研究,仍然存在一些不足之处。通过试验方法,目前已经获得了一些SFRJ自持燃烧条件的规律,但是,针对后台阶结构对SFRJ自持燃烧的影响,并没有进行较深入的分析。在数值仿真研究方面,仍然存在一些问题,在仿真计算中,固体燃料燃速均采用固定值,这与实际情况差别较大,获得的SFRJ自持燃烧条件的可信性有待商榷。

针对目前存在的问题,本文采用地面直连式系统对SFRJ自持燃烧问题进行了试验研究。同时,编制了二维轴对称湍流燃烧程序,对SFRJ的自持燃烧性能进行仿真计算。探索后台阶结构对SFRJ自持燃烧性能的影响规律,并分析影响SFRJ自持燃烧的物理原因。由于目前公开文献中,PE固体燃料的研究成果较多,有较详细的化学反应模型,因此,本文将使用PE作为固体燃料,由于不同固体燃料的差别主要是化学性能的差别,内流场结构基本一致,因此,本文得到的结果对使用其他固体燃料的SFRJ仍然有较大的参考价值。

1 数值计算方法

1.1 控制方程

本文采用有限体积法对Navier-Stokes控制方程进行离散求解。在数值计算中,固体燃料燃速采用阿累尼乌斯定律,其中,影响燃速的主要因素是固体燃料表面温度Tw。因而,为了能够更加准确地计算Tw,本文同时采用数值方法计算了固体温度场,通过保证流固交界面的热流密度以及边界温度相等,实现流固热耦合计算。具体的控制方程、数值方法以及计算方法的验证参考文献[10]。

1.2 化学反应模型

在燃烧室内,固体燃料聚乙烯的分解产物仍然会存在高聚物以及相应的化学反应,由于反应过程非常复杂,很难采用精确模型进行数值仿真计算。简化起见,假设固体燃料分解产物为单质乙烯气体C2H4,固体燃料燃速采用半经验公式,满足阿累尼乌斯定律,表示为[16]

式中:Ru为通用气体常数;Apy为分解反应的指前因子;Eapy为分解反应的活化能;Tw为燃料表面温度;˙r为固体燃料的局部燃速。关于指前因子Apy和活化能Eapy的值,不同文献中差别较大,特别是在SFRJ这种快速升温的环境下的研究结果较少。本文采用文献[17]中的结果,取Apy=8 750 m/s,Eapy=130 kJ/mol。

气相化学反应对火焰稳定有着关键的影响,本文将采用两种化学反应模型对火焰稳定进行预测:一种为层流有限速率(Finite-rate)模型,气相化学反应模型以及化学反应动力学参数见表1[18],表中:A为气相反应的指前因子,n为温度指数,Ea为气相反应的活化能;另外一种为2阶矩湍流燃烧模型(Second-Order Moment Turbulence-Chemistry Model,SOM)[19],该模型考虑了湍流对气相化学反应速率的影响,化学反应速率的计算公式参考文献[19]。

表1 乙烯化学反应模型［18］Table 1 Chemistry reaction model of ethylene［18］

1.3 计算模型

SFRJ计算模型以及边界条件如图2所示,图中没有标明的边界均为壁面边界。其中,阴影部分表示固相计算区域。

仿真中,来流条件保持不变,空气质量流率为0.3 kg/s,来流空气总温为540 K。不同计算工况的燃料内径dp以及入口内径din的数值见表2。表2同时标明了不同计算工况采用的化学反应模型。

仿真中,开始阶段,点火燃气和空气同时从入口处进入,0.5 s之后,点火燃气流量为0,点火过程完成,点火燃气的温度为2 500 K,主要成分为N2、H2O、CO2。大约在1.5 s时刻,流场基本稳定,因此,每个能够稳定燃烧的工况的计算时间为1.5 s。由于计算时间短,燃面退移量很小,因此,仿真中没有考虑燃料表面的网格变形。

图2 SFRJ计算模型Fig.2 Computational model for SFRJ

表2 计算工况Table 2 Computation conditions

2 试验设备及方法

SFRJ直连式系统由空气供应系统、空气加热系统、空气测量控制系统和试验发动机组成。本系统能够模拟海平面马赫数Ma=2.5的工况,空气流量可达到2.0 kg/s。通过航空煤油加热空气的方式,可以使空气的温度达到678 K。燃烧消耗的氧气通过补氧方式补充,使得高温空气的氧气质量百分比始终保持为23%。试验发动机包括入口台阶、燃烧室、补燃室、喷管、点火装置、气流转换装置以及测试系统。

图3为SFRJ试验装置示意图,图4为试验发动机实物照片。气流转换装置主要包括3个电磁阀,用来控制来流空气和氮气的流动。点火装置是利用推进剂的高温喷射气体点燃固体燃料,点火装置位于发动机入口处。在补燃室中,放置一压力传感器,用来测试补燃室内的压强。

图3 SFRJ试验装置示意图Fig.3 Diagram of SFRJ test system

图4 SFRJ试验系统Fig.4 Test system of SFRJ

试验开始阶段,空气和氮气电磁阀关闭,排气阀打开,来流空气通过前端的亚声速壅塞喷管进入试验发动机系统;当来流空气的温度、压强稳定之后,关闭排气阀,打开空气主阀,同时,点火器工作,点燃固体燃料;20 s后,关闭空气主阀,打开排气阀,同时打开氮气阀,吹熄发动机内火焰;最后,关闭供气阀门,试验结束。

本试验模拟了海平面Ma=2.0的工况,因此,来流空气稳定在压强p∞=0.78 MPa,温度T∞=540 K的状态。为了研究入口台阶对SFRJ自持燃烧的影响,试验使用了不同的入口台阶以及固体燃料内径,固体燃料的长度固定为300 mm,空气流量保持0.3 kg/s不变。本研究共进行了9发试验,试验工况如表3所示。

表3 试验工况Table 3 Test conditions

3 SFRJ自持燃烧影响因素分析

目前国外试验研究显示,台阶高度对SFRJ的自持燃烧起着关键作用,但是,对台阶高度影响自持燃烧的物理原因并没有进行更加详细的分析。本文将首先叙述SFRJ自持燃烧性能的仿真与试验结果;随后,根据仿真与试验结果,提出影响SFRJ自持燃烧的主要因素;最后,详细分析台阶高度影响SFRJ自持燃烧的物理原因,并对点火过程对SFRJ自持燃烧的影响进行分析。

3.1 火焰稳定性试验结果

直连式试验结果显示,试验工况1～5、7和8中,固体燃料均能够实现稳定燃烧。由于稳定燃烧情况下的压强变化趋势基本一致,本文仅列举了工况4和5的压强曲线,如图5(a)所示,图中t为试验时间,p为补燃室压强。在前3 s时间内,由于高温点火燃气的存在,固体燃料的燃速会偏高,空燃比降低,使得空燃比更加接近化学当量比,燃烧温度升高,因而,发动机内出现了压力峰值。点火器停止工作之后,燃速下降,温度下降,压强随之下降,随后缓慢上升,在最后5～10 s,发动机基本达到稳定工作状态,压强变化不大。

图5(b)显示了工况6和9中补燃室的压强变化。工况9中,点火器工作结束之后,固体燃料迅速熄火,补燃室压强迅速下降至大约0.3 MPa,即在无后台阶结构情况下,SFRJ无法实现自持燃烧。而工况6中,固体燃料出现了熄火现象,在最后几秒钟出现了复燃现象。

图5 补燃室压强曲线Fig.5 Evolution of pressure in the aft-chamber

3.2 火焰稳定性仿真结果

本文首先分别采用有限速率反应模型和SOM反应模型对工况1～3进行数值仿真。采用有限速率反应模型情况下,工况2、3中,SFRJ能够自持燃烧,而工况1中,火焰无法维持稳定;而采用SOM模型情况下,3种工况中SFRJ均能够自持燃烧,这与直连式试验结果比较吻合。因此,本文随后采用SOM模型对工况4～9进行了仿真计算。仿真结果显示,只有工况6和9中的SFRJ无法自持燃烧。即台阶高度为0时,SFRJ无法自持燃烧,较小的台阶高度依然可以维持火焰的稳定。

对比仿真与试验结果可知,本文建立的仿真模型能够对SFRJ自持燃烧现象进行预测。因此,本文将根据仿真与试验结果,提出影响SFRJ自持燃烧的主要因素,进而分析后台阶影响SFRJ自持燃烧的物理原因。

3.3 SFRJ内流场特性

为了能够获得影响SFRJ自持燃烧的主要因素,首先需要得到SFRJ内的流场结构。图6为din=35 mm,dp=70 mm工况的燃烧室流线图,图中,流线的颜色表示温度T。由图可知,在燃烧室前端的突扩台阶后形成了一回流区,在回流区与来流空气的剪切面上,来流空气与回流区的富燃气体发生燃烧,温度达到最大值;随后,在再发展区,来流空气与分解燃气发生扩散燃烧。

图7为两个轴向位置处,燃烧室内温度以及组分质量分数沿径向的分布,两处分别位于回流区以及再发展区。图中,Δy表示网格点距离固体燃料表面的距离。对比图7(a)和图7(b)可见,两个轴向位置处的组分和温度的分布趋势基本一致,主要差别是,再发展区的火焰面更加靠近固体燃料表面。

由组分和温度的分布可知,在Δy=0.01 m附近,有一层较薄的化学反应区,在该区域,可燃气体与来流空气以接近化学当量比的比例发生扩散燃烧。在化学反应区两侧,存在着富氧和富燃的扩散区域,在靠近壁面的一侧,C2H4从壁面逐渐扩散到化学反应区,而在富氧区域,O2逐渐扩散到化学反应区,化学反应区的反应生成物CO2等组分也逐渐往两侧扩散。但是扩散区有一定的范围,因此,在中心区域基本不存在组分温度的扩散现象。

图6 燃烧室流线图Fig.6 Streamlines in the chamber

图7 燃烧室内温度、组分沿径向的分布Fig.7 Radial distributions of temperature and species in the chamber

3.4 影响SFRJ自持燃烧的主要因素

SFRJ的自持燃烧主要包括两个过程,即固体燃料的分解和分解气体与来流空气的燃烧过程。因而,SFRJ无法自持燃烧的原因可能是由固体燃料分解或者气相化学反应引起的。

表4为试验工况1～5的固体燃料平均燃速˙r,该结果是通过试验前后的称重方法得到的。由表4可见,随着台阶高度的减小,台阶对燃速的影响越来越小,当台阶高度接近0时,固体燃料燃速基本没有变化。由此说明,只要固体燃料在点火阶段能够建立起一定厚度的熔融层,分解反应能够持续进行,固体燃料就不会因为燃速过低导致熄火现象,据此推测,SFRJ的熄火是由气相化学反应导致的。

表4 不同试验工况的平均燃速Table 4 Average regression rates for different test conditions

由3.3节分析可知,在SFRJ内,固体燃料分解燃气与来流空气发生扩散燃烧,燃烧释放的热量反馈给固体燃料,实现固体燃料的自持燃烧。当燃烧无法提供足够的热量反馈给固体燃料时,则会出现熄火现象。一般研究认为,稳定扩散火焰的建立主要由化学反应特征时间与反应物停留时间决定[20]。当空气流速过高时,由于反应物停留时间较短,反应物无法充分反应,没有足够的热量反馈给固体燃料,最终导致发动机熄火。Matsuoka等[20]对小口径PMMA固体燃料管道内的火焰传播问题进行了研究,并定义达姆科勒数(Da)为反应气体混合特征时间与化学反应特征时间的比值。随着Da的减小,化学反应由传热传质物理过程转变为化学反应动力学控制,最终出现熄火现象。

基于上述分析,本文定义一个无量纲变量φ,该变量表示化学反应特征时间τch与反应物停留时间τres的比值。反应物停留时间定义为

τres=L/｜u｜ (2)

式中:L为特征长度,本文取L为药柱长度0.3 m;｜u｜为气流速度值。由3.2节的仿真结果可见,有限速率反应模型预测的化学反应速率较低,当采用考虑湍流影响的SOM模型时,仿真模型预测结果与试验结果比较吻合,由此说明,在本文研究的工况范围内,化学反应速率受到湍流的影响较大。在马格努斯提出的涡耗散湍流燃烧模型中[21],化学反应特征时间τch由k/ε决定,其中,ε为湍动能耗散率,k为湍动能。参考涡耗散模型以及仿真结果,本文采用湍流特征时间来表示化学反应特征时间。由于本文采用k-ωSST湍流模型,同时,比动能耗散率ω与k、ε存在如下转换关系[22]:

ε=βωk (3)式中:常数β=0.09。因此,化学反应特征时间可表示为

τch=1/(βω) (4)

综上所述,无量纲变量φ定义为

根据φ的定义,φ值越小,化学反应特征时间越短,或者反应物停留时间越长,反应越充分,能够为固体燃料提供足够的热反馈,扩散火焰越稳定。需要注意的是,上述定义会受到研究工况的影响,本文中,来流空气总温不变化,主要研究了台阶高度变化对SFRJ自持燃烧的影响。当来流总温有较大幅度提高时,气相化学反应速率会有明显的提高,化学反应特征时间减小,φ减小,火焰更容易稳定,如文献[11]中的研究结果,当来流总温达到600℃时,无台阶情况下,SFRJ依然可以自持燃烧。

3.5 台阶高度变化对φ的影响

3.4节提出了影响SFRJ自持燃烧的主要因素,本节将以此为依据,分析后台阶影响SFRJ自持燃烧的物理原因。

由3.1节和3.2节的仿真与试验结果可见,当台阶高度为0时,SFRJ无法自持燃烧,较小的台阶高度,即5～10 mm的台阶高度即可实现SFRJ自持燃烧。后台阶稳定火焰的原因可能是:后台阶的存在会使固体表面附近产生一个气流低速区,该低速区域的存在使得反应物停留时间增加,而化学反应特征时间变化不大,因而φ值减小,从而保证了固体燃料的分解以及气相化学反应的初步进行,防止高速的来流空气直接终止化学反应过程。随着台阶高度的增加,固体表面附近的低速区增加,另外,台阶高度的增加也使得湍流强度增加,促进了化学反应速率的增加,这两个因素都促进了φ值的减小,因而火焰更加稳定。

为了更加直接明了地了解台阶高度对φ的影响,本文将定量描述不同工况下φ的分布。图8为仿真工况7和9中,不同时刻的φ沿轴向x的分布,其中,每个轴向位置的φ值根据式(5)计算得到,式中,速度等自变量的值均为该轴向位置处火焰面,即温度最大点处的值。由于工况9中,SFRJ不能够自持燃烧,在点火燃气结束的0.1～0.2 s内即出现了完全熄火现象,因此,图8(b)显示了点火燃气刚刚结束之后的3个时刻,即t=0.51,0.53,0.55 s时刻的φ分布。从图8可以看出,同一工况中,3个不同时刻的φ数值差别不大,分布规律基本相似,沿着轴向方向,φ值逐渐增加,火焰稳定性相应降低。

图8 不同时刻的φ沿轴向的分布(仿真结果)Fig.8 Distributions ofφalong axial direction at different time instants(simulation results)

图9为不同仿真工况下的φ沿轴向的分布。对于能够稳定燃烧的工况,φ值为1.5 s时刻的分布。图9(a)为燃料内径为50 mm工况的φ分布,其中,对于工况9,图9(a)仅仅显示了0.55 s时刻的φ值分布规律。由图9(a)可以看到,工况7和8的φ值差别不大,而工况9的φ值大约为工况7和8的2～4倍。φ值越大,意味着反应物停留时间相对于化学反应时间更短,使得化学反应无法充分完成,没有足够的热量反馈到固体燃料表面,最终出现熄火现象。由此可见,φ值的增大是导致SFRJ无法自持燃烧的关键原因。

图9(b)为燃料内径为70 mm工况的φ沿轴向的分布。对于工况6,φ值根据0.55 s时刻的流场计算得到。φ的分布规律与燃料内径为50 mm的工况基本相似,当台阶高度为0时,φ值迅速增加,这是工况6中的SFRJ无法自持燃烧的主要原因。

图9 不同工况下φ沿轴向的分布(仿真结果)Fig.9 Distributions ofφalong axial direction at different conditions(simulation results)

3.6 点火过程对SFRJ自持燃烧的影响

由3.4节的分析可知,影响SFRJ自持燃烧的主要因素是化学反应特征时间与反应物停留时间的比值。但是,该结论的前提是,点火燃气必须能够提供足够的能量,使得固体燃料在点火阶段熔融分解燃烧,固体燃料表面附近建立起一定厚度的熔融分解层,从而保证点火结束后,固体燃料能够持续分解燃烧。

由3.1节图5可知,在试验工况6中,固体燃料在点火器工作结束之后出现了熄火,大约7～8 s之后,固体燃料再次复燃,补燃室压强再次升高。产生这种现象的原因可能是:台阶高度较低的情况下,固体燃料表面附近的流速相对较高,对点火阶段的高温燃气对固体燃料的加热有负面的影响,固体燃料高温分解区较薄,没有形成稳定的温度场。点火器结束工作之后,固体燃料表面附近,因为有低速流动区域的存在,可能仍然存在缓慢的分解放热反应,对固体燃料持续加热,但是该加热过程较慢,大约7～8 s后,固体燃料形成了稳定的温度场,随后,固体燃料稳定持续地分解燃烧,燃烧室压强上升。在试验工况8中,台阶高度同样较小,但是没有出现这种现象,由此说明,点火阶段的高温燃气对固体燃料的加热,以及固体燃料表面区域稳定温度场的形成不仅仅受到台阶高度的影响,还受到发动机结构尺寸,来流条件等参数的影响,不确定性因素较多。

对于dp=70 mm和din=60 mm的工况,仿真结果显示,SFRJ能够自持燃烧,而试验中出现了先熄火后复燃的情况,本文的仿真模型很难对该现象进行预测。这主要有以下几方面原因:① 本文仿真研究中,点火过程与试验中的点火过程存在一些不同,这导致本文仿真模型无法对试验点火阶段的固体燃料的传热过程进行准确的描述;② 由于聚乙烯的分解燃烧过程复杂,很难进行准确的建模,简化的化学反应模型无法对聚乙烯表面的分解以及初始燃烧过程进行准确的描述,而该过程对固体燃料的熄火再复燃现象有着重要的影响;③ 二阶矩湍流燃烧模型对湍流与化学反应的相互作用的描述是不足的,而固体燃料附近的熄火复燃现象比较复杂,需要更加准确的数学模型,比如概率密度函数等,来描述湍流与化学反应的相互作用。

4 结 论

本文采用地面直连式试验系统对SFRJ的自持燃烧现象进行了研究,并编制了二维轴对称湍流燃烧程序,对SFRJ自持燃烧的影响因素进行了分析。获得了以下结论:

1)SFRJ内为典型的扩散燃烧现象,化学反应特征时间与反应物停留时间的比值φ是影响SFRJ自持燃烧的关键因素,该值越小,火焰越稳定。

2)较小的台阶高度即可保证SFRJ能够自持燃烧,后台阶稳定火焰的主要原因是:在固体燃料表面形成一个气流低速区,使得反应物停留时间足够长,保证了固体燃料的分解以及初步气相反应的进行;无台阶情况下,φ值迅速增大,从而造成SFRJ无法自持燃烧。

3)台阶高度较低时,点火过程会影响固体燃料表面的传热分解过程,使得固体燃料无法建立一定的熔融分解层,SFRJ可能出现熄火复燃现象,该过程需要更加准确的计算模型进行描述预测。

[1] 鲍福廷,黄熙君,张振鹏,等.固体火箭冲压发动机[M].北京:中国宇航出版社,2006:25-31.BAO F T,HUANG X J,ZHANG Z P,et al.Integral solid propellant ramjet rocket motor[M].Beijing:China Astronautic Publishing House,2006:25-31(in Chinese).

[2] ELANDS P J M,KORTING P G,WIJCHERS T,et al.Comparison of combustion experiments and theory in polyethylene solid fuel ramjets[J].Journal of Propulsion and Power,1990,6(6):732-739.

[3] FERREIRA J G,CARVALHO J A,SILVA M G.Experimental investigation of polyethylene combustion in a solid fuel ramjet:AIAA-1996-2698[R].Reston:AIAA,1996.

[4] ZVULONI R,LEVY Y,GANY A.Investigation of a small solid fuel ramjet combustor[J].Journal of Propulsion and Power,1989,5(3):269-275.

[5] 郭健.固体燃料冲压发动机工作过程理论与试验研究[D].长沙:国防科学技术大学,2007:42-70.GUO J.Theoretical and experimental investigation on the working process of solid fuel ramjet[D].Changsha:National University of Defense Technology,2007:42-70(in Chinese).

[6] 夏强.固体燃料冲压发动机工作过程研究[D].南京:南京理工大学,2011:82-112.XIA Q.Investigation on the working process of solid fuel ramjet[D].Nanjing:Nanjing University of Science and Technology,2011:82-112(in Chinese).

[7] 成红刚.固体燃料冲压发动机燃烧性能数值仿真及试验研究[D].南京:南京理工大学,2013:38-110.CHENG H G.Numerical simulation and experimental research on combustion characteristics of solid fuel ramjet[D].Nanjing:Nanjing University of Science and Technology,2013:38-110(in Chinese).

[8] 谢爱元,武晓松.外侧面燃烧固体燃料冲压发动机燃烧室流场的数值研究[J].推进技术,2014,35(7):956-964.XIE A Y,WU X S.Numerical simulation on flow field of solid fuel ramjet combustion chamber with outboard burning[J].Journal of Propulsion Technology,2014,35(7):956-964(in Chinese).

[9] 巩伦昆,陈雄,周长省,等.结构尺寸对固体燃料冲压发动机燃速影响的仿真研究[J].兵工学报,2016,37(5):798-807.GONG L K,CHEN X,ZHOU C S,et al.Numerical investigation on the effect of solid fuel ramjet geometry on solid fuel regression rate[J].Acta Armamentarii,2016,37(5):798-807(in Chinese).

[10] 巩伦昆,陈雄,周长省,等.来流条件对SFRJ燃速及自持燃烧性能的影响[J].航空学报,2016,37(5):1428-1439.GONG L K,CHEN X,ZHOU C S,et al.Effect of inlet flow condition on regression rate and self-sustained combustion of solid fuel ramjet[J].Acta Aeronautica et Astronautica Sinica,2016,37(5):1428-1439(in Chinese).

[11] SCHULTE G.Fuel regression and flame stabilization studies of solid-fuel ramjets[J].Journal of Propulsion and Power,1986,2(4):301-304.

[12] NETZER A,GANY A.Burning and flameholding characteristics of a miniature solid fuel ramjet combustor[J].Propellants,Explosives,Pyrotechnics,1991,7(3):357-363.

[13] ELANDS P,DIJSWIJK F,ZANDBERGEN B.Experimental and computational flammability limits in a solid fuel ramjet:AIAA-1990-1964[R].Reston:AIAA,1990.

[14] LIOU T M,LIEN W Y,HWANG P W.Flammability limits and probability density functions in simulated solidfuel ramjet combustors[J].Journal of Propulsion and Power,1997,13(5):643-650.

[15] 朱国强.固体燃料冲压发动机点火过程研究[D].南京:南京理工大学,2014:79-111.ZHU G Q.The research on the ignition process of solid fuel ramjet[D].Nanjing:Nanjing University of Science and Technology,2014:79-111(in Chinese).

[16] STOLIAROV S I,WALTERS R N.Determination of the heats of gasification of polymers using differential scanning calorimetry[J].Polymer Degradation and Stability,2008,93(2):422-427.

[17] GASCOIN N,FAU G,GILLARD P.Flash pyrolysis of high density polyethylene:AIAA-2013-3833[R].Reston:AIAA,2013.

[18] BAURLE R A,MATHUR T,GRUBER M R,et al.A numerical and experimental investigation of a scramjet combustor for hypersonic missile applications:AIAA-1998-3521[R].Reston:AIAA,1998.

[19] ZHOU L X,CHEN X L,ZHENG C G,et al.Second-order moment turbulence-chemistry models for simulating NOxformation in gas combustion[J].Fuel,2000,79(11):1289-1301.

[20] MATSUOKA T,MURAKAMI S,NAGATA H.Transition characteristics of combustion modes for flame spread in solid fuel tube[J].Combustion and Flame,2012,159(7):2466-2473.

[21] MAGNUSSENA B F,HJERTAGERA B H.On mathematical modeling of turbulent combustion with special emphasis on soot formation and combustion[J].Symposium on Combustion,1977,16(1):719-729.

[22] WILCOX D C.Dilatation-dissipation corrections for advanced turbulence models[J].AIAA Journal,1992,30(11):2639-2646.

lnfluencing factors for self-sustained combustion of solid fuel ramjet

GONG Lunkun,CHEN Xiong*,Ll Weixuan,YANG Haitao,ZHOU Changsheng

School of Mechanical Engineering,Nanjing University of Science and Technology,Nanjing 210094,China

By simulation and test,the parameters which influence the self-sustained combustion of the solid fuel ramjet(SFRJ)are investigated.The results show that the combustion in SFRJ is diffusion-controlled,and the ratio of characteristic chemical reaction time to characteristic residence time is the critical parameter influencing the self-sustained combustion of ramjet.The smaller the ratio is,the more stable the flame will be.Under the conditions of the present research,a small step height can guarantee the self-sustained combustion of ramjet.The physical reason for the stabilization of combustion behind a backward-facing step is that existence of a low speed flow area near the solid fuel surface guarantees the pyrolysis and combustion of the solid fuel.When the step height is 0,the flame will go out as the ratio of characteristic chemical reaction time to characteristic residence time increases quickly.When the step height is small,the ignition process may have influence on stable combustion of the SFRJ,and reburning may occur after flameout.

ramjet;solid fuel;self-sustained combustion;backward-facing step;direct-connect test;numerical simulation

2016-09-28;Revised:2016-10-25;Accepted:2016-11-16;Published online:2016-12-14 16:29

URL:www.cnki.net/kcms/detail/11.1929.V.20161214.1629.006.html

General Armaments Department Pre-research Foundation(404040301)

V235

A

1000-6893(2017)07-120821-10

10.7527/S1000-6893.2016.0301

2016-09-28;退修日期:2016-10-25;录用日期:2016-11-16;网络出版时间:2016-12-14 16:29

www.cnki.net/kcms/detail/11.1929.V.20161214.1629.006.html

总装备部预先研究项目(404040301)

*通讯作者.E-mail:chenxiong@njust.edu.cn

巩伦昆,陈雄,李唯暄,等.固体燃料冲压发动机自持燃烧的影响因素[J].航空学报,2017,38(7):120821.GONG L K,CHEN X,Ll W X,et al.lnfluencing factors for self-sustained combustion of solid fuel ramjet[J].Acta Aeronautica et Astronautica Sinica,2017,38(7):120821.

(责任编辑:王娇)

*Corresponding author.E-mail:chenxiong@njust.edu.cn