scholarly journals Investigation of the roughness-induced transition: global stability analyses and direct numerical simulations

2014 ◽  
Vol 760 ◽  
pp. 175-211 ◽  
Author(s):  
Jean-Christophe Loiseau ◽  
Jean-Christophe Robinet ◽  
Stefania Cherubini ◽  
Emmanuel Leriche

AbstractThe linear global instability and resulting transition to turbulence induced by an isolated cylindrical roughness element of height $h$ and diameter $d$ immersed within an incompressible boundary layer flow along a flat plate is investigated using the joint application of direct numerical simulations and fully three-dimensional global stability analyses. For the range of parameters investigated, base flow computations show that the roughness element induces a wake composed of a central low-speed region surrounded by a three-dimensional shear layer and a pair of low- and high-speed streaks on each of its sides. Results from the global stability analyses highlight the unstable nature of the central low-speed region and its crucial importance in the laminar–turbulent transition process. It is able to sustain two different global instabilities: a sinuous and a varicose one. Each of these globally unstable modes is related to a different physical mechanism. While the varicose mode has its root in the instability of the whole three-dimensional shear layer surrounding the central low-speed region, the sinuous instability turns out to be similar to the von Kármán instability in the two-dimensional cylinder wake and has its root in the lateral shear layers of the separated zone. The aspect ratio of the roughness element plays a key role on the selection of the dominant instability: whereas the flow over thin cylindrical roughness elements transitions due to a sinuous instability of the near-wake region, for larger roughness elements the varicose instability of the central low-speed region turns out to be the dominant one. Direct numerical simulations of the flow past an aspect ratio ${\it\eta}=1$ (with ${\it\eta}=d/h$) roughness element sustaining only the sinuous instability have revealed that the bifurcation occurring in this particular case is supercritical. Finally, comparison of the transition thresholds predicted by global linear stability analyses with the von Doenhoff–Braslow transition diagram provides qualitatively good agreement.

Author(s):  
Joshua R. Brinkerhoff ◽  
Metin I. Yaras

This paper describes numerical simulations of the instability mechanisms in a separation bubble subjected to a three-dimensional freestream pressure distribution. Two direct numerical simulations are performed of a separation bubble with laminar separation and turbulent reattachment under low freestream turbulence at flow Reynolds numbers and streamwise pressure distributions that approximate the conditions encountered on the suction side of typical low-pressure gas-turbine blades with blade sweep angles of 0° and 45°. The three-dimensional pressure field in the swept configuration produces a crossflow-velocity component in the laminar boundary layer upstream of the separation point that is unstable to a crossflow instability mode. The simulation results show that crossflow instability does not play a role in the development of the boundary layer upstream of separation. An increase in the amplification rate and most amplified disturbance frequency is observed in the separated-flow region of the swept configuration, and is attributed to boundary-layer conditions at the point of separation that are modified by the spanwise pressure gradient. This results in a slight upstream movement of the location where the shear layer breaks down to small-scale turbulence and modifies the turbulent mixing of the separated shear layer to yield a downstream shift in the time-averaged reattachment location. The results demonstrate that although crossflow instability does not appear to have a noticeable effect on the development of the transitional separation bubble, the 3D pressure field does indirectly alter the separation-bubble development by modifying the flow conditions at separation.


2007 ◽  
Vol 586 ◽  
pp. 259-293 ◽  
Author(s):  
F. X. TRIAS ◽  
M. SORIA ◽  
A. OLIVA ◽  
C. D. PÉREZ-SEGARRA

A set of complete two- and three-dimensional direct numerical simulations (DNS) in a differentially heated air-filled cavity of aspect ratio 4 with adiabatic horizontal walls is presented in this paper. Although the physical phenomenon is three-dimensional, owing to its prohibitive computational costs the majority of the previous DNS of turbulent and transition natural convection flows in enclosed cavities assumed a two-dimensional behaviour. The configurations selected here (Rayleigh number based on the cavity height 6.4 × 108, 2 × 109 and 1010, Pr = 0.71) are an extension to three dimensions of previous two-dimensional problems.An overview of the numerical algorithm and the methodology used to verify the code and the simulations is presented. The main features of the flow, including the time-averaged flow structure, the power spectra and probability density distributions of a set of selected monitoring points, the turbulent statistics, the global kinetic energy balances and the internal waves motion phenomenon are described and discussed.As expected, significant differences are observed between two- and three-dimensional results. For two-dimensional simulations the oscillations at the downstream part of the vertical boundary layer are clearly stronger, ejecting large eddies to the cavity core. In the three-dimensional simulations these large eddies do not persist and their energy is rapidly passed down to smaller scales of motion. It yields on a reduction of the large-scale mixing effect at the hot upper and cold lower regions and consequently the cavity core still remains almost motionless even for the highest Rayleigh number. The boundary layers remain laminar in their upstream parts up to the point where these eddies are ejected. The point where this phenomenon occurs clearly moves upstream for the three-dimensional simulations. It is also shown that, even for the three-dimensional simulations, these eddies are large enough to permanently excite an internal wave motion in the stratified core region. All these differences become more marked for the highest Rayleigh number.


