大型客机涡扇发动机动力特性模拟

机体/发动机干扰问题是现代大型客机设计中必须考虑和解决的核心问题之一。长期以来,商用和in-house计算流体力学（CFD）软件,都是通过在发动机进气口设置质量流量比,排气口和外涵道设置总温比/总压比,建立动力特性模型来等效模拟机体/发动机干扰流场。在现有动力特性模型基础上,借鉴特征边界思想,将外涵道指定为特征边界,建立了一种新的发动机动力特性模型。采用轴对称超高涵道比涡扇发动机模型、轴对称涡轮动力模型以及某型客机带发动机模型对两种动力特性模型进行了系统的验证和确认,结果表明：两种发动机动力特性模型均能很好地模拟涡扇发动机动力效应,且计算结果与试验数据吻合较好,说明了两种动力特性模型的正确性、可靠性以及工程适用性。此外,所建立的新动力特性模型特别适合于大涵道比涡扇发动机的动力效应模拟,且可以在外涵道总温比、总压比等参数未知的情形下,评估发动机的动力特性,在实际工程中应用更广泛。     

利用二、三维嵌套技术数值模拟复杂边界下的流场

通过数值模拟复杂几何区域和边界条件下二、三维嵌套紊流模型,获得了计算区域内任意一点的流场信息。这一方法克服了常规的单纯模拟二维或三维紊流控制方程组而无法获得准确的自由表面和三维紊流场的缺点,通过两者有机的结合,并采用动边界处理等先进的计算技术,得出了精度更高、更为合理的紊流场。该方法曾被用于三峡工程大江截流完成之后导流明渠通航的流场模拟中,计算结果合理可信,部分成果的精度也已得到了实测资料的验证,成果为导流明渠通航路线的选择提供了依据     

粘性非结构网格生成方法及绕流计算

运用非结构网格进行计算时,需要在物体附面层区域内生成高展弦比的三角形网格。本文给出了一种基于物理学中弹簧原理的粘性非结构网格生成新方法－－挤压法,它能有效地实现粘性区域内位伸网格到附面层外一般非拉伸网格的光滑过渡。然后在非结构粘性网格上对Ｎ－Ｓ方程运用有限体积法进行了数值计算,ＮＡＣＡ－００１２翼型以及平板层流的数值试验结果较为满意。     

一个解析三维非结构网格生成方法

给出了一个三维非结构网格生成方法。提出弱区域指示函数描述三维封闭流场形状,并证明了该函数的存在不唯一性。利用该函数基于Delaunay方法依据流场尺度生成网格结构。其中,边界引力有效解决了流场边界恢复和内嵌边界问题,弹簧振子和正四面体趋进技术有效调整网格结构并提高网格质量。最后,弹簧振子和边界引力随着弱区域指示函数的演变与监测函数的指导成功应用到自适应网格生成。     

超声速条件下内埋式武器分离特性的数值分析

采用非结构动网格技术,耦合求解欧拉方程和六自由度(6-DOF)运动方程,对内埋式导弹在亚声速和超声速条件下的分离轨迹和姿态变化进行了比较,分析了不同马赫数下气动力对导弹分离特性的影响,重点讨论了超声速条件下导弹分离过程中相对于载机的位置和导弹的姿态变化情况,分析了超声速飞行环境下导弹安全分离的条件.     

共轭梯度方法在求解不可压N-S方程中的应用

采用非结构网格上的SIMPLE(Semi-implicit Method for Pressure-linked Equations)算法、k-ε湍流模型求解了三维不可压Navier-stokes方程,对低速不可压粘性流场进行数值模拟,在求解压强方程时采用共轭梯度方法,求解其它变量方程时,采用预处理的BICG(biconjugate gradients)算法.对NACA0012翼型绕流流场和飞艇绕流流场进行数值模拟并对结果进行分析,取得较好的结果.     

栅格翼导弹流场混合网格N-S方程数值计算

采用结构和非结构混合网格技术,对栅格翼导弹黏性流场进行数值模拟,预测其气动特性.计算的马赫数为0.7～2.5.用有限体积法和LU-SGS算法求解N-S方程.计算与实验数据以及栅格翼附近流场结构进行了比较.     

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

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

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

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

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

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