高超声速飞行器钝舵缝隙流动数值模拟研究

目前,带缝隙钝舵的缝隙引起的流场结构和气动加热规律变化,还很不明确,需要研究缝隙诱导所形成旋涡的空间分布特征和旋涡运动对物面气动加热的影响规律。通过分析高超声速钝舵缝隙气动加热问题,基于无缝隙钝舵,建立一种带缝隙钝舵简化模型。使用有限体积方法求解可压缩Navier-Stokes方程,通量采用van Leer通量向量分裂方法计算。插值采用MUSCL方法,时间项采用LU-SGS隐式方法。结果表明:无缝隙钝舵流场结构相对简单,带缝隙钝舵流场结构同无缝隙钝舵相比要更为复杂,舵轴上游缝隙内会出现马蹄形涡串结构,相应地在缝隙的上下表面均会出现马蹄形高热流区;受缝隙诱导分离再附流动的影响,在舵轴迎风面以及舵体侧面后部均形成了局部高热流区。     

NNW-HyFLOW高超声速流动模拟软件框架设计

NNW-HyFLOW在国家数值风洞（NNW）工程支持下,由风雷开源软件提供基础框架支撑,将建设成为基于结构/非结构混合网格,面向高超声速应用领域的国产自主工业CFD软件,具备高温气体热化学非平衡效应模拟及其相关气动力、气动热和气动物理特性计算分析等主要特色功能。本文从框架设计、数据结构、耦合方法、并行计算方法以及接口设计等方面对软件设计思想和框架特点进行了介绍,给出了求解器采用的理论模型、核心数值算法及其实现方法,结合HEG风洞试验、RAM-C飞行试验、Electre飞行试验以及航天飞机OV102飞行试验等典型算例开展了数值模拟。研究表明：NNW-HyFLOW具有底层代码复用、功能兼容性好、扩展能力强和接口灵活等优点,其当前测试版本已经具备了较好的高超声速非平衡流动数值模拟能力,在热化学非平衡效应及其影响的气动力特性、气动热环境和等离子体分布特性预测与评估方面,具有较高的数值计算精准度,初步满足了高超声速复杂飞行器高温非平衡流动数值模拟的需求。     

高超声速激波/边界层干扰流动数值模拟研究

在基于Roe格式的全Navier-Stokes方程计算流体力学(CFD)代码中,发展了一种局部熵修正方法,克服了传统熵修正方法在边界层流动模拟中耗散过大的缺点,可用于更加准确的模拟激波/边界层干扰的复杂高超声速流动。对典型高超声速双锥边界层分离与激波干扰的复杂流动进行了模拟,研究了网格密度和熵修正方法耗散性对计算结果的影响。研究表明:高超声速双锥边界层分离与激波干扰流动的数值模拟结果对网格具有很强的敏感性,过稀的网格将产生严重的分离流动预测偏差;低耗散性的局部熵修正方法能更加准确地模拟复杂的高超声速激波与边界层分离流动干扰现象。     

高超声速飞行器流动特征分析

在非流线型构件或突起物的扰动效应、高马赫数和低雷诺数极限效应、低湍流度环境效应和由激波或摩擦导致的气动加热效应等4个方面的影响下,未来高超声速飞行器涉及的流动主要表现出这样的特点:典型流动结构强度高、尺度大,如强激波和厚边界层;局部流动结构数量多;激波、膨胀波和边界层结构之间相互干扰十分严重;转捩、压力脉动和一些流动结构对细微因素非常敏感;压力、摩擦应力和热流峰值现象普遍;升阻比屏障难以突破;流场同时依赖大量无量纲参数和有量纲参数,导致实验模拟难度大。本文在回顾传统高超声速流动主要流动现象的基础上,对上述7个方面涉及的典型流动现象的基础研究现状、问题本质和因果关系进行综合描述,讨论如何更有效地面对基础研究和工程实际问题。该文既可为解决典型流动现象中尚未解决的基础研究提供帮助,也可为如何合理地利用有限的已知知识解决工程应用问题提供指导。     

