高超声速溢流冷却实验研究

高超声速溢流冷却是一种新型的飞行器热防护方法,基本思想为:在高热流区布置溢流孔,控制冷却液以溢流方式流出,之后通过飞行器表面摩阻作用展布为液膜,形成热缓冲层以降低飞行器表面热流. 目前,溢流冷却技术还处于探索阶段,实现工程应用前还需开展大量的实验验证和机理研究工作. 本文首次开展溢流冷却的实验研究工作,采用热流测量、液膜厚度测量及液膜流动特性观测技术,搭建了完善的溢流冷却风洞实验平台,对溢流冷却热防护性能和高超声速条件下液膜流动规律进行了初步研究. 研究表明:(1) 高超声速流场中通过溢流能够在飞行器表面形成液膜并有效隔离外部高温气流,可降低飞行器表面热流率;(2) 楔面上的液膜前缘流动是一个逐渐减速的过程,增加冷却液流量液膜厚度变化不明显,但液膜前缘运动速度增大;(3) 液膜层存在表面波,在时间和空间方向发生演化,导致液膜厚度的微弱扰动;(4) 液膜层存在横向展宽现象,即液膜层宽度大于溢流缝宽度. 原因是液膜层与流场边界层条件不匹配,存在压力梯度,迫使冷却液向低压区流动,从而展宽液膜层,并且流量越高,横向展宽现象越明显.      

高超声速溢流冷却实验研究

高超声速溢流冷却是一种新型的飞行器热防护方法,基本思想为:在高热流区布置溢流孔,控制冷却液以溢流方式流出,之后通过飞行器表面摩阻作用展布为液膜,形成热缓冲层以降低飞行器表面热流.目前,溢流冷却技术还处于探索阶段,实现工程应用前还需开展大量的实验验证和机理研究工作.本文首次开展溢流冷却的实验研究工作,采用热流测量、液膜厚度测量及液膜流动特性观测技术,搭建了完善的溢流冷却风洞实验平台,对溢流冷却热防护性能和高超声速条件下液膜流动规律进行了初步研究.研究表明:(1)高超声速流场中通过溢流能够在飞行器表面形成液膜并有效隔离外部高温气流,可降低飞行器表面热流率;(2)楔面上的液膜前缘流动是一个逐渐减速的过程,增加冷却液流量液膜厚度变化不明显,但液膜前缘运动速度增大;(3)液膜层存在表面波,在时间和空间方向发生演化,导致液膜厚度的微弱扰动;(4)液膜层存在横向展宽现象,即液膜层宽度大于溢流缝宽度.原因是液膜层与流场边界层条件不匹配,存在压力梯度,迫使冷却液向低压区流动,从而展宽液膜层,并且流量越高,横向展宽现象越明显.     

微槽-吸气组合控制平板边界层多模态稳定性研究

边界层转捩会使高超声速飞行器壁面摩阻和热流显著增加,因此在高超声速飞行器设计过程中往往占据重要地位.针对高超声速飞行器多模态转捩控制问题,提出了微槽道(1 mm)与边界层吸气的组合控制方法,并通过直接数值模拟和线性稳定性理论研究了Ma=4.5平板边界层的稳定性及组合控制效果.边界层在无控状态时,同时存在失稳的第一、二模态波,且二维第二模态波最不稳定;单纯施加微槽道控制时,边界层第二模态波会被抑制但第一模态波会被略微激发.对比而言,采用“微槽-吸气”组合控制后,不仅增强了对第二模态波的抑制效果,而且减弱了第一模态波的激发程度;同时随着吸气强度的增加,第二模态波不稳定区域明显收缩、频率显著增高,而第一模态波则变化不明显.相较于单纯的微槽道,吸气增强了“微槽吸收”与“声波散射”作用,因此中等吸气强度下该组合控制方法对第一和第二模态波的增长率分别实现了12.63%和28.02%的抑制效果.以上结果表明“微槽-吸气”组合控制手段具有适用宽频、布置区域灵活的优点,展现出了一定的多模态控制效果.     

高超声速流场激光测速技术研究进展

高超声速气流条件下飞行器内/外部流动中存在强湍流及脉动、边界层转捩、激波-边界层干扰和高温真实气体效应等耦合效应,表征该非定常流动现象对飞行器气动力、气动热以及目标光电特性等产生的影响是高超声速流动研究中的前沿课题.速度作为表征流动过程最重要的参数之一,准确的速度测量对于深入理解上述复杂流动-传输机理以及高超声速飞行器设计具有重要指导意义.文章针对高超声速流场速度测量中几种常用的非接触式激光测试技术进行了综述,主要包括基于空间法的粒子图像测速,基于激光吸收光谱、激光诱导荧光和瑞利散射的多普勒测速,基于飞行时间法的分子标记测速,以及基于流场折射率的聚焦激光差分干涉测速技术.首先简要介绍每种激光测速技术的基本原理,然后进一步介绍该技术在高超声速自由流、层/湍流边界层、激波/边界层干扰、尾流或其他复杂流动区域的速度及其脉动度测量等方面的典型应用,分析各种技术环境适用性及面临的局限性和挑战.最后对基于激光技术的高超声速流场速度测量进行了总结及发展趋势展望.     

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

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

