高超声速下表面凸起干扰气动热实验研究

 对高超声速飞行器表面凸起附近的气流流动和气动加热开展了实验研究和分析。实验在高超声速炮风洞中进行,来流马赫数为8.2、单位雷诺数为9.35×10⁶ m^-1。利用薄膜传热测量方法进行了凸起几何形状和边界层状态对干扰流动加热的影响评估。利用流油图谱和纹影摄像法得到了凸起周围的流动特征:若凸起上游边界层未分离,最大峰值热流发生在凸起侧方附近处;若凸起上游边界层完全分离,最大峰值热流通常发生在凸起的上游表面。实验发现最大峰值热流和来流边界层状态关系不大,原因是流动干扰区表现出较强的三维扰动特性,使得来流层流边界层在干扰区内会转变成过渡甚至完全湍流状态。     

An overview of infrared thermography techniques used in large wind tunnels

Infrared thermography (IRT) is used at Onera in large facilities for boundary layer visualization and for heat flux assessment. Modern IR cameras and insulating paints enable efficient visualization of the laminar/turbulent transition region. This technique is now applied in large transonic test facilities. Heat flux assessment is one of the main purposes of hypersonic tests. It is done mainly with IRT and dedicated softwares, while sensors as thermocouples are used to check the reliability of IRT. A 1D data reduction method has been developed to provide the heat flux through temperature measurements. It takes into account thickness and curvature effects. The method has been recently improved to be used with steel models covered with an insulating paint, which provides a high emissivity. The temperature film is converted into a heat flux film, which is be used to extract the useful information. This requires image processing tools that relate every pixel to a point on the model. A new application of IRT is going on in the Onera's high enthalpy hypersonic wind tunnel F4. The camera is used in single-line scan mode because of the short duration of the run. The main difficulty comes from the flow, which is not transparent. The first trial to cope with this kind of optical pollution is encouraging.     

激波风洞边界层转捩测量技术及应用

高超声速边界层转捩对摩阻、传热等有重要影响。在高超声速飞行器研制中,迫切希望能精确预测和控制边界层转捩。激波风洞作为高超声速气动热环境试验的主要地面模拟设备,是研究高超声速边界层转捩的重要设备。但激波风洞原有测量技术适用于工程型号试验,需要依据高超声速边界层转捩特点进行适应性改造和升级。依据高超声速边界层转捩过程中的热流、压力、密度等物理参数变化,发展了薄膜热流传感器测热技术、温敏热图测量技术、高频脉动压力测量技术、高清晰度纹影显示技术等适用于激波风洞的边界层转捩测量技术。并针对头部钝度0.05 mm的半锥角7°尖锥模型,在中国空气动力研究与发展中心Ø2 m激波风洞（FD-14A）马赫数10、单位雷诺数1.2×10⁷/m的流场条件下开展了边界层转捩试验。采用多种转捩测量技术同时进行测量,获得尖锥模型表面边界层转捩情况、边界层脉动压力频谱特征、边界层内清晰的第2模态波和湍流斑纹影图像,不同测量技术获取的试验结果可相互印证,线性稳定性理论分析结果与试验结果相吻合。     

平板空气舵干扰区烧蚀凹陷对热环境影响

针对高速飞行条件下空气舵干扰区烧蚀产生的局部凹陷对气动加热的影响问题,建立了平板-空气舵流动模型,针对典型高速飞行状态,采用高温热化学非平衡数值模拟,研究了空气舵缝隙区的流动结构和气动加热规律,并对舵缝干扰区的烧蚀外形进行了模化,分析了干扰区烧蚀凹陷对流动结构和气动加热的影响,结果表明:烧蚀凹陷改变了干扰区压力分布规律,降低了沿展向压力梯度,从而抑制了边界层的横向流动和厚度减薄效应,使得干扰区热流降低,且热流降低量值与烧蚀凹陷深度呈正相关,凹陷深度为5 mm时干扰区热流降低量达到28.9%。      

