高超声速条件下7°直圆锥边界层转捩实验研究

在Ma=6低噪声风洞中开展了半锥角7?的直圆锥边界层转捩相关实验研究.利用响应频率达到MHz量级的高频压力传感器对圆锥壁面脉动压力进行了测量,研究了高超声速圆锥边界层中扰动波的发展过程.结果表明:高超声速圆锥边界层中第二模态扰动波产生的位置以及扰动波特征频率和波长等参数受雷诺数影响较大,当单位雷诺数从2×106m~(-1)增加到8×106m~(-1)时,第二模态波的特征频率从55 k Hz增加到226 k Hz;随着单位雷诺数增加,边界层中扰动增长速度加快,第二模态波出现在圆锥表面更靠近上游的位置;相同单位雷诺数条件下,随着第二模态波的向下游传播,其特征频率逐渐减小.通过对比发现自由来流湍流度对边界层中扰动波的发展同样有较大影响,自由来流湍流度降低,边界层中的第二模态波的特征频率明显减小.利用互相关分析得出第二模态扰动波在边界层中的传播速度大约为当地主流速度的0.8—0.9倍.在1?小攻角条件下,圆锥迎风面和背风面边界层发展呈现出明显的差异,背风面边界层中扰动发展提前,第二模态波出现在更靠近上游的位置,而迎风面中扰动发展受到抑制,第二模态波特征频率更大.     

基于温敏漆技术的圆锥高超声速大攻角绕流背风面流动结构实验研究

高超声速流动中, 大攻角下圆锥背风面边界层会存在流动分离与再附、边界层转捩等多种流动现象, 进而对圆锥表面温度分布产生显著的影响。为了对这一复杂流动规律及其对表面温升分布的影响进行讨论, 研究基于温敏漆技术, 得到了在Mach数为6的低湍流度来流条件下, 攻角为10°的圆锥背风面温升分布结果。通过对不同位置、不同方位角处温升分布曲线的分析, 对大攻角下圆锥背风面边界层流动发展过程及不同发展阶段的流动特征进行了讨论。同时, 通过对来流总压的调节, 得到了不同Reynolds数下的圆锥背风面温升分布结果, 总结了Reynolds数对流动的影响规律。研究发现, 高超声速大攻角圆锥背风面边界层流动发展过程中会依次出现层流分离、定常横流涡影响、转捩以及湍流分离与再附等流动特征, 而在不同的Reynolds数下, 各个流动特征产生影响的范围不同, 随着Reynolds数的降低, 层流范围和定常横流涡影响范围均有所增加, 而从观察到横流影响到转捩开始发生的范围基本相同。      

压缩拐角流动中G?rtler涡特性及热流分布实验研究

对高超声速压缩拐角流动中G?rtler涡特性及热流分布进行了实验研究.开发了温敏漆（temperature sensitive paint,TSP）系统,简要介绍了TSP技术的原理、文章所用的TSP涂料的标定曲线、辅助设备参数、实验过程数据后处理过程,采用基于离散Fourier定律的热流算法.研究在Ma=6低噪声风洞中进行,采用TSP技术,得到压缩拐角斜坡板上的热流分布图像,并对高低热流条带现象做出解释,与G?rtler涡有对应关系.通过改变拐角角度及来流参数,获得了不同拐角和单位Reynolds数条件下的热流分布图像,分析得到压缩拐角斜坡上G?rtler涡特性及热流分布在变参数条件下的变化规律.研究发现:当增加拐角角度或增大单位Reynolds数时,G?rtler涡的波长减小,且涡的起始位置更靠近拐角;随单位Reynolds数增加,斜坡上热流值整体增加,热流峰值位置前移;峰值位置后,热流缓慢减小的区域与G?rtler涡位置相对应.      

Experimental study of mean and pulsation flow characteristics in the 2D/3D supersonic boundary layer behind flat roughness elements

