转捩边界层中的上喷下扫运动

采用高精度直接数值模拟方法,研究平板可压缩流边界层转捩后期过程中的上喷和下扫运动.模拟结果清楚地显示了边界层转捩过程中的上喷下扫运动的复杂现象,展示与典型涡系及尖峰结构的紧密联系.详细剖析多层次的上喷下扫运动的形成机制,数值结果证实了实验的新发现,并进行相应的机理分析.     

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

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

低速轴流涡轮内部流动的大涡模拟

在低进口雷诺数下,低速轴流涡轮内部可能存在复杂的边界层转捩和分离流动。准确模拟边界层转捩和流动分离对低速轴流涡轮的气动设计具有重要意义。本文以某单级低速轴流涡轮为研究对象,采用大涡模拟方法对其在进口雷诺数为20000情况下的内部流动进行了数值模拟研究,并与前期采用全层流模型、S-A模型、Abu-GhannamShaw(AGS)转捩模型的模拟结果进行了对比,对比分析发现,大涡模拟结果与实验结果吻合更好,可以准确模拟该涡轮叶片吸力面的流动分离和叶片通道内的二次流动。由大涡模拟结果可以得出,静叶尾迹和分离使尾迹区的流体流动速度降低,但尾迹对流动角的影响较小。动叶入口低速微团在做周向运动的同时沿径向运动;高速微团主要沿周向运动。静叶叶片表面的分离流存在较大的由叶顶向叶根的径向的运动;动叶吸力面叶顶处也存在较大的分离流动.     

湍流边界层内流向涡的数值模拟研究

采用准二维共振三波作为湍流边界层近壁区相干结构初值,用直接数值模拟方法计算了流动从二维结构发展到三维结构并且伴随流向涡生成的整个过程,分析结果显示流向涡对湍流动能和质量传输有着重要作用,是湍流边界层相干结构的重要特征和运动形式.     

壁面换热和压力梯度条件下边界层演化的LES研究

本文采用大涡模拟(LES)方法对壁面换热边界条件下有压力梯度的分离边界层流动问题进行了数值模拟研究,对不同壁面热边界条件下分离边界层进行了统计特性及瞬态流场分析.结果表明:与绝热壁面条件相比,等温冷壁面条件(0.8倍来流总温)下附着边界层内的速度剖面比绝热壁面更加饱满,边界层形状因子降低、分离泡也被显著抑制,其时均分离泡尺寸无论是在流向还是法向均更小;等温冷壁面条件下分离剪切层展向涡脱落的位置提前,展向涡卷起的频率降低,分离剪切层内相干结构尺度整体减小,促使分离边界层转捩加速和再附位置前移。     

边界层内涡旋-颗粒相互作用的直接数值模拟

本文采用直接数值模拟方法结合拉格朗日跟踪方法研究平板边界层内涡旋与重颗粒的相互作用。模拟结果显示了由半球状粗糙元诱导产生的发夹涡沿流向的演化过程。小惯性的重颗粒在低压区上边缘的聚集阻碍了涡量的聚集并破坏了发夹涡的产生。模拟结果还证实颗粒在边界层内的输运主要受扫过和喷出事件的影响,且在边界层近壁区域内,扫过事件输运颗粒的能力强于喷出事件的。     

压缩拐角激波/湍流边界层干扰特性分析

用大涡模拟方法对Mach数3.0下的压缩拐角激波/湍流边界层干扰问题进行数值研究.对拐角上游平板区域边界层转捩及湍流进行模拟,设定平板区域长度,使得转捩过程于平板区域发生并充分完成,从而在拐角处产生激波/湍流边界层相互干扰,研究激波/湍流边界层的作用机理.研究表明:流场能够在非定常扰动激励下迅速转捩,并于平板区发展为完全湍流;湍流边界层与激波相互作用过程中,拐角附近分离区较层流情况明显减小;展向不同区域分离区大小差异较大,局部区域分离现象消失.     

涡发生器结构对翼型绕流场的影响