高超声速边界层流动稳定性实验研究

高超声速条件下边界层转捩相关研究是近年来空气动力学领域的研究热点。通过采用基于纳米粒子示踪平面激光散射(Nano-tracer-based planar laser scattering,NPLS)技术以及温敏漆等测试技术,对高超声速边界层的流动稳定性开展了实验研究。通过NPLS技术,对圆锥边界层中第二模态波精细结构进行了测量,并基于时间相关的测量结果,对第二模态波的波长和频率进行了分析。针对裙锥边界层NPLS结果,计算了特定位置上功率谱空间分布结果,测量得到了高次谐波成分。通过温敏漆和NPLS结果,发现主导三角翼前缘边界层转捩的模态为行进横流模态,分析了该模态的特性,并且与Kulite传感器测量得到的频率进行比较。     

高超声速流动粘性干扰效应研究

基于边界层相似性理论,开展了高超声速流动粘性干扰效应的相关性分析,得到了粘性干扰造成的物面压力系数改变量与粘性干扰因子成正比的结论。通过对升力体简化外形的高超声速绕流流动的模拟,分析了粘性干扰效应对气动力/力矩的影响,计算结果表明:由粘性干扰导致的气动力系数的改变量与粘性干扰因子在特定范围内具有良好的线性关系。     

高超声速流动的气体吹除控制方法研究

为了在超燃冲压发动机直连式试验台上得到真实的非均匀试验气流,采用商业软件Fluent对高超声速流动的气体吹除控制进行了数值研究。在试验台喷管和发动机隔离段之间增加了吹除控制段,试验了4种吹除方案。常温空气从位于控制段顶部、侧部和底部的吹除口注入,吹除口包括长槽、短槽和圆孔等,吹除气体的流量为来流流量的8.3%~13.3%。以超燃冲压发动机自由射流计算结果为基准,比较分析了控制段出口截面马赫数分布。结果表明,"顶槽"方案能够得到较为理想的吹除效果,吹除参数优化后可为试验提供指导。     

高超声速滑翔飞行器表面加热特点研究

高超声速滑翔飞行器具有高机动和远程快速到达能力,是高超声速技术应用的前沿领域。快速准确地预测飞行器表面受热变化特征,对高超声速滑翔飞行器热防护系统设计十分重要。以基于锥导乘波构型的高超声速滑翔飞行器和跳跃飞行弹道为研究对象,对其表面不同特征区域的加热开展了快速预测方法研究;采用选用方法分析驻点辐射平衡温度随飞行弹道的变化规律;在弹道特征位置上针对飞行器前缘及上、下表面的加热情况进行分析,获得高超声速滑翔飞行器受热的整体特征,可为分区域选择热防护措施提供参考。     

均匀磁场中高超声速弱电离气体流动数值模拟

采用7组元化学模型并应用组分公式计算电导率,通过求解黏性MHD(magneto-hydro-dy-namics)方程组,研究了不同强度均匀磁场对三维钝头体高超声速绕流化学非平衡流动的影响.结果表明,随着磁场的增强,激波脱体距离逐渐增加;总阻力系数和壁面温度逐渐减小.在B_y=0.03 T磁场作用下,与无磁场的结果相比,化学非平衡流中的激波脱体距离增加约7%,总阻力系数减小约5%,局部壁面温度最大降低74%;而冻结流中的激波脱体距离增加约43%,总阻力系数减小约6.9%,局部壁面温度最大降低18%.      

高超声速进气道强制转捩流动的大涡模拟

吸气式高超声速飞行器常在进气道边界层内布置粗糙颗粒或涡流发生器强制流动转捩为湍流以确保发动机正常启动。为了清晰认识强制转捩过程,采用隐式大涡模拟方法,对强制转捩问题开展了数值模拟研究。针对钻石形和斜坡形涡流发生器,计算得到涡流发生器诱导的流动结构,显示出强制转捩流动由涡流发生器产生的反向旋转流向涡对主导。扰动沿流向增长和发展,导致流向涡对以偶模式或奇模式失稳,偶模式失稳产生对称形式的涡对破碎,而奇模式失稳则导致非对称（弯曲）形式的涡对破碎。流向涡对破碎后产生一系列发卡涡并最终促使边界层转捩为湍流。最后就计算网格和数值耗散对隐式大涡模拟结果的影响以及计算的收敛性进行了讨论。     

