前缘曲率对三维边界层内被激发出非定常横流模态的影响研究

三维边界层感受性问题是三维边界层层流向湍流转捩的初始阶段,是实现三维边界层转捩预测与控制的关键环节.在高湍流度的环境下,非定常横流模态的失稳是导致三维边界层流动转捩的主要原因;但是,前缘曲率对三维边界层感受性机制作用的研究也是十分重要的课题之一.因此,本文采用直接数值模拟方法研究在自由来流湍流作用下具有不同椭圆形前缘三维（后掠翼平板）边界层内被激发出非定常横流模态的感受性机制;揭示不同椭圆形前缘曲率对三维边界层内被激发出非定常横流模态的扰动波波包传播速度、传播方向、分布规律、感受性系数以及分别提取获得一组扰动波的幅值、色散关系和增长率等关键因素的影响;建立在不同椭圆形前缘曲率情况下,三维边界层内被激发出非定常横流模态的感受性问题与自由来流湍流的强度和运动方向变化之间的内在联系;详细分析了不同强度各向异性的自由来流湍流在激发三维边界层感受性机制的物理过程中起着何种作用等.通过上述研究将有益于拓展和完善流动稳定性理论,为三维边界层内层流向湍流转捩的预测与控制提供依据.     

三维边界层内定常横流涡的感受性研究

层流向湍流转捩的预测与控制一直是研究的前沿热点问题之一,其中感受性阶段是转捩过程中的初始阶段,它决定着湍流产生或形成的物理过程.但是有关三维边界层内感受性问题的数值和理论研究都比较少;实际工程问题中大部分转捩过程都是发生在三维边界层流中,所以研究三维边界层中的感受性问题显得尤为重要.本文以典型的后掠角45?无限长平板为例,数值研究了在三维壁面局部粗糙作用下的三维边界层感受性问题,探讨了三维边界层感受性问题与三维壁面局部粗糙长、宽和高之间的关系;然后,考虑在后掠平板上设计不同的三维壁面局部粗糙的分布状态、几何形状、距离后掠平板前缘的位置以及流向和展向设计多个三维壁面局部粗糙对三维边界层感受性问题有何影响;最后,讨论两两三维壁面局部粗糙中心点之间的距离以及后掠角的改变对三维边界层感受性的物理过程将会发生何种影响等.这一问题的深入研究将为三维边界层流中层流向湍流转捩过程的认识和理解提供理论依据.     

壁面局部吹吸边界层感受性的数值研究

目前理论、实验以及数值模拟主要研究自由来流中的小扰动与壁面局部粗糙相互作用激发边界层感受性问题. 但是, 针对自由来流湍流与壁面局部吹吸相互作用诱导边界层感受性的相关报道甚少. 本文采用直接数值模拟和快速傅里叶变换的方法, 数值研究了二维平板壁面具有局部吹吸的边界层感受性问题. 结果发现, 在二维边界层内能找到一组被激发产生的Tollmien-Schlichting(T-S)波波包的包络序列以及从波包中能够分离出一组稳定的、中性的和不稳定的T-S波, 证明了二维边界层内感受性现象的存在性. 经数值计算获得了T-S波波包传播的群速度; 并建立了自由来流湍流强度、壁面局部吹吸强度和长度与二维边界层感受性之间的关系, 获得了与Dietz感受性实验相类似的结论. 另外, 还发现在自由来流湍流与壁面局部吹、吸相互作用下能诱导二维边界层内产生相位相反的T-S波. 依据这一理论机理来优化设计局部吹吸装置, 不但能促使层流向湍流转捩的提前, 也可以延迟转捩过程的发生, 达到控制湍流运动的目的.     

无限薄平板边界层前缘感受性过程的数值研究

在边界层流中层流向湍流转捩机理的研究一直是人们所关注的重要理论课题之一.感受性阶段是边界层内整个转捩过程中的初始阶段,它在层流向湍流转捩过程中起着关键性的作用.但是,在过去的边界层前缘感受性研究中大多数都是针对外部声波小扰动,很少见到考虑在自由流中普遍存在的自由来流湍流作用下边界层内诱导前缘感受性问题的相关报道.本文采用直接数值模拟的方法,研究自由来流湍流与无限薄平板前缘驻点扰动作用下边界层流中前缘感受性过程的内在机理.数值结果发现,在自由来流湍流与无限薄平板前缘驻点扰动作用下边界层流中能被感受出一组小扰动波,且它们的色散关系、增长率、中性曲线等结果都与流动稳定性中的线性理论获得的理论解相一致,由此可知在边界层内被激发产生的一组小扰动波就是Tollmien-Schlichting波,这也证明自由来流湍流与无限薄平板前缘驻点扰动相互作用是激发边界层流中前缘感受性过程的另一种物理机理;另外,还探讨了自由来流湍流度以及自由来流湍流的运动方向对无限薄平板边界层前缘感受性过程有何影响等.总之,开展边界层前缘感受性过程的深入研究,有益于完善流动稳定性理论,将为层流向湍流转捩过程的预测提供合理的理论依据.     

