Signature of large-scale motions on turbulent/non-turbulent interface in boundary layers

2017 ◽  
Vol 819 ◽  
pp. 165-187 ◽  
Author(s):  
Jin Lee ◽  
Hyung Jin Sung ◽  
Tamer A. Zaki

The effect of large-scale motions (LSMs) on the turbulent/non-turbulent (T/NT) interface is examined in a turbulent boundary layer. Using flow fields from direct numerical simulation, the shape of the interface and near-interface statistics are evaluated conditional on the position of the LSM. The T/NT interface is identified using the vorticity magnitude and a streak detection algorithm is adopted to identify and track the LSMs. Two-point correlation and spectral analysis of variations in the interface height show that the spatial undulation of the interface is longer in streamwise wavelength than the boundary-layer thickness, and grows with the Reynolds number in a similar manner to the LSMs. The average variation in the interface height was evaluated conditional on the position of the LSMs. The result provides statistical evidence that the interface is locally modulated by the LSMs in both the streamwise and spanwise directions. The modulation is different when the coherent structure is high- versus low-speed motion: high-speed structures lead to a wedge-shaped deformation of the T/NT interface, which causes an anti-correlation between the angles of the interface and the internal shear layer. On the other hand, low-speed structures are correlated with crests in the interface. Finally, the sudden changes in turbulence statistics across the interface are in line with the changes in the population of low-speed structures, which consist of slower mean streamwise velocity and stronger turbulence than the high-speed counterparts.

1980 ◽  
Vol 98 (1) ◽  
pp. 65-95 ◽  
Author(s):  
M. Sokolov ◽  
A. K. M. F. Hussain ◽  
S. J. Kleis ◽  
Z. D. Husain

A three-dimensional ‘turbulent spot’ has been induced in the axisymmetric free mixing layer of a 12.7 cm diameter air jet by a spark generated at the nozzle boundary layer upstream of the exit. The spot coherent-structure signature, buried in the large-amplitude random fluctuating signal, has been educed at three downstream stations within the apparent self-preserving region of the mixing layer (i.e. x/D = 1.5, 3.0 and 4.5) at the jet exit speed of 20 ms−1. The eduction has been performed through digital phase averaging of the spot signature from 200 realizations. In order to reduce the effect of the turbulence-induced jitter on the phase average, individual filtered signal arrays were optimally time-aligned through an iterative process of cross-correlation of each realization with the ensemble average. Further signal enhancement was achieved through rejection of realizations requiring excessive time shifts for alignment. The number of iterations required and the fraction of realizations rejected progressively increase with the downstream distance and the radial position.The mixing-layer spot is a large-scale elongated structure spanning the entire width of the layer but does not appear to exhibit a self-similar shape. The dynamics of the mixing-layer spot and its eduction are more complicated than those of the boundary-layer spot. The spot initially moves downstream essentially at a uniform speed across the mixing layer, but further downstream it accelerates on the high-speed side and decelerates on the low-speed side. This paper discusses the data acquisition and processing techniques and the results based on the streamwise velocity signals. Phase average distributions of vorticity, pseudo-streamlines, coherent and background Reynolds stresses and further dynamics of the spot are presented in part 2 (Hussain, Kleis & Sokolov 1980).


2021 ◽  
Vol 932 ◽  
Author(s):  
C.I. Chan ◽  
R.C. Chin

