AP/HTPB推进剂微尺度燃烧特性的数值分析

马龙泽, 余永刚

(南京理工大学能源与动力工程学院, 江苏 南京 210094)

1 引 言

高氯酸铵/端羟基聚丁二烯(AP/HTPB)是一种典型的复合固体推进剂,由氧化剂AP颗粒和粘合剂HTPB混合而成,其燃烧特性和火焰结构依赖于燃烧压力和材料非均匀性质。针对AP/HTPB的燃烧特性,前人做了诸多工作,Price等[1]通过扫描电镜(SEM)观察发现,AP/HTPB在微观尺度上呈准周期性排列,结构表现为三明治分布,因此可以将AP/HTPB复合推进剂宏观燃烧简化为微观三明治单元燃烧。Chorpening[2]利用紫外线透射成像技术观测了不同压力下推进剂燃烧表面上方的火焰结构。Beckstead等[3]对AP/HTPB复合推进剂表面结构进行了大量实验观测后,提出了BDP多火焰模型,能合理地解释推进剂的诸多燃烧特征。Bilger等[4]以Pe数或雷诺数形式研究了推进剂燃烧时对流和扩散在方程和边界条件中应用的重要性。Buckmaster[5]和Jackson[6]采用总包气相反应建立了二维和三维几何模型,对推进剂燃烧的火焰结构进行了模拟,进一步支持了Bilger的研究成果。Hegab等[7]研究了AP/HTPB复合固体推进剂燃烧气相与固相间的耦合问题。Cai[8]分析了AP颗粒尺寸、环境压力和气相反应速率等参数对推进剂燃烧气相火焰结构的影响。Jackson[9]建立了AP颗粒随机分布的三维非定常模型并进行了数值模拟。叶锐和曹永杰[10-11]以Pe数为基础,讨论了不同工况下侧向速度等参数对燃烧的影响。AP/HTPB复合固体推进剂燃烧火焰可以看作是由大量的微尺度火焰单元组成,研究单个火焰单元的燃烧特性,有助于从宏观上深入了解其燃烧机理和性能。

为了进一步研究AP/HTPB微观单元的燃烧特性和火焰结构,本研究基于气相火焰与固相推进剂的耦合,利用计算流体动力学(CFD)方法,从流动燃烧角度,采用简化的两步总包化学反应动力学机理,基于BDP多火焰模型,建立了二维周期性三明治定常燃烧模型,借助FLUENT软件,数值分析了AP体积分数一定(0.75)时,不同燃烧压力下的燃烧参数分布特性以及氧化剂含量(0.7<α<0.95)对燃烧的影响。

2 物理模型

针对本文采用的AP/HTPB二维微尺度周期性三明治燃烧模型,如图1所示,作如下简化假设:

(1)燃气为不可压理想气体,气相各组分Lewis数均为1;

(2)固相热分解和气相反应均满足Arrhenius定律;

(3)固相分解仅发生在燃面层,通过质量、能量通量平衡以及温度连续性关系耦合气固两相;

(4)不考虑燃气辐射及燃面形状;

(5)以源项法描述固相热分解及气相扩散燃烧。

图1AP/HTPB微尺度燃烧模型

Fig.1Microscale combustion model of AP/HTPB

3 数学模型

根据上述物理模型,建立低雷诺数下二维定常流动燃烧基本控制方程,具体如下:

3.1 气相控制方程

(1)质量守恒方程

(1)

(2)动量守恒方程

(2)

(3)组分守恒方程

ρgu·Yi=

(3)

(4)能量守恒方程

ρgcpu·T=·(λgT)+ST

(4)

(1)式～(4)式中,u为速度矢量,m·s-1;uj和uk为速度张量，j和k指标范围是(1,2)；ρg为气相密度,kg·m-3;μg为气相粘性系数,N·s·m-2;λg为气相导热系数,W·m-1·K-1;cp为气相的定压比热容,J·kg-1·K-1;Yi为气相组分i的质量分数;Le为无量纲Lewis数,Sj为气相动量源项,N·m-3;SYi为气相组分源项,为无量纲;ST为气相能量源项,W·m-3。

