收敛-扩张喷管中运用激波矢量控制法的计算研究

采用二阶TVD格式的有限体积法耦合RNG k-ε湍流模型和非平衡壁面函数求解二维守恒型雷诺平均N-S方程，数值模拟了运用激波矢量控制法的二维收敛一扩张喷管中的流动现象，很好的模拟出了由于次流喷射所产生的激波及上下壁面的压力分布，并将计算结果与国外的实验结果进行了对比。     

基于动态嵌套网格的飞行器外挂物投放的数值模拟

运用Euler方程和六自由度动力学方程数值模拟了外挂物投放的气动特性与运行姿态和轨迹。为了模拟外挂物从飞行器分离后的绕流特性和相对运动，采用了动态嵌套网格技术，并发展了一种洞边界建立和外边点的插值单元搜寻蛇行法，提高了计算效率。通过算例，证明了本文方法的有效性和合理性。     

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

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

微射流作动器三维外流场数值模拟

采用有限体积法求解Faver平均N—s方程，对三维、粘性、非定常、不可压微作动器外流场进行数值模拟。计算结果揭示了微射流产生、发展与耗散的过程，和实验中观察到的现象比较吻合。旋涡对沿展向的不稳定性是涡环破裂、旋涡对的相干性消失的主要原因。二次涡的发现为以后研究微射流旋涡对的细微结构打下良好的基础。     

旋翼翼型-桨尖涡干扰的数值模拟

采用Kirchhoff方法计算NACA0012翼型的远场气动噪声,Kirchhoff面为包围整个翼型的一层网格。首先采用欧拉方程计算程序求解得到翼型的Kirchhoff面上在桨尖涡干扰下的近场气动流场解,再根据Kirchhoff方法求解远场声压。从而体现出在计算桨涡干扰噪声时,旋翼CFD技术与Kirchhoff方法相结合的良好效果。      

Euler方程的自适应多重非结构网格计算

介绍一种基于四面体网格的多重网格技术,用于求解三维Euler方程定常解,具有十分明显的加速收敛效果。求解Euler方程的基本算法采用有限体积中心差分格点格式。多重网格技术采用一组几何上相互独立的粗细不同的网格。流动变量、余量和修正量在相邻粗细网格间的相互传递通过线性插值实现。采用自适应技术对网格进行加密,可以获得较优的网格,并作为一套细网格,进行多重网格计算。数值实验表明本文自适应多重网格法十分有效。     

Chaboche热粘塑性损伤模型的应用研究

本文讨论了Chaboche粘塑性模型的材料参数确定问题,提出了一种简单实用的模型参数确定方法,并用本文提出的方法确定了TiAl材料的Chaboche模型参数,并用这些参数进行了不同应变率下的拉伸模拟计算、不同应力下的蠕变模拟计算和CT试样在阶梯载荷下的蠕变模拟计算。     

跨声速HTA涡轮导向器全三维流场计算分析

利用CFD软件CFX-TASCflow对具有试验数据的跨声速HTA涡轮导向器流场进行了计算，结果表明，计算的叶面等熵马赫数分布和出口气流角与试验值吻合良好；而计算的总压恢复系数与试验数据存在一定的差别。     

不同轴长比椭球体的绕流计算与分离形态分析

用拟压缩性方法求解不可压Ｎ－Ｓ方程,对４∶１旋成体椭球和其他轴长比椭球的层流绕流进行数值模拟。分别用二阶中心差分格式和三阶迎风紧致差分格式来逼近无粘项。从计算所得的表面摩擦应力线、空间流线和横截面流场特性,分析了１０°≤α≤７０°迎角范围内的分离和脱体涡系形态,算出了相应的气动特性。计算结果表明,三阶迎风紧致差分比二阶中心差分可更好地模拟分离涡结构。迎角直到４０°时,算出的流态仍然对称,而迎角为７０°时,算出的流态为非对称非定常,与实验流态基本相符。椭球体横截面形状越扁平,开式分离就越提前产生。椭球体横截面形状横向越窄,分离区就越小。由于改变椭球体横截面形状具有改变分离和脱体涡系的效果,本文的研究结果有助于改善三维物体的气动力特性,为改进机体设计提供理论依据     

考虑NND格式修正的非定常Taylor—Galerkin有限元法

以Ｔａｙｌｏｒ－Ｇａｌｅｒｋｉｎ法为基础，吸收了差分格式中ＮＮＤ格式的思想，并在时间方向上采用四阶Ｒｕｎｇｅ－Ｋｕｔｔａ方法，构成了适用于定常和非定常跨声速流动计算的有限元算法。利用文的方法，分别计算了跨声速叶栅流动和工程中一类复杂非定常轴对称跨声速流动。计算表明，方法克服了Ｔａｙｌｏｒ－Ｇａｌｅｒｋｉｎ法的局限性，对激波及流场的分辨率较好，可用于跨声速流动计算，并可用于非定常流场的计算。