高超声速湍流边界层研究现状

一、高超声速湍流边界层研究的重要性随着再入导弹武器从惯性弹道导弹发展到可作机动飞行的多弹头分导弹道导弹,以及航天飞机的出现,高超声速再入飞行器的气动外形变得更复杂了。由于出现了多个激波的相互干扰,激波与边界层的相互影响,以及边界层的分离(入射激波和后台阶产生      

微管内气液两相流中液膜厚度的实验测量

微通道内气液两相流中气柱(plug bubble)与通道壁之间液膜厚度的实验测量,是微热管、微流动、微电子冷却以及气泡雾化等研究中普遍关注的问题.本文利用基于光学干涉和快速傅立叶变换的空间频谱分析方法,实验测量获取了含表面活性剂水中气柱在750μm 通道内运动时其与通道壁面之间的液膜厚度.实验结果表明:表面活性剂对液膜厚度的影响比较明显,即当表面活性剂浓度在一定范围内增大时,液膜厚度会减小;此外,当气柱运动速度在一定范围内增大时,液膜厚度也会减小.     

前体涡发生器对轴对称高超声速进气道激波振荡流动的影响实验

激波振荡是高超声速进气道不起动过程中常见的流动现象,会显著降低进气道气流捕获与压缩效率、产生剧烈的非定常气动力载荷而危害飞行器安全. 从激波振荡的控制出发,实验研究了前体转捩带位置的涡发生器对轴对称高超声速进气道激波振荡流动的影响. 分别在起动和激波振荡两种进气道流态下,选择无、0.5 mm与1 mm高度涡发生器工况进行对比研究. 并采用高速纹影与壁面动态测压同步记录非定常流动特征. 结果表明,1 mm高度内的涡发生器对起动状态的进气道主流流场结构、壁面压强分布影响不显著. 但对于激波振荡流动,涡发生器会明显缩小外压缩面分离区运动范围,缩短振荡周期,提升振荡周期内壁面压强的时均值. 涡发生器的影响程度随其高度的增大而增强,其中振荡周期从无涡发生器的4 ms缩短到1 mm高度涡发生器的3.13 ms. 此外,0.5 mm高度涡发生器会使得进气道内部测点的压强振荡幅值整体下降,相比无涡发生器工况的下降幅度可达23%. 流场结构与壁面压强信号的分析表明,涡流发生器主要通过其产生的流向涡影响激波振荡流动,包含流向涡对下游边界层的扰动以及流向涡与分离区的相互干扰.      

前体涡发生器对轴对称高超声速进气道激波振荡流动的影响实验

激波振荡是高超声速进气道不起动过程中常见的流动现象,会显著降低进气道气流捕获与压缩效率、产生剧烈的非定常气动力载荷而危害飞行器安全.从激波振荡的控制出发,实验研究了前体转捩带位置的涡发生器对轴对称高超声速进气道激波振荡流动的影响.分别在起动和激波振荡两种进气道流态下,选择无、0.5 mm与1 mm高度涡发生器工况进行对比研究.并采用高速纹影与壁面动态测压同步记录非定常流动特征.结果表明,1 mm高度内的涡发生器对起动状态的进气道主流流场结构、壁面压强分布影响不显著.但对于激波振荡流动,涡发生器会明显缩小外压缩面分离区运动范围,缩短振荡周期,提升振荡周期内壁面压强的时均值.涡发生器的影响程度随其高度的增大而增强,其中振荡周期从无涡发生器的4 ms缩短到1 mm高度涡发生器的3.13 ms.此外,0.5 mm高度涡发生器会使得进气道内部测点的压强振荡幅值整体下降,相比无涡发生器工况的下降幅度可达23%.流场结构与壁面压强信号的分析表明,涡流发生器主要通过其产生的流向涡影响激波振荡流动,包含流向涡对下游边界层的扰动以及流向涡与分离区的相互干扰.     

MULTI-OBJECTIVE OPTIMIZATION AND AERODYNAMIC PERFORMANCE ANALYSIS OF THE UPPER SURFACE FOR HYPERSONIC VEHICLES

为分析小攻角巡航条件下吸气式高超声速飞行器上壁面的变化对其气动性能和容积的影响, 以参数化后的飞行器上壁面对称面型线为设计变量, 在飞行马赫数6.5, 飞行高度27 km, 飞行攻角为4°的条件下, 采用计算流体力学为性能分析工具, Pareto多目标遗传算法为优化设计方法, 开展了二维条件下的升阻比/容积双目标优化设计. 在此基础上, 选择典型的二维优化结果, 重构生成对应的三维构型并进行数值分析, 获得了飞行器气动性能和容积间的相互关系. 结果表明在巡航条件下, 尽管二维/三维条件下飞行器的气动参数数值有较大差别, 但在这2种条件下, 飞行器的升阻比和容积间的关系均近似呈线性反比例关系. 同时, 对于三维构型而言, 在给定容积不变的条件下, 通过改变上壁面对称面型线的形状仅能使升阻比获得较小的增量(约0.36%). 相比之下, 当给定升阻比基本不变的条件下, 飞行器容积可调空间相对较大, 约为1.93%. 此外, 计算结果还表明, 在飞行器的容积基本不变情况下, 通过调节上壁面对称面型线, 可使飞行器的俯仰力矩获得5%左右的调节空间, 且其升阻比基本不变.     

微纳米结构壁面对流场及动压润滑的影响研究

固体边界具有的微纳米结构将影响流体在近壁面处的流动行为,进而由于尺度效应改变流体在整个微间隙的流动或润滑规律.将壁面可渗透微纳米结构等效为多孔介质薄膜,采用Brinkman方程来描述流体在近壁面边界渗透层内的流动,并将其与自由流动区域的不可压缩流体Navier-Stokes控制方程耦合,在界面处的连续边界条件下求解和分析了速度分布规律和压力变化规律.针对恒定法向承载力的油膜润滑条件,进一步讨论了静止表面或运动表面的微纳米结构对近壁面流动行为的影响;并揭示了考虑壁面微纳米结构的流体动压润滑的油膜厚度和摩擦系数的变化规律.论文结果为具有可渗透微结构表面的微间隙流动与润滑提供了理论参考.     

Prediction of hypersonic boundary layer transition on sharp wedge flow considering variable specific heat

When the air temperature reaches 600 K or higher, vibration is excited. The specific heat is not a constant but a function of temperature. Under this condition, the transition position of hypersonic sharp wedge boundary layer is predicted by using the improved eN method considering variable specific heat. The transition positions with different Mach numbers of oncoming flow, half wedge angles, and wall conditions are computed condition, the nearer to the Mach number The results show that for the same oncoming flow condition and wall transition positions of hypersonic sharp wedge boundary layer move much leading edge than those of the flat plate. The greater the oncoming flow the closer the transition position to the leading edge.     

Smooth leading edge transition in hypersonic flow

The boundary layer transition along the attachment line of a smooth swept circular cylinder in hypersonic flow is investigated in a blowdown wind tunnel. A wide range of spanwise Mach numbers Me (3.28 to 6.78) is covered with the help of different models at several sweep angles (60°?Λ?80°). The transition is indirectly detected by means of heat flux measurements. The influence of the wall to stagnation temperature ratio is investigated by cooling the model with liquid nitrogen.     

Development of slow disturbances in a hypersonic boundary layer

The mechanisms of development of slow time-dependent disturbances in the wall region of a hypersonic boundary layer are established and a diagram of the disturbed flow patterns is plotted; the corresponding nonlinear boundary value problem is formulated for each of these regimes. It is shown that the main factors that form the disturbed flow are the gas enthalpy near the body surface, the local viscous-inviscid interaction level, and the type, either subsonic or supersonic, of the boundary layer as a whole. Numerical and analytical solutions are obtained in the linear approximation. It is established that enhancement of the local viscous-inviscid interaction or an increased role for the main supersonic region of the boundary layer makes the disturbed flow by and large “supersonic”: the upstream propagation of the disturbances becomes weaker, while their downstream growth is amplified. Contrariwise, local viscous-inviscid interaction attenuation or an increased role for the main subsonic region of the boundary layer has the opposite effect. Surface cooling favors an increased effect of the main region of the boundary layer while heating favors an increased wall region effect. It is also found that in the regimes considered disturbances travel from the turbulent flow region downstream of the disturbed region under consideration counter to the oncoming flow, which may be of considerable significance in constructing the nonlinear stability theory.     

An experimental investigation on the surface water transport process over an airfoil by using a digital image projection technique

In the present study, an experimental investigation was conducted to characterize the transient behavior of the surface water film and rivulet flows driven by boundary layer airflows over a NACA0012 airfoil in order to elucidate underlying physics of the important micro-physical processes pertinent to aircraft icing phenomena. A digital image projection (DIP) technique was developed to quantitatively measure the film thickness distribution of the surface water film/rivulet flows over the airfoil at different test conditions. The time-resolved DIP measurements reveal that micro-sized water droplets carried by the oncoming airflow impinged onto the airfoil surface, mainly in the region near the airfoil leading edge. After impingement, the water droplets formed thin water film that runs back over the airfoil surface, driven by the boundary layer airflow. As the water film advanced downstream, the contact line was found to bugle locally and developed into isolated water rivulets further downstream. The front lobes of the rivulets quickly advanced along the airfoil and then shed from the airfoil trailing edge, resulting in isolated water transport channels over the airfoil surface. The water channels were responsible for transporting the water mass impinging at the airfoil leading edge. Additionally, the transition location of the surface water transport process from film flows to rivulet flows was found to occur further upstream with increasing velocity of the oncoming airflow. The thickness of the water film/rivulet flows was found to increase monotonically with the increasing distance away from the airfoil leading edge. The runback velocity of the water rivulets was found to increase rapidly with the increasing airflow velocity, while the rivulet width and the gap between the neighboring rivulets decreased as the airflow velocity increased.     

Numerical simulations of droplet train and free surface jet impingement

The cooling behavior of the impingement of a droplet train, and free surface jets over a heated and pre-wetted surface is explored employing an Algebraic Volume-of-Fluid methodology. The code is based on a modified version of the two-phase numerical solver interFoam (OpenFOAM) (Trujillo and Lewis, 2012). Two versions of the free surface jet are studied. The first consists of a fully-developed profile exiting the nozzle, and the second is characterized by a uniform velocity distribution. Results show that both jet configurations have higher cooling performance than the droplet train locally and globally, with the fully-developed case being the most effective of the two jet arrangements. Locally, the performance is measured by radial profiles of the boundary-layer-displacement thickness and heat transfer coefficient. Globally, the cooling effectiveness is directly proportional to the surface area that resides within the high-convection region, i.e. before the boundary layer separation point. On a temporal basis, the liquid film within the impingement region of the droplet train exhibits pronounced variations in velocity magnitude and film thickness. This is directly attributed to the nature of continuous droplet impacts affecting the impingement region, and gives rise to an unsteady cooling and heating of the fluid near the wall. In contrast for the jets, the film and the corresponding free surface are nearly steady with only minor perturbations.     

Mechanism of unsteady aerodynamic heating with sudden change in surface temperature

The characteristics and mechanism of unsteady aerodynamic heating of a transient hypersonic boundary layer caused by a sudden change in surface temperature are studied. The complete time history of wall heat flux is presented with both analytical and numerical approaches. With the analytical method, the unsteady compressible boundary layer equation is solved. In the neighborhood of the initial and final steady states, the transient responses can be expressed with a steady-state solution plus a perturbation series. By combining these two solutions, a complete solution in the entire time domain is achieved. In the region in which the analytical approach is applicable, numerical results are in good agreement with the analytical results, showing reliability of the methods. The result shows two distinct features of the unsteady response. In a short period just after a sudden increase in the wall temperature, the direction of the wall heat flux is reverted, and a new inflexion near the wall occurs in the profile of the thermal boundary layer. This is a typical unsteady characteristic. However, these unsteady responses only exist in a very short period in hypersonic flows, meaning that, in a long-term aerodynamic heating process considering only unsteady surface temperature, the unsteady characteristics of the flow can be ignored, and the traditional quasi-steady aerodynamic heating prediction methods are still valid.     

Film cooling by injection into a turbulent boundary layer

Most papers on film cooling concern injection of a homogeneous gas. Stollery et al. [1] examined the case of tangential injection of gas into a boundary layer, the specific heat_63-01 differing little from that of the main flow,_63-02.Here we examine the effectiveness of film cooling of a thermally isolated planar wall by local supply to a turbulent boundary layer.     

Hybrid LES-RANS study of an effusion cooling array with circular holes

A multi-row effusion cooling configuration with scaled gas turbine combustor conditions is studied numerically, using a novel wall-proximity-based hybrid LES-RANS approach. The distribution of the coolant film is examined by surface adiabatic cooling effectiveness (ACE). Simulation results have shown that the accuracy of cooling effectiveness prediction is closely related to the resolution of turbulent flow structures involved in hot-cold flow mixing, especially those close to the plate surface. The formation of the coolant film in the streamwise direction is investigated. It is shown that the plate surface directly downstream the coolant holes are covered well by the coolant jets, while surface regions in between the two columns of the coolant holes could not be protected until the coolant film is developed sufficiently in the spanwise direction in the downstream region. More detailed study has also been carried out to study the time-averaged and time-dependent flow fields. The relation between the turbulent flow structures and coolant film distribution are also examined. The Kelvin–Helmholtz instability in the upper and lower coolant jet shear layer, is found to have the same frequency of around 8000 Hz, and is independent of the coolant hole position. Additionally, it is suggested by the spectral coherence analysis that those unsteady flow structures from the lower shear layer are closely related to the near wall flow temperature, and such effect is also independent of the coolant hole position.