文章快速检索  
  高级检索
基于燃气射流冷凝的氧路系统频率特性研究
张淼1,2, 李斌3, 邢理想1     
1. 西安航天动力研究所, 西安 710100;
2. 液体火箭发动机技术重点实验室, 西安 710100;
3. 航天推进技术研究院, 西安 710100
摘要: 基于热力学不平衡两流体六方程模型模拟了液氧煤油补燃循环发动机燃气射流在氧路系统泵间管路液氧流体中的冷凝过程,获得了燃气射流冷凝特征参数沿管路流向的分布规律。根据Rayleigh-Plesset方程和燃气射流特性建立了泵间管路燃气冷凝过程的传递函数模型,并与氧路管路、泵等组件模型联立求解,分析了发动机氧路系统频率特性。研究了在不同氧路入口压力和液氧温度边界条件下,泵间管路燃气射流冷凝过程对氧路系统频率特性的影响,仿真结果表明高入口压力和过冷液氧改变了射流气体的惯性和柔度,使得氧路系统特征频率增大。不同边界条件下发动机热试车结果表明,提高氧路入口压力或降低液氧温度使得氧路系统频率从8.3 Hz提高至11 Hz,与数值仿真结果一致,验证了泵间管燃气射流冷凝过程是影响氧路系统频率特性的重要环节。
关键词: 氧路系统    频率特性    燃气射流    冷凝    两流体模型    
Study on frequency characteristics of oxygen feed system based on gas jet and condensation
ZHANG Miao1,2, LI Bin3, XING Lixiang1     
1. Xi'an Aerospace Propulsion Institute, Xi'an 710100, China;
2. Science and Technology on Liquid Rocket Engine Laboratory, Xi'an 710100, China;
3. Academy of Aerospace Propulsion Technology, Xi'an 710100, China
Abstract: Based on the one-dimensional thermodynamic disequilibrium two-fluid six-equation model, the condensation of gas jet is simulated. The simulation happens in the pipeline between pumps of the liquid oxygen/kerosene staged combustion cycle rocket engine system, obtaining the characteristic parameters along pipeline flow direction. Firstly, based on Rayleigh-Plesset equation and the characteristics of gas jet and condensation distribution, a transfer function model is established. The condensation model with liquid oxygen pipeline, pump, and other component models are solved simultaneously, and the frequency characteristics of oxygen pipeline system of the engine are analyzed. Secondly, the influence of gas jet and condensation process on the frequency characteristics of the oxygen feed system at different inlet pressure levels and liquid oxygen temperature boundary conditions are studied. The simulation results show that the high inlet pressure and subcooling liquid oxygen changes the inertia and flexibility of jet bubble, and increases the characteristic frequency of the oxygen feed system. The engine firing test results under different boundary conditions show that increasing the oxygen feed system inlet pressure or decreasing the liquid oxygen temperature increases the frequency of oxygen feed system from 8.3 Hz to 11 Hz. The result is consistent with the numerical simulation results, verifying that gas jet and condensation of pump pipeline is an important process affecting the frequency characteristics of the oxygen feed system.
Keywords: oxygen feed system    frequency characteristics    gas jet    condensation    two-fluid model    

液氧煤油补燃循环发动机液氧路系统频率特性具有典型的低频特征[1]。大推力液氧煤油补燃循环发动机采用富氧燃气驱动预压涡轮泵提高发动机性能,富氧燃气从射流孔进入泵间管液氧中发生冷凝形成复杂的两相流动过程,产生剧烈的传热和相变,并存在流体不稳定振荡现象[2]。燃气射流产生大量的气泡以及气液两相流动动态过程对氧路系统频率特性影响很大,特别是在发动机入口边界条件发生变化时,即高入口压力和过冷液氧条件对泵间管两相流动过程的影响未开展深入研究,对氧路系统的频率特性影响程度也不明确。

这种过热蒸汽在液体中通过射流形成直接接触的冷凝过程也出现在电力、制冷、核工业和石油化工等许多工业领域,最早Kerney[3]和Weimer[4]等对单喷嘴水蒸气射流冷凝过程进行了试验研究,总结了蒸汽射流无量纲冷凝特征长度的经验公式。西安交通大学多相流实验室武心壮等[5-7]进行了超声速水蒸气射流冷凝试验,获得了无量纲冷凝特征长度的经验公式。关于蒸汽射流冷凝的研究都是集中在水蒸气在过冷水中的冷凝过程,对低温液氧中氧蒸汽的射流冷凝过程的试验研究和理论研究较少。Pilipenko等[8]在实验室状态下开展了过热氧气射流冷凝试验,获得了冷凝特征长度经验公式。

蒸汽射流冷凝过程的数值仿真研究主要基于流体体积函数(VOF)三维流场仿真。Li等[9]研究了亚声速蒸汽喷入过冷水的冷凝过程,对蒸汽冷凝区域类型分布特性进行验证,蒸汽射流处于振荡冷凝区时存在周期性的低频振荡。Shah等[10]研究了超声速喷嘴下蒸汽喷射入过冷水中的冷凝过程,在蒸汽质量流量密度较大时,随着液体温度升高,蒸汽冷凝边界越来越不清晰,冷凝长度大幅增加。宋纪元和陈听宽[11]基于热不平衡两流体模型研究了一维管流临界两相流动,得出了气液流动热不平衡动力学特性。核工业中基于RELAP5软件拓展了热不平衡两流体在考虑热力过程两相流求解中的应用范围[12-14]。陈二锋等[15]基于两相流压力波速模型研究了泵间管气液两相流压力波速传播特征,并揭示了两相流传播特征。薛帅杰等[16]在研究离心喷嘴自激振荡特性试验中,推测出旋流腔内气涡与液膜两相流耦合过程主导了喷嘴低频自激振荡过程。张国渊等[17]提出了低温高速离心泵的机械密封内低温介质汽化形成两相流会影响离心泵运转稳定性,形成低频自激振荡现象。在特定通道内的气液两相流动动力学过程普遍存在低频自激振荡现象,因此研究蒸汽射流冷凝过程的自激振荡现象是管路系统动态特性研究的重要课题之一。

