Search NASA⌕ Search

SEARCH · Search NASA

Results for “wake effects”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Investigating the interactions between wakes and floating wind turbines using FAST.Farm

As floating offshore wind progresses to commercial maturity, wake and array effects across a farm of floating offshore wind turbines (FOWTs) will become increasingly important. While wakes of land-based and bottom-fixed offshore wind turbines have been extensively studied, only recently has this topic become relevant for floating turbines. This work presents an investigation of the mutual interaction between the motions of floating wind turbines and wakes using FAST.Farm. While FAST.Farm has been extensively validated across a wide range of conditions, it has never been validated for FOWT applications. Hence, in the first part of this work, we validate FAST.Farm by comparing simulations of a single FOWT against high-fidelity results from large-eddy simulations available in the literature. The validation is based on wake meandering, mean wake deflection, and velocity deficit at different downstream locations. This validation showed that the original axisymmetric (polar) wake model of FAST.Farm overpredicts the vertical wake deflection induced by shaft tilt and floater pitch, while the new curled wake model is capable of properly capturing the vertical wake deflection. In the second part, we use FAST.Farm to analyze a small three-unit array of FOWTs with a spacing of 7 diameters across a wide range of environmental conditions. The same National Renewable Energy Laboratory 5 MW reference wind turbine atop the OC4-DeepCwind semisubmersible is adopted for the three FOWTs and for the validation against high-fidelity simulations. To assess the effect of the floating substructure, we compare the power production, tower-base moments, and blade-root moments obtained for the floating turbines with the results obtained in a fixed-bottom configuration. The main differences introduced by the floating substructure are the motions induced by the waves, the change in the natural frequencies of the tower caused by differences in the boundary condition at its base, and the larger vertical deflection of the wake deficit due to the mean pitch of the platform. The impact of these differences, as well as other minor effects, are analyzed in detail.

16 TIDAL AND WAVE POWER↗
DOE OSTI2021

A vortex sheet based analytical model of the curled wake behind yawed wind turbines

Motivated by the need for compact descriptions of the evolution of non-classical wakes behind yawed wind turbines, we develop an analytical model to predict the shape of curled wakes. Interest in such modelling arises due to the potential of wake steering as a strategy for mitigating power reduction and unsteady loading of downstream turbines in wind farms. We first estimate the distribution of the shed vorticity at the wake edge due to both yaw offset and rotating blades. By considering the wake edge as an ideally thin vortex sheet, we describe its evolution in time moving with the flow. Vortex sheet equations are solved using a power series expansion method, and an approximate solution for the wake shape is obtained. The vortex sheet time evolution is then mapped into a spatial evolution by using a convection velocity. Apart from the wake shape, the lateral deflection of the wake including ground effects is modelled. Our results show that there exists a universal solution for the shape of curled wakes if suitable dimensionless variables are employed. For the case of turbulent boundary layer inflow, the decay of vortex sheet circulation due to turbulent diffusion is included. Finally, we modify the Gaussian wake model by incorporating the predicted shape and deflection of the curled wake, so that we can calculate the wake profiles behind yawed turbines. Model predictions are validated against large-eddy simulations and laboratory experiments for turbines with various operating conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Analyzing the Effects of Atmospheric Turbulent Fluctuations on the Wake Structure of Wind Turbines and Their Blade Vibrational Dynamics

