高超声速热流高精度数值模拟方法

通过高超声速下钝锥热流的计算,对影响热流计算精度的因素进行了综合研究.获得了壁面法向网格Re数对热流计算的影响规律.采用边界高阶插值,提高了热流计算精度.用Roe和AUSM(Advection Upstream Splitting Method)+格式搭配minmod和混合limiter两种限制器,探讨了空间离散方法及限制器的耗散性对热流结果的影响.应用五阶WENO(Weighted Essentially Non-Oscillatory)方法,对高阶格式计算热流的性能进行了细致研究,最终确立了网格依赖小,热流精度高的计算方法.研究认为:高阶格式不但可以放宽热流计算对网格Re数的限制,同时使得格式对本身计算方法的依赖性也变小.      

基于闭环控制的带翼导弹虚拟飞行数值模拟

通过耦合求解非定常Navier-Stokes方程、刚体六自由度（6DOF）运动方程及舵偏控制律,发展了一套虚拟飞行数值模拟平台,对基本带翼导弹外形进行了虚拟飞行数值模拟研究。流场求解采用基于结构网格的有限体积法,时间推进采用双时间步法,使用刚性重叠动网格技术模拟舵与弹身的相对运动,采用四阶Runger-Kutta方法求解运动方程组,舵偏控制律采用基于比例-积分-微分（PID）的反馈控制。模拟结果表明：发展的虚拟飞行数值模拟平台能够处理包含复杂运动边界的非定常运动问题,可以为非定常气动力建模提供数据、预示飞行器稳定性及检验飞行控制律,该模拟平台具备一定的工程应用价值。     

转捩模式与转捩准则预测高超声速边界层流动

对原始的k-ω-γ转捩模式和"层流+转捩准则"模型进行了改进,在2种方法中分别增加了横流模态时间尺度和横流转捩准则用于预测横流失稳诱导转捩。通过对网格预处理可并行计算获得边界层外缘信息以及边界层内横流速度。采用不同雷诺数条件下的0°攻角尖锥以及HIFiRE-5外形对2种方法预测高超声速边界层转捩的性能进行了比较分析。研究结果表明,2种方法均能正确反映高超声速边界层转捩起始位置和转捩区长度随雷诺数的变化趋势,但不能捕捉转捩区热流峰值;"层流+转捩准则"模型计算得到的传热系数在全湍流区较k-ω-γ转捩模式偏高。对于同时存在流向不稳定和横流不稳定的HIFiRE-5外形,改进的k-ω-γ转捩模式和改进的"层流+转捩准则"模型相比于原始的模型均能更加准确地预测中心线两侧横流失稳诱导形成的转捩;对于中心线附近因速度剖面拐点引起的边界层转捩,"层流+转捩准则"模型由于与边界层厚度相关,预测得到的转捩位置较试验结果靠前,k-ω-γ转捩模式与试验结果吻合很好。     

适用于混合网格的改进雅可比迭代法及其应用

LU-SGS因有较高的鲁棒性和小的内存需求而得到广泛应用,然而用于混合网格计算前需要进行网格排序和分组来实现算法并行;此外,LU-SGS格式收敛效率不高。针对这些缺点,本文提出了一种改进的适用于复杂混合网格的雅可比迭代方法,无需网格排序和分组就可实现算法的并行化,且有较快的收敛速度。该方法编程实现简单,易于采用OpenMP实现并行。算例研究表明,相比于LU-SGS格式,在各来流条件下,本文提出的方法收敛速度更快,鲁棒性好,并行和串行结果一致,且内存需求增加很少。      

高超声速热流计算湍流模型性能评估

采用计算流体力学方法,通过对高超声速来流下双椭球热流计算,综合分析了影响壁面热流预测精度的因素.针对工程上4种常用的湍流模型:BL(Baldwin-Lomax)模型、SA(Spalart-Allmaras)模型、k-ω模型及SST(Shear-Stress Transport)模型,评估了其在高超声速热流预测中的性能,获得了湍流模型下壁面法向网格雷诺数对热流计算的影响规律.此外,对高马赫数下不同后掠角钝舵模型进行热流计算,获得了后掠角对钝舵前缘热流峰值的影响规律.研究表明,两方程湍流模型相比一方程和零方程模型更适合高超声速气动热计算,其中SST模型对热流预测精度更高,总体性能更为优异;钝舵驻点线最大热流随后掠角增大呈非单调变化,在后掠角约为22°时热流达到峰值.      

