CHINESE JOURNAL OF ENERGETIC MATERIALS
+Advanced Search
en
×

分享给微信好友或者朋友圈

使用微信“扫一扫”功能。
参考文献 1
丁则胜,邱光纯,刘亚飞, 等. 固体燃料底部排气空气动力研究[J]. 空气动力学学报, 1991, 9(3): 300-307
DINGZe‑sheng, QIUGuang‑chun, LIUYa‑fei, et al. An aerodynamic investigation of base bleed by solid fule[J]. Acta Aerodynamic Sinica, 1991, 9(3): 300-307
参考文献 2
JacksonD E,LinnC A B,MoserM D. Ultrasonic burn rate testing of bi-plateau propellants under decreasing pressure conditions[C]//41st AIAA/ASME/SAE/ASEE Joint Propulsion Conference & Exhibit. Tucson,AZ,USA,2005,AIAA2005-4164.
参考文献 3
YuY, ZhouY, LuC, et al. Study on unsteady combustion behaviors of AP/HTPB base bleed propellants under transient depressurization conditions[J]. Propellants ,Explosives, Pyrotechnics, 2014, 39(4): 511-517.
参考文献 4
陆春义, 周彦煌, 余永刚. 高降压速率下复合底排药剂瞬变燃烧特性研究[J]. 含能材料, 2007, 15(6):587-591.
LUChun‑yi, ZHOUYan‑huang,YUYong‑gang. Combustion of composite base bleed charge under rapid depressurization [J]. Chinese Journal of Energetic Materials(Hanneng Cailiao),2007, 15(6) : 587-591.
参考文献 5
陆春义,周彦煌,余永刚.底排点火具在高降压速率下瞬态燃烧特性的实验研究[J].含能材料, 2008, 16(5): 629-632.
LUChun‑yi, ZHOUYan‑huang,YUYong‑gang. Transient combustion performance of igniter for base bleed unit under rapid depressurization[J]. Chinese Journal of Energetic Materials(Hanneng Cailiao) , 2008, 16(5): 629-632.
参考文献 6
陆春义. 底排装置强非稳态燃烧特性研究[D]. 南京:南京理工大学,2009.
LUChun‑yi. Study on characteristics of unstable combustion of base bleed unit[D]. Nanjing: Nanjing University of Science and Technology,2009.
参考文献 7
赵宏立, 刘来东, 靳建伟, 等. 压力下降条件下一种双基发射药的瞬态燃烧特性[J]. 火炸药学报, 2012, 35(1): 64-68.
ZHAOHong‑li, LIULai‑dong, JINJian‑wei, et al. Transient combustion performance of gun propellant under depressurization condition[J]. Chinese Journal of Explosive & Propellants, 2012, 35(1): 64-68.
参考文献 8
YeR, YuY, CaoY. Experimental study of transient combustion characteristics of AP/HTPB base bleed propellant under rapid pressure drop[J]. Combustion Science & Technology, 2015, 187(3):445-457.
参考文献 9
曹永杰, 余永刚, 张明安, 等. 底排推进剂瞬态泄压工况下的燃烧失稳特性[J]. 燃烧科学与技术, 2015(05): 421-427.
CAOYong‑jie, YUYong‑gang, ZHANGMing‑an,et al. Combustion instability characteristics of base bleed propellant under transient depressurization[J]. Journal Combustion Science and Technology, 2015(05): 421-427.
参考文献 10
曹永杰, 余永刚, 叶锐, 等, 底排推进剂瞬态泄压工况下燃烧流场特性的数值模拟[J]. 含能材料, 2013(4): 464-468.
CAOYong‑jie, YUYong‑gang, YERui, et al. Numerical simulation of combustion flow field characteristics of base bleed propellant under transient pressure‑release[J]. Chinese Journal of Energetic Materials(Hanneng Cailiao), 2013(4): 464-468.
参考文献 11
曹颖, 王栋, 武晓松. 底排装置瞬态泄压过程数值模拟[J]. 四川兵工学报, 2014, 35(5): 43-46.
CAOYing, WANGDong, WUXiao‑song. Numerical simulation of transient depressurization process of base bleed[J]. Sichuan Acta Armamentarii, 2014, 35(5): 43-46.
参考文献 12
ZhuoC, FengF, WuX, et al. Numerical simulation of the muzzle flows with base bleed projectile based on dynamic overlapped grids[J]. Computers & Fluids, 2014, 105:307-320.
参考文献 13
ChoiJ Y , ShinE, KimC K. Numerical study of Base-Bleed projectile with external combustion[R]. AIAA 2005-4352, 2005.
参考文献 14
XueX, YuY. An improvement of the base bleed unit on base drag reduction and heat energy addition as well as mass addition[J]. Applied Thermal Engineering, 2016, 109: 238-250.
参考文献 15
余文杰, 余永刚. 底排装置尾部化学非平衡流的数值模拟[J]. 推进技术, 2015, 36(11): 1610-1615.
YUWen‑Jie , YUYong‑gang. Numerical simulation of base flow with chemical non‑equilibrium for base bleed equipment[J]. Journal of Propulsion Technology, 2015, 36(11): 1610-1615.
参考文献 16
丁则胜, 罗荣, 陈少松, 等. 底部燃烧减阻性能的若干参数影响研究[J]. 弹道学报, 1996, 8(4): 79-83.
DINGZe‑sheng, LUORong, CHENShao‑song, et al. A study of some parameters influence on performance of drag reduction by base burning[J]. Journal of Ballistics, 1996, 8(4):79-83.
参考文献 17
余文杰, 余永刚. 二次燃烧对底排装置尾部流场影响的数值模拟[J]. 爆炸与冲击, 2015, 35(1): 94-100.
YUWen‑Jie, YUYong‑Gang. Numerical simulation of secondary combustion affecting base flow of base bleed equipment[J]. Explosion & Shock Waves, 2015, 35(1): 94-100.
参考文献 18
MenterF R. Two‑equation eddy‑viscosity turbulence models for engineering applications[J]. AIAA Journal, 1994, 32(8):1598-1605.
参考文献 19
朱自强. 应用计算流体力学[M]. 北京:北京航空航天大学出版社, 1998: 43-44.
ZHUZi‑qiang. Applied Computational Fluid Dynamics[M]. Beijing University of Aeronautics and Astronautics Press, 1998: 43-44.
参考文献 20
吴子牛. 计算流体力学基本原理[M]. 北京:科学出版社, 2001:48-50.
WUZi‑Niu. Basic Principles of Computational Fluid Dynamics [M]. Bejing:Science China Press, 2001: 48-50.
参考文献 21
武频, 赵润祥, 郭锡福. 弧长网格生成法及其应用[J]. 南京理工大学学报, 2002, 26(5): 482-485.
WUPin, ZHAORun‑xiang, GUOXi‑fu. Arc length method of grid generation and its application[J]. Journal of Nanjing University of Science and Technology, 2002, 26(5): 482-485.
参考文献 22
GibelingH J, BuggelnR C. Reacting flow models for navier‑stokes analysis of projectile Base Combustion[R]. AIAA 91-2077, 1991.
参考文献 23
KuoK K, SummerfieldM. Fundamentals of solid‑propellant combustion[M]. Beijing: China Astronautic Publishing House, 1994.
参考文献 24
刘君,张涵信,高树椿. 一种新型的计算化学非平衡流动的解耦方法[J]. 国防科技大学学报,2000,22(5):19-22
LIUJun, ZHANGXin‑han, GAOShu‑chun. A new uncoupled method for numerical simulation of nonequilibrium flow[J]. Journal of National University of Defense Technology, 2000, 22(5): 19-22.
参考文献 25
梁德旺, 王可. AUSM+格式的改进[J]. 空气动力学学报, 2004, 22(4): 404-409.
LIANGDe‑Wang, WANGKe. Improvement of AUSM+ scheme[J]. Acta AerodynamicSinica, 2004, 22(4): 404-409.
参考文献 26
PulliamT H. Time Accuracy and the use of Implicit Methods[R]. AIAA 93-3360, 1993.
目录 contents

    摘要

    为了研究瞬态泄压过程中二次燃烧对底排弹尾部加能的影响,建立了底排装置瞬态泄压过程中尾部化学非平衡流的数学物理模型。采用对流迎风矢通量分裂格式(AUSM+)的改进格式、切应力运输湍流模型SST kω模型和8组分12步反应的H2─CO反应动力学机理,运用统一算法的思路编程求解二维轴对称Navier‑Stokes方程,对尾部流场进行数值模拟,得到泄压过程中底排燃烧室压力随时间的变化规律,模拟结果和文献实验结果基本吻合。在此基础上对泄压过程中有无二次燃烧两种条件下的尾部流场进行数值预测,分析了泄压过程中底排装置尾部流场特征参数随时间的变化规律。结果表明,泄压过程前期,二次燃烧对尾部流场影响较小。而在泄压中后期,二次燃烧对尾部流场影响较大,此时,尾部流场特征参数更快趋于稳定,且二次燃烧使尾部流场温度提高,加能效果显著,底压增大,底排弹底阻明显减小,减阻率比不含二次燃烧的减阻率提高了75%。

    Abstract

    To investigate the effect of post‑combustion on the energy addition at tail of base bleed projectile in the transient depressurization process, a mathematical and physical model of the chemical non‑equilibrium flow at tail in the transient depressurization process for base bleed equipment was established. Two‑dimensional axial symmetric Navier‑Stokes equations were solved by idea programming of uniform algorithm using the improved format of convective upwind vector flux splitting format(AUSM+), the shear stress transport turbulence model SST kω and the H2─CO reaction kinetic mechanism of 12‑step reaction of 8 components. The numerical simulation of tail flow field was performed, the change rule of pressure in base bleed combustion chamber in the transient depressurization process with time was obtained. The simulated results were basically in agreement with the experimental ones in the literature. Base this basis, the tail flow field under the condition of with and without post‑combustion in the transient depressurization process were numerically predicted, and the change rule of characteristic parameters for the tail flow field of base bleed equipment with time in the transient depressurization process were analyzed. The results show that the post‑combustion has a little effect on the tail flow field in the early stage of transient depressurization process, while in the middle and late stage, the post‑combustion has a great effect on the tail flow field. At this time, the characteristic parameters of tail flow field more quickly tend to be stable, and the post‑combustion makes the temperature of tail flow field increase and the energy addition effect remarkable, the bottom pressure increase, the bottom drag of base bleed projectile decrease significantly and the drag reduction rate increases by 75% compared with that without post‑combustion.

    Graphic Abstract

    图文摘要

    html/hncl/CJEM2018275/media/af25a83d-396b-4531-9e9d-ccb8c93183aa-image009.png

    To research the base drag and energy characteristics of a base bleed projectile under transient depressurization, average base pressure was used to show the difference of with and without post‑combustion. Then, the relative drag reduction rate was calculated to describe the effect of post‑combustion on drag reduction.

  • 1 引 言

    1

    底排弹在离开膛口时会经历一个快速降压过程,燃烧室内压力急剧降低,这种强非稳态过程对底排装置排气造成强烈扰动,从而影响泄压过程中的减阻效率。底排药剂燃烧产生的高温负氧型气体与空气中的氧气接触后发生二次燃烧,二次燃烧释放的能量使减阻率显著提[1]

    目前,在底排弹瞬态泄压方面,国内外已有大量研究。实验方面,Jackson[2]借助超声波技术实验测量了初始压力为16.3 MPa的瞬态泄压条件下,高氯酸铵(AP)/端羟基聚丁二烯(HTPB)复合推进剂的燃速变化;Yu[3]和陆春义[4,5,6]采用半密闭爆发器模拟炮口工况,借助高速录像系统研究了20~90 MPa下高速降压下底排点火具的燃烧特性;赵宏立[7]实验研究了一种双基发射药在瞬态降压下的燃烧行为,给出了不同初始压力段的双基发射药熄火临界特性曲线;Ye[8]实验研究了快速降压条件下AP/HTPB底排推进剂的瞬时燃烧特性和燃烧不稳定性;曹永杰[9]采用可视化半密闭燃烧模拟实验系统研究了3.15~7.24 MPa下底排推进剂在瞬态泄压条件下的燃烧失稳特性。数值模拟方面,曹永杰[10]建立了AP/HTPB底排推进剂燃烧流场的二维轴对称非稳态模型,数值研究了底排模拟装置在瞬态泄压工况下的流场特性;曹颖[11]运用fluent软件对底排弹出膛口过程中近膛口以及底排装置内腔的流场进行了数值模拟;Zhuo[12]采用高分辨率迎风格式AUSMPW+格式编程模拟研究了底排弹出膛口过程。在底排装置尾部流场二次燃烧方面,近些年主要以数值模拟为主。Choi[13]使用改进的两方程k‑ω SST湍流模型,模拟了含二次燃烧条件下的底排弹的尾部流场流动和燃烧过程;Xue[14]通过数值模拟对比研究了超音速下热排气的和有二次燃烧的底排装置底部流场;余文杰[15]通过编程模拟了不同排气温度下的底排装置二次燃烧尾部流场,发现排气温度在大于1100 K时二次燃烧更为充分。

    在底排减阻中,加质的作用仅占20%,而加能的作用达到80%[16]。所加的能量来自于高温燃气的热能和尾部的二次燃烧所释放的热量,而二次燃烧所释放的热量远大于排气所释放的热[17]。有关底排弹瞬态泄压过程中二次燃烧对尾部流场的加能影响目前尚未见报道。为此,本研究在与陆春[6]实验对比验证数值模型的基础上,对泄压过程中有无二次燃烧两种条件下的尾部流场进行了数值模拟,研究了二次燃烧对底排装置泄压过程中尾部流场的影响,为研究底部排气减阻技术以及底部排气弹工程设计提供参考依据

  • 2 控制方程

    2

    假设尾部流场轴向对称,湍流模型选用SST[18]模型,则守恒形式的二维轴对称Navier‑Stokes[19,20]方程为

    U t + F x + G y - F v x - G v y = W + Q + S
    (1)
    U = ρ , ρ u , ρ v , e , ρ k , ρ ω , ρ i
    (2)
    F = ρ u , ρ u 2 + p , ρ u v , ( e + p ) u , ρ k u , ρ ω u , ρ i u
    (3)
    G = ρ v , ρ u v , ρ v 2 + p , ( e + p ) v , ρ k v , ρ ω v , ρ i v
    (4)
    F v = 0 , τ x x , τ x y , u τ x x + v τ x y + q x , ( μ l + μ t σ k ) k x , ( μ l + μ t σ ω ) ω x , ρ D i Y i x
    (5)
    G v = 0 , τ x y , τ y y , u τ x y + v τ x x + q y , ( μ l + μ t σ k ) k y , ( μ l + μ t σ ω ) ω y , ρ D i Y i y
    (6)
    W = 0,0 , 0,0 , S k , S w , 0
    (7)
    Q = - v y ρ , ρ u , ρ v , e + p , ρ k , ρ ω , ρ i
    (8)
    S = 0,0 , 0,0 , 0,0 , σ i
    (9)

    式中, U 为守恒向量;F,G为对流矢通量; F v , G v 为粘性矢通量;Q为轴对称源项;W为湍流源项;S为化学反应源项; x , y 分别代表轴向和径向长度,m; u v 分别为轴向速度和径向速度,m·s-1 ω 为湍动能耗散率,m2·s-3 k 为湍动能,m2·s-2 ρ 为密度,kg·m-3 p 为压力,N·m-2 τ 为粘性应力,N·m-2 ρ i Y i D i h i 分别为 i 组分的密度(kg·m-3)、质量分数、扩散系数(m2·s)、焓(kJ·kg); λ 为热传导系数(W·(m·K)-1); μ l μ t 分别为层流黏性系数和湍流黏性系数,kg·(m·s)-1 q x q y 分别为轴向和径向的导热热流,W·m-2;形式为:

    q x = λ T x + ρ i = 1 i D i h i Y i x q y = λ T y + ρ i = 1 i D i h i Y i y
    (10)

    e 为单位体积总能,kJ·m-3;形式为:

    e = i = 1 i ρ i h i + 1 2 ρ ( u 2 + v 2 ) - p
    (11)

    σ i i 组分的化学反应质量生产率。湍流源项中的 S k S w 的给定参见文献[18]。不含二次燃烧条件下的模型去掉了化学反应源项 s ,排出气体为空气。

  • 3 模拟方法

    3
  • 3.1 计算模型

    3.1

    底排模型的径向截面如图1所示:模型为圆柱体结构(有船尾),x、r分别表示轴向和径向;R为模型最大半径;船尾长和模型最大半径相同,船尾角为5°;Tj、rj分别为排气温度、排气口半径;L为燃烧室长度,Ls为排气口处固壁厚度,vj为排气速度。具体参数值在表1中给出,Ma、p、T、P0、T0分别为来流马赫数、来流静压、来流静温、燃烧室初始压力、燃烧室初始温度。

    图1
                            瞬态泄压底排装置模型

    图1 瞬态泄压底排装置模型

    Fig.1 Model of base bleed equipment under transient depressurization

    表1 瞬态泄压数值模拟参数

    Table 1 Parameters of numerical simulation under transient depressurization

    parametervalueparametervalue
    Ma2.0p/Pa101325
    R/m0.05rj/R0.375
    Ls/m0.005p0/MPa36
    T/K293L/m0.07
    T0/K1812
    表1
                    瞬态泄压数值模拟参数

    NOTE: Ma is the free stream Maher number, Pis the free stream pressure, Tis the free stream temperature, R is the radius of axisymmetric body, rj is the radius of circular jet, L is the length of combustion chamber, Ls is the wall thickness of jet, P0 is the initial pressure of combustion chamber, T0 is the initial temperature of combustion chamber.

    底排模型尾部区域网格采用弧长法生[21]。壁面第一层网格处y+控制在2以内,网格数目35000,远场采用基于一维Riemann不变量的无反射边界条件,固壁采用无滑移、绝热、压力梯度为零的边界条件,中心轴线上采用对称边界条件。燃烧室内初始时刻没有轴向流动,燃面处径向速度vj给定参见文献[22]。

    图2
                            底排装置内部及尾部区域网格图

    图2 底排装置内部及尾部区域网格图

    Fig.2 Grid diagram of the inner and tail regions of base bleed equipment

  • 3.2 燃烧模型

    3.2

    含二次燃烧情况下底排药剂燃面燃气各组分的质量分数参照文献[22],有无二次燃烧两种条件都考虑降压速率对底排推进剂燃速的影响,采用线性燃速修正公[23]

    r ˙ = r ˙ 0 1 + φ n α r ˙ 0 2 p d p d t
    (12)

    式中, r ˙ 0 = b p n bn分别为燃速系数和压力指数; φ 为瞬态燃速修正系数; α 为底排药剂热扩散系数。

    根据Gibeling[22]等建立的8组分12步基元反应的H2─CO反应动力学机理以及有限速率化学模型,模拟底排装置尾部流场的二次燃烧过程。

  • 3.3 数值计算格式

    3.3

    根据刘君[24]的统一算法的思路,采用有限体积法编程求解方程组多组分Navier‑Stokes方程。对流项通过结合改进的AUSM+格式和MUSCL格式,并加入Van Leer限制器离[25],从而提高计算精度。黏性项采用局部坐标变换法来消除中心差分格式中奇偶失连现象。化学反应源项采用二阶隐式梯形公式求解,以消除非平衡流动的刚性问[24]。时间项采用双时间步LU‑SGS方[26]求解。求出守恒变量之后,再通过牛顿迭代法求解温度,压强可由分压定理求出。

  • 4 数值模拟结果与分析

    4
  • 4.1 数值验证

    4.1

    选择陆春义[6]的底排瞬态泄压实验结果与数值计算结果对比。底排瞬态泄压实验模型为圆柱形结构,燃烧室自由容积为550 ml,喷口直径为14 mm,破膜压力为56 MPa。燃烧室压力随时间变化的实验结果如图3所示。针对该实验模型及条件,利用本文数值方法进行模拟,得到底排燃烧室内压力随时间变化曲线,并与实验结果进行对比,结果如图3所示。由图3可见,燃烧室内压力随时间逐渐减小,模拟结果和实验结[6]吻合较好。

    图3
                            燃烧室内压力随时间变化曲线

    图3 燃烧室内压力随时间变化曲线

    Fig.3 Change curve of pressure in combustion chamber with time

  • 4.2 二次燃烧对尾部压力的影响

    4.2

    4为有无二次燃烧两种条件下不同时刻底排装置尾部区域的压力分布云图。单幅云图上半部分为含二次燃烧情况,下半部分为不含二次燃烧情况。由图4a和图4b可见,0.1 ms≤ t ≤0.3 ms时,含二次燃烧和不含二次燃烧两种情况的结果相似,高压燃气喷出底排装置之后都迅速膨胀,在0.3ms时排气口下游均出现明显的低压区。由图4b~图4d所见,0.3 ms< t ≤1.0 ms时,两种情况下燃烧室内压力都逐渐降低。0.7 ms时,排气口附近压力降至来流压力值,底排装置底部低压区都逐渐减小,燃气经过低压区后压力升高,这是因为入射激波在中心轴线上相遇反射,形成反射倾斜激波。1.0 ms时,含二次燃烧的情况下,底排燃烧室内压力降至最低,压力值约为0.75 p,尾部低压区消失,而不含二次燃烧的情况下,尾部压力略有升高,但仅达到0.45 p左右,且低压区有了扩大的趋势。由图4d~图4f所见,1.0 ms< t ≤3.0 ms时,含二次燃烧的情况下,燃烧室内压力先逐渐回升至0.9 p,之后基本不变,尾部压力场逐渐成型。底部压力先升高再略微降低,下游再压缩激波略为增强,并形成一个高压区。不含二次燃烧的情况下,燃烧室压力逐渐回升至0.74 p,之后基本不变,下游再压缩激波显著增强。总的来说,两种情况在泄压前期,压力场变化相似,泄压中后期差距明显,含二次燃烧的条件下底部压力高于不含二次燃烧的情况,减阻效果更为明显。

    html/hncl/CJEM2018275/media/af25a83d-396b-4531-9e9d-ccb8c93183aa-image004.png

    a. t=0.1 ms b. t=0.3 ms

    html/hncl/CJEM2018275/media/af25a83d-396b-4531-9e9d-ccb8c93183aa-image005.png

    c. t=0.7 ms d. t=1.0 ms

    html/hncl/CJEM2018275/media/af25a83d-396b-4531-9e9d-ccb8c93183aa-image006.png

    e. t=1.5 ms f. t=3.0 ms

    图4 不同时刻底排尾部压力分布云图(上半图:有二次燃烧;下半图:无二次燃烧)

    Fig.4 Pressure distribution nephograms at the tail of base bleed at different times (Upper part of the figure: with post‑combustion; lower part of the figure: without post‑combustion)

    5为有无二次燃烧两种条件下不同时刻的装置底部平均压力曲线图。将含二次燃烧情况记为PC(Post‑Combustion),不含二次燃烧记为NPC(Non Post‑Combustion)。由图5可见,含二次燃烧情况下,底压先迅速增加到峰值1.78 p,之后再迅速降低到最低值0.55 p,然后逐渐回升到0.9 p,随后变化较小,3.0 ms时缓慢降至0.82 p。不含二次燃烧情况下,0.1 ms时,底压略有升高,随后迅速降低,0.8 ms时达到最低0.32 p,然后逐渐回升到0.62 p,随后变化很小,在3.0 ms时缓慢升至0.69 p。由此可见,降压过程中,含二次燃烧的底部压力一直高于不含二次燃烧的底部压力,且含二次燃烧的底部压力更快地接近末态压力。

    6为相对减阻率随时间的变化曲线,相对减阻率的计算公式为:

    R B = C 1 - C 2 C 1 = p 2 - p 1 p - p 1
    (13)

    式中,C1C2分别为不含二次燃烧时的底阻和含二次燃烧时的底阻;p1p2p分别为不含二次燃烧时的底压、含二次燃烧时的底压、来流压力。由于0.3 ms内底部压力大于常压,泄压气体主要起推进作用,而本研究主要考虑泄压中底部排气状态,故RB从0.4 ms开始取值。

    由图6可见,泄压前期RB从1.4急剧减小到0.28,然后在1.3 ms时增加到0.75,泄压后期随后逐渐减小到0.43。在整个泄压过程中,相对减阻率都在0.28以上,泄压中后期相对减阻率最高达0.75。

  • 4.3 二次燃烧对尾部温度的影响

    4.3

    7为有无二次燃烧两种条件下不同时刻底排装置尾部区域的温度分布云图。单幅云图上半部分为有二次燃烧情况,下半部分为无二次燃烧情况。由图7a所见,0.1 ms时,两种情况结果类似。燃烧室出口处温度略有降低,底部固壁附近温度还未明显升高。由图7b所见,0.3 ms时,含二次燃烧模型底排装置底部固壁回流区燃气速度较慢,同时因为回流混入了来流空气,燃烧反应更为充分,温度明显升高,达到2400 K以上,剪切层中有部分燃气与空气直接接触发生反应使温度升高;不含二次燃烧模型装置底部固壁处因混入来流空气温度降低。两种情况底排燃烧室出口燃气温度降低,中轴线部分气体未与来流空气直接接触,温度分布相似。由图7a~图7c所见,0.1 ms≤t≤0.9 ms时,两种情况下,除了底排药剂部分,燃烧室内温度都因流速较快而降低。含二次燃烧模型的固壁处由于回流区扩大,高温区扩大,剪切层中也形成明显高温区,轴线部分温度也因此升高。不含二次燃烧的模型燃烧室内温度明显降低,底部固壁处压力略有降低,0.9 ms时因燃烧室内气体流出膨胀从而在出口附近有明显的低温区。由图7c和图7d所见,0.9 ms<t≤1.3 ms,两种条件燃烧室内流速减慢,燃烧室内温度都逐渐升高。含二次燃烧条件下排气装置底部固壁处高温区扩大至中心轴线,下游二次燃烧生成的高温区扩大。不含二次燃烧条件下排气装置底部固壁处温度略有升高,分布逐渐均匀,下游温度降低。由图7d~图7f所见,1.3 ms<t≤3.0 ms时,两种情况下燃烧室温度都逐渐升高。燃烧室内因燃气挤压作用,中心轴线处温度较高,这种现象在含二次燃烧条件下更为明显。含二次燃烧条件下底排装置底部固壁处温度略有降低,但分布逐渐均匀,燃烧室温度升高,由于此时排气质量流率基本稳定,且在回流区内反应殆尽,下游温度逐渐降低达到稳定。不含二次燃烧条件下装置底部温度逐渐降低,下游区域温度逐渐降低。总的来说,在泄压前期,有无二次燃烧对温度场的影响不明显,但随时间推移,两种条件下的温度场差别越来越大,含二次燃烧模型在装置底部出现二次燃烧区域使温度升高,加能效果明显。

    图5
                            底部平均压力随时间变化曲线

    图5 底部平均压力随时间变化曲线

    Fig.5 Change curves of averaged base pressure with time

    图6
                            相对减阻率随时间的变化曲线

    图6 相对减阻率随时间的变化曲线

    Fig.6 Change curve of relative drag reduction rate with times

    html/hncl/CJEM2018275/media/af25a83d-396b-4531-9e9d-ccb8c93183aa-image010.png

    a. t=0.1 ms b. t=0.3 ms

    html/hncl/CJEM2018275/media/af25a83d-396b-4531-9e9d-ccb8c93183aa-image011.png

    c. t=0.9 ms d. t=1.3 ms

    html/hncl/CJEM2018275/media/af25a83d-396b-4531-9e9d-ccb8c93183aa-image012.png

    e. t=1.5 ms f. t=3.0 ms

    图7 不同时刻底排尾部温度分布云图(上半图:有二次燃烧;下半图:无二次燃烧)

    Fig.7 Temperature distribution nephograms at the tail of base bleed at different times (Upper part of the figure: with post‑combustion; lower part of the figure: without post‑combustion)

    8为有无二次燃烧两种条件下不同时刻x=0.01R处径向温度分布曲线。由图8a可知,含二次燃烧条件下在泄压前期和中期,燃烧室出口拐角处(y=0.375R)由于燃气二次燃烧温度显著升高,在泄压后期,燃烧室出口流速较慢,整个底部温度分布均匀。由图8b可知,不含二次燃烧条件下燃烧室出口拐角处(y=0.375R)由于膨胀温度降低,在径向上温度先在0.3R处升到最高高后降低,图中1.5 ms时因为燃烧室内部分燃气贴着壁面流出造成0.4R处温度也较高。整体来看,含二次燃烧条件下底部x=0.01R处径向温度显著高于不含二次燃烧的条件。

    html/hncl/CJEM2018275/media/af25a83d-396b-4531-9e9d-ccb8c93183aa-image007.png

    a. contain post‑combustion

    html/hncl/CJEM2018275/media/af25a83d-396b-4531-9e9d-ccb8c93183aa-image008.png

    b. not contain post‑combustion

    图8 不同时刻x=0.01R处径向温度分布曲线

    Fig.8 Radial temperature distribution curves at different times in x=0.01R

  • 4.4 二次燃烧对尾部马赫数的影响

    4.4

    9有无二次燃烧两种条件下不同时刻底排装置尾部区域的马赫数分布云图。单幅云图上半部分为有二次燃烧情况,下半部分为无二次燃烧情况。由图9a~图9c所见,0.1 ms≤t≤0.7 ms时,底排燃气出排气口后迅速向外膨胀,速度逐渐增大,峰值都达到6 Ma左右。随后由于燃烧室内压力逐渐减小,燃气出排气口后膨胀程度逐渐减小,装置尾部区域的速度峰值逐渐减小,在0.7 ms时两者排气口处的马赫数都维持在当地声速。由图9c~图9e所见,0.7 ms<t≤1.3 ms时,随着燃烧室内压力逐渐降低,1.0 ms时含二次燃烧条件下排气口边缘处有少量气体回流到燃烧室中,1.3 ms时回流消失。底排装置尾部的主回流区逐渐增大,回流区数量先在1.0 ms时增加到四个,然后在1.3 ms时减少到一个。而不含二次燃烧条件下底部主回流区逐渐增大,但燃烧室内压力降低较慢,在1.3 ms时才有气体回流。由于气体从燃烧室出口冲出生成膨胀波,在出口中心轴线处存在高速区,随着燃烧室内压力降低而缩小。由图9e~图9g所见,1.3 ms<t≤3.0 ms时,含二次燃烧条件下马赫数分布基本成型。不含二次燃烧条件下燃烧室内气体回流现象消失,马赫数分布基本成型。总的来说,两种情况下的速度场在泄压前期和泄压后期分布相似,而在泄压中期,不含二次燃烧的速度场变化滞后于含二次燃烧的速度场。

    10给出了有无二次燃烧条件下不同时刻x=0.01R处径向马赫数分布。由图可知含二次燃烧条件下在1.0 ms时速度场就趋于稳定,而不含二次燃烧条件下在1.3 ms速度场才趋于稳定,且在泄压后期,含二次燃烧条件下的径向速度波动更为明显。

    html/hncl/CJEM2018275/media/af25a83d-396b-4531-9e9d-ccb8c93183aa-image015.png

    a. t=0.1 ms b. t=0.3 ms

    html/hncl/CJEM2018275/media/af25a83d-396b-4531-9e9d-ccb8c93183aa-image016.png

    c. t=0.7 ms d. t=1.0 ms

    html/hncl/CJEM2018275/media/af25a83d-396b-4531-9e9d-ccb8c93183aa-image017.png

    e. t=1.3 ms f. t=1.5 ms

    图9 不同时刻底排尾部马赫数分布云图(上半图:有二次燃烧;下半图:无二次燃烧)

    Fig.9 Mach number distribution nephograms at the tail of base bleed at different times (Upper part of the figure: with post‑combustion; lower part of the figure: without post‑combustion)

    html/hncl/CJEM2018275/media/af25a83d-396b-4531-9e9d-ccb8c93183aa-image018.png

    a. contain post‑combustion

    html/hncl/CJEM2018275/media/af25a83d-396b-4531-9e9d-ccb8c93183aa-image019.png

    b. not contain post‑combustion

    图10 不同时刻x=0.01R处径向马赫数分布曲线

    Fig.10 Radial Mach number distribution curves at different times in x=0.01R

  • 5 结 论

    5

    (1) 在底排装置泄压过程中,二次燃烧对燃烧室和尾部流场的影响很大,且主要体现在泄压中后期后期(t≥0.3 ms)。1.0 ms时,含二次燃烧条件下的燃烧室内外压力分布较为均匀,在x=0.01R处径向马赫数分布比较稳定。而此时不含二次燃烧条件下燃烧室内外压力相差较大,在x=0.01R处,径向马赫数依然受到燃烧室出口的影响。说明含二次燃烧条件下泄压过程流场更快趋于稳定。

    (2) 二次燃烧显著改变了泄压过程中温度场分布,底排装置底部温度明显升高,泄压前期与中期,流场下游部分二次燃烧充分,全过程流场整体温度高于不含二次燃烧的情况,加能效果明显。

    (3) 二次燃烧使底排装置底部压力明显升高,压力分布更加均匀,在泄压中后期其减阻率相比不含二次燃烧的减阻率提高了75%。

  • 参考文献

    • 1

      丁则胜,邱光纯,刘亚飞, 等. 固体燃料底部排气空气动力研究[J]. 空气动力学学报, 1991, 9(3): 300-307

      DING Ze‑sheng, QIU Guang‑chun, LIU Ya‑fei, et al. An aerodynamic investigation of base bleed by solid fule[J]. Acta Aerodynamic Sinica, 1991, 9(3): 300-307

    • 2

      Jackson D E,Linn C A B,Moser M D. Ultrasonic burn rate testing of bi-plateau propellants under decreasing pressure conditions[C]//41st AIAA/ASME/SAE/ASEE Joint Propulsion Conference & Exhibit. Tucson,AZ,USA,2005,AIAA2005-4164.

    • 3

      Yu Y, Zhou Y, Lu C, et al. Study on unsteady combustion behaviors of AP/HTPB base bleed propellants under transient depressurization conditions[J]. Propellants ,Explosives, Pyrotechnics, 2014, 39(4): 511-517.

    • 4

      陆春义, 周彦煌, 余永刚. 高降压速率下复合底排药剂瞬变燃烧特性研究[J]. 含能材料, 2007, 15(6):587-591.

      LU Chun‑yi, ZHOU Yan‑huang,YU Yong‑gang. Combustion of composite base bleed charge under rapid depressurization [J]. Chinese Journal of Energetic Materials(Hanneng Cailiao),2007, 15(6) : 587-591.

    • 5

      陆春义,周彦煌,余永刚.底排点火具在高降压速率下瞬态燃烧特性的实验研究[J].含能材料, 2008, 16(5): 629-632.

      LU Chun‑yi, ZHOU Yan‑huang,YU Yong‑gang. Transient combustion performance of igniter for base bleed unit under rapid depressurization[J]. Chinese Journal of Energetic Materials(Hanneng Cailiao) , 2008, 16(5): 629-632.

    • 6

      陆春义. 底排装置强非稳态燃烧特性研究[D]. 南京:南京理工大学,2009.

      LU Chun‑yi. Study on characteristics of unstable combustion of base bleed unit[D]. Nanjing: Nanjing University of Science and Technology,2009.

    • 7

      赵宏立, 刘来东, 靳建伟, 等. 压力下降条件下一种双基发射药的瞬态燃烧特性[J]. 火炸药学报, 2012, 35(1): 64-68.

      ZHAO Hong‑li, LIU Lai‑dong, JIN Jian‑wei, et al. Transient combustion performance of gun propellant under depressurization condition[J]. Chinese Journal of Explosive & Propellants, 2012, 35(1): 64-68.

    • 8

      Ye R, Yu Y, Cao Y. Experimental study of transient combustion characteristics of AP/HTPB base bleed propellant under rapid pressure drop[J]. Combustion Science & Technology, 2015, 187(3):445-457.

    • 9

      曹永杰, 余永刚, 张明安, 等. 底排推进剂瞬态泄压工况下的燃烧失稳特性[J]. 燃烧科学与技术, 2015(05): 421-427.

      CAO Yong‑jie, YU Yong‑gang, ZHANG Ming‑an,et al. Combustion instability characteristics of base bleed propellant under transient depressurization[J]. Journal Combustion Science and Technology, 2015(05): 421-427.

    • 10

      曹永杰, 余永刚, 叶锐, 等, 底排推进剂瞬态泄压工况下燃烧流场特性的数值模拟[J]. 含能材料, 2013(4): 464-468.

      CAO Yong‑jie, YU Yong‑gang, YE Rui, et al. Numerical simulation of combustion flow field characteristics of base bleed propellant under transient pressure‑release[J]. Chinese Journal of Energetic Materials(Hanneng Cailiao), 2013(4): 464-468.

    • 11

      曹颖, 王栋, 武晓松. 底排装置瞬态泄压过程数值模拟[J]. 四川兵工学报, 2014, 35(5): 43-46.

      CAO Ying, WANG Dong, WU Xiao‑song. Numerical simulation of transient depressurization process of base bleed[J]. Sichuan Acta Armamentarii, 2014, 35(5): 43-46.

    • 12

      Zhuo C, Feng F, Wu X, et al. Numerical simulation of the muzzle flows with base bleed projectile based on dynamic overlapped grids[J]. Computers & Fluids, 2014, 105:307-320.

    • 13

      Choi J Y , Shin E, Kim C K. Numerical study of Base-Bleed projectile with external combustion[R]. AIAA 2005-4352, 2005.

    • 14

      Xue X, Yu Y. An improvement of the base bleed unit on base drag reduction and heat energy addition as well as mass addition[J]. Applied Thermal Engineering, 2016, 109: 238-250.

    • 15

      余文杰, 余永刚. 底排装置尾部化学非平衡流的数值模拟[J]. 推进技术, 2015, 36(11): 1610-1615.

      YU Wen‑Jie , YU Yong‑gang. Numerical simulation of base flow with chemical non‑equilibrium for base bleed equipment[J]. Journal of Propulsion Technology, 2015, 36(11): 1610-1615.

    • 16

      丁则胜, 罗荣, 陈少松, 等. 底部燃烧减阻性能的若干参数影响研究[J]. 弹道学报, 1996, 8(4): 79-83.

      DING Ze‑sheng, LUO Rong, CHEN Shao‑song, et al. A study of some parameters influence on performance of drag reduction by base burning[J]. Journal of Ballistics, 1996, 8(4):79-83.

    • 17

      余文杰, 余永刚. 二次燃烧对底排装置尾部流场影响的数值模拟[J]. 爆炸与冲击, 2015, 35(1): 94-100.

      YU Wen‑Jie, YU Yong‑Gang. Numerical simulation of secondary combustion affecting base flow of base bleed equipment[J]. Explosion & Shock Waves, 2015, 35(1): 94-100.

    • 18

      Menter F R. Two‑equation eddy‑viscosity turbulence models for engineering applications[J]. AIAA Journal, 1994, 32(8):1598-1605.

    • 19

      朱自强. 应用计算流体力学[M]. 北京:北京航空航天大学出版社, 1998: 43-44.

      ZHU Zi‑qiang. Applied Computational Fluid Dynamics[M]. Beijing University of Aeronautics and Astronautics Press, 1998: 43-44.

    • 20

      吴子牛. 计算流体力学基本原理[M]. 北京:科学出版社, 2001:48-50.

      WU Zi‑Niu. Basic Principles of Computational Fluid Dynamics [M]. Bejing:Science China Press, 2001: 48-50.

    • 21

      武频, 赵润祥, 郭锡福. 弧长网格生成法及其应用[J]. 南京理工大学学报, 2002, 26(5): 482-485.

      WU Pin, ZHAO Run‑xiang, GUO Xi‑fu. Arc length method of grid generation and its application[J]. Journal of Nanjing University of Science and Technology, 2002, 26(5): 482-485.

    • 22

      Gibeling H J, Buggeln R C. Reacting flow models for navier‑stokes analysis of projectile Base Combustion[R]. AIAA 91-2077, 1991.

    • 23

      Kuo K K, Summerfield M. Fundamentals of solid‑propellant combustion[M]. Beijing: China Astronautic Publishing House, 1994.

    • 24

      刘君,张涵信,高树椿. 一种新型的计算化学非平衡流动的解耦方法[J]. 国防科技大学学报,2000,22(5):19-22

      LIU Jun, ZHANG Xin‑han, GAO Shu‑chun. A new uncoupled method for numerical simulation of nonequilibrium flow[J]. Journal of National University of Defense Technology, 2000, 22(5): 19-22.

    • 25

      梁德旺, 王可. AUSM+格式的改进[J]. 空气动力学学报, 2004, 22(4): 404-409.

      LIANG De‑Wang, WANG Ke. Improvement of AUSM+ scheme[J]. Acta AerodynamicSinica, 2004, 22(4): 404-409.

    • 26

      Pulliam T H. Time Accuracy and the use of Implicit Methods[R]. AIAA 93-3360, 1993.

