作者简介:魏 鑫(1992—),男,硕士,研究领域为高超音速热防护及热管理
(College of Energy and Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China)
precooler; supercritical helium; average heat transfer coefficient; total pressure loss coefficient
备注
作者简介:魏 鑫(1992—),男,硕士,研究领域为高超音速热防护及热管理
为了掌握吸气式火箭发动机(SABRE)空气预冷器的流动换热特性,为设计相应类型的预冷器提供技术基础,针对SABRE预冷器最小周期性单元,以数值方法研究了管间距、管排数、空气入射角度及氦气/空气热容量比对预冷器流动换热的影响。研究结果表明:增大管排数和减小管间距,能够增大预冷器换热功率,降低空气出口温度,但会降低空气侧、氦气侧平均换热系数,减弱对流换热能力,增大空气侧总压损失。空气入射角度对空气侧、氦气侧换热影响微小,但对空气侧总压恢复系数影响显著。增大氦气/空气热容量比能够降低空气侧总压损失,增大空气侧、氦气侧平均换热系数,降低空气出口温度。
In order to study the flow and heat transfer characteristics of air precooler in Synergistic Air-Breathing Rocket Engine(SABRE)and provide the technical basis for designing precooler, the effects of tube pitch, tube row, air incidence angle and helium/air heat capacity ratio on the flow and heat transfer of air precooler were numerically studied.The results show that increasing the number of tube rows and reducing the tube spacing can increase the heat transfer power of the precooler and reduce the air outlet temperature.However, they also reduce the average heat transfer coefficient on the air side and helium side, weaken the convective heat transfer coefficient and increase the total pressure loss coefficient of air side.The air incidence angle has little effect on the heat transfer of both sides, but has a significant effect on the total pressure recovery coefficient of air side.Increasing the helium/air heat capacity ratio can reduce the total pressure loss coefficient on the air side, increase the average heat transfer coefficient on both sides, reduce the air outlet temperature and enhance the pre-cooling effect.
引言
协同吸气式火箭发动机(SABRE)在高马赫数飞行时,因进气滞止温度过高,空气难以压缩,动力系统性能会急剧下降。为了改善吸气式动力系统高飞行马赫数时的性能,通过在进气道中加预冷器,以降低进气温度,减小压气机功率,这对扩展飞行器飞行包线、提高发动机推重比具有重要意义[1-5]。SABRE在氢气和空气之间引入闭式氦气循环,借助低温氦气来冷却高温空气[6-8],在较短的时间内实现空气和低温氦气高效换热,大温差、高效紧凑、轻质的空气预冷器是协同吸气式火箭发动机的核心技术[9-10]。
目前,SABRE型发动机方案已经从SABRE3发展到SABRE4,相比之下,后者优化了循环方式,方案可行性大大增加。SABRE3型预冷器目前存在一些工程性问题,如:空气冷却到-150 ℃会结霜阻塞通道,增设除霜装置会使得预冷器结构复杂、重量增加; 预冷器的长度、体积限制使得预冷器难以达到换热要求。因此,在SABRE3方案基础上,2015年,REL(Reaction Engines Ltd)公司公布了优化后的“佩刀”发动机方案,即SABRE4[11-12],SABRE4预冷换热器在Ma=1.9~2.8仅预冷器HX1工作,在Ma=2.8以上预冷器HX1和HX2两者同时工作[13]。预冷换热器出口温度由SABRE3的-150 ℃提高到0 ℃以上,省去了结霜控制系统。
针对SABRE型发动机空气预冷器,国内外学者做了许多研究。最初,基于“佩刀”发动机先进的预冷概念,在上世纪90年代,布里斯托尔大学研制了由直径0.38 mm微细管组成的实验型预冷器模块JMHX[14-15],研究了JMHX冷却剂为氦气、氮气时在不同流量下的换热性能,并分析了经典理论中摩擦因子和传热系数的公式在微通道换热器中的适用性。Helen Webber[16]实验研究了与SABRE预冷器相似的叉排布局圆柱管束在低速风洞下的对流换热能力,同时测得了空气流过圆柱管束时的流场特性。Huan Lee[17]建立了紧凑式换热器性能试验台,对紧凑式换热器进行了性能测试,初步实验结果表明该换热器具有优良的性能,实验结果与设计值吻合较好。李晨沛[18]针对复合发动机的预冷器,研究管内流体速度、管外流体速度、入射角度以及管间距对于空气侧换热性能的影响机制,结果表明:通过增大管内制冷介质的流量,缩小管间距以及空气垂直入射角度等均能提高预冷器换热效率。程惠尔[19]通过对选定参数的空气预冷换热器作了全工况跟踪计算,给出了空气通道阻力特性的数值结果。
综上,国内外学者通过实验、数值模拟、一维分析的方法针对预冷器做了一些研究,但研究对象大部分是传统的叉排直管束,对于SABRE特殊形式的预冷器,由于结构形式复杂,对其流动换热规律研究很少。为了深入研究空气预冷器的流动换热特性,本文针对SABRE型发动机预冷器,建立预冷器流动换热的最小周期性模型,研究了管间距、管排数、空气入射角度及氦气/空气热容量比对于预冷器流动换热的影响规律,为空气预冷器的工程设计提供参考。
1 物理模型和计算方法
1.1 物理模型和网格划分图1为空气预冷器结构示意图,空气预冷器由弯曲状的预冷片、总管及支管组成,共36个预冷片,预冷片绕圆心旋转形成辐射状排布。低温高压的氦气从内环的总管流入,之后流入每根支管,支管连接着预冷片,每个预冷片由四排毛细管组成叉排管阵,高温空气流经管与管之间的缝隙逐渐被冷却。由于预冷器在圆周方向具有周期性重复的特点,氦气和空气的流动换热也具有周期性重复的特点,对于预冷器的最小周期性单元,其流动换热能够代表整个预冷器的流动换热,因此可以取最小的周期性计算单元(图1黄色扇形区域)进行研究。
图2为预冷器计算单元模型,计算单元内径873 mm,外径1 250 mm。每一排预冷片都是由4排毛细管组成的,毛细管外径1 mm,壁厚0.05 mm。管排数为预冷器计算单元内空气流过的预冷片数目,管阵叉排布局以等边三角形布置,管间距为相邻两个毛细管圆心之间的距离。氦气在1和1',2和2',3和3',…,8和8'的流动状态相同,每根管子的出口都是旋转对应的下一排管子的入口,氦气在计算域中往返流动,计算域内的毛细管可以等效为一根完整的毛细管。计算单元内扇形夹角为10°,由于夹角相对较小,空气在扇形区域上下边界处空气的流动可以近似看成相同的。为了减少计算量,简化计算结构,预冷器轴向选择四根半管作为轴向的最小周期性单元,同时体现了叉排布局结构,对于四根半管的截面,设为对称边界条件,使得半管的流动状态与完整管相似。
1.2 边界条件和计算方法计算模型空气入口采用压力入口边界条件,给定总温总压,空气出口采用压力出口,并给定出口目标流量。氦气入口选择速度入口边界条件,入口静温为340 K,氦气出口采用压力出口边界条件,给定出口静压20 MPa。周向两侧壁面除了氦气入口和氦气出口外均设为周期性边界条件,上下截面设为对称边界条件,管壁设为耦合壁面。由于氦气在SABRE发动机中为闭式循环[20],氦气总流量为定值。根据文献[21]可知,SABRE预冷器空气侧流量随飞行弹道基本保持不变。因此在设置边界条件时考虑到了预冷器空气和氦气流量恒定的工作特性,并且预冷器总长度设定为定值4 m,氦气和空气总的流量保持不变,其中
mHe=(MHe*w)/(2*36*l)(1)
mAir=(MAir*w)/(2*36*l)(2)
式中:mHe为预冷器计算单元氦气流量; MHe为预冷器的氦气总流量; mAir为预冷器计算单元空气流量; MAir为预冷气的空气总流量; w为管间距; l为预冷器长度。
预冷器计算单元的空气流量和氦气流量随管间距增大而增大,但氦气和空气总的流量保持不变。表1为预冷器计算单元的研究变量取值,针对管排数、管间距、氦气/空气热容量比及空气入射角度进行研究,通过控制单一变量,研究各个变量对于空气侧流动换热的影响规律。其中,氦气/空气热容量比定义为(mHecp,He)/(mAircp,Air)。
本文采用Fluent16.0进行数值模拟计算,数值模拟中空气和超临界氦气的物性参数由NIST软件查得,空气设为理想气体,比热容、粘性系数及导热率采用多项式函数进行拟合,由于超临界氦气在毛细管内部流阻很小,压力变化很小,超临界氦气密度设为温度的函数。
1.3 数值方法验证参考文献[14]验证数值计算方法,文献对紧凑型预冷器模块JMHX进行了实验研究,JMHX外形结构为40 mm*40 mm*40 mm大小,共415个直径0.38 mm的微细管束,呈10排83列排布,管内冷却介质为低温高压氦气或者氮气。选择文献中内部介质为氦气、热容量比为1的实验进行对比分析,图3为不同湍流模型下各排管氦气出口温度计算结果与实验结果比较,可以发现,靠近空气入口的前几排管束计算结果与实验结果误差较大,后几排管束的温度计算结果与实验结果误差较小,这可能是因为前面的管排管内外介质温差较大,毛细管的热缩性较大,几何尺寸变化较大,而数值模拟结果没有考虑形变的影响,因此前几排的误差较大,相比来看,SST-kω湍流模型比其他湍流模型更接近实验结果,因此本文数值计算选择SST-kω湍流模型。
2 计算结果与讨论
当氦气压力大于临界压力(0.228 MPa),氦气温度大于临界温度(5.19 K),此时氦气状态为超临界状态[22]。预冷器毛细管内氦气压力远大于临界压力,温度远大于临界温度,因此预冷器内氦气为超临界状态。超临界氦气为单相[23],不会因为温度升高而呈现两相状态,因此本文预冷器管内超临界氦气属于单相换热范畴。
2.1 管间距对于预冷器流动与换热的影响针对管间距对于预冷器流动与换热的影响,本文研究了来流工况不变的条件下,管间距对于空气侧总压损失系数以及空气侧、氦气侧平均换热系数的影响规律。
预冷器空气侧总压损失系数定义为
σ=(p*in-p*out)/(p*in)(3)
式中p*in和p*out 为预冷器空气入口总压和出口总压。
氦气侧平均换热系数基于管内湍流对流换热实验关联式Dittus-Boelter公式可得
Nuf=0.023Re0.8fPr0.4f(4)
空气侧平均换热系数[24],通过以下方法计算。预冷器在传热达到稳态后,换热量
Q=mAircp,AirΔTAir=m<sub>Hecp,HeΔTHe(5)
式中ΔTHe和ΔTAir为氦气、空气进出口温差,预冷器总传热系数为
K=Q/(Aψ(ΔTm)ctf)(6)
式中:A为预冷器的壳侧换热面积;(ΔTm)ctf为将给定的冷热流体的进出口温度布置为逆流时的对数平均温差; ψ为小于1的修正系数,ψ的计算参考文献[25],根据预冷器流体流动形式,选择n个相同的叉流热交换器流动形式计算ψ值,经计算ψ的范围为0.98~1,说明预冷器流体流动形式接近逆流布置。基于壳侧换热面积的总传热系数[26]
K=1/(1/(hHe)+Rw+RS+1/(hAir)(dO)/(di))(7)
式中do和di为管外径和管内径。由于毛细管壁厚很薄,导热热阻Rw以及污垢热阻RS忽略不计,可以计算得空气侧平均换热系数
hAir=[(1/K-1/(hHe))(di)/(do)]-1(8)
图4为预冷器空气侧总压损失系数随管间距变化关系,可以看出,随着管间距从1.5 mm增大到3.5 mm,总压损失系数从2.6%降低到0.5%,且随着管间距增大,总压损失系数下降速率减缓。
横掠管束压降经验关系式[25]Δp=1.5Re-0.2ρV2max((μf)/(μw))0.14N(9)
其中
Re=ρVmax D/μf
Vmax =Vf(s1)/(s1+D)
式中:N为管排数; D为管外径; μf为管壁平均温度; Vmax为管束间的最大速度; Vf为来流速度; s1为管间距。当管间距变化时,空气来流速度基本不变,管束间的最大速度变化较大,而定性温度下的物性参数变化很小,因此压降主要与管束间的最大流速有关,且压降变化趋势与最大流速变化趋势相同,压降近似与管束间最大流速的1.8次方成正比。
图5为预冷器空气侧、氦气侧平均换热系数随管间距变化关系,随着管间距从1.5 mm增大到3.5 mm,空气侧平均换热系数增大了约4.4倍,氦气侧平均换热系数增大了约2倍,管内外对流换热能力增强,可以看出,氦气侧换热系数远大于空气侧换热系数。预冷器氦气和空气总流量不变,在给定预冷器长度不变的约束下,随着管间距增大,氦气入口流速增大,雷诺数增大,而流体普朗特数(Pr)变化很小,根据管内换热关联式(4),氦气侧平均换热系数增大; 而空气侧平均换热系数主要由预冷器总体换热系数决定,随管间距增大,预冷器换热面积与对数平均温差乘积的减小速率大于换热功率的减小速率,根据总体换热系数的表达式(6),总体换热系数增大,空气侧平均换热系数因而增大。
图6为预冷器空气出口温度和预冷器换热功率随管间距变化关系,随着管间距从1.5 mm增大到3.5 mm,空气出口温度升高,氦气出口温度降低,换热功率降低了31%,预冷效果减弱。虽然减小管间距会降低预冷器空气侧、氦气侧换热系数,但预冷器总换热面积增大了,而空气、氦气总流量不变,因此减小管间距能够提高预冷效果,降低空气出口温度。2.2 管排数对于预冷器流动与换热的影响针对管排数对于预冷器流动与换热的影响,本文研究了来流工况不变的条件下,管排数对于空气侧总压损失系数以及空气侧、氦气侧平均换热系数的影响规律。
图7为预冷器空气侧总压损失系数随管排数变化关系。可以看出,随管排数从7排增大到15排,预冷器空气侧总压损失系数从1.5%增大到了7%,且近似呈线性上升。管间距不变,由于空气侧流量为定值,侧管束间的最大速度不变,定性温度下的物性参数也基本不变,由横掠叉排管束压降经验关系式(9)可知,空气侧压降与管排数呈线性关系,因此预冷器空气侧总压损失系数随管排数增大近似线性上升。
图8为预冷器空气侧、氦气侧平均换热系数随管排数变化关系,随管排数从7排增大到15排,空气侧平均换热系数减小了57%,氦气侧换热系数略微减小,减小幅度小于4%。随管排数增大,对数平均温差与换热面积乘积的增大速率大于换热功率增大的速率,根据总体换热系数的表达式(6),总体换热系数减小,空气侧平均换热系数因而减小。氦气侧平均换热系数主要由雷诺数决定,随管排数增大,氦气侧Pr数基本不变,而雷诺数略微减小,根据管内换热关联式(4),氦气侧平均换热系数略微减小。
图9为预冷器空气出口温度和预冷器换热功率随管排数变化关系,随管排数从7排增大到15排,空气出口温度降低,氦气出口温度升高,预冷器换热功率增大了20%,预冷效果增强; 同时可以看出,随管排数增大,预冷器换热功率上升速度减缓,空气出口温度下降速度减缓,说明增大管排数对于空气预冷效果的收益减小,因此管排数不宜过多。2.3 空气入射角度对于预冷器流动与换热的影响针对空气入射角度对于预冷器流动与换热的影响,本文研究了来流总温、总压不变的条件下,空气入射角度对于空气侧总压损失系数以及空气侧、氦气侧平均换热系数的影响规律。
图 10为空气侧总压损失系数随空气入射角度的变化关系,在空气入射角度小于90°时,入射角度对于空气侧总压损失系数的影响较小,总压损失系数变化小于5%; 空气入射角度大于90°时,随着入射角度增大,空气侧总压损失系数逐渐增大,入射角度从90°增加到150°,空气侧总压损失系数增大了20%。
图 10 空气侧总压损失系数随空气入射角度的变化关系
Fig.10 Effects of incidence angle on total pressure loss coefficient of air side2.4 氦气/空气热容量比对于预冷器流动与换热的影响针对氦气/空气热容量比对于预冷器流动与换热的影响,本文研究了来流总温、总压以及空气侧流量不变的条件下,改变氦气总流量,也就是改变氦气/空气热容量比对于空气侧总压损失系数以及空气侧、氦气侧平均换热系数的影响规律。
图 13为空气侧总压损失系数随氦气/空气热容量比的变化关系,随着氦气/空气热容量比从1.06增大到1.64,空气侧总压损失逐渐下降,总压损失下降了13%。这是因为管内氦气流量增加使得空气侧预冷效果增强,空气温度降低,空气密度进而增大,而预冷器空气侧流量不变,因而沿程空气流速降低,使得空气侧总压损失减小。
图 13 空气侧总压损失系数随氦气/空气热容量比的变化关系
Fig.13 Effects of heat capacity ratio on total pressure loss coefficient of air side3 结论
1)随着管间距从1.5 mm增大到3.5 mm,总压损失系数从2.6%降低到0.5%,且随着管间距增大,总压损失系数下降速率减缓; 空气侧平均换热系数增大了约4.4倍,氦气侧增大了约2倍,但换热功率降低了31%,空气出口温度升高,预冷效果减弱。
2)随着管排数从7排增大到15排,预冷器空气侧总压损失系数从1.5%增大到了7%,且近似呈线性上升。空气侧平均换热系数减小了57%,氦气侧减小了4%。预冷器换热功率增大了20%,预冷效果增强。
3)在空气入射角度小于90°时,入射角度对于空气侧总压损失系数的影响小于5%; 空气入射角度大于90°时,随着入射角度增大,空气侧总压损失系数逐渐增大,入射角度从90°增加到150°,空气侧总压损失系数增大了20%; 空气入射角度对于空气侧、氦气侧换热几乎没有影响。
4)随着氦气/空气热容量比从1.06增大到1.64,空气侧总压损失下降了13%,空气侧平均换热系数增大了30%,氦气侧平均换热系数增大了40%,预冷效果增强。
- [1] VERDTR AETE D,HENDRICK P.Hydrogen fueled precooled air breathing engines for hypersonic aircraft and spaceplanes[C]//20th ISABE.Sweden:Chalmers University,2011.
- [2] 张友法, 张文文, 郑日恒, 等.高超声速组合发动机预冷器抗结霜涂层技术研究[J].推进技术, 2017, 38(2): 463-470.
- [3] 聂万胜, 周思引, 雷旭. 协同吸气式火箭发动机研究进展[J]. 装备学院学报, 2016, 27(6): 57-64.
- [4] LONGSTAFF R, BOND A. The SKYLON project[C]//17th AIAA International Space Planes and Hypersonic Systems and Technologies Conference. San Francisco, California. Reston, Virigina: AIAA, 2011.
- [5] 马海波, 张蒙正. 预冷空气类动力系统发展历程浅析[J]. 火箭推进, 2019, 45(2): 1-8.
- [6] 邓帆, 谭慧俊, 董昊, 等.预冷组合动力高超声速空天飞机关键技术研究进展[J].推进技术, 2018, 39(1): 1-13.
- [7] WEBBER H, FEAST S, BOND A. Heat exchanger design in combined cycle engines[J]. Journal of the British Interplanetary Society, 2009, 62(4):122-130
- [8] WEBBER H, TAYLOR N. Tunnel development for heat transfer analysis in compact heat exchangers[C]//27th AIAA Aerodynamic Measurement Technology and Ground Testing Conference. Chicago, Illinois. Reston, Virigina: AIAA, 2010.
- [9] 邹正平, 刘火星, 唐海龙, 等.高超声速航空发动机强预冷技术研究[J].航空学报, 2015, 36(8): 2544-2562.
- [10] 汪元, 王振国. 空气预冷发动机及微小通道流动传热研究综述[J]. 宇航学报, 2016, 37(1): 11-20.
- [11] DISSEL A.SABRE technology development: IAC-16,C4.9.2.[R].[S.l.]:Reaction Engines Inc,2016.
- [12] 张连庆, 刘博, 李浩悦. “佩刀”发动机技术进展分析[C]//中国航天第三专业信息网第三十八届技术交流会暨第二届空天动力联合会议论文集.大连:[s.n.],2017.
- [13] 周建兴, 张浩成, 高启滨, 等. 基于SABRE技术的高超声速预冷飞行器应用分析[J]. 推进技术, 2018, 39(10): 2196-2206.
- [14] MURRAY J J, HEMPSELL C M, BOND A.An experimental precooler for airbreathing rocket engines[J].Journal of the British Interplanetary Society, 2001, 54(5):199-209.
- [15] MURRAY J J, GUHA A, BOND A.Overview of the development of heat exchangers for use in air-breathing propulsion pre-coolers[J].Acta Astronautica, 1997, 41(11): 723-729.
- [16] WEBBER H, TAYLOR N.Tunnel development for heat transfer analysis in compact heat exchangers[C]//27th AIAA Aerodynamic Measurement Technology and Ground Testing Conference.Reston, Virigina: American Institute of Aeronautics and Astronautics, 2010.
- [17] LEE H, MA S, CHEN Y M, et al.Experimental study on compact heat exchanger for hypersonic aero-engine[C]//21st AIAA International Space Planes and Hypersonics Technologies Conference.Reston, Virginia: American Institute of Aeronautics and Astronautics, 2017.
- [18] 李晨沛, 王跃社, 王海军, 等.复合发动机预冷器换热特性研究[J].工程热物理学报, 2017, 38(4): 811-816.
- [19] 程惠尔.空气预冷换热器的性能计算[J].推进技术, 1993, 14(6): 44-50.
- [20] 郭海波, 肖洪, 南向谊, 等.复合预冷吸气式火箭发动机热力循环分析[J].火箭推进, 2013, 39(3): 15-20.
- [21] ZHANG J Q, WANG Z G, LI Q L.Thermodynamic efficiency analysis and cycle optimization of deeply precooled combined cycle engine in the air-breathing mode[J].Acta Astronautica, 2017, 138: 394-406.
- [22] 张鹏, 王如竹.超流氦传热[M].北京: 科学出版社, 2009.
- [23] 张鹏, 黄永华, 陈国邦.氦-4和氦-3及其应用[M].北京:国防工业出版社,2006.
- [24] 王冬兰, 朱卫兵, 平登科, 等.燃-滑油换热器壳侧换热特性研究[J].推进技术, 2016, 37(5): 900-906.
- [25] 钱滨江.简明传热手册[M].北京: 高等教育出版社, 1983.
- [26] CHANG S W.Heat transfer of orthogonal-mode reciprocating tube fitted with twisted-tape[J].Experimental Heat Transfer, 2000, 13(1): 61-86.