导弹与航天运载技术（中英文）

不同设计参数对液氧甲烷发动机喷管效率的影响

PDF下载

刘阳旻 , 田原

导弹与航天运载技术 | 动力系统 2024,47(1): 57-63

收起

导弹与航天运载技术 | 动力系统 2024, 47(1): 57-63

不同设计参数对液氧甲烷发动机喷管效率的影响

全屏

刘阳旻, 田原

作者信息

北京航天动力研究所,北京,100076

刘阳旻(1998—),男,工程师,主要研究方向为液体火箭发动机喷管性能仿真计算。

田原(1979—),男,研究员,主要研究方向为大推力液体火箭发动机推力室设计。

Influence of Different Design Parameters on Nozzle Efficiency of Liquid Oxygen Methane Engine

Yangmin LIU, Yuan TIAN

Affiliations

Beijing Aerospace Propulsion Institute,Beijing,100076

出版时间: 2024-02-25 doi: 10.7654/j.issn.2097-1974.20240111

文章导航

摘要

收起

为了更好更快地计算喷管性能,采用特征线法加附面层修正的方式开发了一种二维喷管化学动力流场计算程序,针对某液氧甲烷发动机计算了不同设计参数对应的喷管流场,研究了不同设计参数对喷管效率的影响,设计参数包括推力、室压、喷管面积比和喷管出口角等。计算喷管效率时考虑了几何损失、附面层损失、化学动力损失。结果表明使用程序计算能够准确快速获得工程上重视的设计指标,大大提高了仿真效率,可用于喷管型面优化。

关键词

液氧甲烷发动机 / 数值仿真 / 喷管效率 / 特征线法 / 附面层修正

Abstract

收起

In order to better and faster calculate nozzle performance, a program for calculating the chemical dynamic flow field of two dimensional nozzles is developed by using the method of characteristic and boundary layer correction. The nozzle field corresponding to different for a liquid oxymethane engine, and the influence of different design parameters on nozzle efficiency are studied, the design parameters include thrust, chamber pressure, nozzle area ratio and nozzle outlet angle, etc. Geometric loss, boundary layer loss and chemical dynamic loss are considered in the calculation of nozzle efficiency. The results show that the program can obtain the number important for engineer quickly and accurately, improve the simulation efficiency greatly, and can be used to optimize the nozzle contour.

Key words

liquid oxygen methane engine / numerical simulation / nozzle efficiencies / method of characteristic / boundary layer correction

引用本文

刘阳旻, 田原. 不同设计参数对液氧甲烷发动机喷管效率的影响. 导弹与航天运载技术, 2024 , 47 (1) : 57 -63 . DOI: 10.7654/j.issn.2097-1974.20240111

Yangmin LIU, Yuan TIAN. Influence of Different Design Parameters on Nozzle Efficiency of Liquid Oxygen Methane Engine[J]. Missiles and Space Vehicles, 2024 , 47 (1) : 57 -63 . DOI: 10.7654/j.issn.2097-1974.20240111

正文

收起

0 引言

收起

火箭发动机中喷管的作用是将燃气的热能转化为动能, 加速后的气流对喷管产生反作用推力。燃气在喷管内流动会有各种损失, 其中以几何损失、附面层损失、化学动力损失为主。喷管效率为相同混合比、相同室压条件下的实际真空比冲与理想真空比冲之比。提高喷管效率是喷管设计的主题, 曹彬彬等[1]使用CFD技术研究了不同设计参数对固液火箭发动机的喷管效率的影响；王长辉等^[2]使用CFD技术计算了不同设计对塞式喷管性能的影响；陈林泉等^[3]使用类似于美国SPP^[4]软件的方法研究了喷管扩张段型面对固体发动机的喷管效率的影响；Manski等^[5]使用TDK^[6]研究了不同设计参数对氢氧发动机喷管效率的影响。而液氧甲烷发动机凭借着其优秀的经济性和可复用性得到了国内外广泛关注, 但目前中国尚未针对液氧甲烷发动机设计参数对喷管效率影响进行过系统研究, 故本文采用一种自研的高效二维喷管化学动力流场计算程序研究了液氧甲烷发动机不同设计参数对喷管效率的影响。

