Sajben跨声速扩压器分离流动中湍流模式数值研究

为研究湍流模式对激波/湍流边界层干扰内流流动的影响,提高数值计算准确度,使用SA,SST k-ω,非线性EASM k-ω,Gao-Yong四个湍流模式对Sajben扩压器内激波/湍流边界层干扰流动进行了数值计算。对流项采用Roe格式离散,扩散项采用二阶中心格式离散,离散后的控制方程用多步Runge-Kutta显示时间推进法求解。文中展示了四个湍流模式计算得到的壁面压力、速度剖面、摩阻系数等分布。计算值与实验值符合很好,四个湍流模式总体上能够较好地模拟扩压器内激波/湍流边界层干扰复杂分离流动。Gao-Yong湍流模式对分离区内的压力、速度型的模拟更加准确,而非线性EASM k-ω模式对分离再附点位置计算最理想。     

绕翼型低速湍流的数值模拟

采用求解低速流动的SIMPLE算法,用标准K-ε两方程模型、非线性K-ε两方程模型和一种改进的K-ε两方程模型来对绕两种翼型（A翼型和GAW-1翼型）在接近失速攻角情况下的低速湍流流动进行了数值模拟。计算结果表明,非线性K-ε两方程模型和改进的K-ε两方程模型较好地模拟出了翼型表面上的分离流动特性,较准确地预测出了分离点位置以及翼型的气动力系数,对翼型与机翼的工程设计具有一定的参考价值。     

二维喷管分离流线性/非线性湍流模式的研究

采用三个线性和七个非线性涡粘性湍流模式，对二维非对称喷管中的不可压缩分离流进行了研究。通过将数值计算结果和实验结果进行比较，对有关的湍流模式进行了评估和分析。计算结果表明，非线性模式的模化系数与平均流动应变不变量以及旋转不变量有关，反映了湍流的各向异性，比线性模式优越得多。另外，所有线性／非线性涡粘性湍流模式不能够捕捉流动的非定常分离过程。     

使用GAO-YONG湍流方程组计算翼型分离流

采用SIMPLE方法求解GAO-YONG不可压湍流方程组, 对不同来流迎角下的NACA0012翼型绕流结构进行了数值模拟, 给出了翼型绕流分离流结构随迎角的变化特征和翼型在分离绕流中的气动力参数.与实验数据以及大涡模拟结果的比较表明, GAO-YONG不可压湍流方程组能够对翼型绕流的分离点、分离涡形态、表面压力分布、升阻特性做出较好的预测, 能够模拟翼型大攻角分离流动, 计算结果优于FLU-ENT软件中的k-ε、k-ω、k-ωsst模型的计算结果.      

翼型动态失速的数值研究

用不可压缩流动的求解算法，结合WilcoxDC提出的k－ω模式和k－ωSST湍流模式，对翼型的动态失速进行了数值模拟。通过对典型的振荡翼型轻失速和深失速算例的计算结果分析可以看出：（1）绕动态失速翼型的流场结构十分复杂，轻失速和深失速在流动特性上有很大区别。计算结果显示：轻失速主要是由于后缘分离引起，分离涡的影响范围主要是在后缘附近。而深失速则首先形成很大的前缘分离涡，该分离涡在翼型表面上运动，并诱发出二次分离涡，引起翼型升、阻力系数的显著变化。（2）对于动态失速的翼型绕流，k－ωSST湍流模式是较为有效的，计算出的气动力系数迟带曲线变化趋势与实验结果符合得比较好。     

使用GAO-YONG湍流方程组对翼型绕流的计算

将GAO-YONG湍流方程组应用到翼型复杂绕流场的计算中,研究了机械能方程中Du/Dt项和各向异性湍流黏性系数分量中的实度系数Cs的影响,给出了控制Du/Dt项作用区域和选择实度系数Cs值的判别函数,对NACA4412翼型最大升力迎角下的低速绕流算例、AS240翼型8°迎角下的转捩算例、AS240翼型19°迎角下的分离流算例进行了模拟,与实验结果比较,得到了比较好的结果.结果表明GAO-YONG湍流方程组能够准确预测翼型前缘层流到湍流的转捩问题,翼型尾迹流动松弛问题和翼型大迎角分离流.     

采用RANS/LES混合方法研究分离流动

采用结合湍流模式和大涡模拟(LES)的混合方法数值模拟Aerospatial A-profile翼型以及细长椭球体大迎角分离流动。该类混合方法比较适于模拟分离流动:近壁区附体流动采用湍流模式,而远离物面分离区域内的大尺度运动则用LES方法。通过计算和实验结果的对比,混合方法切实可行并具有很大的发展潜力。     