周书培

机 构:南京理工大学 能源与动力工程学院,江苏 南京 210094

Affiliation:School of Energy and Power Engineering, Nanjing University of Science & Technology, Nanjing 210094, China

邮 箱:zhou_shupei@163.com

作者简介:周书培(1993-),男,博士研究生,主要从事底排装置尾部流场数值模拟的研究。e‑mail:zhou_shupei@163.com

余永刚

机 构:南京理工大学 能源与动力工程学院,江苏 南京 210094

Affiliation:School of Energy and Power Engineering, Nanjing University of Science & Technology, Nanjing 210094, China

角 色:通讯作者

Role:Corresponding author

邮 箱:yygnjust801@163.com

作者简介:余永刚(1963-),男,教授,主要从事含能材料燃烧推进理论与技术的研究。e‑mail:yygnjust801@163.com

html/hncl/CJEM2018275/media/af25a83d-396b-4531-9e9d-ccb8c93183aa-image001.png
parametervalueparametervalue
Ma2.0p/Pa101325
R/m0.05rj/R0.375
Ls/m0.005p0/MPa36
T/K293L/m0.07
T0/K1812
html/hncl/CJEM2018275/media/af25a83d-396b-4531-9e9d-ccb8c93183aa-image002.png
html/hncl/CJEM2018275/media/af25a83d-396b-4531-9e9d-ccb8c93183aa-image003.png
html/hncl/CJEM2018275/media/af25a83d-396b-4531-9e9d-ccb8c93183aa-image004.png
html/hncl/CJEM2018275/media/af25a83d-396b-4531-9e9d-ccb8c93183aa-image005.png
html/hncl/CJEM2018275/media/af25a83d-396b-4531-9e9d-ccb8c93183aa-image006.png
html/hncl/CJEM2018275/media/af25a83d-396b-4531-9e9d-ccb8c93183aa-image013.png
html/hncl/CJEM2018275/media/af25a83d-396b-4531-9e9d-ccb8c93183aa-image014.png
html/hncl/CJEM2018275/media/af25a83d-396b-4531-9e9d-ccb8c93183aa-image010.png
html/hncl/CJEM2018275/media/af25a83d-396b-4531-9e9d-ccb8c93183aa-image011.png
html/hncl/CJEM2018275/media/af25a83d-396b-4531-9e9d-ccb8c93183aa-image012.png
html/hncl/CJEM2018275/media/af25a83d-396b-4531-9e9d-ccb8c93183aa-image007.png
html/hncl/CJEM2018275/media/af25a83d-396b-4531-9e9d-ccb8c93183aa-image008.png
html/hncl/CJEM2018275/media/af25a83d-396b-4531-9e9d-ccb8c93183aa-image015.png
html/hncl/CJEM2018275/media/af25a83d-396b-4531-9e9d-ccb8c93183aa-image016.png
html/hncl/CJEM2018275/media/af25a83d-396b-4531-9e9d-ccb8c93183aa-image017.png
html/hncl/CJEM2018275/media/af25a83d-396b-4531-9e9d-ccb8c93183aa-image018.png
html/hncl/CJEM2018275/media/af25a83d-396b-4531-9e9d-ccb8c93183aa-image019.png

