单级跨音速压气机三维非定常粘性流动计算

本文以数值求解雷诺平均非定常N- S方程方法完成了单级跨音速压气机动静叶相干三维非定常粘性流动计算,给出了高分辨率NND格式的时间推进有限差分解。对某跨音速单级压气机进行了定常和非定常流的详细数值计算,设计转速下定常粘性流场计算所得到的性能参数与实验结果吻合较好。在此基础上,对最高效率状态进行了动/静叶相互干扰的非定常流动计算,并对结果进行了分析。非定常流计算模式更接近于真实流动,在高性能压气机的设计与分析中具有广阔的应用前景     

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

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

组合抽吸对跨音速压气机稳定性影响分析

为研究组合抽吸对跨音速压气机稳定性的作用，以NAsARotor35为研究对象，应用数值模拟手段深入分析了不同抽吸方案对该压气机稳定性以及流场结构的影响效应，数值计算采用商业软件NUMECA结合S-A模型，求解定常NS方程组。研究结果表明：叶表和端壁组合抽吸不仅对压气机的压比和效率性能有较大的提升，同时对压气机内部的流动失稳机制有抑制作用，是一种可行的流动失稳控制策略；叶表抽吸缝主要是吸除叶片吸力面附面层内的低能高熵流体提高压气机的效率，端壁抽吸主要是吸除叶顶内低能高熵流体，改善叶顶复杂流动结构，进而扩大压气机稳定工作范围。     

压气机叶片排内三维紊流流动的隐式高分辨率数值分析

对跨音速压气机叶片排内具有复杂三维激波的紊流流场所进行的数值分析，采用了空间三阶精度无波动格式以保证对紊流特征的正确模拟，利用ＬＵ－ＳＧＳ算法避免了一般隐式算法的繁琐运算，使之与显式方法推进一步的运算工作量基本相当，从而构成了一种既准确又高效的跨音速压气机紊流流动数值求解系统。计算与实验进行了对比，结果是令人满意的。     

跨音速压气机回转叶片排中三维流场的数值分析

本采用时间分裂显式格式求解Euler方程，获得了跨音速压回转叶片排中的三维流场解。研究了三维扭曲网格的生成，边界条件的处理以及保证稳定和提高流场分辨率的措施等问题。计算结果与实测三维流场的对比表明了本方法的可靠性。对某型高负荷跨音速压气机第一级转子叶片排所作的流场分析计算，为判断该叶片排流场品质，分析叶片设计的合理性提供了依据，证明了本方法的工程实用性。     

用于粘性可压缩流动数值计算的SIMPLE方法

对ＳＩＭＰＬＥ方法加以推广，使之适用于可压缩粘性流动的数值计算。采用非交错网格技术，推出三维任意曲线坐标系下可压缩形式ＳＩＭＰＬＥ方法的计算方程。对亚音速、跨音速及超音速等５个流场进行了数值计算并与有关文献及实验数据进行比较，结果相符很好，且有较快的收敛速度、有较大的应用前景     

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

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

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

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

多级风扇三维粘性流场的数值模拟

 进行了多级跨音速风扇粘性流动的数值模拟研究。求解全三维N-S方程,通过高分辨率格式提高跨音速流场计算的准确性,利用LU-SGS隐式解法加快计算收敛速度。给出了一双级高负荷跨音速风扇设计和非设计工况下的详细计算结果,各工况下的计算与实验数据均吻合较好。     

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

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

跨音速平面叶栅全位势无粘-粘性迭代计算

本文发展了一种全位势与附面层、尾迹流正逆模式积分法迭代计算程序,给出了计算结果与国内外四个不同的叶栅风洞试验结果的对比分析。结果表明,对比是较为满意的,本文在网格生成,无粘跨音流求解和无粘流与附面层、尾迹流迭代等方面所作的发展和改进适用,使整个程序系统具有工程实用价值。     

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

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

跨声速操纵面嗡鸣的数值研究

我国的一些高速歼击机，在飞行中曾多次出现操纵面嗡鸣问题，它严重影响了飞机的性能和飞行安全，本文采用数值计算方法，对此问题做较深入的分析和研究。气动力计算采用Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，并采用ＬＵ－ＴＶＤ混合差分格式求解由此气动力解操纵面的弹性振动方程，求出现嗡鸣的Ｍ数及振动的频率、振幅等参数。计算中采用Ｂ－Ｌ代数湍流模型，网格生成采用代数方法及保角变换方法，文中讨论了隐式算法在非定常气动力计算     

多级轴流压气机非设计性能的数值预估

采用LU-SGS隐式解法以及叶排间掺混面模型数值模拟压气机三维粘性流动,研究了多级轴流压气机的设计与非设计状态性能的预估问题。对一台五级高压压气机进行了详细的流场数值计算,包括了设计转速的100%,90%及80%等多个转速及各种工作状态。结果表明,本文所建立的数值方法与计算机程序对多级轴流压气机各种非设计状态的计算都具有收敛容易和快速的优点,计算得到的压气机总性能和各分级性能与设计参数进行了对比,证明该方法和程序的预估结果是比较准确的。根据设计转速下的计算结果,预估了压气机设计点的喘振裕度。此外,还对不同状态下叶片进口气流角的变化进行了分析。     

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

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

固体火箭喷管两相粘性跨音速流场计算

本文进行了固体火箭喷管两相粘性湍流跨音速流场计算.粘性的气相控制方程用隐式近似因子分解法求解,粒子方程采用跟踪粒子轨迹的特征线法求解,粘性湍流选用代数模型,气相和凝相充分地偶合.CFL数可以取500左右,收敛速度快,使粘性两相跨音速喷管流场计算耗费的机时达到工程计算可以接受的程度.通过计算,获得了在粘性和粒子同时作用下的流场参数分布.这对固体火箭发动机流场需要同时考虑粒子和粘性作用的专题研究大有裨益.     

跨音速压气机动静叶排相干的三维非定常流动数值分析

采用时空二阶精度ＮＮＤ显式差分格式求解三维Ｅｕｌｅｒ方程，对跨音速压气机级中动静叶排相干形成的非定常流动进行了数值分析。除动静叶排流动的控制方程体系、数值方法外，着重讨论了动静叶排计算域、网格以及边界条件的处理。对某压气机第１级动静叶排相干的非定常流场的计算，给出了叶片表面的马赫数分布、压力波动幅值，以及叶片槽道中瞬态密度等值线等，并分析了结果的合理性。     

轴流压气机孤立转子内流场数值模拟

针对跨声速轴流压气机转子叶排中的流动，利用高分辨率的三阶NND格式和LU-SGS隐式推进迭代法，建立了孤立转子内三维N-S方程的高分辨率和高效率的数值分析程序。用该程序分别计算了NASA-37跨声速转子，计算结果表明：所建立的孤立转子程序能够较好的模拟轴流压气机叶排内流动情况。数值分析结果与实验吻合较好，从而证明了该程序的有效性。     

跨声速轴流压气机级全工况三维流场数值分析

以LV－SGS－GE隐式格式和MUSCL TVD迎风格式为基础，结合壁面函数方法和简单的混合长度湍流模型，对三维可压缩雷诺平均N－S方程进行求解。叶列间参数的传递采用混合平面方法并应用了微机网络并行计算技术。计算得到了NASA 35号低展弦比、跨声速轴流压气机组70％和90％设计转速下的全工况性能曲线，并重点分析了其中一些典型工况下的内部流场。计算与实验结果的对比表明，此方法能快速得到三维粘性流场的流动特性且计算精度较高，可用来模拟跨声速轴流压气机组内的全工况三维粘性流动。