含摩擦和加热的跨声速流动数值模拟研究

喷管中跨声速流动的反应流体动力学是实验和理论研究的一个重要领域.本文主要由两部分构成,在第一部分中对文献[2]中的Murman格式[1]的一般形式进行了简化,得到了守恒形式的欧拉方程,并对其求解方法进行了讨论,同时针对Chattot等的工作展开了分析,并给出了带激波的流动计算结果;而文章的第二部分则研究了包含摩擦、加热以及耗散的跨声速流动,并重点关注了一步Arrhenius化学反应模型下的爆震现象.本文还探究了摩擦力和带燃烧的跨声速流动之间的相互作用,同时详细给出了新的欧拉方程数值求解算法并同已有文献算例进行了对比验证.     

高超声速流动分区对接网格算法研究

本文发展了一种能适应强间断流场的满足通量守恒的分区对接网格算法，并给出了保证通量守恒的通量插值系数应满足的条件。采用有限差分NND格式，求解层流NS方程，数值模拟了马赫数6—20绕圆柱的高超声速粘性流场。数值试验表明：方法具有良好的稳定性和收敛性，强激波能自由地穿越内边界。     

高超声速前缘空腔数值模拟研究

介绍了一种前缘空腔高超声速被动流动控制方法并将其应用于轴对称钝锥外形,采用高分辨率的CFD数值模拟方法对不同参数的前缘空腔外形进行了对比研究,阐述了高超声速前缘空腔降低前缘壁面热流的机理,总结了不同空腔参数对空腔外形气动力、热的影响规律。结果表明,合理设计空腔参数可以使空腔外形的前缘壁面热流有较大程度的降低,并且不影响外形的升阻比等气动特性。     

高超声速飞行器气动加热三维数值分析方法研究

高超声速飞行器处于高空高速飞行环境,表面气动加热现象十分严重.有效预测并降低飞行器的表面温度,对防热材料和结构提出要求是高超声速飞行器设计的一个关键问题.采用三维N-S方程、Mac-Cormack中心差分格式对类乘波体构型的高超声速飞行器全机气动加热进行了数值研究.分析了不同马赫数来流对气动加热的影响.研究结果表明,采...     

高超声速壁湍流入口条件生成方法的比较

给定恰当的入口条件是开展壁湍流数值模拟的关键问题。选取基于有限差分的直接数值模拟方法结合高精度的计算格式,详细讨论了自然转捩、波纹壁面结构导致的Bypass转捩和利用时间发展湍流场进行参数回收促发转捩这几种湍流入口生成方法在高超声速条件下的可行性,分析了这几种方法各自存在的优点和不足之处。计算结果表明,Bypass转捩和利用时间发展流场进行参数回收方法相比自然转捩能快速促发转捩,但自然转捩得到的湍流场品质更好。该项研究为高马赫数壁湍流数值模拟入口条件的选取提供参考依据。     

高超声速再入飞行器头部辐射加热特性研究

高超声速再入飞行器的热环境特性是热防护设计的重要基础。它主要包括辐射加热和对流加热,一般情况辐射加热较小可以忽略,但部分再入飞行器在低空还具有极高的速度,辐射加热不能忽略。本文针对高超声速再入飞行器驻点的辐射加热特性开展了研究,采用高温空气辐射加热的四光谱带模型,同时计算方法中考虑了非平衡辐射加热及其"截断效应"。计算结果表明,在低空高速再入阶段,辐射加热比较明显,热防护设计需要考虑辐射加热。辐射加热主要由可见连续光谱辐射构成,红外线谱和紫外线谱的辐射加热也比较明显,紫外连续谱辐射以及非平衡辐射效应均不明显。辐射加热沿球头驻点至球头外边缘逐渐减小,影响辐射加热的主要因素有头部半径、飞行高度和飞行速度。     

