飞行器姿态角解算的全角度双欧法

为了实现飞行器超机动飞行情况下姿态角的连续解算,将传统双欧法进行了改进.通过引入判别函数使得数值模拟过程中正反欧拉角能一一对应,克服了欧拉方程的奇异性,消除了传统双欧法须限定俯仰为 的限制,从而将正反欧拉角取值范围扩展到了全角度,实现了飞行器的任意姿态角连续解算.采用改进的双欧法数值模拟了飞机斤斗飞行动作,时间积分采用了变步长龙格库塔方法,结果表明该方法简单有效,能得到连续的高精度解.      

先进上风格式跨声速CFD计算研究

介绍了空气动力学跨音速计算中应用广泛的几种上风格式,同时引入了JAMSON中心差分格式作为比较,针对激波分辨率,格式精度,计算效率等方面的问题对ONERAM6机翼进行了对比性计算,讨论了ROE的FDS格式,AUSM+格式和AUSMDV格式的计算性能,得出了很有意义的结论。     

一种改进的跨声速旋成体壁面脉动压力经验预测公式

使用大涡模拟（LES）方法对典型整流罩旋成体结构进行非定常数值模拟研究,结果表明引起壁面压力脉动的流场结构随时间和空间变化,不同的流场结构之间相互干扰。同时研究发现,旋成体折角膨胀会造成流动在较短的距离内迅速加减速,同样引起壁面压力脉动。基于上述流动现象,重新划分压力脉动产生区域,通过引入空间修正项,提出一套改进的跨声速旋成体壁面脉动压力经验公式。算例表明,该方法能够显著提高旋成体在跨声速条件下的脉动压力分布的预测精度。     

组合前缘乘波体预压缩性能计算研究

乘波体飞行器是近年兴起的新概念飞行器.现提出了一种组合前缘的乘波体设计方法以改善飞行器前体的预压缩性能,并用Navier-Stokes方程进行了数值模拟,就密度流、压力系数、马赫数、总压等参数与普通乘波体进行了对比,结果表明,该设计方法能显著改善飞行器前体的预压缩性能.     

非定常流动计算的混合时间推进方法研究

结合ADI方法的二阶时间精度特性和高斯-赛德尔迭代法的良好收敛性特性,提出了一种ADI-高斯-赛德尔混合迭代格式,应用该格式,对NACA0012翼型的跨声速俯仰振荡进行了EULER方程求解,空间离散格式采用ROE’S FDS格式,计算结果与实验及成熟的LUSGS双时间步迭代方法计算结果进行了分析比较。计算结果表明,该格式精度能达到LUSGS双时间步迭代方法完全收敛的二阶精度,且该格式精度对内迭代次数的依赖性较低。分析表明,本文提出的方法具有较高的精度和较好的收敛性。     

飞行器大迎角下俯仰静、动导数的数值计算

 采用一种适用于非定常流动计算的逐次超松弛（SOR）时间推进方法耦合求解非定常Navier-Stokes方程和强迫运动方程,数值模拟了NACA0012翼型、弹道外形(HBS)和有翼导弹等飞行器在不同初始迎角下做强迫俯仰振荡的黏性动态流场,给出动态气动力及力矩系数的时间历程;采用积分法获得飞行器的静导数和动导数,将计算结果与实验结果及半经验理论预测方法进行比较并分析飞行器在大迎角下动态特性的非线性特征;计算结果与实验值相吻合,表明本文采用的非定常数值模拟方法在大迎角下能够较准确地模拟飞行器的动态特性,具有工程应用价值。     

一种新型的动态RANS/LES混合方法

为解决RANS/LES混合方法过渡依赖网格疏密的问题,使用熵函数fs作为交界面函数构造一种新型的动态RANS/LES混合方法。同时采用脱体涡模拟方法 (DES97),对平板边界层流动、超声速凹腔-压缩拐角流动以及超声速底部流动进行数值模拟,并与实验数据进行对比。结果表明,无论对于附体流动区域还是分离区域,fs都能够动态地指示出边界层的范围。构造的动态RANS/LES混合方法能够显著缓解DES97由于网格因素导致的雷诺应力模化不足的现象,同时保证分离区域涡结构的求解精度。另外,熵函数fs是从能量耗散的角度进行构造,并不依赖于特定的湍流变量,因此更具有普适性。     

重叠网格技术中割补法的研究与改进

针对计算流体力学中重叠网格生成困难的问题,分析了重叠网格自动生成技术中割补法的主要特点和基本流程,指出其在洞映射单元属性的判断方法和洞边界外推过程中出现的"孤点"现象等方面的问题,提出了合理的改进方案,并将改进后的方法应用到自适应嵌入网格和DRAGON网格的生成中,最后编制出通用程序并给出了三维算例.实际应用表明,改进后的方法具有更稳定、自动化程度更高、适用范围更广的特点.