高焓湍流边界层壁面摩阻产生机制分析

高超飞行器在中低空以极高马赫数飞行时,飞行器表面会遇到湍流与高温非平衡效应耦合作用的新问题.这种高焓湍流边界层壁面摩阻产生机制是新型高超声速飞行器所关注的基础科学问题,厘清此产生机制可以为减阻方法的设计提供指导,具有重要的工程实用价值.本文选取高超声速飞行时楔形体头部斜激波后的高焓流动状态,开展了考虑高温非平衡效应的湍流边界层直接数值模拟研究,并设置同等边界层参数下的低焓完全气体湍流边界层流动作为对比,采用RD (Renard&Deck)分解技术研究了高焓湍流边界层摩阻的主要产生机制,对摩阻产生的主要贡献项积分函数分布进行了详细分析,研究了高温非平衡效应对摩阻产生的影响规律;采用象限分析技术,研究了摩阻分解湍动能生成项的主导流动事件.计算结果表明,高温非平衡效应会使得壁面摩阻脉动条带的流向和展向尺寸均减小.分子黏性耗散项和湍动能生成项是高焓湍流边界层摩阻生成的主要流动过程.分子黏性耗散项主要作用在近壁区,高焓流动的分布与低焓流动存在差异.象限分析表明,上抛和下扫运动是影响摩阻分解中湍动能生成项的主导事件.     

基于多层分块算法的激波干扰流场预测

激波-激波干扰流场预测是超声速乃至高超声速流动中最具挑战性的问题之一. 特别地, 第IV类激波干扰由于其在壁面驻点附近产生极高的热载荷而备受关注. 本文针对圆柱诱导的弓形激波和入射斜激波的干扰问题, 分别基于量热完全气体模型和考虑振动激发的热完全气体模型, 数值求解有黏二维可压缩NS方程, 分析了高温气体效应对激波干扰流场结构, 以及第IV类激波干扰流场状态参数的影响. 接着, 本文基于一种具有广义可分离特性的遗传算法 (多层分块算法), 给出能够预测不同气体模型下第IV类激波干扰流场三波点的坐标位置、超声速射流的几何形状等特征性几何结构的数学模型, 进一步获得高温气体效应对激波干扰类型转变准则影响的定量化评估. 激波干扰类型转变准则面上的多组临界工况的激波干扰流场结构以及壁面压力和壁面热流分布的对比结果表明, 不同气体模型下的激波干扰类型和流场结构差异较为显著, 获得的定量化预测模型对工程中气动热环境的预测具有一定的参考价值.      

气体模型对再入飞行器气动力热环境影响的数值研究

复杂外形再入飞行器的设计,需对气动力热环境进行预测,由于不同的气体模型会对预测的结果产生影响,所以气动设计时就必须考虑这一影响.采用热化学平衡气体模型和双温度热化学非平衡气体模型对复杂外形再入飞行器的气动力热环境进行了数值计算;分析了气体模型对气动力、壁面热流等值线、驻点线平动温度、振动温度、组分质量分数等特征量的影响...     

极高超声速稀薄气体原子辐射效应的p-DSMC方法

极高超声速流动激波层内的高温导致内能模态的激发并伴随热辐射发生, 过高的温度使得空气分子完全解离, 原子组分对辐射热的贡献将达到80%以上. 本文基于优化的原子辐射模型, 提出追踪光子?直接模拟蒙特卡罗(p-DSMC)方法, 研究了稀薄流区不同马赫数下的高超声速二维圆柱绕流的壁面辐射加热, 获得了有无激发辐射效应的壁面压力和热流以及沿驻点线变化的平动、振动和转动温度. 在不考虑激发辐射效应的情况下, 得到的壁面压力和热流与已有的模拟结果符合的非常好, 误差均在5%以内, 尤其是在驻点位置, 误差在1%以内; 获得的平动、振动以及转动温度均与文献结果符合的很好. 在相同的来流条件下, 考虑辐射效应后发现, 来流速度低于10 km/s时, 辐射加热不明显, 在驻点区域, 辐射加热占对流加热比重在7%左右; 来流速度大于10 km/s时, 在驻点区域, 辐射加热占对流加热比重将超过30%. 考虑辐射效应后, 对非平衡区的平动、转动和振动温度的最大值影响不大. 此外, 另一个重要结论是, 流场中原子的浓度是影响壁面辐射热流大小的一个重要因素.      

