大迎角分离流动是现代飞行器气动特性研究中的常见问题,其重要的应用价值不言而喻。然而,迄今为止要对此类流动进行准确地预测是非常困难的,数值模拟的可信度往往难以令人满意。这一方面是由于物理外形过于复杂,难以生成高质量的网格,影响了数值解的精准度。另一方面是由于流动中包含了许多非定常、非线性的物理机制,给数值模拟方法带来严峻挑战。面对这一难题,本文从湍流模型、数值格式、计算网格三个关键因素出发,开展了相关探索和研究,期望能够建立一套更加工程可靠的数值模拟方法,进一步提高复杂飞行器大迎角气动特性数值模拟能力。由于非结构/混合网格对复杂外形具有良好的适应性,因此全文的研究以此类网格为基础。针对大迎角大范围分离湍流模拟问题,本文选用了当前最具工程应用前景的脱体涡模拟(DES)模型。鉴于传统二阶精度迎风格式具有较大的数值耗散和色散,不利于多尺度湍流结构的精细模拟,本文发展了一种基于原始Roe格式的自适应耗散混合格式,较好地控制了数值耗散对DES类模拟的影响,同时还建立了网格自适应方法,可以高效地增加DES类模拟流场的空间分辨率。最后基于这些改进措施进行了三角翼和现代战斗机构型的大迎角分离流动IDDES模拟,并初步开展了飞行器大振幅俯仰运动的IDDES模拟,计算与试验两者数据较为吻合,表明本文算法具有良好的应用前景。全文分为七章。第一章为引言,简要介绍了飞行器大迎角气动特性研究的背景和意义;综述了国内外在该领域的研究进展,包括机理研究、试验技术、数值离散方法、湍流模拟方法、网格技术等方面的发展概况;最后简述了本文的主要工作。第二章介绍了本文采用的数值模拟方法。为了提高二阶精度有限体积算法的湍流数值模拟能力,在原始Roe格式基础上建立了与DES类方法相匹配的自适应耗散混合格式;通过Fourier分析讨论了不同格式的色散、耗散特性,研究表明在二阶精度下混合格式的数值特性有利于DES类模拟;分别基于Spalart-Allmaras(SA)一方程湍流模型和k-ω剪切应力输运(SST)两方程湍流模型,建立了基于非结构/混合网格的DES类方法。第三章对数值模拟方法进行了算例验证。首先利用均匀各向同性衰减湍流算例对整套DES算法的经验参数进行了标定;然后通过雷诺数(Re)为3900的圆柱绕流、NACA0021翼型60°大迎角分离流动、超声速圆柱底部流动和NACA0012翼型俯仰振荡等典型算例,对本文DES类方法进行了检验。计算表明,DES类方法的数值表现总体优于RANS方法,而采用混合格式的DES类算法能够有效减小数值耗散对DES类计算的影响,解析出更小尺度的湍流旋涡结构,计算数据更接近试验值。第四章描述了本文建立的混合网格自适应方法,其目的是期望将DES类方法、混合格式、网格分布三者有机结合在一起,在进一步提高数值解分辨率的同时又尽可能地减少工作量。自适应策略采用了h–type方式;对于自适应过程中产生的悬空点,通过向相关网格单元插入对应节点,来消除其影响,使网格单元之间完全相容,自适应生成的网格可直接用于计算;最后,由一些简单算例验证了程序的正确性。第五章进行了65°后掠三角翼旋涡流动的数值模拟,对自适应耗散混合格式、网格自适应技术和IDDES方法开展了综合应用。计算发现:采用原始Roe格式的IDDES计算仅模拟出了泡型涡破裂,而采用自适应耗散混合格式的IDDES计算模拟出了泡型、螺旋型和双螺旋型等涡破裂形态;基于自适应网格的IDDES计算则捕捉到更为精细的流场结构。第六章开展了复杂飞行器大迎角分离流数值模拟的应用,选取的外形为类F22战斗机模型。首先简要介绍了风洞试验情况。随后进行了静态大迎角分离流的计算。从预测的气动力系数与试验数据的对比来看,IDDES模型相对于RANS模型仍然具有明显的优势。由流动结构的详细比较表明,IDDES模拟流场具有更强的三维非定常特性。最后,初步数值模拟了类F22战斗机模型大振幅俯仰运动,考虑到计算网格的尺度较粗以及动态试验的误差修正等不确定因素,IDDES计算结果与试验测量数据之间的相符程度还是不错的。第七章结束语,总结了本文的工作,指出了其中的不足,并对今后的研究方向进行了展望。