3.2 固相能量方程

(5)

式中,ρc为固相密度,kg·m-3;λc为固相导热系数,W·m-1·K-1;cc为固相比热容,J·kg-1·K-1;r为固相燃速,m·s-1。其中,

(6)

式中,下标AP和B分别代表氧化剂AP和粘结剂HTPB;x为横坐标,m;α为AP体积分数。

3.3 气相化学反应动力学

采用基于BDP多火焰模型的两步总包反应机理[7]:

AP(g)(X)→Ox(Z)

(7)

C2H4(Y)+βOx(Z)→Final products(P)

(8)

式中,X为氧化剂颗粒AP分解反应生成的NH3和HClO4,Z为其分解产生的氧化性气体混合物,Y为粘合剂HTPB的热解产物,β为化学计量系数,为无量纲,其值由AP体积分数α决定,如(9)式所示:

β=ρAPα/(ρB(1-α))

(9)

化学反应速率遵循基于压力的Arrhenius定律[7]:

R1=D1pn1[X]exp(-E1/RuT)

(10)

R2=D2pn2[Y][Z]exp(-E2/RuT)

(11)

式中,D1、D2为化学反应速率常数,[X]、[Y]、[Z]为相应组分的质量分数,n1、n2为压力指数,E1、E2为活化能,Ru为通用气体常数。

3.4 边界条件

针对图1所示计算模型,计算域取周期性半宽L=90 μm,气固两相高度均为H=500 μm,左右边界为周期性边界,上下边界为远场边界,气固耦合面处组分分布为:

(12)

4 数值模拟结果及讨论

采用基于有限体积法的CFD软件FLUENT,对AP/HTPB微观单元的定常燃烧过程进行了数值模拟,模型计算所需的化学反应动力学及物性参数[12]如表1所示。

3.1.3 采取集体备课，案例设计集思广益 基于临床实际的教学案例并不是完全的拿来主义，规范的案例设计应通过教师集体备课完成，将收集的复杂临床病例资料按教材内容和岗位任务进行梳理，使案例经过集体讨论补充与修改简单化、完整化，保证其篇幅、难度与教学相适应，且知识点全面，主次分明，设计的案例问题能导出并围绕当次课程内容。另外，对教学大纲要求和案例中涉及的需要学生实践的操作技术，在集体备课中首先进行教师间的预演，进一步讨论确定需要学生提前预习的知识与技能，估算教学所需时间，评价模拟脚本是否符合逻辑及临床实际，以更有利于教师掌控教学节奏。

表1AP/HTPB的反应动力学及物性参数[12]

Table1Reaction kinetic and physical property parameters of AP/HTPB[12]

parametervalueparametervalueEAP/J·mol-19．145´104ρAP/kg·m-31950EB/J·mol-16．236´104ρB/kg·m-3920E1/J·mol-16．651´104λAP/W·m-1·K-10．405E2/J·mol-19．145´104λB/W·m-1·K-10．276Ru/J·mol-1·K-18．341cp/J·kg-1·K-11255．2

Note:EAPandEBare the decomposition activation energies,E1andE2are activation energies,Ruis the universal gas constant,ρAPandρBare densities of AP and HTPB, respectively,λAPandλBare thermal conductivities of AP and HTPB, respectively,cpis the specific heat.

4.1 火焰结构特性

AP体积分数为0.75时,不同燃烧压力下的温度分布云图如图2所示,可以看出,当0.1 MPa≤p≤6.0 MPa时,随着压力升高,气相火焰温度随之先递增再衰减,火焰结构呈现出不同燃烧特性。由式(10)和式(11)可知,气相化学反应速率与pn相关,随着压力升高,化学反应速率迅速增大。当压力小于0.4 MPa时,火焰温度较低,但分布较为均匀,这是因为AP/HTPB燃烧的火焰温度分布受气相化学反应和扩散的共同作用,此时,扩散混合占据主导位置,固相分解产物能在气相发生反应前混合充分,故而燃烧火焰总体呈预混燃烧特性。当0.4 MPa≤p<2.5 MPa时,气相化学反应速率加快,与扩散速率相当,这时燃烧火焰呈预混-扩散特征,火焰温度分布基本呈现“W”形态,HTPB火焰高于AP火焰,说明HTPB燃速比AP大,形成AP凸出于粘结剂表面,此时当压力为1.5 MPa时,火焰达到最高温度。当2.5 MPa

