Next Article in Journal
Numerical Study on Effects of Air Return Height on Performance of an Underfloor Air Distribution System for Heating and Cooling
Next Article in Special Issue
Numerical Study of the Effects of the Jet Fan Speed, Heat Release Rate and Aspect Ratio on Smoke Movement in Tunnel Fires
Previous Article in Journal
Optimal Inertia Reserve and Inertia Control Strategy for Wind Farms
Previous Article in Special Issue
Surrogate Model with a Deep Neural Network to Evaluate Gas–Liquid Flow in a Horizontal Pipe
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Simulation of Coherent Structures in the Turbulent Boundary Layer under Different Stability Conditions

1
Key Lab of Smart Prevention and Mitigation for Civil Engineering Disasters of the Ministry of Industry and Information, Harbin Institute of Technology, Harbin 150090, China
2
Key Lab of Structures Dynamic Behavior and Control of the Ministry of Education, Harbin Institute of Technology, Harbin 150090, China
3
School of Civil Engineering, Harbin Institute of Technology, Harbin 150090, China
*
Author to whom correspondence should be addressed.
Energies 2020, 13(5), 1068; https://doi.org/10.3390/en13051068
Submission received: 7 February 2020 / Revised: 25 February 2020 / Accepted: 25 February 2020 / Published: 1 March 2020
(This article belongs to the Special Issue Engineering Fluid Dynamics 2019-2020)

Abstract

:
Coherent structures in the turbulent boundary layer were investigated under different stability conditions. Qualitative analyses of the flow field, spatial correlation coefficient field and pre-multiplied wind velocity spectrum showed that the dominant turbulent eddy structure changed from small-scale motions to large- and very-large-scale motions and then to thermal plumes as the stability changed from strong stable to neutral and then to strong unstable. A quantitative analysis of the size characteristics of the three-dimensional turbulent eddy structure based on the spatial correlation coefficient field showed that under near-neutral stability, the streamwise, wall-normal and spanwise extents remained constant at approximately 0.3 δ , 0.1 δ and 0.2 δ ( δ , boundary layer height), respectively, while for other conditions, the extent in each direction varied in a log-linear manner with stability; only the spanwise extent under stable conditions was also independent of stability. The peak wavenumber of the pre-multiplied wind velocity spectrum moves towards small values from stable conditions to neutral condition and then to unstable conditions; thus, for the wind velocity spectrum, another form is needed that takes account the effects of the stability condition.

Graphical Abstract

1. Introduction

Over flat and homogeneous terrain, turbulence is a very important factor affecting turbulent motions near the ground. Consequently, the structure of the turbulent boundary layer depends mainly on vertical turbulence transport which has a great influence on mass, momentum, heat, water vapor and energy exchanges [1] with turbulent coherent structures playing an important role [2,3]. Field measurements, laboratory experiments and numerical simulations of the turbulent boundary layer have revealed the existence of a wide variety of coherent structures such as low-speed streaks, hairpin vortices, large-scale motions and very large-scale motions in canonical wall-bounded flows including pipe flows, channel flows and turbulent boundary layer flows [2,4,5,6,7,8,9,10]. Actually, coherent structures are characterized by a specific scale and lifetime, the recurrence of characteristic and flow-specific organization and regions of concentrated vorticity. Specifically, Theodorsen [4] proposed the existence of hairpin vortices, and such vortices, together with streamwise packets, have been suggested to be the primary elements forming coherent structures in the neutral turbulent boundary logarithmic layer under experimental conditions [6]. However, most studies have focused on neutral conditions. Under neutral conditions, turbulent structures observed in the field are similar to those found in experiments. Furthermore, Marusic et al. [3] indicated that turbulence and even the pressure field in the turbulent boundary layer are also similar to those in canonical wall-bounded flows, and they suggested that the reason for this is the presence of coherent structures.
The prevalent view regarding the formation of coherent structures in the turbulent surface layer is the “bottom-up” mechanism, in which active structures emerge mainly at the surface and then grow up into the outer layer [8,11,12,13]. As is well known, the turbulent boundary layer is usually classified as being stable, neutral or unstable on the basis of the dominant mechanism of turbulence generation and stability [14,15]. Differences in stability will affect energy transport and fluxes between the near-surface layer and upper free layer, through the formation and development of coherent structures [16,17,18,19,20]. The stability depends on its ability to resist three-dimensional turbulent motions. In stable conditions, even slight three-dimensional turbulent motions are difficult or even impossible. By contrast, in unstable conditions, slight three-dimensional turbulent motions tend to become stronger, leading to turbulence and convective activity. Li and Bou-Zeid [18] concluded that coherent structures in the atmosphere must be different from those found in experiments, because in the turbulent boundary layer they can be affected by the surface buoyancy flux under stable (nocturnal) or unstable (daytime) conditions which are not taken into account in most experiments. Specifically, unstable turbulent boundary layers include three kinds of coherent structures: under slightly unstable conditions, large-scale horizontal roll vortices occur, owing to a combination of shear and buoyancy [21,22,23,24]; under moderately unstable conditions, individual hairpins arise from the ground, leading to an increase in the inclination angle of the vortex packets [25,26,27,28]; and under strongly unstable conditions, plume-shaped structures appear, together with streamwise elongated structures, and the magnitude of the vertical turbulent momentum flux increases [17,25,29,30]. Stable turbulent boundary layers can generally be divided into two kinds: weakly stable and strongly stable [31]. Under weakly stable conditions, because turbulence is nearly continuous and quasi-steady, Monin–Obukhov similarity theory applies [32]. Under strongly stable conditions, small eddy motions dominate in the boundary layer, and the layer is much shallower than those under neutral or unstable conditions; however, our understanding of this is still quite limited [31,33,34] partly because of the intrinsic complexity of fluid dynamical phenomena such as intermittency and Kelvin–Helmholtz instability [31,33,35].
In the present study, to obtain a better understanding of turbulent boundary layer flow under different stability conditions, a large eddy simulation (LES) approach was adopted. The LES technique was proposed for meteorological applications by Deardorff [36,37]; since then, it has become an indispensable tool to study the turbulent boundary layer [38,39,40,41]. In LES, large-scale structures are resolved explicitly, while the effect of small-scale structures to large-scale structures is not resolved but is parameterized using a subgrid-scale (SGS) model. Thus, because the unstable turbulent boundary layer is dominated by large-scale structures, LES is an appropriate approach [42,43,44]. Large eddy simulation has also been used successfully for neutral conditions [17,20,42,44,45]. However, simulation of the stable turbulent boundary layer faces a number of challenges. These arise because the way in which the radiative cooling of the ground leads to stable stratification modifying the flow turbulence is still not completely understood [46], and because the characteristic size of eddies decreases with increasing stability. The application of LES to stable conditions has, therefore, been mostly confined to weakly or moderately stable cases [47,48], and this approach has not been considered appropriate for strongly stable conditions [34,49]. However, Churchfield et al. [50] have suggested an approach that is capable of simulating complex terrain and stable conditions using wall model LES (WMLES), and it is this approach that is adopted in the present study.
Turbulent boundary layers have a complex structure and involve some phenomena that are still not completely understood, especially at high Reynolds numbers [9]. Therefore, the objective of this study was to investigate the characteristics of the turbulent boundary layer under different stability conditions and to quantify the variations in the shapes and sizes of eddy structures with stability. The LES model and simulation setup are described in Section 2, and the model is validated in Section 3. Simulation results are discussed in Section 4. Finally, conclusions are given in Section 5.