In recent trends, a rising demand for renewable energy has driven wind turbines to larger proportions, where lighter blade designs are often adopted to reduce the costs associated with logistics and production. This causes modern utility-scale wind turbine blades to be inherently more flexible, and their amplified aeroelastic sensitivity results in complex multi-physics reactions to variant atmospheric conditions, including dynamic patterns of aerodynamic loading at the rotor and vortex structure evolutions within the wake. In this paper, we analyze the influence of inflow variance for wind turbines with large, flexible rotors through simulations of the National Rotor Testbed (NRT) turbine, located at Sandia National Labs’ Scaled Wind Farm Technology (SWiFT) facility in Lubbock, Texas. The Common Ordinary Differential Equation Framework (CODEF) modeling suite is used to simulate wind turbine aeroelastic oscillatory behavior and wind farm vortex wake interactions for a range of flexible NRT blade variations, operating in differing conditions of variant atmospheric flow. CODEF solutions of turbine operation in Steady-In-The-Average (SITA) wind conditions are compared to SITA wind conditions featuring a controlled gust-like pulse overimposed, to isolate the effects of typical wind fluctuations. Finally, simulations of realistic time-varying wind conditions from SWiFT meteorological tower measurements are compared to the solutions of SITA wind conditions. These increasingly complex atmospheric inflow variations are tested to show the differing effects evoked by various patterns of spatiotemporal atmospheric flow fluctuations. An analysis is presented for solutions of wind turbine aeroelastic response and vortex wake evolution, to elucidate the consequences of variant inflow, which pertain to wind turbine dynamics at an individual and farm-collective scale. The comparisons of simulated farm flow for SITA and measured fluctuating wind conditions show that certain regions of the wake contain up to a 12% difference in normalized axial velocity, due to the introduction of wind fluctuations. The findings of this study prove valuable for practical applications in wind farm control and optimization strategies, with particular significance for modern utility-scale wind power plants operating in variant atmospheric conditions.

Farrell, Alayna (ORCID:000000023555720X)↗
DOE OSTI2022

High-Fidelity Simulation Study of the Unsteady Flow Effects on High-Pressure Turbine Blade Performance

Abstract Unsteadiness, in the form of both broadband background disturbances and discrete coherent wakes, can have a strong effect on the performance of turbomachinery blades. The influence of the incoming flow has received much interest as it inevitably affects the blade boundary layers and develops as it passes through the machine. In the present work, we investigate the effect of unsteady flow on high-pressure turbines (HPTs), using high-fidelity datasets produced by wall resolved large-eddy simulation of an HPT stage. The effects of incident wakes from an upstream stator, compounded by the presence of freestream turbulence, on the downstream rotor are investigated. Based on analyzing cases with different turbulence intensities and length scales prescribed at the inlet, we show that changing the freestream turbulence characteristics has a direct effect on the unsteady behavior of the stator wakes. As a result, the performance of the rotor is also significantly affected. By detailing the influence of the wake–turbulence interaction, we aim to distinguish driving forces on rotor performance, be it changes in the incident wakes or direct influence from the freestream turbulence. Furthermore, the aerothermal behaviors of the rotor blades have been extensively investigated, showing that the blade boundary layers on the suction and pressure sides respond differently to external disturbances. The insights gained can provide designers with guidelines in understanding the unsteady flow effects of a given flow state, and how the unsteadiness present, either broadband or deterministic, will affect the performance of downstream blades.

Engineering↗
DOE OSTI2020

Assessing the blockage effect of wind turbines and wind farms using an analytical vortex model

Wind farm blockage effects are currently neglected in the prediction of wind farm energy yield, typically leading to an overestimation of the production. This work presents a novel method to assess wind farm production, while accounting for blockage effects. We apply a vortex model, based on a cylindrical wake, to assess induction effects. We present variations of the model to account for finite wake length, finite tip-speed ratios, and the proximity to the ground. The results are applied to single rotors in aligned and yawed conditions and to different wind farm layouts. We provide far-field approximations for faster estimates of the velocity field. Further, this article includes a new methodology to couple the induction model to engineering wake models, such as the ones present in the FLOw Redirection and Induction in Steady State (FLORIS). We compare the results to actuator disk simulations for various operating conditions of a single turbine and different wind farm layouts. We found that the mean relative error of the model in the induction zone is typically around 0.2% compared with actuator disk simulations. The computational time of the velocity field using the analytical vortex model is three orders of magnitude less than the one obtained with the actuator disk simulation.

17 WIND ENERGY↗
DOE OSTI2023

A Study of Wake Characteristics of Marine Turbine Arrays: Preprint