Well resolved large-eddy simulation data are used to study the physical modulation effects of miniature vortex generators (MVGs) in a moderate Reynolds number zero pressure gradient turbulent boundary layer. Large-scale counter-rotating primary vortex pairs (PVPs) imposed by the MVG contribute to the formation of streamwise streaks by transporting high momentum fluids from the outer regions of the boundary layer towards the wall, giving rise to high-speed regions centred at the PVP. Consequently, low-speed regions are formed along the outer flank of the PVP, resulting in a pronounced alternating high- and low-speed flow pattern. The PVP also relates to regions with skin friction modification, where a local skin friction reduction of up to 15 % is obtained at the low-speed region, but the opposite situation is observed over the high-speed region. The MVG-induced flow feature is further investigated by spectral analysis of the triple decomposition velocity fluctuation. Pre-multiplied energy spectra of the streamwise MVG-induced velocity fluctuation reveal that the large-scale induced modes scale with the spanwise wavelength and the length of the MVG, but the energy peak is eventually repositioned to the size of the near-wall streaks in the streamwise direction. Analysis of the triple decomposition of the kinetic energy transport equations revealed the significance of the mean flow gradient in generating kinetic energy which sustains the secondary motion. There is also an energy transfer between the turbulent and MVG-induced kinetic energy independent of the mean flow.


2018 ◽  
Vol 856 ◽  
pp. 135-168 ◽  
Author(s):  
S. T. Salesky ◽  
W. Anderson

A number of recent studies have demonstrated the existence of so-called large- and very-large-scale motions (LSM, VLSM) that occur in the logarithmic region of inertia-dominated wall-bounded turbulent flows. These regions exhibit significant streamwise coherence, and have been shown to modulate the amplitude and frequency of small-scale inner-layer fluctuations in smooth-wall turbulent boundary layers. In contrast, the extent to which analogous modulation occurs in inertia-dominated flows subjected to convective thermal stratification (low Richardson number) and Coriolis forcing (low Rossby number), has not been considered. And yet, these parameter values encompass a wide range of important environmental flows. In this article, we present evidence of amplitude modulation (AM) phenomena in the unstably stratified (i.e. convective) atmospheric boundary layer, and link changes in AM to changes in the topology of coherent structures with increasing instability. We perform a suite of large eddy simulations spanning weakly ($-z_{i}/L=3.1$) to highly convective ($-z_{i}/L=1082$) conditions (where$-z_{i}/L$is the bulk stability parameter formed from the boundary-layer depth$z_{i}$and the Obukhov length $L$) to investigate how AM is affected by buoyancy. Results demonstrate that as unstable stratification increases, the inclination angle of surface layer structures (as determined from the two-point correlation of streamwise velocity) increases from$\unicode[STIX]{x1D6FE}\approx 15^{\circ }$for weakly convective conditions to nearly vertical for highly convective conditions. As$-z_{i}/L$increases, LSMs in the streamwise velocity field transition from long, linear updrafts (or horizontal convective rolls) to open cellular patterns, analogous to turbulent Rayleigh–Bénard convection. These changes in the instantaneous velocity field are accompanied by a shift in the outer peak in the streamwise and vertical velocity spectra to smaller dimensionless wavelengths until the energy is concentrated at a single peak. The decoupling procedure proposed by Mathiset al.(J. Fluid Mech., vol. 628, 2009a, pp. 311–337) is used to investigate the extent to which amplitude modulation of small-scale turbulence occurs due to large-scale streamwise and vertical velocity fluctuations. As the spatial attributes of flow structures change from streamwise to vertically dominated, modulation by the large-scale streamwise velocity decreases monotonically. However, the modulating influence of the large-scale vertical velocity remains significant across the stability range considered. We report, finally, that amplitude modulation correlations are insensitive to the computational mesh resolution for flows forced by shear, buoyancy and Coriolis accelerations.


1996 ◽  
Vol 326 ◽  
pp. 1-36 ◽  
Author(s):  
FréDÉRic Ducros, Pierre Comte ◽  
Marcel Lesieur