为了研究涡发生器在风力机叶片上的应用,以进一步提高风力机气动效率,本文采用CFD数值模拟方法,分析涡发生器几何形状对其绕流场和翼型边界层特性的影响.涡发生器几何形状为同样高度的矩形、梯形和三角形。翼型为风力机专用翼型DU97-W-300。首先对数值模拟结果与实验值进行了对比,验证了数值方法的可信性。然后详细讨论了各种涡发生器所产生的集中涡涡量、翼型边界层特性、以及绕流场等沿流向的发展演变。总体上看,三角形涡发生器较适合用于风力机翼型的流动控制。     

一维周期与准周期排列沟槽结构的流体减阻特性研究

本文设计了具有相同平均沟槽密度的三种排列类型的一维沟槽结构: 密排列、周期间隔排列和两种准周期间隔排列, 并采用数值模拟和实验验证相结合的方法研究了一维沟槽结构在不同排列下的流体减阻特性. 模拟计算分析流场特征和总阻力, 发现相对于密排列和周期间隔排列的沟槽结构, 准周期间隔排列具有更好的减阻特性, 并且这一结果得到减阻实验的验证. 通过流场分布特性进一步分析沟槽结构的减阻机理. 机理分析发现高速流在经一维准周期结构的扰动波调制后形成了准周期间隔排列的速度条纹相, 这有效地抑制了大涡在流向和展向上的形成, 从而实现较大幅度的减阻. 同时对比分析沟槽排列结构调制展向涡和流向涡各自对流动减阻的贡献, 结果表明, 调制流向涡对减阻的作用更大. 关键词： 流体减阻 沟槽结构 准周期     

超声速平板圆台突起物绕流实验和数值模拟研究

高速飞行器表面不可避免的存在突起物并形成复杂流场, 从而引起飞行器气动特性和热载荷的变化; 同时, 突起物是流动控制的重要方法之一, 合适的突起物形状及安装位置对于改善冲压发动机进气道性能有重要意义. 本文采用基于纳米粒子的平面激光散射技术(NPLS)研究了马赫3.0来流边界层为层流的平板上三个不同高度圆台突起物绕流流场, 主要关注了突起物后方的尾迹边界层, 并采用高精度的显式五阶精度加权紧致非线性格式(WCNS-E-5)离散求解Navier-Stokes方程模拟了该流场. 获得了超声速圆台绕流精细流场结构, 观察到突起物后方尾迹区域边界层发展的过程. 结合实验和数值模拟结果可以发现, 当圆台高度接近或者小于当地边界层厚度时, 突起物对边界层的扰动非常弱, 圆台后方尾迹边界层能够维持较长距离的层流状态, 在边界层转捩阶段也有清晰的发卡涡结构出现; 反之, 边界层受到的扰动明显增大, 在突起物后方很快发展为湍流; 风洞噪声对本文研究圆台引起的边界层扰动有一定影响, 实验获得的边界层转捩位置要比数值结果靠前. 基于NPLS流场图像, 采用间歇性方法分析了圆台突起物后方边界层的特性, 对于高度大于边界层厚度的圆台其间歇性曲线较为接近并且更加饱满, 边界层的脉动也更为强烈.     

Visualization of sinusoidal and varicose instabilities of streaks in a boundary layer

Nonlinear instabilities of boundary layer streaks are investigated experimentally. Extensive measurements visualizing the sinusoidal and varicose instabilities of streaky structures at nonlinear stage of the breakdown process in boundary layer are presented. The flow behaviour in the course of spatial evolution of the streaky structures with a secondary high-frequency disturbance generated on them is discussed. Various scenarios of origination and development of coherent vortex structures examined in physical experiments are considered. Specific features of the development of sinusoidal and varicose cases of destruction of the steady streamwise streaks are demonstrated, such as transverse and streamwise modulations of the streak by the secondary-disturbance frequency, appearance of new streaky structures in the downstream direction, and emergence and evolution of unsteady Λ-shaped structures localized in space in both cases.     

A new kind of hairpin-like vortical structure induced by cross-interaction of sinuous streaks in turbulent channel

