单列叶栅紊流流场有限分析数值研究

本文将有限分析方法用于曲线座标系上紊流 N-S方程数值计算,研究了高雷诺数时叶栅粘性紊流流场。有限分析方法在网格单元上对非线性偏微分程进行线化处理,在解析边界条件下求出线化偏微分方程的解析解,以此解构造离散方程。有限分析方法能够根据对流的方向和大小自动改变格式系数,并具有数值扩散小、精度高和稳定性好等优点。本文以 k-ε紊流模型模化紊流,以壁面函数方法处理近壁区流动参数。      

紊流射流掺混流场的有限差分计算

本文给出了二维紊流射流掺混流场的有限差分数值计算方法。用这种方法可以计算等截面和变截面二维(平面的和轴对称的)可压缩紊流射流掺混的流场分布。本方法已经用FORTRAN算法语言编制成计算程序,在FELIX C-256电子计算机上进行运算。通过算例和实验结果表明,本文提供的方法是可行的。     

一种用于分离流动的网格自适应算法

对于可压缩粘性流动,提出利用流场速度的紊乱度作为指示变量进行网格自适应.Jameson中心格式有限体积法、五步Runge-Kutta时间推进法/双时间推进法求解定常/非定常N-S方程.基于雷诺平均N-S方程模拟紊流,选用SA一方程模型.在数值求解二维静态失速和动态失速问题过程中,加入网格自适应算法,提高数值模拟对流动分离特性的捕捉和分辨能力.算例结果表明在流场发生失速后,运用本文的自适应算法能够在增加少量网格单元的情况下明显提高计算精度.     

N-S方程有限分析数值研究

 本文将有限分析方法用于曲线座标系上紊流N-S方程的数值计算。分别计算了单列和串列叶栅内部流场,计算中采用了k-ε紊流模型和壁面函数。计算结果与试验结果相比较,吻合程度令人满意。有限分析方法在网格单元上对N-S方程进行线化处理,以解析边界条件作为约束,得出解析解,在解析解基础上构造离散代数方程。有限分析方法的最大特点是可以适应对流速度大小和方向,自动调整格式系数,因而具有数值扩散小和稳定性高等优点。     

粘性超声速绕体流动的有限谱法计算

首先推导出粘性可压缩流体动力学基本方程组在广义坐标下矢通量分裂的一般形式，并将有限谱ENO格式用于求解超声速流动问题。通过典型算例，验证了这一格式的精度和稳定性。然后进行了圆柱和机翼超声速绕流的数值模拟，得到了与实验和其他计算相符的结果。这表明有限谱法不单具有精度高，应用灵活等多方面的优点，而且突破了传统谱方法全域性的局限，可以成功地应用到一般形状流场可压缩流动的数值计算中。     

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

对SIMPLE方法加以推广,使之适用于可压缩粘性流动的数值计算。采用非交错网格技术,推出三维任意曲线坐标系下可压缩形式SIMPLE方法的计算方程。对亚音速、跨音速及超音速等5个流场进行了数值计算并与有关文献及实验数据进行比较。      

粘性不可压缩流动三维复杂流场分块耦合求解

提出了一个粘性不可压缩流动三维复杂流场的分块耦合求解方法。可解决科学和工程计算中的各种稳态和瞬态粘性不可压缩流动计算问题。     

亚超音速全机纵向气动系数数值计算方法的改善与发展

本报告是在亚超音速位流解基础上改进发展的数值计算方法。用面分布有限基本解计算作用在物面上的压力,用数值积分方法计算作用在物体上的力与力矩。本文采用联合流场、有效剖面吸力概念以及顺流向片条作可压缩紊流附面层的积分解,来改进阻力特性的数值计算。 数值方法已编制成计算程序。 本文提供了用这个计算程序计算的若干算例,并与实验结果作了对比。比较表明,计算结果与实验值有良好的一致性。     

