基于非结构网格的不可压N-S方程多矩有限体积法

提出了一种基于三角形及四面体非结构网格的有限体积法(FVM),用以鲁棒且精确地求解不可压粘性流动问题.与传统的FVM方法仅将体积分平均值(VIA)作为计算变量的做法不同,本文提出的方法将VIA及点值(PV)同时作为计算变量并在每个迭代步进行计算更新.VIA以通量形式进行计算以确保数值守恒,PV可以通过控制方程的不同形式进行求解更新,无需守恒,因此可以采用非常高效的方法进行求解.将PV作为增加的变量使得紧致网格模板得以实现更高阶精度的重构,而且由此获得的数值模型对于非结构网格变得更鲁棒.本文针对二维/三维的三角形/四面体非结构网格提出了数值格式,给出了几个基准测试算例,验证了本文提出的数值方法在采用非结构网格求解不可压粘性流动问题时的精确性和鲁棒性.     

跨声速串列转子失速机制的数值研究

为了解决跨声速串列转子的低裕度问题,就必须了解跨声速串列转子的流场结构与失速机制.设计了叶尖切线速度为450m/s,负荷系数为0.56的高负荷跨声速串列转子.基于数值模拟的结果,分析了该串列转子在0.5mm叶尖间隙下的叶尖流场结构与失速机制,并在此基础上分别探讨了叶尖间隙和前、后排叶片周向位置对串列转子特性的影响和失速机制的变化.结果表明:前排叶片的叶尖区域是 影响串列转子稳定性的关键;随着叶尖间隙的增加,串列转子的失速机制也发生变化,从前排叶片叶尖区域的尾迹与径向潜流堵塞后排叶片通道转变为前排叶片叶尖泄漏流堵塞;在较大周向相对位置(后排叶片压力面周向远离前排吸力面)的情况下,串列转子获得最好的效果,随着周向相对位置(PP)的增加,失速部位从后排叶片转移至前排叶片.      

蜂窝式轴心通风器油气分离性能计算

为对航空发动机蜂窝式轴心通风器油气分离效率进行研究,建立考虑油气双向耦合的流场计算方法及油滴/壁面相互作用模型,在验证通风阻力及油气分离效率可靠性的基础上,对不同转速、通风流量和环境温度下蜂窝式轴心通风器的油气分离效果进行计算和分析.结果表明:转速的增加会使油气分离效率得到提升,而通风流量和环境温度的增加则导致油气分离效率的降低.蜂窝孔结构的加入对通风阻力影响不大,却对通风器的滑油分离过程起主要作用,计算表明其对滑油分离贡献率在80%以上.      

New constitutive relation for eddy viscosity models

A new turbulent constitutive relation was directly derived from Boussinesq's hypothesis and mixing length theory, and then implemented in the standard k-ε model. The performance of this constitutive relation was validated in zero pressure gradient flat-plate boundary layer flow, fully-developed turbulent channel flow and separated flow in a plane asymmetric diffuser. The investigation demonstrated that, this new constitutive relation gave very accurate results in the former two basic cases and provided significant improvement in prediction of separated and reattachment points in the plane asymmetric diffuser. Separation and reattachment points at x/H=7.5 and 29 were calculated accurately in comparison to experimental results, and the static pressure coefficient of 0.82 was very close to large eddy simulation calculation. These results are very encouraging but further verification and extensive application of the new constitutive relation to other two-equation eddy viscosity model are needed.      

温度对飞机燃油射流泵工作特性影响

建立了射流泵数学模型,采用数值计算方法对其工作特性进行了模拟,搭建了实验台,采用RP-3燃油,通过实验数据验证了模型的正确性.在此基础上,分别选择了RP-3和JET-A两种燃油,计算了从-40℃~40℃温度下,不同操作压力对射流泵流量比的影响.结果表明:燃油温度高于0℃时,温度对射流泵工作特性影响并不明显,随着温度的降低,射流泵性能与标准温度下（20℃）偏差越来越大,且压比越高此偏差越明显.当温度高于0℃时,采用RP-3和JET-A燃油性能相差不大,RP-3燃油性能略优于JET-A燃油,但是在低温下,RP-3燃油性能远高于JET-A燃油.因此当燃油温度较低且压比较高时,必须充分考虑温度对射流泵特性的影响.      

Effect of blade tip winglet on the performance of a highly loaded transonic compressor rotor

