高超声速圆锥边界层转捩数值研究

以类似于处理湍流脉动的方式来模型化非湍流脉动,考虑湍流的间歇性特征及引起超声速流动失稳机制,研究包括转捩起始位置及转捩区范围的高超声速流动转捩问题,发展可反映扰动模态和可压缩性影响的两方程SST转捩/湍流模式.针对尖锥模型,在马赫数6、攻角O°～2°时,采用第二模态(Second-mode)机制,预测并比较完全层流态、完全湍流态,以及自然转捩发展过程中的表面热流、表面摩阻等气动特性.     

高超声速边界层的转捩问题

取40km高空处的气体参数,对来流马赫数为8、10和12等三种不同情况的平板边界层和楔角为20°的楔型体边界层,做了转捩发生位置的预测。所用方法为改进的eN方法。结果发现,壁面条件对是第一还是第二模态波决定转捩位置有很大关系。如对于传统的eN方法,如果取N为10作为转捩判断标准,对马赫数为12的楔体边界层,等温壁条件下由第二模态波决定的转捩位置距前缘13m左右,而平板边界层中两模态波所决定的转捩位置距前缘均超过50m。而采用改进的eN方法,平板边界层中两模态波所决定的转捩位置距前缘均超过10m,而对马赫数为12的楔体边界层,两模态波所决定的转捩位置距前缘均出现小于5m的情况。     

高超声速进气道边界层强制转捩试验

在FL-31高超声速风洞分别开展了进气道的自然转捩和强制转捩风洞试验,试验Ma数为5、6和7,迎角为1°。通过红外热图得到了壁面的热流分布,从中得到了转捩区域。强制转捩装置为钻石型涡流发生器。随着涡流发生器高度的增加,强制转捩区域逐渐前移,得到了涡流发生器的有效高度,实现了强制转捩的目的。     

高超声速三维边界层横流转捩的数值研究

横流失稳是高超声速三维边界层转捩的主要机制。然而到目前为止,由横流失稳导致高超声速三维边界层转捩的数值研究还不多见。本文采用直接数值模拟方法计算了马赫数6、后掠角45°钝板边界层中定常横流涡的演化,在此基础上引入高频二次失稳模态,计算了二次失稳模态的非线性演化,直至湍流发生。计算结果表明:横流定常涡的非线性作用引起平均流修正,可使壁面摩擦系数曲线有一定程度的抬升。而高频二次失稳波的增长导致低频及定常扰动快速增长,促使壁面摩擦系数急剧抬升,同时饱和横流涡结构破碎,最终触发转捩发生。     

高超声速钝头体边界层转捩试验

为研究高超声钝头体边界层转捩以及头部钝度对转捩的影响,在FD-14和FD-14A两座激波风洞中开展了热流、压力扰动和高速纹影显示等综合测量。试验结果表明,转捩雷诺数关于钝度雷诺数的变化显示出转捩反转的趋势。压力扰动的功率谱密度（PSD）分析结果以流向离散分布云图形式显示,边界层高速纹影图像显示了第二模态波的发展、湍流的生成和熵层对边界层结构的显著影响。大头部钝度带来的强熵梯度熵层流动对边界层压力扰动频谱特性和流动结构影响显著,在转捩反转机理中起到重要作用。此外,马赫数对转捩的影响不容忽视。     

从总体设计角度透视高超声速飞行器边界层转捩问题

边界层转捩对飞行器总体性能影响很大,按照弹道式再入、再入机动、高超声速滑翔、高超声速巡航四类飞行器进行分类,与边界层转捩相关的典型问题可归结为气动力、气动热及掺混效率三类,相应带来飞行稳定性与配平能力、落点散布、防热风险、飞行器减重、推进系统优化等问题.针对边界层转捩控制问题,从总体设计方面可采取弹道设计、翼载荷控制、...     

基于边界层转捩的高超声速进气道特性研究

为了探索边界层非强迫转捩对进气道性能的影响,采用数值计算的方法开展了边界层转捩对轴对称混压式高超声速进气道流场特性的研究。研究表明:随着进气道中心锥锥尖钝化半径增大,边界层转捩先推迟。当锥尖钝度大到一定程度时,边界层转捩位置前移。随着钝化半径进一步增大,边界层转捩再次推迟,转捩位置逐渐后移。来流湍流度越大,边界层越不稳定,边界层转捩越易发生。与湍流边界层相比,考虑边界层转捩时进气道的总压恢复系数及流量系数较高、热载荷及阻力系数较小,Ma=6.5时喉道处总压恢复系数最高上升17.3%,进气道阻力最大下降17.4%。边界层转捩对壁面热流密度分布影响较大,但对壁面压力分布影响较小。钝化影响进气道的自起动性能,随着钝化半径增大,自起动马赫数升高,而边界层转捩对进气道自起动性能影响较小。     

