跨音速压气机转子中三维湍流流场计算及涡系分析

从叶轮机械转动坐标系中的三维、可压、湍流、Ｒｅｙｎｏｌｄｓ平均Ｎ－Ｓ方程组出发,用Ｂａｌｄ－ｗｉｎ－Ｌｏｍａｘ湍流模式使方程组闭合;结合有限体积离散并采用隐式矢通量分裂技术;求解了西德宇航院（ＤＦＶＬＲ）的单级跨音速压气机转子内的三维湍流流场。计算结果可分辨出该转子内三维流场的细微结构、激波的空间曲面;得到马蹄涡、通道涡、尾迹涡和角涡的形态及其发展,流场的三维性非常明显。计算得到的壁面极限流线图可反映出压气机中分离流动的拓扑结构     

跨音速压气机转子内三维流场数值模拟

采用多步Runge-Kutta与有限差分相结合的格式,离散叶轮机械内三维欧拉方程,应用代数法生成了压气机转子三维网格。用此方法模拟了NASA Rotor37跨音速压气机转子流场,计算得到的流场与试验结果进行了对比。结果表明,该方法收敛速度快,实用可行,计算简单省时。并且,该方法为研究叶轮机械内部流场计算提供了很好的基础。在此基础上进行更深层次的研究会变得很容易。     

高速微型离心压气机的改型设计

近年来微型涡喷发动机的研究受到国内外的高度重视。微型叶轮机作为微型涡喷发动机的核心,其内部流动状况及特性是研究的重点。本文首先采用S-A湍流模型,对某微型涡喷发动机的离心压气机(转速118 000 r/min)的地面状态和设计转速下的粘性流场进行了三维数值模拟计算,得到了压气机特性曲线。接着对该离心压气机进行了改型。通过与原型三维流场的对比分析,实施了改型流道设计,重点在于调整转子流道的形状以及转子出口与径向扩压器之间的流道匹配。改型离心压气机的三维流场计算结果表明,基本达到改型设计要求。     

跨声速轴流风扇转子流场三维粘性数值分析

采用一种快速求解三维粘性流场的计算方法求解跨声速压气机风扇转子内部流场。该方法以ＬＵ型隐式格式和ＭＵＳＣＬ ＴＶＤ迎风格式为基础,结合壁面函数方法和简单的混合长度湍流模型,使用多重网格迭代加速收敛技术,对三维可压缩雷诺平均Ｎ－Ｓ方程进行求解。计算得到了ＮＡＳＡ Ｌｅｗｉｓ ６７跨声速、低展弦比轴流风扇转子的性能曲线,并重点分析了最高效率点工况和近失速工况下的内部流场。计算与实验结果的对比表明,此方     

压气机转子三维流场计算中叶尖间隙网格的影响

对一带叶尖间隙的低速轴流压气机转子三维流场进行了数值模拟,通过与实验结果对比,分析了叶尖网格匹配、网格质量对压气机转子三维流场计算的影响,发现网格正交性和网格疏密过渡会影响二阶计算精度的收敛,研究还显示压升特性线不能准确的反映网格数目对计算结果的影响。      

吸附式压气机转子叶片气动优化设计

进行吸附式压气机转子叶片气动设计方法研究.提出通过迭代设计,应用三维流场计算结果提供S₂流面通流计算损失模型,提高了S₂流面流场计算精度;将最优化方法与数值模拟技术相结合,建立吸气与型面耦合的回转面二维叶型优化设计和三维叶片优化设计软件;应用所构建的设计软件,进行吸附式压气机转子设计,数值模拟结果表明:该转子在0.86的叶尖载荷下,设计点总压比为1.631、等熵效率为0.965,实现了高气动性能.      

跨音速低展弦比轴流风扇转子内部流场的数值研究

采用三维雷诺平均Ｎ－Ｓ方程程序求解跨音速压气机风扇转子内部流场。这个程序采用有限体积显式时间推进方法求解控制方程,方程自变量选取在控制体顶点,采用简单的Ｈ型网格离散方程,离散后的代数方程用两步Ｒｕｎｇｅ－Ｋｕｔａ法进行时间推进求解,使用多重网格和当地时间步长方法加速计算收敛。程序中使用了Ｂａｌｄｗｉｎ－Ｌｏｍａｘ模型模拟湍流影响。用该方法计算了ＮＡＳＡ６７＃跨音速风扇转子的性能曲线,并重点计算最高效率点工况下内部流场。计算结果与实验结果的比较,证明该法可以用来模拟跨音速风扇转子内部流动。     

