Numerical study on evolution of subharmonic varicose low-speed streaks in turbulent channel flow

The evolution of two spanwise-aligned low-speed streaks in a wall turbulent flow, triggered by the instability of the subharmonic varicose (SV) mode, is studied by a direct numerical simulation (DNS) method in a small spatial-periodic channel. The results show that the SV low-speed streaks are self-sustained at the early stage, and then transform into subharmonic sinuous (SS) low-speed streaks. Initially, the streamwise vortex sheets are formed by shearing, and then evolve into zigzag vortex sheets due to the mutual induction. As the intensification of the SV low-speed streaks becomes prominent, the tilted streamwise vortex tubes and the V-like streamwise vortex tubes can be formed simultaneously by increasing \( + \frac{{\partial u}}{{\partial x}}\). When the SV low-speed streaks break down, new zigzag streamwise vortices will be generated, thus giving birth to the next sustaining cycle of the SV low-speed streaks. When the second breakdown happens, new secondary V-like streamwise vortices instead of zigzag streamwise vortices will be generated. Because of the sweep motion of the fluid induced by the secondary V-like streamwise vortices, each decayed low-speed streak can be divided into two parts, and each part combines with the part of another streak, finally leading to the formation of SS low-speed streaks.     

纵向涡对近壁湍流结构的影响

在封闭水槽中用氢气泡方法观察了纵向涡对近壁湍流结构的影响，涡的下洗侧出现了展向距离较宽、流动较稳定的快斑区，流向速度快；涡的上洗侧出现了含有慢斑的区域，流动结构复杂，流向速度慢，纵向涡使下洗侧速度较快的流体流向壁面，使上洗侧速度较慢的流体远离壁面。     

A visual study on complex flow structures and flow breakdown in a boundary layer transition

A visualization study is conducted on the excited laminar-turbulent transition within a flat plate boundary layer flow in a water tunnel. The hydrogen bubble technique is employed to investigate the complex characteristics of the flow structure and its breakdown in the later stages of the transition. A new flow structure is observed, which involves two secondary hairpin vortices outboard of both legs of a primary hairpin vortex. This complex structure is argued to be a precursor of a turbulent spot in this K-type transition. Also reported in the paper is the evolution of the flow structure and its subsequent breakdown, manifested by the emergence of dark spots, low-speed fluid bumps, and near-wall hairpin vortex groups. The results indicate that the near-wall flow breakdown is the result of instability of a local three-dimensional high-shear layer between the low-speed fluid bump and the outer higher-speed region.     

Dynamics of low-speed streak evolution and interaction in laminar boundary layer

The present paper presents an experimental effort on the regeneration process of two low-speed laminar streaks in a zero-pressure-gradient laminar boundary layer. Two vertical thin wires separated by a spanwise distance of 30 mm are used to introduce disturbances of two rolls of transitional Karmain vortex street to the downstream boundary layer. Both hydrogen bubble visualization and particle image velocimetry （PIV） measurement show that two lowspeed streaks are induced through leading-edge receptivity process. As these streaks develop in the downstream, two additional low-speed streaks begin to appear outboard of the flank of the original two, together with complex dynamics of streak splitting and merging. A flow pattern of four streaks aligned along the spanwise direction occurs finally in the far downstream. It is found that besides the mechanisms of streak breakdown, the streak interaction is also an important factor characterizing the instability of low speed streaks and their regeneration process.     

On the etiology of shear layer vortices

In an effort to isolate the mechanism by which streamwise structures form in turbulent wall layers, evolution equations were derived for the streamwise velocity and vorticity perturbations about a mean turbulent fully developed channel flow. The stability of these equations, which take their most concise form when derived from the Generalized Lagrangian mean equations of Andrews and McIntyre, are studied assuming normal modes and infinitesimal disturbances. The resulting stability diagram yields, inter alia, the spanwise periodicity of the resulting structures, which we term shear layer vortices. If streaks are thought of as the footprints of these vortices, we then have a formal way of determining the spacing of streaks. The first three modes of instability are determined; at the first not just two vortices form per period, but four. It is also evident that an intense local shear layer forms about the plane in which the convection velocity equals the mean Eulerian velocity.Dedicated to Professor J.L. Lumley on the occasion of his 60th birthday.This work was supported in part by the U.S. Office of Naval Research under SRO IV Grant No. N00014-85-K-0172 and in part by the National Science Foundation Grant CTS-9008477.     