本文首先对推力室中各种损失的定义进行了介绍, 随后简要介绍了自研程序的基本工作流程和方法原理, 并给出针对SSME喷管的计算结果, 与文献[5]和[7]使用TDK计算的结果进行了对比, 结果基本一致, 在此基础上对某液氧甲烷发动机喷管流场进行了仿真, 研究了不同设计参数对喷管效率的影响, 并且其计算速度极快, 能够有效提高工程仿真效率。

1 火箭喷管效率定义

收起

本文将喷管损失作为火箭设计参数的函数来进行研究。化学动力损失、几何损失及附面层损失^[8]是研究的主题。

为了衡量化学动力损失, 化学动力效率被定义为

(1)${\eta }_{{I}_{{S}_{\text{kin }}}}= {I}_{{S}_{\text{ODK }}}/{I}_{{S}_{\text{ODE }}}$

比冲${I}_{{S}_{\text{ODE }}}$由一维平衡分析计算得出,比冲${I}_{{S}_{\text{ODK }}}$是由一维动力学流动分析计算得出。

为了衡量几何损失,几何效率${\eta }_{{I}_{{s}_{\mathrm{{dw}}}}}$被定义为

(2)${\eta }_{{I}_{{S}_{\text{div }}}}= {I}_{{S}_{2\mathrm{D}}}/{I}_{{S}_{\text{ODK }}}$

比冲${I}_{{S}_{20}}$使用包括化学动力效应的二维流动分析计算得出。

为了衡量附面层损失,附面层效率${\eta }_{{I}_{{S}_{\min }}}$被定义为

(3)${\eta }_{{I}_{{S}_{\text{fric }}}}= {I}_{{S}_{\mathrm{{BL}}}}/{I}_{{S}_{2\mathrm{D}}}$

比冲${I}_{{S}_{\mathrm{{BL}}}}$使用包括化学动力效应、附面层及传热的二维流动分析计算得出。

喷管效率${\eta }_{{I}_{{s}_{\max }}}$由下式得到:

(4)${\eta }_{{I}_{{S}_{\max }}}= {\eta }_{{I}_{{S}_{\text{kin }}}}{\eta }_{{I}_{{S}_{\text{div }}}}{\eta }_{{I}_{{S}_{\text{div }}}}$

随着外界气压的不同, 计算的比冲值也会不同, 而本文均使用真空比冲进行计算分析, 而上述比冲的具体公式可参考文献[6]。

2 计算方法

收起

本文使用自研程序进行计算, 它是一款无黏流场结合附面层修正的二维喷管化学动力流场求解器, 其中无黏主流部分使用特征线法^[9]计算考虑化学反应的喷管流场, 化学反应以有限速率模型建模, 附面层修正部分采用的是T.Cebeci^[10]的二维/轴对称可压缩附面层方程求解方法, 并且附面层计算也考虑了化学效应, 程序结构和计算框架示意见图1和图2。

程序基本流程是先进行一维平衡计算再进行一维非平衡计算, 一维平衡计算主要使用的是最小吉布斯自由能法,其原理可参考美国CEA^[11]软件,一维平衡计算得到的轴向压力分布和收敛段入口的燃气组份及物性参数将作为一维非平衡计算的输入, 一维非平衡计算的方程为

(5)$\frac{\mathrm{d}{c}_{i}}{\mathrm{\;d}x}= \frac{{\omega }_{i}{r}^{* }}{\rho V}$

(6)$\frac{\mathrm{d}V}{\mathrm{\;d}x}= -\frac{1}{\rho V}\frac{\mathrm{d}P}{\mathrm{\;d}x}$

