CT-1标模大迎角静态气动特性数值模拟

 采用雷诺平均的Navier-Stokes方程，对接网格技术和TRIP（TRI sonic Platform）软件，数值模拟了CT-1标模大迎角静态气动特性。来流马赫数0.5，迎角0°～100°，侧滑角0°和5°，在与1.2 m跨声速风洞试验结果对比的基础上，重点讨论了对流项离散方法、湍流模型、远场距离和不同尾部处理形式对CT-1标模大迎角静态气动特性影响。研究结果表明，不同尾部处理形式对大迎角静态气动特性数值模拟结果影响显著，模拟试验的尾部支杆是对比大迎角静态测力试验结果与计算结果的基本要求。     

三角翼大迎角风洞试验支架干扰数值模拟研究

现代战争要求战斗机能够在大迎角(AOA)状态下进行过失速飞行,对飞机大迎角绕流流场的研究主要的方法有风洞试验和数值模拟。在大迎角风洞试验中,常用的是尾支撑方法,支架的存在会对模型的试验结果产生一定的影响,本文通过数值模拟来对这个影响进行研究。以开源计算流体力学软件OpenFOAM 2.3为平台,采用PIMPLE算法求解Navier-Stokes(N-S)方程,PIMPLE算法是SIMPLE(Semi-Implicit Method for Pressure-linked Equations)算法和PISO(Pressure Implicit with Splitting of Operator)算法的结合体;采用基于有限体积的空间离散方法和空间二阶精度的线性插值方法,时间离散采用后向差分方法,湍流模型采用SA-DDES(Spalart-Allmaras-Delayed Detached Eddy Simulation)模型。为了验证方法的可靠性,首先对0°、10°、30°、50°、70°以及90°迎角下的有支架三角翼绕流流场进行计算,并将计算结果与试验结果进行对比,两者吻合较好。在此基础上,数值模拟了无支架的三角翼绕流流场,对比有/无支架情况下数值模拟结果,得到支架对三角翼绕流流场、背风面压强分布和气动力的影响。计算结果表明:大迎角情况下,有支架与无支架时相比,支架的存在会影响三角翼附近的流场(但是不会改变涡系等流动结构)、改变翼表面压强分布,从而导致三角翼的法向力系数和俯仰力矩系数发生明显变化。     

使用GAO-YONG湍流方程组计算翼型分离流

采用SIMPLE方法求解GAO-YONG不可压湍流方程组,对不同来流迎角下的NACA0012翼型绕流结构进行了数值模拟,给出了翼型绕流分离流结构随迎角的变化特征和翼型在分离绕流中的气动力参数.与实验数据以及大涡模拟结果的比较表明,GAO-YONG不可压湍流方程组能够对翼型绕流的分离点、分离涡形态、表面压力分布、升阻特性做出较好的预测,能够模拟翼型大攻角分离流动,计算结果优于FLUENT软件中的k-ε、k-ω、k-ωsst模型的计算结果.     

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

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

Gamma-Theta转捩模型在绕翼型流动问题中的应用

在非结构RANS求解器上实现了完全基于流场当地变量和转捩经验关系式的γ-Reθ转捩模型,并利用该模型对S809和NLF0416两种翼型在中等迎角以及大迎角下的低湍流度绕流问题进行了计算.计算结果表明该模型可以准确地预测自然转捩、分离泡转捩和再层流化等各种转捩现象,中等迎角下预测的转捩发生位置和升阻力系数与实验吻合较好,大迎角下该模型也捕捉到了前缘分离泡.同时计算结果也表明该模型预测结果受物面流向网格密度和对流项数值离散格式影响较大,计算时须加以注意.     

BGK-NS格式在复杂流动中的应用研究

作为构造能够同时模拟稀薄和连续流动的混合算法的第一步,本文在前期工作的基础上,考察了气体动理学BGK-NS格式在复杂流动数值模拟中的能力,选取了高超声速再入飞行器大迎角侧向喷流及带栅格翼导弹超声速绕流两个算例.数值模拟结果与实验,以及文献中的计算结果都吻合很好.这充分说明了BGK-NS格式对复杂流动具有良好的模拟能力,也为下一步混合算法的建立打下了良好的基础.     

亚声速大迎角模型试验洞壁干扰修正方法研究

通过测量洞壁附近的压力分布来模拟透气壁试验段的洞壁边界条件，采用数值求解Ｅｕｌｅｒ方程的方法模拟模型在风洞中的绕流场，然后将洞壁边界条件用远场边界条件替换，用同样方法算出模型在自由流中的绕流场，从而计算出洞壁干扰对模型的绕流场和气动力的影响。针对大迎角模型的洞壁干扰问题的特殊性，采用重叠网格技术，对飞机大迎角标模和ＳＢ－０３模型进行了亚声速大迎角的洞壁干扰计算与修正。     