For the first time, detailed experiments aimed at evaluation of the influence of flat microroughness elements (adhesive tape stickers attached to the streamlined surface) on the mean and pulsation flow characteristics in flat-plate and swept-wing boundary layers were carried out under comparable dimensional conditions. At Mach number M = 2, the sticker shape was found to have almost no influence on the amplitude of mass-flow pulsations in the wake of sticker. The relative receptivity of 3D boundary layer to stationary disturbances was obtained to be approximately 1.5 times higher than the receptivity of 2D boundary layer. It was found that the presence of a sticker at the line of the source of controlled disturbance changes the pulsation spectra in 2D and 3D boundary layers in comparison with unperturbed flow. In the wake behind the center of a sticker, the flow can be most unstable.     

Direct numerical simulation of a supersonic turbulent boundary layer on a flat plate and its analysis

In the past three decades, considerable progress has been made in the investigation of incompressible turbulent boundary layer through experiments, DNS and theoretical works, including: (1) the statistics characteristic and structure of turbulence; (2) the co-herent structures in turbulent flows; (3) turbulence modeling and the large eddy simula-tion (LES). In contrast, the progress was very slow for the compressible, in particular, the super-sonic turbulent boundary layer. Recent works on d…     

Experimental and numerical investigation of the recovery ratio of a wedge-shaped hot-film probe

The recovery ratio of a wedge-shaped hot-film probe was determined in an experimental as well as numerical study, since this information is still unpublished and essential for using the probe in hot-film anemometry. The experiments were conducted at the Khristianovich Institute of Theoretical and Applied Mechanics (ITAM) in Novosibirsk, Russia, and the simulations were performed with StarCCM+, a commercial 2nd order finite volume code. In the analysis, the Mach number was varied between M = 2 and M = 4, and the unit Reynolds number ranged from Re₁ = 3.8?10⁶ to Re₁ = 26.1?10⁶ m^?1, depending on the Mach number. During the experiment, the stagnation temperature was kept constant for each Mach number at a separate value in the range of T ₀ = 289 ± 7 K. Three different stagnation temperatures were used in the simulations: T ₀ = 259 K, T ₀ = 289 K, and T ₀ = 319 K. The difference between the experimental and the numerical results is ≤ 0.5 %, and, therefore, both are in very good accordance. The influence of the Mach number, of the unit Reynolds number, and of the stagnation temperature was analysed, and three different fitting functions for the recovery ratio were established. In general, the recovery ratio shows small variations with all three tested parameters. These dependencies are of the same order of magnitude.     

Influence of blunt-cone tip temperature on the laminar-turbulent transition in hypersonic boundary layers

In the present paper, we report on the results of a combined experimetal and numerical study of the laminar-turbulent transition at Mach number 5.95 in the boundary layers of 7-deg cones with a small nose-tip bluntness radius (down to 1.5 mm). The tip temperature of the model was varied in the range from 90 to 440 K. It was confirmed that the small nose-tip bluntness substantially shifts the position of the transition in downstream direction. This effect was also retained in testing the models with local heating/cooling of the cone tip. It is shown that for the experimental conditions implemented in the Transit-M wind tunnel, ITAM SB RAS, the local heating of blunt cone tips exerts almost no influence on the position of the transition. The local cooling of the cone tip with R = 1.5 mm leads to a shift of the transition position in upstream direction.     

Effect of the surface roughness of blunt cone forebody on the position of laminar-turbulent transition

In the present paper, data on the effect of the surface roughness of blunt cone forebody on the position of laminar-turbulent transition are reported. The study was carried out under freestream Mach 5.95. It was found that the roughness position plays a substantial role in the transition process. Critical Reynolds numbers at which the laminar-turbulent transition occurs on the nose-tip of the model were identified. For the first time, hysteresis in transition position was observed.     

基于转捩SST模型凸起圆柱绕流数值研究