2. Simulation Details

The governing equations of the LES model are:
u ¯ i t + x j ( u ¯ j u ¯ i ) = x i p ¯ x j ( τ i j D ) + f i T
u ¯ i x i = 0
For potential temperature transport, the governing equation is:
θ ¯ t + x j ( u ¯ j θ ¯ ) = x j ( q j )
where x i represent Cartesian coordinates, u i represents the corresponding velocity components, p represents the pressure which, divided by the fluid density and τ i j D , represents the SGS stress, and f i T represents the density-normalized forces. The Coriolis force was not considered, the same as in Fang and Porté-Agel [20]. A bar on a variable indicates the mean grid-filtered quantity. θ ¯ is the resolved potential temperature, and q j represents the flux of temperature due to the viscous and sub-filter-scale effects. The details of how the SGS momentum flux τ i j D was modelled is shown in Ren et al. [51]. The heat flux q j at the surface was modelled using a similar method [38]:
q j = [ 0 0 q 3 t o t ]
where, in the simulation, the total average temperature flux is specified, and the fluctuating temperature flux q 3 t o t is created by the wall model.
For numerical discretization, the finite volume method was used; for the viscous and convective terms, a second-order central difference scheme was adopted; and for the unsteady term, a second-order backward scheme was adopted.
Schumann’s model [52] was used as the wall model to model the near wall shear stress, and Moeng’s model [38] was adopted for the total temperature flux at the surface. The subgrid stress (SGS) model is the one-equation eddy viscosity model [53]; this is because the one-equation eddy viscosity model does not require the assumption in the algebraic eddy viscosity models, that is, local balance between the production and dissipation of SGS energy. The SGS stress can be modelled as:
τ i j , s g s D = 2 3 k s g s δ i j 2 υ s g s S ¯ i j
where τ i j , s g s D is the SGS stress part of the total stress, τ i j D . S ¯ i j is the resolved-scale strain rate tensor, δ i j is the Kronecker delta function, k s g s is the SGS kinetic energy; and υ s g s is the SGS stress eddy viscosity.
Accounting for the historic effect of k s g s , because of diffusion, production and dissipation, a transportation equation was derived:
k s g s t + u ¯ j k s g s x j = τ i j , s g s D S ¯ i j C ε k s g s 3 / 2 Δ + x j ( υ s g s k s g s x j )
where Δ represents the filter scale, and C k and C ε are the model constants.
In the present study, the size of the computational domain was 5 × 1 × 0.3 km for all cases. When the hexahedral mesh was applied, the mesh grid size was 5 m in each direction from the ground surface to 0.2 km, and above 0.2 km, the mesh grid size was 10 m. For the streamwise and spanwise directions, periodic boundary conditions were applied; this agrees well with Fang and Porté-Agel [20] of LES of turbulent boundary layer flow and with Lee et al. [54] of DNS of channel flow; for the upper boundary, the slip boundary condition was adopted; and for the solid wall, the wall model was used. For the initial conditions, the logarithmic mean wind speed profile was used in all cases (gradient wind height of 250 m and a gradient wind speed of 8 m/s).
The computation timestep was 0.5 s which ensures the stability and viscous stability by a given Courant number. First, to obtain a fully developed flow field and quasi-steady state, a spin-up simulation was performed. For the statistical analysis, to ensure statistical convergence, over 40,000 timesteps were simulated. For the perturbation component, the initial peak was set at the height of 0.015% of the computational domain, and for the horizontal direction, the maximum perturbation was 0.25 m/s; at the beginning of the simulation, for the streamwise and spanwise direction, the number of periods of the perturbations were 4.0 and 20.0, respectively.
The aim of this numerical simulation was to calculate the forms of coherent structures under specifically assumed stability conditions and to quantify the variations in the shapes and sizes of eddy structures with stability. The simulation was conducted within a pre-defined stability range. As is well known, stability refers to the air’s tendency to either resist vertical movement (stability) or to rise and create storms (instability) which can be expressed in terms of a dimensionless stability parameter, z / L , where z represents the height above the surface ground, and L represents the Obukhov length, given by L = u 3 T v / ( k g Q v 0 ) , where k represents the von Kármán constant, u represents the friction velocity, Q v 0 is the kinematic virtual temperature flux at the surface, T v is a reference virtual temperature, and g is the gravitational acceleration. L expresses the relative roles of shear and buoyancy in the production/consumption of turbulence kinetic energy with positive and negative values indicating stable and unstable conditions, respectively, and it approaches zero in the limit of neutral stratification. For z / L , positive and negative values indicate unstable and stable conditions. In this study, 17 different turbulent boundary layer cases were analysed, from case 1 to 17 which range from stable to unstable, with z / L = −98.27 (case 1), −16.19 (case 2), −2.29 (case 3), −1.15 (case 4), −0.84 (case 5), −0.63 (case 6), −0.35 (case 7), −0.05 (case 8), 0 (case 9), 0.54 (case 10), 1.01 (case 11), 1.79 (case 12), 3.16 (case 13), 4.48 (case 14), 8.13 (case 15), 11.31 (case 16), and 31.46 (case 17). It should be noted that to represent the average stability of the turbulent boundary layer, the z used for calculating z / L here is the boundary layer height δ ; this is different from Figure 1 and Figure 2 in which z represents the level height above the ground.