喷口布局对导弹侧向喷流控制作用的数值模拟

为保持飞行器在稀薄大气中的机动性,通常采用喷流反控制作用(RCS),但在超声速来流中,这会导致飞行器表面出现复杂的喷流干扰流场,对飞行控制造成了巨大影响。为提高对超声速条件下的侧向喷流控制作用的规律性认识,应用数值模拟方法,研究了超声速条件下的无舵光滑弹体和带尾舵的弹-翼组合体上的声速侧向喷流控制问题。开展了关于喷口布局对侧向喷流控制效果影响规律的研究工作,并通过引入法向干扰力沿程增长系数从定量角度加以分析。计算结果表明:在有尾舵的情况下,喷口位置的后移和马赫数的增加能够显著增强侧向喷流控制效果;当喷口位置位于舵面之前时,喷流干扰力放大系数随迎角增大而增大,随来流静压增大而减小;当喷流位置后移至舵面之后时,规律相反;在某些喷口位置和来流条件下,弹-翼组合体的侧向控制效果与无舵光滑弹体相比并不具备优势。     

DDES延迟函数在超声速底部流动中的性能分析

DDES是广泛应用的一类RANS/LES混合方法,其通过引入延迟函数保证近壁区的RANS模化,对分离流动十分有效。目前DDES已发展了多种不同的延迟函数,但对各延迟函数的性能特点认识尚不够充分,尤其缺乏超声速流动中的相关研究。围绕DDES方法中不同延迟函数开展研究工作,选取超声速底部流动作为测试算例,通过与实验数据的系统对比分析,考察不同延迟函数在超声速分离流动中的分布规律、作用效果及模型求解能力。研究表明,不同延迟函数作用范围与求解能力存在差异,其中DDES-F1能够在起到保护作用的同时不损害模型的求解精度,对该流动较为有效,所得结果与实验数据吻合较好。     

叶片式涡流发生器对压缩拐角流动分离的控制

采用基于k-ω湍流模型的雷诺平均Navier-Stokes（RANS）方程方法,研究了叶片式涡流发生器（VG）对于马赫数Ma_∞=2.9时24°压缩拐角边界层分离的控制作用。计算结果表明：叶片式涡流发生器诱发的流向涡,是控制拐角处边界层分离的主要因素,流向涡强度越大控制效果越好。流向涡增大了主流与边界层内的动量输运,沿壁面法向速度型更加饱满,并使得压缩拐角处的二维分离转变为三维分离,改变了激波边界层干扰的结构,分离区长度减小了39.68%。相比于相向旋转,同向旋转叶片式涡流发生器改善了分离区内的压力分布,分离区总长度减小量相当,但分离点距转折点处的长度更短,且系统阻力增量更小。对于相向旋转叶片式涡流发生器,后缘高度增大,分离区总长度减小,系统阻力增量先减小后增大;相向旋转叶片间距越大,分离区总长度越小,系统阻力增量越大;同向旋转叶片间距越大,分离区总长度越大,系统阻力增量越小。高度对叶片式涡流发生器诱发的流向涡强度起主要作用,异向与同向叶片间距的影响较小。     

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

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

高超声速热化学非平衡空间格式的扩展与改进

扩展、改进和分析了ROE,MAPS,LDFSS,HLLC,HCUSP和ECUSP6种高分辨率空间格式于热化学非平衡流的计算中.采用5组分17反应的Dunn-Kang化学反应模型和Park的双温模型计算了三维球头算例,得到如下结论:ROE和MAPS格式精度高于LDFSS,HLLC和CUSP类格式;基于ROE平均的界面音速可提高格式的鲁棒性及斜激波捕捉能力,同时对流场参量的计算无影响;ROE平均音速的引用改进了MAPS和ECUSP格式.为热化学非平衡流场计算格式的选择提供了参考.