高超声速椭圆锥边界层横流转捩特性大涡模拟

横流效应显著影响高超声速飞行器的三维边界层转捩过程, 深化对该流动机制的认识有助于提升和改善飞行器气动性能及热力学环境. 针对HIFiRE5椭圆锥绕流问题, 采用大涡模拟方法计算分析了超声速边界层横流转捩特性, 并揭示其中的流动机理. 参考HIFiRE5风洞模型试验条件, 数值模拟中椭圆锥来流入口处施加人工速度扰动以激发边界层内不稳定扰动波, 进而预测了高超声速边界层流动横流失稳、转捩过程等基本流动特征, 并基于转捩热流分布形态对比, 获得了与试验数据基本吻合的计算结果. 研究发现, 椭圆锥中心线流动汇聚形成的流向涡结构非常容易失稳, 另外在中心线及侧缘之间的中部区域存在较强的横流不稳定性, 两种机制共同作用影响边界层转捩过程. 此外, 分析了来流扰动幅值对边界层横流失稳转捩的影响, 并发现静来流条件下, 横流区域出现两组独立的定常横流涡结构, 而强噪声来流条件下, 中心线主涡和中部横流涡均发生失稳转捩, 且在椭圆锥表面形成多峰状的转捩阵面. 最后, 深入分析流场的压力脉动动力学特性, 揭示了三维边界层发生失稳转捩的非线性演化机制.      

前缘曲率变化对平板边界层感受性问题的影响

边界层感受性问题是层流向湍流转捩的初始阶段,是实现边界层转捩预测和控制的关键环节.目前已有的研究成果显示,在声波扰动或涡波扰动作用下前缘曲率变化对边界层感受性机制有着显著的影响.本文采用直接数值模拟方法,研究了在自由来流湍流作用下具有不同椭圆形前缘平板边界层感受性问题,揭示椭圆形前缘曲率变化对平板边界层内被激发出Tollmien-Schlichting (T-S)波波包的感受性机制以及波包向前传播群速度的影响;通过快速傅里叶分析方法从波包中提取获得了不同频率的T-S波,详细分析了前缘曲率变化对不同频率的T-S波的幅值、色散关系、增长率、相速度以及形状函数的作用;确定了前缘曲率在平板边界层内激发T-S波的感受性过程中所占据的地位.通过上述研究能够进一步认识和理解边界层感受性机制,从而丰富和完善了流动稳定性理论.     

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

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

展向凹槽及泄流孔对高超声速平板边界层转捩影响的试验研究

针对展向凹槽和泄流孔对高超声速钝平板边界层转捩的影响,在中国空气动力研究与发展中心F2 m激波风洞(FD-14A)开展了试验及初步的计算与理论研究.试验的来流马赫数为6、单位雷诺数为3.3×107/m,平板的前缘半径为1 mm,攻角为–4°.在距平板前缘110 mm处布置三组不同的二维展向凹槽,凹槽的宽度与深度分别为凹槽1(2.5 mm,1 mm)、凹槽2(3.75 mm,1.5 mm)、凹槽3(5 mm,2 mm),同时凹槽1的两端可以打开泄流孔,记为凹槽4,不含凹槽时的光滑平板情况记为凹槽5或平板.采用热流传感器测量了不同情况下平板中心线的热流分布,测量结果显示,光滑平板情况在x≈340 mm处开始转捩,在x≈425 mm处转捩接近完成.凹槽导致平板边界层的转捩位置提前,且随着凹槽宽度及深度的增加,对转捩的促进作用增强,转捩位置向上游移动.凹槽1增加泄流孔后(凹槽4)其热流分布及转捩位置与光滑平板情况基本一致.边界层流动完全转捩为湍流后,各情况下的热流差别较小,表明不同规格的凹槽只影响转捩过程中的热流分布,对转捩完成后的湍流壁面热流影响较小.数值计算(CFD)结果显示,泄流孔导致了被动抽吸,试验结果显示凹槽两端的泄流孔抽吸效应抵消了凹槽对平板中心线边界层转捩的促进作用.采用线性稳定性理论(LST)及最优扰动方法分析了光滑钝平板情况的流动失稳机制.LST结果显示,本文平板流动不存在Mack第一模态、第二模态失稳,因此传统的模态失稳机制无法解释试验中观测到的转捩现象.最优扰动计算显示,平板流动存在较强的非模态失稳,可以定性解释观测到的转捩现象.     

