关于三维Navier－Stokes方程的粘性项计算

借助于张量分析和张量计算，在贴体曲线坐标系下本文讨论了不同求解变量导致了粘性项个数上的重大差异和不同大小的计算量，并提出了便于粘性计算的最佳形式。文中借助于有限体积离散技术，通过引进两个对称辅助矩阵［Ａ］和［Ｂ］，使粘性项的计算量大大减少，这对完成三维粘性流的数值计算具有重要的指导意义。借助于上述方法，本文完成了某型真实进气道两种工况的三维粘性Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程计算（即Ｍ∞＝３．０，α＝０°和设计工况Ｍ∞＝２．６５，α＝０°），获得了满意的全场结果；对于Ｍ∞＝２．６５的设计工况，同实验数据作了比较，符合良好。由于本文的方法明显的减少了粘性项的计算量且节省了大量内存，以致于使三维流场的Ｎ－Ｓ求解能在普通微机上进行。     

虚拟飞行试验技术介绍

飞行器在大迎角下机动飞行，流动状态的非定常性是其重要特性。对飞行器的非定常气动特性研究预测，主要还是以风洞试验为主。以前风洞试验主要是动导数试验和大迎角非定常试验。由于这两种试验的局限性，AEDC最近开始研究一种新的风洞试验技术一虚拟飞行试验技术。本文主要是对这种新的风洞试验技术研究情况做了介绍。     

转轴位置对俯仰动态特性的影响

研究了转轴位置对机翼、机身和三个组合体上翼面的动态特性的影响，对于机翼和机身，当转轴位置后移时其动态迟滞效应更为强烈。当翼面处于转轴之前，上仰时会出现更为的动态迟滞效应。因此鸭翼达到失速状态较晚。而当翼面处于转轴之后时，如尾翼，将会更早地出现动失速，转轴位置对动态迟滞特性的影响可以认为主要是由于俯仰运动时的动态洗流影响所导致的实际有效迎角的变化。由此影响了动态分离涡及其发展。     

预估校正的高分辨率格式及其应用

本文从基于通量分裂的一阶迎风格式出发,构造了一类校正型具有高分辨瘁的二阶全变差递减(TVD)格式,并用此格式进行了一维激波管问题及二维问题的计算。结果表明格式的解为二阶精度,激波处无明显的振荡,激波过渡区为1～2个点。     

再入飞行器高超声速绕流的一步和二步有限元计算

本文针对再入飞行器的空气动力问题，采用一步和二步有限元格式求解可缩守恒型Ｅｕｌｅｒ方程，并结合不同的人工粘性模型来捕捉激波，给出了高超声速二维钝头绕流的结果。     

超声速主流中横向喷流干扰流场数值模拟

本文用二阶迎风 TVD格式求解层流 N- S方程 ,数值模拟了超声速主流与横向喷流相互作用复杂干扰流场 ,得到了合理的涡系结构 ,计算波系结构及沿物面压力分布也与实验结果一致。计算还表明 ,TVD格式不但具有较强的捕捉激波的能力 ,而且用来模拟复杂涡系结构也是有效的     

三维可压缩混合层中扰动演化的研究

为了利用被动控制增强可压缩混合层中混合效率,本文用数值模拟的方法研究了三维可压缩混合层中扰动向下游的演化.通过在入口处引入不同的不稳定T-S波,给出了扰动随空间的演化过程,及马赫数对扰动演化的影响.计算结果表明在马赫数较小时,扰动波增长很快,当马赫数增大时,扰动所形成的流向涡的强度会减弱.这一结论与流动稳定性分析结果是一致的.     

突扩后突缩的三维流场计算

根据有关固体火箭发动机装药后端到喷管间的空间特征,对三维突扩后突缩内部流场计算进行了探索和研究。计算中成功地应用了SIMPLE算法解决复杂的三维流场计算问题,并在初场构成、源项离散及松弛系数的选取等方面都进行了调试改进。通过建立计算模型、划分网格、对方程及边界条件离散化和非线性处理,以及线性方程组求解等,计算得到了符合物理规律的数值解,指出了空间涡及二次涡的存在。     

大攻角非对称流动的非定常弱扰动控制

 研制了一种新的大攻角细长旋成体非对称涡的主动控制技术, 即在细长旋成体头部施加非定常弱扰动来控制头部非对称背涡。应用七孔探针测量的空间截面流场揭示了非定常控制下非对称涡变成对称涡的流态特征。测力试验研究结果表明该方法不仅能完全消除背涡的非对称性及其产生的侧向力, 并且有效控制攻角范围从30b直到80b。     

跨音速压气机动静叶排相干的三维非定常流动数值分析

采用时空二阶精度ＮＮＤ显式差分格式求解三维Ｅｕｌｅｒ方程，对跨音速压气机级中动静叶排相干形成的非定常流动进行了数值分析。除动静叶排流动的控制方程体系、数值方法外，着重讨论了动静叶排计算域、网格以及边界条件的处理。对某压气机第１级动静叶排相干的非定常流场的计算，给出了叶片表面的马赫数分布、压力波动幅值，以及叶片槽道中瞬态密度等值线等，并分析了结果的合理性。