自然循环预冷是低温液体火箭发动机预冷的重要发展趋势,其是利用低温推进剂在两垂直管路中的密度差驱动自发性流动实现对火箭发动机及其增压输送系统进行预冷,具有节省推进剂消耗量、简化工艺流程等特点。然而,低温液体管路预冷过程内在机理极其复杂,预冷管路内可能产生剧烈气化,引起较大压力波动。同时,伴随预冷过程不断进行,循环驱动力逐渐减弱,预冷效率降低。因此,研究低温液体管路充填及预冷特性具有重要的现实意义。

目前,国内外对低温火箭发动机及低温液体的管路预冷研究已有相关的文献报道。文献[1]～[3]主要针对低温管路内部传热机理及流型转化进行了理论分析和数值研究; 文献[4]～[7]采用实验方法研究了低温管路预冷过程的流动特性; 文献[8]分析了低温流体预冷竖直管路的高速再淹没传热特征,研究了液体火箭发动机自然循环预冷过程中的温降特性及其循环驱动力变化趋势。结合国内外文献发现,针对低温液体火箭发动机自然循环预冷管路特性的研究相对较少。本文通过构建一维均相平衡态流体数学模型,建立低温流体管路传热计算模型,重点研究液氢燃料火箭发动机循环预冷管路非稳态过程的变化规律。

以芯一级低温氢氧液体火箭发动机液氢预冷管路为研究对象,重点进行液氢加注充填过程管路内部流动及传热的非稳态特性分析。图2为输送管进口速度及循环流量变化图,由图2可知,循环管路中输送管进口速度先快速增至最大值,而后迅速下降,并逐步趋于稳定值。这是由于循环管路中的流动压降主要受流速和含气率两因素影响,在预冷过程开始阶段,其出口为全气相,循环管路中流动阻力很大,使得输送管进口速度较小。由于回流管出口全气相及大部分气相的时间很短,在这一段阻力损失变化较大,使得速度变化比较陡峭。在此之后,流动阻力主要受流速影响,随着预冷介质产生气化的位置逐渐后移,流动产生的动力及阻力逐渐减小。

图2 输送管进口速度及循环流量变化
Fig.2 Velocity and mass flux of feeding pipe inlet

图3给出了回流管出口质量/截面含气率随预冷时间的变化曲线。由于液氢过冷度很小,在预冷开始阶段,回流管出口为全气相,质量含气率、截面含气率均为1,然后迅速减小。非稳态过程中回流管出口质量/截面含气率变化主要受输送管进口速度和预冷管路壁面温度的影响。对于非稳态传热,在开始预冷阶段管路温度变化很快,随后逐渐减慢,基本成指数曲线关系变化,这使得回流管出口质量含气率、截面含气率具有类似的特征。边界条件一定时,输送管进口速度增加时,回流管出口质量含气率、截面含气率将会减小,反之亦然。

图4给出了不同时刻时循环管路沿程各位置点的截面含气率变化。从图中可以看出,气化点位置随时间的推移逐渐往后推移,约5 000 s以后截面含气率开始趋于相对稳定状态。

图3 预冷介质回流管出口含气率变化
Fig.3 Variation of void fraction at return pipe outlet

图4 循环管路中各点截面含气率变化
Fig.4 Variation of void fraction at each point in circulation pipe

图5给出了输送管内壁面一特征点温度变化曲线,该特征点依次经历四个典型阶段:膜态沸腾、过渡沸腾、核态沸腾以及单相液体强制对流。循环管路上其他各点与该特征点具有相类似的温度变化曲线,只是相对应的各转折点温度及相应的时间略有区别而已。由图可知,由于T_FL和T_ONB很低,加上过渡沸腾阶段与核态沸腾阶段涉及的温度区间太小,因而壁面温度的降低主要集中在膜态沸腾阶段。

图5 输送管特征点预冷过程传热工况过渡
Fig.5 Heat transfer transition during the precooling process of the characteristic points in feeding pipe

图6给出了泵管道内壁面一特征点处的温度变化曲线,达到稳定状态时该点壁温高于核态沸腾起始点温度,处于两相流动区,故没有单相对流阶段。初始阶段,管壁温度下降很快,而后为一定值。这是由于初始阶段管内全部为气相,预冷介质处于过热状态,当管内很快进入两相状态后,由于管路压力为定值,使得该特征点处温度此后为一定值。

图6 泵管道特征点预冷过程传热工况过渡
Fig.6 Heat transfer transition during the precooling process of the characteristic points in pump pipe

通过构建液体火箭发动机自然循环预冷的一维均相平衡态流体数学模型,分析了非稳定状态时液氢循环预冷管路内介质流速、流量、含气率和壁面温度的变化规律,获得了不同管段特征点的温度曲线,为提高低温发动机预冷效率和优化管路绝热设计提供了重要的理论依据。