高超声速高焓边界层稳定性与转捩研究进展

边界层由层流向湍流的转捩是高超声速飞行器设计面临的重大空气动力学问题. 随着飞行速域与空域的不断拓展, 高超声速高焓边界层中的高温气体效应会使得量热完全气体假设失效, 从而深刻影响流动转捩过程. 相关研究涉及多个学科, 是典型的多物理场耦合问题. 近年来, 随着相关飞行器技术的快速发展, 高超声速高焓边界层转捩问题的重要性越来越得到体现, 相关研究已成为国际上的热点领域. 本文综述相关研究进展, 首先介绍目前常用的高温气体物理模型, 尤其关注热化学非平衡模型, 并介绍激波捕捉、激波装配和边界层方程解等常用的高焓流动求解方法, 以及相关风洞和飞行试验技术的进展. 然后综述高温气体效应对转捩过程中的感受性、模态增长、瞬态增长和非线性作用等的影响的相关研究, 其中流向不稳定性中出现较大增长率的第三模态和超声速模态引起了广泛的研究兴趣. 最后进行总结, 并对未来发展略作展望.      

基于虚功率原理的鲹科鱼类波状推进分析

自然界中的大多数鱼类通过波状摆动的方式实现推进,这是动态变形的鱼体和周围流体相互作用的结果,研究推进中流体对鱼体变形的响应不仅可以增强对波状推进的认识,还可以为流动控制提供依据.以鲹科模式推进的仿生二维模型为研究对象,通过数值计算获得鱼体波动产生的流场以及鱼体受到的流体力数据.基于虚功率原理,将鱼体受力分解为4部分,分别是鱼体边界变速运动的瞬时贡献、流场中流体旋转和应变速率相对大小的贡献、壁面剪切应力的类摩阻分量和壁面摩阻分量.结果表明,当鱼体波状摆动产生推力时,鱼体边界变速运动是主要的正推力来源,并且该项80%的推力贡献来源于20%的鱼尾部分的边界变速运动.鱼尾两侧边界层中的流体旋转和应变速率的相对大小和壁面摩阻对推力都是负贡献.对于低雷诺数的情况,流体旋转和应变速率的相对大小的负贡献低于壁面摩阻的负贡献,而在高雷诺数的情况下,流体旋转和应变速率的相对大小的负贡献强于壁面摩阻的负贡献.壁面类摩阻分量相对于其他3项总是较小的.结合标度律分析,在摆动推进的标度关系中,与雷诺数无关的推力部分是由边界的变速运动、流场中流体旋转和应变速率共同提供,且流体旋转和应变速率也贡献了摆动推力中与雷诺数...     

尖锐前缘气动加热机理与激波脱体距离研究

针对尖锐前缘所处的气动环境,分析了由于前缘尖化所带来的稀薄效应影响。应用理论方法推导了尖锐前缘驻点气动热的计算方法和计算公式,并结合CFD方法、DSMC方法和Fay-Riddel工程方法对其前缘钝头的流场、驻点热流进行了计算和对比分析。计算结果表明,由于前缘的尖化使得稀薄气体效应提前出现,从而影响到前缘的激波厚度、激波形状和激波脱体距离等流动现象,导致激波结构复杂化,可能会对进气道唇口的斜激波带来不利的影响。另外稀薄效应的影响降低了驻点热流,缓解了尖化前缘的气动热环境。     

管道凸起激波/湍流边界层干扰的数值模拟