As the marine energy industry continues to develop and installations grow from single turbines to small arrays and larger farms, understanding and predicting wake behavior using modeling tools is necessary for array design and optimization. The wake characteristics of upstream turbines, including velocity deficit, wake swirling, tip vortices, wake shape and direction, and flow recovery, influence the inflow conditions experienced by downstream turbines in an array and, subsequently, affect the individual turbine loading and power generated by the array. These wake characteristics are influenced by many factors, including site geography, inflow shear profile, turbulence, waves, the strength of the current resource, gravitational forces, storm surges, and seasonal effects like temperature and wind conditions. Efficient array modeling tools that capture these wake effects and resulting turbine-to-turbine interactions need to be developed to optimize farm layouts and design individual turbines for operation within arrays. Unlike the wind industry, the marine energy research community is largely limited to the use of either numerical tools or tank/flume testing to understand marine turbine wake interactions. This paper will investigate recent advances in numerical and experimental marine turbine wake modeling, with a focus on recent publications and on both axial-flow and cross-flow marine turbines. A summary of recent research on these topics is provided. Additionally, a demonstration of turbine wake interactions using computational fluid dynamics is presented.

marine hydrokinetic turbine↗
DOE OSTI2020

How Does Inow Veer Affect the Veer of a Wind-Turbine Wake?

Stably stratified flow conditions often exhibit wind veer, or a change of wind direction with height. When wind turbines experience this veered flow, the resulting wake structure tends to exhibit a stretching into an ellipsoid, rather than a symmetric shape or a curled shape. Observational studies suggest that the magnitude of wake veer is less than the veer of the inflow, whereas large-eddy simulations with actuator disk models and actuator line models suggest a range of relationships between inflow veer and wake veer. Here we present a series of large-eddy simulations with a range of veer shapes, a range of magnitudes of veer, a range of wind speeds, and both rotational directions of the wind-turbine rotor investigating the effect on the wake deflection angle. These results can guide the application of wake steering in stably stratified flow.

17 WIND ENERGY↗
DOE OSTI2025

Modeling and Analysis of Wind Turbine Wake Vortex Evolution Due to Time-Constant Spatial Variations in Atmospheric Flow

Modern utility-scale wind turbines are evolving toward larger, lighter, and more flexible designs to meet the growing demand for renewable energy while minimizing logistical costs. However, these advancements in lightweight design result in heightened aeroelastic sensitivity, leading to complex interactions which affect the rotor’s capacity to withstand aerodynamic loading and the cascading effects that manifest in the wake’s vortex-structure evolution under variable atmospheric conditions. In this paper, we analyze the influence of stream-wise fluctuating atmospheric flow conditions on wind turbines with large, flexible rotors through simulations of the National Rotor Testbed (NRT) turbine, located at Sandia National Labs’ Scaled Wind Farm Technology (SWiFT) facility in Lubbock, Texas. The Common Ordinary Differential Equation Framework (CODEF) modeling suite is used to simulate wind turbine aeroelastic oscillatory behavior and wind farm vortex–wake interactions for a range of conditions with spatially variant atmospheric flow. CODEF solutions for turbine operation in wind conditions featuring only one parameter fluctuation are compared to wind conditions with several wind parameter variations in combination. By isolating individual inflow variations and comparing them to multi-parameter scenarios, we determine the contributions of each atmospheric factor to rotor dynamics, wake evolution, and downstream wind farm interactions. The purpose of this paper is to analyze the effects of spatial variations in atmospheric flow on the topological evolution of wind turbine vortex wakes, which constitutes a gap in the current understanding of wind turbine wake dynamics. The insights gained from this study are particularly valuable for the development of wind farm control strategies aimed at mitigating the adverse effects of wake interactions, enhancing energy capture, and improving the overall stability of wind farm operations. With these insights, we aim to contribute to the development of modeling and simulation tools to optimize utility-scale wind power plants operating in diverse atmospheric environments.

Farrell, Alayna (ORCID:000000023555720X)↗
DOE OSTI2026

Idealized simulations of wind farm interactions with intermittent turbulence in stable boundary layer conditions