(7)$\frac{\mathrm{d}\rho }{\mathrm{d}x}= \left\lbrack {\frac{1}{\gamma P}\frac{\mathrm{d}P}{\mathrm{\;d}x}- A}\right\rbrack \rho $

(8)$\frac{\mathrm{d}T}{\mathrm{\;d}x}= \left\lbrack {\frac{\gamma - 1}{\gamma }\frac{1}{P}\frac{\mathrm{d}P}{\mathrm{\;d}x}- B}\right\rbrack T $

马赫数大于1.04时, 采用下述方程:

(9)$\frac{\mathrm{d}V}{\mathrm{\;d}x}= \left\lbrack {\frac{1}{\bar{a}}\frac{\mathrm{d}\bar{a}}{\mathrm{\;d}x}- A}\right\rbrack \frac{V}{{M}^{2}-1}$

(10)$\frac{\mathrm{d}\rho }{\mathrm{d}x}= -\left\lbrack {\left({\frac{1}{\bar{a}}\frac{\mathrm{d}\bar{a}}{\mathrm{\;d}x}- A}\right)\frac{{M}^{2}}{{M}^{2}-1}+ A}\right\rbrack \rho $

(11)$\frac{\mathrm{d}T}{\mathrm{\;d}x}= -\left\lbrack {\left({\frac{1}{\bar{a}}\frac{\mathrm{d}\bar{a}}{\mathrm{\;d}x}- A}\right)\frac{\left({\gamma - 1}\right){M}^{2}}{{M}^{2}-1}+ B}\right\rbrack T $

(12)$ A =\frac{{r}^{* }}{PV}\left\lbrack {\mathop{\sum }\limits_{{i = 1}}^{N}{\omega }_{i}{R}_{i}T -\frac{\gamma - 1}{\gamma }\mathop{\sum }\limits_{{i = 1}}^{N}{\omega }_{i}{h}_{i}}\right\rbrack $

(13)$ B =\frac{\gamma - 1}{\gamma }\frac{{r}^{* }}{PV}\mathop{\sum }\limits_{{i = 1}}^{N}{\omega }_{i}{h}_{i}$

式中${c}_{i}$为$i$组分的浓度;${\omega }_{i}$为$i$组分的净反应生成速率;${R}_{i}$为$i$组分的气体常数;$x$为喷管轴向距离/喉部半径或者曲线坐标系流向坐标;${r}^{* }$为喉部半径;$V$为总速度；$P$为压强；$\rho$为密度；$T$为温度；$\gamma$为比热比；$\bar{a}$为喷管面积比;$A$为联系流动和化学反应过程的绝热额外项;$B$为联系流动和化学反应过程的能量交换项;${h}_{i}$为$i$组分的静焓;$M$为马赫数或者第三体。

单组分的净生成速率${\omega }_{i}$使用有限速率模型建模求解,采用碳、氢、氧、氮11组元$\left({\mathrm{H}\text{、}{\mathrm{H}}_{2}\text{、}{\mathrm{H}}_{2}\mathrm{O}}\right.$、$\mathrm{{OH}}\text{、}\mathrm{O}\text{、}{\mathrm{O}}_{2}\text{、}\mathrm{{CO}}\text{、}{\mathrm{{CO}}}_{2}\text{、}\mathrm{\;N}\text{、}\mathrm{{NO}}\text{、}{\mathrm{N}}_{2}$),16步化学反应模型^[7]:

${\mathrm{O}}_{2}+ \mathrm{H}\rightleftharpoons \mathrm{O}+ {\mathrm{{OHH}}}_{2}+ \mathrm{O}\rightleftharpoons \mathrm{H}+ {\mathrm{{OHH}}}_{2}+ \mathrm{{OH}}\rightleftharpoons {\mathrm{H}}_{2}\mathrm{O}+ \mathrm{H}\\\mathrm{{OH}}+ \mathrm{{OH}}\rightleftharpoons {\mathrm{H}}_{2}\mathrm{O}+ \mathrm{{OCO}}+ \mathrm{{OH}}\rightleftharpoons {\mathrm{{CO}}}_{2}+ {\mathrm{{HN}}}_{2}+ \mathrm{O}\rightleftharpoons \mathrm{{NO}}+ \mathrm{N}\\{\mathrm{O}}_{2}+ \mathrm{N}\rightleftharpoons \mathrm{{NO}}+ \mathrm{{OCO}}+ \mathrm{O}\rightleftharpoons {\mathrm{{CO}}}_{2}{\mathrm{{CO}}}_{2}+ \mathrm{O}\rightleftharpoons \mathrm{{CO}}+ {\mathrm{O}}_{2}\mathrm{H}+ \\\mathrm{H}+ \mathrm{M}\rightleftharpoons {\mathrm{H}}_{2}+ \mathrm{{MH}}+ \mathrm{{OH}}+ \mathrm{M}\rightleftharpoons {\mathrm{H}}_{2}\mathrm{O}+ \mathrm{{MO}}+ \mathrm{O}+ \\\mathrm{M}\rightleftharpoons {\mathrm{O}}_{2}+ \mathrm{{MN}}+ \mathrm{O}+ \mathrm{M}\rightleftharpoons \mathrm{{NO}}+ \mathrm{{MN}}+ \mathrm{N}+ \mathrm{M}\rightleftharpoons {\mathrm{N}}_{2}+ \mathrm{M}\\\mathrm{{CO}}+ \mathrm{O}+ \mathrm{M}\rightleftharpoons {\mathrm{{CO}}}_{2}+ \mathrm{{MO}}+ \mathrm{H}+ \mathrm{M}\rightleftharpoons \mathrm{{OH}}+ \mathrm{M}$

这里$\mathrm{M}$代表第3体,有限速率模型原理可参考文献[12], 一维非平衡求解方法可参考文献[7]和文献[13]。

喉部跨声速区的控制方程为

(14)$\frac{\partial {v}^{\prime }}{\partial x}- \frac{\partial {u}^{\prime }}{\partial r}= 0 $

(15)$\left({\gamma + 1}\right){u}^{\prime }\frac{\partial {u}^{\prime }}{\partial x}- \frac{\partial {v}^{\prime }}{\partial r}- \frac{{v}^{\prime }}{r}= 0 $

式中${u}^{\prime },{v}^{\prime }$为小扰动量。

利用改良的Sauer方法对其进行求解, 得到特征线法初始线, 求解结果与参考文献[7]一致。

接着进行二维特征线法计算, 由于特征线法只适用于双曲型方程的性质, 故该模块只计算喷管扩张段流场, 采用柱坐标系, 求解的方程沿流线为

(16)$\frac{\mathrm{d}r}{\mathrm{\;d}x}= \tan \theta $

(17)$\mathrm{d}\frac{{V}^{2}}{2}+ \frac{\mathrm{d}P}{\rho }= 0 $

(18)$\frac{\mathrm{d}P}{\gamma P}- \frac{\mathrm{d}\rho }{\rho }= \frac{A}{\cos \theta }\mathrm{d}x $

(19)$\frac{\gamma - 1}{\gamma }\frac{\mathrm{d}P}{P}- \frac{\mathrm{d}T}{T}= \frac{B\mathrm{\;d}x}{\cos \theta }$

(20)$\mathrm{d}{c}_{i}= \frac{{\omega }_{i}}{{\rho V}\cos \theta }\mathrm{d}x $

沿左特征线:

(21)$\frac{\mathrm{d}x}{\mathrm{\;d}r}= \cot \left({\theta +\alpha }\right)$

(22)$\frac{\mathrm{d}P}{P}= G\left\lbrack {\left({A -\frac{\sin \theta }{r}}\right) F\mathrm{\;d}r -\mathrm{d}\theta }\right\rbrack $

沿右特征线:

(23)$\frac{\mathrm{d}r}{\mathrm{\;d}x}= \tan \left({\theta -\alpha }\right)$

(24)$\frac{\mathrm{d}P}{P}= - G\left\lbrack {\left({A -\frac{\sin \theta }{r}}\right) H\mathrm{\;d}x -\mathrm{d}\theta }\right\rbrack $

这里的A和B同前文。

(25)$ F =\cos \theta -\sin \theta \cot \left({\theta +\alpha }\right)$

(26)$ G =\frac{\gamma }{\sin \alpha \cos \alpha }$

(27)$ H =\cos \theta \tan \left({\theta -\alpha }\right)- \sin \theta $

式中$r$为从对称轴开始的径向距离/喉部半径；$\theta$为轴向与径向速度夹角；$\alpha$为马赫角；$H$为总焓。

使用特征线法具体的求解过程可参照文献[9]。

收缩段一维和扩张段二维非平衡计算的结果将作为附面层模块的输入进行附面层修正计算, 附面层模块求解的是曲线坐标系下的轴对称可压缩薄边界层方程, 如下:

(28)$\frac{\partial }{\partial x}\left({{\rho u}{r}^{k}}\right)+ \frac{\partial }{\partial x}\left({\overline{\rho v}{r}^{k}}\right)= 0 $

(29)${\rho u}\frac{\partial u}{\partial x}+ \overline{\rho v}\frac{\partial u}{\partial y}= -\frac{\mathrm{d}p}{\mathrm{\;d}x}+ \frac{1}{{r}^{k}}\frac{\partial }{\partial y}\left\lbrack {{r}^{k}\left({\mu \frac{\partial u}{\partial y}- \rho \overline{{u}^{\prime }{v}^{\prime }}}\right)}\right\rbrack $

(30)${\rho u}\frac{\partial H}{\partial x}+ \overline{\rho v}\frac{\partial H}{\partial y}= \frac{1}{{r}^{k}}\frac{\partial }{\partial y}\left\lbrack {{r}^{k}\left({k\frac{\partial u}{\partial y}- {c}_{p}\rho \overline{{T}^{\prime }{v}^{\prime }}+ u\left({\mu \frac{\partial u}{\partial y}- \rho \overline{{u}^{\prime }{v}^{\prime }}}\right)}\right)}\right\rbrack $

式中$u$为轴向速度或者曲线坐标系流向速度;$v$为径向速度或者曲线坐标系纵向速度;$y$为曲线坐标系纵向坐标;${c}_{p}$为定压比热容;$k$为流动指标数,$k = 0$时为二维流动,$k = 1$时为轴对称流动,雷诺应力使用C-S零方程模型建模, 采用盒式格式差分求解, 具体方法可参考文献[14]。附面层模块没有组分方程但是在物性设置上使用了近似化学平衡的物性插值方法, 具体方法可参见文献[7]。

附面层模块只运行一次,其输出的位移厚度随后被向内叠加到原型面上, 得到新型面, 并将计算得到的再生冷却热增加到推进剂焓中(初始推进剂焓采用未吸收再生冷却热的值), 然后再进行一次无黏计算, 这便是本文所用程序的全部流程(不考虑传热时, 可不进行二次无黏计算)。

图3和图4是使用自研程序对SSME比冲、喷管效率、位移厚度分布及动量厚度分布计算的结果。图3中不同颜色的散点对应文献[5]中的结果。图4中实线为程序计算结果, 散点为文献结果。

由图3和图4可知, 本文所使用的的自研程序在SSME喷管上的计算结果与Manski^[5]使用TDK/BLM对SSME喷管比冲与喷管效率计算的结果基本一致, 数值上存在误差主要是由于程序差异, 并且低室压时误差有所增加, 原因可能为Manski使用厚附面层方程组进行附面层修正计算, 低室压时喷管附面层厚度增加, 本文使用的附面层修正模块所采用的薄附面层假设适用性降低, 但这并不影响后面对甲烷发动机进行性能分析, 同时SSME标准工况下自研程序计算得到的位移厚度与动量厚度分布与Nickerson^[6]使用TDK/MABL得到的结果几乎完全一致, 因此总的来说自研程序的计算结果与文献符合很好, 程序正确性得到了验证。