可压缩粘性流绕翼型的N─S方程自适应网格解

从Ｎ－Ｓ方程出发,采用ＬＵ－ＡＤＩ隐式分解方法求解绕翼型的可压缩粘性流动。湍流模式采用Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模式和ｑ－ω二方程微分模式以研究湍流边界层的非平衡效应;为更好地捕捉激波,对网格进行了自适应调整。通过对ＮＡＣＡ００１２和ＲＡＥ２８２２两种翼型在亚音速和跨音速不同马赫数和雷诺数下进行大量计算,表明该方法对改进湍流计算、提高激波分辨率有较好的效果     

不可压缩实际粘性流体绕流孤立翼型的数值研究

本文提出了基于微分-积分方程组求解N-S方程的有限差分法求解不可压缩实际粘性流体绕孤立翼型流动。计算中采用了K-ε湍流模型结合考虑逆压梯度的壁面函数法,并导出任意非正交曲线坐标系下的压力修正方法和适用于大流通量的混合格式。分别计算了层流和湍流情况下带分离翼型绕流,与流动显示及实验数据比较表明,在预测分离点位置,翼型表面压力、平均速度分布、雷诺切应力uv方面,两者吻合,结果是满意的。     

一种基于湍流模式的转捩预测方法

通过对平板边界层进行计算分析,基于湍流模式来预测流动转捩是可行的,但预测出的转捩位置不准,比实际情况靠前。因此,本文从转捩流动的物理特征出发对湍流模式中的混合尺度作出了修正,提出了一种考虑间歇函数影响的k-ωSST两方程湍流模式。采用该模式来对平板湍流边界层和不同攻角下的NACA0012翼型绕流进行了转捩预测,计算结果表明,改进后的模式对转捩位置具有较好的预测能力。     

塞式喷管气动特性与性能研究

对6单元瓦状塞式喷管、6单元圆转方内喷管塞式喷管和1单元直排塞式喷管进行了实验和数值模拟研究.实验表明,3个塞式喷管具有良好的高度补偿特性和较高的效率,6单元瓦状塞式喷管和1单元直排塞式喷管在设计高度的效率接近100%;由于内喷管型面不够理想和加工时存在的缺陷,6单元圆转方内喷管塞式喷管的设计点效率在95%左右.数值模拟结果与实验数据吻合较好,结合实验和数值模拟分析了外界反压对气流,进而对底部气动特性的影响.实验研究了底部二次流流量变化对性能的影响.     

塞式喷管底部气动特性及其对性能的影响

在内喷管倾角为10°,20°,30°,40°和底部二次流为0.0%,0.4%,1.0%,1.4%,2.0%,3.0%的工况下,实验测量了两单元瓦状塞式喷管的底部压强。结果表明在不同的外界反压作用下,底部具有开放、闭合和开闭过渡三种气动状态,底部的这种特性不随内喷管倾角和底部二次流的变化而变化。实验研究了80%,40%,30%和20%截短两单元直排塞式喷管的高度特性,数值模拟计算了实验喷管的性能和塞锥表面的压强分布,数值模拟结果与实验测量数据吻合较好。结合实验数据和数值模拟结果分析了塞锥截短和底部气动特性对塞式喷管性能的影响。      

塞式喷管设计和性能验证

提出了简化的塞式喷管型面设计和优化方法,并以气氢/气氧为推进剂,对圆转方内喷管一单元直排塞式喷管性能进行了热试实验验证.用圆弧和抛物线近似内喷管型面,用抛物线和三次曲线近似塞锥型面,以从海平面到设计高度的飞行总冲最大为目标函数,进行塞式喷管型面设计和优化.介绍了实验系统及实验发动机主要零部件的结构和设计参数,给出了实验参数测量结果、实验照片和数据分析.实验表明,塞式喷管具有良好的高度补偿能力和较高的效率.在三个不同高度下,实验喷管效率在93%～98%之间,预计设计点效率不低于98%.实验结果表明,所提出的塞式喷管型面设计和优化方法是合理可行的.     

高度补偿喷管的氢氧热试研究

采用气氧作氧化剂、气氢作燃料, 对具有高度补偿特性的塞式喷管和双钟型喷管进行了点火热试.介绍了气氢/气氧试验系统, 以及试验喷管的结构形式、设计参数和装配照片, 喷管试验件采用耐烧蚀的钨渗铜材料加工, 成功进行了多次短时间点火热试.给出了试验测量参数曲线、点火热试照片和数据结果, 获得了不同高度下塞式喷管和双钟型喷管的热试性能数据, 和当量钟型喷管相比两者都具有较好的高度补偿特性.以推力系数效率为例, 塞式喷管在低空压强比下达到92%-95%, 双钟型喷管低空下为96%-98%, 高空下则都在95%左右.      