This work is motivated by previous experimental and numerical studies which reveal that the hairpin vortex could be formed by the interaction between spanwise adjacent low-speed streaks. To prove that such an interaction mechanism is still applicable in the normal direction, two sinuous low-speed streaks with the same streamwise phase are set to be in the upper half and bottom half of a small size channel, respectively, and their evolution and interaction are investigated by direct numerical simulation. A new kind of hairpin-like vortical structure, distributed in the normal direction and straddled across both halves of the channel, is found during the cross-interaction process of the low-speed streaks. The influence of such a normal-distributed hairpin-like vortex(NHV) on the turbulent statistical regularity is also revealed. It is observed that the NHV can lead to a sudden surge of wall skin friction, but the value of the normal velocity as well as the streamwise and spanwise vorticity sharply decrease to zero in the center of the channel.     

Origin of effectiveness degradation in active drag reduction control of turbulent channel flow at Reτ = 1000

Direct numerical simulations of turbulent channel flows are performed with opposition control at Re_τ = 180 and 1000. The drag reduction rate at the higher Reynolds number is reduced by 25% compared with that at the lower Reynolds number. In order to investigate the reason for the degradation of the control effectiveness, we examine the response of Reynolds stresses and coherent structures in both the outer and inner regions to the control and the role that large-scale motions play therein. In the outer region, the Reynolds stresses at different length scales are reduced at the same rate as the drag reduction rate, and conditionally averaged large-scale motions with spanwise scale larger than half channel width are still large-scale low-speed streaks flanked by a pair of large-scale counter-rotating streamwise vortices but with reduced velocity amplitudes. In the inner region, the effectiveness of the control in suppressing the turbulence deteriorates at the higher Reynolds number. In response to the superimposition effect of large-scale motions, the contribution to near-wall wall-parallel velocity fluctuations from large-scale motions becomes larger at the higher Reynolds number, while the suppression of large-scale motions by the control is weaker than that of near-wall coherent structures. In both controlled and uncontrolled cases, large-scale motions can modulate the amplitudes of near-wall coherent structures, and the attenuation of streamwise vortices by the control under large-scale high-speed streaks is significantly less effective than that under large-scale low-speed streaks. As a result, the effectiveness of control in suppressing near-wall coherent structures and Reynolds shear stresses becomes weaker at the higher Reynolds number. The quantitative analysis of the contributions to the drag reduction rate from outer and inner regions shows that the effectiveness of the control is mainly determined by the suppression degree of near-wall motions. Furthermore, budgets of streamwise enstrophy are analysed to reveal the interaction of large-scale motions with near-wall streamwise vorticity. The titling effect induced by large-scale motions is positive under large-scale high-speed streaks, but negative under large-scale low-speed streaks, which could be a possible way of large-scale motion to modulate streamwise vorticity. In the controlled cases, the positive titling effect induced by large-scale motions under large-scale high-speed streaks is even enhanced, while other terms in the budgets are reduced, which could explain the degradation of control effectiveness in suppressing near-wall streamwise vortices under large-scale high-speed streaks. Therefore, the loss in the drag reduction rate at the higher Reynolds number is due to the weakened control effectiveness on near-wall coherent structures, which are exposed to the modulation effect of large-scale motions.     

高聚物减阻溶液对壁湍流输运过程的影响

基于相同雷诺数下清水和高分子聚合物溶液壁湍流的高时间分辨率粒子图像测速技术(time-resolved particle image velocimetry, TRPIV)的对比实验, 从高聚物溶液对湍流边界层动量能量输运影响的角度分析其减阻的机理. 对比两者的雷诺应力发现高聚物的存在抑制了湍流输运过程. 这一影响与高聚物对壁湍流中占主导地位的涡旋运动和低速条带等相干结构的作用密切相关. 运用条件相位平均、相关函数和线性随机估计(linear stochastic estimation, LSE)等方法, 分析提取了高聚物溶液流场中的发卡涡和发卡涡包等典型相干结构的空间拓扑形态. 相比于清水, 高聚物溶液中相干结构的流向尺度增大, 涡旋运动的发展及低速流体喷射的强度受到削弱, 表明了添加的高聚物阻碍了湍流原有的能量传递和自维持的机理. 正是通过影响相干结构, 高聚物抑制了湍流边界层中近壁区与外区之间的动量和能量输运, 使得湍流的无序性降低, 从而减小了湍流流动的阻力.     