李斌[1]、邢理想[18]和刘上[19]等对液氧煤油补燃循环发动机的氧路系统频率特性开展了深入的理论研究,提出并改进了系统传递函数模型,研究主要集中在热力组件建模及燃烧过程对系统频率的影响。张青松[20]和朱平平[21]等基于动力学系统模型研究了蓄压器模型对氧路输送系统频率特性的影响,研究了大型运载火箭纵向耦合振动(POGO)特性。上述研究对不同入口压力和液氧温度条件下的氧路系统频率特性影响都未开展深入研究。本文的研究重点是基于热不平衡两流体模型研究泵间管燃气射流冷凝过程的分布特性,以及基于气泡动力学方程建立动态特性传递函数模型,研究不同入口压力和液氧温度条件下泵间管燃气射流冷凝过程对氧路系统频率的影响。

1 氧路系统传递函数模型 1.1 氧路系统

液氧煤油补燃循环发动机氧路系统示意图如图 1所示,液氧从贮箱开始经过输送管路进入发动机预压泵,驱动氧化剂预压涡轮后的高温燃气从射流孔进入主路液氧中直接冷凝,泵间管路混合流体形成热不平衡两相流动,燃气完全冷凝后单相液氧进入氧主泵增压后进入富氧燃气发生器形成富氧燃气驱动主涡轮,最后在燃烧室再次燃烧形成推力。本文中主要研究氧入口边界条件对频率特性的影响,而主泵后参数受边界条件影响较小,因此将氧路系统传递函数模型简化为贮箱到氧主泵出口的液路系统频率特性仿真模型的假设是合理的。

图 1 发动机氧路系统示意图 Fig. 1 Schematic diagram of engine oxygen feed system

泵间管路混合流体两相流动过程复杂,燃气射流结构如图 2所示,燃气与液氧工况条件见表 1。泵间管路燃气射流后在主路液氧中形成离散的气泡,并在冷凝过程中气泡直径逐渐变小直至完全消失。由于燃气中含质量分数为96%的气氧,其余为水蒸气和CO2,为了简化计算,采用气氧物性参数代替燃气进行仿真计算。

图 2 泵间管路燃气射流冷凝示意图 Fig. 2 Schematic diagram of gas jet and condensation in pump pipeline
表 1 燃气和液氧的参数 Table 1 Parameters of gas and liquid oxygen
参数 数值
液氧入口流量/(kg·s-1) 288
液氧入口温度/K 82~97
贮箱入口压力/MPa 0.25~1.2
泵间管入口压力/MPa 1.2~2.1
燃气入口流量/(kg·s-1) 11
燃气入口温度/K 260
1.2 液体管路模型

基于一维流动的波动方程、连续和状态方程所建立的液路管路分布参数传递矩阵为

$ \left[ {\begin{array}{*{20}{c}} {{\rm{ \mathsf{ δ} }}{P_2}}\\ {{\rm{ \mathsf{ δ} }}{Q_2}} \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} {\rm{ch}\left( {\frac{{sl}}{a}} \right)}&{ - {Z_{\rm{c}}} \cdot \rm{sh}\left( {\frac{{sl}}{a}} \right)}\\ { - \frac{1}{{{Z_{\rm{c}}}}}\rm{sh}\left( {\frac{{sl}}{a}} \right)}&{\rm{ch}\left( {\frac{{sl}}{a}} \right)} \end{array}} \right]\left[ {\begin{array}{*{20}{l}} {{\rm{ \mathsf{ δ} }}{P_1}}\\ {{\rm{ \mathsf{ δ} }}{Q_1}} \end{array}} \right] $ (1)
 

式中:δP、δQ分别表示压力和质量流量的无量纲变化率,下标1表示入口,下标2表示出口;s为拉普拉斯算子;l为管路长度;a为管路流体声速;Zc为无量纲特征阻抗;sh和ch分别为双曲正弦和双曲余弦函数。

液路局部阻力元件的传递矩阵为

$ \left[ {\begin{array}{*{20}{c}} {{\rm{ \mathsf{ δ} }}{P_2}}\\ {{\rm{ \mathsf{ δ} }}{Q_2}} \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} {\frac{{{p_1}}}{{{p_2}}}}&{ - \frac{{2\left( {{p_1} - {p_2}} \right)}}{{{p_2}}}}\\ 0&1 \end{array}} \right]\left[ {\begin{array}{*{20}{l}} {{\rm{ \mathsf{ δ} }}{P_1}}\\ {{\rm{ \mathsf{ δ} }}{Q_1}} \end{array}} \right] $ (2)
 

式中:p1p2分别为入口和出口压力。

液路分支传递矩阵为

$ \left[ {\begin{array}{*{20}{c}} {{\rm{ \mathsf{ δ} }}{P_2}}\\ {{\rm{ \mathsf{ δ} }}{Q_2}} \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} 1&0\\ { - \frac{{\bar p}}{{\bar q}} \cdot \frac{1}{{{Z_k}}}}&1 \end{array}} \right]\left[ {\begin{array}{*{20}{l}} {{\rm{ \mathsf{ δ} }}{P_1}}\\ {{\rm{ \mathsf{ δ} }}{Q_1}} \end{array}} \right] $ (3)
 

式中:pq分别为平均压力和平均质量流量;Zk为分支管输入阻抗。

1.3 气蚀泵模型

考虑泵气蚀过程,基于试验获得的动态增益m+1、气蚀柔性Cb、质量增益Mb、泵的惯性Lp和泵的阻力Rp等参数建立气蚀泵传递函数,其表达式为