CRM-WB风洞模型高阶精度数值模拟

通过求解雷诺平均Navier-Stokes(RANS)方程,采用5阶空间离散精度的WCNS和多块对接结构网格技术,开展了CRM翼身组合体风洞试验模型的高阶精度数值模拟,计算构型和计算状态来自AIAA第6届阻力预测研讨会(DPWⅥ)。主要目的是采用高阶精度方法,评估静气动弹性变形和模型支撑装置对CRM翼身组合体数值模拟结果的影响。通过与刚性外形的计算结果和NASA NTF风洞试验结果的对比,高阶精度数值模拟结果表明:迎角为3.0°时,静气动弹性变形使得机翼上表面激波位置前移并显著降低了外侧机翼上表面激波位置前的负压;迎角为3.0°时,模型支撑装置使得机翼上表面激波位置进一步前移,并导致翼梢处机翼上表面流动结构发生变化;迎角为4.0°时,计算模型中没有包含模型支撑装置是导致升力系数下降的主要原因;计算模型中包含机翼静气动弹性变形和模型支撑装置的数值模拟结果更加接近试验结果。     

平流层飞艇大迎角气动特性分析

为了获得大迎角下的平流层飞艇的气动特性,基于N-S方程及k-ωSST湍流模型,对某平流层飞艇的大迎角气动特性进行数值模拟分析.结果表明：平流层飞艇的大迎角失速特性及力矩特性与飞机具有较大的区别,飞艇的失速迎角较大（45°左右）;俯仰力矩特性在失速前呈现完全不同的两种趋势,小迎角范围内其静不稳定,迎角增加到一定程度后变为静稳定且维持该状态直至失速.研究结果对平流层飞艇的大迎角操稳设计、飞控设计具有重要的参考意义.     

飞机机翼表面霜冰的三维数值模拟

基于欧拉两相流理论对三维情况下飞机机翼表面的霜冰进行了数值模拟.根据水滴拟流体模型建立三维水滴控制方程;提出一套水滴控制方程的数值求解方法;由三维水滴流场的求解结果计算机翼表面的水滴收集特性,提出一种三维积冰外形的生成方法,完成了对飞机机翼表面霜冰的三维数值模拟.对ONERA M6机翼在不同迎角下霜冰的积冰情况进行了数值预测,并分析了结冰条件对积冰的影响.     

一种前体加宽型高超声速进气道试验方案研究

 根据矩形截面高超声速进气道前体的流动特征,对一种前体加宽型高超声速进气道试验方案开展了数值仿真及高焓风洞试验研究。首先,对不同前体宽度的高超声速进气道开展了三维数值仿真研究,结果显示:随着前体宽度的增加,进气道的流量系数和静压比逐渐增加,而总压恢复系数和隔离段出口马赫数逐渐减小,表现为先急后缓,且当来流马赫数和来流攻角变化时依旧保持上述变化规律。其次,对前体加宽型高超声速进气道试验方案开展了高焓风洞试验研究,结果表明:加宽前体可有效地提高进气道的流量系数,较为真实地反映此类进气道的流动特征,试验结果与数值仿真结果吻合较好。考虑到进气道性能参数随前体宽度变化规律表现为先急后缓,建议在试验条件下前体宽度比取0.5~0.8之间较为适宜。     

有间隙涡轮叶栅壁面流谱实验与数值模拟研究

对具有 3 .6 %相对叶顶间隙的涡轮叶栅进行了壁面流动显示和三维流场数值模拟 ,分析了大叶顶间隙涡轮叶栅的壁面流动特点 ,从数值模拟上特别对尾缘附近的壁面流动结构进行了详细讨论。实验与计算结果表明叶栅尾缘附近的流动是非常复杂的 ,同时由于叶顶间隙的存在 ,上下尾缘附近的壁面流谱明显不同     

细长体大攻角绕流的数值模拟

 用数值模拟方法研究了细长体大攻角的绕流特性（求解在周向３６０°范围内进行）。数值模拟的出发方程是三维可压流的全Ｎ－Ｓ方程,格式为二阶迎风ＴＶＤ格式。通过数值模拟研究了附近,攻角１９°、３２．５°和５０°状态下,细长体绕流场中的激波、分离流动及表面压力分布等重要特征。     

船用燃气轮机间冷器流路的数值计算分析

