含偏心孔复杂盘体三维有限元分析

一、复杂三维盘体网格划分的基本策略　　图 1为含偏心孔的压气机盘实体 ,略去其与盘强度分析无关的若干细节 (如榫槽和花边 )。图 1　偏心孔的压气机盘实体对于图 1所示的压气机盘 ,根据其周期对称性取 30°扇区为分析对象 ,含半个偏心孔。这是一个非常复杂的几何域 ,其分网时遇到的主要困难是形体不规则和三维各向尺寸差异过大。广义上说用自由网格生成技术是可行的 ,但是实际上是行不通的。主要原因有 :( 1 )为了使单元自身的几何性质保持良好 ,对四面体、五面体单元畸变值要大于 0 .4,六面体单元畸变值要大于 0 .7。这势必导致节点数剧…     

大型运输机动力增升喷流效应数值模拟

参照C-17运输机发动机安装位置,考虑内、外涵道分开排气,建立了外吹式襟翼动力增升全机几何分析模型以及相应的巡航构型.采用结构化多块网格技术,基于雷诺平均Navier-Stokes方法,分别对全机增升构型和单独发动机动力喷流进行数值模拟验证,在此基础上对外吹式襟翼动力喷流效应展开研究.对于低速动力增升构型,发动机喷流大部分直接冲刷襟翼下表面而后向下偏转,部分高速气流经襟翼缝道引射并加速后吹向襟翼上表面,两部分气流在襟翼后缘汇合并向下游延伸,喷流冲刷襟翼时存在明显展向横流特征.在动力喷流影响下,不仅襟翼环量大幅增加,缝翼和主翼上的环量也均有所增加,全机可用升力系数和最大升力系数均突破了机械式增升装置的极限,达到4.0以上.同时,全机低头力矩大幅增加,为纵向配平带来额外的压力.对于相应的高速巡航构型,发动机喷流主要影响机翼下表面的压力分布,使得全机升力减小,阻力明显增大.动力增升构型在基本翼设计过程中应充分考虑喷流的影响.      

稳定性导数计算方法研究

稳定性导数的准确计算对于飞行器的操稳特性具有重要意义.应用计算流体力学(CFD)方法中的非结构化动网格技术建立能够模拟飞行器做周期性俯仰运动的强迫振荡法,以国际动导数标模Finner导弹为验证算例,获得不同马赫数下俯仰力矩系数的迟滞环曲线,进而计算Finner导弹在不同马赫数下的静稳定性导数和动稳定性导数.结果表明:本文计算得到的俯仰静稳定性导数与试验数据非常接近;在亚音速和超音速范围内,动稳定性导数计算结果与试验值很接近,但在跨音速范围内,本文计算结果与试验曲线的规律性一致,但误差较大.     

非结构变形网格和离散几何守恒律

数值模拟流固耦合问题或多体分离问题的非定常流动时,常采用基于任意拉格朗日-欧拉（ALE）方程的有限体积法,涉及到变形网格和离散几何守恒律。在对变形网格算法进行综述时,按照构造思想分为物理比拟、椭圆光顺、插值、运动子网格（MSA）及其混合法共5类,分别介绍了基本原理、研究现状和适用范围,通过算例比较表明,径向基函数（RBFs）和运动子网格相结合的混合方法既有很好的变形能力,也有较高的计算效率,值得进一步发展和推广。在介绍了离散几何守恒律（DGCL）概念之后,采用二维几何模型进行分析,指出其机理是离散过程中体积增量与网格面元扫过的体积不相等造成的,把目前国内外应用的算法分为面积修正法、给定速度的面积修正法、速度修正法和体积修正法共4类,对其应用范围和存在的问题进行讨论,认为提出的体积修正算法既可以保证流固界面条件,也可以用于时间多层格式。     

一种三维并行自适应计算方法

基于串行网格划分软件METIS与并行化消息传递编程接口（ MPICH2）对现有串行自适应程序进行简单的并行化改造,给出了一种三维可压缩无粘流数值模拟的并行自适应方法。首先利用单个进程调用METIS,串行划分网格;然后对所有进程并行计算以获得初始网格下的流场解;再次利用单个进程对整个流场运用自适应方法进行局部网格加密并调用METIS串行划分网格;最后全部进程在流场初始解的基础上继续并行计算,以获得自适应网格下的流场解。数值模拟算例验证了此方法的可靠性与高效性。     

