用于压气机流动计算的3种模型比较

为了更好地对压气机流动进行模拟,在课题组自行开发的结构化有限体积解算器上实现了用于压气机流场计算的混合平面法、谐波平衡法及相滞后法.以NASA Stage 35为例,对3种方法的计算结果进行了比较分析.结果表明:相滞后法的计算精度最高,混合平面法的计算精度最低;相滞后法与半环的双时间推进法结果相近,计算速度提高了20倍;相比混合平面法,谐波平衡法能准确地模拟动静叶间的非定常干涉及进出口参数变化;在谐波阶数达到5阶后,谐波平衡法计算得到的结果不随阶数变化,且与相滞后法的结果基本吻合;混合平面法的计算效率远高于另外两种方法,相滞后法与谐波平衡法在谐波阶数为5阶时的计算效率相当.      

不同雷诺数下翼型气动特性及层流分离现象演化

低雷诺数下空气黏性效应突出,翼型表面普遍存在层流分离现象,相比常规雷诺数情况气动特性显著恶化。采用带预处理的Roe方法求解非定常可压缩Navier-Stokes方程的数值模拟技术和低雷诺数低湍流度风洞油流显示试验技术,对FX63-137翼型不同雷诺数下气动特性和流动结构展开深入研究。通过风洞油流显示试验可以清晰获得低雷诺数层流分离流动的两道油流汇集线。数值模拟结果表明其分别为时均化主分离线和二次分离线,两种结果定性定量均吻合较好,证明了本文的研究方法有效可靠;雷诺数从500 000降至20 000,翼型气动特性和层流分离流动结构均发生显著的变化,伴随阻力系数剧增和升力系数剧降,时均化流动结构从附体至出现经典的长层流分离泡,并最终演化为后缘层流分离泡,相应的两种分离泡的非定常流动结构也存在显著差异;对于阻力系数和升力系数而言,存在不同的临界雷诺数,因为导致阻力系数剧增的机理在于经典长层流分离泡的产生使翼型压差阻力大增,而造成升力系数剧降的主要原因在于后缘层流分离泡使得等效翼型后部弯度减小;非定常结果显示正是由于翼型表面漩涡周期性的生成与脱落,才造成了低雷诺数下升力系数的周期性波动。翼型上表面主分离涡即将脱落时,流线在后缘附近再附,升力系数达到峰值;而当流体从下表面向上卷起二次分离涡时,尾部流线大尺度分离,升力系数降至谷值。     

螺旋桨飞机滑流非定常数值模拟研究

采用非定常数值模拟方法,基于滑移网格技术,分析研究螺旋桨滑流对全机气动特性的干扰影响。用多块结构网格数值模拟软件,从搭接边界建立、搭接网格生成、流场插值传递3个方面对螺旋桨滑移搭接网格的数值模拟方法进行全面详细的介绍,并且阐述基于该技术的非定常数值模拟方法。在此基础上,开展螺旋桨飞机模型验证计算,分析研究螺旋桨滑流对全机流场和气动性能的影响,通过与风洞实验数据的对比分析,验证滑移搭接网格技术的可行性与准确性。研究结果表明,基于滑移搭接网格技术的螺旋桨滑流数值模拟,可以真实反映螺旋桨滑流流场特性,为螺旋桨类飞机的设计评估提供一种有效的数值模拟手段,在工程方面具有重要的实用价值。     

粒子图象电影测速技术

本文介绍一种粒子图象电影测速技术,该技术采用机械调制装置及光学透镜,将连续输出的氩离子激光调制成周期性脉冲片光源序列,并与快速照相枪匹配,可以连续拍摄流场一个截面上随时间发展过程中不同瞬态的一系列两次或多次曝光的粒子图象底片,因而不仅可以取得流动一个截面上的瞬时二维速度场,并且可以获得该流动截面速度场发展的时间历程,为定量研究非定常流动随时间的发展特性提供了一种有效途径。     

亚格子模型对泵喷推进器宽带非定常力预测的影响研究

为了研究亚格子模型对泵喷推进器非定常流动与宽带非定常力预报结果的影响,采用分块结构化网格建立了模型尺度下艇后泵喷推进器的计算模型,并进行了大涡模拟数值计算。从艇尾非定常流场特征量和推进器转子脉动载荷两个方面对比了三种不同亚格子模型计算结果的差异,并分析了泵喷内部流动与转子非定常力间的内在联系。研究结果表明：三种亚格子模型均能得到含有叶频宽带峰的轴向推力谱,且整体趋势相近。泵喷转子上游的湍流强度和尺度的分布对亚格子模型较敏感,其中Smargrinsky-Lilly模型得到的湍流强度较强,尺度较大,该模型下的湍流谱在非平衡区量级较大,但由低频向高频的衰减较快,并且预测到的分离区范围大,导叶尾缘脱落涡分散,导致转子上游来流空间分布不均程度较强。对于转子叶片上的载荷脉动,Smargrinsky-Lilly模型预测的推力谱中线谱成分明显,并且叶频处宽带谱峰“陡峭”,而WALE模型和KET模型的结果宽带特性较强,对比非定常推力测试结果可知,WALE模型和KET模型更适于该问题宽带非定常力预测。     

绕三角翼流动中的非定常现象研究

利用流动显示技术、片光技术及PIV技术,对三角翼面上的各种非定常特性做了研究.分别对涡层中的小涡结构,二次涡同主涡相互作用所引起的二次涡的喷射效应,螺旋破裂扰动的非定常特性,涡破裂点沿轴方向振动的非定常特性,以及完全分离流中的非定常特征进行了较为细致的研究,给出了不同扰动的频率特性.      

导弹滚转通道气动伺服弹性稳定性分析

弹性飞行器可能由于结构、气动及控制的不利耦合引发气动伺服弹性失稳,为此分析了带有差动舵面及控制系统的细长体导弹滚转通道气动伺服弹性稳定性.导弹的广义非定常气动力采用气动导数方法计算,动力学方程考虑了弹体的刚体滚转运动和一阶弹性扭转振动.由导弹动力学方程求解得到滚转通道的开环传递函数,并应用Nyquist方法进行闭环系统稳定性分析.算例表明,若控制系统设计不当,系统有可能在弹性振动频率处发生气动伺服弹性失稳.      

圆柱绕流尾迹中涡的特征分析

在低雷诺数圆柱绕流中,卡门涡街中涡的形状随着其向下游的发展是不断变化的,而且它也不是理想轴对称的.针对该问题,利用有限元法进行数值模拟,研究了由涡量等值线表示的涡的形状特征.由流场中几种等值线的对比可以得到,主要是涡量等值线与压力等值线之间的差异造成了涡的形状变化.利用广义的茹柯夫斯基变换描述涡相对于轴对称的变形状态.结果表明,描述涡的主要3个参数随半径的变化规律不同,即椭圆率的抛物形分布,偏心率的线性分布和弯度的近似常量的分布.在此基础上给出了3个参数各自对应的运动学量表达式,两者对比的结果还是比较相符的.     

对转螺旋桨流场气动干扰数值模拟

基于结构网格动态面搭接技术,通过求解三维非定常雷诺平均纳维-斯托克斯(RANS)方程,对单独螺旋桨(SRP)滑流流场进行了数值模拟,结果与风洞试验数据吻合良好,在此基础上深入研究了对转螺旋桨（CRP）的非定常滑流效应及前后桨之间的气动干扰现象。计算结果表明:所采用的非定常RANS方法能够良好刻画对转螺旋桨滑流流场的发展变化特征,适用于分析前后桨气动干扰特性。前后桨桨尖涡的复杂相互作用是导致对转螺旋桨气动干扰的直接原因,后桨桨尖涡对于前桨桨尖涡具有一定耗散作用。在气动干扰影响下,对转螺旋桨气动力会产生周期性波动,波动周期与螺旋桨桨叶数相关, 6×6对转螺旋桨在1个旋转周期内产生了12次波动。气流穿过对转螺旋桨会发生两次加速,因此在相同工况下,对转螺旋桨的拉力系数和功率系数相比于单独螺旋桨分别增大约1倍,效率提升约1.5%。      

摆线桨非定常动态失速的三维效应

对摆线桨三维全尺寸模型进行了非定常数值模拟。验证了两种嵌套网格方法的可靠性与适用性,针对二维/三维情况下摆线桨非定常涡流动特性及诱导分离特性开展了数值分析,重点研究了展向气动力随方位角的分布与全局流场的诱导速度分布,并分析了尾迹捕捉精度与涡量耗散特征。结果表明:受摆线桨切向速度、展向诱导速度等三维效应影响,两种计算结果的瞬态气动力极值差值达到44%。三维动态失速涡的产生、脱落与再附明显弱于二维翼型情况,这对于摆线桨非定常气动特性有较大影响。