4.2 火焰热释放及其对燃面的热反馈

氧化剂AP体积分数为0.75时,不同压力下的容积释放速率分布如图4所示,总体来看,气相火焰热释放速率随压力升高而递增,火焰对燃面的热反馈随压力增大而加强。压力较低(<0.4 MPa,图4a)时,预混火焰热释放区域为一个整体,纵向分布较大,容积放热速率较小,在AP燃面上方各有一个较大的放热核心,

a. 0.1 MPab. 0.4 MPa

c. 1.5 MPad. 2.5 MPa

e. 4.0 MPaf. 6.0 MPa

图2不同压力下的气相火焰温度分布

Fig.2Gas phase flame temperature distribution at different pressure

a. 0.2 MPab. 0.4 MPa

c. 1.5 MPad. 3.2 MPa

图3不同压力下的火焰紫外发射图[2]

Fig.3Flame UV emission images at different pressure[2]

随着压力升高,火焰热释放区域轴向缩减变窄,面积减小,两个放热核心收缩且逐渐分离,火焰由一个整体逐渐分裂形成两个扩散火焰带,并且压力越高,火焰放热核心越贴近燃烧表面。

图5为不同压力下燃面各处的温度分布。由图5可以看出,随着压力增大,氧化剂AP燃面温度比粘结剂燃面温度低,燃面各处温度升高,但HTPB燃面最高温度变化幅度大于AP燃面处的温度增幅,燃面最高温度位置处于粘结剂一侧。由于燃面退移速率与温度成正相关,可推断出不同压力下燃面各处的燃速,如图6所示,燃速整体随着压力升高而增大,但AP燃速随压力变化增幅更大,所以压力较低时,HTPB燃速高于AP燃速,压力较高时,HTPB燃速反而低于AP燃速。图7为Price[1]扫描电镜观测到的低压和高压下的燃面轮廓,可知低压下AP相对HTPB凸起,高压下AP相对于HTPB凹陷,计算结果与该实验相吻合。

a. 0.1 MPab. 0.4 MPa

c. 1.5 MPad. 2.5 MPa

e. 4.0 MPaf. 6.0 MPa

图4不同压力下的容积放热速率分布

Fig.4Volume heat release rate distribution at different pressure

4.3 氧化剂含量对燃烧的影响

当燃烧压力不变时,以2.5 MPa为例,不同氧化剂AP体积分数下的容积放热速率云图,如图8所示。可以看出,AP体积分数为0.7～0.95时,随着AP含量增大,则粘结剂相对越窄,容积放热速率越小,气相热释放区域由两个合并连接形成一个整体,AP/HTPB交界处上方的火焰带逐渐变短消失。图9为AP高体积分数和低体积分数下的火焰温度分布云图,结合图8可知,α=0.7时,气相火焰热释放最大,此时微观单元火焰呈现出富燃特性,有两个狭长带状火焰面,最大容积放热速率为2.54×1013W·m-3,燃烧火焰最高温度为3376 K,随着氧化剂AP体积分数增大,α=0.95时,气相热释放区域为一个整体,热释放量最小,火焰单元呈贫燃特性,此时最大容积放热速率为7.53×1012W·m-3,燃烧火焰最高温度降低为2085 K。可见,随着氧化剂AP体积分数增大,容积放热速率减小,气相对固相燃面的热反馈降低,导致燃面温度衰减,燃速降低。

图5不同压力下的燃面温度分布

Fig.5Temperature distribution of the burning surface at different pressure

图6不同压力下的燃面燃速分布

Fig.6Burning rate distribution of the burning surface at different pressure

a. 1.4 MPab. 6.9 MPa

图7低压和高压下熄火推进剂的燃面[1]