为了研究局部凸起对边界层转捩的影响,采用转捩SST模型分别对亚临界、临界和超临界状态下带突起的圆柱绕流问题进行了数值模拟,分析了不同Reynolds数下带突起的圆柱绕流问题的近壁面流动特征以及表面时均压力与摩擦力系数的分布和凸起对圆柱表面流动分离以及转捩的影响,对比了有无凸起两侧圆柱表面时均压力、摩擦力系数的不同. 结果表明:当来流Reynolds数处于临界区时,气流在圆柱上表面凸起处形成了3个反向旋转的漩涡,之后随着θ的增大,发生了流动分离和流动转捩现象;对于不同Reynolds数下的来流,圆柱上表面的凸起可以使气流发生转捩的位置提前;圆柱上表面的凸起使流速增大、压强降低,从而导致圆柱产生升力,随着来流Reynolds数的增大,其升力逐渐变大.      

Effects of mean flow convection,quadrupole sources and vortex shedding on airfoil overall sound pressure level

This paper presents a further analysis of results of airfoil self-noise prediction obtained in the previous work using large eddy simulation and acoustic analogy. The physical mechanisms responsible for airfoil noise generation in the aerodynamic flows analyzed are a combination of turbulent and laminar boundary layers, as well as vortex shedding (VS) originated due to trailing edge bluntness. The primary interest here consists of evaluating the effects of mean flow convection, quadrupole sources and vortex shedding tonal noise on the overall sound pressure level (OASPL) of a NACA0012 airfoil at low and moderate freestream Mach numbers. The overall sound pressure level is the measured quantity which eventually would be the main concern in terms of noise generation for aircraft and wind energy companies, and regulating agencies. The Reynolds number based on the airfoil chord is fixed at _Rec=408,000${\mathit{Re}}_c=408,000$ for all flow configurations studied. The results demonstrate that, for moderate Mach numbers, mean flow effects and quadrupole sources considerably increase OASPL and, therefore, should be taken into account in the acoustic prediction. For a low Mach number flow with vortex shedding, it is observed that OASPL is higher when laminar boundary layer separation is the VS driving mechanism compared to trailing edge bluntness.     

An experimental study of the natural noise in the Transit-M hypersonic wind tunnel

Tests were performed at Mach number M = 6. Observations of flow pulsations were performed with the help of high-frequency measurements of pressure and heat fluxes and with the help of hot-wire anemometry. Measurements were carried out in a broad range of Re₁ numbers at various orientations and locations of sensors with respect to the nozzle exit plane. It is shown that the Transit-M noise varies weakly over the longitudinal coordinate, whereas on moving along the radial coordinate, some increase of the fluctuation level is observed due to approaching the mixing layer. Also, the distribution of axial asymmetry of the noise level was obtained.     

Influence of viscosity on aerodynamic sound emission in free space

The theory of aerodynamic sound is reformulated with account taken of the influence of viscosity in the source flow on the sound emission in free space. This is based on the Ffowcs Williams form of the Lighthill equation. The source flow is assumed to be determined by a localized vorticity field and characterized by low Mach numbers and high Reynolds numbers. It is found that the acoustic pressure emitted by the viscous vortex motion is composed of a quadrupole and two kinds of monopole. The Reynolds stress in the source flow is decomposed into isotropic and non-isotropic parts. The non-isotropic part leads to the quadrupole wave derived by Möhring, while the isotropic part is related to one of the monopole-like waves radiated when the total kinetic energy changes. The other monopole wave is associated with entropy production by viscous dissipation of the kinetic energy. All three components are influenced by viscosity.     

强梯度复杂流场中的粒子动力学响应试验研究

示踪粒子在(高)超声速流场中的动力学响应是粒子成像测速等粒子示踪测量技术的关键问题之一.现有文献对粒子动力学响应的试验测量往往是通过单个斜激波响应的测量方法.然而,当示踪粒子用于测量高速飞行器发动机内部复杂的激波串流场时,粒子将经历由多道激波导致的速度、压力、黏性等剧烈变化.本文结合目前(高)超声速飞行器的研究热潮,重点关注示踪粒子在应用于发动机内部具有连续激波的复杂流场测量中存在的跟随性评估方面,开展了一系列的相关试验研究.包括测量超声速风洞的喷管出口速度分布以验证测试系统的性能,在马赫4.2和3.0流场中测量了粒子对二维10°和15°单斜劈绕流中的斜激波动力响应,并测量了模拟发动机内部连续梯度的双斜劈粒子斜激波动力响应.结合粒子动力学的理论模型,得到了各状态的粒子弛豫时间、弛豫距离、Stokes数.基于图像方法、统计学规律分析了激波非定常抖动对测量结果的影响,并对测量结果进行了修正.结果显示,相同斜劈角度下,马赫数越高,粒子的弛豫时间、弛豫距离就越大.但是在相同的来流马赫数下,斜劈角度越大,粒子的弛豫时间、弛豫距离反而减小.在强梯度之后由于流场的雷诺数和黏性系数变化剧烈,粒子的跟随性降低了大约5.7%,stokes数增加了约1%.虽然在本文条件下Stokes数仍满足超声速流场对粒子跟随性的要求,但粒子响应的降低无疑是值得关注的,尤其是当其被应用于具有更多连续梯度的复杂流场测量中.     