轴流压气机子午流面计算中三维剪切作用的湍流模型

本文发展了一个轴流压气机子午流面二维计算中模拟三维剪切结构对子午流场作用的湍流模型。试验与这一模型进行了初步的比较。模型表明,在多级轴流压气机子午流场计算中,必须计入实际流场三维剪切结构对子午流场的湍流扩散作用。     

三维C-H网格下跨音压气机转子Euler流场数值模拟

介绍了三维C-H网格的生成和全三维转子流场的数值解。以法国宇航院设计的一个高负荷跨音转子为算例, 并与他们的计算结果作了比较。应用这套数值模拟技术得以计算带进口预旋的某跨音压气机转子流场, 并对进口预旋对其下游的转子流场的影响作了数值模拟, 得出了若干有意义的结论。      

压气机转子出口流场的发展及三维紊流特性

用单斜热丝和高频压力探针仔细测量了单级压气机转子出口及下游的三维紊流流场, 揭示转子出口流场的变化, 分析不同流量状态下单级压气机转子出口的三维紊流特性。      

考虑顶隙的压气机单转子内三维跨音速粘性流场结构的数值模拟

借助于考虑叶片顶部间隙的压气机转子内跨音速湍流流场的三维数值解,得到了该转子流场中的复杂涡系与波系结构。给出了叶片表面及轮毂面上的壁面流型,分析并揭示了转子叶道内的马蹄涡、通道涡、角涡、尾涡等复杂涡系及其与激波的相互干涉;特别是,分析并揭示了转子叶道内顶隙射流的特点,得到了顶隙附近激波与泄漏流动的干涉形态及泄漏流动的发展,并着重给出了其与激波的干涉形态。     

单转子风扇的三维反问题气动设计

采用准三维流函数反问题设计方法与三维 N-S求解方法的相互迭代对单转子风扇进行气动设计。在叶片初步计算中得到叶栅进出口气流角沿径向分布 ,并将它作为本文的目标函数。采用准三维反问题求解方法 ,依次构造出各个 S1流面上的叶片几何形状和气流角分布。然后再采用 N-S方程的求解方法 ,对叶片进行全三维流场的数值计算。通过 N-S方程计算结果与目标函数的对比 ,重新修正叶片出口气流角分布 ,并作为下一次反问题设计的目标函数。经过反问题与 N-S方程求解的反复迭代 ,最终得到满足设计要求的叶型。     

CFD simulation of helicopter rotor flow based on unsteady actuator disk model

Actuator Disks (AD) can provide characterizations of rotor wakes while reducing computational expense associated with modeling the fully resolved blades. This work presents an unsteady actuator disk method based on surface circulation distribution combined with empirical data, blade element theory and rotor momentum theory. The nonuniform circulation distribution accounts for 3D blade load effects, and in particular, tip loses. Numerical simulations were conducted for the isolated pressure sensitive paint model rotor blade in hover and forward flight using the HMB3 CFD solver of Glasgow University. Validation of CFD results in comparison with published numerical data was performed in hover, for a range of blade pitch angles using fully turbulent flow and the k-ω SST model. In forward flight, the vortex structures predicted using the unsteady actuator disk model showed configurations similar to the ones obtained using fully resolved rotor blades. Despite the reduced grid cells number, the CFD results for AD models captured well the main vortical structures around the rotor disk in comparison to the fully resolved cases.     

Kirchhoff方法在旋翼前飞噪声预测中的应用研究

将CFD(计算流体力学)技术与噪声预测的Kirchhoff方法相结合,发展了前飞状态下直升机旋翼远场气动噪声的定量预测方法,得到了旋翼噪声的时域解。先用具有解析解的“点源”流动代替旋翼流动,通过比较计算值与理论值,验证了该预测方法的可行性和可靠性;尔后针对AH—1／OLS旋翼模型的两种前飞状态进行计算,通过比较计算值与实验值,进一步验证了本文方法及其所发展的程序的正确性。     

轴流压气机转静二维非定常薄层N-S数值模拟

采用非定常薄层Navier-Stokes紊流运动方程来描述一跨声轴流压气机中径处的二维转-静非定常流场，速度脉动的相关分析表明转子尾迹对静子的热流流场几乎没有什么影响，而粘性的影响取决于叶片排之间的轴向距离。     

叶轮机转子内流动的三维跨声速Navier-Stokes方程组数值计算及试验验证

