The significance of air flow within dense canopies situated on hilly terrain is not in dispute given its relevance to a plethora of applications in meteorology, wind energy, air pollution, atmospheric chemistry and ecology. While the mathematical description of such flows is complex, progress has proceeded through an interplay between experiments, mathematical modelling, and more recently large‐eddy simulations (LESs). In this contribution, LES is used to investigate the topography‐induced changes in the flow field and how these changes propagate to scalar transport within the canopy. The LES runs are conducted for a neutral atmospheric boundary layer above a tall dense forested canopy situated on a train of two‐dimensional sinusoidal hills. The foliage distribution is specified using leaf area density measurements collected in an Amazon rain forest. A series of LES runs with increasing hill amplitude are conducted to disturb the flow from its flat‐terrain state. The LES runs successfully reproduce the recirculation region and the flow separation on the lee‐side of the hill within the canopy region in agreement with prior laboratory and LES studies. Simulation results show that air parcels released inside the canopy have two preferential pathways to escape the canopy region: a “local” pathway similar to that encountered in flat terrain and an “advective” pathway near the flow‐separation region. Further analysis shows that the preferential escape location over the flow‐separation region leads to a “chimney”‐like effect that becomes amplified for air parcel releases near the forest floor. The work here demonstrates that shear‐layer turbulence is the main mechanism exporting air parcels out the canopy for both pathways. However, compared to flat terrain, the mean updraught at the flow separation induced by topography significantly shortens the in‐canopy residence time for air parcels released in the lower canopy, thus enhancing the export fraction of reactive gases.
more »
« less
Wind-Tunnel Experiments of Turbulent Wind Fields over a Two-dimensional (2D) Steep Hill: Effects of the Stable Boundary Layer
Flow separation caused by steep topography remains a significant obstacle in accurately predicting turbulent boundary-layer flows over complex terrain, despite the utilization of sophisticated numerical models. The addition of atmospheric thermal stability, in conjunction with steep topography, further complicates the determination of disrupted turbulent wind patterns. The turbulent separated flows over a two-dimensional (2D) steep hill under thermal stratification has not been extensively addressed in previous experimental studies. Such measurements are crucial for enhancing our comprehension of flow physics and validating numerical models. We measured the turbulent wind flows over a 2D steep hill immersed in a stable boundary layer (of the bulk Richardson Number = 0.256) in a thermally-stratified boundary-layer wind tunnel. The flow separation, re-circulation zone and flow reattachment were characterized by the planar particle image velocimetry technique. Vertical profiles of mean air temperature and its fluctuations are also quantified at representative locations above the 2D steep hill and in the near wake region. Results indicate that the separated shear layer, initiated near the crest of the 2D steep hill, dominates the physical process leading to high turbulence levels and the turbulent kinetic energy production in the wake region for both stable and neutral thermal stability. Although the stable boundary layer does not dramatically change the turbulent flow pattern around the hill, the mean separation bubble is elongated by 13%, and its vertical extent is decreased by approximately 20%. Furthermore, the reduced turbulence intensities and turbulent kinetic energy of the near wake flow are attributed to the relatively low turbulence intensity and low momentum of the stable boundary layer due to buoyancy damping, compared to the neutral boundary layer. Additionally, a distinct low-temperature region—a cold pool—is extended beyond the separation bubble, reflecting the significant sheltering effect of the 2D steep hill on the downwind flow and temperature field.
more »
« less
- Award ID(s):
- 1944776
- PAR ID:
- 10492597
- Publisher / Repository:
- Springer
- Date Published:
- Journal Name:
- Boundary-Layer Meteorology
- Volume:
- 188
- Issue:
- 3
- ISSN:
- 0006-8314
- Page Range / eLocation ID:
- 441 to 461
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
As turbines continue to grow in hub height and rotor diameter and wind farms grow larger, consideration of stratified atmospheric boundary layer (ABL) processes in wind power models becomes increasingly important. Atmospheric stratification can considerably alter the boundary layer structure and flow characteristics through buoyant forcing. Variations in buoyancy, and corresponding ABL stability, in both space and time impact ABL wind speed shear, wind direction shear, boundary layer height, turbulence kinetic energy, and turbulence intensity. In addition, the presence of stratification will result in a direct buoyant forcing within the wake region. These ABL mechanisms affect turbine power production, the momentum and kinetic energy deficit wakes generated by turbines, and the turbulent mixing and kinetic energy entrainment in wind farms. Presently, state-of-practice engineering models of mean wake momentum utilize highly empirical turbulence models that do not explicitly account for ABL stability. Models also often neglect the interaction between the wake momentum deficit and the turbulence kinetic energy added by the wake, which depends on stratification. In this work, we develop a turbulence model that models the wake-added turbulence kinetic energy, and we couple it with a wake model based on the parabolized Reynolds-averaged Navier–Stokes equations. Comparing the model predictions to large eddy simulations across stabilities (Obukhov lengths) and surface roughness lengths, we find lower prediction error in both power production and the wake velocity field across the ABL conditions and error metrics investigated.more » « less
-
To achieve decarbonization targets, wind turbines are growing in hub height and rotor diameter, and they are being deployed in new locations with diverse atmospheric conditions not previously seen, such as offshore. Physics-based analytical wake models commonly used for design and control of wind farms simplify atmospheric boundary layer (ABL) and wake physics to achieve computational efficiency. This is accomplished primarily through a simplified model form that neglects certain flow processes, such as atmospheric stability, and through the parametrization of ABL and wake turbulence through a wake spreading rate. In this study, we systematically analyze the physical mechanisms that govern momentum and turbulence within a wind turbine wake in the stratified ABL. We use large-eddy simulation and analysis of the streamwise momentum deficit and wake-added turbulence kinetic energy (TKE) budgets to study wind turbine wakes under neutral and stable conditions. To parse the turbulence in the wake from the turbulent, incident ABL flow, we decompose the flow into the base ABL flow and the deficit flow produced by the presence of a turbine. We then analyze the decomposed flow field budgets to study the effects of changing stability on the streamwise momentum deficit and wake-added TKE. The results demonstrate that stability changes the relative balance of turbulence and advection for both the streamwise momentum deficit and wake-added TKE primarily through the nonlinear interactions of the base flow with the deficit flow. The stable cases are most affected by increased shear and veer in the base flow and the neutral case is most affected by the increased ambient turbulence intensity. These differences in the base flow that arise from stratification are relatively more important than the buoyancy forcing terms in the wake-added TKE budget. The wake-added TKE depends on the ABL stability. An existing wake-added TKE model that neglects the effects of ABL stability yields error compared to large-eddy simulation, with errors that are higher in stable conditions than neutral. These results motivate future research to develop fast-running models of wake-added TKE that account for stability effects. Published by the American Physical Society2024more » « less
-
Turbulent wake flows behind helical- and straight-bladed vertical axis wind turbines (VAWTs) in boundary layer turbulence are numerically studied using the large-eddy simulation (LES) method combined with the actuator line model. Based on the LES data, systematic statistical analyses are performed to explore the effects of blade geometry on the characteristics of the turbine wake. The time-averaged velocity fields show that the helical-bladed VAWT generates a mean vertical velocity along the center of the turbine wake, which causes a vertical inclination of the turbine wake and alters the vertical gradient of the mean streamwise velocity. Consequently, the intensities of the turbulent fluctuations and Reynolds shear stresses are also affected by the helical-shaped blades when compared with those in the straight-bladed VAWT case. The LES results also show that reversing the twist direction of the helical-bladed VAWT causes the spatial patterns of the turbulent wake flow statistics to be reversed in the vertical direction. Moreover, the mass and kinetic energy transports in the turbine wakes are directly visualized using the transport tube method, and the comparison between the helical- and straight-bladed VAWT cases show significant differences in the downstream evolution of the transport tubes.more » « less
-
An incoming canonical spatially developing turbulent boundary layer (SDTBL) over a 2-D curved hill is numerically investigated via the Reynolds-averaged Navier–Stokes (RANS) equations plus two eddy-viscosity models: the K−ω SST (henceforth SST) and the Spalart–Allmaras (henceforth SA) turbulence models. A spatially evolving thermal boundary layer has also been included, assuming temperature as a passive scalar (Pr = 0.71) and a turbulent Prandtl number, Prt, of 0.90 for wall-normal turbulent heat flux modeling. The complex flow with a combined strong adverse/favorable streamline curvature-driven pressure gradient caused by concave/convex surface curvatures has been replicated from wind-tunnel experiments from the literature, and the measured velocity and pressure fields have been used for validation purposes (the thermal field was not experimentally measured). Furthermore, direct numerical simulation (DNS) databases from the literature were also employed for the incoming turbulent flow assessment. Concerning first-order statistics, the SA model demonstrated a better agreement with experiments where the turbulent boundary layer remained attached, for instance, in Cp, Cf, and Us predictions. Conversely, the SST model has shown a slightly better match with experiments over the flow separation zone (in terms of Cp and Cf) and in Us profiles just upstream of the bubble. The Reynolds analogy, based on the St/(Cf/2) ratio, holds in zero-pressure gradient (ZPG) zones; however, it is significantly deteriorated by the presence of streamline curvature-driven pressure gradient, particularly due to concave wall curvature or adverse-pressure gradient (APG). In terms of second-order statistics, the SST model has better captured the positively correlated characteristics of u′ and v′ or positive Reynolds shear stresses ( > 0) inside the recirculating zone. Very strong APG induced outer secondary peaks in and turbulence production as well as an evident negative slope on the constant shear layer.more » « less
An official website of the United States government