基于深度强化学习的近壁圆柱绕流控制

基于深度强化学习方法（deep reinforcement learning,DRL）,采用一对施加在圆柱表面的质量流量为零的射流,对Re=200,400,间隙比G/D=0.5,0.7,1.0,1.5,2.0的近壁圆柱进行主动流动控制研究。通过DRL方法获得不同参数下的射流控制策略与相应控制效果,并讨论了不同射流位置（90°,270°）、（90°,320°）、（90°,360°）对控制效果的影响。研究发现Reynolds数、间隙比和射流位置都对控制效果有重要的影响。射流位置为（90°,270°）时通过DRL控制可以有效降低阻力系数及其波动,控制后的圆柱尾迹被拉长且圆柱前后压差降低。射流位置为（90°,320°）和（90°,360°）的控制效果相似,都能使平均阻力系数有所降低,但控制位置的不对称性导致控制后的阻力系数波动较大。Reynolds数和间隙比的增大会增加控制射流的质量流量水平,在相同条件下,使用（90°,270°）的射流位置可以用相同的质量流量得到更好的控制效果。      

激光聚焦扰动作用下高超声速边界层稳定性实验研究

在马赫数6、单位雷诺数3.1×10⁶/m的条件下对半锥角7°直圆锥边界层稳定性开展了实验研究.以激光聚焦于流场中局部空间而产生的膨胀冲击波作为人工添加的小扰动,分析了该扰动对高超声速圆锥边界层流动稳定性的影响.实验中利用响应频率达到兆赫兹量级的高频压力传感器对圆锥壁面脉动压力进行测量,通过对压力数据进行短时傅里叶分析和功率谱分析发现,相比于不添加激光聚焦扰动的结果,添加激光聚焦扰动使边界层中第二模态波的出现位置提前,且扰动波的幅值大幅度地增加,在相同的流向范围内,激光聚焦扰动将边界层中的扰动波从线性发展阶段推进到非线性发展阶段,其对边界层中扰动波发展的促进效果明显.同时,激光聚焦位置的不同对边界层中扰动波的发展也具有不同的影响.当激光直接聚焦于圆锥壁面X=100 mm位置时,边界层中频率为90 kHz的扰动波幅值增长最快,在X=500 mm的位置处其幅值放大倍数为3.81,相比而言当激光聚焦位置位于圆锥前方自由来流中时,边界层幅值增长最快的扰动波频率大幅减小为73 kHz,相同范围内,其幅值放大倍数为4.51倍.由此可见,当激光聚焦位置位于圆锥上游的自由来流中时,其对边界层中扰动波的影响更为显著.     

Auto-oscillations in supersonic boundary layer

The production of periodic oscillations in a supersonic boundary layer at the moderate and high Mach numbers (M = 2 and 5.35) is investigated within the framework of the weakly nonlinear stability theory of the second order in nonlinearity. The model includes the effects of self-action, such as the generation of stationary secondary harmonics and the disturbances of double frequencies. It is shown that for two-dimensional vortex disturbances, the character of the excitation of vortex disturbances changes from the mild one to the stiff one with the increasing Mach number, which leads to a reduction of the critical Reynolds number Re_c. For three-dimensional disturbances of low azimuthal wave numbers, a supercritical auto-oscillatory regime sets in. A complex regime realizes for two-dimensional acoustic disturbances at M = 5.35 with a stiff excitation in the region of Re_c.     