借助计算流体动力学(CFD)数值模拟技术,实现了船用燃气轮机间冷器整个流路的总压损失分析.数值模拟计算中针对间冷器模块采用了等效流通面积和等效换热效果的方法,有效解决了整个流道计算间冷器建模问题.数值模拟计算结果表明,CFD数值模拟技术可以有效的分析包括间冷器在内的整个流路的流场细节,改善间冷器的流路设计.CFD数值模拟技术与间冷换热器设计技术的结合,可以有效提高船用燃气轮机间冷器整个流路的优化设计水平.     

不同壁温及差分格式下超燃冲压发动机的仿真

针对圆截面超燃冲压发动机相关实验进行数值仿真，采用了不同的壁面温度条件和差分格式，并探究其对仿真结果的影响规律。结果表明采用常用的绝热壁面边界条件会导致激波边界层交互位置处的回流区变大，进而导致流道收缩，壁面压力升高。采用高精度格式不仅不能更为精确地计算回流区的大小，反而会使回流区面积变得更大。通过改变壁面边界条件，采用300 K的等温壁面，可以使激波边界层交互位置处的回流区变得极小，使仿真结果与实验结果十分吻合。计算结果表明，针对当前构型与实验条件，盲目提高数值格式的精度并不一定会得到更好的结果，相反可能会使仿真与实验更加偏离，而适当修改壁面温度条件，即使数值格式精度较低，仍然可以得到很好的仿真结果。最后，针对与实验结果吻合的仿真结果，分析了圆截面超燃冲压发动机的流场特性，重点研究了凹腔剪切层及其质量交换特性。     

绕钝体分离流动的非定常数值模拟研究

利用非定常RNGk-ε模型对三种典型的钝体绕流问题进行了研究,这三种流动分别为:方柱绕流、圆柱绕流和绕立方柱流动。研究表明了采用非定常雷诺平均方法可以研究钝体绕流的非定常流动,但所获得的结果各不相同。对于分离角固定的方柱绕流和绕立方柱流动,数值模拟结果和实验以及大涡模拟结果吻合较好,对于分离角不固定的圆柱绕流,随着雷诺数的增加,数值模拟结果和实验结果相差越大,分析认为这是由于壁面模型无法准确预测流动中不固定的分离角所致。改进了绕立方柱流动的模拟方法,采用非定常模拟方法后获得的结果要比定常模拟方法得到的结果有了根本改善。     

超声速和高超声速横向喷流数值模拟

为了解横向喷流的干扰影响和喷口附近的流场结构,采用分块对接网格和“0”型网格技术,数值求解N—S方程来模拟超声速和高超声速流场中横向喷流的干扰流场。对两种尖拱弹身外形的超声速和高超声速喷流干扰流场进行了数值计算,计算结果与风洞实验数据吻合一致。在此基础上,开展了某型导弹多喷构型干扰流场的数值模拟,得出一些横向喷流数值模拟的结论。     

绕翼型分离流结构的数值研究

 利用大涡模拟技术，对不同来流迎角下的NACA0012翼型绕流结构进行了数值模拟，详细地给出了翼型绕流分离流结构随迎角的变化特征和翼型在分离绕流中的气动力参数。在数值模拟中，采用了弱压缩流的控制方程，用贴体坐标技术进行了网格生成。     

高超声速飞行器高温流场数值模拟面临的问题

随着高超声速飞行器目标光辐射和电磁散射特性研究的发展和深入,高温流场特性日益引起人们的关注。由于高温流场特性研究中涉及到非常多的复杂气动现象,如气动加热、烧蚀、辐射、燃烧、化学反应以及湍流等,因此其数值模拟面临着诸多挑战。这里基于连续流计算流体力学(CFD)技术和稀薄气体蒙特卡罗直接仿真(DSMC)方法,从化学物理模型建模、方法稳定性与数值求解效率出发,分析了高超声速飞行器外部绕流、尾迹和发动机喷焰三方面的流场特性数值模拟在不同弹道、热防护手段和飞行流域环境下所面临的问题。在此基础上提出了数值求解技术和化学物理模型建模今后需要发展的方向,为有效提高高超声速高温流场特性数值模拟效率、增加流场特性预测精度提供了指导,从而为研究流场对高超声速飞行器目标光辐射和电磁散射特性影响提供有效的基础数据。     

轻质复合材料发射箱盖结构设计与参数优化

 针对国内现阶段复合材料发射箱盖研制领域中所存在的一些缺陷,通过设计一种方形轻质复合材料发射箱盖结构提出解决方案。首先提出了结构中薄弱区与抛出部分的初步设计方案;然后对导弹发射时燃气流冲击箱盖过程进行了瞬态数值模拟,并根据数值模拟结果对箱盖结构参数进行了优化;最后针对优化结构冲破过程进行了数值模拟,结果表明,按优化后方案设计的箱盖结构既能实现结构的轻质化,又能在开盖过程中保证较好的冲破性能。