粗糙壁热交换实验研究

 本文介绍了来流马赫数5的条件下,典型球锥模型的粗糙壁热交换实验结果。模型头部半径R为27.4毫米,底部直径D为60毫米,对五个不同粗糙度的模型进行了实验。模型表面粗糙颗粒直径d分别为0、0.3、0.5、0.7、0.9毫米。风洞前室总压P_t为10～45公斤/厘米。,相应的来流雷诺数R_eD为（O.8～3.6）×10⁶。 实验结果表明:光滑壁模型表面是层流加热,驻点热流与层流理论计算值较一致。粗糙度的影响,在低总压条件下（10公斤/厘米）主要在于促使边界层的转捩和发展。随着风洞总压的提高,物面静压和局部雷诺数的相应增大,粗糙度对热流的影响才明显增强,而严重的区域在端头。在实验最大粗糙度和最大总压条件下（d=O.9毫米、p_t=45公斤/厘米。）,除驻点值外,热流与光滑壁层流驻点值相比(q_i/q_so)的峰值在音点区域且接近4,而在驻点,此模型有别于其它模型,较为特殊,比热流最大值接近6,看来这可能与驻点局部外形变化有关。     

隐式紧耦合SST和TNT湍流模型的高速流动数值模拟

将SST(shear stress transport)和TNT(turbulent/non-turbulent)湍流模型输运方程与平均流场控制方程进行隐式紧耦合求解,结合当地时间步长方法和湍流源项隐式处理确保求解过程的快速和稳定.采用AUSMPW+(AUSM by pressure-based weight functions)格式和LU-SGS(lower-upper symmetric Gauss-Seidel)隐式紧耦合方法对高超声速压缩拐角流动、锥柱裙流动和超声速非对称激波/边界层干扰问题进行了数值模拟.计算结果与实验值的对比表明:SST模型和TNT湍流模型可以很好地预测15°压缩拐角流动的壁面压力和热流密度;随着压缩拐角的增大,计算结果与实验值偏差增大;可压缩性修正对压缩拐角流动的压力和热流密度分布有很大影响,对超声速非对称激波/边界层干扰流动影响很小;隐式紧耦合方法比显式紧耦合方法具有更好的收敛特性.      

压气机叶栅叶片表面附面层流态变化影响因素探讨

以平面叶栅中的二元叶栅模型为试验对象,测量了在不同来流条件下叶片表面流场分布情况及栅后气流参数,分析了不同来流条件下叶片表面附面层流动状态的变化。并借助数值模拟手段重点研究了在不同来流马赫数和冲角下,叶片表面压力梯度对层流附面层向紊流附面层转捩过程的影响,通过利用实验数据分析研究来流条件对转捩过程的影响,为从机理上更深刻地认识叶片表面粘性附面层转捩机制提供了科学参考依据。     

高超声速飞行器-进气道一体化热流数值计算

采用CFD(计算流体动力学)技术, 开展了飞行器前体/发动机一体化气动热环境分析.对层流区、转捩区和湍流区分别采用计算模型, 在湍流区利用压缩性修正的SSGZ-Jk-ε湍流模型, 在转捩区引入代数型转捩因子模型描述边界层由层流逐渐过渡为完全湍流的流动过程.计算了前体和内通道的表面热流, 并与实验结果进行了对比.结果表明所采用的计算方法可以较好地预测前体及发动机内通道热流率, 流动状态、几何结构及激波入射对热流值影响较大.      

高超声速圆锥边界层失稳条纹结构实验研究

边界层转捩的准确预测是高超声速飞行面临的关键气动问题之一。为研究高超声速边界层失稳和转捩机理,以前缘半径1.6mm、半锥角7°的圆锥模型为研究对象,在FD-07高超声速风洞中采用红外热图技术开展边界层转捩实验测量。通过与工程计算结果对比,确认模型表面边界层流态。实验结果表明:有迎角条件下,模型表面中后段出现条纹结构,条纹结构的起始位置随着周向角的增加而向上游移动;随着迎角的增加,条纹起始位置向上游移动,条纹强度差异和条纹与模型中心线的夹角越来越大。实验获得的条纹结构与不同频率扰动波相互作用直接数值模拟获得的条纹结构现象一致。通过对比分析,认为边界层内不同频率扰动波相互作用是产生条纹结构的一种机制。     

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

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

高超声速气动热辨识技术实验研究