混合网格方法研究栅格翼亚跨超声速气动特性

采用基于结构/非结构混合网格的CFD方法对栅格翼的亚跨声速气动特性进行了研究,计算方法经过试验数据的验证,可达到工程精度。对单独栅格翼的研究表明,跨声速壅塞时,流动通过栅格前方亚声速气流减速从栅格外侧溢流而实现流量调节;超声速壅塞时,流动通过在栅格翼上产生网状脱体激波使来流减速至亚声速来进行流量调节,且计算结果与理论估算上下临界马赫数范围一致。对某栅格翼气动布局导弹的计算结果表明,在跨声速壅塞区间内,壅塞有自我调节作用,通过溢流使栅格通道内维持相似流动,从而使栅格气动性能保持基本平稳;同时由于溢流使流量下降,栅格翼升力效率下降,导弹静稳定性减弱。     

非结构网格限制器研究

研究了6种限制器形式及其参数计算方法,其中新的自适应VenkatakrishnanM限制器能够根据当地流场特征自动调节参数从而实现不同区域不同特性.数值试验表明:非结构网格限制器对邻居集合选择方法较敏感,其中共面法比共点法精度高约15%;高马赫数时压缩型限制器比耗散型限制器计算精度高约10%~40%;参变量算法中,方向梯度法得到的流场比增量法平滑;自适应VenkatakrishnanM限制器能够模拟复杂超声速流动.      

超声速内埋武器不同分离方式分析

为给高空高速无人机内埋武器分离方式的选择提供权衡依据,建立了弹舱及弹体模型,采用chimera嵌套网格方法与Menter SST k-ω湍流模式,对舱内重力投放、舱内弹射投放及舱外重力投放三种内埋武器分离方式在高空高速条件下的内埋武器分离过程进行了仿真分析,对比了三种分离方式分离过程的流场特性与弹体运动参数。结果表明,在Ma=3.7,高度为2.5 km的条件下,舱内弹射投放与舱内重力投放能够使内埋武器安全分离,舱外重力投放在高空高速条件下无法完成内埋武器的安全分离,弹体无法快速下落并伴有大幅振荡,威胁无人机的飞行安全。     

基于三维非结构混合网格的离散伴随优化方法

伴随优化方法在优化过程中不受设计变量个数的限制,有巨大的工程应用前景。基于非结构混合网格技术,采用雷诺平均Navier-Stokes方程、离散伴随优化方法以及自由形面变形（FFD）技术,建立了飞行器气动外形优化设计系统。离散伴随方程采用公式推导的方法直接获取,并运用LU-SGS迭代求解。通过与差分结果进行比较的方式,对通量的雅可比矩阵和目标函数敏感导数开展了确认,并对文中黏性通量项、限制器以及湍流黏性系数的近似处理方法做了分析对比,验证了近似处理方法的合理性和可行性。采用建立的优化系统,完成了ONERA M6机翼在跨声速条件下的减阻优化,并比较了有无容积约束下优化结果的差别。模拟验证结果表明,建立的飞行器气动外形优化设计系统具有较高的可靠性和有效性,可以用于三维飞行器外形的减阻优化。     

几何尺寸对高超声速进气道气动性能的影响

为了探索模型缩尺比对高超声速进气道气动性能的影响,对不同缩尺比的二元高超声速进气道开展了数值模拟研究,结果表明：随着缩尺比的增大,进气道流量系数、隔离段出口总压恢复系数和马赫数均逐渐增大,而静压比逐渐减小,且来流马赫数越高,上述参数变化幅度越大。由理论与数值模拟分析可知,上述现象主要是由于不同缩尺比下,进气道当地雷诺数不同,导致进气道附面层相对厚度变化,进而影响进气道气动性能。理论分析了进气道总压恢复系数与缩尺比的定量关系,就进气道而言,进气道进口处附面层相对厚度减小1%,隔离段出口总压恢复系数提高约0.7%。