It is well known that subgrid models such as Smagorinsky's cannot be used for the spatially growing simulation of the transition to turbulence of flat-plate boundary layers, unless large-amplitude perturbations are introduced at the upstream boundary: they are over-dissipative, and the flow simulated remains laminar. This is also the case for the structure-function model (SF) of Métais & Lesieur (1992). In the present paper we present a sequel to this model, the filtered-structure-function (FSF) model. It consists of removing the large-scale fluctuations of the field before computing its second-order structure function. Analytical arguments confirm the superiority of the FSF model over the SF model for large-eddy simulations of weakly unstable transitional flows. The FSF model is therefore used for the simulation of a quasi-incompressible (M∞ = 0.5) boundary layer developing spatially over an adiabatic flat plate, with a low level of upstream forcing. With the minimal resolution 650 × 32 × 20 grid points covering a range of streamwise Reynolds numbers Rex1 ε [3.4 × 105, 1.1 × 106], transition is obtained for 80 hours of time-processing on a CRAY 2 (whereas DNS of the whole transition takes about ten times longer). Statistics of the LES are found to be in acceptable agreement with experiments and empirical laws, in the laminar, transitional and turbulent parts of the domain. The dynamics of low-pressure and high-vorticity distributions is examined during transition, with particular emphasis on the neighbourhood of the critical layer (defined here as the height of the fluid travelling at a speed equal to the phase speed of the incoming Tollmien–Schlichting waves). Evidence is given that a subharmonic-type secondary instability grows, followed by a purely spanwise (i.e. time-independent) mode which yields peak-and-valley splitting and transition to turbulence. In the turbulent region, flow visualizations and local instantaneous profiles are provided. They confirm the presence of low- and high-speed streaks at the wall, weak hairpins stretched by the flow and bursting events. It is found that most of the vorticity is produced in the spanwise direction, at the wall, below the high-speed streaks. Isosurfaces of eddy viscosity confirm that the FSF model does not perturb transition much, and acts mostly in the vicinity of the hairpins.


2009 ◽  
Vol 622 ◽  
pp. 33-62 ◽  
Author(s):  
R. A. HUMBLE ◽  
G. E. ELSINGA ◽  
F. SCARANO ◽  
B. W. van OUDHEUSDEN

An experimental study is carried out to investigate the three-dimensional instantaneous structure of an incident shock wave/turbulent boundary layer interaction at Mach 2.1 using tomographic particle image velocimetry. Large-scale coherent motions within the incoming boundary layer are observed, in the form of three-dimensional streamwise-elongated regions of relatively low- and high-speed fluid, similar to what has been reported in other supersonic boundary layers. Three-dimensional vortical structures are found to be associated with the low-speed regions, in a way that can be explained by the hairpin packet model. The instantaneous reflected shock wave pattern is observed to conform to the low- and high-speed regions as they enter the interaction, and its organization may be qualitatively decomposed into streamwise translation and spanwise rippling patterns, in agreement with what has been observed in direct numerical simulations. The results are used to construct a conceptual model of the three-dimensional unsteady flow organization of the interaction.


1994 ◽  
Vol 259 ◽  
pp. 345-373 ◽  
Author(s):  
ROY Y. Myose ◽  
Ron F. Blackwelder

The dynamics and interaction of turbulent-boundary-layer eddy structures was experimentally emulated. Counter-rotating streamwise vortices and low-speed streaks emulating turbulent-boundary-layer wall eddies were generated by a Görtler instability mechanism. Large-scale motions associated with the outer region of turbulent boundary layer were emulated with — ωzspanwise vortical eddies shed by a periodic non-sinusoidal oscillation of an airfoil. The scales of the resulting eddy structures were comparable to a moderate-Reynolds-number turbulent boundary layer. Results show that the emulated wall-eddy breakdown was triggered by streamwise acceleration associated with the outer region of turbulent boundary layer. This breakdown involved violent mixing between low-speed fluid from the wall eddy and accelerated fluid associated with the outer structure. Although wall eddies can break down autonomously, the presence of and interaction with outer-region — ωzeddies hastened their breakdown. Increasing the — ωzeddy strength resulted in further hastening of the breakdown. Conversely, + ωzeddies were found to delay wall-eddy breakdown locally, with further delays resulting from stronger + ωzeddies. This suggests that the outer region of turbulent boundary layers plays a role in the bursting process.


2020 ◽  
Author(s):  
Benjamin Schumacher ◽  
Marwan Katurji ◽  
Jiawei Zhang