小攻角高超声速钝锥边界层失稳特性

利用高阶紧致格式,采用直接数值模拟(DNS)和线性稳定性(LST)分析方法,对高超声速边界层的失稳机制和转捩特点进行了研究.通过对马赫数为6的2?攻角高超声速钝锥边界层的稳定性分析发现:小攻角高超声速钝锥边界层存在多枝不稳定模态;周向速度使钝锥的稳定性特征与不考虑周向速度时有本质的差别;转捩线在接近背风面处出现拐折现象是由失稳模态发生转换引起的.      

马赫数为6的高超声速钝锥湍流边界层空间演化的直接数值模拟

首先求解了来流马赫数为6的零攻角高超声速钝锥边界层的层流基本流场,在选定的计算域入口引入一组有限幅值的T-S波扰动,用高精度差分格式对流动进行了直接数值模拟.引入的扰动触发了转捩,从而得到了空间模式下的湍流边界层.研究了湍流平均场与脉动场的统计特性,给出了相干结构的流动显示图,并对强雷诺比拟的结论进行了检验.     

高超声速边界层转捩对旋转钝锥自由飞运动的影响

通过在钝锥模型表面上布置人工绊线促使边界层强迫转捩,采用运动自由度不受约束的风洞模型自由飞试验技术研究边界层转捩对高超声速旋转钝锥自由飞行运动特性和气动特性的影响规律,并与自然转捩的旋转钝锥风洞模型自由飞试验结果作对比分析,试验马赫数为5.0,以模型长为特征尺寸的自由流雷诺数为1.68×106。研究结果表明:有人工绊线的旋转钝锥在自由飞行过程中有"激励稳定"的绕流流场,产生动态稳定的自由飞运动(动稳定导数系数小于0),而无转捩绊线的旋转钝锥在自由飞行中则有"激励不稳定"的绕流流场,产生动态不稳定的自由飞运动(动稳定导数系数大于0)。     

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

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

粗糙元诱导的高超声速边界层转捩

基于有限体积方法,直接数值模拟了高超声速边界层内不同形状粗糙元导致的强制转捩现象;为了能够深入探究强制转捩机理,解析小尺度运动,同时又能够较好地捕捉激波,采用高阶色散最小耗散可调（MDCD）格式对Navier-Stokes方程组对流项进行重构。计算结果表明,数值结果与对应的实验值吻合较好;该方法能解析小尺度的流动结构以及规则结构的破碎与失稳过程,可揭示粗糙元引起的强制转捩机理,即此类强制转捩主要由粗糙元顶部的三维剪切层失稳导致。对多种粗糙元的转捩效果进行了定量研究,影响因素包括粗糙元形状、几何参数等。     

高超声速二维钝楔边界层参数分析

为了指导二元高超声速进气道型面设计,基于两类高超声速二维钝楔外形,采用高精度CFD方法计算了四种典型工况下二维钝楔边界层参数,分析了前缘半径、壁面温度、压力梯度对边界层参数的影响,得出了边界层位移厚度、边界层总压恢复及壁面压力梯度沿流向的变化规律。研究表明,前缘半径及壁面温度对边界层参数影响较大,在两者一定的情况下,逆压梯度使边界层位移厚度增速放缓,顺压梯度会使边界层位移厚度快速增长。     

高超声速锥柱裙模型边界层转捩的弹道靶实验

为研究高超声速边界层转捩现象、给边界层计算提供可靠的对比数据,在中国空气动力研究与发展中心趟高速弹道靶上开展了锥柱裙模型高超声速边界层转捩的自由飞实验。所采用的锥柱裙模型全长105mm,飞行速度1．94km／s（Ma=5．65）,单位雷诺数4．32×10^7～1．20×10^8m-1。使用激光阴影成像技术,获得了锥柱裙模型边界层转捩和湍流边界层发展的图像,测得的湍流边界层厚度在0．6～2．2mm之间,湍流涡的流向尺寸与边界层厚度的比值介于0．3～0．8之间且沿流向呈下降趋势。实验结果表明：弹道靶实验能够获得给定飞行环境下的高超声速边界层转捩图像,从图像中可以清晰判断转捩位置或区域、测量边界层厚度和分析湍流涡的尺寸。     

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

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