3. Validation

The basic statistical properties of the velocity under neutral stability conditions were compared with those from previous studies, and good agreement was found when dimensioned both by inner and outer variables; the details are shown in previous studies by Ren et al. [51] and Ren et al. [55].
The ratio of the heat fluxes in the horizontal and vertical directions under unstable conditions is shown in Figure 2. Wyngaard and Cote [56] proposed the following formula for this ratio:
u θ ¯ w θ = a φ m φ h
where φ m = ( k z / u * ) ( u / z ) represents the non-dimensional wind shear function, u is the surface friction velocity, φ h = ( k z / u * ) ( θ / z ) is the non-dimensional potential temperature gradient function, and a is a constant. In the present study, φ m = ( 1 28 z / L ) 1 / 4 , φ h = ( 1 20 z / L ) 1 / 2 [57,58], and a = 5 [57,58,59]. In Figure 1, the ratio of the horizontal and vertical heat fluxes decreases with increasing z / L away from the near-surface region, and u θ ¯ / w θ ¯ tends towards 0 which indicates that at higher elevations the shear stress is very small compared with the thermal plume force. The simulation results show good agreement with theory under strongly unstable conditions, but under moderately or slightly unstable conditions, they deviate from the theoretical values. This deviation occurs because under these conditions, the shear stress driving the horizontal motions near the surface is greater than under strongly unstable conditions, while the vertical transport is weaker than under strongly unstable conditions. Therefore, Equation (7) is appropriate for strongly unstable conditions only, and perhaps, after the forms of a, φ m and φ h are changed in the present study, it will also prove to be appropriate for moderately or slightly unstable conditions. Another phenomenon revealed in Figure 1 is a shift in the peak of w θ ¯ / w θ ¯ towards larger z / L ; this is because as conditions become more unstable, thermal plumes reach higher. Similarly, the peak of u θ ¯ / w θ ¯ shifts to smaller z / L as the stability increases under stable conditions. This suggests that as conditions become more unstable, z / L becomes larger which is related to the upper location of thermal plumes.
Figure 2 shows the normalized standard deviation of the vertical velocity component, σ w / u * , which is important for describing stability. The curves in Figure 2 represent the formula proposed by Panofsky et al. [59] in the case of the unstable condition:
σ w / u * = 1.3 ( 1 3 z / L ) 1 / 3
where u * is the friction velocity, L represents the Obukhov length, and z is the height above the ground. It can be seen that σ w / u * remains constant under slightly or moderately unstable conditions but increases with increasing height under strongly unstable conditions. The simulation results show good agreement with Equation (8) which means that σ w / u * obeys Monin–Obukhov similarity theory under unstable conditions, fitting well with the 1/3 power law. Under neutral conditions, σ w / u * is constant, independent of height and equal to 1.32 (shown in Figure 2b); indeed, a large number of observations show that under neutral conditions, σ w / u * should be constant [58,60,61,62,63,64]. In particular, Bowne and Ball [60] and Counihan [61] both concluded that this constant value is 1.3. Roth [64] summarized the values of different researchers, which ranged from 1.1 to 1.7. When the turbulent boundary layer is under stable conditions, the turbulent motions are weakly connected to the ground, and the height is no longer a similarity parameter of this motion, except in the region near the ground. As can be seen from Figure 2c, the normalized standard deviation of the vertical velocity component behaves as expected for stable conditions. Also, similar to the ratio of horizontal and vertical heat fluxes, the peak of σ w / u * shifts to smaller z / L as stability increases, except for neutral conditions; this is related to the upper location of thermal plumes.

4. Results

The above verification analysis has illustrated the feasibility of the WMLES approach and the accuracy of simulation of turbulent boundary layer flow under different stability conditions. In the following, we investigate the characteristics of turbulent boundary layers under different stability conditions and quantify the variations in the shape and size of eddy structures with stability.

4.1. Flow Field

We now turn to an analysis of the flow field with the aim of providing a more intuitive view of its basic characteristics and those of the turbulent eddy scale. We consider a horizontal plane (whole domain) at the height of 0 . 2 δ . Figure 3 shows the flow field for all 17 cases. On the whole, under stable conditions, the dominant turbulent eddy structure was small which is consistent with many previous results [31,33,34]. This is because, under stable conditions, the ground temperature is lower than that of the near-surface flow, and the vertical development of the turbulent eddy structure is limited which suppresses the development of three-dimensional turbulent eddy structure thus small-scale motion is prevalent. The dominant turbulent eddy structure scale increases as the stability decreases, and for the weakly stable condition in case 8, very-large-scale motion structures [8,10,65] even occur, just as under neutral conditions. It is obvious under neutral conditions that both large-scale and very-large-scale motion structures are visualized, based on the streamwise-elongated volumes associated with negative streamwise velocity, which is consistent with the results of Dennis and Nickels [66,67], Lee and Sung [65], Baltzer et al. [68] and Lee et al. [54]. Under slightly and moderately unstable conditions, obvious large-scale horizontal roll vortices form which are usually attributed to the combined effects of shear and buoyancy [21,22,24]. These can be seen as corresponding to the extended turbulent eddy structure that appears under neutral stability conditions. Under strongly unstable conditions, the effects of buoyancy dominate those of shear, so there is an obvious distinction between low-velocity and high-velocity regions which extends over the whole spanwise direction. This is consistent with the experimental results of Hommema and Adrian [25].
In Figure 3, going from the strongly stable condition to the neutral stable condition to the strongly unstable condition, some interesting phenomena can be seen. First, the form and shape of the dominant turbulent eddy structure changed from small-scale motions to large-scale and very-large-scale motions and then to regional motions as the turbulent boundary layer became unstable. Second, the dominant size of the turbulent eddy structure varied as the stability changed; a quantitative analysis of this is presented in Section 4.2 based on the three-dimensional spatial correlation coefficient field.

4.2. Spatial Correlation Coefficient Field

