半正交曲线坐标系下三维不可压边界层的一种差分解法

本文提出一种三维不可压边界层的差分解法。此法中采用了半正交曲线坐标系和涡粘性各向异性的湍流模型,数值分析中采用了变换域和Keller分块差分格式。用本法计算了一个三维边界层实例,计算结果同实验数据吻合得很好。     

前缘剖面形状对80°／60°后掠双三角翼上涡流运动影响的数值模拟研究

采用上风高分辨率格式———通量差分分裂格式离散求解三维可压非定常薄层N S方程 ,数值模拟前缘剖面形状对 80°/60°后掠双三角翼上涡流运动的影响。计算模型包括尖前缘、菱形前缘、圆形前缘等三种不同前缘形状的机翼。计算结果表明 :三种不同剖面形状的前缘可以诱导产生不同的前缘分离 ,形成的各前缘剪切层的特点也不同。尖前缘机翼的边条翼涡和外翼涡合并点最靠前 ,其次为菱形前缘和圆形前缘的机翼。在大迎角情况下 ,三种机翼上的内翼涡发生破裂。圆形前缘机翼上的内翼涡破裂点比其它两种情况下破裂点的位置较靠前。涡的合并对二次涡的结构和特点也有显著的影响。     

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

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

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

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

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

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

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

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

扰流片快速张开诱发的非定常升力过冲

通过数值求解二维非定常可压N-S方程,研究了二维翼型上扰流片快速张开诱发的复杂非定常流动问题.利用动态重叠网格方法处理扰流片和主翼相对移动,分析了扰流片快速张开诱发的升力过冲现象及其形成机理.着重探讨了来流攻角和扰流片高度对过冲量的影响.     

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

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

后缘喷流对三角翼前缘涡影响的数值模拟

为了研究三角翼后缘对称喷流对前缘涡破裂位置和旋涡结构的影响,通过高分辨率的N-S方程数值模拟方法,对60°后掠角三角翼后缘有对称喷流及无喷流情况下的绕流进行了研究.结果表明,后缘喷流速度与自由来流速度之比影响前缘涡破裂的位置.与无喷流情况相比,喷流与来流速度比的不同造成了涡破裂的提前或推后.而与传统结论有所不同的是,并非所有的后缘喷流都能延迟涡的破裂.另外,后缘对称喷流对涡轴位置的影响很小.后缘对称喷流不改变三角翼前缘涡横截面流动拓扑结构的变化规律,但影响极限环的扩张速度.      

战术导弹横向喷流数值模拟

为了研究喷流对导弹不同升力面的干扰影响,通过数值求解N-S方程来模拟超声速外流场中横向喷流的干扰流场,采用分块对接网格和"O"型网格技术,精确模拟喷口截面及弹翼形状,生成高质量的贴体计算网格.通过几种升力面和弹身组合的超声速横向喷流干扰流场的数值模拟,计算的马赫数为2和4.5,计算结果显示具有升力面的喷流干扰比光弹身的喷流干扰大.研究和分析了喷口附近流场的涡系结构和波系结构,并将导弹喷流干扰力放大因子及喷流引起的压心变化的计算结果与目前国外公布的试验值进行比较分析,二者显示了良好的一致性.