Fig.7Burning surface of the propellant quenched at low and high pressure[1]

a.α=0.7b.α=0.8

c.α=0.9d.α=0.95

图8不同AP体积分数下的体积热释放速率分布

Fig.8Volume heat release rate distribution under the different volume fraction of AP

a.α=0.7b.α=0.95

图9不同AP体积分数下的气相火焰温度分布

Fig.9Gas phase flame temperature distribution under the different volume fraction of AP

5 结 论

建立了AP/HTPB微观单元的二维周期性三明治定常燃烧模型,对不同压力和不同氧化剂含量下的扩散燃烧过程进行了数值模拟,可得如下结论:

(2) 当AP体积分数为0.75,燃烧压力为0.1～6.0 MPa时,压力越大,气相燃烧热释放及其对固相的热反馈越强,燃面的温度越高,对颗粒燃烧速度起更大的推进作用,但AP燃速随压力变化更大,故低压下AP相对HTPB凸起,高压下AP相对于HTPB凹陷。

(3) 当压力不变时,即在2.5 MPa下,AP体积分数为0.7～0.95时,随着AP含量增大,粘结剂相对越窄,气相放热量变小,火焰温度降低,两个分裂的独立放热区域逐渐合并形成一个整体,且在AP/HTPB交界处上方的两个狭长带状火焰面逐渐消失。

参考文献:

[1] Price E W, Sambamurthi J K, Sigman R K, et al. Combustion of ammonium perchlorate-polymer sandwiches[J].CombustionandFlame, 1986, 63(3): 381-413.

[2] Chorpening B T, Brewster M Q. Flame structure and burning rate of ammonium perchlorate/hydroxyl-terminated polybutadiene[J].ProceedingsoftheCombustionInstitute, 2000, 28(1): 847-853.

[3] Beckstead M W, Derr R L, Price C F. A model of composite solid-propellant combustion based on multiple flames[J].JournalofAIAA, 1970, 8(12): 2200-2207.

[4] Bilger W R, Jia X.The burke-schumann diffusion flame with zero net flux boundary condition[J].CombustionScienceandTechnology, 1994, 99(4): 371-376.

[5] Buckmaster J, Jackson T L, Yao J. An elementary discussion of propellant flame geometry[J].CombustionandFlame, 1999, 117(3): 541-552.

[6] Jackson T L, Buckmaster J. Nonpremixed periodic flames supported by heterogeneous propellants[J].JournalofPropulsionandPowder, 2000, 16(3): 498-504.

[7] Hegab A, Jackson T L,Buckmaster J, et al. Nonsteady burning of periodic sandwich propellants with complete coupling between the solid and gas phases[J].CombustionandFlame, 2001, 125(1): 1055-1070.

[8] Cai W, Thakre P, Yang V. A model of AP/HTPB composite propellant combustion in rocket-motor environments[J].CombustionScienceandTechnology, 2008, 180(12): 2143-2169.

[9] Jackson T L. Modeling of heterogeneous propellant combustion: a survey[J].JournalofAIAA, 2012, 50(5): 993-1006.

[10] 叶锐, 余永刚, 曹永杰. AP/HTPB二维火焰结构和燃速数值分析[J]. 工程热物理学报, 2013, 34(3): 576-580.

YE Rui, YU Yong-gang, CAO Yong-jie. Numerical analysis of flame structure and burning rate of two dimensional AP/HTPB sandwich combustion[J].JournalofEngineeringThermophysics, 2013, 34(3): 576-580.

[11] 曹永杰, 余永刚, 叶锐. AP/HTPB复合推进剂微尺度燃烧模型及数值分析[J]. 推进技术, 2013, 34(11): 1567-1574.

CAO Yong-jie, YU Yong-gang, YE Rui. A microscale combustion model and numerical analysis of AP/HTPB composite propellant[J].JournalofPropulsionTechnology, 2013, 34(11): 1567-1574.

[12] Jackson L T, Buckmaster J. Heterogeneous propellant combustion[J].JournalofAIAA, 2002, 40(6): 1122-1130.