$ \begin{array}{l} \left[ {\begin{array}{*{20}{c}} {{\rm{ \mathsf{ δ} }}{P_2}}\\ {{\rm{ \mathsf{ δ} }}{Q_2}} \end{array}} \right] = \\ \;\;\;\;\;\;\;\;\left[ {\begin{array}{*{20}{c}} {m + 1 + s{C_{\rm{b}}}{Z_{\rm{p}}}}&{ - {Z_{\rm{p}}}\left( {1 - s{M_{\rm{b}}}} \right)}\\ { - s{C_{\rm{b}}}}&{\left( {1 - s{M_{\rm{b}}}} \right)} \end{array}} \right]\left[ {\begin{array}{*{20}{l}} {{\rm{ \mathsf{ δ} }}{P_1}}\\ {{\rm{ \mathsf{ δ} }}{Q_1}} \end{array}} \right] \end{array} $ (4)
 

式中:Zp=Rp+sLp为泵特征阻抗。

1.4 泵间管路燃气射流冷凝模型

泵间管燃气射流在主流液氧中形成离散气泡,假设离散气泡为球形气泡,并不考虑气泡的聚合与破碎对射流冷凝过程的影响,气泡的动态过程由气泡动力学方程Rayleigh-Plesset方程求解[22],其表达式为

$ \begin{array}{*{20}{c}} {R\frac{{{{\rm{d}}^2}R}}{{{\rm{d}}{t^2}}} + \frac{3}{2}{{\left( {\frac{{{\rm{d}}R}}{{{\rm{d}}t}}} \right)}^2} + \frac{{4{\nu _{\rm{L}}}}}{R} \cdot \frac{{{\rm{d}}R}}{{{\rm{d}}t}} = }\\ {\frac{{{p_{{\rm{g}}0}}}}{{{\rho _{\rm{L}}}}}{{\left( {\frac{{{R_0}}}{R}} \right)}^{3\gamma }} - \frac{{2S}}{{{\rho _{\rm{L}}}R}} - \frac{{{p_{\rm{L}}}\left( t \right)}}{{{\rho _{\rm{L}}}}}} \end{array} $ (5)
 

式中:R为气泡半径;νL为液体运动黏度;pg0为气泡初始压力;R0为气泡初始直径;ρL为液体密度;S为表面张力;pL为液体当地压力;γ为气体多变指数。

考虑整个计算域气泡核密度N0,将气泡半径R与气泡总体积Vg联系起来,则有$\frac{{{\rm{d}}{V_{\rm{g}}}}}{{{\rm{d}}t}} = {N_0}4{\rm{ \mathsf{ π} }}{R^2}{V_{\sum {} }}\frac{{{\rm{d}}R}}{{{\rm{d}}t}}$,其中VΣ为泵间管路总体积,代入式(5),通过小偏差线性化得出泵间气泡总体积动力学方程,其表达式为

$ {m_{{\rm{eq}}}}\frac{{{{\rm{d}}^2}\left( {{\rm{ \mathsf{ δ} }}{V_{\rm{g}}}} \right)}}{{{\rm{d}}{t^2}}} + b\frac{{{\rm{d}}\left( {{\rm{ \mathsf{ δ} }}{V_{\rm{g}}}} \right)}}{{{\rm{d}}t}} + \frac{1}{{{C_{\rm{g}}}}}\delta {V_{\rm{g}}} = - \frac{{{\rm{ \mathsf{ δ} }}{p_{\rm{L}}}\left( t \right)}}{{{\rho _{\rm{L}}}}} $ (6)
 

式中:${m_{{\rm{eq}}}} = \sum\limits_{k = 0}^N {\frac{{{\rm{d}}Z}}{{4{\rm{ \mathsf{ π} }}{R_k}{l_x}}}} $为气泡等效惯性质量;$b = \sum\limits_{k = 0}^N {\frac{{{\rm{d}}Z}}{{{l_x}}}\left( {\frac{{{v_{\rm{L}}}}}{{{\rm{ \mathsf{ π} }}{R_k}^3}} - \frac{{\left| {{U_{\rm{g}}} - {U_{\rm{L}}}} \right|}}{{4{\rm{ \mathsf{ π} }}{R_k}^2}}} \right)} $为线性阻尼;Cg为等效气泡柔度,忽略表面张力后${C_{\rm{g}}} = {\rho _{\rm{L}}}\frac{{{V_{\rm{g}}}}}{{\gamma {p_{\rm{g}}}}}$, pg为气泡压力。随着气泡冷凝,流体流动方向的气泡半径Rk逐渐减小直至消失。因此,将泵间管分成N段,通过每一段的特征参数积分计算,获得等效惯性质量和线性阻尼,其中,dZ为每段管路长度,lx为燃气冷凝长度,Ug为气泡速度,UL为液体速度。另外,气泡冷凝过程既不是绝热过程也不是等温过程,而是气泡与液氧传质传热过程,Cg与气泡温度相关,则有${C_{\rm{g}}} = {\rho _{\rm{L}}}\frac{{{V_{{\rm{g0}}}}}}{{\gamma {p_{\rm{g}}}}} \cdot \frac{{{T_{{\rm{g}}k}}}}{{{T_{{\rm{g0}}}}}}$Vg0为气泡初始总体积,Tgk为每一段气泡温度,Tg0为气泡初始温度。

主路流体的连续方程为

$ - {\rho _{\rm{L}}}\frac{{{\rm{d}}\left( {{\rm{ \mathsf{ δ} }}{V_{\rm{g}}}} \right)}}{{{\rm{d}}t}} = {\rm{ \mathsf{ δ} }}{q_2} - {\rm{ \mathsf{ δ} }}{q_1} $ (7)
 

