Cross-Stream Vorticity Field Induced by Streamwise Vortices in Unbounded and Wall-Bounded Shear Flows

This computational study examines the unsteady cross-stream vorticity structures that form when one or more streamwise vortices are immersed in homogeneous and boundary-layer shear flows. A quasi-two-dimensional limit is considered in which the velocity and vorticity fields, while still possessing three nonzero components, have vanishing gradient in the streamwise direction. This idealization is suitable to applications such as streamwise vortices that occur along a ship hull or airplane fuselage and it can be used as an idealized representation of the quasi-streamwise vortices in the near-wall region of a turbulent boundary layer. In this quasi-two-dimensional idealization, the streamwise velocity has no effect on the cross-stream velocity associated with the vortex. However, the vortex acts to modify the cross-stream vorticity component, resulting in regions of the flow with strong deviations in streamwise velocity. This paper examines the complex structures that form as the cross-stream vorticity field is wrapped up by the vortex and the effect of these structures on the streamwise velocity field, first for vortices immersed in homogeneous shear flow and then for vortices immersed in a boundary layer along a flat wall. Received 2 January 2002 and accepted 13 August 2002 Published online 3 December 2002 RID="*" ID="*" This project was supported by the Office of Naval Research under Grant Number N00014-01-1-0015. Dr. Thomas Swain is the program manager. Communicated by T.B. Gatski     

Maintenance of Oscillations in Three-Dimensional Traveling Waves in Plane Poiseuille Flows

Two solutions of three-dimensional Navier–Stokes equations are studied numerically. These solutions describe the fluid motion in a plane channel, are of the traveling-wave form, and are periodic in the streamwise and spanwise directions. It is shown that, in each solution, the oscillations arise as a result of linear instability in the streamwise averaged velocity field. This instability is due to the existence of streamwise streaks known as the regions where the velocity is higher or lower than the mean velocity. A mechanism for the maintenance of streamwise vortices causing the formation of streaks is revealed. The obtained results confirm and extend the existing knowledge about the mechanism for the formation of near-wall turbulent structures.     

An experimental study on turbulent coherent structures near a sheared air-water interface

The turbulence structures near a sheared air-water interface were experimentally investigated with the hydrogen bubble visualization technique. Surface shear was imposed by an airflow over the water flow which was kept free from surface waves. Results show that the wind shear has the main influence on coherent structures under air-water interfaces. Low- and high- speed streaks form in the region close to the interface as a result of the imposed shear stress. When a certain airflow velocity is reached, “turbulent spots” appear randomly at low-speed streaks with some characteristics of hairpin vortices. At even higher shear rates, the flow near the interface is dominated primarily by intermittent bursting events. The coherent structures observed near sheared air-water interfaces show qualitative similarities with those occurring in near-wall turbulence. However, a few distinctive phenomena were also observed, including the fluctuating thickness of the instantaneous boundary layer and vertical vortices in bursting processes, which appear to be associated with the characteristics of air-water interfaces. The project supported by the National Natural Science Foundation of China (Grant No.19672070)     

Tomographic-PIV measurement of the flow around a zigzag boundary layer trip

Tomographic-PIV was used to measure the boundary layer transition forced by a zigzag trip. The resulting instantaneous three-dimensional velocity distributions are used to quantitatively visualize the flow structures. They reveal undulating spanwise vortices directly behind the trip, which break up into individual arches and then develop into the hairpin-like structures typical of wall-bounded turbulence. Compared to the instantaneous flow structure, the structure of the average velocity field is very different showing streamwise vortices. Such streamwise vortices are often associated with the low-speed streaks occurring in bypass transition flows, but in this case clearly are an artifact of the averaging. Rather, the present streaks in the separated flow region directly behind the trip are resulting from the waviness in the spanwise vortices as introduced by the zigzag trip. Furthermore, these streaks and the separated flow region are observed to be related to a large-scale, spanwise uniform unsteadiness in the flow that contributes significantly to the velocity fluctuations over large downstream distances (up to at least the edge of the present measurement domain).     

