Study on multiple ring-like vortex formation and small vortex generation in late flow transition on a flat plate

A high-order direct numerical simulation (DNS) of flow transition over a flat plate at a free stream Mach number 0.5 has been carried out. During the simulation, we cannot find, according to the classical theory of boundary layer transition, the “hairpin vortex breakdown to smaller structures” in the last stage of flow transition on a flat plate. However, we did discover the so-called spikes as a result of a multibridge or multiring formation. This indicated a large and stable vortex structure which can travel for a long distance. We believe that this is a result of ring heads that are located in an inviscid region. These heads of the “turbulence spot” never seem to break down and persist as a stable structure. In addition, we discovered that the U-shaped vortex is a part of an existing coherent structure instead of a newly generated one. The U-shaped vortex also provides an additional channel to transfer vorticity to the ring from the wall. During travel, the leading primary ring in the front of the spot is sloped and skewed, causing disappearance of the second sweep. As a consequence, no energy is brought down to the lower boundary layer near the vortex ring head of the spot. Thus, the small length scale structures become damped and the existing U-shaped vortex structure becomes distinguishable. So, the question is where do the turbulent small length scale vortices come from? We will address this with a new theory which states that all of small length scales (turbulence) are generated by high-shear (HS) layers rather than being produced by “vortex breakdown.” The new DNS shows that the HS layers are produced by strong positive spikes surrounded by low-speed fluids and by the interaction between the secondary and higher-level vortices and the wall surface especially near the ring neck. This multiple ring-like vortex generation also follows the first Helmholtz vortex conservation law. Furthermore, the Λ-shaped vortex is formed and rolling up, and the Λ-vortex is stretched and narrowed, and a new bridge is generated after the neck. The bridge will further become a second ring and so on. Besides the original vortex legs, there are also U-shaped vortex tubes. Finally, the multiple ring vortex structure is formed. From this process, no evidence is found to support that two consequent multiring circles are mixed to generate small vortices. The connection of downdraft/updraft motions is also studied.     

A study in transitional flat plate boundary layers: measurement and visualization

 This study is concerned with transition in flat plate boundary layer flow. Sets of results are obtained as follows: (1) Very clear pictures of the formation and the development of the butterfly-like structures rather than ∧-structures in the K-regime of boundary layer transition are obtained. (2) A chain of ring like vortices, which generate the high-frequency spikes on the time traces of velocity and still present periodical behaviour, at the tip of each ∧-vortex, which is the part of the butterfly-like structure, are visualized for the first time. (3) A wave-like structure is observed to occupy the whole boundary layer, extending from the near-wall region to the outer edge of the boundary layer. Received: 24 September 1998/Accepted: 24 April 1999     

On the formation of the chain of ring-like vortices in a transitional boundary layer

One mechanism by which laminar flow becomes turbulence is through the instability of vortex interaction. The nonlinear evolution of the interactions between the Λ-vortex and the secondary closed vortex has been investigated. This nonlinear instability process leads to some fundamental changes in the topology of flows. Experimental observations are presented here showing the formation of a chain of ring-like vortices, which generate high-frequency spikes on time traces in a transitional boundary layer. Received: 13 March 2000/Accepted: 5 July 2000     

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.     

Numerical investigation of turbulent channel flow controlled by spatially oscillating spanwise Lorentz force

A formulation of the skin-friction drag related to the Reynolds shear stress in a turbulent channel flow is derived. A direct numerical simulation (DNS) of the turbulent control is performed by imposing the spatially oscillating spanwise Lorentz force. Under the action of the Lorentz force with several proper control parameters, only the periodically well-organized streamwise vortices are finally observed in the near-wall region. The Reynolds shear stress decreases dramatically, especially in the near-wall area, resulting in a drag reduction.     

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

On the mechanism of turbulence suppression by spanwise surface oscillations

The attenuation of turbulent pulsations in near-wall flows by means of spanwise periodic surface oscillation is examined. A direct numerical simulation of the flow in a circular pipe with imposed rotational oscillations has shown that for Re=4000 and the optimal oscillation frequency, the degree of turbulence attenuation increases with increase in the oscillation amplitude until the flow relaminarizes. The estimated optimal frequency ω⁺=0.06. The results of applying the theory of the development of near-wall coherent structures agree qualitatively with those of numerical simulation. It is concluded that the intensity of the pulsations is reduced because the spanwise movements weaken the longitudinal vortices which cause turbulent bursts in near-wall flows. Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 2, pp. 37–44, March–April, 2000. The research was carried out with financial support from the Russian Foundation for Basic Research (project No. 99-01-01095).     

Numerical studies on evolution of secondary streamwise vortices in transitional boundary layers

