Abstract. Blade element momentum (BEM) codes traditionally neglect the radial induced velocity ur, although for the strongly deflected or prebent blades of current multi-MW turbines it produces an additional in-plane driving force. Two engineering models proposed by Li et al. (2022) — the smooth correction attributed to Madsen, driven by the area-averaged thrust coefficient, and the superposition of discrete semiinfinite vortex cylinders formed from the converged annulus inductions — were implemented independently in an open-source blade design and analysis code. For the NREL 5 MW reference rotor with an operational tip deflection of 6 m both models change thrust and power by less than 1 %, and their integral effects on the power coefficient agree to within 0.05 %. This agreement conceals a pronounced local disagreement: at mid-span the vortex-cylinder ur changes sign, following the non-monotonic radial derivative of the axial induction, and at the outermost section it grows logarithmically with the number of blade sections, ur ≃ 0.78 ln Nsec + 3.05 (m s−1), in close agreement with the analytic slope |γt tip|/2π = 0.82, reaching 6.6 m s−1 at Nsec = 70 against 2.3 m s−1 for the Madsen correction. We show that the integral agreement is systematic rather than coincidental, and that local load predictions near the tip of non-planar rotors are model- and grid-dependent and require an explicit, resolution-independent regularisation.
Laminar–turbulent transition on the suction surface of the LM45.3p blade ( $20\,\%$ thickness) was investigated using wall-resolved large eddy simulation (LES) at a chord Reynolds number of $Re_c=10^6$ and angle of attack $4.6^\circ$ . The effects of anisotropic free stream turbulence (FST) with intensities $TI=0\,\%$ – $7\,\%$ were examined, with integral length scales scaled down from atmospheric measurements. At $TI=0\,\%$ , a laminar separation bubble (LSB) forms and transition is initiated by Kelvin–Helmholtz vortices. At low FST levels ( $0\,\%\lt TI \leqslant 2.4\,\%$ ), robust streak growth via the lift-up mechanism suppresses the LSB, while transition dynamics shifts from two-dimensional Tollmien–Schlichting (TS) waves ( $TI=0.6\,\%$ ) to predominantly varicose inner and outer instabilities ( $TI=1.2\,\%$ and $2.4\,\%$ ) induced by the wall-normal shear and inflectional velocity profiles. The critical disturbance kinetic energy scales with $TI^{-1.80\pm 0.11}$ , compared with $TI^{-2.40}$ from Mack’s correlation. For $TI\geqslant 4.5\,\%$ , bypass transition dominates, driven by high-frequency boundary layer perturbations and streak breakdown via outer sinuous modes induced by the spanwise shear and inflectional velocity profiles. The scaling of streak amplitudes with $TI$ becomes sub-linear and spanwise non-uniformity characterises the turbulent breakdown. The critical disturbance kinetic energy reduces to $TI^{-0.90\pm 0.16}$ , marking a transition regime distinct from modal mechanisms. The onset of bypass transition ( $TI\approx 2.4\,\%{-}4.5\,\%$ ) aligns with prior studies of separated and flat-plate flows. A proposed turbulence spectrum cutoff links atmospheric measurements to wind tunnel data and Mack’s correlation, offering a framework for effective $TI$ estimation in practical environments.
This article reports about a wind-tunnel experiment carried out in the ONERA F2 low-speed wind tunnel on a model of the DU 97-W-300Mod airfoil designed for wind turbine application. The wind tunnel, the airfoil model, and experimental techniques used are presented, with special emphasis on the data processing and corrections required to derive airfoil forces and pressure distribution. To better document the flow physics at play, the results are illustrated by infrared thermography and surface oil flow visualization. The test allowed investigating Reynolds number effects between 1 and 3.8 millions. To ameliorate the understanding of the benefits and limitations of such airfoil testing, one section is devoted to the comparison of present results with previous experiments in other wind tunnels. Some of the difficulties arising in airfoil testing are evidenced and discussed to contribute to the improvement of test methods.
We designed 60% thick airfoil to improve the aerodynamic performance in the root region of wind turbine rotor blades, taking into account current constraints. After an extensive literature review and patent research, a design methodology (including the considerations of simple manufacturing) was set up, and extensive 2D- and 3D-CFD investigations with four codes (Xfoil, MSES, ANSYS fluent, and DLR-tau) were performed, including implementation inside a generic 10 MW test-blade (CIG10MW). Comparison with results from Blade Element Momentum (BEM) methods and the estimation of 3D effects due to the rotating blade were undertaken. One specific shape (with a pronounced flat-back) was selected and tested in the Deutsche WindGuard aeroacoustic Wind Tunnel (DWAA), in Bremerhaven, Germany. A total of 34 polars were measured, included two trailing edge shapes and aerodynamic devices such as vortex generators, gurney flaps, zig-zag tape, and a splitter plate. Considerable changes in lift and drag characteristics were observed due to the use of aerodynamic add-ons. With the studies presented here, we believe we have closed an important technological gap.
Within the framework of the fourth phase of the International Energy Agency (IEA) Wind Task 29, a large comparison exercise between measurements and aeroelastic simulations has been carried out featuring three simulation cases in axial, sheared and yawed inflow conditions. Results were obtained from more than 19 simulation tools originating from 12 institutes, ranging in fidelity from blade element momentum (BEM) to computational fluid dynamics (CFDs) and compared to state-of-the-art field measurements from the 2 MW DanAero turbine. More than 15 different variable types ranging from lifting-line variables to blade surface pressures, loads and velocities have been compared for the different conditions, resulting in over 250 comparison plots. The result is a unique insight into the current status and accuracy of rotor aerodynamic modeling. For axial flow conditions, a good agreement was found between the various code types, where a dedicated grid sensitivity study was necessary for the CFD simulations. However, compared to wind tunnel experiments on rotors featuring controlled conditions, it remains a challenge to achieve good agreement of absolute levels between simulations and measurements in the field. For sheared inflow conditions, uncertainties due to rotational and unsteady effects on airfoil data result in the CFD predictions standing out above the codes that need input of sectional airfoil data. However, it was demonstrated that using CFD-synthesized airfoil data is an effective means to bypass this shortcoming. For yawed flow conditions, it was observed that modeling of the skewed wake effect is still problematic for BEM codes where CFD and free vortex wake codes inherently model the underlying physics correctly. The next step is a comparison in turbulent inflow conditions, which is featured in IEA Wind Task 47. Doing this analysis in cooperation under the auspices of the IEA Wind Technology Collaboration Program (TCP) has led to many mutual benefits for the participants. The large size of the consortium brought ample manpower for the analysis where the learning process by combining several complementary experiences and modeling techniques gave valuable insights that could not be found when the analysis is carried out individually.
The interaction of several instabilities and the influence of free-stream turbulence on laminar-turbulent transition on a 20% thick wind-turbine blade section with a laminar separation bubble (LSB) are investigated with wall-resolved large-eddy simulations (LES). Turbulence intensities (TI) of 0%, 2.2%, 4.5%, 8.6%, and 15.6% at chord Reynolds number 105 are considered. Linear receptivity occurs for the most energetic disturbances; high-frequency perturbations are excited via non-linear mechanisms for TI≥8.6%. Unstable Tollmien–Schlichting (TS) waves appear in the inflectional flow region for TI≤4.5%, shifting to inviscid Kelvin–Helmholtz (KH) modes upon separation and forming spanwise rolls. Sub-harmonic secondary instability occurs for TI=0%, with rolls intertwining before transition. Streaks spanwise modulate the rolls and increase their growth rates with TI for TI≤4.5%, reducing separation and shifting transition upstream. The TI=4.5% case presents the highest perturbations, leading to the smallest LSB and most upstream transition. Earlier inception of TS/KH modes occurs on low-speed streaks, inducing premature transition. However, for TI=8.6%, the effect of the streaks is to stabilize the attached mean flow and front part of the LSB. This occurs due to the near-wall momentum deficit alleviation, leading to the transition delay and larger LSB than TI=4.5%. This also suppresses separation and completely stabilizes TS/KH modes for TI=15.6%. Linear stability theory predicts well the modal evolution for TI≤8.6%. Optimal perturbation analysis accurately computes the streak development upstream of the inflectional flow region but indicates higher amplification than LES downstream due to the capture of low-frequency, oblique modal instabilities from the LSB. Only low-amplitude [O(1%)] streaks displayed exponential growth in the LES since non-linearity precludes the appearance of these modes.
Among a few field experiments on wind turbines for analyzing laminar–turbulent boundary layer transition, the results obtained from the DAN-AERO and aerodynamic glove projects provide significant findings. The effect of inflow turbulence on boundary layer transition and the possible transition mechanisms on wind turbine blades are discussed and compared to CFD (computational fluid dynamics) simulations of increasing fidelity (Reynolds-averaged Navier–Stokes, RANS; unsteady Reynolds-averaged Navier–Stokes, URANS; and large-eddy simulations, LESs). From the experiments, it is found that the transition scenario changes even over a single revolution with bypass transition taking place under the influence of enhanced upstream turbulence, for example, such as that from wakes, while natural transition is observed in other instances under relatively low inflow turbulence conditions. This change from bypass to natural transition takes place at azimuthal angles directly outside the influence of the wake indicating a quick boundary layer recovery. The importance of a suitable choice of the amplification factor to be used within the eN method of transition detection is evident from both the RANS and URANS simulations. The URANS simulations which simultaneously check for natural and bypass transition match very well with the experiment. The LES predictions with anisotropic inflow turbulence show the shear-sheltering effect and a good agreement between the power spectral density plots from the experiment and simulation is found in case of bypass transition. A condition to easily distinguish the region of transition to turbulence based on the Reynolds shear stress is also observed. Overall, useful insights into the flow phenomena are obtained and a remarkably consistent set of conclusions can be drawn.
Wind energy plays a leading role as a renewable energy source, especially in countries with a low annual solar yield.
Free-stream turbulence (FST) alters the boundary layer of wind turbine blades, changing the hydrodynamic stability and near wake. Large-eddy simulations (LES) of a blade section with a laminar separation bubble for several turbulence intensities ( TI ) and a Reynolds number of 10 5 are performed. The effects of boundary-layer streaks generated by FST on Tollmien-Schlichting (TS) and Kelvin-Helmholtz (KH) instabilities are analyzed with a model based on the parabolized stability equations (PSE). Two competing effects on flow stability are identified. The spanwise-averaged mean-flow distortion stabilizes primary TS/KH modes for increasing TI . However, this contribution seems dominant only for TI ≥ 8.6%. For lower TI , the spanwise-oscillating distortion caused by streaks destabilizes the flow, and the growth rates of secondary modal instabilities increase with the streak amplitude. The destabilization occurs mainly at spanwise locations with negative streaks since the inflection point shifts away from the wall, enhancing inviscid instabilities. Inflection points in the spanwise direction formed by the streaks also contribute to the destabilization. The modal structures from PSE and LES agree. Finally, the trailing-edge near-wake coherent structures are more energetic for TI ≥ 8.6% due to the partial stabilization of modal instabilities, delaying the turbulent breakdown.
An approach to aerodynamically optimizing cycling posture and reducing drag in an Ironman (IM) event was elaborated. Therefore, four commonly used positions in cycling were investigated and simulated for a flow velocity of 10 m/s and yaw angles of 0–20° using OpenFoam-based Nabla Flow CFD simulation software software. A cyclist was scanned using an IPhone 12, and a special-purpose meshing software BLENDER was used. Significant differences were observed by changing and optimizing the cyclist’s posture. Aerodynamic drag coefficient (CdA) varies by more than a factor of 2, ranging from 0.214 to 0.450. Within a position, the CdA tends to increase slightly at yaw angles of 5–10° and decrease at higher yaw angles compared to a straight head wind, except for the time trial (TT) position. The results were applied to the IM Hawaii bike course (180 km), estimating a constant power output of 300 W. Including the wind distributions, two different bike split models for performance prediction were applied. Significant time saving of roughly 1 h was found. Finally, a machine learning approach to deduce 3D triangulation for specific body shapes from 2D pictures was tested.
You have accessEinführung in die WindenergietechnikApr 2022Einführung in die WindenergietechnikAlois P. SchaffarczykAlois P. Schaffarczykhttps://doi.org/10.3139/9783446473225.fmSectionsAboutPDF ToolsAdd to FavoritesDownload CitationTrack CitationsCopy LTI LinkPDF ShareFacebookTwitterEmailLinkedIn next chapter FiguresReferencesRelatedDetails 2022Pages: I-XXPrint ISBN: 978-3-446-46576-3eISBN: 978-3-446-47322-5 Copyright & Permissions© 2022 Carl Hanser Verlag GmbH & Co. KGPDF DownloadLoading ...