Using digital holographic microscopy for simultaneous measurements of 3D near wall velocity and wall shear stress in a turbulent boundary layer

A digital holographic microscope is used to simultaneously measure the instantaneous 3D flow structure in the inner part of a turbulent boundary layer over a smooth wall, and the spatial distribution of wall shear stresses. The measurements are performed in a fully developed turbulent channel flow within square duct, at a moderately high Reynolds number. The sample volume size is 90 × 145 × 90 wall units, and the spatial resolution of the measurements is 3–8 wall units in streamwise and spanwise directions and one wall unit in the wall-normal direction. The paper describes the data acquisition and analysis procedures, including the particle tracking method and associated method for matching of particle pairs. The uncertainty in velocity is estimated to be better than 1 mm/s, less than 0.05% of the free stream velocity, by comparing the statistics of the normalized velocity divergence to divergence obtained by randomly adding an error of 1 mm/s to the data. Spatial distributions of wall shear stresses are approximated with the least square fit of velocity measurements in the viscous sublayer. Mean flow profiles and statistics of velocity fluctuations agree very well with expectations. Joint probability density distributions of instantaneous spanwise and streamwise wall shear stresses demonstrate the significance of near-wall coherent structures. The near wall 3D flow structures are classified into three groups, the first containing a pair of counter-rotating, quasi streamwise vortices and high streak-like shear stresses; the second group is characterized by multiple streamwise vortices and little variations in wall stress; and the third group has no buffer layer structures.     

Streamwise streaks and secondary instability in a two-dimensional bent channel

Streamwise streaks generated from a pair of oblique waves and secondary instability of the streaks are studied in a two-dimensional bent channel. Nonlinear parabolized stability equations (NPSE) are employed to investigate streamwise streaks and vortices. A pair of oblique waves from linear stability analysis is imposed as initial disturbances. Generation of streamwise streaks and vortices and subsequent development are described in detail. The case of plane channel is also studied to provide comparable data. Through comparison, the effect of bent region is clearly highlighted. Results of parametric studies to examine the effect of Reynolds number, radius of curvature, and bent angle are also given and discussed in detail. Secondary instability analysis for the modified mean flow due to the streamwise streaks is carried out by solving a two-dimensional eigenvalue problem. Several unstable modes which can be classified into fundamental and subharmonic mode of secondary instability are identified. Among several unstable modes, two modes are turned out to be dominant modes. Details on these two modes including generation mechanism, typical pattern, and dependency on wave number and streak amplitude are discussed. It is found that the presence of bent channel can lead to early oblique-mode breakdown via strong growth of the streamwise streaks due to the curved section. Such large amplitude of streaks and its secondary instability eventually could trigger transition even for small amplitude oblique waves at subcritical channel Reynolds numbers.     

Numerical simulation on evolution of subharmonic low-speed streaks in minimal channel turbulent flow

The evolution of low-speed streaks in the turbulent boundary layer of the minimum channel flow unit at a low Reynolds number is simulated by the direct numer- ical simulation （DNS） based on the standard Fourier-Chebyshev spectral method. The subharmonic sinuous （SS） mode for two spanwise-aligned low-speed streaks is excited by imposing the initial perturbations. The possibilities and the physical realities of the turbulent sustaining in the minimal channel unit are examined. Based on such a flow field environment, the evolution of the low-speed streaks during a cycle of turbulent sus- taining, including lift-up, oscillation, and breakdown, is investigated. The development of streamwise vortices and the dynamics of vortex structures are examined. The results show that the vortices generated from the same streak are staggered along the streamwise direction, while the vortices induced by different streaks tilt toward the normal direction due to the mutual induction effect. It is the spatial variations of the streamwise vortices that cause the lift-up of the streaks. By resolving the transport dynamics of enstrophy, the strength of the vortices is found to continuously grow in the logarithmic layer through the vortex stretching mechanism during the evolution of streaks. The enhancement of the vortices contributes to the spanwise oscillation and the following breakdown of the low-speed streaks.     

Near-field measurements and development of a new boundary layer over a flat plate with localized suction