采用基于部分平均的湍流方程组,模拟了管道凸起流动中两种不同强度的激波-湍流边界层干扰现象,对流场的时均参数与实验值进行了比较,计算得到的壁面压力分布、摩阻系数分布、边界层厚度和速度型与实验值比较吻合很好.结果表明基于部分平均的湍流方程组既能够准确模拟小分离的激波/湍流边界层干扰流动又能够准确模拟大分离的激波/湍流边界层干扰流动,较好地预测了典型的λ激波结构.     

翼型跨声速气动特性的不确定性及全局灵敏度分析

针对马赫数和仰角的随机不确定性会导致气动性能波动的现象, 采用非嵌入式的混沌多项式方法对绕NACA0012 翼型跨声速随机气动特性进行不确定性及全局灵敏度分析. 具体分析了飞行状态的不确定性对气动载荷分布、流场及气动力系数的影响并通过全局灵敏度分析找出重要因素. 不确定性分析结果表明翼型上表面的激波以及激波后分离泡是造成气动性能剧烈波动的主要原因. 灵敏度分析结果表明在跨声速区域马赫数对激波处气动性能影响最大, 此外, 虽然马赫数和仰角相互耦合作用对气动力系数贡献比较小, 但对于激波位置处的流场, 这种互耦合作用不可忽略.      

触摸高温气体动力学

回顾了高温气体动力学与高超声速科技相关的一些重要研究进展,探讨几个具有基础性研究意义的方向：即高超声速流动模拟;高温气体热化学反应机制;高超声速流动滞止区预测;高超声速边界层转捩和激波/激波相互作用诱导的气动热问题．这些研究方向与高温气体效应和强激波密切相关,对高超声速科技关键技术的突破起着重要作用．     

Laminar-turbulent flow over wedges mounted on sharp and blunt plates

The flow pattern and the heat transfer on sharp and blunt flat plates near a wedge in a Mach 6 stream are experimentally investigated for two Reynolds numbers corresponding to the laminar and transitional states of the undisturbed boundary layer ahead of the wedge. It is shown that, as in a two-dimensional flow, plate bluntness leads to the attenuation of the heat transfer in the boundary layer/shock wave interference zone. However, when a certain threshold value of the bluntness is exceeded, a further increase in the bluntness has almost no effect on the heat transfer. For the first time, an experiment conducted in an intermittent (blow-down) wind tunnel has been based on the comprehensive use of panoramic (global) techniques for measuring the heat transfer and pressure coefficients and a method for visualizing the surface friction employing the luminescence effect after UV irradiation.     

A numerical study of oblique shock wave reflections over wedges in an ideal quantum gas

The various oblique shock wave reflection patterns generated by a moving incident shock on a planar wedge using an ideal quantum gas model are numerically studied using a novel high resolution quantum kinetic flux splitting scheme. With different incident shock Mach numbers and wedge angles as flow parameters, four different types of reflection patterns, namely, the regular reflection, simple Mach reflection, complex Mach reflection and the double Mach reflection as in the classical gas can be classified and observed. Both Bose–Einstein and Fermi–Dirac gases are considered.      

Self-similar hypervelocity shock interactions with oblique contact discontinuities

The two-dimensional inviscid refraction of a shock wave at an oblique contact discontinuity is self-similar; i.e. it depends only upon the variables , . We transform the compressible Euler equations into the self-similar () coordinates and solve the resulting boundary-value problem by an implicit equilibrium flux method. We present results for strong shock ( where M is the incident shock Mach number) interactions with an oblique contact discontinuity separating an inert gas from a gas which exhibits high-temperature gas chemistry effects. To model high-temperature effects we employ Lighthill's ideal dissociating gas (IDG) model. Comparison between the frozen and equilibrium limits, both of which are self-similar, indicate large changes in peak density and temperatures. Significant differences in the overall flow pattern between the frozen and equilibrium limits are observed for interfaces with low negative Atwood ratio and high positive Atwood ratio. Results from a local analysis are presented when the shock refraction is regular. The critical angle at which the transition from regular to irregular refraction occurs is slightly larger for the equilibrium chemistry case. Received 6 September 1996 / Accepted 20 May 1997     