2013 ◽  
Vol 135 (4) ◽  
Author(s):  
Joshua R. Brinkerhoff ◽  
Metin I. Yaras

This paper describes numerical simulations of the instability mechanisms in a separation bubble subjected to a three-dimensional freestream pressure distribution. Two direct numerical simulations are performed of a separation bubble with laminar separation and turbulent reattachment under low freestream turbulence at flow Reynolds numbers and streamwise pressure distributions that approximate the conditions encountered on the suction side of typical low-pressure gas-turbine blades with blade sweep angles of 0 deg and 45 deg. The three-dimensional (3D) pressure field in the swept configuration produces a crossflow-velocity component in the laminar boundary layer upstream of the separation point that is unstable to a crossflow instability mode. The simulation results show that crossflow instability does not play a role in the development of the boundary layer upstream of separation. An increase in the amplification rate and the most amplified disturbance frequency is observed in the separated-flow region of the swept configuration and is attributed to boundary-layer conditions at the point of separation that are modified by the spanwise pressure gradient. This results in a slight upstream movement of the location where the shear layer breaks down to small-scale turbulence and modifies the turbulent mixing of the separated shear layer to yield a downstream shift in the time-averaged reattachment location. The results demonstrate that although crossflow instability does not appear to have a noticeable effect on the development of the transitional separation bubble, the 3D pressure field does indirectly alter the separation-bubble development by modifying the flow conditions at separation.


Author(s):  
Stepan Tolkachev ◽  
Victor Kozlov ◽  
Valeriya Kaprilevskaya

In this article, the results of research about stationary and secondary disturbances development behind the localized and two-dimensional roughness elements are presented. It is shown that the two-dimensional roughness element has a destabilizing effect on the disturbances induced by the three-dimensional roughness element lying upstream. In this case, the two-dimensional roughness element causes the appearance of stationary structures, and then secondary perturbations, whose frequency range lies lower than in the case of the stationary vortices excited by a three-dimensional roughness element.


2008 ◽  
Vol 602 ◽  
pp. 175-207 ◽  
Author(s):  
L. E. JONES ◽  
R. D. SANDBERG ◽  
N. D. SANDHAM

Direct numerical simulations (DNS) of laminar separation bubbles on a NACA-0012 airfoil at Rec=5×104 and incidence 5° are presented. Initially volume forcing is introduced in order to promote transition to turbulence. After obtaining sufficient data from this forced case, the explicitly added disturbances are removed and the simulation run further. With no forcing the turbulence is observed to self-sustain, with increased turbulence intensity in the reattachment region. A comparison of the forced and unforced cases shows that the forcing improves the aerodynamic performance whilst requiring little energy input. Classical linear stability analysis is performed upon the time-averaged flow field; however no absolute instability is observed that could explain the presence of self-sustaining turbulence. Finally, a series of simplified DNS are presented that illustrate a three-dimensional absolute instability of the two-dimensional vortex shedding that occurs naturally. Three-dimensional perturbations are amplified in the braid region of developing vortices, and subsequently convected upstream by local regions of reverse flow, within which the upstream velocity magnitude greatly exceeds that of the time-average. The perturbations are convected into the braid region of the next developing vortex, where they are amplified further, hence the cycle repeats with increasing amplitude. The fact that this transition process is independent of upstream disturbances has implications for modelling separation bubbles.


2019 ◽  
Vol 870 ◽  
pp. 680-697
Author(s):  
Dominik K. Puckert ◽  
Ulrich Rist

The interaction of disturbance modes behind an isolated cylindrical roughness element in a laminar boundary layer is investigated by means of hot-film anemometry and particle image velocimetry in a low-turbulence laminar water channel. Both sinuous and varicose disturbance modes are found in the wake of a roughness with unit aspect ratio (diameter/height $=$ 1). Interestingly, the frequency of the varicose mode synchronizes with the first harmonic of the sinuous mode when the critical Reynolds number from three-dimensional global linear stability theory is exceeded. The coupled motion of sinuous and varicose modes is explained by frequency lock-in. This mechanism is of great importance in many aspects of nature, but has not yet received sufficient attention in the field of boundary-layer theory. A Fourier mode decomposition provides detailed analyses of sinuous and varicose modes. The observation is confirmed by a second experiment with the same aspect ratio at a different position in the laminar boundary layer. When the aspect ratio is increased, the flow is fully governed by the varicose mode. Thus, no frequency lock-in can be observed in this case. The significance of this work is to explain how sinuous and varicose modes can co-exist behind a roughness and to propose a mechanism which is well established in physics but not encountered often in boundary-layer theory.


Sign in / Sign up

Export Citation Format

Share Document