轴流压气机引气后流场的数值模拟

选取了NASA Stage 35作为算例进行压气机引气数值模拟研究.通过对比验证了网格无关性和计算的准确性,应用两种数值实现方法分别进行圆孔和槽引气计算,两种数值方法的结果基本吻合.进行了两种引气方式（即圆孔和槽）的对比.结果表明:槽引气转子出口总压损失小,效果好于圆孔引气.应用NUMECA软件自带引气模块操作简单,容易实现.绘制引气管路网格增加了流道网格绘制的难度和计算边界条件设置的复杂度,但更加符合压气机实际的引气情况.引气后压气机效率 流量曲线整体上移,总压比 流量曲线下降,失速裕度有较大提升.      

适用于倾转旋翼机气动分析的新型嵌套网格方法

为克服围绕复杂倾转旋翼机生成高质量一体化结构网格的困难,将分块网格策略与嵌套网格方法融合,提出并建立了一套适用于倾转旋翼机气动干扰分析的新型嵌套网格生成方法。该方法中,将整体网格按照倾转旋翼机结构特性进行了分块生成,在保证插值精度的基础上降低了网格规模。为高效连接不同部件网格,首先,引入Inverse Map方法来辅助部件网格定位,采用新颖的多向投影方法保证了部件间洞边界的连续性;其次,基于独立挖洞嵌套策略完成机身与背景网格的组合,并使用了相应的辅助计算坐标进行背景边界网格标识。在贴体网格区域采用基于S A（Spalart Allmaras）湍流模型可压缩RANS(雷诺平均Navier Stokes）方程求解流场,背景网格区域采用欧拉方程,并运用了SPMD(single program, multiple data)模式的并行加速技术,建立了适合倾转旋翼机气动干扰特性分析的高效混合CFD方法。采用Robin直升机机身作为数值算例,验证了CFD方法的有效性。在此基础上着重对倾转旋翼机在悬停、过渡和前飞巡航模式下的干扰流场进行了数值分析研究,分别得出了悬停状态下倾转旋翼机典型的“喷泉效应”干扰现象、过渡状态下的非定常气动干扰特性以及巡航状态下的部件气动变化特性,计算结果表明建立的新型嵌套网格方法能够较好地表征倾转旋翼机外形特性,并能够有效地用于倾转旋翼机的气动特性分析,加速比能超过5.0。      

三种T700级碳纤维及其复合材料性能比较

对MT700、T700-A及T700-B三种碳纤维拉伸性能、表面形貌、单向板力学性能及网格加筋圆筒轴压稳定性进行逐级对比研究。结果表明:MT700碳纤维拉伸性能达到同级别进口碳纤维水平且具有高模量特征;MT700碳纤维表面均布沟槽的结构特点使得MT700/603复合材料体系表现出良好的界面性能和拉伸-压缩匹配性,单向板压缩强度、层剪强度及弯曲强度均明显高于T700-A/603和T700-B/603;MT700/603网格加筋圆筒轴压破坏强度及模量分别达到870 k N和108.2 GPa,相比于T700-B/603分别提高11.5%和33.1%。MT700碳纤维更适用于制备航天领域结构复杂承力构件。     

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

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

颗粒无序分布的堆积床内部流动与传热分析

为了详细分析单组元发动机催化床内部的流动与传热,针对颗粒无序分布的堆积床,提出了一种网格生成途径,适合于对大量颗粒自然堆积的管道流动进行数值模拟。根据该方法,对不同颗粒直径、不同进口速度的管道堆积床进行了液体介质的流动与非定常传热模拟,并与实验结果和Ergun方程计算进行了对比。结果表明,本文提出的方法可以有效模拟颗粒无序分布的堆积床内部流动与传热,适合于大量颗粒的堆积床仿真;网格尺度小于1/20颗粒直径时,模拟结果与实验结果符合很好;在相同雷诺数下,摩擦系数随颗粒直径增大而减小,Ergun方程在高雷诺数下计算的流阻偏大;非定常传热时,壁面效应会使壁面附近温度先接近流体温度。     

化学非平衡流动解耦算法的回顾与新进展

超声速化学非平衡流动的数值模拟一直是计算流体力学领域的难点,主要体现在如下几个方面：物理过程非常复杂,存在着激波、燃烧波等各种复杂波系的相互作用;超声速化学非平衡流动属于典型的时空多尺度物理问题,其控制方程存在严重的刚性,给数值求解带来了很大困难。对国内外的解耦算法研究现状简单综述后,主要介绍1993年刘君提出的解耦算法的理论基础,流动方程采用冻结流模型,源项方程模拟流体微团在当地绝热、定容的热力学系统内发生的化学反应过程。通过引入两个中间变量,即等效内能和等效比热比,将与温度无关的生成焓从流动方程组能量项中分离出去,源项方程组中包含等效内能,使用不同算子对流动方程和源项方程解耦求解。与传统解耦算法相比,源项方程的求解过程中包含状态参数和组元同时变化。结合刘君解耦算法机理和有限体积法空间平均特性,介绍近期在提高算法计算效率方面的研究进展,包括流动方程优化算法和耦合过程优化。采用优化算法对经典的激波诱导燃烧算例进行数值模拟,与不同文献结果进行对比,验证了优化算法的时空精度。通过总结经验发现化学非平衡反应仅发生在流场局部区域,提出质量生成率判据,结合相应模拟结果验证了方法的可行性,可以进一步提高化学反应算子的计算效率。     

基于惯性补偿的三分量测力天平瞬态特性的研究

高超声速一体化试验模型在脉冲燃烧风洞中进行测力试验时,模型振动导致天平无法准确测量模型所受到的气动力载荷。为研究脉冲风洞瞬态测力问题,采用了以下方法:根据测力天平的结构特点建立了动力学方程;对其进行了虚拟标定和模态分析;对测力天平进行了瞬态分析并对输出结果进行了惯性补偿。结果表明,单分量阶跃载荷加载时,补偿后输出结果超调量大大下降,振动衰减时间缩短,俯仰力矩会对轴向力输出产生干扰,且三分量阶跃载荷加载时,干扰降低,各分量的超调量分别为162.6N,574.4N和38N·m;单分量正弦加载时,加载分量上输入输出基本一致,俯仰力矩对轴向力的干扰作用仍然存在,3分量正弦载荷加载时,输入输出结果具有相同的周期特性,俯仰力矩结果与输入结果一致,轴向和法向输入输出间存在一定偏差,3分量的超调量分别为24.1N,375.7N和70.8N·m。     

叶轮行波振动下的气动阻尼计算

为展开运动边界下离心叶轮流场的数值分析,独立开发了网格变形程序和非定常流动模拟程序,实现了流场中振动离心叶轮的气动阻尼计算。采用紧支径向基函数法进行结构到气动表面变形的数据传递,应用二叉树技术进行壁面距离的计算,大幅提高了网格变形和流场分析中距离搜索的计算效率。通过振动叶栅和径向叶轮的算例,验证了程序应用于运动边界流场计算和叶轮流场模拟的正确性。以半开式叶轮为对象,展开其在行波振动条件下气动阻尼的计算分析。结果表明,叶轮在前、后行波振动下的气动阻尼均为正,后行波振动引起的气动推力和气动阻尼较大,叶轮表面的平均气动功率密度呈现循环对称的分布特征。     

航空发动机两级反向旋流燃烧室燃烧流场大涡模拟研究

为深入了解真实航空发动机内燃烧流场,采用动态亚网格模型结合单步快速化学反应(FC)、火焰面(FLM)及反应进度变量(FPV)等三种燃烧模型对径向两级反向旋流燃烧室单头部构型进行了大涡模拟。为避免模型简化误差,对燃烧室包括全部气膜冷却孔在内的精细结构进行了完全仿真。在已达到统计定常状态的冷态流场基础上首先模拟了燃料喷注和掺混过程,约2.6ms后燃料到达真实的点火位置,随后采用FPV模型在半径3mm的球形区域数值模拟了点火,展示了在主燃孔横向射流作用下初始火焰沿化学恰当比混合分数等值线传播并充满整个火焰筒的发展过程,结果显示火焰到达燃烧室出口的耗时约为6～7ms。不同模型算法预测的平均温度场与CARS测量结果作了对比,LES-FPV,RANS-FPV,LES-FLM以及参考文献中RANS-FLM计算平均误差分别为3.47%,4.17%,7.76%和11.22%,表明LES改进了模拟精度,且FPV模型显著优于FLM模型。RANS-FPV预测的出口存在严重热斑,导致其给出的出口温度分布因子(OTDF)及最大径向温度分布因子(RTDF)值分别达到0.593和0.313;LES-FPV结果均匀性最好,其预测值分别为0.284和0.193。     

一种部分森林覆盖环境电波传播特性预测方法

针对部分森林覆盖环境中传播媒介的非均匀性,提出一种基于动态网格的抛物方程(PE)预测方法,该方法采用宽角PE法进行大区域电波传播特性预测。为平衡PE求解的准确性与时效性,引入动态网格技术。通过仿真算例,探讨了基于动态网格的PE法在部分森林覆盖环境电波传播特性预测中的有效性。结果表明:在计算精度相似的情况下,相对于均匀细网格,动态网格会大大减少PE所需的计算时间;而相比均匀粗网格,动态网格具有更高的计算精度。