基于径向基函数与混合背景网格的动态网格变形方法

发展了一种基于径向基函数(RBF)插值和混合背景网格映射的动态网格变形方法。混合背景网格由运动物面附近的各向异性网格和偏离运动物面的各向同性网格组成。物面变形或运动后,包含在各向异性背景网格中的所有物面网格点的位移由指定的物面运动规律得到,而各向同性背景网格点的位移通过径向基函数插值获得。然后变形后的计算网格利用计算网格与混合背景网格之间固定的映射关系代数插值得到。各向同性网格距离物面一定的距离,消除了物面保形对径向基函数插值的要求,从而可以减少径向基函数基点的数量,以提高插值的效率和可靠性。最后,通过二维和三维算例对该方法的变形网格质量、参数影响、网格拓扑依赖性、复杂外形变形能力及网格变形效率进行了测试。测试结果表明,该方法能够高效生成大变形下高质量的变形网格。     

高超声速导弹等离子体合成射流控制数值研究

快响应控制技术已成为高超声速飞行器发展的关键技术之一,具有极快响应、零质量特性的等离子体合成射流(PSJ)已在超声速流动控制方面初步显示出优异的控制能力,极有潜力应用于高超声速飞行器的快响应控制。基于等离子体合成射流的快响应特性,提出了高超声速飞行器等离子体合成射流快响应控制技术,并通过建立简化的高超声速导弹流场控制模型,对等离子体合成射流控制高超声速导弹进行数值研究。首先,理论分析了高超声速导弹流场的典型结构特征,导弹流场中存在3个特征流场结构。在此基础上,在导弹3个特征位置前面安装等离子体合成射流激励器,研究等离子体合成射流对高超声速流场结构的控制作用,分析由此导致的导弹表面压力分布、升阻特性以及俯仰力矩特性变化。数值仿真结果表明:等离子体合成射流对高超声速导弹外流场中膨胀波和斜激波都具有控制作用,使得波的强度均变弱,且对斜激波的控制效果更为显著;导弹流场结构及气动特性变化具有很强的射流跟随性,射流作用下的导弹流场变化响应时间非常短,仅为0.2 ms;通过合理布置等离子合成射流激励器的位置,可以使得导弹表面压力分布快速改变,从而实现高超声速导弹姿态的快速控制。     

一种适用于三维混合网格的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求解器的嵌套网格计算能力有待进一步挖掘和验证。计算还发现,不同湍流模型得到的表面压力分布、分离区大小以及最大升力系数存在差别,而风洞模型支撑和机翼静气弹变形对飞行器力矩特性预测影响较大。     

基于DQN的反舰导弹火力分配方法研究

反舰导弹作为海上作战的主战武器,由于其精度高、射程远、威力大等特性长期以来一直被当作舰艇编队的主要防御对象。针对反舰导弹打击舰艇编队的火力分配问题,我们提出了一种基于深度Q值网络求解反舰导弹火力分配策略的算法。不同于现有的基于领域知识的方法,深度Q值网络无需依赖任何先验信息,就能够通过与模拟器的交互自动求解最佳的攻击策略。该算法使用深度神经网络拟合Q值函数,解决了传统强化学习中的状态空间过大无法遍历的问题。实验结果表明,在各种不同的舰队防御配置下,深度Q值网络求解得到的攻击策略均获得了最佳的毁伤效果。     

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

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

吸气式空空导弹外形多学科一体化优化设计

针对采用整体式固冲发动机的吸气式空空导弹外形气动与推进耦合的推阻匹配设计难题,引入多学科优化设计方法提出了一种综合考虑气动/推进/质量/弹道的导弹外形多学科一体化优化设计技术。其中,气动性能预测采用代理模型技术,主要基于外形参数化建模、非结构网格技术和流场精细数值计算来自动构建气动数据库,据此建立了包含外形几何信息的气动预测代理模型,并对其预测精度进行了验证;推进性能预测采用推进求解模型,该模型根据固冲发动机理论建立,精度满足工程要求。对所建立的学科预测模型完成一体化集成后,以质点弹道仿真评估的战技指标为优化目标,对一款吸气式空空导弹进气道和翼面外形进行了优化设计,取得了推阻匹配的优化外形,优化后导弹动力射程提高10%。所提出的一体化优化设计技术,有助于吸气式空空导弹外形气动与推进耦合推阻匹配设计和提高导弹动力射程。