滑行过程中火箭发动机泵系统各构件温度数值模拟

利用数值计算的方法，研究了某型液体火箭发动饥泵系统各构件及泵腔残余氧化剂在二次启动前滑行过程的温度变化。基于集总参数法建立了温度变化计算模型，并建立了泵腔残余氧化剂排空计算模型。计算结果与试车数据吻合良好。滑行过程中．涡轮温度下降，而泵壳体及进口管的温度升高；泵腔中残余氧化剂间歇排放。研究结果表明，滑行结束时泵系统的温度高于100℃，泵腔中残余氧化剂含汽率为1，不满足二次点火启动条件。     

液氢液氧发动机喷管化学反应流数值模拟

用时间相关的半隐格式求解液氢液氧发动机喷管中的化学反应流动，考虑了６种燃烧产物，８个基元反应，得到了流动参数和产物的质量分数在流场中的分布，为发动机化学动力学损失计算和喷管热结构设计提供了必要数据。计算结果表明，采用半隐格式解决组分连续方程刚性引起的数值不稳定是有效的。     

在火箭发动机燃烧室中气相湍流流动的理论计算

本文用单流体湍流模型对液体火箭发动机突扩燃烧室中冷态气相流动的速度场和湍流强度分布进行了数值计算并与实验数据进行了对比和分析.     

塞锥侧板对直排塞式喷管性能的影响

为了解塞锥侧板对直排塞式喷管性能的影响,提出了对应着同一个塞锥的不同形状和高度的侧板。采用数值模拟的方法,得到了带有不同侧板的直排塞式喷管在不同工况下的流场和性能。比较了不同侧板形状和高度对直排塞式喷管性能的影响。结果表明:带有曲边侧板的塞式喷管的性能略高于带有直边侧板的塞式喷管,但不超过1%;没有侧板时塞式喷管性能要下降4%～6%;塞锥侧板的高度应不小于内喷管出口的高度,低于内喷管出口高度的侧板会降低塞式喷管性能;当侧板高于内喷管出口时,直排塞式喷管性能不随侧板高度的增加而增大。最后定性分析了带有不同侧板塞锥的流场结构,找出了产生这种结果的原因。     

双钟型喷管高度补偿特性的数值分析

针对双钟型喷管流场中存在的激波、激波与附面层相互作用、大尺度分离流动等复杂的物理现象，应用RNGκ-ε湍流模型封闭二维粘性可压缩N-S方程组，采用二阶迎风格式及隐式Gauss-Seidel迭代方法进行耦合求解，对双钟型喷管流场进行数值模拟研究。结果表明，双钟型喷管在低空出现壁面可控流动分离，避免了过膨胀现象，性能优于大面积比钟型喷管；在高空下为完全附着流动，性能优于小面积比钟型喷管，实现了高度补偿特性。     

旋转对内冷通道气膜孔流量系数的影响

为了获得涡轮叶片旋转时对气膜孔流量系数的影响,采用非结构化网格及标准k-ε紊流模型,求解三维N-S方程,对带90°肋和气膜孔出流的旋转矩形通道内的三维流场进行了数值模拟,气动参数变化范围是:通道入口雷诺数Re=60000,150000,罗斯贝数Ro=0,0.11,0.22,气膜孔总出流比为0.22和0.09。Ro≠0时,旋转效应对气膜孔流量系数有明显影响,通道顺时针旋转时,科氏力由上壁面指向带气膜孔的下壁面,引起流量系数增加;通道逆时针旋转时,情况相反。计算结果还表明,在通道内同一径向位置处的两个气膜孔的流量系数是不同的,通道顺时针旋转时,进入左侧孔的流体和气膜孔轴线的夹角小于进入右侧孔的流体和气膜孔轴线的夹角,导致左侧孔的流量系数大于右侧孔的流量系数;通道逆时针旋转时,情况相反。Ro=0时的计算结果与实验数据符合很好。     

电动气阀动态特性及反力因素的影响

建立了电动气阀吸动过程的数学模型并进行了动态特性的数值计算，详细分析了反力因素对电动气阀动态特性的影响。计算结果表明：该阀的衔铁在被触动后能很快地吸到最小气隙处；适当地减小压缩气的入口压力和弹簧预紧力有利于提高电动气阀的响应能力，有利于延长电动气阀的使用寿命。     

空腔倾向对泄漏影响的流动分析

给出了前倾斜腔，矩形空腔和后倾斜腔三种结构的四空腔迷密封实验结果。通过数值求解非定常Ｎ－Ｓ方程，预测上述三种结构的单空腔迷宫通道的内部流动，计算结果表明：后倾斜腔结构内部是稳态流动，而对于前斜腔和矩形空腔结构，流动表现为总体结构的基本稳定和局部结构的不稳定，局部不稳定的强弱顺序依次是前倾斜腔，矩形腔和后倾斜腔的迷宫结构。通过分析局部结构不稳定对减小泄漏的作用机制，说明不稳定效应的强弱是造成密封效果     

直接的矩阵分裂计算跨音速喷管流场

给出了分离系数矩阵求解欧拉方程组时系数矩阵的直接分离公式，减少和简化了应用该方法时的矩阵运算，参照Ｂｅａｍ－Ｗａｒｍｉｎｇ显式格式构造与ＳＣＭ方法一致的内点差分格式，并分析了它的数值特性，对有斜激波的跨音速内流场进行了数值仿真。计算中采用变系数的当地时间步长，加快了流场达到的收敛速度，计算无需引入任何人工参数，计算值与实验结果吻合。     

高超音速化学非平衡钝劈绕流数值计算

 自1970年Davis提出粘性激波层方法以来,用数值方法求解高超音速轴对称钝体绕流问题国内外已做了大量工作,但至今仍未见到关于平面问题的计算结果。对于平面问题,虽然方程在形式上比轴对称简单,但由于二维效应,激波层较厚,用文献[2]的方法向下游区推进有困难。另外,在驻点线上采用极限关系式虽能克服方程的奇性,但驻点解对流向步长Δξ有依赖。