Supersonic flow past a wedge on a flat plate. Laminar boundary layer separation and reattachment

Supersonic laminar flow past a two-dimensional “flat-plate/wedge“ configuration is numerically investigated. The pressures at the boundary layer separation and reattachment points are calculated over wide Mach and Reynolds number ranges. The minimum angles of the wedge surface inclination at which a return flow occurs are determined. The results are presented in the form of generalized Mach-number-dependences of the theoretical pressure on the wedge surface initiating boundary layer separation and the pressure at the boundary layer reattachment point.     

The interaction of a shock wave with the boundary layer in a reflected shock tunnel

The influence of a nontotal reflection on the interaction of a reflected shock wave with the boundary layer in a reflected shock tunnel has been investigated. The calculating method of the velocity, the temperature and the Mach number profiles in the boundary layer in reflected shock fixed coordinates has been obtained. To account for equilibrium real gas effects of nitrogen, the numerical results show that the minimum Mach number in the boundary layer has been moved from the wall into the boundary layer with the increasing of the incident shock Mach number. The minimum Mach number, the shock angle in the bifurcated foot and the jet velocity along the wall to the end plate are reduced owing to the increasing of the area of nozzle throat. The numerical results are in good agreement with measurements.     

Shock boundary layer interaction under the influence of a normal suction slot

The phenomenon of shock boundary layer interaction of a shock train under the influence of a normal suction slot is studied. In previous work, it was found that a normal, circumferential suction slot is sufficient to stabilize the primary shock of a shock train in as much as that the back pressure of the shock train can be increased until the shock train gradually changes into a single normal shock. Based on the experimental and numerical results, a flow model was derived which explains the transition of a shock train into a single shock under the influence of boundary layer suction. In this work, the normal shock boundary layer interaction model is validated against flow cases with different upstream Mach and Reynolds numbers. For that purpose three different nozzle flows are investigated at various total pressure levels. In a second step, the flow model is extended to the oblique shock case, correlating the suction mass flow with the total pressure distribution of the incoming boundary layer and the static pressure downstream of the oblique shock. Finally, the influence of the suction cavity pressure onto the shock boundary layer interaction is considered.     

Effect of boundary layer on Mach reflection over a wedge surface

Abstract. In this paper, we consider the phenomenon of unsteady Mach reflection generated by a plane shock wave advancing over a straight wedge surface, with particular attention to the deviation of the flow field from the self-similar nature. We examine the observed change in angle between incident and reflected shocks, which is in contrast to the fact that the angle should remain constant with time in a self-similar flow. The effect of the boundary layer behind the advancing shock wave over the surface of the wedge is considered to cause this, and boundary layer theory is utilized to estimate the thickness of the layer. It is found that the thickness increases as to the time t compared with t by the overall expansion in the self-similar flow. Assuming that the thicker boundary layer is effectively equivalent to a change in wedge angle, the effect of the boundary layer on the flow field should be less in later stages with larger t values in accordance with the observation above. Received 6 March 2000 / Accepted 23 April 2001     

Double Mach reflection of strong shock waves

The Mach reflection of shock waves in those cases in which the gas ideality condition is satisfied with high accuracy is well-known. The effects associated with the excitation of the internal degrees of freedom for the molecules lead to a qualitative change in the reflection pattern. The present study is an extension of [1, 2], devoted to the study of the Mach reflection of shock waves from a wedge under conditions in which the physical and chemical transformations in the gas heated by the shock wave play a significant role.     

Oblique shock wave with sweep

An attached planar, oblique shock with sweep is investigated for the inviscid flow of a perfect gas. The ratio of specific heats, freestream Mach number, and wedge angle in the plane of the freestream velocity are prescribed, with the sweep angle as a free parameter. Explicit relations are provided for jump and detachment conditions. A number of trends, some non-intuitive, are discussed, e.g., the downstream Mach number may increase with sweep.