Stable atmospheric boundary layer conditions typically correspond to weak turbulence levels, but intermittent periods of elevated turbulence can occur during otherwise quiescent conditions. The interaction between intermittent turbulence and wind turbines is not well understood because of sparse observations, as well as the difficulty in realistically resolving small-scale turbulence during strongly stable conditions with numerical simulations. In this study, an explicit filtering and reconstruction approach for large-eddy simulation (LES) is used to simulate weakly and strongly stable conditions, with surface cooling rates of −0.2 and −2.0 K h −1 , respectively. This approach can sustain resolved background turbulence at relatively coarse grid spacing and stronger stratification compared to conventional closures, permitting more realistic intermittent stable boundary layer (SBL) turbulence. The idealized LES capability of the Weather Research and Forecasting model is employed with turbine rotors parameterized using generalized actuator disks to examine (1) how the presence of turbine wakes affects SBL evolution and (2) the effect of intermittent turbulence on power production and wake recovery. Wakes increase mixing and deepen the SBL, with a stronger effect under strongly stable conditions, primarily because the SBL is shallower and closer to the top of the wind turbine rotor layer. Intermittent turbulence does not have a significant impact on mean power generation and wake recovery because the relevant intermittent turbulent structures in this study only affect the bottom half of the rotor disk. Power production is, however, more variable during periods of elevated turbulence, demonstrating the impact of SBL intermittency. This study uses an idealized configuration, focusing on LES model performance and physical understanding, with the goal of informing future simulations of the conditions observed during the American Wake Experiment.

Energy - Wind↗
DOE OSTI2023

Increased power gains from wake steering control using preview wind direction information

Abstract. Yaw controllers typically rely on measurements taken at the wind turbine, resulting in a slow reaction to wind direction changes and subsequent power losses due to misalignments. Delayed yaw action is especially problematic in wake steering operation because it can result in power losses when the yaw misalignment angle deviates from the intended one due to a changing wind direction. This study explores the use of preview wind direction information for wake steering control in a two-turbine setup with a wind speed in the partial load range. For these conditions and a simple yaw controller, results from an engineering model identify an optimum preview time of 90 s. These results are validated by forcing wind direction changes in a large-eddy simulation model. For a set of six simulations with large wind direction changes, the average power gain from wake steering increases from only 0.44 % to 1.32 %. For a second set of six simulations with smaller wind direction changes, the average power gain from wake steering increases from 1.24 % to 1.85 %. Low-frequency fluctuations are shown to have a larger impact on the performance of wake steering and the effectiveness of preview control, in particular, than high-frequency fluctuations. From these results, it is concluded that the benefit of preview wind direction control for wake steering is substantial, making it a topic worth pursuing in future work.

17 WIND ENERGY↗
DOE OSTI2024

Near wake evolution of a tidal stream turbine due to asymmetric sheared turbulent inflow with different integral length scales

Tidal stream turbines deployed at highly energetic open water sites are subjected to sheared inflow in the rotor plane. The inflow shear is expected to cause asymmetric loading on the rotor blades and affect the downstream wake. In the current study, two different turbulent inflow conditions, static-high shear and dynamic shear, were generated via an active-grid turbulence generator. A 1:20 scaled three-bladed horizontal axis tidal turbine model was tested in those conditions. The results were compared to a quasi-laminar case with no imposed turbulence or shear. The results show that the high shear reduces the average performance, with a drop of up to 16% in the optimal power coefficient. Besides, the shear profiles increase torque fluctuations and induce significant differences in wake hydrodynamics between the high-speed (upper) and low-speed (lower) regions. The large integral length scales further enhance the load fluctuations perceived by the rotor but have a negligible effect on the mean wake field quantities and the wake recovery. The lower half region featured a faster breakdown of tip vortex structure and a rapid drop of swirl number, a phenomenon conjectured to be a consequence of the strong turbulence intensities and Reynolds stresses in the lower half region. Furthermore, the sheared turbulent inflow also results in a very intensive energy redistribution process towards large-scale, low-frequency motions, which is important to the downstream turbines.

16 TIDAL AND WAVE POWER↗
DOE OSTI2020

Effect of optical field ionisation on the generation of wake fields by femtosecond laser pulses in an inhomogeneous plasma