合成冷/热射流控制超声速边界层流动稳定性

为了探究超声速边界层流动稳定性及其转捩控制机理,提出基于合成冷/热射流的边界层速度-温度耦合控制方法,并通过数值模拟研究了Ma=4.5超声速平板边界层不稳定波的传播,采用线性稳定性理论中的时间模式分析了壁面吹吸、射流温度、扰动频率、扰动振幅等对不稳定波控制效果的影响.结果表明:无射流控制时,边界层内同时存在不稳定的第一模态扰动波和第二模态扰动波,且二维波形式的第二模态占主导地位;壁面吹吸作用下,仅出现更加不稳定的第二模态,第一模态被抑制;速度-温度耦合控制下,射流温度对扰动模态的不稳定区域大小及扰动增长率影响显著,射流温度与来流温度不同时,温度的脉动使得流动转捩为湍流的速度加快,边界层速度型更加饱满,抗干扰能力增强,流动稳定性提高;高频的吹吸扰动对流场的控制效果优于低频扰动,扰动频率超过400 Hz时,第二模态扰动波时间增长率降低,扰动分量对边界层速度剖面和温度剖面的修正加快,第二模态更加稳定;扰动振幅减小为主流速度的1%时,仅出现时间增长率较小的第二模态,控制效果较好,进一步减小时,第一模态重新出现,并且波数范围与第二模态先重合后分离,对应的时间增长率先增加后减小.研究结果为边界层转捩控制技术提供了新的思路.     

高超声速飞行器复杂外形转捩预测

采用k-ω-γ 转捩模式对某新型飞行器外形的典型流动特征和边界层失稳特性进行了分析.研究结果表明,横流是影响飞行器大面积转捩的主要因素.随着高度增加,来流Reynolds数减小,迎风面和背风面的转捩起始位置均向下游移动.随着攻角增加,头部附近背风面的展向压力梯度增大,横流效应增强,转捩起始位置向上游移动;另一方面攻角增...     

Distributed two-dimensional boundary-layer receptivity to non-stationary vortical disturbances in the presence of surface roughness

The distributed (over the longitudinal coordinate) excitation of two-dimensional (2D) Tollmien — Schlichting (TS) waves by weak non-stationary free-stream vortices propagating along the edge of the laminar boundary layer developing over a surface with small-amplitude 2D roughness is examined. The vorticity vector of the free-stream vortices was oriented over the span of the model, i. e., did not depend on the transverse coordinate. The theoretical analysis of the excitation mechanism reported in [1] was refined to develop, around it, a procedure making it possible to experimentally determine the coefficients of distributed vortical receptivity of the flow and the coefficients of “vortex-roughness” receptivity by fitting the experimental distributions with analytical solutions. Under conditions with controllable excitation of disturbances, a detailed hot-wire study of free-stream disturbances and boundary-layer disturbances at several vortex frequencies and at several surface-roughness periods was performed, and the shape of the controllable surface roughness was measured. The point-source method was employed to experimentally examine the characteristics of linear three-dimensional (3D) stability of the flow to TS waves, necessary for determination of the coefficients of distributed receptivity. It was found that the free-stream vortices with transverse orientation of the vorticity vector excited boundary-layer TS waves via two receptivity mechanisms: (a) on the smooth surface (due to natural non-uniformity of the flow) and (b) during interaction of the vortices with the surface roughness. The developed approach was used to experimentally estimate the amplitudes and phases of the coefficients of both types of distributed vortical receptivity as dependent on problem parameters. The absolute values of both types of receptivity coefficients were found to rapidly grow in value with increasing vortex frequency. It is shown that the most efficient excitation of TS waves is observed in the situation with satisfied resonance conditions for streamwise vortex, surface-roughness, and TS-wave streamwise wavenumbers, resulting in strong deviation of the increments of the TS waves from the linear-stability increments. Under no-resonance conditions, only amplitude beats of boundary-layer disturbances were observed. This work was financialy supported by the Russian Foundation for Basic Research (Grant No. 03-01-00299).     