<p>The evolution of micrometeorological measurements has been recently manifested by developments in methodological and analytical techniques using spatial surface brightness temperature captured by infrared cameras (Schumacher et al. 2019, Katurji and Zawar-Reza 2016). The Thermal Image Velocimetry (TIV) method can now produce accurate 2D advection-velocities using high speed (>20Hz) infrared imagery (Inagaki 2013, Schumacher 2019). However, to further develop TIV methods and achieve a novel micrometeorological measurement technique, all scales of motion within the boundary layer need to be captured.</p><p>Spatial observations of multi-frequency and multi-scale temperature perturbations are a result from the turbulent interaction of the overlying atmosphere and the surface. However, these surface signatures are connected to the larger scales of the atmospheric boundary layer (McNaughton 2002, Träumner 2015). When longer periods (a few hours to a few days) of spatial surface brightness temperatures are observed, the larger scale information needs to be accounted for to build a comprehensive understanding of surface-atmospheric spatial turbulent interactions. Additionally, the time-frequency decomposition of brightness temperature perturbations shows longer periods of 4-15 minutes superimposed over shorter periods of ~ 4–30 seconds. This suggests that that boundary layer dynamic scales (of longer periods) can influence brightness temperature perturbations on the local turbulent scale. An accurate TIV algorithm needs to account for all scales of motion when analysing the time-space variability of locally observed spatial brightness temperature patterns.</p><p>To analyse these propositions temporally high resolved geostationary satellite infrared data from the Himawari 8 satellite was compared to near-surface and high speed (20 Hz) measured air and brightness temperature using thermocouple measurements and infrared cameras. The satellite provides a temporal resolution of 10-minutes and a horizontal resolution of 2 by 2 km per pixel and therefore captures the atmospheric meso γ and micro α scale which signals are usually active for ~10 minutes to < 12 hours. Moreover, the Himawari 8 brightness temperature was used to create the near-surface mean velocity field using TIV. Afterwards, the velocity field was compared to the in-situ measured wind velocity over several days during January 2019.</p><p>The results show that the atmospheric forcing from the micro α scale to lower atmospheric scales has a major impact on the near-surface temperature over several minutes. A significant (p-value: 0.02) positive covariance between the Himawari 8 measurement and the local measured temperature 1.5 cm above the ground on a 10 minute average, specifically concerning cooling and heating patterns, has been found.</p><p>Further analysis demonstrates that the retrieved near-surface 2-D velocity field calculated from the Himawari 8 brightness temperature perturbations is correctly representing the mean velocity. This finding allows the classification of meso-scale atmospheric forcing and its direct connection to local scale turbulent 2-D velocity measurements. This extends the TIV algorithm by a multi-scale component which allows to address inter-scale boundary layer analysis from a new point of view. In respect to the current findings a new experiment will focus on the repeated induced local velocity patterns from large scale forcing which will be measured through the surface brightness temperature.</p>


2014 ◽  
Vol 748 ◽  
pp. 848-878 ◽  
Author(s):  
Pramod K. Subbareddy ◽  
Matthew D. Bartkowicz ◽  
Graham V. Candler

AbstractWe study the transition of a Mach 6 laminar boundary layer due to an isolated cylindrical roughness element using large-scale direct numerical simulations (DNS). Three flow conditions, corresponding to experiments conducted at the Purdue Mach 6 quiet wind tunnel are simulated. Solutions are obtained using a high-order, low-dissipation scheme for the convection terms in the Navier–Stokes equations. The lowest Reynolds number ($Re$) case is steady, whereas the two higher $Re$ cases break down to a quasi-turbulent state. Statistics from the highest $Re$ case show the presence of a wedge of fully developed turbulent flow towards the end of the domain. The simulations do not employ forcing of any kind, apart from the roughness element itself, and the results suggest a self-sustaining mechanism that causes the flow to transition at a sufficiently large Reynolds number. Statistics, including spectra, are compared with available experimental data. Visualizations of the flow explore the dominant and dynamically significant flow structures: the upstream shock system, the horseshoe vortices formed in the upstream separated boundary layer and the shear layer that separates from the top and sides of the cylindrical roughness element. Streamwise and spanwise planes of data were used to perform a dynamic mode decomposition (DMD) (Rowley et al., J. Fluid Mech., vol. 641, 2009, pp. 115–127; Schmid, J. Fluid Mech., vol. 656, 2010, pp. 5–28).