图1 瞬态泄压底排装置模型

Fig.1 Model of base bleed equipment under transient depressurization

表1 瞬态泄压数值模拟参数

Table 1 Parameters of numerical simulation under transient depressurization

图2 底排装置内部及尾部区域网格图

Fig.2 Grid diagram of the inner and tail regions of base bleed equipment

图3 燃烧室内压力随时间变化曲线

Fig.3 Change curve of pressure in combustion chamber with time

图4 不同时刻底排尾部压力分布云图(上半图:有二次燃烧;下半图:无二次燃烧) -- a. t=0.1 ms b. t=0.3 ms

Fig.4 Pressure distribution nephograms at the tail of base bleed at different times (Upper part of the figure: with post‑combustion; lower part of the figure: without post‑combustion) -- a. t=0.1 ms b. t=0.3 ms

图4 不同时刻底排尾部压力分布云图(上半图:有二次燃烧;下半图:无二次燃烧) -- c. t=0.7 ms d. t=1.0 ms

Fig.4 Pressure distribution nephograms at the tail of base bleed at different times (Upper part of the figure: with post‑combustion; lower part of the figure: without post‑combustion) -- c. t=0.7 ms d. t=1.0 ms

图4 不同时刻底排尾部压力分布云图(上半图:有二次燃烧;下半图:无二次燃烧) -- e. t=1.5 ms f. t=3.0 ms

