Hypersonic shock wave and turbulent boundary layer interaction in a sharp cone/flare model

A direct numerical simulation of hypersonic Shock wave and Turbulent Boundary Layer Interaction(STBLI) at Mach 6.0 on a sharp 7° half-angle circular cone/flare configuration at zero angle of attack is performed. The flare angle is 34° and the momentum thickness Reynolds number based on the incoming turbulent boundary layer on the sharp circular cone is Re_θ = 2506. It is found that the mean flow is separated and the separation bubble occurring near the corner exhibits unsteadiness. The Reynolds analogy factor changes dramatically across the interaction, and varies between 1.06 and 1.27 in the downstream region, while the QP85 scaling factor has a nearly constant value of 0.5 across the interaction. The evolution of the reattached boundary layer is characterized in terms of the mean profiles, the Reynolds stress components, the anisotropy tensor and the turbulence kinetic energy. It is argued that the recovery is incomplete and the near-wall asymptotic behavior does not occur for the hypersonic interaction. In addition, mean skin friction decomposition in an axisymmetric turbulent boundary layer is carried out for the first time. Downstream of the interaction, the contributions of transverse curvature and body divergence are negligible, whereas the positive contribution associated with the turbulence kinetic energy production and the negative spatial-growth contribution are dominant. Based on scale decomposition, the positive contribution is further divided into terms with different spanwise length scales. The negative contribution is analyzed by comparing the convective term, the streamwise-heterogeneity term and the pressure gradient term.     

Numerical simulation of compression corner flows at Mach number 9

A hypersonic vehicle encounters a wide range of conditions during its complete flight regime. These flight conditions may vary from low to high Mach numbers with varying angles of attack. The near-wall viscous dissipation associated with flows at combined high Mach and Reynolds numbers leads to significant wall heat transfer rates and shear stresses. The shock wave/boundary-layer interaction results in a flow separation region, which commonly augments total pressure losses in the flow and lowers the efficiency of aerodynamic control surfaces such as fins installed on a vehicle. The standard turbulence models, when used to resolve such flows, result in incorrect separation bubble size for large separated flows. Therefore, it results in an inaccurate aerodynamic load, such as the wall pressures, skin friction distribution, and heat transfer rate. In previous studies, the application of the shock-unsteadiness correction to the standard two-equation k-ω turbulence model improved the separation bubble size leading to an accurate pressure prediction and shock definition with the assumption of constant Prandtl number. In the present work, the new shock-unsteadiness modification to the k-ω turbulence model is applied to the hypersonic compression corner flows. This new model with variable Prandtl number is based on the model parameter, which depends upon the local density ratio. The computed wall pressures, heat flux and flow field are compared to the experimental data. A parametric study is carried out by varying compression deflection angles, free stream Reynolds number and wall temperatures to compute the flow field and wall data accurately, particularly in the shock boundary layer interaction region. The new shock-unsteadiness modified k-ω model with variable Prandtl number shows an accurate prediction of initial pressure rise location, pressure distribution in the plateau region and heat flux in comparison to the standard k-ω model.     

Behaviour of hypersonic wind tunnels diffusers at low Reynolds numbers

This paper deals with the behaviour of hypersonic wind tunnels diffusers at low Reynolds numbers. The phenomena that occur in the diffusers at relatively low Reynolds numbers are particularly critical due to the interactions between shock waves and the large boundary layers present in the long throat, which can reduce the maximum recovery pressure and therefore the efficiency of the diffuser. The aim of the present paper is to correlate the results obtained by numerical simulations of the diffusers with the experimental results, obtained using a high enthalpy blow down arc jet facility. The numerical and experimental results are in agreement and show that, for the range of Reynolds number investigated, the diffuser efficiency is smaller compared to the values determined for similar diffusers operating at larger Reynolds numbers.     