Suction on a turbulent boundary layer is applied through a narrow strip in order to understand the effects suction can have on the boundary layer development and turbulent structures in the flow. Detailed two-component laser Doppler velocimetry (LDV) and laser-induced fluorescence (LIF) based measurements have been undertaken in regions close to the suction strip and further downstream. The region close to the strip involves a flow reversal accompanied by a change in sign for the Reynolds shear stress and strong gradients in the flow variables. The mean streamwise velocity after suction remains larger than its corresponding no-suction value. Relative to the no-suction case, the velocity fluctuations first decrease with suction followed by a slow recovery which may involve a slight overshoot. LIF visualizations indicate that compared to the no-suction case, the low-speeds streaks stay closer to the wall and exhibit a smaller amount of spanwise and wall-normal oscillations with suction. The visualization results are consistent with two-point velocity correlation measurements. The streamwise and spanwise correlation measurements indicate that the structures are disrupted or removed from the boundary layer due to suction suggesting that the original boundary layer has been strongly influenced by suction. The results are explained by the development of a new inner layer that forms downstream of the suction strip.     

微型射流涡流器对湍流边界层控制的数值研究

运用数值方法,模拟出展向分布的同向倾斜微型射流列与平板湍流边界层相互作用形成流向涡列的流场结构,验证了利用其来对湍流边界层进行控制的可能性.随射流间距减小,流向涡列控制作用流向渗透能力增强,但作用区域减小;随射流速度提高,流向涡列控制作用增强,但过大的射流速度反而会导致流向涡列在局部区域内控制作用的下降;随射流俯仰角减小、倾斜角增大,流向涡列初始控制作用增强,但过小的俯仰角、过大的倾斜角会导致流向涡列流向控制区域明显缩小.要保证流向涡列具有较强的湍流边界层控制作用,必须通过合理配置射流列各主要参数,在保证各流向涡具有一定强度的同时,还要确保各流向涡在形成时部分嵌入边界层内部.     

Three-dimensional instability analysis of boundary layers perturbed by streamwise vortices

A parametric study is presented for the incompressible, zero-pressure-gradient flat-plate boundary layer perturbed by streamwise vortices. The vortices are placed near the leading edge and model the vortices induced by miniature vortex generators (MVGs), which consist in a spanwise-periodic array of small winglet pairs. The introduction of MVGs has been experimentally proved to be a successful passive flow control strategy for delaying laminar-turbulent transition caused by Tollmien–Schlichting (TS) waves. The counter-rotating vortex pairs induce non-modal, transient growth that leads to a streaky boundary layer flow. The initial intensity of the vortices and their wall-normal distances to the plate wall are varied with the aim of finding the most effective location for streak generation and the effect on the instability characteristics of the perturbed flow. The study includes the solution of the three-dimensional, stationary, streaky boundary layer flows by using the boundary region equations, and the three-dimensional instability analysis of the resulting basic flows by using the plane-marching parabolized stability equations. Depending on the initial circulation and positioning of the vortices, planar TS waves are stabilized by the presence of the streaks, resulting in a reduction in the region of instability and shrink of the neutral stability curve. For a fixed maximum streak amplitude below the threshold for secondary instability (SI), the most effective wall-normal distance for the formation of the streaks is found to also offer the most stabilization of TS waves. By setting a maximum streak amplitude above the threshold for SI, sinuous shear layer modes become unstable, as well as another instability mode that is amplified in a narrow region near the vortex inlet position.     

The Relaminarization Mechanisms of Turbulent Channel Flow at Low Reynolds Numbers

The mechanisms of laminarization in wall-bounded flows have been investigated by performing direct numerical simulations (DNS) of turbulent channel flows. By decreasing Reynolds numbers systematically, the effects of the low Reynolds number are studied in connection with the near-wall turbulent structure and turbulent statistics. At approximately the critical Reynolds number, the turbulent skin friction is reduced, and the turbulent structure changes qualitatively in the very near-wall region. Instantaneous turbulent structures reveal that streamwise vortices, the cores of which are at y+ 10, disappear, although low speed streaks and Reynolds shear stress are still produced by larger streamwise vortices located in the buffer region y+ > 10. Sweep motions induced by these vortical structures are shifted toward the center of a channel and also significantly deterred, which may heighten the effects of the viscous sublayer over most of the channel section and suppress the regeneration mechanisms of new streamwise vortices in the very near-wall region. To investigate the details of how large-scale coherent vortices affect the viscous sublayer and the relevant small-scale streamwise vortices, a body force is virtually imposed in the wall-normal direction to enhance the large streamwise vortices. As a result, it is found that when they are sufficiently enhanced, the small-scale vortices reappear, and the sweep events are again dominant in the viscous sublayer.     