To quantify the size characteristics of the three-dimensional turbulent eddy structure under different stability conditions, a flow field correlation analysis was performed with reference location ( x , y , z ) = ( 2500 m , 500 m , 35 m ) . Figure 4 presents contours of the spatial correlation coefficient for all 17 cases. The plots in the left column are for the x z plane, and those in the right column are for the y z plane; thus, it is possible to visualize the three-dimensional turbulent eddy structure. The conclusion is similar to that from the flow field analysis, namely, that with increasing stability, the dominant size of the turbulent eddy decreases. This clearly shows that under weakly stable, neutral and slightly and moderately unstable conditions, very-large-scale motion structures appeared, while strongly unstable conditions led to thermal plumes. Meanwhile, an interesting phenomenon was that the spanwise size of the dominant vortices remained almost constant, except under strongly unstable conditions which may be due to the fact that the effect of the shear force was mainly in the streamwise direction and that of the buoyancy force was mainly in the wall-normal direction; however, under strongly unstable conditions, the temperature of the ground was higher than that of the adjacent near-surface air, and so the buoyancy force had a positive effect on the upward motion of air near the ground in the whole spanwise direction. The spanwise and wall-normal velocity correlations were different from the streamwise velocity correlation. For the spanwise velocity correlation under stable conditions, the spanwise size of the correlation structure was smaller than the streamwise velocity correlation; moreover, the degree of tilt was smaller than for the streamwise velocity correlation. Under unstable conditions, even light unstable conditions, the plume-shaped structure occurred, meaning that this structure form occurred earlier than the streamwise velocity correlation. For the wall-normal velocity correlation, from stable to unstable conditions, the spanwise and wall-normal sizes of the correlation structure became wide and deep which is very different to the other two direction velocity correlations. These results are not shown here.
Figure 5 shows the size of the turbulent eddy structure under different stability conditions at a height of 0 . 2 δ . The distances L x u u , L y u u and L z u u are defined in Christensen and Wu [69] and Volino et al. [70] as twice the distances from the self-correlation peak to the most distant locations on the 0.5 contour of the two-point correlation of the fluctuating velocity in the downstream, spanwise and wall-normal directions, respectively. We can see that under strongly and moderately stable conditions, the streamwise extent L x u u increased in a log-linear manner as z / L increased, and that it behaved similarly under strongly unstable conditions. However, under near-neutral stability conditions, L x u u was independent of z / L and remained constant at approximately 0.3 δ . Similar behaviour was found for the wall-normal extent L z u u which remained constant at approximately 0.1 δ . From strongly stable through neutral to slightly unstable conditions, the spanwise extent L y u u remained constant at approximately 0.2 δ , while under strongly unstable conditions, L y u u also increased in a log-linear manner with increasing z / L . This corresponds to the behaviour of the spatial correlation coefficient field.

4.3. Pre-Multiplied Spectral Analysis

Pre-multiplied spectral analysis is a popular approach for determining the distribution characteristics of vortices of different scales in canonical wall-bounded flows [8,10,71,72]. Figure 6 shows the multipoint averaged [73] pre-multiplied wind velocity spectra for all 17 cases. As we know, the most energetic point in the pre-multiplied wind velocity spectrum represents the dominant turbulent eddy structure. Therefore, under stable conditions, the peak wavenumber of the wind velocity spectrum is located in the high-wavenumber region which indicates that the dominant vortices are small. With decreasing stability, the peak wavenumber moved to smaller values which means that the size of the dominant turbulent eddy structure increased and very-large-scale motion structures occurred under the slightly stable conditions of case 8. Here, it should be noted that very-large-scale motion structures are those with turbulent eddy size larger than 3 δ . Under near-neutral stability conditions (slightly stable, neutral and slightly unstable conditions), the stationary stage of the pre-multiplied wind velocity spectrum has a large width, because the stationary stage corresponds to the very-large-scale motion structure [6,8,10,12]. Under strongly unstable conditions, the peak wavenumber is located in the small wavenumber region which means that the dominant turbulent eddy structure is very large, and the shape of the wind velocity spectrum is similar to that of a typhoon wind velocity spectrum near the eye-wall region. This may be because the eye-wall region of a typhoon is also characterized by strongly convective conditions, accompanied by a strong updraft. For spanwise and wall-normal pre-multiplied wind velocity spectra, the maximum energy points are consistent with the streamwise pre-multiplied wind velocity spectrum, and as the instability increases, it moves to a small wavenumber. The difference is that the steady stage is more obvious in the streamwise pre-multiplied wind velocity spectrum than the other two directions, and that the spanwise and wall-normal pre-multiplied wind velocity spectra have obvious cross-over phenomenon at different heights, which does not occur for streamwise pre-multiplied wind velocity spectrum. These results are not shown here.
It can be seen from Figure 6 that as the stability conditions changed from stable to neutral to unstable, the peak wavenumber of the pre-multiplied wind velocity spectrum moved towards smaller values which is consistent with the above analyses of the flow field and the spatial correlation coefficient field, and the stationary stage of the wind velocity spectrum existed only under near-neutral conditions. Under strongly stable or unstable conditions, the small-scale and large-scale turbulent eddy structures were the dominant structures, respectively, and only one dominant peak existed in the pre-multiplied wind velocity spectrum.
The canonical wind velocity spectral styles are based on the neutral stability condition, but in this paper, we found that the different stability conditions gave different wind velocity spectral forms. Therefore, in future studies, we should try to find a stability parameter that represents the differences in the wind velocity spectra and allows us to obtain a universal wind velocity spectrum that is suitable for all stability conditions.

5. Conclusions