直升机旋翼悬停流场的粘性数值模拟

在旋转坐标系下求解可压缩NS方程,对Caradona-Tung旋翼的跨声速悬停流场进行了数值模拟,计算中采用嵌套网格方法,并且应用了多重网格技术加速流场收敛。计算结果表明,多重网格技术显著改善了流场的收敛效率,计算量仅为单重网格的四分之一。文中同时还对无粘和粘性流场进行了对比,NS方程考虑了物面的粘性作用,能正确模拟出激波和边界层的干扰现象,计算得到的桨叶表面压力分布与实验值更加吻合。     

带圆筒头部火焰筒两相燃烧流场及性能计算

在三维任意曲线坐标系下采用代数雷诺应力紊流模型模拟紊流粘性,Arrhen inus-EBU紊流模型模拟燃烧室内化学反应速率,采用随机轨道模型模拟气液两相之间的相互作用,数值分析不同进口气流参数对带双圆筒头部环形燃烧室两相紊流燃烧流场的影响,同时采用多维经验分析法预估燃烧室性能。通过两种工况计算表明:在最大工况下所得的燃烧流场和燃烧室性能均优于慢车工况。计算结果与试验数据符合较好,表明燃烧过程的数值模拟以及多维经验分析法可为燃烧室优化设计提供有用的设计依据,适合于工程应用。     

考虑叶片径向和垂直于壁面方向导热的涡轮叶片对流冷却模型研究

为了提高涡轮叶片对流冷却模型预测精度,提出了一种在叶片固壁内同时考虑叶片径向和垂直于壁面方向(法向)导热的二维对流冷却模型。该模型在弦长方向划分多个元素,忽略元素内弦长方向叶片温度变化,在元素内的径向和法向建立二维导热方程作为叶片固壁温度场的控制方程,其边界条件包括叶表燃气绝热温度、燃气侧对流换热系数和叶片叶根、叶顶热流密度等。给出了该模型二维导热方程和边界条件的差分求解方法。以E~3涡轮高压导叶为例,将模型与CFD计算的叶片外壁面温度分布进行了对比。结果表明,该模型在给定冷气量下预测的叶片温度分布变化趋势与CFD相近,最大温度误差不超过6.5%,计算时间与CFD相比缩短了95%,能够快速、准确预测涡轮对流冷却叶片的冷气需求量。     

基于优化技术的叶型反方法改进设计

在势流函数反方法设计叶型的过程中,通过分析叶片表面速度分布对叶型形状的影响规律,发现采用局部区域速度分布调整,可有效的克服以往反方法设计叶型中易出现的非物理解现象,并针对叶型前缘和尾缘的几何封闭问题分别建立了优化目标函数,借助流场求解中的优化迭代技术获得了具有良好气动性能并满足几何约束的叶型设计结果。对典型叶型的三个设计算例表明:所研究的叶型反问题优化设计方法,具有设计精度较高,适用范围较宽,工程应用性较好等优点。      

一种叶轮机三维叶型黏性反问题方法

提出了一种应用于叶轮机三维叶片气动设计的黏性反问题方法,该方法将叶片表面边界视为是运动的壁面,通过建立并求解壁面法向动量方程,从而可根据给定的压力分布或压力面与吸力面之间的压力差分布直接求解叶型几何.结果表明:该方法可与定常雷诺平均Navier-Stokes(RANS)模拟配合进行求解,鲁棒性强、计算耗时短,且符合真实黏性流动环境,无需对叶型进行黏性修正.算例验证结果显示,该方法在单排和多排叶片环境下都能获得反问题解,计算稳定性和收敛性良好.      

透平机械叶片的遗传优化设计

提出一个基于遗传优化理论的透平机械叶片设计方法。该方法以叶片形状作为优化对象,利用遗传算法通过使其表面边界层中流动损失的极小化来搜索最佳的叶片形状。叶片形状被参数化表示。已知叶片形状的流场分析由一个叶栅正命题CFD程序完成。方法应用于一个离心压缩机扩压器叶片设计。数值计算显示,该方法可成功地求得具有最小流动损失的扩压器叶片形状。     