Fig.4 Pressure distribution nephograms at the tail of base bleed at different times (Upper part of the figure: with post‑combustion; lower part of the figure: without post‑combustion) -- e. t=1.5 ms f. t=3.0 ms

图5 底部平均压力随时间变化曲线

Fig.5 Change curves of averaged base pressure with time

图6 相对减阻率随时间的变化曲线

Fig.6 Change curve of relative drag reduction rate with times

图7 不同时刻底排尾部温度分布云图(上半图:有二次燃烧;下半图:无二次燃烧) -- a. t=0.1 ms b. t=0.3 ms

Fig.7 Temperature distribution nephograms at the tail of base bleed at different times (Upper part of the figure: with post‑combustion; lower part of the figure: without post‑combustion) -- a. t=0.1 ms b. t=0.3 ms

图7 不同时刻底排尾部温度分布云图(上半图:有二次燃烧;下半图:无二次燃烧) -- c. t=0.9 ms d. t=1.3 ms

Fig.7 Temperature distribution nephograms at the tail of base bleed at different times (Upper part of the figure: with post‑combustion; lower part of the figure: without post‑combustion) -- c. t=0.9 ms d. t=1.3 ms

图7 不同时刻底排尾部温度分布云图(上半图:有二次燃烧;下半图:无二次燃烧) -- e. t=1.5 ms f. t=3.0 ms