Quantitative visualization of the leading-edge vortices on a delta wing by using pressure-sensitive paint

Leading-edge vortices on a simple delta wing were visualized by using pressure-sensitive paint (PSP). PSP is an optical pressure measurement technique based on oxygen quenching of luminescent molecules. In the present study, we used PSP composed of platinum octaethylporphyrine (PtOEP) and fluoropolymer (poly-IBM-co-TFEM [Poly (isobutylmethacrtlate-co-trifluoroethylate)]). This new paint has higher sensitivity to pressure and lower sensitivity to temperature than previous ones, reducing an error due to temperature variation during a wind tunnel test. A thin coating of PSP was applied to a delta wing model with 70-degree leading-edge sweep. The coating was excited by Xenon light and emission from the coating was detected by a high-resolution CCD camera. Tests were done at subsonic speeds in the 0.2-m Supersonic Wind Tunnel at the National Aerospace Laboratory in Japan. Complicated flow structures on the delta wing including primary and secondary vortices were clearly visualized using pressure-sensitive paint. An a priori calibration technique was used to convert measured luminescent intensity into pressure. The obtained pressure distributions were in good agreement with pressure tap data. Pressure maps were obtained for various Mach numbers, Reynolds numbers and angles of attack. It was found that an increase in Mach number delayed vortex breakdown while Reynolds number had little effect on the vortex formation.     

基于线性抛物化稳定性方程的后掠翼边界层内横流稳定性研究

采用有限差分法求解考虑模型曲率的线性抛物化稳定性方程(LPSE),分析无限展长后掠翼边界层内的横流驻波不稳定,并与实验结果进行对比,研究LPSE方法的模拟效果及其适用范围.研究表明,在横流驻波扰动增长的初期,LPSE能够准确的预测扰动的e^N曲线,较好地描述边界层内的流动结构和扰动形态;当扰动增长到足够大时,扰动的高阶项不能再被忽略,LPSE的线性假设不再成立,需要采用非线性的方法(NPSE)来分析该状态.计算分析发现,模型曲率和边界层非平行性对后掠翼边界层内横流驻波的稳定性分析影响很大,影响程度与雷诺数无关.对于本文研究的模型,曲率对边界层内的横流扰动起着稳定的作用,而非平行性对扰动起不稳定的影响.     

Mixing in a T-shaped micromixer at moderate Reynolds numbers

In the present work, the regimes of the flow and mixing of fluids in a T-shaped micromixer in the range of the Reynolds numbers from 1 to 1000 are investigated systematically with the aid of numerical modeling. The flow and mixing regimes are shown to alter substantially with increasing Reynolds numbers. Five different flow regimes have been identified in the total. The dependencies of the friction coefficient and mixing efficiency on the Reynolds number are obtained. A sharp increase in the mixing efficiency at a flow transition from the symmetric to asymmetric steady regime is shown. On the other hand, the mixing efficiency slightly drops in the laminar-turbulent transition region. A substantial influence of the slip presence on walls on flow structure in the channel and mixing efficiency has been revealed.     

DNS of compressible turbulent boundary layer around a sharp cone

Direct numerical simulation of the turbulent boundary layer over a sharp cone with 20° cone angle (or 10° half-cone angle) is performed by using the mixed seventh-order up-wind biased finite difference scheme and sixth-order central difference scheme. The free stream Mach number is 0.7 and free stream unit Reynolds number is 250000/inch. The characteristics of transition and turbulence of the sharp cone boundary layer are compared with those of the flat plate boundary layer. Statistics of fully developed turbulent flow agree well with the experimental and theoretical data for the turbulent flat-plate boundary layer flow. The near wall streak-like structure is shown and the average space between streaks (normalized by the local wall unit) keeps approximately invariable at different streamwise locations. The turbulent energy equation in the cylindrical coordinate is given and turbulent energy budget is studied. The computed results show that the effect of circumferential curvature on turbulence characteristics is not obvious.