The tip leakage flow has an important influence on the performance of transonic com-pressor. Blade tip winglet has been proved to be an effective method to control the tip leakage flow in compressor, while the physical mechanisms of blade tip winglet have been poorly understood. A numerical study for a highly loaded transonic compressor rotor has been conducted to understand the effect of varying the location of blade tip winglet on the performance of the rotor. Two kinds of tip winglet were designed and investigated. The effects of blade tip winglet on the compressor over-all performance, stability and tip flow structure were presented and discussed. It is found that the interaction of the tip winglet with the flow in the tip region is different when the winglet is located at suction-side or pressure-side of the blade tip. Results indicate that the suction-side winglet (SW) is ineffective to improve the performance of compressor rotor. In addition, a significant stall range extension equivalent to 33.74% with a very small penalty in efficiency can be obtained by the pressure-side winglet (PW). An attempt has been made to explain the fundamental mechanisms of blade tip winglet in detail.     

水下超声速射流对上浮水雷受力特性影响研究

固体火箭发动机具有功率密度大、推力大等优势,常被用于上浮水雷的推进器。然而发动机在水下工作时燃气射流使流场压力发生剧烈脉动,进而影响发动机的推力性能及上浮水雷的受力特性。基于VOF多相流模型和理想气体模型,建立了上浮水雷在不同工作状态(欠膨胀、完全膨胀、过膨胀)下的数值模型,研究了水下燃气喷射流对上浮水雷的受力特性影响。结果显示,欠膨胀工况时,发动机推力平缓,大小为12.2KN,上浮水雷受力未出现负值;完全膨胀时,射流发生颈缩、胀鼓现象位置距离雷体较远,导致发动机推力及水雷壳体受力振荡不剧烈;过膨胀工况时,射流发生颈缩、胀鼓现象距离喷管较近,发动机推力发生剧烈脉动,产生21.37%的振荡幅度,胀鼓现象发生时,流场压力显著降低使得上浮水雷后体受力减小,壳体阻力增大,上浮水雷最大产生27KN的负推力。     

基于自适应网格的超声速液体射流破碎过程研究

为研究超声速气流中液体横向射流的破碎过程,采用脉冲背景光方法和VOF方法开展了实验和数值研究。为提高液体横向射流中气液界面和气流场特征捕捉的精确性,采用自适应网格技术对于气液界面、激波出现位置进行网格细化,计算得到了较为精细的气液界面、激波特征及涡系结构。研究结果表明:在低成本仿真模拟条件下,利用自适应网格计算得到的射流轨迹和轮廓与实验吻合较好,射流轨迹的最大误差为10%;射流初始段在超声速气流条件下,仍然存在一段高度约为1.9倍喷孔直径且圆柱形态保持较好的连续光滑液柱。随着喷注压降的升高,液柱的长度逐渐增大;主流气体流经液柱发生三维绕流,在射流附近和近壁面区域形成不断演化的反转涡对,反转涡对的形成加速了液体射流一次破碎过程。     

基于功能脆弱性的空中交通相依网络流量分配

依据空中交通管理与航班运行规则,采用复杂网络理论构建由机场、航路与管制扇区组成的相依网络模型,建立不同扰动策略的影响规则,提出以网络流量熵和交通流损失比变化率为指标识别网络功能脆弱性。并以网络总流量熵最小为目标,建立基于改进遗传算法的网络流量协调分配策略,以降低空中交通相依网络的脆弱性。以民航华北地区空域为原型,发现了其相依网络脆弱性表现规律和脆弱源,采用遗传算法求解网络流量分配方案,优化结果降低了网络总熵值和功能脆弱性,其中机场网络流量分配后效果最为显著,验证了方法的有效性,研究结果可为空中流量管理决策提供一定的理论支撑。     

基于RANS/LES混合方法的分离流动模拟

飞行器在大迎角、快速俯仰机动时,流场中含有大尺度、非定常的涡结构,传统雷诺平均Navier-Stokes （RANS）模型不能准确模拟流场结构,根据国际上相关研究的发展趋势,需要采用混合RANS/大涡模拟（LES）模型来对复杂分离流动进行准确模拟。本文对基于分区混合与湍流尺度混合的双重RANS/LES混合计算模型进行发展与应用。通过典型简化模型的静、动态湍流大分离流动,测试和验证所采用的脱体涡模拟（DES）类方法,重点研究改进的延迟DES （IDDES）模型在动态问题应用中的正确性和有效性,并对所采用的数值模拟方法和相应的计算软件的可靠性、鲁棒性以及精度进行了考核验证。典型算例包括超声速圆柱底部流动、跨声速方腔流动、NACA0015机翼深失速分离涡模拟等。计算表明：发展的IDDES类混合计算模型可有效解决对数层不匹配的问题;对于定态非定常分离流动,DES、DDES、IDDES等模型计算结果差别不大,随着流动的非定常特性增强,IDDES模型的优势逐渐显现;对于动态非定常分离流动,则需要采用IDDES类模型。