Formation and evolution of secondary streamwise vortices in the compressible transitional boundary layers over a flat plate are studied using a direct numerical simulation method with high-order accuracy and highly effective non-reflecting characteristic boundary conditions. Generation and development processes of the secondary streamwise vortices in the complicated transitional boundary flow are clearly analyzed based on the of numerical results, and the effects on the formation of the ring-like vortex that is vital to the boundary layer transition are explored. A new mechanism forming the ring-like vortex through the mutual effect of the primary and secondary streamwise vortices is expressed.     

On the role of the smallest scales of a passive scalar field in a near-wall turbulent flow

Role of the smallest diffusive scales of a passive scalar field in the near-wall turbulent flow was examined with pseudo-spectral numerical simulations. Temperature fields were analyzed at friction Reynolds number Re _τ=171 and at Prandtl numbers, Pr=1 and Pr=5.4. Results of direct numerical simulations (DNS) were compared with the under-resolved simulations where the velocity field was still resolved with the DNS accuracy, while a coarser grid was used to describe the temperature fields. Since the smallest temperature scales remained unresolved in these simulations, an appropriate spectral turbulent thermal diffusivity was applied to avoid pile-up at the higher wave numbers. In spite of coarser numerical grids, the temperature fields are still highly correlated with the DNS results, including instantaneous temperature fields. Results point to practically negligible role of the diffusive temperature scales on the macroscopic behavior of the turbulent heat transfer.     

Spatial relation between fluctuating wall pressure and near-wall streamwise vortices in wall bounded turbulent flow

A new view of the spatial relation between fluctuating wall pressure and near-wall streamwise vortices (NWSV) is proposed for wall bounded turbulent flow by use of the direct numerical simulation (DNS) database. The results show that the wall region with low pressure forms just below the strong NWSV, which is mostly associated with the overhead NWSV. The wall region with high pressure forms downstream of the NWSV, which has a good correspondence with the downwash of the fluids induced by the upstream NWSV. The results provide a significant basis for the detection of NWSV.     

Flow Structure of Swirling Turbulent Propane Flames

Flow structure of premixed propane–air swirling jet flames at various combustion regimes was studied experimentally by stereo PIV, CH* chemiluminescence imaging, and pressure probe. For the non-swirling conditions, a nonlinear feedback mechanism of the flame front interaction with ring-like vortices, developing in the jet shear layer, was found to play important role in the stabilisation of the premixed lifted flame. For the studied swirl rates (S = 0.41, 0.7, and 1.0) the determined domain of stable combustion can be divided into three main groups of flame types: attached flames, quasi-tubular flames, and lifted flames. These regimes were studied in details for the case of S = 1.0, and the difference in the flow structure of the vortex breakdown is described. For the quasi-tubular flames an increase of flow precessing above the recirculation zone was observed when increased the stoichiometric coefficient from 0.7 to 1.4. This precessing motion was supposed to be responsible for the observed increase of acoustic noise generation and could drive the transition from the quasi-tubular to the lifted flame regime.     

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.     

Direct numerical simulation of particle transport by hairpin vortices in a laminar boundary layer

The transport of solid particles by coherent wall structures is studied here. This phenomenon is present in numerous environmental and engineering flows. The flow above a wall-mounted hemisphere is used for generating hairpin vortices in a laminar boundary layer in a controlled way. By means of direct numerical simulation (DNS) of the fluid flow and simultaneous Lagrangian tracking of particles, the influence of hairpin vortices on solid particles released in the wake of the obstacle is analyzed.     

Large eddy simulation of gas-particle turbulent channel flow with momentum exchange between the phases

This paper presents results of a large eddy simulation (LES) combined with Lagrangian particle tracking and a point-force approximation for the feedback effect of particles on the downward turbulent gaseous flow in a vertical channel. The LES predictions are compared with the results obtained by direct numerical simulation (DNS) of a finer computational mesh. A parametric study is conducted for particles with two response times in simulations with and without streamwise gravitational settling and elastic, binary interparticle collisions. It is shown that the classical and the dynamic Smagorinsky turbulence models adequately predict the particle-induced changes in the mean streamwise velocity and the Reynolds stresses of the carrier phase for the range of parameters studied. However, the largest discrepancies between the LES and DNS results are found in the cases of particle-laden flows. Conditional sampling of the instantaneous resolved flow fields indicates that the mechanisms by which particles directly oppose the production of momentum and vorticity of the organized fluid motions are also observed in the LES results. However, the geometric features of the near-wall quasistreamwise vortices are overestimated by the use of both turbulence models compared to the DNS predictions.     

Vortex stretching and filaments