应用有限体积法及Baldwin-Lomax代数湍流模型求解叶轮机转子内流动的三维可压缩非定常薄层Navier-Stokes方程组,计算一个绕X轴等转速旋转的转子内流场问题,用当地时间步长方法来加速收敛.采用隐式算法对对流通量进行求解,而耗散通量采用显式算法.采用高阶TVD下的Roe格式计算对流通量,采用中心差分格式计算耗散通量.关于守恒变量Q的非线性方程组,采用牛顿迭代方法求解,线化后的方程用对称的Gauss-Seidle方法求其渐近解.应用该算法模拟了转子Rotor67流场,与试验做了大量的对比研究,包括性能、等马赫线以及总压和总温等.在最高效率点附近计算与试验比较接近,但在靠近失速点时计算与试验有较大的差别.     

Application of immersed boundary method in turbomachines

Simulating unsteady turbulent flow in turbomachines is still challenging due to the complexity of blade geometry and relative motion between rotor and stator. This study presents an Immersed Boundary Method (IBM) for high-Reynolds turbomachinery internal flows, and shows the advantage of the automatic grid generation techniques and flexible moving boundary treatments. The wall functions are used in the present method to alleviate the wall resolution restriction of turbulence simulation. The Two-Dimensional (2-D) IBM solver, which was previously developed and tested for a low-speed compressor, is further validated for a well-documented Low-Pressure Turbine (LPT) cascade. Both the blade loading and the total pressure losses in the wake are well captured by the present 2-D solver. The complex Three-Dimensional (3-D) effects in turbomachines motivate the further development of an extended 3-D IBM solver by using a curvilinear-coordinate system that facilitates the hub and casing boundary treatment. The good performance of the 3-D solver is demonstrated through comparison with CFX solver solutions for the rotor configuration of Advanced Noise Control Fan (ANCF). Further effects of the grid resolution on capturing the blade wake are discussed. The results indicate that the present 3-D solver is capable of reproducing the evolution of the blade wake with suitable computational grid.     

旋翼非定常动态失速的三维效应分析

采用建立的高精度计算流体动力学(CFD)方法,针对旋翼非定常动态失速的三维（3D）效应特性进行研究。以Helishape 7AD旋翼为基准,开展三维效应对旋翼非定常气动特性的影响分析。研究了来流马赫数对旋翼翼型动态失速特性影响。在此基础上,针对三维情形旋翼动态失速非定常涡流动特性及诱导分离特征开展了数值分析,通过与二维情形对比表明:受旋翼旋转、轴向诱导速度等三维效应影响,旋翼桨叶剖面动态失速涡的产生、对流和脱落明显滞后于二维翼型情形,并且涡强度也更弱。越靠近桨叶内段,桨叶剖面非定常动态失速特性与二维旋翼翼型情形的差距越为明显。      

先进直升机旋翼悬停状态气动性能计算

 为更好地模拟先进直升机旋翼流场,同时准确计算其悬停状态的气动性能,建立了一套基于Navier-Stokes/Euler方程的混合CFD方法。该方法的求解域由两部分组成：一是围绕旋翼桨叶周围的黏性区域,采用可压缩Navier-Stokes方程来模拟旋翼附近的黏性流动和近场尾涡的捕捉;二是离桨叶较远、黏性可以忽略的远场区域,用Euler方程来描述其流动。在该方法中,将三阶逆风格式（MUSCL）与通量差分分裂方法相结合,无需添加人工黏性,因而可有效地减少旋翼尾迹数值耗散。?阌诹鞒》智蠼庖约爸芷谛员呓缣跫氖凳?采用了嵌套网格方法,并给出旋翼网格与背景网格交界面的信息传递方式。应用所建立的计算方法,首先对二维翼型、三维M6机翼的流场进行了数值模拟,以验证计算方法;然后,着重计算了有实验结果可供对比的具有先进气动外形的HELISHAPE 7A模型旋翼和UH-60A直升机旋翼,通过计算得到旋翼表面压力分布、桨叶展向拉力系数分布、桨叶表面细节流动以及气动性能等,进一步验证了该方法的有效性。     

旋翼近地飞行时尾迹及地面涡的N─S方程模拟

用有限体积法、三维Ｎ－Ｓ方程对旋翼近地飞行时的尾迹及地面涡进行了数值模拟。由分布在特定区域内的动量源项模拟旋翼对流场的作用,动量源项的强度由旋翼的几何尺寸、叶片剖面的空气动力特性及当地流场的性质迭代定出。进而,研究了旋翼近地飞行时,在不同前进比及不同离地高度情况下,尾迹与地面的干扰问题及地面涡的现象