We investigated the coherent structures in the turbulent boundary layer under different stability conditions and performed qualitative and quantitative analyses of the shape and size of eddy structures.
First, the WMLES approach and the basic variables of the simulation results were validated by comparing these results with theoretical predictions, observed data, experimental data and other simulation results.
Then, based on a flow field analysis, it was found that the form and shape of the dominant turbulent eddy structure changed from small-scale motions to large-scale and very-large-scale motions and then to thermal plumes as the stability condition changed from strong stable conditions to neutral stability conditions and then, to strongly unstable conditions. A spatial correlation analysis was performed, and the results led to similar conclusions as those obtained from the flow field analysis. An interesting phenomenon was also found, namely, that the spanwise size of the dominant vortices remained almost constant, except under strongly unstable conditions which may be because the main effect of shear is in the streamwise direction while that of buoyancy is in the wall-normal direction. From an analysis of the pre-multiplied wind velocity spectrum, it was found that as the stability condition goes from stable to neutral to unstable, the peak wavenumber of the pre-multiplied wind velocity spectrum moves towards smaller values, and the stationary stage of wind velocity spectrum only exists under near-neutral stability, whereas under strongly stable or unstable conditions, the small- and large-scale turbulent eddy structures have an overwhelming advantage, so only one dominant peak exists in the pre-multiplied wind velocity spectrum. In future work, we need to consider a new style for the wind velocity spectrum that takes account of the effects of stability.
Finally, a quantitative analysis of the size characteristics of the three-dimensional turbulent eddy structure under different stability conditions based on the spatial correlation coefficient field showed that for near-neutral stability, the streamwise, wall-normal and spanwise extents remain constant at approximately 0.3 δ , 0.1 δ and 0.2 δ , respectively. For other conditions, it was found that the streamwise, wall-normal and spanwise extents increase in a log-linear manner with increasing stability, except for the spanwise extent under stable conditions which is independent of stability.

Author Contributions

Conceptualization, H.R., S.L., H.L. and J.O.; Methodology, H.R. and S.L.; Software, H.R.; Validation, H.R.; Formal Analysis, H.R.; Investigation, H.R. and S.L.; Resources, S.L.; Data Curation, H.R.; Writing–Original Draft Preparation, H.R.; Writing–Review and Editing, S.L.; Visualization, H.R.; Supervision, H.L. and J.O.; Project Administration, H.L.; Funding Acquisition, H.L. All authors have read and agreed to the published version of the manuscript.

Funding