Near-wall structure of three-dimensional turbulent boundary layers

 Modifications to near-wall turbulent boundary layer structure with increased three-dimensionality have been investigated through the use of hydrogen bubble wire flow visualization. Results indicate that three-dimensionality does not influence the strength or sign of near-wall streamwise vortices. Increased three-dimensionality does stabilize the near-wall structure resulting in less ejection type activity. The spanwise spacing between low-speed streaks also decreased slightly with increased cross-flow. Received: 15 October 1996/Accepted: 2 April 1997     

Vortex ring/viscous wall layer interaction model of the turbulence production process near walls

An experimental simulation of the interaction of vortex ring-like eddies with the sublayer of a turbulent boundary layer is investigated. An artificially generated vortex ring interacting with a Stokes' layer enables investigation of the interaction with reproducible initial conditions and in the absence of background turbulence. All of the observed features in the turbulent boundary layer production process such as the streaky structure, the pockets, the hairpin vortices, streak lift-up, oscillation, and breakup, have been observed to form. The model shows us that hairpin vortices can pinchoff and reconnect forming new vortex ring-like eddies. Interestingly, the model includes interactions that occur with low probability in the turbulent boundary layer, but which contribute significantly to transport, and may be the events most readily controllable.     

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).     

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.     

Investigation of turbulence structures in a drag-reduced turbulent channel flow with surfactant additive by stereoscopic particle image velocimetry

In the present study, we employed stereoscopic particle image velocimetry (PIV) to investigate the characteristics of turbulence structures in a drag-reduced turbulent channel flow with addition of surfactant. The tested drag-reducing fluid was a CTAC/NaSal/Water (CTAC: cetyltrimethyl ammonium chloride; NaSal: sodium salicylate) system at 25°C. The weight concentration of CTAC was 30 ppm. Stereoscopic PIV measurement was performed for a water flow (Re=1.1×10⁴) and a CTAC solution flow (Re=1.5×10⁴ with 54% drag reduction) in both the streamwise–spanwise and wall-normal-spanwise planes, respectively. The three-dimensionality of hairpin vortex structures in the near-wall region for wall-bounded turbulent flow was reproduced by conditionally averaging the stereoscopic two-dimensional-three-component velocity fields. A series of wall-normal vortex cores were found to align with the near-wall low-speed streaks with opposite vorticity signals at both sides of the streaks and with the vorticity decreased on average by about one order of magnitude in CTAC solution flow compared with water flow; the spanwise spacing between the near-wall low-speed streaks in the solution flow is increased by about 46%. The streamwise vorticity of the vortex cores appearing in the wall-normal-spanwise plane was also decreased by the use of drag-reducing surfactant additives.     

Dynamics of vortical structures in cylinder/wall interaction with moderate gap ratio

The interaction between the wake of a transverse circular cylinder and the underlying flat-plate boundary layer with a moderate gap ratio G/D=1.0 is investigated using both hydrogen-bubble-based and PIV-based visualization techniques. The spanwise rollers in the cylinder wake are found to be capable of inducing secondary vortices in the near-wall region. The mutual induction from the counter-clockwise rollers, which are closer to the wall, plays a primary role, so that these secondary vortices present linear lift-up motion at first. Their subsequent evolution dominantly determines the characteristics of the wake/boundary-layer interaction. Two different vortex interaction scenarios are observed: the secondary vortices can be either entrained into the rollers or pushed down towards the wall. This leads to a rapid three-dimensional destabilization process, through which streamwise vortices are generated. And it is suggested that these streamwise vortices are the dominant structures to promote the following boundary layer transition.     

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     

Bypass transition receptivity modes

