高过载条件下绝热层烧蚀实验方法研究 (Ⅰ)方案论证及数值模拟

分析了高过载对固体火箭发动机流场和绝热层烧蚀的影响规律，提出一种新的研究思路：采用数值模拟方法来预示发动机三维两相流场，建立高过载流场模拟实验装置，开展绝热层烧蚀实验，建立高过载条件下的绝热层烧蚀模型，在此基础上发展高过载发动机绝热层设计和烧蚀预示方法。其中关键技术是高过载流场的模拟，对粒子加入法和弯管分离法两种方案进行了论证，认为弯管分离法原理上是可行的。为了验证这种方案的可行性，开展了弯管通道两相流的数值模拟研究，计算结果表明弯管装置具有使凝相粒子聚集形成高浓度粒子流的功能。     

固体粒子在压缩机叶栅内的部分沉积模型研究

针对轴流压缩机叶栅内固体微粒沉积的问题,发展了新的粒子沉积模型(部分沉积模型)。采用此模型和普遍采用的完全沉积模型进行了计算,计算结果表明两种模型对叶栅前缘处的粒子沉积影响明显。采用部分沉积模型时,叶片前缘处不会出现短时间内不合实际的大量沉积现象,而表现为一种比较均匀的粒子沉积。     

湍流参数对复杂形状涡轮盘腔流场的数值研究

以罗罗公司未来某型双轴发动机的高压涡轮旋转盘腔为计算模型，目的是要对复杂的内部旋转空腔流场有一个详细的了解，通过CFD(FLUENT6．0．12)的研究为即将进行的放大比例的旋转涡轮盘腔试验提供指导。计算中采用了Spalart-Allmaras湍流模型，结果表明湍流参数相同的流动，其流体结构是非常相似的；流场主要由湍流参数λT决定涡旋的大小和数量；对于湍流参数小于0．199的旋转流动，盘腔内存在着自由涡流。     

涡轮效率改变对发动机加速特性的影响

采用零维闭环数值仿真平台，对某型单轴涡轮喷气发动机进行启动加速过程的零维闭环仿真，通过选择三种不同的调节方式，研究了涡轮效率的改变对发动机的转速、推力和耗油率等发动机特性的影响，其中某些结果与实验数据进行了比较。     

固体火箭燃气射流近场形成与发展的数值模拟

固体火箭武器燃气射流的近场区域存在有复杂的波系，膨胀波、压缩波的混合作用是机-弹、弹-舰相容性研究的重要内容，复杂波系压力、温度场的分布是燃气射流动力冲击和热冲击作用的直接因素。因为这一区域流动的复杂性，使得实验研究非常困难。为此应用数值分析的方法，研究了固体火箭高度欠膨胀燃气射流流场起始冲击波、近场波系的形成及发展。计算结果揭示了近场区域的一些特性，与文献提供的数据较为吻合，为工程应用提供了帮助。     

H型叶栅通道流计算网格的改进

提出一种壁面正交、尾迹区网格节点分布逐渐均匀化叶栅通道网格生成方法。为了研究壁面处网格正交性对数值计算结果的影响，推导出壁面附近流动参数一阶导数离散误差与网格线交角的数学关系。采用N-S方程分别对平板紊流附面层流和叶栅通道流进行计算分析。结果表明：壁面处网格正交性对附面层内流动计算预测精度影响较大，而对势流区流动影响很小；尾迹区切向网格节点逐渐均匀化分布能更准确计算尾迹区流动。     

二维过渡区喷流干扰流场的DSMC/EPS仿真

应用DSMC／EPS混合算法，数值模拟过渡区有RCS横向喷流干扰的二维有限平板高超声速流场，分析了流场结构、分离形态及物面气动系数分布等特性。结果表明，不同：Knudsen数下，喷流干扰流场存在差异，过渡区中的稀薄气体效应得以体现，证明了DSMC／EPS混合算法是研究稀薄，连续混合区域高超声速流动的有效方法。     

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

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

浮升力对双旋转轴间流动换热影响的数值模拟

利用Boussinesq假设，采用混合长度模型，对具有轴向通流的双旋转轴间的流动换热进行了数值研究，得到了轴间气体充分发展段的速度分布、温度分布和加热边的对流换热系数，其结果在实验模型参数范围内与PfitzerH的实验研究符合较好。在此基础上，进一步数值模拟研究了同向等速旋转双轴间不可压流的离心浮升力对流动换热的影响，结果表明：离心浮升力对流动换热的影响取决于浮升力与惯性力之比GrB/Rex^2。当GrB/Rex^2≥24时，浮升力的影响十分显著，不可忽略。     

论定常不可压粘流的积分方程解法

本文对定常不可压粘性绕流积分解法中的涡量迭代过程、涡量边界条件的提法以及压强系数的计算等方面都作了改进,使数值解更易于稳定收敛并提高精度。为了验证本文的方法,对低雷诺数下的圆柱绕流进行了计算,结果表明计算效率是高的,迭代过程是稳定的。将本文方法推广到三维流动没有原则上的困难。