Temporal development of turbulent boundary layers with embedded streamwise vortices

The interaction of streamwise vortices with turbulent boundary layer has been investigated using large-eddy simulation. The initial conditions are a pair of counterrotating Oseen vortices with flow between them directed toward the wall (common-flow-down), superimposed on various instantaneous realizations of a turbulent boundary layer. The time development of the vortices and their interaction with the boundary layer are studied by integrating the filtered Navier-Stokes equations in time. The most important effects of the vortices on the boundary layer are the thinning of the boundary layer between vortices (downwash region) and the thickening of the boundary layer in the upwash region. The vortices first move toward the wall as a result of the self-induced velocity, and then apart from each other because of the image vortices due to the solid wall. The Reynolds stress profiles highlight the highly three-dimensional structure of the turbulent boundary layer modified by the vortices. The presence of significant turbulent activity near the vortex center and in the upwash region suggests that localized instability mechanisms in addition to the convection of turbulent energy by the secondary flow are responsible for this effect. High levels of turbulent kinetic energy and secondary stresses in the vicinity of the vortex center are also observed. The numerical results show good agreement with experimental results.This work was supported by the Office of Naval Research under Grant N00014-89-J-1638. Computer time was supplied by the San Diego Supercomputing Center.     

The Effect of Pressure Gradient on Boundary Layer Receptivity

This paper presents numerical results for the receptivity of three laminar boundary layers with zero (ZPG), adverse (APG) and favourable (FPG) pressure gradients. Each boundary layer is subjected to a series of simple freestream waveforms which can be considered as constituent parts of either an isotropic or a non-isotropic turbulent freestream. Each freestream waveform has a single frequency in each spatial direction and is divided into two mutually perpendicular components. The first component has a zero spanwise velocity and hence lies in the streamwise normal plane whereas the second component lies in a plane which is perpendicular both to this plane and the spatial frequency vector. High boundary layer receptivities are only obtained for a minority of these waveforms and so only the resulting flow structures for these waveforms are considered in detail. The dominant flow structures are identified as either Tollmien Schlichting (T-S) waves or streaky structures. The streaky structures can be induced by both freestream components, but the response to the second component, which results in streamwise vortices in the freestream, is considerably stronger and occurs over a much larger streamwise frequency range. The boundary layer is only receptive to a relatively narrow band of spanwise wavelengths ranging from approximately one to four times the local boundary layer thickness. The APG leads to receptivities which are more than double those for the FPG case. The ratio of the freestream fluctuation streamwise wavelength to the distance from the plate leading edge is identified as an important influential parameter for receptivity leading to streaks. Significant T-S activity is only observed for APG, but is also detected for ZPG.     

The effect of streamwise vortices on the turbulence structure of a separating boundary layer

A high Reynolds number flat plate turbulent boundary layer is investigated in a wind-tunnel experiment. The flow is subjected to an adverse pressure gradient which is strong enough to generate a weak separation bubble. This experimental study attempts to shed some new light on separation control by means of streamwise vortices with emphasize on the change in the boundary layer turbulence structure. In the present case, counter-rotating and initially non-equidistant streamwise vortices become and remain equidistant and confined within the boundary layer, contradictory to the prediction by inviscid theory. The viscous diffusion cause the vortices to grow, the swirling velocity component to decrease and the boundary layer to develop towards a two-dimensional state. At the position of the eliminated separation bubble the following changes in the turbulence structure were observed. The anisotropy state in the near-wall region is unchanged, which indicates that it is determined by the presence of the wall rather than the large scale vortices. However, the turbulence in the outer part of the boundary layer becomes overall more isotropic due to an increased wall-normal mixing and a significantly decreased production of streamwise fluctuations. The turbulent kinetic energy is decreased as a consequence of the latter. Despite the complete change in mean flow, the spatial turbulence structure and the anisotropy state, the process of transfer of turbulent kinetic energy to the spanwise fluctuating component seems to be unchanged. Local regions of anisotropy are strongly connected to maxima in the turbulent production. For example, at spanwise positions in between those of symmetry, the spanwise gradient of the streamwise velocity cause significant production of turbulent fluctuations. Transport of turbulence in the spanwise direction occurs in the same direction as the rotation of the vortices.     