式中:δq1和δq2分别表示泵间管入口和出口的质量流量变化率。联立式(6)与式(7),通过拉普拉斯变换将泵间液体流动连续方程写成复数域传递函数,并考虑液氧主流受气泡挤压的流动变形,得出泵间管传递函数为

$ \begin{array}{l} \left[ {\begin{array}{*{20}{c}} {{\rm{ \mathsf{ δ} }}{P_2}}\\ {{\rm{ \mathsf{ δ} }}{Q_2}} \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} 1&0\\ { - \frac{s}{{{s^2}m + sb + \frac{1}{{{C_g}}}}}}&1 \end{array}} \right] \cdot \\ \left[ {\begin{array}{*{20}{c}} {ch\left( {\frac{{s{l_x}}}{a}} \right)}&{ - \frac{{{Z_{\rm{c}}}}}{{1 - \alpha }} \cdot sh\left( {\frac{{s{l_x}}}{a}} \right)}\\ { - \frac{{1 - \alpha }}{{{Z_{\rm{c}}}}}sh\left( {\frac{{s{l_x}}}{a}} \right)}&{ch\left( {\frac{{s{l_x}}}{a}} \right)} \end{array}} \right]\left[ {\begin{array}{*{20}{l}} {{\rm{ \mathsf{ δ} }}{P_1}}\\ {{\rm{ \mathsf{ δ} }}{Q_1}} \end{array}} \right] \end{array} $ (8)
 

式中:α为空泡份额,表示气液两流体通过管路截面的体积比。

为了准确获得泵间管传递函数各特征参数,需要获得燃气射流气泡沿流动方向冷凝参数分布特性。燃气以高速射入主路液氧中冷凝,采用均相流模型无法模拟射流燃气与液氧因速度差和温度差而形成的质量、动量和能量传递过程。因此,本文中采用热不平衡两流体六方程模型[11]计算燃气射流冷凝过程的分布特性。简化计算射流装置结构,认为燃气从射流孔沿轴向射入主路液氧,两相流体在管路中为一维流动,并作如下假设:①两相存在清晰的相界面,质量、动量和能量的传递在相界面上完成,相内部参数均匀;②不考虑气泡的聚合与破碎作用,相界面按球形表面积确定,相界面假设为存在不计厚度的饱和液体薄膜,能量和动量传递按球面对流过程计算;③气氧为过热蒸汽,根据气氧温度确定气体状态方程,液氧为过冷液体,液体物性参数由当地压力和温度关系获得。忽略两相间的黏性应力以及黏性耗散。

热不平衡两流体模型为

气相质量方程:

$ {\rho _{\rm{g}}}{U_{\rm{g}}}\frac{{{\rm{d}}\alpha }}{{{\rm{d}}Z}} + {\rho _{\rm{g}}}\alpha \frac{{{\rm{d}}{U_{\rm{g}}}}}{{{\rm{d}}Z}} + \alpha {U_{\rm{g}}}\frac{{{\rm{d}}{\rho _{\rm{g}}}}}{{{\rm{d}}Z}} + {\rho _{\rm{g}}}\alpha {U_{\rm{g}}}\frac{1}{A} \cdot \frac{{{\rm{d}}A}}{{{\rm{d}}Z}} = \mathit{\Gamma } $ (9)
 

液相质量方程:

$ \begin{array}{l} - {\rho _{\rm{L}}}{U_{\rm{L}}}\frac{{{\rm{d}}\alpha }}{{{\rm{d}}Z}} + {\rho _{\rm{L}}}\left( {1 - \alpha } \right)\frac{{{\rm{d}}{U_{\rm{L}}}}}{{{\rm{d}}Z}} + \left( {1 - \alpha } \right){U_{\rm{L}}}\frac{{{\rm{d}}{\rho _{\rm{L}}}}}{{{\rm{d}}Z}} + \\ \;\;\;{\rho _{\rm{L}}}\left( {1 - \alpha } \right){U_{\rm{L}}}\frac{1}{A} \cdot \frac{{{\rm{d}}A}}{{{\rm{d}}Z}} = - \mathit{\Gamma } \end{array} $ (10)
 

气相动量方程:

$ \begin{array}{l} {\rho _{\rm{g}}}{U_{\rm{g}}}\alpha \frac{{{\rm{d}}{U_{\rm{g}}}}}{{{\rm{d}}Z}} + \alpha \frac{{{\rm{d}}p}}{{{\rm{d}}Z}} = - {F_{{\rm{wg}}}} - {F_{{\rm{vm}}}} - {F_{\rm{D}}} + {F_{{\rm{gi}}}} - \\ \;\;\;\;\;\alpha {\rho _{\rm{g}}}g\sin \theta \end{array} $ (11)
 

液相动量方程:

$ \begin{array}{l} {\rho _{\rm{L}}}{U_{\rm{L}}}\left( {1 - \alpha } \right)\frac{{{\rm{d}}{U_{\rm{L}}}}}{{{\rm{d}}Z}} + \left( {1 - \alpha } \right)\frac{{{\rm{d}}p}}{{{\rm{d}}Z}} = - {F_{{\rm{wL}}}} + {F_{{\rm{vm}}}} + \\ \;\;\;\;\;\;\;{F_{\rm{D}}} + {F_{{\rm{Li}}}} - \left( {1 - \alpha } \right){\rho _{\rm{L}}}g\sin \theta \end{array} $ (12)
 

总能量方程:

$ \begin{array}{l} \left( {{h_{\rm{g}}} - {h_{\rm{L}}} + \frac{{U_{\rm{g}}^2 - U_{\rm{L}}^2}}{2}} \right)\frac{{{\rm{d}}x}}{{{\rm{d}}Z}} + x\left( {\frac{{{\rm{d}}{h_{\rm{g}}}}}{{{\rm{d}}Z}} + {U_{\rm{g}}}\frac{{{\rm{d}}{U_{\rm{g}}}}}{{{\rm{d}}Z}}} \right) + \\ \;\;\;\;\;\left( {1 - x} \right)\left( {\frac{{{\rm{d}}{h_{\rm{L}}}}}{{{\rm{d}}Z}} + {U_{\rm{L}}}\frac{{{\rm{d}}{U_{\rm{L}}}}}{{{\rm{d}}Z}}} \right) + g\sin \theta = 0 \end{array} $ (13)
 

