作者简介:刘昊(1984—),男,博士,研究员,研究领域为组合推进技术。
(Xi'an Aerospace Propulsion Institute,Xi'an 710100,China)
RBCC engine; kerosene; methane; hydrogen; engine performance
备注
作者简介:刘昊(1984—),男,博士,研究员,研究领域为组合推进技术。
为获得不同燃料对RBCC发动机性能影响,基于部件及发动机试验数据,建立了采用控制体法,考虑热完全气体效应、化学平衡流动效应、黏性损失及热损失等影响的发动机性能分析模型。完成模型发动机分别采用煤油、CH4、H2燃料性能仿真,获得了给定飞行剖面不同燃料对发动机推力、比冲性能的影响。在给定研究条件下,结果表明:H2推力、比冲最高,火箭引射、火箭冲压及冲压模态推力分别是煤油的1.40、1.28、1.14倍,比冲分别是煤油的1.40、1.43、2.73倍; CH4推力、比冲与煤油相当,火箭引射、火箭冲压及冲压模态推力分别是煤油的1.08、1.04、0.97倍,比冲分别是煤油的1.08、1.08、1.17倍。
In order to obtain the influence of different fuel on engine performance,based on the components and engine experimental data,an engine characteristic analysis model was established,which adopt the control body method and taked into account the effects of heat complete gas,chemical equilibrium flow,viscosity loss and heat loss. Model engine was used to complete the simulation with kerosene,methane,and hydrogen. The effects of different fuel types on the engine force and specific impulse under given ballistic were obtained as the same time. The results show that:H2 achieved maximum performance,the engine force on ejector mode,rocket ramjet mode,and ramjet mode is 1.40,1.28,1.14 times to kerosene respectively,and 1.40,1.43,2.73 times respectively for engine specific impulse; methane and kerosene had similar properties,the engine force on ejector mode,rocket ramjet mode,and ramjet mode is 1.08,1.04,0.97 times as much as kerosene respectively,and 1.08,1.08,1.17 times respectively for engine specific impulse.
引言
0 引言
截至目前,国内外学者针对RBCC(rocket based combined cycle)发动机已持续开展近60年的研究[1-8],针对可重复使用天地往返,先后提出ERJ、SERJ、A5、Strutjet、GTX、ISTAR等多种发动机方案[9-14],并开展了大量实验研究。然而,由于技术复杂性,RBCC发动机至今仍未实现工程应用。掌握宽飞行马赫数范围内不同工作模态发动机特性,是发动机方案设计、总体应用论证及技术指标确定的前提,因此RBCC发动机性能分析一直为国内外学者所关注。
Olds等[15-17]较早建立了RBCC发动机全模态性能分析软件SCCREAM,并应用于H2/O2单级入轨飞行器概念方案设计[18],软件对引射掺混段、喷管、主火箭及冲压燃烧室部件简单引入效率系数进行模拟,难以反映部件真实特性。Mckamey等开发了EPSURBCC软件,软件仅支持H2燃料,并且未考虑黏性影响及部件损失[19]。黄生洪等通过变步长半隐式多步龙格库塔方法求解考虑化学反应源项的一维流动方程组,建立了发动机性能预估模型[20],同SCCREAM类似,通过引入效率系数求解燃烧室及引射段出口参数,并对文献[21]给出的H2/O2发动机引射模态完成了性能仿真。吕翔等则通过采用MacCormack格式求解带有化学反应源项的一维流动方程组,发展了发动机准一维性能分析方法[22-23],模型通过引入修正因子对前体边界层效应进行修正,完成采用H2/O2推进剂的Hyperion计划RBCC发动机[24]性能预测,并将比冲预测结果与SCCREAM软件进行了比对。安佳宁建立了RBCC发动机引射模态性能计算准一维模型,模型通过采用飞行马赫数单变量多项式拟合求解进气道总压恢复系数,采用引射器理论建立了等截面混合及等压混合两种构型计算模型,并完成H2/O2模型发动机引射模态性能分析[25]。
然而,现有发动机性能分析软件存在如下问题:①进气道、燃烧室等关键部件特性多采用效率系数假设,甚至不予考虑,难以反映发动机真实特性; ②学者大多集中于H2/O2推进剂发动机研究,缺乏其他推进剂组合发动机性能数据,特别是不同推进剂组合发动机性能对比数据。
针对上述问题,本文基于部件及发动机试验数据,建立了发动机关键部件特性数学模型,开发了RBCC发动机一体化性能计算平台,并完成飞行马赫数Ma0=0~8范围分别采用煤油/O2、CH4/O2、H2/O2推进剂组合发动机性能仿真,获得了不同工作模态下推进剂组合类型对发动机推力、比冲性能的影响,为发动机总体方案论证及燃料体系确定提供参考。
1 数学模型及验证
2 计算结果及分析
2.1 火箭引射模态图4给出了火箭引射模态发动机推力、比冲曲线。3种燃料发动机推力、比冲特性一致,即随着飞行马赫数的增加,发动机推力、比冲均先略微减小后增加,并在Ma0=0.3达到最小值。相同工况下,H2推力、比冲最高,CH4次之,煤油最低。
图5给出了以煤油性能为基准计算获得的发动机相对推力及相对比冲曲线。从图中可以看出,H2推力、比冲约是煤油的1.40倍,CH4推力、比冲约是煤油的1.08倍。
2.2 火箭冲压模态图6给出了火箭冲压模态发动机推力、比冲曲线。随着飞行马赫数的增加,发动机推力、比冲均先增加后减小,并于Ma0=3.7达到最大值。相同工况下,H2推力、比冲最高,CH4次之,煤油最低。
图7给出了以煤油性能为基准计算获得的发动机相对推力、相对比冲曲线。图中可见,H2相对推力及相对比冲随着飞行马赫数增加先增加后减小,CH4相对推力及相对比冲基本保持不变。整体而言,H2推力约是煤油的1.28倍,比冲约是煤油的1.43倍; CH4推力约是煤油的1.04倍,比冲约是煤油的1.08倍。
2.3 冲压模态图8给出了冲压模态发动机推力、比冲曲线。随着飞行马赫数的增加,发动机推力、比冲先增加后减小,并于Ma0=3.7达到最大值。相同工况下,H2推力最高,煤油次之,CH4最低; H2比冲最高,CH4次之,煤油最低。
图9给出了以煤油性能为基准计算获得的相对推力、相对比冲曲线。对于H2,随着飞行马赫数的增加,相对推力及相对比冲先增加后减小,并于Ma0=6.6达到最大值,整体而言,H2相对推力约为煤油1.14倍,相对比冲约为煤油的2.73倍。对于CH4,随着飞行马赫数增加,相对推力先略增后减小再略增加,相对比冲基本保持不变,整体而言,CH4相对推力约为煤油的0.97倍,相对比冲约为煤油的1.17倍。
3 结论
基于部件及发动机试验数据,建立了RBCC发动机部件特性模型,采用控制体法,发展了考虑热完全气体效应、化学平衡流动效应、黏性损失及热损失的RBCC发动机性能分析模型。针对几何结构可调尾喷管模型发动机,完成飞行Ma0=0~8范围采用煤油/O2、CH4/O2、H2/O2推进剂组合发动机仿真,获得了不同工作模态下发动机性能,分析了燃料类型对推力、比冲的影响。在本文给定研究条件下,获得结论如下:
1)火箭引射模态,随着飞行马赫数的增加,推力、比冲均先略微减小后增加。Ma0=0~2范围内,H2推力、比冲约是煤油燃料的1.40倍,CH4推力、比冲约是煤油燃料的1.08倍。
2)火箭冲压模态,随着飞行马赫数的增加,推力、比冲均先增加后减小。Ma0=2~8范围内,H2推力约是煤油的1.28倍,比冲约是煤油的1.43倍; CH4推力约是煤油的1.04倍,比冲约是煤油的1.08倍。
3)冲压模态,随着飞行马赫数的增加,推力、比冲均先增加后减小。Ma0=2~8范围内,H2推力约是煤油的1.14倍,比冲约是煤油的2.73倍; CH4推力约是煤油的0.97倍,比冲约是煤油的1.17倍。
1.1 计算模型本文采用文献[26]中的模型RBCC发动机方案,其结构示意及特征截面定义如图1所示。发动机由二元多楔进气道、二元单侧扩张燃烧室、火箭推力室及尾喷管组成,火箭推力室采用中心支板布局方式。图1中0截面为自由来流,2截面为进气道出口(燃烧室入口),4截面为燃烧室出口(尾喷管入口),6截面为尾喷管出口。模型发动机进气道捕获面积A0=0.08 m2,进气道、燃烧室固定几何结构,进气道收缩比5.8,燃烧室扩张比2.5,尾喷管根据具体工作条件进行调节。
1.2 计算方法及验证采用控制体法完成发动机各特征截面参数计算,计算模型考虑了各部组件黏性损失及热损失,具体计算过程如下:
1.2.1 0截面参数参考文献[27]完成计算。
1.2.2 2截面参数引射模态根据式(1)给出的引射比模型完成捕获空气流量计算,火箭冲压及冲压模态根据式(2)给出的流量系数模型完成捕获空气流量计算,结合式(3)给出的进气道总压恢复系数模型完成2截面参数计算。
n=cej·Pt2(1)
φs=b0+b1·Ma0+b2·α+b3·Ma02+
b4·Ma0·α+b5·α2(2)
σ2=a0+a1·Ma0+a2·α+a3·Ma02+
a4·Ma0·α+a5·α2(3)
式中:n为引射系数,定义为一次流流量与二次流流量之比; 常系数cej与发动机具体构型、火箭布局方式、火箭工作参数等有关,可根据试验或仿真获得; Pt2为进气道出口总压; φs为进气道流量系数; σ2为进气道总压恢复系数; Ma0为飞行马赫数; α为飞行攻角; 系数a0~a5、b0~b5由试验确定。
1.2.3 4截面参数冲压模态,根据式(4)、式(5)给出的燃烧效率及燃烧室总压恢复系数模型完成4截面参数计算; 火箭引射及火箭冲压模态,采用冲压模态与火箭冲压模态燃烧室出口马赫数相似假设完成4截面参数计算。燃烧室内热力计算调用CEA[28]完成,考虑了化学平衡流动影响。
ηc=c0+c1·Ma0+c2·φ+c3·Ma02+
c4·Ma0·φ+c5·φ2(4)
σ4=d0+d1·Ma0+d2·φ+d3·Ma02+
d4·Ma0·φ+d5·φ2(5)
式中:ηc为燃烧效率; σ4为燃烧室总压恢复系数; φ为当量比; 系数c0~c5、d0~d5由试验确定。
1.2.4 6截面参数根据式(6)给出的尾喷管总压恢复系数模型及完全膨胀条件完成6截面尾喷管出口参数计算。
σ6=e0+e1·Ma4+e2·Ma42(6)
式中:σ6为尾喷管总压恢复系数; 系数e0~e2由CFD仿真结果确定。
完成发动机进出口截面参数计算后,采用式(7)、式(8)计算获得发动机推力、比冲。
F=(m·6v6-m·0v0)+(p6-p0)A6(7)
Isp=F/(m·6-m·0)(8)
式中:F为发动机推力; Isp为发动机比冲; m·0、v0、p0分别为发动机捕获空气流量、自由流速度、静压; m·6、v6、p6、A6分别为尾喷管出口燃气流量、速度、静压、面积。
采用上述发动机性能分析方法,对文献[8]中给出的发动机构型进行了性能计算,图2、图3分别给出了Ma0=4、Ma0=6工况计算与自由射流试验结果对比。图中可见,采用本文所建立的发动机性能计算方法,冲压模态及火箭冲压模态下发动机推力、比冲计算误差均在10%以内。
图2 Ma0=4工况计算与自由射流试验结果对比
Fig.2 Comparison between simulation and test data at Ma0=4图3 Ma0=6工况计算与自由射流试验结果对比
Fig.3 Comparison between simulation and test data at Ma0=61.3 计算方案冲压分别采用煤油、CH4、H2燃料; 火箭推力室燃料与冲压相同,火箭推力室氧化剂均为O2,火箭推力室混合比采用化学恰当比,火箭额定流量为1 kg/s,火箭推力室室压3 MPa。火箭引射模态冲压不供油,火箭冲压及冲压模态发动机当量比φ=0.8。
发动机火箭引射模态从Ma0=0、H=0 km弹道点开始工作,按照H=6.27Ma0弹道爬升至Ma0=2、H=12.54 km弹道点(Q=50 kPa); 而后以火箭冲压模态或冲压模态等动压爬升至Ma0=8、H=30.48 km弹道点。
- [1] ESCHER W,SCHNURSTEIN R. A retrospective on early cryogenic primary rocket subsystem designs as integrated into rocket-based combined-cycle(RBCC)engines[C]//29th Joint Propulsion Conference and Exhibit. Reston,Virginia:AIAA,1993.
- [2] EHRLICH C. Early studies of RBCC applications and lessons learned for today[C]//36th AIAA/ASME/SAE/ASEE Joint Propulsion Conference and Exhibit. Reston,Virginia:AIAA,2000.
- [3] CZYSZ P,LITTLE M. Rocket based combined cycle engine(RBCC):A propulsion system for the 21st century[C]//5th International Aerospace Planes and Hypersonics Technologies Conference. Reston,Virginia:AIAA,1993.
- [4] ESCHER W J D. A US history of airbreathing/rocket combined-cycle(RBCC)propulsion for powering future aerospace transports,with a look ahead to the year 2020[EB/OL].https://www.researchgate.net/publication/4690052_A_US_History_of_AirbreathingRocket_Combined-Cycle_RBCC_Propulsion_for_Powering_Future_Aerospace_Transports_with_a_Look_Ahead_to_the_Year_2020,1999.
- [5] MCDANIEL J,GOYNE C,EDWARDS J,et al.US national center for hypersonic combined cycle propulsion:an overview[C]//16th AIAA/DLR/DGLR International Space Planes and Hypersonic Systems and Technologies Conference. Reston,Virginia:AIAA,2009.
- [6] KOUICHIRO T,SADATAKE T,KANENORI K,et al.Current status of researches of the combined cycle engine at JAXA[R].ISABE-2011-1334.
- [7] ZHANG T T,WANG Z G,HUANG W,et al.The overall layout of rocket-based combined-cycle engines:a review[J].Journal of Zhejiang University-SCIENCE A,2019,20(3):163-183.
- [8] 张玫,张蒙正,刘昊. 火箭基组合循环动力研究进展[J].科技导报,2020,38(12):54-68.
- [9] ESCHER W,PRYOR D,HYDE E. Development of a 12-thrust chamber kerosene/oxygen primary rocket subsystem for an early(1964)air-augmented rocket ground test system[C]//9th International Space Planes and Hypersonic Systems and Technologies Conference. Reston,Virginia:AIAA,1999.
- [10] HYDE E,ESCHER W,RODDY J.Marquardt's Mach 4.5 supercharged ejector ramjet high-performance aircraft engine project[C]//36th AIAA/ASME/SAE/ASEE Joint Propulsion Conference and Exhibit. Reston,Virigina:AIAA,2000.
- [11] RATEKIN G. Rocketdyne RBCC engine concept development[R].ISABE-99-7179.
- [12] BULMAN M,SIEBENHAAR A. The strujet engine:exploding the myths surrounding high speed airbreathing propulsion[C]//31st Joint Propulsion Conference and Exhibit. Reston,Virigina:AIAA,1995.
- [13] THOMAS S R,PALAC D T,TREFNY C J,et al.Performance evaluation of the NASA GTX RBCC flowpath[EB/OL].https://www.researchgate.net/publication/24301438_Performance_Evaluation_of_the_NASA_GTX_RBCC_Flowpath,2001.
- [14] QUINN J. ISTAR:Project status and ground test engine design[C]//39th AIAA/ASME/SAE/ASEE Joint Propulsion Conference and Exhibit. Reston,Virginia:AIAA,2003.
- [15] OLDS J R,BRADFORD J E. SCCREAM(simulated combined-cycle rocket engine analysis module):A conceptual RBCC engine design tool[C]//33rd AIAA/ASME/SAE/ASEE Joint Propulsion Conference and Exhibit. Reston,Virginia:AIAA,1997.
- [16] BRADFORD J E,OLDS J R. Improvements and Enhancements to SCCREAM,conceptual RBCC engine analysis tool[C]//34th AIAA/ASME/SAE/ASEE Joint Propulsion Conference and Exhibit.Reston,Virginia:AIAA,1998.
- [17] OLDS J R,BRADFORD J E. SCCREAM:A conceptual rocket-based combined-cycle engine performance analysis tool[J].Journal of Propulsion and Power,2001,17(2):333-339.
- [18] OLDS J R. Launch vehicle systems analysis[EB/OL].https://www.researchgate.net/publication/27525887_Launch_vehicle_systems_analysis,1999.
- [19] MCKAMEY R,LANDRUM D. A one-dimensional engineering model for the evaluation of rocket-based combined cycle engine performance[C]//37th Joint Propulsion Conference and Exhibit. Reston,Virginia:AIAA,2001.
- [20] 黄生洪,何洪庆,何国强,等. 火箭基组合循环(RBCC)推进系统概念设计模型[J].推进技术,2003,24(1):1-5.
- [21] LEHMAN M,PAL S,SANTORO R. Experimental investigation of the RBCC rocket-ejector mode[C]//36th AIAA/ASME/SAE/ASEE Joint Propulsion Conference and Exhibit. Reston,Virginia:AIAA,2000.
- [22] 吕翔,刘佩进,何国强. RBCC发动机性能分析方法研究[J].固体火箭技术,2007,30(2):120-123.
- [23] 吕翔,何国强,刘佩进,等. RBCC发动机性能分析模型改进方法研究[J].固体火箭技术,2010,33(4):387-390.
- [24] OLDS J,BRADFORD J,CHARANIA A,et al.Hyperion:an SSTO vision vehicle concept utilizing rocket-based combined cycle propulsion[C]//9th International Space Planes and Hypersonic Systems and Technologies Conference. Reston,Virigina:AIAA,1999.
- [25] 安佳宁. RBCC(火箭基组合循环)引射模态研究[D].长沙:国防科学技术大学,2011.
- [26] 刘昊,王君. RBCC发动机火箭推力增益之探讨[J].火箭推进,2017,43(1):18-23.
- [27] 梁德旺,李博,容伟. 热完全气体的热力学特性及其N-S方程的求解[J].南京航空航天大学学报,2003,35(4):424-429.
- [28] GORDON S,MCBRIDE B J. Computer program for calculation of complex chemical equilibrium compositions and applications[EB/OL].https://www.docin.com/p-1743998778.html,1994.