Generation of nonstationary Görtler vortices by localized surface nonuniformities. Receptivity coefficients

The mechanism of production of nonstationary Görtler vortices in a boundary layer on concave wall by surface nonuniformities (vibrations and roughness) has been experimentally examined. The nonuniformities were produced by a specially developed disturbance source. They were controlled, localized along the streamwise coordinate, and periodic over the span of the experimental model. Tests in a low-turbulence wind tunnel have proved that the disturbance source is an efficient means of experimental study of the receptivity and stability problem for boundary layers dominated by Görtler instability. The operation of the disturbance source leads to the production of small-amplitude nonstationary Görtler vortices (tenth or hundredth fractions of a per cent of the free-stream velocity) with predefined characteristics (frequency and spanwise wavelength). In our experiments, we quantitatively examined the problem of linear receptivity of boundary layer to surface nonuniformities in a broad range of frequencies for the most dangerous spanwise scales of Görtler vortices. The values of the amplitudes and phases of the receptivity coefficients were determined. The amplitudes proved to be much smaller in magnitude in comparison with the excitation of modes of hydrodynamic instabilities of other types (Tollmien-Schlichting waves and cross-flow-instability modes). It was found that, with increasing the frequency, the amplitudes of the receptivity coefficients showed a distinct growth while for high frequencies those amplitudes also exhibited a growth with the spanwise scale of perturbations, although for stationary surface roughness no effect due to this scale was observed. It was found that the dependences on frequency of the efficiency of the mechanisms of stability and receptivity showed opposing behaviors, were in competition, and could partially compensate each other, promoting, thus, the production of boundary-layer Görtler vortices in a broad range of frequencies.     

Viscous near-wall flow in a wake of circular cylinder at moderate Reynolds numbers

Here we present the results of experimental investigation of a cross flow around a circular cylinder mounted near the wall of a channel with rectangular cross section. The experiments were carried out in the range of Reynolds numbers corresponding to the transition to turbulence in a wake of the cylinder. Flow visualization and SIV-measurements of instantaneous velocity fields were carried out. Evolution of the flow pattern behind the cylinder and formation of the regular vortex structures were analyzed. It is shown that in the case of flow around the cylinder, there is no spiral motion of fluid from the side walls of the channel towards its symmetry plane, typical of the flow around a spanwise rib located on the channel wall. The laminar-turbulent transition in the wake of the cylinder is caused by the shear layer instability.     

High-order non-uniform grid schemes for numerical simulation of hypersonic boundary-layer stability and transition

The direct numerical simulation of receptivity, instability and transition of hypersonic boundary layers requires high-order accurate schemes because lower-order schemes do not have an adequate accuracy level to compute the large range of time and length scales in such flow fields. The main limiting factor in the application of high-order schemes to practical boundary-layer flow problems is the numerical instability of high-order boundary closure schemes on the wall. This paper presents a family of high-order non-uniform grid finite difference schemes with stable boundary closures for the direct numerical simulation of hypersonic boundary-layer transition. By using an appropriate grid stretching, and clustering grid points near the boundary, high-order schemes with stable boundary closures can be obtained. The order of the schemes ranges from first-order at the lowest, to the global spectral collocation method at the highest. The accuracy and stability of the new high-order numerical schemes is tested by numerical simulations of the linear wave equation and two-dimensional incompressible flat plate boundary layer flows. The high-order non-uniform-grid schemes (up to the 11th-order) are subsequently applied for the simulation of the receptivity of a hypersonic boundary layer to free stream disturbances over a blunt leading edge. The steady and unsteady results show that the new high-order schemes are stable and are able to produce high accuracy for computations of the nonlinear two-dimensional Navier–Stokes equations for the wall bounded supersonic flow.     

Onset of turbulence from the receptivity stage of fluid flows

The traditional viewpoint of fluid flow considers the transition to turbulence to occur by the secondary and nonlinear instability of wave packets, which have been created experimentally by localized harmonic excitation. The boundary layer has been shown theoretically to support spatiotemporal growing wave fronts by Sengupta, Rao, and Venkatasubbaiah [Phys. Rev. Lett. 96, 224504 (2006)] by a linear mechanism, which is shown here to grow continuously, causing the transition to turbulence. Here, we track spatiotemporal wave fronts to a nonlinear turbulent state by solving the full 2D Navier-Stokes equation, without any limiting assumptions. Thus, this is the only demonstration of deterministic disturbances evolving from a receptivity stage to the full turbulent flow. This is despite the prevalent competing conjectures of the event being three-dimensional and/or stochastic in nature.