端区流动对跨声速压气机失速影响的数值研究

为了揭示近失速工况下动叶间隙和静叶角区的复杂流动结构,探索诱发压气机失速的主要因素,对高负荷跨声速压气机开展数值研究,获得了整级压气机和单转子的特性曲线,进而对压气机的流场进行了详细分析。结果表明:对于研究对象,失速最先起始于静叶上角区。上角区完成由开式分离到闭式分离的转变过程,可以认为是静叶发生失速的标志。在近失速工况下,动叶叶尖泄漏涡发生"泡式破裂"堵塞流道,泄漏涡破裂引起的堵塞作用一直从叶顶延伸至70%叶高左右。     

冲击距与气膜孔方位角对旋流气膜冷却性能影响

为了获得冲击距Hi与气膜孔方位角α对旋流气膜冷却性能的影响规律,以六边形供气腔圆形气膜孔平板气膜冷却结构为研究对象,对五种冲击距参数(0.74D,1.14D,1.54D,1.94D,2.34D)(D为气膜孔直径)和五种气膜孔方位角参数(0°,10°,15°,20°,25°)进行了三维数值计算研究,得到了绝热壁面气膜冷却效率、展向平均气膜冷却效率、流场空间无量纲浓度分布等随冲击距与气膜孔方位角的变化规律,分析了肾形涡对旋流气膜冷却性能的影响机理。结果表明:冲击距对绝热壁面气膜冷却效率展向分布规律影响不大,而方位角增加能够明显提高绝热壁面气膜冷却效率及展向气膜覆盖面积,方位角0°模型展向气膜冷却效率最大值为0.42,方位角25°模型的最大值为0.48,相比前者增加14.3%;绝热壁面同一流向位置的展向平均气膜冷却效率随冲击距的增加而增大,随方位角的增加而增大,Hi=2.34D时的展向平均气膜冷却效率最佳,α=20°时的展向平均气膜冷却效率曲线最佳。方位角的增加能够明显破坏流场中存在的肾形涡结构。     

单转子轴流压气机涡动力学失稳机理研究

为了研究单转子轴流压气机的涡动力学失稳机理,采用基于Shear Stress Transport(SST)湍流模型的尺度自适应雷诺平均/大涡(RANS/LES)混合模拟的方法对低速单转子轴流压气机进行了非定常数值模拟。研究结果表明:在设计转速3kr/min条件下,叶顶泄漏涡、二次泄漏涡以及诱导涡破碎引起的叶顶区域的堵塞是触发单转子轴流压气机内部流动失稳的主要因素。压气机由近堵塞工况点向小流量工况点逼近的过程中,叶顶泄漏涡轨迹与轴向的夹角由70°增加到76°,二次泄漏涡起始点位置前移加速叶顶泄漏涡向转子前缘移动。近失速工况点叶顶泄漏涡的轴向动量与主流的轴向动量之间存在一种平衡,叶顶泄漏涡稳定在转子前缘。压气机进一步节流主流的轴向动量减小,对叶顶泄漏涡轴向动量的抑制能力减弱,叶顶泄漏涡的位置不再稳定,诱发尖脉冲型失速先兆。     

基于 LES 方法的增升装置气动噪声特性分析

在气动噪声数值计算中,流场的求解精度对涡流扰动的细节计算以及声学的求解结果有着重要的影响。本文应用 LES 方法对增升装置的流场进行数值模拟,采用可穿透积分面的 Ffcows Wil1iams-Hawkings(FW-H)积分方法进行远场噪声计算。采用圆柱绕流算例对本文的数值计算方法进行了验证,验证结果表明：本文所使用的LES 方法能准确地捕捉到涡脱落、流动分离等非定常流动现象,可为远场气动噪声的计算提供精确的近场流动的数值解;基于 FW-H 的声类比方法能够精确高效求解远场气动噪声。在此基础上,对增升装置噪声产生的流动特性、远场特性、风速影响等进行了数值模拟研究。结果表明：缝翼产生气动噪声的主要原因是,流动在缝翼和主翼之间的凹槽形成的不稳定波以及缝翼钝后缘的小脱落涡;襟翼产生气动噪声的主要原因是,襟翼附近由于流动分离产生的高频的小尺度不稳定涡和低频的大尺度涡。     

含激波的湍流流动高精度大涡数值模拟方法

针对采用亚格子模型进行含激波的湍流流动模拟时会面临激波附近的精度损失问题,考虑从通过亚格子模型以及数值模拟方法两方面的改进来实现湍流流动大涡模拟的精度提高.大涡模拟采用了Yee及Sj(o)green (2009)提出的高阶低耗散方法.该方法采用自适应的流场探测器以控制计算中所需区域的数值耗散,并考虑对动力学模型采用在激波位置使用Sagaut和Germano(2005)提出的单边亚格子过滤器和(或)直接禁用亚格子项等方法加以改进.对于标准的马赫数1.5和3条件下的激波-湍流干扰问题,上述新方法相较于全区域采用亚格子模型的方法均表现出了相似的精度提升.同时实现的数值精度改进方案采用了Harten的亚单元分辨过程来定位和锐化激波,并在精确激波位置附近的网格点处采用了单边测试滤波.     

带自由表面水流的分离涡模拟

目前水动力学问题的数值求解仍然广泛采用RANS模型,特别对于大尺度的地表水流的数值模拟,该类模型计算效率较高,但其难以给出湍流场的小尺度涡结构.深入研究复杂地形条件下的水流运动,需要发展诸如LES等的高级数值模型.天然地表水流,地形、边界复杂,雷诺数通常较高,高级数值模拟受限于计算能力的限制,还难以工程应用.RANS和LES混合模型是目前具有工程应用的一种混合数值模型,其已经在CFD领域获得了长足发展.然而对于地表水流运动的数值模拟,相关工作还较少.本文将分离涡模型(DES),即一种RANS和LES的混合模型,应用于带自由表面的地表水流运动,建立了一套数值仿真模型.模型基于有限体积法,水平面内采用非结构计算网格,垂向为结构化网格,对流项离散格式采用二阶TVD格式,并行基于OpenMP语言库.算例表明DES模型有助于揭示复杂地形条件下带自由表面水流的大涡拟序结构.     

不同钝体宽度下的驻涡燃烧室排放试验

为了研究空气流量分配对驻涡燃烧室对排放特性的影响,了解对驻涡燃烧室内污染物生成的过程及其影响因素,设计了一个能够改变中心钝体宽度、仅凹腔供油的驻涡燃烧室.在常压下对该驻涡燃烧室进行了排放特性试验,进口温度保持200℃.试验中,燃烧室进口马赫数为0.15~0.3.影响排放的因素主要包括雾化质量、凹腔当量比以及与进口马赫数相关的驻留时间等.总体来说雾化质量、凹腔当量比的提高对降低CO和HC的排放是有利的,但是这会使NOx排放增加.在低凹腔当量比时,CO排放曲线变化下降比较平缓,甚至出现上升趋势,而HC排放曲线比较陡峭.这是由于HC的消耗速度比CO消耗速度快,随着凹腔当量比的增加,供油压力提高,燃油雾化粒径变小,燃油蒸发时间缩短,使HC排放快速减少,中间产物CO大量产生而来不及消耗.凹腔当量比进一步上升时,由于燃烧温度的提高,使得CO排放快速减少.在燃烧室内燃烧过程中,NOx的形成和消耗是非常复杂的过程,目前只能作一些定性的分析,而CO和HC的反应过程相对简单.通过对不同钝体槽宽下,具有相似凹腔前壁流量的工况的比较,发现CO和HC的形成主要受凹腔内工作状况影响,而NOx的形成过程更复杂,主流也对其产生着重要的影响.      

垂直/短矩起降飞机机翼内埋式风扇布局气动特性分析

以垂直/短距起降飞机过渡飞行状态为背景,针对机翼内埋式风扇布局的自由来流/风扇喷流混合型流动,基于结构/非结构混合网格使用CFD方法进行了非定常数值模拟和分析.首先使用滑移网格技术对NASA涵道螺旋桨进行算例验证,其时均计算结果与实验值的误差为5.3%,证明了计算方法的可靠性和准确性,然后数值模拟了机翼内埋式风扇布局在不同迎角下的气动性能.结果表明:风扇喷流在机翼上产生了特有的“抽吸”和“堵塞”效应,引起了机翼总升阻力的显著增加,升力最大增量达到干净机翼升力的2.6倍,阻力最大增量为干净机翼阻力的3.2倍,混合流场在机翼后缘引起了升力损失并卷起对涡.      

尾缘锯齿降低叶栅噪声的数值模拟

采用大涡模拟与声类比的方法研究了尾缘锯齿对涡轮叶栅噪声的影响.设计了两种不同的尾缘锯齿,对比了Re=3.3×105(基于叶片弦长与叶栅出口速度)下两种不同结构锯齿尾缘叶栅与直尾缘叶栅的声功率.结果表明:尾缘锯齿可以降低叶片吸力面边界层分离噪声约5dB,降低尾缘涡脱落噪声约10dB.进一步的研究表明,尾缘锯齿可以降低叶片尾缘附近表面的压力脉动幅值约50%,将展向相关尺度较大的涡破碎成展向相关尺度较小的涡,并消除尾缘脱落涡,这三者的综合作用使噪声得到降低.      

边界层厚度对腔体气动声学特性影响数值模拟

为了研究来流边界层厚度对开式腔体气动声学特性的影响,基于分离涡模拟方法,计算了来流马赫数为2.0条件下,不同来流边界层厚度与腔体深度比时,长深比为5.88的腔体流动特性,得到了该腔体声压级的频谱特性.计算结果表明:随着来流边界层厚度增加,形成的剪切层稳定性增强,失稳后上下摆动幅度减少,失稳生成的大尺度涡与超声速主流的相互作用减弱,使得大尺度涡发展到腔体后缘时所具有的平动动能和转动动能降低.大尺度涡撞击腔体后缘在腔体内形成的气动噪声的声压级降低,最大减小幅度达7.5dB.同时各阶模态的频率也发生偏移,偏移值在100Hz左右.基于新的假设重新推导了Rossiter公式,明确了经验常数的物理意义,并以此解释了频率偏移现象.