湍流模型对压气机数值模拟精度的影响

开发了一个求解叶轮机全三维黏性流场的程序,采用了目前叶轮机数值模拟中常用的且模拟性能较好的三种湍流模型:k-ε模型,k-ω模型,Spalart-Allmaras(S-A)模型,计算了跨声压气机NASA(National Aeronautics and Space Administration) Rotor 37的流场,并与实验进行了性能参数、三维流场和角区流动的对比分析.考察在相同的数值计算平台上进行比较这几个湍流模型对压气机的数值模拟性能的影响.最终结果表明:对于具有强剪切、存在分离流的复杂的叶轮机三维流场来说,k-ω模型数值模拟精度更高,相比其他两个模型具有一定的优势.      

立楔诱导的高超声速流动分离数值分析

高超声速飞行器进气道等关键部件引起的激波与边界层相互作用将导致流动分离,从而改变当地压力分布与局部受热情况,影响飞行稳定性与飞行安全,因此需要对高超声速流动的分离现象进行细致研究。采用高精度5阶特征型WENO格式与3阶TVD型Runge-Kutta方法,求解三维Navier-Stokes方程,对立楔诱导的高超声速激波与边界层相互作用引起的分离流动流场结构进行了细致的数值模拟与分析。结果表明,5阶特征型WENO格式分辨率远高于类TVD格式;Ma=6时得到清晰的激波结构、分离涡结构及其演化过程和壁面极限流线的拓扑结构,证明了WENO格式应用于高超声速分离流动的可行性与高分辨率;对不同来流Mach数的对比证明Mach数的增大抑制流动分离,导致分离涡减小。     

基于无网格算法的化学非平衡流数值模拟

为了验证AUSMPW+格式在无网格算法的精度等以及无网格算法能否用于复杂化学非平衡流数值模拟中,首先将AUSMPW+格式推广到无网格算法中,给出了AUSMPW+在无网格算法中的具体形式。其次,将无网格算法应用于带化学反应的多组分气体Euler方程的求解,给出无网格条件下求解化学非平衡流的具体过程。最后,对NACA0012跨声速绕流场、高超声速弹丸绕流场、氢气/氧气诱导燃烧流场、球锥高超声速绕流场、楔体斜爆轰流场进行了数值模拟。数值模拟结果与相关文献计算结果和实验结果吻合较好,表明了AUSMPW+格式在无网格算法中表现较高精度。同时,采用无网格算法能较好地模拟复杂化学非平衡流场,正确分辨复杂的物理现象,为化学非平衡流数值模拟提供一种新的算法,拓宽了无网格算法的应用范围。     

多重网格技术在复杂粘性流场计算中的应用及研究

流体力学控制方程的数值求解过程中,当网格加密、粘性效应强或来流马赫数低的时候,流场收敛非常缓慢。本文在多块网格、LUSGS隐式格式基础上采用多重网格技术,计算了NLR7301case#1两段翼型的低速绕流及DLRF4跨声速流场,在复杂二维、三维粘性流场求解中获得了满意的加速收敛效果。文中还比较了V循环、W循环的加速效率,为多重网格的工程应用奠下基础。     

基于多重网格方法的跨声速颤振数值模拟研究

使用多重网格的方法对二维翼型和三维机翼的跨声速颤振进行了数值求解.流场控制方程为N-S方程,湍流模型采用SA模型网格.计算网格采用结构网格,空间离散采用中心格式,使用双时间法进行时间推进.流场与结构之间的数据通过径向基函数(RBF)插值方法进行交换.通过耦合求解流场控制方程和结构运动方程得到Isogai(case A)机翼模型和Agard445.6机翼的跨声速颤振边界.通过与使用单重网格时得到的结果相比较得出,多重网格方法能够节省计算时间,提高颤振分析计算效率.     

飞机亚跨超绕流的数值模拟研究

本文以Euler方程为数学模型,采用一种高精度的TVD(Total Variation Dimishing)离散格式及一种含近似因式分解的推进迭代方法,求解亚跨超绕流’流场。通过若干算例的试算,证明方法是可行的,可以用来模拟飞机的复杂流场。     

超声速湍流流场的RANS/LES混合计算方法研究

采用对接/拼接网格技术,建立了基于分区混合和基于湍流尺度混合的双重RANS/LES混合计算模型,并对环翼低速绕流、翼型跨声速绕流和球锥带凹窗外形二维超声速绕流进行了初步的数值模拟.环翼和翼型绕流计算表明,该混合模型可给出较合理的湍流宏观平均量;球锥带凹窗外形二维超声速绕流计算表明,该混合模型可得到超声速瞬态湍流脉动流场,凹窗处存在复杂的旋涡结构和波系结构,呈现较大尺度的脉动.但该模型还需要进一步的考核验证.     

一种适用于热流计算的改进WENO格式

数值格式是影响高超声速热流数值计算准确度的主要因素之一.根据高超声速流场的特性合理的控制了数值黏性,对WENO(weighted essentially non-oscillatory)格式进行了改进,使其适用于高超声速热流计算.研究结果表明:改进后的格式在捕捉激波时具有良好的激波辨率,在壁面附近尽量减小了数值黏性,提高了格式的黏性分辨率.数值试验也表明了这种改进算法的有效性.      

基于人工黏性的二维激波小波多尺度数值计算

基于自适应小波配点法和人工黏性，构造了二维激波问题的小波数值计算格式。利用最细尺度的小波系数构造两组激波定位函数，分别用以控制两个正交方向人工黏性的大小和分布。对三类二维激波问题进行计算和验证，结果表明：自适应小波配点法利用小波阈值滤波删除变化平缓区域大量网格点，而保留变化急剧区域的网格点，比传统方法计算效率高，且分辨率越高，计算效率更高；利用最细层小波系数构建的幂函数形式的激波定位函数能准确判断激波位置和控制X、Y方向人工黏性的大小和分布，从而捕捉不同方向和不同强度的激波。     

