机翼-机身-挂架-短舱复杂流场数值模拟

利用超级计算机资源求解雷诺平均N-S方程（RANS）,计算评估大量外形方案性能,在可接受的时间周期内完成民机飞发集成构型下机翼多目标优化设计。搭建了集成机翼CST参数化、复杂外形网格变形、快速CFD流场解算和自动后处理等关键环节,用遗传算法全局寻优的优化系统。合理的流程设计使该系统可在“天河2号”超级计算机上同时对数百个方案实施计算评估。在并行计算、加速收敛等技术的综合运用下,使用包含800万单元的多块结构化网格,对NASA CRM（Common Research Model）机翼/机身/短舱/吊挂构型的计算分析可在15 min内结束。在3点3目标优化案例中,用90个设计变量表达CRM机翼9个控制剖面的中弧线和扭转角,60 h内完成了超过10000个外形方案的计算分析,遗传进化40代。与初始外形相比,PARETO前缘上选择的最优解的各设计点取得了2~10 count（1 count=阻力系数0.0001）的减阻效果。     

物理质团法关于可压流计算的程序验证

为了评估无网格物理质团法(PECM)对可压流计算的可信性,结合物理认识对该方法的设计思想和数值格式进行了先验分析,基于9个一维模型数值解与参照解的对比信息,通过后验分析方法对PECM进行了程序验证,参照解包括无网格再生核粒子法(RKPM)的数值解以及精确解。验证结果表明:各种物理参数以及状态方程的强烈间断会导致RKPM严重的失稳或失真,而PECM在包括极大密度比在内的各种强间断条件下仍保持稳定和收敛并具有较高的准确度。      

气动弹性计算中网格变形方法研究进展

网格变形是气动弹性计算中实现计算网格随运动边界变形的主要方法。在总结目前网格变形方法发展现状的基础上,对近几年常用的网格变形方法,即弹簧法、弹性体法、超限插值法、Delaunay背景网格法、径向基函数插值法和温度体法等做了简要的总结。根据各方法构建模型的不同,将它们分成物理模型法、数学插值法和混合方法3类,简要介绍了各方法的基本思想和研究进展,重点比较了各方法的网格变形特性（变形能力、变形质量和变形效率）和优缺点。总结了当前流场边界与结构边界的数据传递方法,对当前气动弹性计算中遇到的网格变形的难点问题作了简要评述并对未来网格变形方法的发展方向进行了探讨。     

一种基于有限元法和弹性接触理论的齿轮啮合刚度改进算法

提出了一种将有限元法和弹性接触理论相结合的齿轮啮合刚度计算方法.该方法利用子结构法提取齿面原始柔度矩阵并分离出接触点弯曲变形,根据线弹接触变形解析公式计算接触变形,通过求解非线性变形协调方程得到齿轮时变啮合刚度和齿面载荷分布.以一对齿轮副为例,计算的啮合刚度与航空标准计算结果相差在6%以内.该方法发挥了有限元法在预测物体整体变形方面的优势,同时结合弹性接触理论能够准确计算局部接触变形的优点,与常规有限元法相比,能够有效地提升计算效率.由于接触变形问题的非线性,啮合刚度随总啮合力增加呈现非线性增大的趋势.      

分区对接网格在跨流域气体运动论统一算法的应用研究

为了模拟跨流域复杂物形绕流或多体流场,需要在位置空间建立能描述整个流场特性的网格系统,相应的数值计算方法也需要有所改进。本文采用分区对接网格处理复杂物形绕流或多体干扰流场,得到能准确描述流场特征的贴体网格系统。在此基础上数值求解考虑转动非平衡影响的Boltzmann-Rykov模型方程,在速度空间应用离散速度坐标法、位置空间采用数值差分NND格式。网格对接面上考虑分布函数的传递,包含网格点信息和网格点上速度的信息。通过计算得到复杂物形的流动特性。相关算例表明网格分区之间的流场参数光滑过渡。与文献及DSMC结果的对比证实本文基于分区对接网格的气体运动论统一算法是可靠的。     

流体力学方程的间断有限元法

本文对近三十年来,国内外对于高精度数值方法研究中的热点——间断伽辽金方法在可压缩流数值模拟方面的应用研究进行了综述。首先对间断伽辽金方法的基本概念和特点作了简单介绍,然后对应用该方法解决双曲型及椭圆型问题的发展历程进行了回顾,并重点梳理了其在计算流体力学领域可压缩流数值模拟方面的应用发展以及研究现状,之后对该方法在对应的网格技术、激波捕捉方法、湍流流动模拟以及计算量需求方面目前仍然存在的研究难点和可能的发展趋势做出了总结和分析。最后给出了间断伽辽金方法在可压缩流数值模拟中的若干应用实例。     

倾转旋翼机多部件对机翼气动干扰的分析及优化

为研究倾转旋翼机各部件对机翼的气动干扰效应,建立了一套适用于倾转旋翼机流场CFD求解方法,提出并生成了一套由三棱柱/四面体组成的适用于倾转旋翼机多种飞行状态的非结构混合网格系统.以三维Navier-Stokes方程为主控方程,采用动量源方法进行倾转旋翼的模拟,并引入了高效的OpenMP并行加速技术等,提出了一套新型的动量源项添加及搜索方法.计算分析了倾转旋翼机旋翼/机身/短舱对于机翼气动特性的干扰影响,得出了一些对设计有指导意义的结论.采用基于径向基函数（RBF）的代理模型方法,根据机翼不同展向位置的干扰程度对各段翼型配置进行气动优化.优化结果表明:考虑气动干扰作用,在巡航速度下优化后的全机升阻比增长达到了36.78%.      

高雷诺数下的翼型绕流LBM数值模拟

采用格子Boltzmann方法(LBM)的二维9速度(D2Q9)模型和贴体网格,通过引入非均匀网格插值方法和非平衡态外推边界处理,分别结合Baldwin-Lomax(B-L)湍流模型和Spalart-Allmaras(S-A)湍流模型,对高雷诺数Re≥5×105下的NACA0012翼型绕流进行了数值模拟和对比研究,两者的结果与CFL3D的结果和实验结果均吻合的很好,相比之下,采用S-A模型能更好地预测失速迎角,其处理分离流动的能力要强于B-L模型。改进后的LBM适用于非均匀贴体网格,曲边边界,计算简单,并可应用于更复杂的高雷诺数流动中。     

基于动态边界控制的优化方法在高升阻比气动外形优化中的初步应用

本文是利用动网格技术和有限体积法来研究基于动态边界控制的气动优化的方法.该优化方法是通过在求解非定常流动的过程中计算优化的目标函数随设计变量变化的梯度.而利用参数化的结构动网格技术能有效、快速地生成与优化过程出现的外形对应的网格并进行计算,具有较高的效率.本文以典型的高升阻比外形优化为例,用非定常NS方程为流动控制方程进行优化,取得了初步的结果.     

自力式调压阀变开度流场及特性分析

以研究自力式调压阀性能为目的,利用动态网格方法对阀门内流场进行数值模拟。通过对不同入口流量和弹性系数的多个算例分析得到:稳定压力和入口流量呈直线关系;稳定压力、阀门开度和弹性系数呈二次曲线关系。调压阀在入口质量流量变化时,阀门可以很快稳定,实现对压力和流量进行的调节。