Fig.7 Temperature distribution nephograms at the tail of base bleed at different times (Upper part of the figure: with post‑combustion; lower part of the figure: without post‑combustion) -- e. t=1.5 ms f. t=3.0 ms

图8 不同时刻x=0.01R处径向温度分布曲线 -- a. contain post‑combustion

Fig.8 Radial temperature distribution curves at different times in x=0.01R -- a. contain post‑combustion

图8 不同时刻x=0.01R处径向温度分布曲线 -- b. not contain post‑combustion

Fig.8 Radial temperature distribution curves at different times in x=0.01R -- b. not contain post‑combustion

图9 不同时刻底排尾部马赫数分布云图(上半图:有二次燃烧;下半图:无二次燃烧) -- a. t=0.1 ms b. t=0.3 ms

Fig.9 Mach number distribution nephograms at the tail of base bleed at different times (Upper part of the figure: with post‑combustion; lower part of the figure: without post‑combustion) -- a. t=0.1 ms b. t=0.3 ms

图9 不同时刻底排尾部马赫数分布云图(上半图:有二次燃烧;下半图:无二次燃烧) -- c. t=0.7 ms d. t=1.0 ms

Fig.9 Mach number distribution nephograms at the tail of base bleed at different times (Upper part of the figure: with post‑combustion; lower part of the figure: without post‑combustion) -- c. t=0.7 ms d. t=1.0 ms