为估计高超声速飞行中的气动热流,采用热流辨识技术研究了平板表面的气动热流。首先采用仿真数据对辨识方法进行验证,结果表明,在测量噪声较低时热流辨识结果误差较小。此后,开展了平板的高超声速气动热风洞实验,通过内埋热电偶测量平板内部温度,采用热流辨识法获得平板表面的气动热流,采用辅助测点对比验证了内埋测点辨识结果的可信度。最后将辨识结果与平板边界层热流理论估算结果进行对比,分析了两者产生误差的原因。本文的研究为热流辨识技术的实际应用奠定了基础。     

基于雷诺平均Navier-Stokes方程的表面传热系数计算

采用有限体积法数值求解控制二维绕流的雷诺平均Navier-Stokes(RANS)方程组,计算了光滑和粗糙NACA0012翼型以及圆柱表面的局部表面传热系数.分析了近壁面网格间距、湍流模式和表面粗糙度模型对数值计算结果的影响.结果表明:切应力输运(SST)湍流模型能够区分层流和湍流边界层的对流传热特性,并能预测转捩的发生;采用Spalart-Allmaras(S-A) 扩展模型能够计算粗糙壁面的对流传热系数,但采用忽略转捩函数的S-A模型不能有效计算层流边界层的传热系数.当近壁面网格间距接近10-5量级的黏性子层时,在光滑和粗糙壁面都能得到准确的传热系数分布.结合合适的近壁面网格间距,湍流模式和表面粗糙度模型可以得到与实验数据十分接近的表面传热系数曲线.通过与求解不可压缩RANS方程得到的结果比较后发现,不可压缩RANS方程主要忽略了压缩和黏性耗散效应,这种效应可以通过绝热升温项的形式并入总体热分析.      

(高)超声速流动试验技术及研究进展

近年来,与高速飞行器相关的(高)超声速流动受到了极大的关注。这类流动所具有的非定常性、强梯度和可压缩性对试验方法和风洞设计技术提出了挑战。超声速纳米示踪平面激光散射(NPLS)技术是由作者所在团队研发的非接触光学测试技术。它能够以较高的空间分辨率来揭示超声速三维流场的一个瞬态剖面的时间解析的流动结构。介绍了NPLS技术以及基于NPLS开发的密度场测量、雷诺应力测量和气动光学波前测量等方法,并回顾了这些技术在超声速边界层、超声速混合层、超声速压缩拐角、激波/边界层相互作用和光学头罩绕流等流动中的应用,清晰地再现了边界层、混合层、激波等典型流场结构及其时空演化特性。另外,为了模拟和研究高空大气条件下边界层自然转捩和超声速混合层的转捩特性,介绍了高超声速静风洞、超-超混合层风洞的设计技术以及层流化喷管的设计方法。     

A swept fin-induced flow field with different height mounting gaps

In order to apply the air fin successfully and ensure the maneuverability of hypersonic vehicle, a key problem to be studied urgently is the heat flux brought by the fin mounting gap. The appearance of mounting gap and fin shaft can induce many complex flow structures which need more attentions to be investigated. Under Ma 6, Nano-tracer-based Planar Laser Scattering (NPLS) and Temperature Sensitive Paints (TSP) were applied to visualize and measure transient flow structures and heat flux distribution of a swept fin-induced flow field with different height mounting gaps. Complementarily, Reynolds-averaged N-S equations were solved with k-ω SST turbulent model. The heat flux distribution results of numerical simulation and TSP observed the change of high heat flux region with different mounting gap, both in position and magnitude. The streamlines based on Computational Fluid Dynamics (CFD) and flow visualization results obtained by NPLS revealed the cause of high heat flux region. The high heat flux region in this flow field is mainly related to the reattachment of vortex and flow stagnation. The increase of gap height can lead to stronger gap overflow and shaft-induced horseshoe vortex, which are source of the high heat flux around the fin. The case with the highest mounting gap (4 mm) en-counters the most severe aerodynamic heating, both on the surface of fin and plate. Thus, under the premise of ensuring the flexibility of the fin, the gap should be set as small as possible.     

Flap contour optimization for highly integrated SERN nozzles