基于Reynolds平均的高超声速二方程湍流模型

 在Reynolds平均的框架下推导了可压缩湍流Reynolds应力方程和湍动能方程。根据一定的假设和尺度分析简化并封闭了所推导的湍动能方程。在均匀湍流假设下,湍动能耗散率可分解成可压缩性耗散和旋度耗散,并对其中的可压缩性耗散进行了封闭;同时认为旋度耗散不受可压缩性影响,直接引用传统的Reynolds平均不可压缩湍动能耗散率模型方程。由此构造了适用于高马赫数的二方程可压缩湍流模型。应用所发展的模型计算了高超声速平板绕流,并与若干现有模型的计算、实验与半经验公式的计算结果进行了对比,验证了所发展的模型。在此基础上,通过对压缩拐角的高超声速湍流的数值模拟,对所发展的湍流模型,以及若干现有模型进行了对比,显示了不同湍流模型及可压缩性修正在计算壁面压力分布和热流分布上的特点,说明了湍流模型可压缩修正的必要性,得到了所发展模型的计算结果最接近实验结果的结论。     

A scaling procedure for measuring thermal structural vibration generated by wall pressure fluctuation

This paper attempts to develop a scaling procedure to measure structural vibration caused simultaneously by wall pressure fluctuations and the thermal load of hypersonic flow by a wind tunnel test. However, simulating the effect of thermal load is difficult with a scaled model in a wind tunnel due to the nonlinear effect of thermal load on a structure. In this work, the temperature variation of a structure is proposed to indicate the nonlinear effect of the thermal load,which provides a means to simulate both the thermal load and wall pressure fluctuations of a hypersonic Turbulent Boundary Layer(TBL) in a wind tunnel test. To validate the scaling procedure,both numerical computations and measurements are performed in this work. Theoretical results show that the scaling procedure can also be adapted to the buckling temperature of a structure even though the scaling procedure is derived from a reference temperature below the critical temperature of the structure. For the measurement, wall pressure fluctuations and thermal environment are simulated by creating hypersonic flow in a wind tunnel. Some encouraging results demonstrate the effectiveness of the scaling procedure for assessing structural vibration generated by hypersonic flow. The scaling procedure developed in this study will provide theoretical support to develop a new measurement technology to evaluate vibration of aircraft due to hypersonic flow.     

SST湍流模型在高超声速绕流中的改进

为模拟高超声速湍流问题,对剪切应力输运(SST)湍流模型系数进行了修正。数值格式采用改进的总变差递减(TVD)格式,并对湍流模型的负值强制项进行了隐式处理。在此基础上计算了绕平板以及具有分离、再附、激波/边界层干扰等复杂流动结构的压缩拐角的高超声速流动。计算结果与试验数据及半经验公式的对比表明:SST湍流模型引入的雷诺剪切应力与湍动能之比为常数(Bradshaw数)在高超声速绕流中并不成立。Bradshaw数修正后的SST湍流模型与原模型相比,所计算的壁面压力、摩擦阻力和壁面热流分布更接近试验结果。     

高超声速风洞变雷诺数试验技术研究

为满足高超声速飞行器气动力雷诺数效应研究需求,在CARDC的Φ1米高超声速风洞中开展了变雷诺数试验技术研究.该项试验技术是利用Φ1米高超声速风洞采用高压下吹-真空抽吸驱动运行方式、风洞运行参数范围宽的特点,通过宽范围内调节风洞运行总压而大幅改变模拟雷诺数.研究采用了单点变雷诺数试验技术和连续变雷诺数试验技术两种手段来开展高超声速飞行器气动力雷诺数效应模拟.单点变雷诺数试验是通过一系列不同雷诺数条件、不同试验车次的试验结果,获得气动特性随雷诺数的变化规律;连续变雷诺数试验时,控制风洞总压从高到低连续变化,测量获取模型处于某一姿态角条件时气动力随雷诺数的变化规律.本文介绍了变雷诺数试验的风洞开车方式、试验及数据处理方法等,并开展了某升力体飞行器和某弹头模型雷诺数效应试验研究.研究结果表明:采用单点和连续变雷诺数试验技术相结合的方式,能较为完整、准确地获得飞行器模型气动力随雷诺数的变化规律.     

一种二元高超声速进气道起动特性的尺度效应研究