基于全局气动优化方法的跨声速叶栅气动优化

提出了适用于叶栅三维气动设计优化的全局自动气动优化方法.对NASA Rotor 37转子叶栅进行了气动优化设计.利用该叶栅的试验数据校核了计算流体(CFD)程序的可靠性.以等熵效率最高为目标函数,在满足流量约束和总压比约束的条件下,完成了跨声速叶栅的气动优化设计.优化叶栅的等熵效率提高了1.66%,具有优秀的气动性能和变工况性能.优化结果表明,通过优化三维跨声速叶栅的型线和径向基迭方式,可以有效的减小跨声速叶栅的激波损失.      

跨声速轴流压气机间隙泄漏流触发旋转失速

通过对跨声速轴流压气机NASA转子37进行单通道定常及多通道定常、非定常数值模拟,单通道定常数值模拟结果与实验结果能较好吻合。多通道非定常数值模拟结果显示,间隙泄漏流及其与激波干涉的非定常振荡,触发突尖型旋转失速先兆,具体表现为叶顶前缘间隙泄漏流溢出。失速团首先在叶顶处形成,且速度约为80%转速。随着流量的下降,失速团进一步发展,在失速通道内,激波与叶片前缘完全分离,且在叶片尾缘出现回流。当转子完全数值失速时,失速团周向尺度约为4个通道,且径向占据约半个叶高。     

跨声速压气机数值模拟格式效应研究

针对跨声速叶轮机械复杂内部粘性流动,开发了全三维高精度RANS数值求解器。分别采用FVS格式与AUSM+格式结合多种限制器对Rotor 67跨声速压气机转子进行数值模拟,分析和比较了各格式的计算效果并与试验结果对比。结果表明:与FVS格式相比,AUSM+格式的数值粘性更小,边界层的模拟精度更高;Hemker限制器的综合表现最优,Van Albada限制器的粘性分辨率略低于Hemker限制器,Minmod限制器对流动分离现象的捕捉能力较差,Van Leer限制器容易引入色散误差。     

跨声速风扇叶片变形对气动性能的影响

为了研究跨声速风扇叶片变形对气动性能的影响,发展了流固耦合方法计算风扇叶片在气动力作用下变形的数值程序,耦合求解了三维可压缩雷诺平均的Navier-Stokes (N-S)方程和描述叶片振动的结构动力学方程,并模拟了跨声速风扇Rotor 67转子叶片在气动载荷下的变形,重点分析了由于叶片变形带来的气动性能的变化.      

低雷诺数下跨声速转子流动失稳及周向槽处理机匣扩稳

采用数值方法模拟了低雷诺数条件下NASA Rotor 37跨声速压气机转子内部流场。结果表明,附面层径向涡是该压气机转子流动失稳的一个很重要的原因。根据该压气机转子低雷诺数条件下流动失稳的特点,研究了周向槽处理机匣结构对其性能的影响。结果表明,引入处理机匣后,附面层径向涡得到一定程度的抑制,由附面层径向涡所引发的叶顶阻塞区有所减小,提高了压气机转子的失速裕度。     

Multivariate interpolation for fluid-structure-interaction problems using radial basis functions

A multivariate interpolation scheme for coupling fluid (CFD) and structural models (FE) in three-dimensional space is presented using radial basis functions. For the purpose of numerical aeroelastic computations, a selection of applicable functions is chosen: a classical without compact support, and some recently presented smooth compactly supported radial basis functions. The scheme is applied to a typical static aeroelastic problem, the prediction of the equilibrium of an elastic wing model in transonic fluid flow. The resulting coupled field problem containing the fluid and the structural state equations is solved by applying a partitioned solution procedure. The structure is represented by finite elements and the related equations are solved using a commercial FE analysis code. The transonic fluid flow is described by the three-dimensional Euler equations, solved by an upwind scheme procedure.     

附面层抽吸对低雷诺数下压气机稳定性的影响

数值模拟了低雷诺数下附面层抽吸流动控制对跨声速压气机稳定性的影响。低雷诺数下附面层径向涡及其诱发的叶顶分离阻塞触发压气机流动失稳,通过在NASA Rotor 37跨声速轴流压气机转子叶片吸力面上设计抽吸槽,探讨了抽吸流动控制对低雷诺数下压气机性能和稳定性的影响．并对比分析了抽吸前后压气机流场特性的变化和抽吸流动控制提高低雷诺数下跨声速轴流压气机稳定性的作用机理。     

跨音速压气机附面层组合抽吸方案对比分析

以NASA Rotor35为研究对象,应用数值模拟方法深入研究了叶表抽吸和端壁抽吸的不同抽吸组合布局对该跨音速压气机流场结构以及其性能影响规律.研究结果表明:四种抽吸组合在不同的抽吸量下均使压气机的压比和效率有所提高,相对叶表抽吸或端壁抽吸,其在较大的抽吸量时对压气机性能的提升较大,在抽吸量Sc=2.0%时最佳抽吸组合...     

跨音速透平叶栅多目标优化设计

本文在应用二维Euler方程及边界层方程相结合的跨音速粘流的计算方法基础上,以叶栅损失和做功能力为目标函数,采用无量纲化的多目标最小偏差法构造统一函数,然后采用可变容差法进行优化求得较为满意的解,从而形成了一种带有多混合变量、多约束以及多目标的跨音速叶栅优化设计方法。