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

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

二元矢量喷管流场的N－S方程计算

对二元矢量喷管在不同工况下的流场进行了数值模拟。求解的控制方程为三维非定常Ｎ－Ｓ方程。空间方向用有限体积离散，时间方向用Ｒｕｎｇｅ－Ｋｕｔｔａ法推进。粘性项的计算上采用了一种有限体积框架下新的离散化方法。湍流模型采用了广泛使用的Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数模型。计算表明，本文方法能准确地计算二元矢量喷管、非矢量喷管在设计与大设计落压比下的流场。     

火箭冲压发动机反应流场数值模拟

用偶合算法求解三维定常椭圆型Ｎ－Ｓ方程组,计算了火箭冲压发动机内部的湍流反应流场。用有限体积法离散控制方程,用块隐式法求解离散后的代数方程组。湍流流动用ｋ－ε模型描述,湍流燃烧采用掺混控制的局部瞬时不混合模型。实践表明,尽管密度场变化剧烈,仍能快速收敛。算例中网格点数为３９×１３×１８。计算出的燃烧效率符合实验结果。     

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

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

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

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

姿控发动机喷管内流场及推力计算研究

运用时间相关法数值求解定常Ｎ－Ｓ方程和采用有限体积方法离散控制方程,对姿控发动机喷管内流场及其产生的推力进行了计算,得到了较为精细的流场图象。计算结果表明,不同壁面温度条件可以显著改变喷管壁面附近的流场结构及其产生的推力。     

双时间法求解大迎角翼型绕流非定常N-S方程

本文采用双时间推进法、顶点格式有限体积法求解非定常ＮＳ方程、应用多层网格技术计算了ＮＡＣＡ００１２翼型大迎角振荡绕流,并对伪时间推进收敛精度,物理时间步长,多层网格迭代方式等对计算结果与计算ＣＰＵ时间的影响进行了讨论,另外还讨论了采用Ｂａｌｄｗｉｎ－Ｌｏｍａｘ湍流模型时用固定转捩或全湍流模式对计算结果的影响,并与实验结果进行了比较。     

具有复杂边界的火箭冲压发动机流场数值计算

针对实际复杂的补燃结构和计算边界,以块隐式法计算流场,用代数法生成计算网格,通过改变离散方程来处理三维定常Ｎ－Ｓ方程边界条件。算例表明,该方法正确可靠,迭代过程能够快速收敛。     

超燃冲压发动机燃烧室三维化学反应流场的数值模拟

采用ＬＵ隐式方法强耦合求解有限体积法离散的完全Ｎ－Ｓ方程及组分方程，数值模拟超燃冲压发动机机燃烧室机的化学反应流流场。计算中和了ＡＵＳＭ通量分裂格式，两方程的化学反应模型和Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代化学反应流场。     

用双时间法数值模拟悬停旋翼流场(英文)

本文用改进的双时间法和重叠网格通过求解可压缩的欧拉和Ｎ－Ｓ方程计算了悬停旋翼流场。在双时间法中,我们采用牛顿松弛以加速每个真实时间步中的子迭代。高效的双时间法加上重叠网格使在微机上计算复杂的悬停旋翼流场成为可能,计算结果和实验符合得很好,计算得到的下洗流和载荷证明了我们的方法是可行的。     

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

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

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

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

两类对转风扇的设计与气动特征数值研究

采用一维设计程序分析了前后转子设计转速比的影响,研究了平均半径处的增压比、绝热效率、扩散损失和激波损失随转速比的变化规律.用计算流体力学分析了设计点与非设计状态的两个对转级流场,研究了其详细物理现象.结果表明两个对转级的设计与非设计性能均良好.发现低速风扇的两个转子均为常规跨声速转子,而高速风扇的前转子常规,后转子则为前缘激波和通道激波均贯穿全叶展的全超声速转子.同时发现,均带有与常规风扇级相当的失速裕度,低速对转级是两个转子同时达到失速点并且激波被推出叶栅,而高速对转级则是后转子先达到失速点并激波推出,从而后转子决定着级失速裕度.      

一类轴流叶轮机转/静干扰问题的 3D Euler 数值研究方法

本文将作者发展的定常N-S解程序,改编成时间精确Euler解程序,用于研究转/静干扰问题。对一台高负荷风扇级无粘流场进行时间精确模拟的结果表明,该方法可被应用于轴流叶轮机内非定常流动机理研究。      

掠弯叶片前缘曲线同流场结构的关联

在全三维反问题还不成熟的条件下,由二维反问题和三维正问题组成的准三维设计体系中,叶片的前缘是由给定的积迭轴和基元叶型特征决定,具有任意性,难于对其加以控制。因此本文试图把掠弯叶片的前缘空间曲线作为设计自由度,增加到设计体系中,实现设计中前缘空间曲线同流场结构的关联,利用前缘形状控制流场中的激波结构和二次流的迁移,改善流场结构。并以一个高负荷跨音单级风扇的气动设计为例,对这一改进进行了验证。      

风扇转子叶片防颤振设计技术

以某两级风扇为对象开展防颤振设计技术研究。风扇转子叶片设计中，选择了合适的展弦比，既考虑到气动性能水平，还考虑到结构质量和颤振稳定性。采用流固耦合能量法评估所设计风扇的各排转子叶片的颤振稳定性问题，并通过叶片厚度、三维造型、根尖弦长比等设计参数的调整消除颤振风险。研究结果表明:展弦比不应是方案设计阶段防颤设计唯一的关注参数；在叶尖跨声速的转速更容易发生颤振现象；较强的叶尖前缘激波会造成较强的流固耦合作用，形成复杂的气动功分布结构；叶片厚度和尖根弦长比等参数是改善叶片颤振风险的有效参数。      

跨声速工况下流体诱发叶片振动研究

对同一跨声速风扇的非定常流场进行了跨声与亚声两种工况下的数值模拟.对两个工况下流场计算得到的转子叶片表面非定常气动激振力进行频谱分析,并对其进行对比分析,从而得到跨声工况下气动激振力特性.应用ANSYS软件的瞬态分析方法进行了叶片的强迫响应分析,对跨声和亚声两种工况下的叶片强迫响应进行对比分析,从而得到跨声工况下叶片强迫响应特性.      

喷管跨声速流场计算

本文用时间相关法求解定常跨声速喷管流场。对守恒型和非守恒型方程采用全场统一时间步长、当地时间步长和修正超阻尼法进行了计算。计算表明:采用当地时间步长可大大提高收敛速度,采用全场统一时间步长时,数值解非守恒型方程收敛较快;采用修正超阻尼法也能提高收敛速度,但与阻尼系数和时间间隔有关。     

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

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

提高跨声速风扇转子气动负荷系数途径

回顾了跨声速风扇/压气机转子气动负荷的研究历程,并对基准风扇转子和两个方案转子进行了三维流场分析。通过对195个研究方案的对比研究,目前阶段数值模拟的结果中风扇转子气动负荷系数达到0.65。在来流条件相同的情况下:合理组织激波结构可以大幅提高叶尖气动负荷;叶根大转角设计是提高根区气动负荷的可行途径;静子跨声速来流是高气动负荷跨声速风扇/压气机的必然结果。