We consider the interaction of a high-intensity short laser pulse with an argon gas during optical field ionisation. Modelling in a three-dimensional cylindrical symmetric geometry is performed at various moderately relativistic intensities and positions of the focal plane of the laser beam in the case of both optical gas ionisation and a pre-ionised plasma. The influence of the processes occurring during optical field ionisation of the gas on the generation of wake waves is investigated, and the conditions are found under which intense wake fields are produced in a plasma formed from an inhomogeneous argon gas jet. The possibility of using a gas with a large number of electrons on the outer shell (argon) to excite intense wake waves and accelerate electrons is demonstrated. Despite significant ionisation refraction of the laser pulse on a radially inhomogeneous density profile of plasma electrons formed during optical field ionisation, the region of the parameters of the laser pulse and the gas target is determined at which ionisation refraction leads to the excitation of an intense wake wave. It is found that in a pre-ionised plasma a wake wave is not generated even at the same laser radiation intensity. (paper)

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Impacts to FeRAM design arising from interfacial dielectric layers and wake up modulation in ferroelectric hafnium zirconium oxide

As ferroelectric hafnium zirconium oxide (HZO) becomes more widely utilized in ferroelectric microelectronics, integration impacts of intentional and non-intentional dielectric interfaces and their effects upon the ferroelectric film wake up and circuit parameters become important to understand. In this work, the effect of the addition of a linear dielectric aluminum oxide, Al 2 O 3 , below a ferroelectric Hf 0.58 Zr 0.42 O 2 film in a capacitor structure for FeRAM applications with NbN electrodes was measured. Depolarization fields resulting from the linear dielectric is observed to induce a reduction of the remanent polarization of the ferroelectric. Addition of the aluminum oxide also impacts the wake up of the HZO with respect to the cycling voltage applied. Intricately linked to the design of a FeRAM 1C/1T cell, the metal-ferroelectric-insulator-metal (MFIM) devices are observed to significantly shift charge related to the read states based on aluminum oxide thickness and wake up cycling voltage. As a result, a 33% reduction in the separation of read states is measured, which complicates how a memory cell is designed and illustrates the importance of clean interfaces in devices.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2024

Recommendations on setup in simulating atmospheric gravity waves under conventionally neutral boundary layer conditions

Wind farm-induced atmospheric gravity waves have been the subject of recent research as they can impact wind farm performance. Pressure variations associated with gravity waves can contribute to the global blockage effect and wind farm wake recovery. Therefore, accurate numerical simulation of flow fields, where wind-farm-induced gravity waves may be produced, is important. Three main considerations in such simulations are the overall domain size, the use of Rayleigh damping near domain boundaries to dampen gravity waves, and advection damping at the inlet to prevent spurious oscillations. Often these considerations are treated ad hoc rather than systematically. This work aims to test and extend the systematic modelling of internal gravity waves proposed in a preliminary investigation to modelling of both internal and trapped gravity waves. The preliminary study identifies the length scales to set the domain and damping layer sizes and the time scale to configure the Rayleigh damping coefficient but under linearly stratified conditions. Large eddy simulations of flow through a wind farm canopy are performed under conventionally neutral boundary layer (CNBL) conditions to test the validity of proposed setups for CNBL conditions. Background atmospheric parameters, such as Froude number (Fr), inversion height (H i ), and inversion layer Froude number (Fr i ) control most of the atmospheric gravity wave characteristics. We validated for CBNL conditions that the effective wavelengths of the internal gravity waves are the correct length scale to configure the domain size and damping layer thickness. Likewise, the optimum damping coefficient to dampen the internal gravity waves relates to the free atmosphere's buoyancy frequency or buoyant perturbations' time scale. We infer that the damping coefficient in the inversion layer may relate to the inversion buoyancy frequency to effectively dampen the trapped gravity waves. Moreover, the advection damping length is linked to the horizontal wavelength of the trapped gravity waves in the inversion layer to prevent spurious waves at the inlet by retaining wave energy accumulation.

17 WIND ENERGY↗
DOE OSTI2023

Influence of the Ozone Dose Time during Atomic Layer Deposition on the Ferroelectric and Pyroelectric Properties of 45 nm-Thick ZrO 2 Films