为了研究设计参数对甲烷发动机喷管效率的影响, 本文选择了室压、喉部半径、面积比和喷管出口角度作为自变量。本程序在喷管型面构建上使用参数化建模的方法, 只需输入喉部半径、收敛段半径、上下游圆弧段半径、面积比、喷管出口角、扩张段入口角、收缩比等关键几何特征便可自动构建出喷管型线, 收敛段使用双圆弧构造, 扩张段为下圆弧接上拟抛物线型面, 其中拟抛物线公式如下:

(31)$ y = d +\sqrt{c +{bx}}+ {ax}$

计算中室压${P}_{\mathrm{C}}$依次取$2\mathrm{{MPa}}$、$3\mathrm{{MPa}}$、$4\mathrm{{MPa}}$、$5\mathrm{{MPa}}$、$6\mathrm{{MPa}}$、$7\mathrm{{MPa}}$、$8\mathrm{{MPa}}$、$9\mathrm{{MPa}}$、${10}\mathrm{{MPa}}$、 11 MPa、11.5 MPa、13.5 MPa、20 MPa、30 MPa；喉部半径${r}_{\mathrm{t}}$取${55}\mathrm{\;{mm}}$、${110.1}\mathrm{\;{mm}}$、${311}\mathrm{\;{mm}}$,改变喉部半径是在几何相似的情况下, 相当于将喷管整体进行放大或者缩小,见图5。面积比$\varepsilon$取${11.7}\text{、}{20}\text{、}{33}$,改变面积比而不改变喉部半径只是放大了喷管的扩张段, 见图6。出口角度$\alpha$取${2}^{\circ }$、${3}^{\circ }$、${4}^{\circ }$、${6}^{\circ }$、${8}^{\circ }$、${8.436}^{\circ }$(该值是一个喷管性能较优的值)、${9}^{\circ }$、${10}^{\circ }$、${12}^{\circ }$、${14}^{\circ }$。由于该型甲烷发动机拟采用全再生冷却设计, 可认为燃气没有向壁面传热, 故壁面边界条件取绝热, 附面层计算使用全湍流计算, 所有计算都没有进行附面层修正型面后的二次无粘计算。

3 计算结果

收起

3.1 各类效率随室压变化的趋势

图7为不同的效率随室压的变化, 其中混合比为3.5,喉部半径为${110}\mathrm{\;{mm}}$,面积比为20,出口角度为${8.436}^{\circ }$。从图7可以看出随着室压增加化学动力效率增加。这可解释为由于室压增加,密度增加,导致化学平衡向重组反应倾斜, 使得燃烧室中离解反应程度降低, 喷管中的化学动力损失降低, 因此化学动力效率增加。但从图7中可知在其他条件不变时, 改变室压对附面层效率和几何效率影响不大。

3.2 化学动力效率随室压、喉部半径和面积比的变化

图8显示了化学动力效率随室压、喉部半径、喷管面积比变化而变化的趋势, 其中混合比为3.5, 出口角度为${8.436}^{\circ }$。从图8中可看出,喷管喉部半径减少, 绝对长度减少, 导致流体在喷管中的停留时间减少, 这意味着离解产物发生重组的时间减少, 使得化学动力损失增大, 化学动力效率降低。而更高的面积比使得离解产物经历的冻结流动时间增加, 化学能向动能的转化减少, 故化学动力效率降低。但总的来说, 喉部半径和面积比变化对化学动力效率的影响均不及前文讨论过的室压, 而且从图8中可看出喷管尺寸越小, 室压对化学动力效率的影响越大。

3.3 几何效率随室压、喉部半径和面积比的变化