This study was financially supported by the National Key Research and Development Program of China under grant no. 2018YFC0705605, NSFC under grant no. 51878230, the China Postdoctoral Science Foundations (2015M580266 and 2016T90295), and the Heilongjiang Postdoctoral Fund (LBH-Z15087).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Monin, A.S.; Obukhov, A.M. Basic laws of turbulent mixing in the surface layer of the atmosphere. Contrib. Geophys. Inst. Acad. Sci. USSR 1954, 24, 163–187. [Google Scholar]
  2. Robinson, S.K. Coherent motions in the turbulent boundary-layer. Annu. Rev. Fluid Mech. 1991, 23, 601–639. [Google Scholar] [CrossRef]
  3. Marusic, I.; Mathis, R.; Hutchins, N. Predictive model for wall-bounded turbulent flow. Science 2010, 329, 193–196. [Google Scholar] [CrossRef] [PubMed]
  4. Theodorsen, T. Mechanism of turbulence. In Proceedings of the Second Midwestern Conference on Fluid Mechanics, Ohio State University, Columbus, OH, USA, 17–19 March 1952; pp. 1–19. [Google Scholar]
  5. Kline, S.J.; Reynolds, W.C.; Schraub, F.A.; Runstadl, P.W. Structure of turbulent boundary layers. J. Fluid Mech. 1967, 30, 741–773. [Google Scholar] [CrossRef] [Green Version]
  6. Adrian, R.J. Hairpin vortex organization in wall turbulence. Phys Fluids 2007, 19, 041301. [Google Scholar] [CrossRef] [Green Version]
  7. Zhou, J.; Adrian, R.J.; Balachandar, S.; Kendall, T.M. Mechanisms for generating coherent packets of hairpin vortices in channel flow. J. Fluid Mech. 1999, 387, 353–396. [Google Scholar] [CrossRef]
  8. Kim, K.C.; Adrian, R.J. Very large-scale motion in the outer layer. Phys Fluids 1999, 11, 417–422. [Google Scholar] [CrossRef]
  9. Hunt, J.C.R.; Morrison, J.F. Eddy structure in turbulent boundary layers. Eur. J. Mech. B 2000, 19, 673–694. [Google Scholar] [CrossRef]
  10. Wang, G.H.; Zheng, X.J. Very large scale motions in the atmospheric surface layer: A field investigation. J. Fluid Mech. 2016, 802, 464–489. [Google Scholar] [CrossRef]
  11. Adrian, R.J.; Meinhart, C.D.; Tomkins, C.D. Vortex organization in the outer region of the turbulent boundary layer. J. Fluid Mech. 2000, 422, 1–54. [Google Scholar] [CrossRef] [Green Version]
  12. Guala, M.; Hommema, S.E.; Adrian, R.J. Large-scale and very-large-scale motions in turbulent pipe flow. J. Fluid Mech. 2006, 554, 521–542. [Google Scholar] [CrossRef]
  13. Balakumar, B.J.; Adrian, R.J. Large-and very-large-scale motions in channel and boundary-layer flows. Philos. Trans. R. Soc. Lond. A 2007, 365, 665–681. [Google Scholar] [CrossRef] [PubMed]
  14. Stull, R.B. An Introduction to Boundary Layer Meteorology; Kluwer Academic: Dodrecht, The Netherlands, 1988. [Google Scholar]
  15. Arya, S.P. Introduction to Micrometeorology; Academic Press: New York, NY, USA, 2001. [Google Scholar]
  16. Kumar, V.; Kleissl, J.; Meneveau, C.; Parlange, M.B. Large-eddy simulation of a diurnal cycle of the atmospheric boundary layer: Atmospheric stability and scaling issues. Water Resour. Res. 2006, 42, 650–664. [Google Scholar] [CrossRef] [Green Version]
  17. Park, S.B.; Baik, J.J. A large-eddy simulation study of thermal effects on turbulence coherent structures in and above a building array. J. Appl. Meteorol. Climatol. 2013, 52, 1348–1365. [Google Scholar] [CrossRef]
  18. Li, D.; Bou-Zeid, E. Coherent structures and the dissimilarity of turbulent transport of momentum and scalars in the unstable atmospheric surface layer. Bound. Layer Meteorol. 2011, 140, 243–262. [Google Scholar] [CrossRef]
  19. Basu, S.; Porté-Agel, F. Large-eddy simulation of stably stratified atmospheric boundary layer turbulence: A scale-dependent dynamic modeling approach. J. Atmos. Sci. 2005, 63, 2074–2091. [Google Scholar] [CrossRef] [Green Version]
  20. Fang, J.N.; Porté-Agel, F. Large-eddy simulation of very-large-scale motions in the neutrally stratified atmospheric boundary layer. Bound. Layer Meteorol. 2015, 155, 397–416. [Google Scholar] [CrossRef] [Green Version]
  21. Mason, P.J.; Sykes, R.I. A two-dimensional numerical study of horizontal roll vortices in the neutral atmospheric boundary-layer. Q. J. R. Meteorol. Soc. 1980, 106, 351–366. [Google Scholar] [CrossRef]
  22. Etling, D.; Brown, R.A. Roll vortices in the planetary boundary-layer—A review. Bound. Layer Meteorol. 1993, 65, 215–248. [Google Scholar] [CrossRef]
  23. Khanna, S.; Brasseur, J.G. Three-dimensional buoyancy-and shear-induced local structure of the atmospheric boundary layer. J. Atmos. Sci. 1998, 55, 710–743. [Google Scholar] [CrossRef]
  24. Boppe, R.S.; Neu, W.L.; Shuai, H. Large-scale motions in the marine atmospheric surface layer. Bound. Layer Meteorol. 1999, 92, 165–183. [Google Scholar] [CrossRef]
  25. Hommema, S.E.; Adrian, R.J. Packet structure of surface eddies in the atmospheric boundary layer. Bound. Layer Meteorol. 2003, 106, 147–170. [Google Scholar] [CrossRef]
  26. Carper, M.A.; Porté-Agel, F. The role of coherent structures in subfilter-scale dissipation of turbulence measured in the atmospheric surface layer. J. Turbul. 2004, 5, 32. [Google Scholar] [CrossRef] [Green Version]
  27. Chauhan, K.; Hutchins, N.; Monty, J.; Marusic, I. Structure Inclination Angles in the Convective Atmospheric Surface Layer. Bound. Layer Meteorol. 2013, 147, 41–50. [Google Scholar] [CrossRef]
  28. Salesky, S.T.; Anderson, W. Buoyancy effects on large-scale motions in convective atmospheric boundary layers: Implications for modulation of near-wall processes. J. Fluid Mech. 2018, 856, 135–168. [Google Scholar] [CrossRef]
  29. Cortese, T.; Balachandar, S. Vortical nature of thermal plumes in turbulent convection. Phys. Fluids 1993, 5, 3226–3232. [Google Scholar] [CrossRef]
  30. Zhou, S.Q.; Xia, K.Q. Plume statistics in thermal turbulence: Mixing of an active scalar. Phys. Rev. Lett. 2002, 89, 184502. [Google Scholar] [CrossRef] [Green Version]
  31. Mahrt, L. Stratified atmospheric boundary layers and breakdown of models. Theor. Comput. Fluid Dyn. 1998, 11, 263–279. [Google Scholar] [CrossRef]
  32. Monin, A.S. Atmospheric boundary layer. Annu. Rev. Fluid Mech. 1970, 2, 225. [Google Scholar] [CrossRef]
  33. Derbyshire, S.H. Stable boundary-layer modelling: Established approaches and beyond. Bound. Layer Meteorol. 1999, 90, 423–446. [Google Scholar] [CrossRef]
  34. Holtslag, A.A.M. GABLS Initiates Intercomparison for Stable Boundary Layer Case; GEWEX News, No. 13.; International GEWEX Project Office: Silver Spring, MD, USA, 2003; pp. 7–8. [Google Scholar]
  35. Hunt, J.C.R.; Shutts, G.J.; Derbyshire, S.H. Stably stratified flows in meteorology. Dyn. Atmos. Oceans 1996, 23, 63–79. [Google Scholar] [CrossRef]
  36. Deardorff, J.W. Theoretical expression for the countergradient vertical heat flux. J. Geophys. Res. 1972, 77, 5900–5904. [Google Scholar] [CrossRef]
  37. Deardorff, J.W. Three-dimensional numerical study of turbulence in an entraining mixed layer. Bound. Layer Meteorol. 1974, 7, 199–226. [Google Scholar] [CrossRef]
  38. Moeng, C.H. A large-eddy simulation model for the study of planetary boundary-layer turbulence. J. Atmos. Sci. 1984, 41, 2052–2062. [Google Scholar] [CrossRef]
  39. Sullivan, P.P.; McWilliams, J.C.; Moeng, C.H. A subgrid-scale model for large-eddy simulation of planetary boundary-layer flows. Bound. Layer Meteorol. 1994, 71, 247–276. [Google Scholar] [CrossRef] [Green Version]
  40. Porté-Agel, F.; Meneveau, C.; Parlange, M.B. A scale-dependent dynamic model for large-eddy simulations: Application to a neutral atmospheric boundary layer. J. Fluid Mech. 2000, 415, 261–284. [Google Scholar] [CrossRef] [Green Version]
  41. Beare, R.J.; Macvean, M.K.; Holtslag, A.A.M.; Cuxart, J.; Esau, I.; Golaz, J.C.; Jimenez, M.A.; Khairoutdinov, M.; Kosovic, B.; Lewellen, D.; et al. An intercomparison of large-eddy simulations of the stable boundary layer. Bound. Layer Meteorol. 2006, 118, 247–272. [Google Scholar] [CrossRef] [Green Version]
  42. Schmidt, H.; Schumann, U. Coherent structure of the convective boundary layer derived from large-eddy simulations. J. Fluid Mech. 1989, 200, 511–562. [Google Scholar] [CrossRef] [Green Version]
  43. Nieuwstadt, F.T.M.; Mason, P.J.; Moeng, C.H.; Schumann, U. Large-eddy simulation of the convective boundary layer: A comparison of four computer codes. In Turbulent Shear Flows 8; Durst, F., Ed.; Springer: Berlin/Heidelberg, Germany, 1991; pp. 343–367. [Google Scholar]
  44. Churchfield, M.J.; Moriarty, P.J.; Vijayakumar, G.; Brasseur, J. Wind Energy-Related Atmospheric Boundary-Layer Large-Eddy Simulation Using OpenFOAM; NREL/CP-500-48905; National Renewable Energy Lab.: Golden, CO, USA, 2010.
  45. Andren, H. Effects of Habitat Fragmentation on Birds and Mammals in Landscapes with Different Proportions of Suitable Habitat: A Review. Oikos 1994, 71, 355–366. [Google Scholar] [CrossRef] [Green Version]
  46. Flores, O.; Riley, J.J. Analysis of turbulence collapse in the stably stratified surface layer using direct numerical simulation. Bound. Layer Meteorol. 2011, 139, 241–259. [Google Scholar] [CrossRef]
  47. Mason, P.J.; Derbyshire, S.H. Large-eddy simulation of the stably stratified atmospheric boundary layer. Bound. Layer Meteorol. 1990, 53, 117–162. [Google Scholar] [CrossRef]
  48. Brown, A.; Derbyshire, S.H.; Mason, P.J. Large-eddy simulation of stable atmospheric boundary layers with a revised stochastic subgrid model. Q. J. R. Meteorol. Soc. 1994, 120, 1485–1512. [Google Scholar] [CrossRef]
  49. Saiki, E.M.; Moeng, C.H.; Sullivan, P.P. Large-eddy simulation of the stably stratified planetary boundary layer. Bound. Layer Meteorol. 2000, 95, 1–30. [Google Scholar] [CrossRef] [Green Version]
  50. Churchfield, M.J.; Sang, L.; Moriarty, P.J. Adding Complex Terrain and Stable Atmospheric Condition Capability to the OpenFOAM-Based Flow Solver of the Simulator for on/Offshore wind Farm Applications (SOWFA); National Renewable Energy Lab.: Golden, CO, USA, 2013; Volume 2.
  51. Ren, H.; Laima, S.; Li, H. Numerical study of amplitude modulation in the atmospheric boundary layer at very high Reynolds number. AIP Adv. 2019, 9, 105112. [Google Scholar] [CrossRef]
  52. Schumann, U. Subgrid-scale model for finite-difference simulations of turbulent flow in plane channels and annuli. J. Comput. Phys. 1975, 18, 376–404. [Google Scholar] [CrossRef] [Green Version]
  53. Yoshizawa, A.; Horiuti, K. A statistically-derived subgrid-scale kinetic energy model for the large-eddy simulation of turbulent flows. J. Phys. Soc. Jpn. 1985, 54, 2834–2839. [Google Scholar] [CrossRef]
  54. Lee, J.; Lee, J.H.; Choi, J.; Sung, H.J. Spatial organization of large- and very-large-scale motions in a turbulent channel flow. J. Fluid Mech. 2014, 749, 818–840. [Google Scholar] [CrossRef]
  55. Ren, H.; Laima, S.; Li, H. Numerical Investigation of Very-Large-Scale Motions in a Turbulent Boundary Layer for Different Roughness. Energies 2020, 13, 659. [Google Scholar] [CrossRef] [Green Version]
  56. Wyngaard, J.C.; Cote, O.R. Cospectral similarity in the atmospheric surface layer. Q. J. R. Meteorol. Soc. 1972, 98, 590–603. [Google Scholar] [CrossRef]
  57. Zhang, H.S.; Chen, J.Y.; Zhang, A.C.; Wang, J.M.; Mitsuta, Y. An experiment and the results on flux-gradient relationships in the atmospheric surface over Gobi Desert Surface. In Proceedings of the International Symposium on HEIFE, Kyoto, Japan, 1–5 March 1993; Volume IV-12, pp. 349–362. [Google Scholar]
  58. Zhang, H.; Chen, J.Y.; Park, S.U. Turbulence structure in unstable conditions over various surfaces. Bound. Layer Meteorol. 2001, 100, 243–261. [Google Scholar] [CrossRef]
  59. Panofsky, H.A.; Tennekes, H.; Lenschow, D.H.; Wyngaard, J.C. The characteristics of turbulent velocity components in the surface layer under convective conditions. Bound. Layer Meteorol. 1977, 11, 355–361. [Google Scholar] [CrossRef]
  60. Bowne, N.E.; Ball, J.T. Observational comparison of rural and urban boundary-layer turbulence. J. Appl. Meteorol. 1970, 9, 862–873. [Google Scholar] [CrossRef] [Green Version]
  61. Counihan, J. Adiabatic atmospheric boundary layers: A review and analysis of data from the period 1880–1972. Atmos. Environ. 1975, 19, 871–905. [Google Scholar] [CrossRef]
  62. Steyn, D.G. Turbulence in an unstable surface layer over suburban terrain. Bound. Layer Meteorol. 1982, 22, 183–191. [Google Scholar] [CrossRef]
  63. Rotach, M.W. Turbulence Within and Above an Urban Canopy; ETH: Zurich, Switzerland, 1991. [Google Scholar]
  64. Roth, M. urbulent transfer relationships over an urban surface. II: Integral statistics. Q. J. R. Meteorol. Soc. 1993, 119, 1105–1120. [Google Scholar] [CrossRef]
  65. Lee, J.H.; Sung, H.J. Very-large-scale motions in a turbulent boundary layer. J. Fluid Mech. 2011, 673, 80–120. [Google Scholar] [CrossRef]
  66. Dennis, D.J.C.; Nickels, T.B. Experimental measurement of large-scale three-dimensional structures in a turbulent boundary layer. Part 1. Vortex packets. J. Fluid Mech. 2011, 673, 180–217. [Google Scholar] [CrossRef]
  67. Dennis, D.J.C.; Nickels, T.B. Experimental measurement of large-scale three-dimensional structures in a turbulent boundary layer. Part 2. Long structures. J. Fluid Mech. 2011, 673, 218–244. [Google Scholar] [CrossRef]
  68. Baltzer, J.R.; Adrian, R.J.; Wu, X. Structural organization of large and very large scales in turbulent pipe flow simulation. J. Fluid Mech. 2013, 720, 236–279. [Google Scholar] [CrossRef]
  69. Christensen, K.T.; Wu, Y. Characteristics of vortex organization in the outer layer of wall turbulence. In Proceedings of the Fourth International Symposium on Turbulence and Shear Flow Phenomena, Williamsburg, VA, USA, 27–29 June 2005; Volume 3, pp. 1025–1030. [Google Scholar]
  70. Volino, R.J.; Schultz, M.P.; Flack, K.A. Turbulence structure in rough- and smooth-wall boundary layers. J. Fluid Mech. 2007, 592, 263–293. [Google Scholar] [CrossRef]
  71. Hutchins, N.; Marusic, I. Evidence of very long meandering features in the logarithmic region of turbulent boundary layers. J. Fluid Mech. 2007, 579, 1–28. [Google Scholar] [CrossRef] [Green Version]
  72. Mathis, R.; Hutchins, N.; Marusic, I. Large-scale amplitude modulation of the small-scale structures in turbulent boundary layers. J. Fluid Mech. 2009, 628, 311–337. [Google Scholar] [CrossRef] [Green Version]
  73. Balasubramaniam, B.J. Nature of turbulence in wall-bounded flows. Ph.D. Thesis, University of Illinois at Urbana-Champaign, Urbana, IL, USA, 2005. [Google Scholar]
