-
摘要:
针对激波诱导推力矢量控制研究仅限于主流和次流气体为同种气体的研究现状,开展了不同次流气体分子质量对推力矢量性能的影响规律研究。首先采用二阶精度AUSM+格式和
k -ω SST两方程湍流模型求解三维Favre平均Navier-Stokes方程,模拟了喷管复杂干扰内流场。然后计算了He、N2和CO2等次流气体在不同注气角度、注气压力和主流落压比下的矢量偏转角度和推力系数。计算结果表明:平均分子质量越小的次流气体矢量偏转角度越大,推力损失越小。因此可选用平均分子质量小的气体作为次流气源,或者将从燃烧室引出的高温燃气与分子质量小的气体混合。Abstract:Aimed at the research status of the shock thrust vector control, which is limited to the mainstream and secondary flow gas as the same gas, the influence of different secondary flow gas molecular mass on thrust vector performance is investigated. First, a turbulence model described by two equations (AUSM+ scheme and
k-ω SST) at two-order accuracy was utilized to solve the Favre averaged three-dimensional Navier-Stokes equations, which simulated the complex interference inner flow field of the nozzle, and the vectorial deflection angles and thrust coefficients were calculated under different gas injection angles, injection pressures and nozzle pressure ratios when the secondary flow gas of He, N2 and CO2 were selected. The calculation results show that the smaller the mean molecular mass of the secondary flow gas is, the larger the vectorial deflection angle is, the less the thrust loss is. Therefore, the gas with smaller mean molecular mass could be used as the gas source of the secondary flow, or the high temperature gas derived from the combustion chamber could be mixed with the gas with a smaller mean molecular mass. -
表 1 空气流和氮气流物性参数
Table 1. Physical property parameters for air stream and nitrogen stream
参数 空气流 氮气流 Ma 3.71 1.0 p/MPa 1.236 0.606 T/ K 301 301 YN2 0.766 4 1 YO2 0.233 6 0 注:Ma—马赫数;YN2—氮气组分质量分数;YO2—氧气组分质量分数。 表 2 不同注气角度下矢量喷管性能参数对比
Table 2. Comparison of vector nozzle performance parameter under different gas injection angles
气体 θ/(°) α/(°) CR, i ω He 90 4.91 0.912 0.026 105 5.08 0.906 0.026 N2 90 4.89 0.911 0.035 105 4.96 0.905 0.035 CO2 90 4.80 0.910 0.044 105 4.88 0.904 0.044 表 3 不同NPR下矢量喷管性能参数对比
Table 3. Comparison of vector nozzle performance parameter under different NPRs
气体 NPR α/(°) CR, i He 10 5.08 0.912 15 4.73 0.886 N2 10 4.96 0.911 15 4.66 0.885 CO2 10 4.88 0.910 15 4.61 0.884 表 4 不同SPR下矢量喷管性能参数对比
Table 4. Comparison of vector nozzle performance parameter under different SPRs
气体 SPR α/(°) CR, i He 0.8 4.09 0.917 1.0 5.08 0.912 N2 0.8 3.97 0.917 1.0 4.96 0.911 CO2 0.8 3.93 0.917 1.0 4.88 0.910 -
[1] DENG R Y, SETOGUCHI T, KIM H D.Large eddy simulation of shock vector control using bypass flow passage[J].International Journal of Heat and Fluid Flow, 2016, 62:474-481. doi: 10.1016/j.ijheatfluidflow.2016.08.011 [2] SHI J W, ZHOU L, WANG Z X, et al.Investigation on flowfield characteristics and performance of shock vector control nozzle based on confined transverse injection[J].Journal of Engineering for Gas Turbines and Power-Transactions of the ASME, 2016, 138(10):101502. doi: 10.1115/1.4033140 [3] 刘辉, 邢玉明.燃气喷射推力矢量喷管气固两相流数值模拟[J].航空动力学报, 2013, 28(1):151-157. http://d.old.wanfangdata.com.cn/Periodical/hkdlxb201301021LIU H, XING Y M.Numerical simulation of gas-particle flow in hot gas injection thrust vector nozzle[J].Journal of Aerospace Power, 2013, 28(1):151-157(in Chinese). http://d.old.wanfangdata.com.cn/Periodical/hkdlxb201301021 [4] 马伟, 杜刚, 金捷, 等.激波诱导矢量喷管动态数值模拟[J].燃气涡轮试验与研究, 2014, 27(5):30-37. doi: 10.3969/j.issn.1672-2620.2014.05.006MA W, DU G, JIN J, et al.Dynamic numerical simulation of shock vector control exhaust nozzle[J].Gas Turbine Experiment and Research, 2014, 27(5):30-37(in Chinese). doi: 10.3969/j.issn.1672-2620.2014.05.006 [5] DEERE K A.Summary of fluidic thrust vectoring research conducted at NASA Langley Research Center: AIAA-2003-3800[R].Reston: AIAA, 2003. [6] WAITHE K A, DEERE K A.Experimental and computational investigation of multiple injection ports in a convergent-divergent nozzle for fluidic thrust vectoring: AIAA-2003-3802[R].Reston: AIAA, 2003. [7] DEERE K A, BERRIER B L.Computational study of fluidic thrust vectoring using separation control in a nozzle: AIAA-2003-3803[R].Reston: AIAA, 2003. [8] 林飞, 王根斌.固体火箭发动机推力向量控制[M].北京:国防工业出版社, 1981:121-125.LIN F, WANG G B.Thrust vectoring control for solid rocket motor[M].Beijing:National Defense and Industry Press, 1981:121-125(in Chinese). [9] FAVRE A.Equations des gaz turbulents compressibles[J].Journal de Mechanique, 1965, 4:361-390. [10] MENTER F R.Two equation eddy viscosity turbulence models for engineering applications[J].AIAA Journal, 1994, 32(8):1598-1605. doi: 10.2514/3.12149 [11] LIOU M S.Progress towards an improved CFD method: AUSM+: AIAA-95-1701[R].Reston: AIAA, 1995. [12] VAN LEER B.Towards the ultimate conservative difference scheme V:A second order sequel to Godunov's method[J].Journal of Computational Physics, 1979, 32(1):101-136. http://www.sciencedirect.com/science/article/pii/0021999179901451 [13] YOON S, JAMESON A.Lower-upper symmetric Gauss-Seidel method for the Euler and Navier-Stoker equations[J].AIAA Journal, 1988, 26(9):1025-1026. doi: 10.2514/3.10007 [14] ASO S, OKUYAMA S, KAWAT M, et al.Experimental study on mixing phenomena in supersonic flows with slot injection: AIAA-91-0016[R].Reston: AIAA, 1991. [15] 刘辉, 邢玉明, 额日其太.气体二次喷射矢量喷管三维流场计算[J].北京航空航天大学学报2009, 35(10):1174-1178. http://bhxb.buaa.edu.cn/CN/abstract/abstract8720.shtmlLIU H, XING Y M, ERIQITAI.Computation of three-dimensional flow field in secondary gas injection for thrust vectoring nozzle[J].Journal of Beijing University of Aeronautics and Astronautics, 2009, 35(10):1174-1178(in Chinese). http://bhxb.buaa.edu.cn/CN/abstract/abstract8720.shtml [16] DEERE K A.Computational investigation of the aerodynamic effects on fluidic thrust vectoring: AIAA-2000-3598[R].Reston: AIAA, 2000.