图9 不同时刻底排尾部马赫数分布云图(上半图:有二次燃烧;下半图:无二次燃烧) -- e. t=1.3 ms f. t=1.5 ms

Fig.9 Mach number distribution nephograms at the tail of base bleed at different times (Upper part of the figure: with post‑combustion; lower part of the figure: without post‑combustion) -- e. t=1.3 ms f. t=1.5 ms

图10 不同时刻x=0.01R处径向马赫数分布曲线 -- a. contain post‑combustion

Fig.10 Radial Mach number distribution curves at different times in x=0.01R -- a. contain post‑combustion

图10 不同时刻x=0.01R处径向马赫数分布曲线 -- b. not contain post‑combustion

Fig.10 Radial Mach number distribution curves at different times in x=0.01R -- b. not contain post‑combustion

image /

无注解

Ma is the free stream Maher number, Pis the free stream pressure, Tis the free stream temperature, R is the radius of axisymmetric body, rj is the radius of circular jet, L is the length of combustion chamber, Ls is the wall thickness of jet, P0 is the initial pressure of combustion chamber, T0 is the initial temperature of combustion chamber.

无注解

无注解

无注解

无注解

无注解

无注解

无注解

无注解

无注解

无注解

无注解

无注解

无注解

无注解

无注解

无注解

