高阶精度方法下的湍流生成项对跨声速流动数值模拟的影响研究

利用五阶空间离散精度的WCNS格式和多块结构网格技术,通过求解雷诺平均NS方程,开展了SST两方程模型不同湍流生成项组合方式对跨声速流动数值模拟影响的计算分析。研究的主要目的是为高阶精度格式在复杂外形上的工程应用提供技术支撑。计算模型采用了RAE2822超临界翼型和DLR-F6翼身组合体构型。研究内容主要包括不同湍流生成项对残差收敛历程、边界层湍流粘性系数分布、边界层速度分布、压力系数分布以及模型整体气动力特性的影响。不同湍流生成项组合方式的流场计算结果还与风洞试验数据进行了对比。研究结果表明:对于小迎角不存在明显分离的跨声速流动,不同湍流生成项对流场的高精度计算结果的影响很小,可以不用考虑。     

两方程湍流模型的可压缩性修正及其应用

本文在不可压两方程湍流模型的基础上进行可压缩性修正,给出了ｋ－ε、ｋ－ε两方程模型的可压缩性修正方程,并对亚声平板湍流流动和跨声喷管流动进行了数值模拟,得到了较好的结果。从而证实了两方程湍流模型经过可压缩性修正,能较好地模拟可压缩流动。     

基于Van Leer＋AUSM混合格式求解三维N-S方程

基于多块粘性结构网格,开展了三维N-S方程数值算法的研究。控制方程的空间离散采用有限体积法,在前人工作的基础上,发展了Van Leer＋AUSM混合格式并应用到对流通量的离散中,粘性项采用中心格式离散并利用格林定理计算粘性通量中的导数项,时间推进采用五步R-K法,湍流模型为S-A一方程模型。最后,以M6机翼和某超声速弹丸的粘性流场作为数值算例,计算表明：发展的数值算法对跨声速、超声速流场均具有较高的分辨率,适用于跨、超声速流场的数值模拟。     

高阶精度格式WCNS在三角翼大攻角模拟中的应用研究

采用5阶精度的加权紧致非线性格式(WCNS-E-5)数值模拟了65°三角翼的大攻角绕流流场,主要目的是考核高阶精度格式WCNS在大攻角旋涡流动方面以及跨声速流场的激波附面层干扰、涡破裂位置的模拟能力,重点研究不同网格规模和湍流模型对尖前缘三角翼涡系之间的相互作用的影响。通过求解任意坐标系下的雷诺平均N-S方程,采用5阶精度的加权紧致非线性格式(WCNS-E-5)和多块对接结构网格技术,两种湍流模型分别是一方程SA和两方程SST湍流模型,在与相应试验结果对比的基础上,详细研究了WCNS-E-5格式在跨声速大攻角旋涡流动中的表现,以及不同网格规模、两种湍流模型对主涡二次涡相互作用、涡破裂位置和表面压力分布的影响。本文的研究结果表明,高阶精度格式WCNS-E-5能成功应用于三角翼的跨声速大攻角流动,网格规模的增加进一步提高流场分辨率,SST湍流模型相对SA湍流模型在三角翼大攻角流动中具有更好的适用性。     

考虑转捩的跨声速气冷涡轮叶片气热耦合计算

为了研究转捩对气热耦合计算的影响,在B-L代数模型与SST(shear stress transport)k-ω二方程模型的基础上,增加了两类基于间歇因子的转捩模型:代数AGS(Abu-Gharmam Shaw)模型与一方程间歇因子输运方程.选取NASA-MARKⅡ跨声速叶片为算例,分别采用全湍流模型与加入转捩的模型进行气热耦合计算.数值计算结果与试验对比表明由于能够预测附面层中的转捩过程,采用转捩模型的耦合计算得到的结果与试验吻合最好,由于在叶片壁面附近的网格较粗,采用间歇因子输运方程的转捩模型计算的结果要逊于采用代数转捩模型的结果.      

采用TRIP2.0软件计算DLR-F6构型的阻力

采用"亚跨超CFD软件平台"(TRIP2.0)数值模拟了DLR-F6构型,主要目的是通过计算DLR-F6构型的安装阻力考察TRIP2.0软件的数值模拟精度,并为运输机构型的气动特性计算积累经验.本文数值模拟采用的多块对接网格,测压和测力的试验结果均来自AIAA CFD Drag Prediction Workshop II(DPWII),对比计算结果采用了CFL3D的结果.本文详细研究了网格密度、湍流模型对DLR-F6翼身组合体和翼/身/架/舱复杂组合体两种构型的的总体气动特性和压力分布的影响,计算结果与相应的试验结果取得了较好的一致.本文采用SST两方程模型计算两种构型均得到了网格收敛结果;不同的湍流模型对压差阻力影响较小,对摩擦阻力影响较大;不同的网格密度和湍流模型对压力分布影响较小.     

预旋进气的盘腔系统内流场的数值研究

用CFD软件(FLUENT6.0.12)对预旋进气的盘腔进行了数值模拟，计算中采用了多种湍流模型，包括Spalart—Allmaras单方程、K—epsilon和K-omega多种双方程模型，结果表明Spalart-Allmaras模型与结果符合得较好；通过计算研究了流场结构，流场主要由两个参数旋转比βp和湍流参数λT控制，λT主要决定源区的大小，计算结果验证了自由涡流的存在，并且盘腔压力分布的计算值与理论曲线符合得很好。     

显式代数雷诺应力模型在超声速/高超声速流场计算中的应用

为了提高显式代数雷诺应力模型在高超声速流场中的计算精度,以修正后的CG k-ε模型作为湍流时间尺度决定方程,将Wallin和Johansson显式代数雷诺应力模型(WJ2000模型)应用于超声速/高超声速流场数值计算,并与线性湍流模型的计算结果进行了对比分析。结果表明WJ2000模型对高超声速压缩拐角附近壁面压力和热流分布的预测要优于线性湍流模型,同时网格依赖性明显低于后者;对CG k-ε湍流时间尺度决定方程所做修正不影响WJ2000模型对超声速流场的预测结果,但使其对高超声速流场的预测精度明显提高。     

超声速横喷干扰湍流场数值模拟

本文探讨了超声速平板横向喷流干扰湍流流场的数值模拟方法。分别采用Roe、Harten-Yee及AUSM 格式对三维雷诺平均Navier-Stokes方程进行离散,湍流模拟应用Baldwin-Lomax代数模型,讨论了B-L模型中Fy最大值的不同选取方式。分析了数值格式、网格划分对计算结果的影响。将来流马赫数为5、喷流马赫数为3的实验结果与计算结果进行比较,选择了与实验符合较好的数值模拟方法,为进一步利用数值计算方法研究喷流干扰流场特性奠定基础。     

增程制导航弹绕流场数值模拟

研究某增程制导航弹绕流场的结构,为气动特性分析和外形改进提供依据。以三维N-S方程为主控方程,引入S-A一方程湍流模型,对某增程制导航弹的绕流场进行了数值模拟,给出了流场图像和气动特性计算结果,并将CFD计算结果与试验结果进行了比较分析。CFD计算结果与试验结果吻合良好,验证了所采用CFD方法的准确性和有效性,数值模拟结果已得到应用。