对一种二元混压式进气道三维流场进行了数值和实验考察,研究了不同尺度进气道模型自起动性能的变化。结果表明,在相同来流单位雷诺数条件下,随着模型尺度的减小,进气道自起动马赫数有所提高,起动性能有所降低。同时对不同尺度模型进行雷诺数匹配,发现在相同雷诺数下,不同尺度模型的起动性能相近,表明雷诺数是影响不同缩尺模型起动性能不同的主要原因。在可获得的实验结果范围内,数值模拟所得到的自起动结果基本与之相符。此外,对实验中发现在低雷诺数下进气道反而呈现出自起动特征的异常现象进行了初步分析,通过数值模拟比较指出了低雷诺数下来流偏向层流流态,可能会导致进气道呈现一种“起动”状态。     

宽高比对二元高超声速进气道性能的影响

对相同迎风面积、不同宽高比的二元高超声速进气道在设计马赫数6.0和非设计马赫数下的三维流场进行了数值模拟,研究了宽高比对进气道流场特征及性能参数的影响。结果表明,随着宽高比的增加,由于进气道长度和气流浸润面积的变化,内压段进口总压恢复系数、进气道流量系数和内部阻力系数逐渐降低;由于侧壁附近三维流动区域占整个流场的比例不同,当宽高比较小时,侧壁附近三维流动效应对进气道性能影响显著,进气道的总压恢复系数相对降低、增压比升高、温升比升高、出口马赫数降低,小宽高比进气道的低马赫数起动性能趋于恶化;设计马赫数下,宽高比的增加使二元高超声速进气道的反压承受能力降低。     

粗糙壁热交换实验研究

 本文介绍了来流马赫数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,看来这可能与驻点局部外形变化有关。     

几何尺寸对高超声速进气道气动性能的影响

为了探索模型缩尺比对高超声速进气道气动性能的影响,对不同缩尺比的二元高超声速进气道开展了数值模拟研究,结果表明：随着缩尺比的增大,进气道流量系数、隔离段出口总压恢复系数和马赫数均逐渐增大,而静压比逐渐减小,且来流马赫数越高,上述参数变化幅度越大。由理论与数值模拟分析可知,上述现象主要是由于不同缩尺比下,进气道当地雷诺数不同,导致进气道附面层相对厚度变化,进而影响进气道气动性能。理论分析了进气道总压恢复系数与缩尺比的定量关系,就进气道而言,进气道进口处附面层相对厚度减小1%,隔离段出口总压恢复系数提高约0.7%。     

非对称来流隔离段流动特性研究

1引言隔离段是超燃冲压发动机的一个重要气动部件,在进气道与燃烧室之间构建一气动热力缓冲区域,为进气道提供一个较宽的连续工作范围。通常超燃冲压发动机位于高超声速飞行器下表面的中后部,这样自由来流需要经过一段较长的前体后进入进气道,这就造成进气道进口靠近机体一侧存     

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

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

基于磁流体控制的高超声速进气道黏性效应

建立引入电磁源项的二维低磁雷诺数磁流体动力学(MHD)方程组,对高超声速二维前体/进气道黏性流场进行了数值模拟.在给出了进气道高于设计马赫数的非设计工况下黏性流场的基本特征基础上,进一步分析了施加MHD控制对进气道黏性效应的影响.结果表明:施加MHD控制可以有效抑制非设计工况下内进气道表面的附面层分离,改善上壁面的热状况,平衡上、下壁面之间的热负担;黏性作用下,进气道流场及性能参数随磁感强度的变化规律与无黏模型计算结果存在较大差别,对磁流体控制的高超声速进气道研究不可忽略黏性的影响.      

高温化学非平衡湍流边界层脉动量象限分析

高超声速飞行器在较低空域以极高马赫数飞行时，表面会同时存在湍流与化学非平衡流动，但目前针对此类高温化学非平衡湍流边界层流动特性的研究工作还比较有限，对不同湍流特征的主导流动机制的认识还有待于进一步深入。选取高超声速楔形体头部斜激波后的流动状态，设置3种不同的壁面温度，通过直接数值模拟对比了不同壁温条件下的边界层参数分布特性，并采用象限分解技术分析了边界层不同象限流动事件对雷诺剪切应力、湍流热流、湍流质量扩散的贡献。结果显示：在整个边界层中上抛和下扫运动对雷诺剪切应力的贡献占优；冷壁效应会使得流向和法向湍流热流通量的主导流动事件在温度峰值两侧发生改变。O原子组分流向湍流组分扩散主要受到高质量分数流体慢速运动事件和低质量分数流体快速运动事件的影响，而法向湍流组分扩散则主要受到高质量分数流体向上运动事件和低质量分数流体向下运动事件的影响。     