气相能量方程:

$ x\left( {\frac{{{\rm{d}}{h_{\rm{g}}}}}{{{\rm{d}}Z}} + {U_{\rm{g}}}\frac{{{\rm{d}}{U_{\rm{g}}}}}{{{\rm{d}}Z}}} \right) + {h_{{\rm{gL}}}}\frac{{{\rm{d}}x}}{{{\rm{d}}Z}} + xg\sin \theta = \frac{{{q_{\rm{i}}}}}{G} $ (14)
 

式中:ρ为密度;p为压力;x为干度;U为速度;A为管路截面积;Γ为相界面传质率;h为流体焓值;g为重力加速度;θ为流体运动方向与重力加速度方向的夹角;hgL为汽化潜热;qi为界面换热率;G为两流体总质量流速;Fvm为虚拟质量力,表示两相速度差对应的加速度引起的虚拟力;FD为界面拖曳力,表示气泡在流体中受到的流体阻力;FLiFgi为动量交换,表示气体冷凝导致的质量变化产生的动量变换;FwLFwg为壁面摩擦力,由于气泡位于流体中,只考虑液体作用;下标g表示气体,L表示液体,w表示壁面,i表示相界面。

源项模型如表 2所示,其中,CD(1-α)为拖曳力系数;dg为气泡直径;Cvm为质量力系数;η为分布指数;f为摩擦阻力系数;D为管路直径;φLO2为摩擦乘子;Re为雷诺数;kL为液体导热系数;Nu为努塞尔数;Pr普朗特数;μL为液体动力黏度;HL为液体换热系数;qL为单位面积液体换热率;qg为单位面积气体换热率。通过仿真计算获得燃气射流冷凝过程沿主路的特征参数分布特性,得出等效气泡直径、空泡份额、气体速度和温度的沿程分布特性。由这些参数得出特征气泡半径Rk、气泡总体积Vg、气泡当地温度Tgk以给定泵间管传递函数式(8)中的气泡惯性质量、线性阻尼以及气泡柔度等状态参数。