Over a decade ago, ferroelectricity was discovered in doped HfO 2 thin films. The HfO 2 -based thin films have attracted much attention due to their remarkable scalability and CMOS compatibility. Other than the HfO 2 -based thin films, the undoped ZrO 2 thin films are understudied despite their commonly reported antiferroelectric behavior. However, being of the same fluorite structure as HfO 2 -based thin films, the undoped ZrO 2 also displayed considerable ferroelectricity as demonstrated in recent studies. Here, 45 nm-thick polycrystalline undoped ZrO 2 films are synthesized using atomic layer deposition with different ozone dose times. The ZrO 2 films are crystallized after atomic layer deposition at 350 °C without anneals. In general, the longer ozone dose time causes a lower in-plane tensile stress and oxygen vacancy content, which help facilitate an irreversible non-polar tetragonal to polar orthorhombic phase transition with electric-field cycling. However, the lower in-plane tensile stress and oxygen vacancy content also stabilize the monoclinic phase so that a long ozone dose time (>17.5 s) reduces the ferroelectric behavior. After wake-up cycles, the ZrO 2 thin film with an ozone dose time of 17.5 s exhibits a remanent polarization of 6 μC·cm –2 and a pyroelectric coefficient of -35 μC·K –1 ·m –2 . Moreover, the wake-up behavior is consistent between the ferroelectric and pyroelectric response. As essential factors in optimizing the growth of fluorite-structure thin films for ferroelectric applications, the in-plane tensile stress and oxygen vacancy content significantly influence the ferroelectric and pyroelectric properties. Additionally, the low thermal budget for processing ferroelectric ZrO 2 thin films is valuable for semiconductor back-end-of-line processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Thermal evolution of ferroelectric behavior in epitaxial Hf 0.5 Zr 0.5 O 2

In this paper, we report a cryogenic-temperature study on the evolution of the ferroelectric properties of epitaxial Hf 0.5 Zr 0.5 O 2 thin films on silicon. Wake-up, endurance, and fatigue of these films are found to be intricately correlated, strongly hysteretic, and dependent on available thermal energy. Field-dependent measurements reveal a decrease in polarization with temperature, which has been determined not to be an intrinsic change of the material property, rather a demonstration of the increase in the coercive bias of the material. Our findings suggest that a deficiency in thermal energy suppresses the mobility of defects presumed to be oxygen vacancies during wake-up and trapped injected charge during fatigue, which is responsible for polarization evolution during cycling. This permits accelerated wake-up and fatigue effects at high temperatures where thermal energy is abundant but delays these effects at cryogenic temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Assessment of the MACCS Code Applicability for Nearfield Consequence Analysis

Three codes are used for comparisons in this report to evaluate the adequacy of MACCS in the nearfield, AERMOD, ARCON96 and QUIC. Test cases were developed to give a broad range of weather conditions, building dimensions and plume buoyancy. Based on the comparisons of MACCS with AERMOD, ARCON96 and QUIC across the test cases, the following observations are made: MACCS calculations configured with point-source, ground-level, nonbuoyant plumes provide nearfield results that bound the centerline, ground-level air concentrations from AERMOD, ARCON96, and QUIC. MACCS calculations with ground-level, nonbuoyant plumes that include the effects of the building wake (area source) provide nearfield results that bound the results from AERMOD and QUIC and the results from ARCON96 at distances >250 m. If using a point-source is too conservative and it is desired to bound the results from all three codes, another alternative is to use area source parameters in MACCS that are less than the standard values, i.e., an area source intermediate between the standard recommendation and a point source. All these options provide results from MACCS that are bounding for the test cases evaluated. Based on these observations, it appears that MACCS is adequate for use in nearfield calculations, given the correct parameterization.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Investigation of Farm-to-Farm Interactions and Blockage Effects from AWAKEN using Large-Scale Numerical Simulations

A large-scale numerical computation of five wind farms was performed as a part of the American Wake Experiment (AWAKEN). This high-fidelity computation used the ExaWind/AMR-Wind LES solver to simulate a 100km × 100km domain containing 541 turbines under unstable atmospheric conditions matching previous measurements. The turbines were represented by Joukowski and OpenFAST coupled actuator disk models. Results of this qualitative comparison illustrate the interactions with wind farms with large scale ABL structures in the flow, as well as the extent of downstream wake penetration in the flow and blockage effects around wind farms.

Cheung, Lawrence C.↗