共轴双旋翼非定常流场干扰特性

直升机共轴双旋翼相遇过程的非定常气动干扰在其周期性干扰中最为强烈。为了研究该过程中的干扰特征,建立了一个基于非定常雷诺平均Navier-Stokes方程的气动干扰数值方法,双旋翼的反转运动采用运动嵌套网格进行模拟。为了探究双旋翼相遇时的气动干扰机理,以两个相对运动的双翼型系统来模拟双旋翼特定展向截面的非定常相遇过程。分析了双旋翼和双翼型相遇时的气动特性和流场特征,并对双翼型系统进行了参数影响研究,结果表明:双旋翼上下桨叶相遇时,上下旋翼拉力均会出现先增后减的变化趋势且波动幅值分别为其对应总拉力的30%和22%,双翼型系统的升力波动趋势与之相似;间距增大使上下翼型间的气动干扰减小,且上下翼型升力波动对间距变化的敏感时间不同;翼型相对厚度增大使双翼型升力波动幅值增大,并出现二次波动;有附加来流时,相对来流速度较小的翼型升力波动的幅值与范围更大。      

非结构嵌套网格中的一种改进型径向基函数插值方法

非结构混合网格拓扑关系相对复杂,在嵌套边界区域进行流场信息插值传递容易有精度的损失。利用已知的解析函数作为插值样本,对线性插值(LINE)、距离倒数权重插值(WA)、逆二次径向基函数插值(IQ)、Wendland暞sC2径向基函数插值四种方法的精度进行对比,筛选出精度较高的C2径向基函数;提出一种改进的C2径向基函数插值方法,通过调整径向基函数的作用半径来控制插值矩阵条件数,进而消除使用径向基函数插值在嵌套边界区域出现的数值奇性问题;选用MD30P/30N 多段翼与三维AEDC典型外挂物分离模型对该插值方法进行验证。结果表明:改进后的径向基函数插值方法能够有效消除嵌套边界区域数值奇性,同时计算收敛速率更快,插值精度更高。     

共轴刚性旋翼气动干扰及操纵特性分析

建立了一套基于高效配平策略和嵌套网格技术的共轴刚性旋翼气动数值分析方法。流场求解采用雷诺平均Navier Stokes方程,空间离散采用Roe monotone upstream centered schemes for conservation laws（MUSCL）格式,湍流模型为Spalart Allmaras模型;配平计算中引入“差量修正”策略。分析了共轴刚性旋翼模型在不同飞行状态下的气动干扰现象和操纵特性。结果表明:所建立的配平方法在同等精度下,比传统方法可显著节省计算时间;共轴刚性旋翼间距较小,其气动干扰明显强于普通共轴旋翼;悬停状态下,由于受上旋翼下洗流影响,下旋翼总距较大但拉力较小;小速度前飞时气动干扰引起下旋翼气流不对称性增大,双旋翼周期变距差也变大;采用适当的升力偏置量,可以提高共轴刚性旋翼气动效率。      

矩形/非结构混合网格技术及在二维/三维复杂无粘流场数值模拟中的应用

建立了一套模拟复杂无粘流场的矩形／非结构混合网格技术,其中非结构网格仅限于物面附近,发挥非结构网格的几何灵活性,以少量的网格模拟复杂外形;同时在以外的区域采用矩形结构网格,发挥矩形网格计算简单快速的优势,有效地克服全非结构网格计算方法需要较大内存量和较长ＣＰＵ时间的不足。混合网格系统由修正的四分树（２Ｄ）／八分树（３Ｄ）方法生成。本文将ＮＮＤ有限差分与ＮＮＤ有限体积格式有机地融合于混合网格计算,消除了全矩形网格模拟曲面边界的台阶效应,同时保证了网格间的通量守恒。该算法被成功地运用于ＮＡＣＡ００１２翼型的跨、低超声速绕流,二维弯管内运动激波的非定常流,以及三维捆绑火箭的超声速绕流的数值模拟。计算结果表明本方法在模拟复杂无粘流场方面的灵活性和高效性。     

