新一代战斗机非定常流动数值研究综述

新一代战斗机强调超机动能力和强隐身性，其中大攻角下的静态失速、动态失速及内埋弹仓绕流是与高机动和强隐身密切相关的、极具挑战的几类典型的非定常流动，它们对数值方法提出了极高的要求。为了高精度地仿真流场、清楚地揭示流动机理，有效地控制非定常流动，非常有必要发展高精度且高效率的RANS-LES混合方法体系，包含RANS-LES混合方法本身、与RANS-LES混合方法匹配的高精度自适应耗散格式、基准湍流模式、高质量计算网格、高精度时间推进方法、非定常量的统计方法等，具有极强的紧迫性。提出、发展、验证并应用该类方法数值仿真新一代战斗机（包括单独部件、组合体、甚至全机）的非定常流动，数值预测结果与风洞实验数据吻合良好；此类方法可为新型战斗机设计提供理论依据和分析手段。     

正反顺序加载下安装边螺栓预紧力分散度的数值模拟与试验验证

在航空发动机机匣安装边装配过程中,不同位置螺栓之间会产生弹性交互作用,不同加载方案会使装配后的螺栓预紧力存在一定的分散度。针对这一问题,建立了3层带止口机匣安装边3维接触有限元模型,研究不同加载方案下螺栓预紧力的施加方法。基于有限元模型数值模拟了3层机匣安装边分别在顺序加载方案、交叉加载方案和正反顺序2轮加载方案下螺栓预紧力的分散度。搭建了相应的试验装置,分别对3种加载方案的预紧力分散度进行试验验证。结果表明：采用正反顺序2轮加载后的螺栓预紧力分散度最小,在3层止口机匣安装边上采用正反顺序2轮加载方案后的螺栓预紧力分散度较小,在有限元仿真中螺栓预紧力分散度为2.43%,在试验中螺栓预紧力分散度为4.98%。采用此方案加载后,止口机匣安装边能达到更好的密封效果和较好的刚度对称性。     

高速旋转冲压增程弹用进气道复杂流场数值模拟

采用块结构网格与二阶精度流场分区求解技术,对高速旋转、含侧向支柱冲压增程弹丸进气道内外复杂流场进行了数值模拟,得到了高速旋转工况下对应于不同来流马赫数和攻角,以及临界工况时超音速进气道内外流场复杂的波系结构。在旋转工况下,进气道流场结构与不旋转时相似,但旋转速度以及攻角的存在对冲压发动机进气道的总体性能产生了负面影响,进气道总压恢复系数和动能系数均有所降低,而流场畸变指数则显著增大;冲压发动机进气道在较低马赫数工作时的综合性能优于在较高马赫数工作时的综合性能。     

二维复值Ginzburg-Landau方程的一个高阶紧致ADI差分格式

对二维复值金兹堡朗道(Ginzburg-Landau,GL)方程提出一个基于时间分裂的高阶紧致交替方向隐式有限差分格式。本文通过时间分裂法将GL方程分裂成一个非线性子问题及两个线性子问题,对非线性子问题以及其中一个线性子问题均通过精确积分进行计算,并对另一线性子问题构造紧致交替方向隐式差分格式进行数值计算。实际计算中,在每一时间步,利用追赶法求解一族常系数三对角线性代数方程组,从而使得算法既具有较高精度又拥有较快的计算速度。数值实验表明该算法在时间和空间方向分别具有二阶和四阶精度,并模拟了方程的一些动力学行为。     

飞行器最优机动策略研究方法及进展

回顾了国内外最优机动策略研究的理论和方法,着重论述了近几年国内外的最优轨道机动策略问题的研究状况,比较了国内外的发展差别,讨论了轨道最优机动策略的发展趋势,为我国飞行器最优轨道机动策略研究方向提供了有意义的参考。     

超声速空气介质自由扩散混合流涡的演化

在不同谐波扰动促发方式下,观察二维超声速空气/空气自由扩散混合流中涡的形成及演化.考虑含空气组分扩散效应的可压缩Navier-Stokes方程,对流项采用3阶迎风紧致格式离散,输运项采用6阶对称紧致格式离散,非定常时间推进采用3阶紧致存储显式Runge-Kutta方法.在无相差情形下,获得了大尺度基频涡结构的饱和、一次对并、二次对并、三涡对并等现象.在Ma_c=0.3低对流Mach数下,基频涡较饱满,但流向尺度较小.受空气真实气体特性影响,在Ma_c=0.8高对流Mach数下未发现涡/小激波结构.      