This paper presents a new technique to produce controlled stretched vortices. Intense elliptical vortices are created by stretching of an initial vorticity sheet. The initial vorticity comes from a laminar boundary layer flow and the stretching is parallel to the vorticity vectors. This low velocity flow enables direct observation of the formation and destabilization of vortices. Visualizations are combined with quasi-instantaneous measurements of a full velocity profile. The velocity profile is obtained with an ultrasonic pulsed Doppler velocimeter. The evolution of the central diameter of the vortices is related to the stretching. It is observed that destabilization occurs by pairing of two vortices, by hairpin deformation, and by breakdown of vortices into a “coil shape”.     

上海陆域古河道分布及对工程建设影响研究

采用直接数值模拟方法, 对槽道湍流中展向振动流向传播的波动壁面的流动 控制和减阻问题进行了研究, 讨论了流向参数k_{x}对Stokes层、湍流拟序结构、湍流猝 发事件以及壁面阻力的影响, 并对此类波动壁面的湍流控制和减阻机理进行了讨论. 结果表 明, 当此类波动壁面被用来调制近壁流动时, 仅低频波对湍流流场具有显著影响, 可导致湍 流猝发事件的频率和强度的显著变化; 波数k_{x}的增大对于湍流猝发事件的频率和强度增 减的影响并不同步, 存在一个最优的波数k_{x}, 在其调制下, 固有流场对诱导流场的影响 最弱, 而诱导流场对固有流场的影响显著, 减阻效果最好.     

Direct Numerical Simulation of Turbulent Flow over a Rectangular Trailing Edge

This paper describes a direct numerical simulation (DNS) study of turbulent flow over a rectangular trailing edge at a Reynolds number of 1000, based on the freestream quantities and the trailing edge thickness h; the incoming boundary layer displacement thickness δ^* is approximately equal to h. The time-dependent inflow boundary condition is provided by a separate turbulent boundary layer simulation which is in good agreement with existing computational and experimental data. The turbulent trailing edge flow simulation is carried out using a parallel multi-block code based on finite difference methods and using a multi-grid Poisson solver. The turbulent flow in the near-wake region of the trailing edge has been studied first for the effects of domain size and grid resolution. Then two simulations with a total of 256 × 512 × 64 (∼ 8.4×10⁶) and 512 × 1024 × 128 (∼ 6.7×10⁷) grid points in the computational domain are carried out to investigate the key flow features. Visualization of the instantaneous flow field is used to investigate the complex fluid dynamics taking place in the near-wake region; of particular importance is the interaction between the large-scale spanwise, or Kármán, vortices and the small-scale quasi-streamwise vortices contained within the inflow boundary layer. Comparisons of turbulence statistics including the mean flow quantities are presented, as well as the pressure distributions over the trailing edge. A spectral analysis applied to the force coefficient in the wall normal direction shows that the main shedding frequency is characterized by a Strouhal number based on h of approximately 0.118. Finally, the turbulence kinetic energy budget is analysed. Received 4 March 1999 and accepted 27 October 2000     

A New Low-Reynolds-Number One-Equation Model of Turbulence

In this study, we propose a new Low-Reynolds-Number (LRN)one-equation model, which is derived from an LRN two-equation(k-ε) model. The derivation of the transport equation, in principle, is based on the assumption that the turbulent structure parameter remains constant. However, the relation for the turbulent structure parameter a ₁(=|− |/k) is modified to account for near-wall turbulence. As a result, the present one-equation model contains a term which takes the near-wall limiting behavior explicitly into account. Thus, the present model provides the correct wall-limiting behavior of turbulence in the vicinity of the wall and can be applied to the analysis of heat transfer. The validity of the present model is tested in channel flows, boundary layer flows with and without pressure gradient, plane wall jet, and flow with separation and reattachment. The calculated results showed good agreement with the direct numerical simulation (DNS) and experimental data. This revised version was published online in July 2006 with corrections to the Cover Date.     

Spatial resolution correction for hot-wire anemometry in wall turbulence

We investigate spatial resolution issues in hot-wire anemometry measurements of turbulence intensity and energy spectra. Single normal hot-wire measurements are simulated by means of filtering direct numerical simulation (DNS) of turbulent channel flow at Re_t = 934Re_\tau = 934. Through analysis of the two-dimensional energy spectra from the DNS, the attenuation of the small-scale energy levels is documented, especially in the near-wall region. The missing energy displays anisotropic characteristics, and an attempt is made to model this using an empirical equation, thus providing a correction scheme for all wall normal locations. The empirical model is assessed using experimental boundary layer data and shown to effectively correct both the streamwise one-dimensional energy spectra and turbulence intensity at a Reynolds number significantly above that of the DNS.