表 2 源项模型 Table 2 Source term model
项目 模型
界面拖曳力 ${F_{\rm{D}}}\frac{3}{4} \cdot \frac{{{C_{D\left( {1 - \alpha } \right)}}}}{{{d_{\rm{g}}}}}\alpha {\left( {1 - \alpha } \right)^3}{\rho _{\rm{L}}}\left( {{U_{\rm{g}}} - {U_{\rm{L}}}} \right)\left| {{U_{\rm{g}}} - {U_{\rm{L}}}} \right|$
虚拟质量力 ${F_{{\rm{vm}}}} = {C_{{\rm{vm}}}}{\rho _{\rm{L}}}{U_{\rm{g}}}\alpha \frac{{{\rm{d}}\left( {{U_{\rm{g}}} - {U_{\rm{L}}}} \right)}}{{{{\rm{d}}_Z}}}$
动量交换 ${F_{{\rm{Li}}}} = \left( {1 - \eta } \right)\left( {{U_{\rm{L}}} - {U_{\rm{g}}}} \right)\Gamma , {F_{gi}} = \eta \left( {{U_{\rm{g}}} - {U_{\rm{L}}}} \right)\Gamma $
壁面摩擦力 ${F_{{\rm{wL}}}} = \varphi _{{\rm{LO}}}^2{\left( {\frac{{{\rm{d}}p}}{{{\rm{d}}Z}}} \right)_{\rm{L}}}$, φLO2按Chisholm[23]方法确定, ${\left( {\frac{{{\rm{d}}p}}{{{\rm{d}}Z}}} \right)_{\rm{L}}} = \frac{f}{D} \cdot \frac{{G{{\left( {1 - x} \right)}^2}}}{{2{\rho _{\rm{L}}}}}$
相界面传热系数 ${H_{\rm{L}}} = \frac{{{k_{\rm{L}}}Nu}}{{{d_{\rm{g}}}}}, {\mathop{\rm Re}\nolimits} = \frac{{{\rho _{\rm{L}}}\left| {{U_{\rm{g}}} - {U_{\rm{L}}}} \right|{d_{\rm{g}}}}}{{{\mu _{\rm{L}}}}}$
$Nu = \left\{ {\begin{array}{*{20}{c}} \begin{array}{l} 2 + 0.6{{\mathop{\rm Re}\nolimits} ^{0.5}}{\Pr ^{0.33}}\;\;\;0 < {\mathop{\rm Re}\nolimits} \le 776.06, \\ \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;0 < \Pr \le 250 \end{array}\\ \begin{array}{l} 2 + 0.27{{\mathop{\rm Re}\nolimits} ^{0.62}}{\Pr ^{0.33}}\;\;\;\;\;\;\;{\mathop{\rm Re}\nolimits} > 776.06, \\ \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;0 < \Pr \le 250 \end{array} \end{array}} \right.$
相界面传质方程 $\Gamma = G\frac{{{\rm{d}}x}}{{{\rm{d}}Z}} = \frac{{{q_{\rm{L}}} - {q_{\rm{g}}}}}{{{h_{{\rm{gL}}}}}}$
2 泵间管路燃气射流冷凝过程的数值仿真

结合表 1工况参数,研究不同液氧入口压力和液氧温度对泵间管路燃气射流冷凝过程的影响,获得额定工况条件下特征参数沿管路的分布特性。将沿程空泡份额积分可以求出泵间管路内停留的气体总体积Vg,从而计算出等效气泡柔度Cg,通过对气泡特征半径Rk积分得出气泡群等效惯性。

2.1 额定工况燃气射流冷凝特征参数分布特性

发动机在额定工况工作时,发动机液氧入口压力为0.4 MPa,液氧温度为90 K,经过氧化剂预压泵的增压后,泵间管入口压力增加到1.35 MPa,液氧温度升高至92 K。以发动机额定工况条件进行仿真计算,图 3~图 6为燃气射流在额定工况条件下特征参数沿管路的分布特性,包括空泡份额α、燃气流速Ug和液氧流速UL、燃气温度Tg和液氧温度TL以及气泡半径Rk。受流体压力和气体温度影响,使得气泡密度较小同时射流速度远大于流体速度,燃气射流受到界面拖曳力较大,气体速度迅速降低,并表现为气泡在射流孔下游0.15 m位置迅速堆积,使得空泡份额迅速增大至峰值。此时气体换热面积大、冷凝效果强,气泡半径开始迅速减小,使得空泡份额开始变小。在0.3 m后由于气体温度降低至饱和温度,流速也与液氧接近,换热效果变差,气泡半径缩小速度变慢。气泡继续随流体流动,最后当气泡半径小于临界值时,由于惯性作用气泡溃灭消失,气体全部冷凝。

图 3 额定工况下空泡份额沿管路的变化 Fig. 3 Variation of void fraction with pipeline at rated condition
图 4 额定工况下气氧速度和液氧速度沿管路的变化 Fig. 4 Variation of gas oxygen velocity and liquid oxygen velocity with pipeline at rated condition
图 5 额定工况下气氧温度和液氧温度沿管路的变化 Fig. 5 Variation of gas oxygen temperature and liquid oxygen temperature with pipeline at rated condition
图 6 额定工况下气泡半径沿管路的变化 Fig. 6 Variation of bubble radius with pipeline at rated condition

通过仿真计算得出燃气与液氧等各参数沿流动方向变化大,气液界面发生强烈的质量、动量以及能量交换。气泡在其温度降低至当地压力下的饱和温度后仍需要较长行程才能完全冷凝,在这段行程中(0.3 m至完全冷凝),气泡温度和速度都已接近液体,使得气液换热效率低、冷凝效果差,且此时空泡份额仍然较大,主路液氧的压力或流量的波动都可能造成气泡冷凝状态发生较大变化形成不稳定振荡源,导致整个系统发生振荡。

2.2 边界条件对燃气射流冷凝过程的影响

比较不同泵间管入口压力和液氧温度对燃气射流冷凝过程的影响。图 7展示了入口压力和液氧温度对燃气冷凝长度的影响。随着入口压力升高,燃气密度升高,射流速度降低,对应的饱和温度也升高,冷凝长度不断增加,入口压力2.1 MPa对应的气氧冷凝长度比1.2 MPa对应的冷凝长度长0.15 m。这说明随着压力升高,气泡变小导致表面积缩小,引起气体的换热面积较低压时小,则其换热能力变差,导致了高压下冷凝长度较低压更长。仿真结果表明在不同入口压力条件下,燃气都能在泵间管内完全冷凝,确保进入主泵时为单相液氧不会引起气蚀现象。对于液氧温度的影响,由于液氧过冷度与对流换热率直接相关,温度越低换热效率越高,则燃气冷凝速率越快。

图 7 入口压力和液氧温度对燃气冷凝长度的影响 Fig. 7 Effect of inlet pressure and liquid oxygen temperature on gas condensation length

图 8展示了入口压力和液氧温度对气泡柔度Cg的影响。Cg是关于气泡总体积和压力的函数,是影响泵间管传递函数的重要参数。燃气冷凝速率越快,气泡柔度就越低。仿真结果表明高入口压力和过冷液氧都使得气泡柔度减小,且压力为1.95 MPa(对应贮箱压力为1.1 MPa)和温度为82 K的过冷液氧条件下,气泡柔度Cg相当。

图 8 入口压力和液氧温度对气泡柔度的影响 Fig. 8 Effect of inlet pressure and liquid oxygen temperature on bubble flexibility
3 泵间管路燃气射流冷凝对氧路系统频率特性的影响

基于串联系统流体网络理论[1]将发动机氧路传递矩阵关系式联立求解,并设定贮箱为系统开端,泵出口为系统闭端,考虑泵间管路燃气射流冷凝过程对氧路系统影响,气泡冷凝过程存在的扰动形式在系统出口以流量扰动的形式对整个氧路系统产生影响,因此通过求解氧路系统各位置参数对出口流量扰动响应特性得出氧路系统频率特性。

图 9展示了在额定工况下试车台与发动机组成的氧路输送系统频率特性仿真曲线,其一阶频率和二阶频率分别为3 Hz和8.3 Hz,仿真结果与试车数据分频值一致,发动机在额定工况工作时氧路系统稳定。图 10展示了在高入口压力(1.2 MPa)下氧路输送系统频率特性仿真曲线与试车压力脉动测点数据分频对比,一阶频率略微升高,而二阶频率从8.3 Hz提高至11 Hz,说明高入口压力对泵间管燃气射流冷凝过程产生了较大影响,使得系统频率升高。另外,二阶频率响应幅值迅速增大表明氧路系统在这一频率下发生谐振,稳定性降低。这一特征与发动机在高入口压力条件工作氧路系统发生频率为11 Hz的低频振荡现象一致。图 11展示了过冷液氧(82 K)条件下氧路输送系统频率特性仿真曲线与试车压力脉动测点数据分频对比,同样二阶频率从8.3 Hz提高至11 Hz,说明过冷液氧对泵间管燃气射流冷凝过程也产生了较大的影响,使得频率升高。不同的是,其二阶频率响应幅值并没有增大很多。在过冷液氧试车过程中也出现了11 Hz的低频振荡频率,但其幅值要远小于高入口压力所对应频率的振荡幅值。高入口压力和过冷液氧条件下发动机氧路系统频率表现出相同二阶频率11 Hz,但高入口压力下系统表现的响应幅值要远远高于过冷液氧下的幅值,这一现象也在热试车过程中得到验证。

图 9 额定工况下氧路系统的频率特性 Fig. 9 Frequency characteristic of oxygen feed system at rated condition
图 10 高入口压力条件下氧路系统的频率特性 Fig. 10 Frequency characteristic of oxygen feed system at high inlet pressure condition
图 11 过冷液氧条件下氧路系统的频率特性 Fig. 11 Frequency characteristic of oxygen feed system at subcooled liquid oxygen condition
4 结论

基于两相流热不平衡两流体模型研究了泵间管燃气射流冷凝过程,建立了泵间管传递函数矩阵,通过建立氧路输送管路系统传递函数模型研究了不同边界条件下氧路系统频率特性。通过仿真计算与分析得出了以下结论:

1) 基于气泡动力学方程建立了泵间管传递函数模型,并通过两流体模型方法获得两相流体各状态参数沿流动方向的分布特性,由此得出的等效气泡惯性质量、等效气泡柔度等参数提高了泵间管传递函数模型的准确性。