图9为几何效率随室压、喉部半径、喷管面积比的变化,其中混合比为3.5,出口角度为${8.436}^{\circ }$。由图9可知几何效率随着室压增加而降低, 这是因为室压越高,出口气流速度越大,出口流场愈不均匀。在室压达到${14.5}\mathrm{{MPa}}$之前,喉部半径增加使得几何效率增加,超过${14.5}\mathrm{{MPa}}$后,喉部半径增加使得几何效率降低。面积比增加,喷管长度增加,同样导致气流速度增大, 出口流场更易不均匀, 使得几何效率降低。

3.4 附面层效率随室压、喉部半径和面积比的变化

附面层效率随喷管设计参数变化如图10所示, 其中混合比为3.5,出口角度为${8.436}^{\circ }$。更高的室压产生了更小的位移厚度 (如图11所示, 其中混合比为3.5,出口角度为${8.436}^{\circ }$),因此提高了附面层效率。更大的面积比意味着更长的壁面,这会增加附面层损失,因此附面层效率随着面积比增大而减小。更大的喉部半径增加了喷管长度, 位移厚度增大, 但是相对位移厚度减小 (如图11所示), 附面层损失减少, 所以最终附面层效率随着喉部半径增大而增大。

3.5 喷管效率随室压、喉部半径和面积比的变化

总的喷管效率变化规律如图12所示, 其中混合比为3.5,出口角度为${8.436}^{\circ }$,从前面的讨论可看出, 对于该型甲烷发动机, 随着室压增加, 附面层效率和化学动力效率增加, 而几何效率变化不大, 故室压增加,喷管效率增加。同样的, 喉部半径增加, 附面层效率和化学动力效率增加, 几何效率也变化不大, 故喉部半径增大, 喷管效率增大。而随着面积比增大, 附面层效率、化学动力效率及几何效率都减小, 故面积比越大,喷管效率越低。

3.6 喷管出口角度变化对各类效率的影响

图13展示了各类喷管效率随壁面出口角度的变化,其中混合比为3.5,室压为${11.5}\mathrm{{MPa}}$,喉部半径为${110}\mathrm{\;{mm}}$。对于该型甲烷发动机来说,标准出口角度为${8.436}^{\circ }$,不同的出口角度意味着抛物线型面的小变化(如图14所示,其中无量纲坐标为实际单位坐标与喉部半径的比值)。化学动力效率和附面层效率随喷管出口角度变化较小,几何效率在${8}^{\circ }$左右达到最大值,之后开始减少,喷管效率在${8.436}^{\circ }$左右达到最大值。

4 结论

收起

本文用化学动力效率、几何效率、附面层效率以及总喷管效率来衡量喷管性能, 采用特征法加附面层修正的方式进行计算, 得到了以下结论:

a)增加室压、喉部半径或降低面积比有利于增加化学动力效率。降低面积比、室压有利于增加几何效率, 喷管半径对几何效率的影响不是单调的。增加室压、喉部半径,降低面积比有利于增加附面层效率。而总的喷管效率随室压、喉部半径的增加而增加, 随面积比的增加而降低。

b)改变出口角度, 对化学动力效率和附面层效率影响较小, 主要影响几何效率, 并且存在一个最佳角度。

c)采用本程序,所有工况的计算时间均在${60}\mathrm{\;s}$左右, 直接输出关键结果, 无需再处理, 相比现在的CFD技术使用的全N-S方程计算,效率大大提高,可用于优化型面。

参考文献

收起

文献

收起

参考文献引证文献

排序方式：

[1]

曹彬彬, 蔡国飙, 田辉, 等. 固液火箭发动机矢量控制与喷管效率研究[J]. 宇航学报, 2017,38(10): 1124-1130.

CAO

Binbin

, CAI

Guobiao

, TIAN

Hui

, et al. Study on thrust vector control and nozzle efficiency for hybrid rocket motor[J]. Journal of Astronautics, 2017,38(10): 1124-1130.