On the nature of the organized structures in turbulent near-wall flows

A physical mechanism of onset of large-scale organized structures in turbulent flows along a plane wall which are the cause of intensification of turbulent fluctuations is formulated. The structures take the form of high-speed and low-speed streaks caused by streamwise vortices, i.e., motions in the plane of the transverse cross-section. The streamwise vortices are excited as a result of instability under the action of the anisotropy of the normal components of the Reynolds stress tensor. A model for describing these vortices that gives characteristics in qualitative and quantitative agreement with the experimental data is proposed. In particular, the most probable and mean distances between neighboring vortices are correctly reproduced. The theory makes it possible to explain certain methods of turbulent flow control for the purpose of drag reduction. Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 1, pp. 24–30, January–February, 1997. The work was carried out with financial support from the Russian Foundation for Fundamental Research (project No. 96-01-00602).     

Numerical studies of the instability and breakdown of a boundary-layer low-speed streak

The experimental configuration in [M. Asai, M. Minagawa, M. Nishioka, The instability and breakdown of a near-wall low-speed streak, J. Fluid Mech. 455 (2002) 289–314] is numerically reproduced in order to examine the instability of a single low-speed streak in a laminar boundary layer and to investigate the resulting generation of coherent structures. Such a configuration is chosen since the experimental data show that the two instability modes, varicose and sinuous, are of comparable strength. The instability characteristics are retrieved from the simulation of the flow impulse response. The varicose instability is associated to higher frequencies and lower group velocities than those of the sinuous modes. The latter are less affected by the diffusion of the streak mean shear and are amplified for a longer streamwise distance. Analysis of the perturbation kinetic energy production reveals that both the varicose and the sinuous instability are driven by the work of the Reynolds stress against the wall-normal shear of the streak. The base flow considered here therefore presents an exception to the common knowledge, supported by several previous studies, that the sinuous instability is associated to the streak spanwise shear. The vortical structures at the late stage of the varicose breakdown are identified from the numerical data. By comparing them with those pertaining to other transition scenarios, it is confirmed that streaks and streamwise vortices are universal features of boundary layer transition.     

Vortex dynamics mechanisms of separated boundary layer in a highly loaded low pressure turbine cascade

This paper investigates the vortex dynamics in the suction-side boundary layer on an aero-engine low pressure turbine blade at two different Reynolds numbers at which short and long laminar separation bubbles occur. Different vortical patterns are observed and investigated through large eddy simulation (LES). The results show that at the higher Reynolds number, streamwise streaks exist upstream of separation line. These streaks initiate spanwise undulation in the form of vortex tubes, which roll-up and shed from the shear layer due to the Kelvin–Helmholtz instability. The vortex tubes alternately pair together and eventually distort and break down to small-scale turbulence structures near the mean reattachment location and convect into a fully turbulent boundary layer. At the lower Reynolds number, streamwise streaks are strong and the separated flow is unable to reattach to the blade surface immediately after transition to turbulence. Therefore, bursting of short bubbles into long bubbles can occur, and vortex tubes have larger diameters and cover a part of the blade span. In this case vortex pairing does not occur and vortex shedding process is promoted mainly by flapping phenomenon. Moreover, the results of dynamic mode decomposition (DMD) analysis show a breathing motion as a source of unsteadiness in the separation location, which is accompanied by the flapping phenomenon.     

基于Navier-Stokes方程残差的隐式大涡模拟有限元模型