Author(s):  
J. Stephen Hu ◽  
Jian Sheng ◽  
Michele Guala ◽  
Leonardo Chamorro

The focus of this paper is to characterize the upstream wake of a three bladed Horizontal Axis Wind Turbine (HAWT) and its interaction with the native structures within a turbulent boundary layer (TBL). The overarching question is the most prevailing length and time scales of coherent structures that would interact with a HAWT and how they would be affected. The implications include wall flow and structure interaction and flow induced noise generation in large scale turbo machineries. The experiments are performed on a turbine that has a 0.128 m rotor diameter, a hub height of 0.104 m and a tip speed ratio of 4. The HAWT model is placed in a large scale wind tunnel in a boundary layer with a thickness δ of ∼0.6 m. The boundary layer is generated by a 60 mm picket fence trip and developed over a smooth wall under thermally neutral conditions. Measurements are performed under ReD of 4 × 105 and 6 × 105. Both turbine geometries and flow conditions are scaled from operating conditions in the field. High speed Particle Image Velocimetry (PIV), turbine voltage output, and angular velocity measurements are conducted simultaneously, by which one could relate the upwind flow structures with the power output of the turbine. High speed PIV offer details in spatial and temporal characteristics of the impinging flow structures, whilst the voltage anemometer and tachometer provide instantaneous measurement of angular velocity of the turbine. PIV measurements are taken at a rate of 1500 image pairs per second with a 100 μs delay between laser pulses. Each sample area is 0.15 × 0.15 cm. Two locations up to two rotor diameters upwind are measured. Instantaneous voltage is taken at a sampling rate of 30 kHz and a sampling time of 60s to ensure sufficient temporal resolution and coverage. Ongoing analysis using conditional averaging based on extreme power output events will provide insights in assessing a HAWT performance in unsteady flow conditions.


2009 ◽  
Vol 48 (10) ◽  
pp. 2069-2085 ◽  
Author(s):  
Cesar Azorin-Molina ◽  
Bernadette H. Connell ◽  
Rafael Baena-Calatrava

Abstract The aim of this study was to identify clear air boundaries and to obtain spatial distribution of convective areas associated with the sea breeze over the Iberian Mediterranean zone and the isle of Mallorca, both in Spain. Daytime Advanced Very High Resolution Radiometer (AVHRR) data from National Oceanic and Atmospheric Administration (NOAA) polar-orbiting satellites were collected for May–October 2004. A cloud detection algorithm was used to identify clouds to derive daytime sea-breeze cloud frequency composites over land. The high-resolution composites aided in identifying the location of five preferential sea-breeze convergence zones (SBCZ) in relation to the shape of coastline and orographic effects. Additionally, eight regimes were designated using mean boundary layer wind speed and direction to provide statistics about the effect of prevailing large-scale flows on sea-breeze convection over the five SBCZ. The offshore SW to W and the NW to N regimes were characterized by high cloud frequencies parallel to the coast. Small differences in mean cloud frequency values from morning to afternoon composites were detected with these regimes because sea-breeze fronts tended to form early and persist into the afternoon. Just the opposite occurred under the onshore NE to E and SE to S regimes. It was found that light to moderate (≤5.1 m s−1) winds aloft result in more clouds at the leading edge of sea breezes. In contrast, strong synoptic-scale (>5.1 m s−1) flows weaken boundary layer convergence. The results from this satellite meteorology study could have practical applications for many people including those that forecast the weather and those that use the forecast for making decisions related to energy use, fishing, recreation, or agriculture activities, as well as for estimating pollution or issuing warnings for heavy rain or flash flooding.


Sign in / Sign up

Export Citation Format

Share Document