单元塞式喷管结构对性能影响的数值研究

从N－S方程出发，采用二阶精度的NND格式对塞式喷管的流场进行了数值模拟，重点研究了单元喷管的倾角的塞锥的截断长度在偏离设计高度时对性能的影响。结果表明，不管是在低于还是高于设计高度，随着倾角的增加，内喷管对推力的贡献在减小，塞锥和底部对推力的贡献却增大，而侧喷管的最佳倾角却基本相同，都在全长型塞式喷管设计角度附近。截断后的塞式喷管虽然减小了重量，但不管是在低空还是高空，性能均低于全长型塞式喷管，因此增加底部的推力是提高性能的一个重要途径。     

三单元气氢/气氧塞式喷管发动机实验

以气氢/气氧为推进剂, 对三单元直排塞式喷管模型发动机(内喷管和燃烧室内衬采用耐烧蚀材料钨渗铜, 塞锥为碳钢)进行了热试实验研究.给出了实验结果, 进行了数据分析.在两个压比(燃烧室压强与环境反压之比)下, 获得了塞式喷管性能数据和塞锥表面5个测点处的压强分布.实验表明, 塞式喷管具有良好的高度补偿能力和较高的效率.在压比值为50附近, 效率达到92%-93.5%;在压比值为350附近, 效率达到95%-96%.数值模拟的预示结果与实验数据吻合较好.      

瓦"状塞式喷管的实验研究

为了了解塞式喷管的性能特征和更好地设计塞式喷管,本文对两种不同结构的"瓦"状塞式喷管进行了试验研究。实验采用气氧/酒精作为推进剂。实验发动机是一个四单元和一个十四单元的"瓦"状塞式喷管。文中介绍了实验系统和实验发动机的主要结构。给出了典型的试验曲线,并分析了"瓦"状塞式喷管的高度特性和底部二次流对性能的影响情况。结果显示,该种喷管具有很好的高度补偿特性,在底部加入适当的二次流能使喷管性能稍有提高。      

矩形内喷管塞式喷管的数值计算与实验研究

为了了解内喷管为二维矩型的塞式喷管性能,设计了一个二单元的实验塞式喷管,并对模型进行了数值模拟和实验研究。数值模拟采用无波动、无自由函数耗散(NND)差分格式求解三维NS方程,利用空气冷流实验方法评价了喷管性能。研究模型的内喷管喉部面积为4×60mm2,内喷管面积比为4,总面积比为24.05,设计压力比为500。计算得到了正确的流场结构和塞锥表面压强分布,结果与实验数据吻合很好,效率数值最大相差1%。模型的性能也比较理想:最大的推力系数效率为0.995,同钟型喷管相比,具有很好的高度补偿能力:从地面到高空,推力系数效率在0.97～0.995之间变化。不同压强比下全锥塞式喷管的塞锥表面压强分布规律,可以作为研究截短型塞式喷管塞锥压强分布的基础。      

多单元圆转方塞式喷管的性能计算和分析

采用NND差分格式求解三维平均雷诺N-S方程的数值方法, 对多单元圆转方塞式喷管的性能进行了计算.首先比较了钟型喷管、二维喷管、圆转方喷管、方形喉部方形出口喷管的流动特点和性能, 然后研究了转方位置、转方后型面和出口圆角对圆转方内喷管性能的影响以及圆转方结构型式、内喷管倾角、塞锥型面变化对塞式喷管性能的影响, 并且给出了具有较高性能的多单元圆转方塞式喷管设计方案.研究结论对于多单元圆转方塞式喷管的优化设计有一定的参考作用.      

外流对塞式喷管性能的影响

为了了解外流对塞式喷管性能的影响。本文从N-S方程出发，采用NND格式对塞式喷管有外流的流场进行了数值模拟，重点研究了外流马赫数和外流迎角对塞式喷管的影响，研究表明，在一定飞行高度下，外流会使塞式喷管的性能降低，并且外流为亚声速时对塞式喷管的影响较小，在跨声速和超声速时影响较大；当外流有迎角时，会改变塞式喷管推力矢量，正迎角会使塞式喷管性能得到补偿，负迎角则会降低塞式喷管的性能，随着背压减小，外流对塞式喷管的影响也会减小。