Visualizing shear stress in Görtler vortex flow

Shear stress distributions were obtained from velocity measurements in a concave surface boundary layer flow in the presence of Görtler vortices by means of a single hot-wire probe for several streamwise (x) locations. A set of vertical wires of 0.20 mm diameter were positioned at a distance of 10 mm upstream from the leading edge of a concave surface of radius of curvature R=1.0 m to pre-set the wavelength of the vortices so to obtain the most amplified wavelength Görtler vortices. Consequently, the wavelength of the vortices was set equal to the wire spacing and preserved downstream. In addition to the high shear regions near the wall, one positive peak at the head of the mushroom-like structures and two relatively weak negative peaks at the vicinity of the low-speed streaks are found in the iso-?u/?y contours. They are believed to be related to the formation of the inflectional point in the velocity profile across boundary layer. The occurrence of the inflection points in the spanwise distributions of streamwise velocity component u is associated with the appearance of the second peak of the ?u/?z shear near the boundary layer edge. The nonlinear effect of Görtler instability is to increase the wall shear stress, and further enhancement beyond the turbulent values is due to the presence of secondary instability.     

关于转捩边界层中流向涡的产生

给出了转捩边界中流向涡产生的物理过程,这个过程是A－涡和二次涡环相互作用产生旋涡轴向失稳断裂造成的。     

Experimental investigation on drag reduction in a turbulent boundary layer with a submerged synthetic jet

This work investigates the active control of a fully developed turbulent boundary layer by a submerged synthetic jet actuator.The impacts of the control are explored by measuring the streamwise velocities using particle image velocimetry,and reduction of the skin-friction drag is observed in a certain range downstream of the orifice.The coherent structure is defined and extracted using a spatial two-point correlation function,and it is found that the synthetic jet can efficiently reduce the streamwise scale of the coherent structure.Proper orthogonal decomposition analysis reveals that large-scale turbulent kinetic energy is significantly attenuated with the introduction of a synthetic jet.The conditional averaging results show that the induction effect of the prograde vortex on the low-speed fluid in a large-scale fluctuation velocity field is deadened,thereby suppressing the bursting process near the wall.     

ILES for mechanism of ramp-type MVG reducing shock induced flow separation

A high order implicit large eddy simulation (ILES) is carried out to study the mechanism of shock induced flow separation reduction under ramp-type MVG control. The mechanism was originally considered as that MVG can generate streamwise vortices which strongly mix boundary layer and the boundary layer becomes more capable to resist strong adverse pressure gradient caused by shock and to keep the boundary layer attached. However, according to our ILES, a chain of ring-like vortices is generated behind the ramp-type MVG and goes further to interact with the shock. When the ring-like vortices pass through the shock, the shock wave is weakened and altered while the vortex structures are quite stable. The instantaneous simulation shows that the spanwise ring-like vortex, not the streamwise vortex, plays a key role to weaken the shock and reduce the shock-induced separation. Detailed investigation on ring-like vortices and shock interaction will be given in this paper.     

The investigation of coherent structures in the wall region of a supersonic turbulent boundary layer based on DNS database

Through temporal mode direct numerical simulation, flow field database of a fully developed turbulent boundary layer on a flat plate with Mach number 4.5 and Reynolds number Reθ =1094 has been obtained. Commonly used detection meth- ods in experiments are applied to detecting coherent structures in the flow field, and it is found that coherent structures do exist in the wall region of a supersonic turbulent boundary layer. The detected results show that a low-speed streak is de- tected by using the Mu-level method, the rising parts of this streak are detected by using the second quadrant method, and the crossing regions from a low-speed streak to the high-speed one are detected by using the VITA method respectively. Notwithstanding that different regions are detected by different methods, they are all accompanied by quasi-stream-wise vortex structures.