一种适用于三维混合网格的GMRES加速收敛新方法

为提高流场计算收敛效率,发展了一套适用于三维混合网格Naiver-Stokes方程求解的并行广义最小残差（GMRES）隐式时间推进方法。该方法由科学计算可移植扩展工具包（PETSc）中的Krylov子空间求解器实现,线性方程系统中的系数矩阵直接以显式给出以提高算法的稳定性。为进一步提高GMRES方法的收敛速度,对非结构网格的序号进行了重排序,使得系数矩阵的非零元素尽量向主对角线靠近。利用所发展的GMRES方法,完成了对ONERA-M6机翼、AIAA阻力预测会议通用研究模型（CRM）等算例的计算,计算结果与试验结果吻合良好。通过与其他隐式推进方法进行比较,对算法的收敛特性进行了研究。结果表明,所发展的GMRES方法计算更加稳定,残差下降速度相对LU-SGS（Lower-Upper Symmetric Gauss-Seidel）方法更快,尤其是气动力系数向着收敛解逼近的速度更加明显,提高了计算效率。     

基于非结构网格流场超大规模并行计算

大规模并行的计算流体力学已成为现代航空工业研发的核心手段之一。基于非结构混合网格和有限体积法,发展了适用于工业级复杂外形气动计算的并行流动数值模拟方法。文中首先介绍了紧致数值离散格式、基于Metis的分布式多核系统网格分区技术、并行边界虚拟单元技术和MPI并行实现等相关算法。采用网格量相对较小的旋成体构型绕流模型对比分析多核并行计算结果与单核计算结果以验证并行计算的正确性,比较了不同并行规模下并行效率和残差收敛情况。然后通过对上亿网格单元的运输机复杂构型绕流模拟,开展并行效率的测试,结果表明,本文方法并行加速性能高,直到多达18816核并行效率都保持在80%以上。     

基于非结构混合网格的CHN-T1标模气动特性预测

本文使用自行研制的基于非结构混合网格的亚跨超声速流场解算器MFlow,对AeCW-1提供的客机标模CHN-T1进行了数值模拟研究。介绍了非结构混合网格的生成情况,重点分析了网格收敛特性、压力分布、气动特性曲线、流动分离等。计算得到了近似线性的网格收敛特性。随着网格加密,对激波和分离气泡的模拟更精细。尾支撑对气动特性的影响非常明显,特别是对平尾气动特性有很大影响。机翼静气动弹性变形的影响主要是使升力系数和阻力系数减小。湍流模型的QCR修正对大迎角计算结果有较大的影响。计算结果表明,MFlow程序能够准确地预测客机标模的气动特性。     

基于嵌套网格的CHN-T1标模数值模拟

为响应第一届航空CFD可信度研讨会(AeCW-1),对CHN-T1标模翼身-尾翼组合体构型采用两种RANS求解器进行了基于结构化嵌套网格的网格收敛性研究和抖振特性计算,并与对接网格结果进行了对比分析。参考AeCW-1网格生成指南自主研发了一族以约3.3倍规模增长的稀-中等-密-特密嵌套网格系列。相比于分区点对接网格,融合了贴体网格和笛卡尔网格优势的嵌套网格从根本上降低了网格生成难度,单块网格拓扑形式更为合理且正交性较好,同时能够有效地平衡近场和远场的网格量。对于跨音速运输机构型,同等规模下嵌套网格的物面网格密度较对接网格更大,能够对激波和分离等复杂流动进行更理想的预测。对比计算过程中,自研求解器OFS3D表现出了较高的嵌套网格计算效率和可靠性,而NASA CFL3D求解器的嵌套网格计算能力有待进一步挖掘和验证。计算还发现,不同湍流模型得到的表面压力分布、分离区大小以及最大升力系数存在差别,而风洞模型支撑和机翼静气弹变形对飞行器力矩特性预测影响较大。     

无人机类脑吸引子神经网络导航技术

当前无人机在非结构化或未知环境下飞行主要采用SLAM进行导航与定位,存在如下突出问题：依赖高精度昂贵激光雷达等环境感知传感器;需要建立准确世界和无人机物理模型;受环境影响较大;自主智能水平较低,无法较好地满足无人机对导航系统的要求,需要发展自主智能的导航方式。基于吸引子神经网络的类脑导航技术,无需训练模型参数,不依赖高精度传感器,无需精确建模,且复杂环境下鲁棒性较强,具有解决上述问题的潜力。简要阐述了动物大脑导航机理,分析了吸引子神经网络和基于吸引子神经网络的类脑导航关键技术,最后讨论了吸引子类脑导航技术在无人机应用中的挑战。     