生成非结构网络的改进阵面推进法

通过对非结构化网格生成的阵面推进方法进行分析,对该算法进行了优化设计,采用笛卡尔直角坐标网格来代替背景非结构化网格,其上网格尺度参数通过控制体有限元法求解Possion方程的Dirchlet问题得到,提高了网格生成的效率和质量。同时利用Laplacian滤波算子和对角交换方法来改善三角形网格的质量,尤其是对角交换技术对消除畸形三角形单元效果明显。对固壁边界采用三次样条拟合后重新离散的方式形成初始阵面,改进了网格的贴体性,提高了网格的光滑性。在编制通用程序时利用链表结构和派生类型的数据变量,节省了计算内存,提高了网格生成效率。     

Euler方程的自适应多重非结构网格计算

介绍一种基于四面体网格的多重网格技术,用于求解三维Euler方程定常解,具有十分明显的加速收敛效果。求解Euler方程的基本算法采用有限体积中心差分格点格式。多重网格技术采用一组几何上相互独立的粗细不同的网格。流动变量、余量和修正量在相邻粗细网格间的相互传递通过线性插值实现。采用自适应技术对网格进行加密,可以获得较优的网格,并作为一套细网格,进行多重网格计算。数值实验表明本文自适应多重网格法十分有效。     

非结构动网格在三维可动边界问题中的应用研究

对用于非结构动网格生成的弹簧近似方法进行了研究。通过采用顶点弹簧方法,分析研究了弹簧倔强系数的取值,同时通过引入挤压倔强系数和边界修正,对标准弹簧近似方法进行了改进。改进后的方法可以大大提高网格变形能力和网格质量。应用本文发展的非结构动网格生成方法并通过耦合求解基于(Arb itrary Lagrangian-Eu lerian ALE)描述的Eu ler方程,模拟了做谐和振动的刚性无后掠机翼及柔性M6机翼的跨音速绕流,计算结果与实验结果吻合良好。     

基于混合网格的三维Navier-Stokes方程并行算法

提出了一种基于混合网格的三维Navier-Stokes方程的并行计算方法。Navier-Stokes的求解采用了基于面的有限体积方法,该方法适用于任何网格类型。采用一方程Spalart-Allmaras模型来计算紊流黏性。并行计算采用区域分裂的方法,利用METIS网格分区系统实现了各节点的加载平衡。节点间的数据交换通过调用MPI库函数来实现,采用非阻断通讯的方式来减少数据交换时间。充分利用FORTRAN90的动态存储特性来减少对内存的需求。最后,通过对绕DLR-F6外形(翼身组合体+挂架+发动机短舱)黏性流动的数值模拟,验证了该并行程序的准确性,高性能并行计算以及处理复杂几何外形的能力。     

三维非结构贴体直角网格的N-S方程数值模拟

基于八叉树数据结构,采用直接分层光顺的方法,解决了非结构直角网格在物面不正交贴体的难点,发展了任意复杂外形的三维非结构贴体直角网格,并且可以根据流场特征进行自适应处理。在此基础上采用最小二乘法对N-S方程进行了MUSCL有限体积法的非结构网格的粘性求解,提出了带约束边界条件的最小二乘法的处理办法。计算表明,三维非结构贴体直角网格的粘性求解,可以适应各种高超声速的复杂外形,计算结果与试验吻合较好。     

HLLC方法在三维非结构网格上的应用

以三维非结构网格的显式有限体积法为基础，采用格心方法以及基于近似黎曼解的二阶精度Godonov扩展方法求解三维Euler方程，使用HLLC(Harten，Lax，van Leer,Contact)近似黎曼解的方法计算网格单元边界处的守恒量通量。为了验证方法的可行性，利用三维解算器对NACA0012翼型绕流进行了数值计算，结果令人满意。在此基础上对某型弹的流场进行了数值模拟，取得了比较好的结果。