The prediction of bypass transition remains an important problem in many engineering applications. This is largely because there is no suitable theoretical model for bypass transition and predictions are made using empirical models. This paper presents numerical results for the receptivity of a zero pressure gradient boundary layer subjected to simple freestream waveforms which are the constituent parts of a turbulent flow field. Significant receptivities are only obtained for a minority of freestream waveforms and these lead to two types of flow structure in the boundary layer. The first type of flow structure is essentially two dimensional in nature and consists of two rows of counter-rotating spanwise vortices and is induced by freestream waves of large normal and spanwise wavelength and streamwise wavelengths approximately equal to the boundary layer thickness. The second type of flow structure are the streamwise streaks frequently observed in flow visualisation experiments. These streaks are induced by freestream waves of long streamwise and normal wavelength and spanwise wavelengths in the range of 14.5-46 θ (1.7-5.4δ). The freestream waves can be formed of velocity components in any direction, however the boundary layer is most receptive to fluctuations that lie in a plane perpendicular to the streamwise direction. The overall receptivity to a full spectrum of waves typical of freestream turbulence is considered and is shown to have similar characteristics to those from experiments.     

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.     

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).     

Population,characteristics and kinematics of vortices in a confined rectangular jet with a co-flow

Vortex behavior and characteristics in a confined rectangular jet with a co-flow were examined using vortex swirling strength as a defining characteristic. On the left side of the jet, the positively (counterclockwise) rotating vortices are dominant, while negatively rotating vortices are dominant on the right side of the jet. The characteristics of vortices, such as population density, average size and strength, and deviation velocity, were calculated and analyzed in both the cross-stream direction and the streamwise direction. In the near-field of the jet, the population density, average size and strength of the dominant direction vortices show high values on both sides of the center stream with a small number of counter-rotating vortices produced in the small wake regions close to jet outlet. As the flow develops, the wake regions disappear, these count-rotating vortices also disappear, and the population of the dominant direction vortices increase and spread in the jet. The mean size and strength of the vortices decrease monotonically with streamwise coordinate. The signs of vortex deviation velocity indicate the vortices transfer low momentum to high-velocity region and high momentum to the low velocity region. The developing trends of these characteristics were also identified by tracing vortices using time-resolved particle image velocimetry data. Both the mean tracked vortex strength and size decrease with increasing downstream distance overall. At the locations of the left peak of turbulent kinetic energy, the two-point spatial cross-correlation of swirling strength with velocity fluctuation and concentration fluctuation were calculated. All the correlation fields contain one positively correlated region and one negatively correlated region although the orientations of the correlation fields varied, due to the flow transitioning from wake, to jet, to channel flow. Finally, linear stochastic estimation was used to calculate conditional structures. The large-scale structures in the velocity field revealed by linear stochastic estimation are spindle-shaped with a titling stream-wise major axis.     

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 vorticity stretching diagnostic for turbulent and transitional flows

Vorticity stretching in wall-bounded turbulent and transitional flows has been investigated by means of a new diagnostic measure, denoted by Γ, designed to pick up regions with large amounts of vorticity stretching. It is based on the maximum vorticity stretching component in every spatial point, thus yielding a three-dimensional scalar field. The measure was applied in four different flows with increasing complexity: (a) the near-wall cycle in an asymptotic suction boundary layer (ASBL), (b) K-type transition in a plane channel flow, (c) fully turbulent channel flow at Re _τ = 180 and (d) a complex turbulent three-dimensional separated flow. Instantaneous data show that the coherent structures associated with intense vorticity stretching in all four cases have the shape of flat ‘pancake’ structures in the vicinity of high-speed streaks, here denoted ‘h-type’ events. The other event found is of ‘l-type’, present on top of an unstable low-speed streak. These events (l-type) are further thought to be associated with the exponential growth of streamwise vorticity in the turbulent near-wall cycle. It was found that the largest occurrence of vorticity stretching in the fully turbulent wall-bounded flows is present at a wall-normal distance of y ⁺?=?6.5, i.e. in the transition between the viscous sublayer and buffer layer. The associated structures have a streamwise length of ~200–300 wall units. In K-type transition, the Γ-measure accurately locates the regions of interest, in particular the formation of high-speed streaks near the wall (h-type) and the appearance of the hairpin vortex (l-type). In the turbulent separated flow, the structures containing large amounts of vorticity stretching increase in size and magnitude in the shear layer upstream of the separation bubble but vanish in the backflow region itself. Overall, the measure proved to be useful in showing growing instabilities before they develop into structures, highlighting the mechanisms creating high shear region on a wall and showing turbulence creation associated with instantaneous separations.