共轴刚性双旋翼/机身干扰流场数值模拟

采用滑移网格技术求解Navier-Stokes (RANS)方程的方法,研究了共轴刚性双旋翼/机身的干扰问题。通过Caradonna-Tung旋翼、Robin直升机、Maryland直升机旋翼/机身干扰和Harrington 2共轴双旋翼等算例,验证了所提出的旋翼流场数值模拟方法的正确性。在此基础上,以Maryland机身为原型,分析了不同桨距的共轴刚性双旋翼与机身之间的干扰特性。结果表明:所提出的数值模拟方法能够很好地模拟共轴刚性双旋翼/机身的气动干扰特性;由于机身对于共轴刚性双旋翼下洗流场的阻滞作用,旋翼的悬停效率增加5%左右,并且随着拉力系数的增大使得悬停效率的增量更加明显;旋翼的悬停效率增加主要来源于下旋翼0°方位角附近的桨叶升力系数的增大,并且拉力系数的增量由桨根向桨尖方向逐渐减小。      

涡桨飞机螺旋桨滑流气动干扰效应及流动机理

螺旋桨滑流对飞机各气动部件的干扰是涡桨飞机气动设计中面临的难点之一。以某双发涡桨支线客机为对象,采用数值模拟手段分析了滑流对机翼、平尾、垂尾的影响及其流动机理。计算采用基于动态面搭接网格技术的非定常方法,气动力与表面压力分布的计算结果均与实验值吻合良好。研究结果表明：在滑流影响下,全机升力和阻力有所增加,升阻比和纵向静稳定度有所降低,并在无侧滑状态下产生了滚转力矩和偏航力矩;在不同展向位置,滑流对机翼表面流动分离的影响存在显著的差异。在螺旋桨下行运动一侧,滑流的旋转减小了当地迎角,同时桨后气流速度较高,翼面流动分离被有效抑制。在螺旋桨上行运动一侧,滑流的旋转增大了当地迎角,而且桨后气流速度较低,因而翼面流动分离并未得到明显改善;在中小迎角下,滑流对背景飞机平尾当地的动压没有产生影响,但增加了下洗角变化率,进而导致平尾效率降低;垂尾的侧力与偏航力矩是由滑流的侧洗引起的,而滑流的侧洗又与左右两侧机翼升力分布的不对称性有关。     

镍铝涂层晶粒细化与金属氧化物掺杂改性研究

包埋渗铝获得的镍铝涂层是一种最早使用的Al_2O_3膜热生长型高温涂层。自20世纪50年代应用于航空发动机热端部构件的高温防护以来,进一步提高其抗高温氧化性能的机理和技术研究延续至今。基于对合金氧化及Al_2O_3膜热生长机制的理解,提出了晶粒细化与特定金属氧化物掺杂可提升镍铝涂层抗氧化性能的观点,介绍了涂层晶粒细化与金属氧化物弥撒掺杂方法,讨论了这些结构和成分改性影响涂层抗高温氧化性能的关键因素:包括Al_2O_3膜生长速度、亚稳态相向稳态相转变、涂层的黏附性以及涂层与合金基体的互扩散。这些新的研究结果有望为进一步挖掘渗铝涂层的应用潜力、延长其服役寿命提供理论和试验基础。     

流向强迫作用下的液体初始雾化机制及动力学特征

使用Volume of fluid (VOF)方法和基于树形数据结构的自适应算法来研究射流雾化的破碎过程以及扰动对射流破碎机理产生的影响。在无扰动情况下,液体射流的头部、液丝和液滴随着射流时间的发展不断发生演变。射流头部先呈现蘑菇状外形,随后液丝生成,并慢慢转变成网兜状,直至断裂形成小液滴。在周期性流向强迫的作用下,射流液柱的表面会形成周期波,其液丝破裂形成液滴的时机与稳定射流情形相比会有所提前,射流形成的头部更趋于扁平,最终生成的液滴数量更多。低中频阶段,随着扰动频率的增大,射流未扰动液柱长度(L)逐渐缩短,液滴直径的概率密度分布(PDF)趋于尖锐,液滴平均直径(SMD)增大。在高频阶段,随着扰动频率的增大,L会随之增大,液滴直径的PDF分布变得平缓,SMD会减小。     

涡轮叶片晶粒细化及热等静压技术的应用

论述了涡轮工作叶片晶粒细化方法的机理,温度参数的影响、铸造工艺及其试验效果,简要阐述了热等静压和热处理的应用技术及工艺参数等对叶片晶粒度、贫铬层、显微组织、力学性能及其尺寸的影响。