Figure 1. Ratio of horizontal and vertical heat fluxes u θ ¯ / w θ ¯ versus z / L . In (a), for unstable conditions, the blue circles, red crosses, green plus signs, cyan stars, magenta squares, yellow asterisks, black diamonds and blue rightward-pointing triangles represent cases 17, 16, 15, 14, 13, 12, 11 and 10, respectively, from very to slightly unstable. The curve corresponds to Equation (7). In (b), for stable conditions, the blue circles, red crosses, green plus signs, cyan stars, magenta squares, yellow asterisks, black diamonds and blue rightward-pointing triangles represents cases 8, 7, 6, 5, 4, 3, 2, and 1, respectively, going from slightly to very stable.
Figure 1. Ratio of horizontal and vertical heat fluxes u θ ¯ / w θ ¯ versus z / L . In (a), for unstable conditions, the blue circles, red crosses, green plus signs, cyan stars, magenta squares, yellow asterisks, black diamonds and blue rightward-pointing triangles represent cases 17, 16, 15, 14, 13, 12, 11 and 10, respectively, from very to slightly unstable. The curve corresponds to Equation (7). In (b), for stable conditions, the blue circles, red crosses, green plus signs, cyan stars, magenta squares, yellow asterisks, black diamonds and blue rightward-pointing triangles represents cases 8, 7, 6, 5, 4, 3, 2, and 1, respectively, going from slightly to very stable.
Energies 13 01068 g001
Figure 2. Normalized standard deviation of w ( σ w / u * ) versus z / L : (a) unstable conditions; (b) neutral stability conditions; and (c) stable conditions. The symbols are the same as in Figure 1, except that the curve in (a) now corresponds to Equation (8).
Figure 2. Normalized standard deviation of w ( σ w / u * ) versus z / L : (a) unstable conditions; (b) neutral stability conditions; and (c) stable conditions. The symbols are the same as in Figure 1, except that the curve in (a) now corresponds to Equation (8).
Energies 13 01068 g002
Figure 3. Streamwise velocity in the x y plane (whole domain) at a height of 0 . 2 δ for each case. (a) to (q) represent cases 1 to 17, respectively, going from very stable, to neutral and, finally, very unstable conditions.
Figure 3. Streamwise velocity in the x y plane (whole domain) at a height of 0 . 2 δ for each case. (a) to (q) represent cases 1 to 17, respectively, going from very stable, to neutral and, finally, very unstable conditions.
Energies 13 01068 g003aEnergies 13 01068 g003b
Figure 4. Contours of the two-point correlation of streamwise velocity fluctuations. Subfigures (a1) to (q2) means: each row represents case 1 to 17, from very stable, to neutral, and then very unstable conditions, with the left and right columns showing the contours in the x z and y z planes, respectively.
Figure 4. Contours of the two-point correlation of streamwise velocity fluctuations. Subfigures (a1) to (q2) means: each row represents case 1 to 17, from very stable, to neutral, and then very unstable conditions, with the left and right columns showing the contours in the x z and y z planes, respectively.
Energies 13 01068 g004aEnergies 13 01068 g004bEnergies 13 01068 g004cEnergies 13 01068 g004d
Figure 5. Extents of the 0.5 contours of the two-point correlations of the fluctuating velocity as functions of z / L at a height 0.2 δ : (a) streamwise extent; (b) spanwise extent; and (c) wall-normal extent. The left and right columns are for stable and unstable conditions, respectively. For the neutral stability condition, L x u u = 0 . 32 δ , L y u u = 0 . 20 δ and L z u u = 0 . 12 δ .
Figure 5. Extents of the 0.5 contours of the two-point correlations of the fluctuating velocity as functions of z / L at a height 0.2 δ : (a) streamwise extent; (b) spanwise extent; and (c) wall-normal extent. The left and right columns are for stable and unstable conditions, respectively. For the neutral stability condition, L x u u = 0 . 32 δ , L y u u = 0 . 20 δ and L z u u = 0 . 12 δ .
Energies 13 01068 g005aEnergies 13 01068 g005b
Figure 6. Pre-multiplied wind velocity spectra of the streamwise velocity fluctuation for each case at different heights: (a) to (q) are for cases 1 to 17, respectively. The dashed and solid black lines represent 0 and –2/3 slope, respectively.
Figure 6. Pre-multiplied wind velocity spectra of the streamwise velocity fluctuation for each case at different heights: (a) to (q) are for cases 1 to 17, respectively. The dashed and solid black lines represent 0 and –2/3 slope, respectively.
Energies 13 01068 g006aEnergies 13 01068 g006b

Share and Cite

MDPI and ACS Style

Laima, S.; Ren, H.; Li, H.; Ou, J. Numerical Simulation of Coherent Structures in the Turbulent Boundary Layer under Different Stability Conditions. Energies 2020, 13, 1068. https://doi.org/10.3390/en13051068

AMA Style

Laima S, Ren H, Li H, Ou J. Numerical Simulation of Coherent Structures in the Turbulent Boundary Layer under Different Stability Conditions. Energies. 2020; 13(5):1068. https://doi.org/10.3390/en13051068

Chicago/Turabian Style

Laima, Shujin, Hehe Ren, Hui Li, and Jinping Ou. 2020. "Numerical Simulation of Coherent Structures in the Turbulent Boundary Layer under Different Stability Conditions" Energies 13, no. 5: 1068. https://doi.org/10.3390/en13051068

APA Style

Laima, S., Ren, H., Li, H., & Ou, J. (2020). Numerical Simulation of Coherent Structures in the Turbulent Boundary Layer under Different Stability Conditions. Energies, 13(5), 1068. https://doi.org/10.3390/en13051068

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop