基于Ansys Fluent的近场翼尖涡数值模拟与分析

为了进一步研究飞机远场尾涡,提供网格分配及湍流模型的参考,并为整机模拟提供必要的参考依据,通过基于Ansys Fluent的数值模拟方法,研究了NACA0012机翼的近场翼尖涡流场,采用有限体积法求解不可压缩雷诺平均Navier-Stokes方程,其中雷诺应力项分别以S-A和Realizable k-ε模型封闭,模拟了近场翼尖涡卷起的过程,分析了机翼表面压力以及涡核参数,包括轴向涡量、涡核位置、涡核粘性等,并与风洞实验结果进行了对比。结果分析表明:基于局部O-网的六面体网格,RKE模型要优于S-A模型,与实验值更为吻合。     

DES与DDES在湍流分离中的原理与性能研究

随着工程上流动结构的日益复杂，兼具雷诺平均Navier-Stokes(RANS)方法高效率与大涡模拟（LES）高精度的分离涡模拟（DES）类混合方法成为现阶段工程中最有效的湍流模拟方法之一。围绕DES类混合方法中的DES与延迟分离涡模拟（DDES）方法开展工作，分析二者开关函数构造上的不同，研究延迟因子作用机理，并考察DES与DDES方法的求解能力。研究表明：DES与DDES方法在模拟表现上存在一定差异，DDES方法通过引入延迟因子，保护RANS求解区域，改善模化应力不足，降低了DDES方法对交界面系数CDES敏感程度；DDES方法在计算过程中容易出现过度保护，导致求解瞬时涡结构能力不如DES方法，分析与延迟因子引入比重及开关函数构造形式有关。     

旋翼桨尖涡生成及演化机理的高精度数值研究

为了细致捕捉直升机旋翼桨尖涡的生成和演化过程,建立了一个基于高精度网格和高阶通量计算格式的旋翼桨尖涡计算流体力学(CFD)求解方法。在该方法中,流场求解选取旋转坐标系下的Navier-Stokes方程为控制方程;空间离散采用迎风Roe格式,并采用低耗散的5阶WENO(Weighted Essentially Non-Osciltatory)格式进行对流通量的计算;时间推进则采用双时间法,在伪时间步上使用隐式LU-SGS(Lower-Upper Symmetric Gauss-Seidel)推进格式;应用嵌套网格方法实现桨叶网格和背景网格的数据交换。应用所建立的方法对悬停状态的旋翼桨尖涡流场进行了高精度模拟,在桨叶网格上精细地捕捉到了桨尖涡的具体生成过程,在背景网格上捕捉到了更多圈数的桨尖涡尾迹,并对桨尖涡的演化机理进行了深入研究。结果表明:建立的高精度数值方法能够有效地对旋翼桨尖涡的生成和演化过程进行细致模拟;悬停状态下旋翼桨尖气流在上下表面压力梯度的作用下经历了边界层增厚、逐渐卷起形成涡核以及最终脱离桨叶形成桨尖涡的过程。     

气液针栓式喷注器液膜破碎过程和结构的数值研究

为了更好地认识针栓式喷注器雾化场的结构，基于网格自适应加密技术以及VOF (volume of fraction)方法追踪气液的分界面，采用realizable k-ε湍流模型模拟整个流动过程，还原了不同时刻气/液撞击的初次破碎过程，数值模拟结果与高速摄影试验结果定性定量对比均吻合较好，验证了数值方法的准确性。进一步对针栓式喷注器气/液撞击的初次破碎过程、内部流场涡结构、速度场进行分析，研究了初次破碎雾化的动力学过程和机理。研究结果表明：液桥的形成主要是由液洞的扩展和拉伸、合并而形成，而液滴主要是由中心液膜拉伸、液丝断裂以及液桥断裂而形成，液膜破碎阶段形成的涡结构是造成液膜断裂的主要原因。     

涡轴发动机双回路PI控制器多发功率匹配

为了应对多个发动机共同驱动同一直升机时,单个发动机性能衰退所引发的输出功率不平衡问题,以控制2个转子转速的双回路结构PI控制器为基础,将外回路改为直接功率控制回路,搭配能够计算旋翼在一定转速下需求功率的机载模型,构建了1种涡轴发动机多发功率平衡匹配控制系统.内回路分别采用燃气发生器转子转速控制回路和动力涡轮转子转速控制回路,得到参数不同的控制器并进行了仿真验证和对比.结果表明:所设计的双回路PI控制器能够在保证涡轴发动机动力涡轮转速恒定的同时,使性能衰退程度不同的2台发动机输出相同的功率.     

三维超声速不等密度平面混合层的小激波形成与演变

基于大涡模拟方法与高精度混合格式,对两种密度比(1.0与0.138)条件下的高对流马赫数(Mac=1.1)三维空间发展超声速平面混合层进行了数值模拟。数值结果清晰描述了扰动激励下超声速平面混合层的演变过程,讨论了三维混合层中小激波的结构及其形成过程。分析了混合层两侧密度变化对混合层失稳及小激波形成的影响,发现当上、下层密度不等时,混合层失稳加速,小激波强度减弱,而密度较低侧的小激波会更弱。     

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

为了揭示近失速工况下动叶间隙和静叶角区的复杂流动结构,探索诱发压气机失速的主要因素,对高负荷跨声速压气机开展数值研究,获得了整级压气机和单转子的特性曲线,进而对压气机的流场进行了详细分析。结果表明:对于研究对象,失速最先起始于静叶上角区。上角区完成由开式分离到闭式分离的转变过程,可以认为是静叶发生失速的标志。在近失速工况下,动叶叶尖泄漏涡发生"泡式破裂"堵塞流道,泄漏涡破裂引起的堵塞作用一直从叶顶延伸至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 的声类比方法能够精确高效求解远场气动噪声。在此基础上,对增升装置噪声产生的流动特性、远场特性、风速影响等进行了数值模拟研究。结果表明：缝翼产生气动噪声的主要原因是,流动在缝翼和主翼之间的凹槽形成的不稳定波以及缝翼钝后缘的小脱落涡;襟翼产生气动噪声的主要原因是,襟翼附近由于流动分离产生的高频的小尺度不稳定涡和低频的大尺度涡。