2) 气泡在其温度降低至当地压力下的饱和温度后仍需要较长行程才能完全冷凝,在这段行程中气液换热效率低且空泡份额大,受主路液氧的压力或流量变化影响可能造成气泡冷凝状态发生较大变化形成不稳定振荡源,导致整个系统发生振荡。

3) 高入口压力和过冷液氧条件下,氧系统的二阶频率从8.3 Hz提高至11 Hz,与泵间管燃气射流冷凝过程密切相关。高入口压力下二阶频率响应幅值迅速增大表明氧路系统在这一频率下发生谐振,稳定性降低。仿真结果与热试车结果一致。

参考文献
[1] 李斌, 杜大华, 张贵田, 等. 液氧/煤油补燃发动机低频频率特性研究[J]. 航空动力学报, 2009, 24(5): 1187-1191.
LI B, DU D H, ZHANG G T, et al. Research on the low frequency characteristics of LOX/kerosene staged combustion cycle engine[J]. Journal of Aerospace Power, 2009, 24(5): 1187-1191. (in Chinese)
Cited By in Cnki (7) | Click to display the text
[2] 张贵田. 高压补燃循环发动机[M]. 北京: 国防工业出版社, 2005: 250-251.
ZHANG G T. High pressure staged combustion LOX/kerosene rocket engine[M]. Beijing: National Defense Industry Press, 2005: 250-251. (in Chinese)
[3] KERNEY P J, FAETH G M, OLSON D R. Penetration characteristics of a submerged steam jet[J]. AIChE Journal, 1972, 18(5): 548-553.
Click to display the text
[4] WEIMER J C, FAETH G M, OLSON D R. Penetration of vapor jets submerged in subcooled liquids[J]. AIChE Journal, 1973, 19(3): 552-558.
Click to display the text
[5] WU X Z, YAN J J, SHAO S F, et al. Experimental study on the condensation of supersonic steam jet submerged in quiescent subcooled water:Steam plume shape and heat transfer[J]. International Journal of Multiphase Flow, 2007, 33(12): 1296-1307.
Click to display the text
[6] 武心壮, 邱斌斌, 种道彤, 等. 单喷嘴蒸汽射流凝结引起的压力振荡研究[J]. 西安交通大学学报, 2014, 48(1): 48-52.
WU X Z, QIU B B, CHONG D T, et al. Pressure oscillation induced by steam jet condensation in water through a nozzle[J]. Journal of Xi'an Jiaotong University, 2014, 48(1): 48-52. (in Chinese)
Cited By in Cnki (5) | Click to display the text
[7] 武心壮, 严俊杰, 潘冬冬, 等. 过膨胀超音速蒸汽射流的一种流形及其换热研究[J]. 工程热物理学报, 2010, 31(12): 2031-2034.
WU X Z, YAN J J, PAN D D, et al. Experimental study of over-expanded supersonic steam jet in subcooled water[J]. Journal of Engineering Thermophysics, 2010, 31(12): 2031-2034. (in Chinese)
Cited By in Cnki (2) | Click to display the text
[8] PILIPENKO V V, DOROSH I L, MANKO I K.气氧束流进入液氧主流时的气体凝结实验[J].李向阳, 译.研制动态, 2014, 335(8): 1-3.
PILIPENKO V V, DOROSH I L, MANKO I K. Experimental research on vapor condensation of oxygen vapor inject into LOX mainstream[J]. LI X Y, translated. Development Information, 2014, 335(8): 1-3(in Chinese).
[9] LI S Q, WANG P, LU T. Numerical simulation of direct contact condensation of subsonic steam injected in a water pool using VOF method and LES turbulence model[J]. Progress in Nuclear Energy, 2015, 78: 201-215.
Click to display the text
[10] SHAH A, CHUGHTAI I R, INAYAT M H. Numerical simulation of direct-contact condensation from a supersonic steam jet in subcooled water[J]. Chinese Journal of Chemical Engineering, 2010, 18(4): 577-587.
Click to display the text
[11] 宋纪元, 陈听宽. 汽液两相临界流动的热力学非平衡两流体模型[J]. 核科学与工程, 1997, 17(3): 193-201.
SONG J Y, CHEN T K. A non-equilibrium two-fluid model for prediction two-phase critical flow[J]. Chinese Journal of Nuclear Science and Engineering, 1997, 17(3): 193-201. (in Chinese)
Cited By in Cnki (4) | Click to display the text
[12] 张小英, 丁斐, 陈佳跃. 反应堆一维两流体模型二阶精度数值解法研究[J]. 核动力工程, 2013, 34(4): 27-32.
ZHANG X Y, DING F, CHEN J Y. Study on second-order scheme for one-dimensional two-fluid model of reactor[J]. Nuclear Power Engineering, 2013, 34(4): 27-32. (in Chinese)
Cited By in Cnki (2) | Click to display the text
[13] 李伟, 苏光辉, 秋穗正, 等. 一种稳定性增强及高精度数值方法在RELAP5中的实现与评价[J]. 核技术, 2016, 39(11): 110601.
LI W, SU G H, QIU S Z, et al. Implementation and assessment of a stability-enhancing and high-resolution numerical scheme in RELAP5[J]. Nuclear Techniques, 2016, 39(11): 110601. (in Chinese)
Cited By in Cnki | Click to display the text
[14] 巢飞, 单建强, 张勇, 等. 两流体双压力模型半隐数值算法研究[J]. 核动力工程, 2018, 39(2): 142-148.
CHAO F, SHAN J Q, ZHANG Y, et al. Study on semi-implicit scheme for two-fluid seven-equation two-pressure model[J]. Nuclear Power Engineering, 2018, 39(2): 142-148. (in Chinese)
Cited By in Cnki | Click to display the text
[15] 陈二锋, 厉彦忠, 应媛媛. 泵间管气液两相流压力波传播速度数值研究[J]. 航空动力学报, 2010, 25(4): 754-760.
CHEN E F, LI Y Z, YING Y Y. Numerical investigation on pressure wave propagation speed of gas-liquid two-phase flow in pump pipeline[J]. Journal of Aerospace Power, 2010, 25(4): 754-760. (in Chinese)
Cited By in Cnki (21) | Click to display the text
[16] 薛帅杰, 刘红军, 洪流, 等. 厚液膜敞口型离心喷嘴自激振荡特性试验[J]. 航空学报, 2018, 39(9): 122189.
XUE S J, LIU H J, HONG L, et al. Test on self-excited oscillation characteristics of an open-end swirl injector with thick liquid film[J]. Acta Aeronautica et Astronautica Sinica, 2018, 39(9): 122189. (in Chinese)
Cited By in Cnki | Click to display the text
[17] 张国渊, 党佳琦, 赵伟刚, 等. 高速水润滑机械密封的两相流热振动现象[J]. 航空学报, 2019, 40(3): 422532.
ZHANG G Y, DANG J Q, ZHAO W G, et al. Two-phase flow thermal vibration phenomenon of high-speed water-lubricated mechanical seal[J]. Acta Aeronautica et Astronautica Sinica, 2019, 40(3): 422532. (in Chinese)
Cited By in Cnki | Click to display the text
[18] 邢理想, 杜大华, 李斌. 液氧/煤油补燃循环发动机氧路低频动特性分析[J]. 火箭推进, 2009, 35(5): 24-28.
XING L X, DU D H, LI B. Low frequency characteristics analysis of LOX/kerosene staged combustion cycle rocket engine[J]. Journal of Rocker Propulsion, 2009, 35(5): 24-28. (in Chinese)
Cited By in Cnki | Click to display the text
[19] 刘上, 刘红军, 陈宏玉. 液体火箭发动机热力组件动力学模型[J]. 宇航学报, 2012, 33(10): 1512-1518.
LIU S, LIU H J, CHEN H Y. Dynamics models for the combustor component in liquid rocket engine[J]. Journal of Astronautics, 2012, 33(10): 1512-1518. (in Chinese)
Cited By in Cnki (3) | Click to display the text
[20] 张青松, 张兵. 大型液体运载火箭POGO动力学模型研究[J]. 中国科学:技术科学, 2014, 44(5): 525-531.
ZHANG Q S, ZHANG B. POGO dynamic model research for liquid launch vehicles[J]. Scientia Sinica(Technologica), 2014, 44(5): 525-531. (in Chinese)
Cited By in Cnki (1) | Click to display the text
[21] 朱平平, 潘辉, 黄辉, 等. 氧输送管路动态水试仿真分析[J]. 导弹与航天运载技术, 2015(1): 66-68.
ZHU P P, PAN H, HUANG H, et al. Simulating analysis of liquid oxygen pipe dynamic test with water[J]. Missiles and Space Vehicles, 2015(1): 66-68. (in Chinese)
Cited By in Cnki | Click to display the text
[22] 卢义玉, 葛兆龙, 李晓红. 空化空泡发育和溃灭过程的数值分析[J]. 中国矿业大学学报, 2009, 38(4): 582-585.
LU Y Y, GE Z L, LI X H. Numerical analysis on growth and collapse of cavitation bubble[J]. Journal of China University of Mining & Technology, 2009, 38(4): 582-585. (in Chinese)
Cited By in Cnki (10) | Click to display the text
[23] 徐济鋆. 沸腾传热和气液两相流[M]. 北京: 原子能出版社, 2000: 120-121.
XU J Y. Boiling heat transfer and gas-liquid two-phase flow[M]. Beijing: Atomic Energy Press, 2000: 120-121. (in Chinese)
http://dx.doi.org/10.7527/S1000-6893.2019.23393
中国航空学会和北京航空航天大学主办。
0

文章信息

张淼, 李斌, 邢理想
ZHANG Miao, LI Bin, XING Lixiang
基于燃气射流冷凝的氧路系统频率特性研究
Study on frequency characteristics of oxygen feed system based on gas jet and condensation
航空学报, 2020, 41(2): 123393.
Acta Aeronautica et Astronautica Sinica, 2020, 41(2): 123393.
http://dx.doi.org/10.7527/S1000-6893.2019.23393

文章历史

收稿日期: 2019-08-20
退修日期: 2019-09-25
录用日期: 2019-10-12
网络出版时间: 2019-10-18 13:48

相关文章

工作空间