The performance of a Single Expansion Ramp Nozzle (SERN) for hypersonic cruise vehicles is examined with focus on the flow at the nozzle flap and the corresponding separation of the boundary layer. Two wind tunnel models are used to analyze the interaction of the free stream and the nozzle flow and the flow at the flap in detail. A finite element flow code is employed to support the design of the flap model and the analysis of the flow.Comparative contour studies show that the size and the location of the separation at the flap cowl and therefore the aerodynamic quality can be influenced by optimizing the flap cowl contour. The position of the separation is measured for different contours and varying Reynolds numbers at laminar, transitional and turbulent flow conditions. Numerical studies of a nozzle/afterbody configuration suggest an improvement of up to 3.5% in axial thrust for optimized flap contours in the examined Mach range between 1.64 and 6. At the same time the rotation of the thrust vector is reduced by up to 5° at lower Mach numbers.     

壁面温度对斜爆震发动机进气道流动特性影响的数值研究

为研究壁面温度条件对层流、转捩、湍流状态下斜爆震发动机进气道流场结构、流场参数的影响,选取Ma10级、具有曲面压缩段的斜爆震发动机进气道为研究对象进行数值模拟,对进气道壁面附近激波诱导分离区、热边界层的变化进行了深入探讨。数值模拟结果表明,进气道肩部圆弧过渡段出现的再层流化现象,壁面冷却对其起抑制作用,绝热壁面条件下再层流化程度最为严重。壁面温度的增加有利于延缓流动转捩,同时也导致了分离区尺寸的增加以及转捩、湍流状态下分离区主体位置逐渐前移,进气道内通道的转捩为分离诱导转捩,转捩位置主要受到分离点位置的影响,整体表现为壁面温度增加转捩位置前移。进气道出口顶板侧热边界层厚度随着壁面温度的增加逐渐变厚,转捩状态下热边界层厚度变化可达5%,温度峰值也随着壁面温度的增加逐渐增加,且峰值位置逐渐靠近壁面。壁面温度条件相同时,层流状态下热流、热边界层厚度均较小。转捩、湍流状态下进气道出口顶板侧热边界层较厚,约为层流状态3倍,同时转捩、湍流状态下热边界层厚度相差可达2%。     

电弧加热流场湍流度对尖锥边界层转捩影响的研究

电弧加热流场的热环境特性直接影响热防护系统的地面试验数据,由于电弧加热器高温气流和参数波动的原因,直接测量湍流度非常困难。为研究电弧加热流场湍流度对于边界层转捩的影响,采用红外热图热像仪,在电弧加热流场中进行了5°尖锥模型边界层转捩研究。结合数值计算,将试验结果与常规风洞的尖锥边界层转捩结果进行了比较。结果表明：马赫数影响的雷诺数转捩判别准则可以用于计算电弧加热流场的转捩雷诺数;电弧加热流场的尖锥边界层转捩雷诺数显著小于常规风洞的转捩雷诺数,表明在该试验条件下,电弧加热流场的湍流度显著大于常规风洞。     

翼型风洞侧壁干扰的数值模拟研究

运用Navier-Stokes数值模拟对翼型模拟试验对风洞侧壁干扰进行模拟，将简单代数湍流模型扩展用于机翼／风洞侧壁结合区流动，分析了风洞实验侧壁干扰问题的形成机理和影响翼型实验侧壁干扰的各种因素，如翼型展长、风洞侧壁边界层厚度及侧壁边界层抽吸等，对实验结果的影响，得出了一些有用的结论。计算格式空间采用中心有限体积离散，时间采用多步Runge-Kutta时间步长格式进行积分。计算结果证明了该方法的可行性和优越性。     

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

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

基于LES方法的平板非定常激波/湍流边界层干扰研究

以高超声速发动机进气道湍流分离控制为应用背景,采用大涡模拟(LES)方法进行马赫数为3.0(唇口附近马赫数约为3.0)的激波/湍流边界层干扰(SWTBLI)流场机理研究.利用扰动循环引入的方法,先得到充分发展湍流场,然后根据斜激波关系式引入激波的方法进行激波/湍流干扰模拟.研究结果显示:充分发展湍流场在激波作用下产生逆...