适于黏性计算的自适应笛卡儿网格生成及其应用

复杂外形的黏性网格生成已经成为计算流体力学(CFD)在工程应用中的主要瓶颈,高效而鲁棒的网格生成技术研究显得尤为迫切。使用基于物面外形自适应加密的笛卡儿网格方法生成计算网格,投影方法拟合壁面边界,将投影得到的柱形单元沿法向分层获得适于黏性计算的网格。应用特征恢复技术恢复凹面特征并改善凸面处的网格质量。流场解算中应用最小二乘法重构获得2阶精度,Roe格式计算无黏通量,Holmes-Connell方法计算黏性通量。加入Spalart-Allmaras一方程湍流模型以获得湍流解。上下对称高斯 赛德尔迭代(LU-SGS)格式隐式时间推进,对网格重排序保证格式的稳定性和平衡性。算例显示网格生成方法是快速的和鲁棒的,流场计算精度满足工程需求,因此该方法非常适合应用于实际工程计算。     

旋翼非定常气动特性CFD模拟的通用运动嵌套网格方法

针对直升机旋翼非定常气动特性CFD模拟中的网格生成难题,提出了一套高效、通用的运动嵌套网格生成方法.首先,基于Poisson方程求解和翻折法生成旋翼桨叶的正交贴体网格.其次,针对旋翼桨叶的扭转分布及变距、挥舞等复杂运动,建立了一套通用的洞单元识别的扰动衍射法;为保证洞包络面的封闭性,完善了挖洞过程中网格加密策略;在洞边界确立基础上,提出了一种高效、鲁棒的最小距离法贡献单元搜索的改进方法.在此基础上,建立了基于RANS(Reynolds-averaged Navier-Stokes)方程的旋翼非定常流场CFD模拟方法.最后,采用所建立的方法分别对悬停和前飞状态下的C-T(Caradonna-Tung)旋翼和7A(Helishape 7A)旋翼的气动特性、桨尖涡的位置进行了计算,计算结果与试验值误差小于5%,验证了该运动嵌套网格生成方法在旋翼非定常气动特性CFD模拟中的有效性.      

两种空气/煤油燃气发生器富油燃烧组织方案对比

根据空气/煤油富油燃烧的特点,提出了两种空气/煤油燃气发生器富油燃烧组织方案,设计了采用钝体稳定火焰和旋流空气、二次喷注空气稳定火焰的两种燃气发生器.为了对比两种方案的点火和燃烧特性,对两种燃气发生器进行了一系列热试,结果表明余氧系数是燃气发生器最重要的工况参数.随着余氧系数的增加,燃气发生器的状态逐渐从启动失败变为中途熄火,最终呈正常启动状态.采用钝体稳定火焰的燃气发生器稳定工作的余氧系数边界为0.518,采用旋流空气和二次喷注空气稳定火焰可将该边界延伸到0.237,极大地扩大了燃气发生器的工作范围.与钝体稳定火焰的燃气发生器相比,旋流空气和二次喷注空气稳定火焰的燃气发生器的富油燃烧的燃烧效率提高了20%.两种方案结构复杂性相当,旋流空气和二次喷注空气稳定火焰的燃气发生器不需要冷却火焰稳定器,可提高燃气发生器的工作时间.      

含激波的湍流流动高精度大涡数值模拟方法

针对采用亚格子模型进行含激波的湍流流动模拟时会面临激波附近的精度损失问题,考虑从通过亚格子模型以及数值模拟方法两方面的改进来实现湍流流动大涡模拟的精度提高.大涡模拟采用了Yee及Sj(o)green (2009)提出的高阶低耗散方法.该方法采用自适应的流场探测器以控制计算中所需区域的数值耗散,并考虑对动力学模型采用在激波位置使用Sagaut和Germano(2005)提出的单边亚格子过滤器和(或)直接禁用亚格子项等方法加以改进.对于标准的马赫数1.5和3条件下的激波-湍流干扰问题,上述新方法相较于全区域采用亚格子模型的方法均表现出了相似的精度提升.同时实现的数值精度改进方案采用了Harten的亚单元分辨过程来定位和锐化激波,并在精确激波位置附近的网格点处采用了单边测试滤波.