考虑空腔的高超声速多流动区域同步数值模拟

先进的高超声速飞行器具有薄壁空腔结构,在飞行过程中受热会产生空腔内气体流动现象,从而影响流场和结构的温度分布。采用数值方法准确模拟高超声速流场、结构温度场和空腔内流动对热结构分析是很有必要的。以研究空腔流动对结构温度分布影响为目的,发展了一种适用于多流动区域流场/结构温度场耦合问题的同步计算方法,并以高超声速带空腔结构物体为例,数值研究了其外部气动热/结构热传导引起的空腔热对流问题。以已发展的高超声速外流场/结构温度场同步计算方法为基础,为了进一步考虑空腔内低速流场,采用了预处理矩阵方法。在流场与结构温度场的交界面两侧分别引入虚拟单元,从而高效地实现相邻场之间物理信息交换。首先通过标准算例验证了方法在求解单独气动热/结构热传导问题以及空腔自然对流问题中的准确性。进而对封闭和带有开孔的两种高超声速运动圆环分别进行多流动区域同步数值模拟。计算结果表明,由于结构温度不均匀引起的空腔内热对流反之也会对结构温度场分布产生轻微的影响。在空腔内气体流动的影响下,封闭圆环的前缘温度在35s内最多下降0.8%左右。对于带开孔空腔的圆环,其孔壁周边温度在0.5s内能够超过外流前缘驻点温度。     

绕HERMES外形高超声速无粘流动的数值模拟

本文采用时间相关法和隐式推进迭代方法求解了绕ＨＥＲＭＥＳ外形高超声速无粘流动。将抛物型方程和代数方程相结合生成计算网格，既节省了网格生成的时间，又能得到较合理的网格分布。在时间相关法计算中，采用了显式ＮＮＤ格式；在隐式推进迭代计算中，采用了推广的隐式ＮＮＤ格式。     

有局部稀薄气体效应的高超声速流动数值模拟

近空间高超声速飞行器当飞行高度和速度足够高时,其流场计算可能要考虑稀薄气体效应,传统的计算流体力学(CFD)方法预测的阻力和升阻比将不够准确。而现有的模拟稀薄气体流动的计算方法由于其计算量巨大,难以在工程实际中应用。因此需要发展能用于近空间高超声速飞行器流场的可行、可靠的计算方法。陈杰和赵磊在文献[1]中针对边界层中既有强剪切而气体分子自由程又相对较大的情况进行分析,提出了刻画此类局部稀薄效应的无量纲参数Zh,并提出了在传统CFD中通过采用依赖于Zh参数的等效黏性系数考虑局部稀薄效应对阻力计算影响的研究思路。因此,本文尝试将此等效黏性系数纳入CFD模型中,以在70km高空,以马赫数15飞行的小迎角钝平板为例,来检验计算方法是否合理可行。结果表明:和传统的CFD方法所得结果相比,新模型计算的阻力减小,升阻比增加,其改进的方向与现有飞行试验结果定性相符,且所增加的计算时间非常有限,可方便地应用于现有的计算空气动力学中。     

高超声速飞行器一体化方法研究

针对高超声速飞行器一体化设计方法现状的分析,阐明了吸气式发动机与乘波体飞行器之间高效的一体化对于高超声速飞行的重要作用,并从理论、原理、设计方法 3方面进行介绍。在激波理论方面,通过从直线激波的求解拓展到二次曲面激波的求解,为3维曲面激波的研究提供了帮助;在乘波原理方面,将乘波原理从外流乘波拓展到内流乘波,继而提出1种兼顾内外流需求的双乘波原理,深化了乘波原理的内涵;在设计方法方面,对于基本流场的气动设计问题,提出更加高效的一体化气动反设计方法。综上分析并归纳出准3维内外流一体化乘波理论与方法,从而在现有"准3维"研究体系上,构建并完善了全3维内外流一体化乘波理论与方法,对于复杂3维超声速内外流一体化设计技术的发展具有一定借鉴作用。