某小型离心压气机气动设计

针对无人机动力装置的需求,为某小推力的涡轮喷气发动机进行了单级离心压气机的气动设计.所用的压气机设计方法为:从一维的热力计算出发来进行流道结构设计,以二维通流计算为基础,通过调整各排叶片的环量分布,来获得所要求的气动性能,并通过叶片任意造型最终生成与通流计算结果较为一致的三维叶型,最后通过三维计算流体动力学(CFD)计算来分析气动设计的性能结果.最终设计出的离心压气机在设计点的各项性能参数均达到了设计要求,并且在各个转速下压气机的特性线比较平缓,压气机有较为宽广的稳定工作范围.      

多级轴流压气机端壁边界层中叶片力亏损的研究

本文基于对有经向间隙和无任向间隙两种情况下叶栅二次流造成边界层内叶片力的偏转的分析建立了计算叶片力亏损的模型。文中给出了由瑞壁边界层叶片力亏损计算压气机效率的方法。用上述方法对两台轮流压气机和一个进口导向叶片排计算了端壁边界层厚度,得到了环面堵塞系数,计算了一台压气机的效率。对计算结果的分析表明,用这种方法可见得到较合理的计算结果。      

求解跨声速压气机叶栅粘性流动反问题的数值解

阐述了一种以有限体积的时间推进方法为基础求解跨声速压气机叶栅设计的反问题方法。应用这一技术设计叶型,规定叶片表面无量纲目标速度分布,通过比较目标速度分布和计算获得的速度分布来修改叶片压力面和吸力面的坐标,最终获得要求的叶片形状。计算程序使用分布体力方法模拟粘性,用局部时间步长和多重网格方法加速收敛。还介绍了这种方法求解的基本方程系统,并给出了采用这种方法设计的无激波超临界跨声速压气机叶栅的设计结果。     

涡轮冷却叶片参数化建模及多学科设计优化

建立了一个涉及结构、气动、传热、振动、强度和寿命等学科的涡轮冷却叶片多学科设计优化系统, 进行了单孔薄壁冷却叶片的多学科设计优化.提出了单孔薄壁冷却叶片的参数化造型方法, 叶片叶型采用5次多项式构造, 气动与传热为三维耦合分析;叶片体积平均温度与最高温度为优化目标, 强度、振动和寿命等学科相关参数为约束, 模拟退火与序列二次规划组合算法进行叶片参数空间寻优, 在保持冷却气体流量不变的条件下, 优化提高了冷却效果, 降低了叶片材料的性能要求.      

高亚声速压气机叶片优化设计

为实现压气机叶片的优化设计,采用Hicks Henne函数进行叶型参数化,N S方程流场计算与混合遗传算法结合构成设计软件。以给定叶片表面压力分布为目标,以损失小而扩压度大和给定压比损失最小作为目标,所得优化叶片吸力面等熵马赫数分布合理、符合控制扩散规律,具有较好的压比和损失指标。采用几何方法与椭圆型方程方法结合生成壁面正交H型网格,可提高计算精度和便于采用代数紊流模型的流场计算。     

透平叶顶泄漏能量损失的数值计算

商用N S求解器已经被用于透平动叶顶部空隙的非定常流场的模拟和研究。本文介绍了一种新的用于透平动叶片顶部空隙区域流场计算的能量损失方法 ,可以方便地用于各种流场。发现高黏度效应区域位于上端壁附近而不是在叶片顶部。研究表明对于不同的动叶片其顶部泄流的能量时均损失是不一样的。该结果提供了一个有用的线索 ,可以通过叶片设计者设计出不同几何和空气动力学参数的非统一的动叶片 ,从而减小能量损失