对应于湍流的大尺度与小尺度流场信息,本文在有限元的框架下,假设Navier-Stokes方程的解的形函数由大尺度和不可解尺度形函数叠加组成,引入对应的权函数,将Navier-Stokes方程的有限元变分形式分解为大尺度和不可解尺度系统.根据不可解尺度系统,构建基于Navier-Stokes大尺度方程残差的不可解尺度模型,将其代入Navier-Stokes方程的大尺度系统,进而数值求解大尺度系统得到Navier-Stokes方程的大尺度解.该方法无需像传统的大涡模拟方法那样对方程的解进行过滤,通过对形函数进行尺度分解实现解的尺度分解.本文使用该方法的自编程序代码开展了槽道湍流的数值模拟.通过与有限差分大涡模拟、DNS计算结果的比较,发现在使用较少网格情况下该方法预测的平均流向速度在近壁区与DNS数据吻合,在黏性外层略偏高;该方法对雷诺应力预测偏低导致从流向向垂向方向上湍动能输运略偏低.流向速度等值面图显示该方法有效捕捉到了大尺度旋涡结构;同时在近壁区可以观察到明显的低速条带结构.     

A RESIDUAL-BASED UNRESOLVED-SCALE FINITE ELEMENT MODELLING FOR IMPLICT LARGE EDDY SIMULATION 1)

对应于湍流的大尺度与小尺度流场信息, 本文在有限元的框架下, 假设Navier-Stokes方程的解的形函数由大尺度和不可解尺度形函数叠加组成, 引入对应的权函数, 将Navier-Stokes方程的有限元变分形式分解为大尺度和不可解尺度系统. 根据不可解尺度系统, 构建基于Navier-Stokes大尺度方程残差的不可解尺度模型, 将其代入Navier-Stokes方程的大尺度系统, 进而数值求解大尺度系统得到Navier-Stokes方程的大尺度解. 该方法无需像传统的大涡模拟方法那样对方程的解进行过滤, 通过对形函数进行尺度分解实现解的尺度分解. 本文使用该方法的自编程序代码开展了槽道湍流的数值模拟. 通过与有限差分大涡模拟、DNS计算结果的比较, 发现在使用较少网格情况下该方法预测的平均流向速度在近壁区与DNS数据吻合, 在黏性外层略偏高; 该方法对雷诺应力预测偏低导致从流向向垂向方向上湍动能输运略偏低. 流向速度等值面图显示该方法有效捕捉到了大尺度旋涡结构; 同时在近壁区可以观察到明显的低速条带结构.     

Skin-friction drag reduction in turbulent channel flow based on streamwise shear control

It is known that stretching and intensification of a hairpin vortex by mean shear play an important role to create a hairpin vortex packet, which generates the large Reynolds shear stress associated with skin-friction drag in wall-bounded turbulent flows. In order to suppress the mean shear at the wall for high efficient drag reduction (DR), in the present study, we explore an active flow control concept using streamwise shear control (SSC) at the wall. The longitudinal control surface is periodically spanwise-arranged with no-control surface while varying the structural spacing, and an amplitude parameter for imposing the strength of the actuating streamwise velocity at the wall is introduced to further enhance the skin-friction DR. Significant DR is observed with an increase in the two parameters with an accompanying reduction of the Reynolds stresses and vorticity fluctuations, although a further increase in the parameters amplifies the turbulence activity in the near-wall region. In order to study the direct relationship between turbulent vortical structures and DR under the SSC, temporal evolution with initial eddies extracted by conditional averages for Reynolds-stress-maximizing Q2 events are examined. It is shown that the generation of new vortices is dramatically inhibited with an increase in the parameters throughout the flow, causing fewer vortices to be generated under the control. However, when the structural spacing is sufficiently large, the generation of new vortex is not suppressed over the no-control surface in the near-wall region, resulting in an increase of the second- and fourth-quadrant Reynolds shear stresses. Although strong actuating velocity intensifies the near-wall turbulence, the increase in the turbulence activity is attributed to the generation of counter-clockwise near-wall vortices by the increased vortex transport.