[2]

王长辉, 刘宇. 主要结构参数对直排塞式喷管性能的影响[J]. 推进技术, 2004(5): 439-443.

WANG

Changhui

, LIU

. Effects of primary configurative parameters on the performance of linear aerospike nozzles[J]. Journal of Propulsion Technology, 2004(5): 439-443.

[3]

陈林泉, 李岩芳, 王建儒, 等. 喷管扩散段型面对固体发动机性能的影响[J]. 固体火箭技术, 2004(1): 9-11.

CHEN

Linquan

, LI

Yanfang

, WANG

Jianru

, et al. Effect of the nozzle divergent cone contour on solid rocket motor performance[J]. Journal of Solid Rocket Technology, 2004(1): 9-11.

[4]

Solid performance program (SPP)[R]. 1987.

[5]

MANSKI

, HAGEMANN

. Influence of rocket design parameters on en-gine nozzle efficiencies[J]. Journal of Propulsion and Power, 1996,12(1): 41-47.

[6]

Two-dimensional reference computer program[R]. 1985.

[7]

Additional support for the TDK/MABL computer program[R]. 1993.

[8]

刘国球. 液体火箭发动机原理[M]. 北京: 中国宇航出版社, 2005.

LIU

Guoqiu

. Principle of liquid rocket engine[M]. Beijing: Chinese Aerospace Publishing House, 2005.

[9]

ZUCROW

M J

, HOFFMAN

J D

. Gas dynamics[M]. New York: Wiley, 1977.

[10]

BRADSHAW

, CEBECI

, WHITELAW

L H

. Engineering calculation methods for turbulent flows[M]. London: Academic, 1981.

[11]

Computer program for calculation of complex chemical equilibrium compositions and applications[R]. 1994.

[12]

JOHN

, ANDERSON

. Hypersonic and high-temperature gas dynamics[M]. Reston: AIAA, 2006.

[13]

李宗杰, 于康元. 火箭发动机化学非平衡流的组合计算方法[J]. 北京联合大学学报, 1995(1): 25-30.

Zongjie

, YU

Kangyuan

. The continent method for nonequilibrium flow of chemical reaction inside rocket nozzle[J]. Journal of Beijing Union University, 1995(1): 25-30.

[14]

CEBECI

, BRADSHAW

. Physical and computational aspects of convective heat transfer[M]. New York: Sping Verlag, 1994.

2024年第47卷第1期

PDF下载

453

195

引用本文

BibTeX

文章信息

doi: 10.7654/j.issn.2097-1974.20240111

接收时间：2020-08-27
首发时间：2025-07-04
出版时间：2024-02-25

补充材料

相关文章

文章信息

作者

出版历史

收稿日期：2020-08-27
修回日期：2021-05-25

基金

作者信息

北京航天动力研究所,北京,100076

参考文献

分享链接

https://castjournals.cast.org.cn/joweb/ddyht/CN/10.7654/j.issn.2097-1974.20240111

分享至

全文二维码

扫描看全文

引用本文

BibTeX

本文的引用情况

2种不同金属材料的力学参数

科 Family	属数 Number of genus	种数 Number of species	占总种数比例 Percentage of total species (%)	属 Genus	种数 Number of species	占总种数比例 Percentage of total species (%)
鹅膏菌科Amanitaceae	2	11	5.26	鹅膏菌属 Amanita	10	4.78
小菇科 Mycenaceae	2	12	5.74	丝盖伞属 Inocybe	5	2.39
多孔菌科 Polyporaceae	8	14	6.70	蜡蘑属 Laccaria	5	2.39
红菇科 Russulaceae	3	23	11.00	小皮伞属 Marasmius	6	2.87
				小菇属 Mycena	11	5.26
				光柄菇属 Pluteus	5	2.39
				红菇属 Russula	17	8.13
				栓菌属 Trametes	5	2.39

关闭全屏

BibTeX
EndNote
RefWorks
TxT