无注解

  • 参考文献

    • 1

      丁则胜,邱光纯,刘亚飞, 等. 固体燃料底部排气空气动力研究[J]. 空气动力学学报, 1991, 9(3): 300-307

      DING Ze‑sheng, QIU Guang‑chun, LIU Ya‑fei, et al. An aerodynamic investigation of base bleed by solid fule[J]. Acta Aerodynamic Sinica, 1991, 9(3): 300-307

    • 2

      Jackson D E,Linn C A B,Moser M D. Ultrasonic burn rate testing of bi-plateau propellants under decreasing pressure conditions[C]//41st AIAA/ASME/SAE/ASEE Joint Propulsion Conference & Exhibit. Tucson,AZ,USA,2005,AIAA2005-4164.

    • 3

      Yu Y, Zhou Y, Lu C, et al. Study on unsteady combustion behaviors of AP/HTPB base bleed propellants under transient depressurization conditions[J]. Propellants ,Explosives, Pyrotechnics, 2014, 39(4): 511-517.

    • 4

      陆春义, 周彦煌, 余永刚. 高降压速率下复合底排药剂瞬变燃烧特性研究[J]. 含能材料, 2007, 15(6):587-591.

      LU Chun‑yi, ZHOU Yan‑huang,YU Yong‑gang. Combustion of composite base bleed charge under rapid depressurization [J]. Chinese Journal of Energetic Materials(Hanneng Cailiao),2007, 15(6) : 587-591.

    • 5

      陆春义,周彦煌,余永刚.底排点火具在高降压速率下瞬态燃烧特性的实验研究[J].含能材料, 2008, 16(5): 629-632.

      LU Chun‑yi, ZHOU Yan‑huang,YU Yong‑gang. Transient combustion performance of igniter for base bleed unit under rapid depressurization[J]. Chinese Journal of Energetic Materials(Hanneng Cailiao) , 2008, 16(5): 629-632.

    • 6

      陆春义. 底排装置强非稳态燃烧特性研究[D]. 南京:南京理工大学,2009.

      LU Chun‑yi. Study on characteristics of unstable combustion of base bleed unit[D]. Nanjing: Nanjing University of Science and Technology,2009.

    • 7

      赵宏立, 刘来东, 靳建伟, 等. 压力下降条件下一种双基发射药的瞬态燃烧特性[J]. 火炸药学报, 2012, 35(1): 64-68.

      ZHAO Hong‑li, LIU Lai‑dong, JIN Jian‑wei, et al. Transient combustion performance of gun propellant under depressurization condition[J]. Chinese Journal of Explosive & Propellants, 2012, 35(1): 64-68.

    • 8

      Ye R, Yu Y, Cao Y. Experimental study of transient combustion characteristics of AP/HTPB base bleed propellant under rapid pressure drop[J]. Combustion Science & Technology, 2015, 187(3):445-457.

    • 9

      曹永杰, 余永刚, 张明安, 等. 底排推进剂瞬态泄压工况下的燃烧失稳特性[J]. 燃烧科学与技术, 2015(05): 421-427.

      CAO Yong‑jie, YU Yong‑gang, ZHANG Ming‑an,et al. Combustion instability characteristics of base bleed propellant under transient depressurization[J]. Journal Combustion Science and Technology, 2015(05): 421-427.

    • 10

      曹永杰, 余永刚, 叶锐, 等, 底排推进剂瞬态泄压工况下燃烧流场特性的数值模拟[J]. 含能材料, 2013(4): 464-468.

      CAO Yong‑jie, YU Yong‑gang, YE Rui, et al. Numerical simulation of combustion flow field characteristics of base bleed propellant under transient pressure‑release[J]. Chinese Journal of Energetic Materials(Hanneng Cailiao), 2013(4): 464-468.

    • 11

      曹颖, 王栋, 武晓松. 底排装置瞬态泄压过程数值模拟[J]. 四川兵工学报, 2014, 35(5): 43-46.

      CAO Ying, WANG Dong, WU Xiao‑song. Numerical simulation of transient depressurization process of base bleed[J]. Sichuan Acta Armamentarii, 2014, 35(5): 43-46.

    • 12

      Zhuo C, Feng F, Wu X, et al. Numerical simulation of the muzzle flows with base bleed projectile based on dynamic overlapped grids[J]. Computers & Fluids, 2014, 105:307-320.

    • 13

      Choi J Y , Shin E, Kim C K. Numerical study of Base-Bleed projectile with external combustion[R]. AIAA 2005-4352, 2005.

    • 14

      Xue X, Yu Y. An improvement of the base bleed unit on base drag reduction and heat energy addition as well as mass addition[J]. Applied Thermal Engineering, 2016, 109: 238-250.

    • 15

      余文杰, 余永刚. 底排装置尾部化学非平衡流的数值模拟[J]. 推进技术, 2015, 36(11): 1610-1615.

      YU Wen‑Jie , YU Yong‑gang. Numerical simulation of base flow with chemical non‑equilibrium for base bleed equipment[J]. Journal of Propulsion Technology, 2015, 36(11): 1610-1615.

    • 16

      丁则胜, 罗荣, 陈少松, 等. 底部燃烧减阻性能的若干参数影响研究[J]. 弹道学报, 1996, 8(4): 79-83.

      DING Ze‑sheng, LUO Rong, CHEN Shao‑song, et al. A study of some parameters influence on performance of drag reduction by base burning[J]. Journal of Ballistics, 1996, 8(4):79-83.

    • 17

      余文杰, 余永刚. 二次燃烧对底排装置尾部流场影响的数值模拟[J]. 爆炸与冲击, 2015, 35(1): 94-100.

      YU Wen‑Jie, YU Yong‑Gang. Numerical simulation of secondary combustion affecting base flow of base bleed equipment[J]. Explosion & Shock Waves, 2015, 35(1): 94-100.

    • 18

      Menter F R. Two‑equation eddy‑viscosity turbulence models for engineering applications[J]. AIAA Journal, 1994, 32(8):1598-1605.

    • 19

      朱自强. 应用计算流体力学[M]. 北京:北京航空航天大学出版社, 1998: 43-44.

      ZHU Zi‑qiang. Applied Computational Fluid Dynamics[M]. Beijing University of Aeronautics and Astronautics Press, 1998: 43-44.

    • 20

      吴子牛. 计算流体力学基本原理[M]. 北京:科学出版社, 2001:48-50.

      WU Zi‑Niu. Basic Principles of Computational Fluid Dynamics [M]. Bejing:Science China Press, 2001: 48-50.

    • 21

      武频, 赵润祥, 郭锡福. 弧长网格生成法及其应用[J]. 南京理工大学学报, 2002, 26(5): 482-485.

      WU Pin, ZHAO Run‑xiang, GUO Xi‑fu. Arc length method of grid generation and its application[J]. Journal of Nanjing University of Science and Technology, 2002, 26(5): 482-485.

    • 22

      Gibeling H J, Buggeln R C. Reacting flow models for navier‑stokes analysis of projectile Base Combustion[R]. AIAA 91-2077, 1991.

    • 23

      Kuo K K, Summerfield M. Fundamentals of solid‑propellant combustion[M]. Beijing: China Astronautic Publishing House, 1994.

    • 24

      刘君,张涵信,高树椿. 一种新型的计算化学非平衡流动的解耦方法[J]. 国防科技大学学报,2000,22(5):19-22

      LIU Jun, ZHANG Xin‑han, GAO Shu‑chun. A new uncoupled method for numerical simulation of nonequilibrium flow[J]. Journal of National University of Defense Technology, 2000, 22(5): 19-22.

    • 25

      梁德旺, 王可. AUSM+格式的改进[J]. 空气动力学学报, 2004, 22(4): 404-409.

      LIANG De‑Wang, WANG Ke. Improvement of AUSM+ scheme[J]. Acta AerodynamicSinica, 2004, 22(4): 404-409.

    • 26

      Pulliam T H. Time Accuracy and the use of Implicit Methods[R]. AIAA 93-3360, 1993.