三棱柱的低雷诺数绕流研究

数值模拟了截面为等腰三角形的棱柱在不可压牛顿流体中的二维绕流问题,分析了不同的截面尺度比以及三角形顶点迎风和背风安置时的尾流稳定性与阻力系数.研究发现对于较大截面尺度比的棱柱,在相同雷诺数的情况下顶点迎风相较于顶点背风的情况更易失稳从而产生侧向力,但阻力系数相对较小,当雷诺数为160时,最多能够减少19.4%.进一步地对比研究了顶点迎风与顶点背风在不同截面尺度比时升、阻力系数及涡脱落频率随雷诺数的变化关系.      

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.     

Accurately predicting hypersonic transitional flow on cone via a symmetry approach

A new algebraic transition model is proposed based on a Structural Ensemble Dynamics (SED) theory of wall turbulence, for accurately predicting the hypersonic flow heat transfer on cone. The model defines the eddy viscosity in terms of a two-dimensional multi-regime distribution of a Stress Length (SL) function, and hence is named as SED-SL. This paper presents clear evidence of precise predictions of transition onset location and peak heat flux of a wide range of hypersonic Transitional Boundary Layers (TrBL) around straight cone at zero incidence, to an unprecedented accuracy as validated by over 70 measurements for varying five crucial influential factors (Mach number, temperature ratio, cone half angle, nose Reynolds number and noise level). The results demonstrate the universality of the postulated multi-regime similarity structure, in characterizing not only the spatial non-uniform distribution of the eddy viscosity in hypersonic TrBL on cone, but also the dependence of the transition onset location on the five influential factors. The latter yields a novel correlation formula for transition center Reynolds number which takes similar functional form as the SL function within the symmetry approach. It is concluded that the SED-SL model simulates TrBL around cone with uniformly high accuracy, and then points out to an optimistic alternative way to construct hypersonic transition model.     

Numerical simulations of shock/boundary layer interactions using time-dependent modeling techniques: A survey of recent results

Supersonic or hypersonic flows within and around flight vehicles inevitably involve interactions of strong shock waves with boundary layers. Flows within inlet/isolator configurations, and flows induced by control surface deflections are some examples. Such interactions are time dependent in nature and are often subject to low-frequency, large-scale motion that induces local pressure and heating loads. With recent increases in available computer power, it has now become possible to simulate such interactions at experimentally relevant Reynolds numbers using time-dependent techniques, such as direct numerical simulation (DNS), large-eddy simulation (LES), and hybrid large-eddy simulation/Reynolds-averaged Navier–Stokes (LES–RANS) methods. This paper will survey some recent work in this area and will describe insights in shock/boundary layer interaction physics gained by using these high-fidelity methods. Attention will be focused on studies that compare directly with experimental data at the same (or nearly the same) Reynolds number. Challenges in the application of these techniques to even more complicated high-speed flow fields are also outlined.     

高超声速化学非平衡粘性激波层方程的数值解

本文讨论了高超声速粘性激波层方程数值计算时差分格式引起的物理失真问题。具体分析了全隐格式格式粘性的影响,并作了数值试验。为了验证隐式结果的可靠性。在超声速激波风洞中测量了钝锥的表面压力分布,并与计算结果作了比较,两者基本一致。 本文采用隐式有限差分法数值计算了高超声速化学非平衡粘性激波层绕细长球锥的流动。计算时采用连续方程和法向动量方程耦合求解的方法以解决细长体远后身区计算中的问题。应用网格技术和加强系数矩阵主对角元素优势的方法提高了化学非平衡流计算的雷诺数范围。文中给出了高超声速化学非平衡流的计算结果,并与其它文献的结果作了比较。