The main objective of the presented work is the validation of the simulation tool FAST.Farm for the calculation of power and structural loads in
single wake situations; the basis for the validation is the measurement database of the operating offshore wind farm alpha ventus. The approach is
described in detail and covers the calibration of the aeroelastic turbine model, transfer of environmental conditions to simulations, and comparison
between simulations and adequately filtered measurements. It is shown that FAST.Farm accurately predicts power and structural load distributions
over wind direction with discrepancies of less than 10 % for most of the cases compared to the measurements. Additionally, the frequency
response of the structure is investigated, and it is calculated by FAST.Farm in good agreement with the measurements. In general, the calculation of
fatigue loads is improved with a wake-added turbulence model added to FAST.Farm in the course of this study.
Introduction
Wind conditions inside a wind farm are strongly influenced by the interaction of individual turbines with the atmospheric boundary layer. In
particular, wake effects of upstream-located turbines affect the inflow conditions of downstream turbines. FAST.Farm is a new numerical tool developed
by the National Renewable Energy Laboratory (NREL), which simulates wake effects and predicts power output, as well as structural loads of turbines
within wind farms. It implements the dynamic wake meandering (DWM) model originally described by but with further model
advancements.
The original DWM model and modified versions of it were verified and validated in previous studies. published results of a
validation with respect to loads and power production for single and multiple wake scenarios. In an additional validation study,
took a closer look at tower loads and showed dependencies of wake loads on turbine spacing. Improvements to the DWM
model were suggested by with respect to atmospheric stability and regarding atmospheric shear and turbulence
build-up in a wind farm. In both studies, they validated their improvements by comparing them to power production measurements, as well as results from a
large eddy simulation (LES) code in terms of velocity deficit and turbulence intensity (TI) profiles. Recently, validation efforts of the DWM model were
conducted by incorporating flow field measurements from a lidar device (). derived new
calibration factors by using lidar measurements and validated their results in terms of power production and structural loads
().
FAST.Farm was validated with LES simulations in the prediction of wake characteristics, turbine power, and structural loads for a three turbine case by
. They showed that FAST.Farm calculates results in good agreement with the reference LES simulations for most analyzed
quantities. However, in low-ambient-turbulence conditions, higher differences were observed which were attributed to a missing wake-added turbulence
feature in FAST.Farm in this study. performed a validation of FAST.Farm against full-scale data of the turbine's
supervisory control and data acquisition (SCADA) system (generator power, rotor speed, blade pitch) of a five turbine configuration. Despite problems
with the used generic controller, FAST.Farm captured the trends of the measurements with good accuracy. In the Scaled Wind Farm Technology (SWiFT)
benchmark study by , FAST.Farm calculated flow characteristics of a single wake in good agreement with the other
considered simulation tools. Underperformance in capturing wakes compared to flow measurements was predominantly traced to inaccuracies in the inflow
modeling.
In the present study, we compare results from FAST.Farm against full-scale measurement data from the offshore wind farm alpha ventus. The overall
objectives are summarized as follows: (1) validating FAST.Farm for power and structural load predictions in single wake situations, (2) providing a detailed
path on how to perform load validation with field data of an operational wind farm, and (3) providing insights into load characteristics of an
offshore wind turbine subjected to single wake conditions. A one-to-one approach is followed for the validation concept. Here, environmental
conditions measured at the meteorological mast FINO1 are aggregated and directly fed as inputs into the simulations; this means that each measured
10 min event is represented in the simulations with its unique environmental parameter combination.
Methods and dataFAST.Farm overview
FAST.Farm is a multiphysics engineering software tool that accounts for wake interaction effects on turbine performance and structural loading within
wind farms. FAST.Farm is an extension of the NREL software OpenFAST, which solves the aero-hydro-servo-elasto dynamics of individual
turbines. FAST.Farm extends this analysis to include wake effects in wind farms. As in OpenFAST, rotor aerodynamics in FAST.Farm are modeled using the
blade-element-momentum (BEM) theory with options for advanced corrections, such as the inclusion of unsteady aerodynamics. Wake aerodynamics in
FAST.Farm are based on the DWM model but expand on it to address many limitations of past DWM implementations.
Using the DWM method, the wake deficit behind each turbine in the wind farm is computed quasi-steadily, with the wake evolution solved via the thin shear-layer approximation of the Navier–Stokes equations in axisymmetric coordinates. Turbulence closure is captured through an eddy-viscosity model,
including the influence of ambient turbulence and the wake shear layer. The wake deficit at the rotor is based on the low-pass time-filtered and
azimuth-averaged radially dependent thrust coefficient calculated by OpenFAST. Wake expansion in the pressure gradient zone is solved with a near-wake
correction. Each wake meanders due to large-scale turbulent structures using a three-dimensional passive tracer, with a meandering velocity based on
spatial averaging of the disturbed wind field. The disturbed wind field is calculated by superimposing the ambient turbulence and the individual wake
deficits from each rotor. When multiple wakes overlap, the superposition of axial wake deficits is based on a root sum square method.
Some of the unique innovations of FAST.Farm relative to DWM implementations in other simulation tools include
the improvement of wake advection, deflection, and merging;
the calibration of wake-related model parameters against results from high-fidelity LES simulations;
the ability to solve all wind turbines and the farm-wide disturbed wind field in parallel;
the optional inclusion of a wind-farm-wide super controller (not used in this paper); and
the optional inclusion of LES-generated ambient wind data (not used in this paper).
More information on the implementation and theory behind FAST.Farm is provided by , , and
. derived the calibration parameters of FAST.Farm through optimization with regard to LES simulations. These constants were directly used in this study.
Implementation of wake-added turbulence
Besides modeling wake deficit and wake meandering, the DWM model by includes modeling of wake-added turbulence. This term
describes the generation of turbulence behind a wind turbine rotor due to shear forces in the wake, as well as the breakdown of mainly the tip and
root vortices. The contribution of wake-added turbulence to the total turbulence level inside a turbine's wake is higher for low-ambient-turbulence
conditions (). Therefore, the inclusion of wake-added turbulence is especially important for offshore conditions
in which ambient turbulence levels are often low (i.e., less than 10 %). Preliminary results in the course of this work (see
Fig. ) supported this observation and led to the implementation of a wake-added turbulence feature in FAST.Farm which was not
present in previous versions. In summary, it is seen that the tower loads are strongly influenced by wake-added turbulence, and the corresponding
improvement of FAST.Farm leads to an increase in accuracy. In contrast, Fig. indicates that the blade loads are not very
sensitive to wake-added turbulence.
The implementation of wake-added turbulence in FAST.Farm presented herein follows mainly the approach by and
, which is included in the IEC 61400-1 standard (). In addition to the ambient turbulence domain,
it uses a new wake-added turbulence domain defined in the meandering frame of reference. This domain is generated with Mann's spectral turbulence
model (), defining turbulence as homogeneous and isotropic with a length scale that equals the rotor diameter. The domain
should have a fine spatial discretization to resolve the smaller turbulent scales of the wake-added turbulence.
The wake-added turbulence velocity components are scaled with the factor kmt defined by Eq. (). It consists of two
terms, which are influenced by (1) the quasi-steady wake deficit Udef(x,r̃)=U(x,r̃)/VDiskAvg, expressed with
the wake velocity U(x,r̃) that is normalized by the rotor disk-averaged ambient wind speed normal to the disk VDiskAvg, and
(2) the radial wake deficit gradient. The contributions of those terms are controlled with the empirical coefficients km1 and
km2. The factor kmt is dependent on the axial distance to the rotor x and the radial location r̃ normalized by the
rotor radius.
kmt(x,r̃)=km11-U(x,r̃)VDiskAvg+km2VDiskAvgδU(x,r̃)δr̃
Exemplary results from the calculation of wake-added turbulence from left to right: velocity deficit distribution, radial gradient of the velocity deficit, scaling factor kmt, and resulting distribution of the turbulence intensity from the wake-added turbulence TIWAT. Calculations are made for two ambient turbulence levels at 6 D behind the rotor. Simulations are performed with FAST.Farm in neutral atmospheric conditions with a mean wind speed of 6.5 ms-1 at hub height.
Exemplary distributions of the TI from the wake-added turbulence TIWAT and the corresponding velocity deficit profiles are
displayed in Fig. . The empirical constants involved in Eq. () were re-calibrated to
km1=1.44 and km2=0.84 by using the measurements. For the calibration, events with low ambient TI and close to full-wake
conditions (± 5∘) were utilized from the measurements as reference; then, the damage equivalent loads (DELs), as well as the frequency response
at the tower base, were compared. The recommended values by the IEC 61400-1 standard are km1=0.6 and km2=0.35. The reason for the
discrepancy between the IEC values and the re-calibrated values is not clear. However, the resulting turbulence levels from the wake-added turbulence
model seen in Fig. are similar to turbulence values found in the literature, e.g., from and
.
The implementation of wake-added turbulence includes three major additions to the FAST.Farm code:
A new instance of the FAST.Farm module InflowWind is initialized for the wake-added turbulence domain; it is reused for each turbine in the
simulation domain to ensure computational efficiency. The turbulent wind box that is created beforehand with the Technical University of Denmark (DTU) Mann turbulence generator
() is loaded into this instance of InflowWind. In the course of a FAST.Farm simulation, the wake-added turbulence wind
field is propagated with the ambient wind speed at hub height. The usage of a Mann turbulence field is mainly motivated by practical reasons with
regard to the implementation in InflowWind; when using a Mann turbulence field, a reference wind speed can be directly defined inside InflowWind,
which is eventually used to propagate the wake-added turbulence box downstream.
The scaling factor kmt is calculated in the meandering frame of reference inside FAST.Farm's “Wake Dynamics (WD)” module. The
calculation is based on the quasi-steady velocity deficit and its radial gradient that are already available inside the WD module.
In the module “Ambient Wind and Array Effects (AWAE)” of FAST.Farm, the velocities of the wake-added turbulence field are interpolated based
on their spatial location and scaled with the spatially interpolated value of the factor kmt. The resulting velocities are added to the
ambient wind vector via vector addition in the low- and high-resolution domains of FAST.Farm and transformed from the meandering frame of reference
to the fixed frame of reference.
Wind farm layout of alpha ventus including FINO1 met mast. AV1-6: type Senvion 5M. AV7-12: type Adwen AD 5-116. North is pointing up. The shaded triangle indicates the wind direction sector 257–287∘ seen from AV5.
Alpha ventus measurement database
The wind farm alpha ventus is located 45 km north of the German island of Borkum in the North Sea. It consists of 12 turbines with a rated
power of 5 MW, which is shown in Fig. . This study focuses on the turbines AV4 and AV5, which are Senvion 5M turbines with a rotor
diameter of 126 m and a hub height of 92 m. They are mounted on a jacket substructure and are located approximately 6.7 D
apart. Within the initiative Research at alpha ventus (RAVE, ) measurement data from both turbines have been acquired since
2011. For example, these data were used in load validation studies for freestream conditions by and
. For this work, we used data from the period January 2016–July 2018 because of good
availability and quality. In front of turbine AV4, the FINO1 met mast is located at a distance of approximately 3.2 D providing environmental
data.
Turbine measurements
The turbines are equipped with load sensors at various locations. Additionally, data from the turbine's SCADA system are available. The time
resolution of all sensors is 50 Hz. The following list explains the sensors used in this study and their calibration.
Generator power, generator speed, and blade-pitch angle measurements are directly taken from the SCADA system.
Nacelle yaw position data are also available through the SCADA system. Over longer time periods, a drift was observed in the data. This
was corrected by using nacelle rotation events and correlating the known tower-base strain gauge positions with the nacelle yaw signal. In this way,
sensor offset values were derived to make the data consistent over time.
Tower-base bending moments were calculated from strain gauges located above the transition piece. The strain gauges are placed at four locations
separated by 90∘ around the tower cylinder. By combining the strain gauge measurements with the nacelle's yaw position, fore–aft (FA) and
side–side (SS) bending moments were derived. Nacelle rotation events during turbine shutdown and calm wind conditions were used to determine
calibration factors in terms of slope and offset.
Blade-root bending moments in edgewise and flapwise directions are measured via four strain gauges placed near the blade root. They were
calibrated with rotor idling events during calm winds, as well as 10 min mean operational data. The measurements were made consistent over time by
adjusting slope and offset. Strain gauge signals were combined to reduce cross-talk effects.
Environmental data
Meteorological and sea conditions are measured at the FINO1 met mast and are available as 10 min statistics. Wind speed is measured with cup
anemometers at seven locations starting at 41.5 mheighta.s.l. (above sea level) and increasing in approx. 10 m
increments to 100 mheighta.s.l. Wind speed data are corrected for met mast shadow effects as explained by
. Vertical wind shear is described in terms of the power law exponent α, which is derived by fitting the
power law on the wind speed measurements from all available heights of the met mast.
Wind direction is taken from the wind vane located at 91.5 m height, with the correction given by . An additional
offset of +3∘ is applied to wind direction, which was derived by correlating the wake deficit of turbine AV4 with the measured wind
direction at FINO1.
Atmospheric stability is estimated by using the power law shear exponent α and applying the limits given in
Table . This simplified approach is motivated by and has the advantage of good sensor
availability. It is considered to be sufficient for this study because it only serves as refinement for environmental conditions, and specific analyses
of atmospheric stability are excluded.
Sea state is measured in terms of significant wave height and peak wave period with a directional waverider buoy (DWR-MkIII by Datawell).
Classification of atmospheric stability with power law exponent α. The given wind field parameters σk are dependent on stability and based on findings by , which is explained in Sect. .
Around alpha ventus, new wind farms have been commissioned over the years. This changes inflow conditions at alpha ventus compared to situations in which
no other wind farms were in close vicinity, as shown by . In this work, these effects are partly taken into account because
measured environmental conditions are directly transferred to simulations on a 10 min event basis.
Filtering approach
The measurement data are clustered in 10 min events for which the statistics are calculated to allow appropriate filtering. In particular, the
following filter criteria were applied:
Only events which have a data availability of more than 99 % are considered.
Wind direction is constrained from 257 to 287∘ to ensure that only wake effects of turbine AV4 affect AV5. At the boundaries of this
wind direction sector, nearly freestream conditions exist for both turbines; effects from wind farm blockage are ignored.
Both turbines are operated under normal conditions. For example, down-regulation, startup, or shutdown events are omitted.
No yaw action takes place during an event. Additionally, both turbines operate at similar yaw angles, allowing only events for which the
difference in the mean yaw position is less than 6∘.
Transfer of environmental conditions to simulations
An important part of a validation procedure is to ensure high-quality simulation inputs. In this study, turbulent wind fields were created with the
software TurbSim () fed by measured meteorological conditions at FINO1. The wind field generation with TurbSim applies
Taylor's frozen-turbulence assumption and propagates turbulent eddies with a constant velocity while neglecting the evolution of the turbulent eddies
in the downstream direction. According to , who compare results from FAST.Farm with results from LES, this assumption is
considered a reasonable simplification.
The von Karman turbulence spectrum is adapted with a modified wind profile: at heights where measurements at FINO1 exist, mean wind speed and TI were
directly used. For heights exceeding the met mast, a power law is utilized to extrapolate mean wind speed; for TI, the standard deviation value at the
top met mast position is kept constant over higher heights. It is assumed that the cup anemometer only measures turbulent fluctuations of the wind
speed component u, defined as the standard deviation σu. The relationship of turbulent fluctuations of the three wind speed components u,v, and w is then adjusted according to atmospheric stability (see Table ). It is based on findings by , who
investigated, among other things, the anisotropy of turbulence with respect to atmospheric stability. Their values of the anisotropy
parameter Γ are translated to standard deviations by using the relationships provided by .
Turbulence length scale L is approximated with Eq. (), as suggested by . Here, the TI measurement at
z= 90 m height and the power law shear coefficient α serve as inputs.
L≈zTIα
Coherence γ of the wind speed components k=u,v,w between points i and j is defined according to the model
in Eq. ():
γi,j=exp-ckfδu‾,
where f is the frequency, δ is the separation distance between points i and j, and u‾ is the mean wind speed of both
points. The decay coefficient ck is chosen according to and is dependent on
atmospheric stability. Note that for stable conditions, the value for coefficient ck that is dedicated to neutral conditions is used. The
coherence model is applied to the wind speed components individually and ignores correlations between different velocity components (e.g., u–v correlation).
Calibration of aeroelastic simulation model
The aeroelastic simulation model was generated with structural and aerodynamic information provided by the manufacturer. Simulations were performed
with the original turbine controller. In addition, a thorough calibration of the simulation model of both turbines AV4 and AV5 was performed to match
the turbine's measured load characteristics in the field as closely as possible. This involves the determination of structural damping of the first
tower eigenmodes in the FA and SS direction by analyzing turbine shutdown events.
Furthermore, imbalances in the rotating system were identified by looking at the frequency response during freestream events with low TI
conditions. The imbalances were introduced for each blade i as variation in blade mass BM,i, variation in blade flapwise BFK,i and
edgewise BEK,i stiffness, and blade-pitch offset BPO,i. A summary of the introduced imbalances in the simulation model is
listed in Table . Using these values, we were able to reproduce the turbines' frequency response from the field at the
desired sensor locations with satisfactory accuracy. It is noted that these imbalances do not necessarily reflect the real existing imbalances;
however, this approach was regarded as the best solution to overcome missing exact blade calibration measurements such as static blade deflection
tests.
Identified simulation model imbalances. Values without units are factors that are multiplied with the given parameter.
BM,1BM,2BM,3BFK,1BFK,2BFK,3BEK,1BEK,2BEK,3BPO,1BPO,2BPO,3AV410.9961110.95110.950∘-0.5∘+0.5∘AV51110.95110.9511-0.5∘+0.5∘+0.5∘Simulation setup and presentation of results
The information presented in the previous sections was used to set up the FAST.Farm simulations. One-to-one simulations were performed, in which the
measured environmental conditions of 10 min events are directly used as simulation inputs. Here, six random and uncorrelated turbulence
realizations of the wind field were created. Similarly, the sea state is modeled with six random realizations of the JONSWAP spectrum using the
measured significant wave height and peak wave period. The wake analysis focuses on two wind speed bins in the below-rated conditions: (1) wind speed
6.5–7.5 ms-1 labeled in the following with “I” and (2) wind speed 9.0–10.0 ms-1 labeled with ”II”. These bins were chosen
because they imply high rotor thrust values and hence strong wake effects. In total, 1014 FAST.Farm simulations were run for wind speed bin I and 1072
for wind speed bin II.
The numerical setup of FAST.Farm is guided by the recommendations given in . The domain size in the longitudinal direction
is X= 1890 m. In the lateral direction, the largest shift of both turbines is 220 m, and it occurs at the considered wind
directions of 257 and 287∘. By considering the rotor radius and reserving additional space (2 D) for the wake meandering in both
directions, the domain is set to Y= 850 m. The domain size in the vertical direction is set to Z= 300 m, which includes
a safety margin of ≈ 1.2 D above the rotor tip for vertical wake meandering. The wakes are modeled up to a downstream distance of
8 D by using wake planes with a spatial resolution of 5 m in the radial direction.
The TurbSim wind fields for the ambient turbulence and the low-resolution domain for resolving the wake meandering in FAST.Farm have a spatial
resolution of dy=dz=6m in the lateral and vertical directions. The spatial resolution for the wake-added turbulence
domain and the high-resolution domain around the turbines is dy=dz=3m. The time step for the low-resolution domain equals
dtlow=2s and dthigh=0.2s for the high-resolution domain, which is also used for the
TurbSim wind fields for the ambient turbulence and the Mann wind fields for the wake-added turbulence. At the beginning of each simulation, a
transient period of 400 s is removed to allow the wakes to develop and to damp initial oscillations. The TurbSim wind fields are 600 s
long and set to be periodic in order to achieve a total simulation time of 1000 s. The Mann wind fields for the wake-added turbulence cover
3 D in the longitudinal direction; turbulence is reused when points outside this box are requested.
The results in Sects. and are presented in terms of statistics of measured quantities, which are plotted
against the wind speed and wind direction. The data are clustered in wind speed bins of the size 0.5 ms-1 and wind direction bins of the
size 2.5∘. For each bin, the mean value plus the 15th–85th percentile range is shown. DELs are calculated according to Eq. ():
DEL=∑i=1nSimNeq1m,
where Si refers to the range of a load cycle and Neq is the number of equivalent cycles (in this work: Neq=600 for
obtaining the 1 Hz DEL). The Wöhler exponents are defined as m=4 in the case of steel (tower) and m=10 for composite materials
(blades). The DEL is not corrected for mean load effects.
To protect proprietary data, all results are normalized by either user-defined values or by values of the freestream turbine, which is indicated in
the figure legend. Measurements are labeled with “AV4” and “AV5”, whereas simulations are named “FFarm4” and “FFarm5” making reference to the
turbine numbers 4 and 5.
Initial results
A check of the turbine characteristics during freestream conditions, as well as an analysis of the environmental conditions, is presented before taking
wake effects into account in Sect. .
Comparison of freestream (wind direction 240–257∘) characteristics between measurements (AV4, AV5) and simulations (OpenFAST: results are shown with calibrations based on AV4). Statistics are shown as mean values per wind bin. The shaded area indicates the 15th and 85th percentiles. Patches in grey show the wake bins I and II analyzed in Sect. .
Investigation of turbine characteristics in freestream conditions
Figure illustrates results of the comparison between simulations and field measurements for close-to-freestream
conditions. This freestream sector of 240–257∘ was detected beforehand by correlating loads of turbine AV5 with the wind direction. In total,
2980 simulations were conducted with the software OpenFAST, which employs the same aeroelastic model as described for FAST.Farm but without wind-farm-wide effects, enabling faster simulation times. Although an individual calibration of both turbines AV4 and AV5 was conducted, the differences in
the considered aggregated load quantities is negligible; hence, only the results of the AV4 calibration are shown in Fig. .
For the operational sensors of generator power (Fig. a) and blade-pitch angle (Fig. b), a
discrepancy of less than 5 % is found between the data of turbines AV4, AV5, and simulations described as OpenFAST. The DEL computed for the FA
bending moment at the tower base (Fig. c) shows a close match with a difference of less than 7 % in below-rated conditions
(4–12 ms-1). In the above-rated conditions, simulations agree well with measurements from turbine AV4; the loads of turbine AV5 are up
to 20 % higher compared to turbine AV4. This difference is most likely related to an imbalance in the rotor system during blade-pitch
actuation. Loads at the blade root are compared by means of DEL of the bending moment in the flapwise direction in Fig. d. It
is observed that measurements of both turbines and simulations match each other well with differences of less than 10 % for wind speeds up
to 16 ms-1.
Overall, the simulations can predict the measured load quantities in freestream conditions with high accuracy. This indicates that the aeroelastic
simulation model is set up appropriately and that the meteorological conditions are transferred into realistic wind fields. The exact representation
of hydrodynamic excitation in terms of wave loads is considered of less importance because the substructure is quite rigid and only sensors above the
seawater level are taken into account in this study.
Turbulence intensity measured at 90 m height at FINO1. Statistics are shown as mean values per wind bin. The shaded area indicates the 15th and 85th percentiles. Scatter shows TI of single 10 min events. AS = atmospheric stability.
Power law shear exponent α derived from FINO1 measurements. Statistics are shown as mean values per wind bin. The shaded area indicates the 15th and 85th percentiles. Scatter shows α of single 10 min events. AS = atmospheric stability.
Distribution of environmental conditions
For the analyzed wind speed bins, the TI distribution over wind direction is shown in Fig. . The mean TI for wind speed bin I
is calculated as ≈ 6.4 % and for wind speed bin II as ≈ 5.9 %. Highest TI values are found in unstable atmospheric conditions,
whereas the lowest TI values occur in stable conditions. Especially in wind speed bin II, more events are found towards wind directions from the
southwest, which is the predominant wind direction. A more uniform event distribution with regard to the wind direction is found for wind speed bin
I. Figure depicts the distribution of the power law shear exponent α over wind direction. It can be seen that wind
speed bin II contains higher α values on average compared to wind speed bin I.
Comparison of generator power in wake conditions between measurements (AV4, AV5) and simulations (FFarm4, FFarm5). Statistics are shown as mean values per wind direction bin. The shaded area indicates the 15th and 85th percentiles.
Results of wake validationTurbine performance: statistics
Results of the mean generator power prediction are presented in Fig. . They show a maximum power loss of ≈ 48 % in
full-wake conditions for wind speed bin I. At nearly 43 %, power loss is slightly less for wind speed bin II due to a decreased rotor thrust
coefficient compared to wind speed bin I. Wake effects are visible in the turbine power generation over a wind direction sector of
25∘. FAST.Farm is able to predict the width and depth of the power deficit with high accuracy compared to the measurements for wind speed
bin I (Fig. a and b). For wind speed bin II, deviations of 5–10 % are observed in Fig. c. Relative
comparisons are shown in Fig. b and d, where results of the waked turbine are divided by the freestream results of simulations and
measurements correspondingly. It is observed that relative plotting produces a better match of FAST.Farm with the measurements. This is due to the
decreased uncertainties potentially arising from the turbine model. Similar levels in the scatter of events indicated by the error range are found for
simulations and measurements.
Probability density function (PDF) of generator speed for measurements (AV4, AV5) and simulations (FFarm4, FFarm5) in the wind direction sector 270.75–273.25∘. Thick lines show the mean of the PDFs of all considered events. The shaded area indicates the 15th and 85th percentiles.
Turbine performance: detailed results
A more detailed analysis of the turbine performance is shown in Fig. by plotting the probability density function (PDF) of the
generator speed. The PDF was calculated for each 10 min event in the wind sector corresponding to full-wake conditions. Afterwards, mean values
with the error range expressed as the 15th and 85th percentiles across the PDFs of all events were derived. In both wind speed bins, a reduction in
generator speed is found for the waked turbine, caused by the wind speed deficit from the upstream turbine. In wind speed bin I, the downstream
turbine operates near the cut-in wind speed, which is indicated by the peak around the normalized generator speed of 0.6 in
Fig. a. From Fig. b, it can be seen that a wider range of generator speeds is covered by the waked turbine
compared to the freestream turbine. This can be related to a varied operational point due to the wake deficit. Additionally, it can be partly
attributed to the increased turbulence in the wake originating from wake-added turbulence, as well as wake meandering. In both wind speed bins,
FAST.Farm predicts the distributions from the measurements with high accuracy.
Comparison of the blade-root bending moment in the flapwise direction in waked conditions between measurements (AV4, AV5) and simulations (FFarm4, FFarm5). Statistics are shown as mean values per wind direction bin. The shaded area indicates the 15th and 85th percentiles.
Structural loads: statistics
Figure shows results of the fatigue loads expressed as DEL of the blade-root bending moment in the flapwise direction. By
comparing the two wind speed bins, different load distributions over wind direction for the waked turbine are observed. For wind speed bin I
(Fig. a), a dip in the DEL for the downstream turbine occurs around full-wake conditions at 272∘; this is not visible
for wind speed bin II (Fig. b). Influencing factors on the load distribution for the waked turbine are the wind direction and
connected mean wake position, magnitude of wake meandering, ambient wind conditions, and operational point of the turbine. Varying combinations of
these effects lead to different load distributions. With the chosen one-to-one simulation approach, the aim was to reduce the uncertainty arising from
the different combinations. Hence, it is seen for both wind speed bins that FAST.Farm agrees well with the measurements and predicts the increase in
loads and trends over wind direction with good accuracy. Overall, a load increase of approximately a factor of 1.7 (wind speed bin I) and factor of 2.0 (wind
speed bin II) are identified for the waked turbine.
Comparison of the tower-base bending moment in the fore–aft direction in waked conditions between measurements (AV4, AV5) and simulations (FFarm4, FFarm5). Statistics are shown as mean values per wind direction bin. The shaded area indicates the 15th and 85th percentiles.
Figure displays the DELs of the tower-base bending moment in the FA direction. Contrary to the fatigue loads at the
blade root, a higher increase in loads for the waked turbine is observed for wind speed bin I compared to wind speed bin II. In particular, DELs are
increased by a factor of ≈ 2.4 (wind speed bin I) and factor of ≈ 2.0 (wind speed bin II) for the waked turbine compared to the DELs of the
freestream turbine. FAST.Farm produces results in good agreement with the measurements in terms of magnitude and wind direction dependency. For most
of the considered wind directions, the discrepancy in the mean value per bin between FAST.Farm and the measurements is less than 10 %; in some
wind directions, the difference is increased to ≈ 25 %. The uncertainty range per bin indicated by the percentile range is predicted by
FAST.Farm in good agreement with the measurements.
Comparison of power spectral densities (PSDs) between measurements (AV4, AV5) and simulations (FFarm4, FFarm5) for wind speed bin I (wind direction 270.75–273.25∘). Thick lines show the mean of the PSDs for all considered events. The shaded area indicates the 15th and 85th percentiles.
Structural loads: frequency response
Figure depicts the frequency response of the structure at the tower base and blade root. Only events during nearly
full-wake conditions (wind direction 270.75–273.25∘) in wind speed bin I were considered. For each 10 min event, the power spectral
density (PSD) was calculated. Then mean values for each frequency with the corresponding uncertainty range expressed as the 15th and 85th percentiles were
determined. They are shown for the blade-root bending moment in the flapwise direction (Fig. a) and the tower-base
bending moment in the FA direction (Fig. b).
In the case of the blade (Fig. a), an increase in the energy at the blade passing frequency 1P is observed for the waked
turbine compared to the freestream turbine. It comes with a reduction in the blade passing frequency (freestream:
1P ≈ 0.15 Hz; wake: 1P ≈ 0.13 Hz) due to the reduced wind speed inside the wake and consequently a reduction in
rotor speed. The magnitude of the first blade passing frequency of turbine AV5 is predicted by FAST.Farm to be a factor of 3 higher compared to the
measurements. In contrast, the excitation of the second blade passing frequency 2P indicated by the peaks between 0.2–0.3 Hz is higher in the
measurements. A possible explanation is the modeling of the wake, which has a Gaussian shape in FAST.Farm. In reality, the wake is more likely to be
distorted, leading to smoother transitions to undisturbed winds. Similar observations are made by , who compare FAST.Farm
with LES simulations. However, more detailed analyses are required, e.g., using LES simulations, to derive a conclusive explanation for the observed characteristics in the
frequency response.
The signal at the tower base (Fig. b) reveals a strong excitation of the first global mode at around 0.3 Hz
for the waked turbine. This is observed in both the measurements and FAST.Farm, whose prediction of the peak is 15 % below the measurements. It
was found that especially the inclusion of wake-added turbulence in the FAST.Farm simulations has an intensifying effect on the peak of the first
global mode. Higher energy content is also seen in the frequency range 0.35–0.4 Hz for the waked turbine compared to the freestream
turbine. Although FAST.Farm captures an increase in the energy content compared to the freestream turbine, it underestimates the energy level detected
in the measurements by 50 % in this frequency range.
Discussion
The present investigation concentrates on two wind speed bins in below-rated conditions. This choice is motivated by the connected high rotor thrust
conditions and hence strong wake effects. Moreover, analyses of quantities over wind direction are enabled. Another reason is that the behavior of
both turbines AV4 and AV5 is comparable in the measurements, whereas in above-rated conditions, differences occur even under ambient inflow for both
turbines (see Fig. c). For the analysis of structural loads, we focus on sensors and directions that are mainly affected by
the change in turbulence characteristics in the wake, i.e., tower-base FA bending moment and blade-root flapwise bending moment.
A crucial part of this load validation is the generation of adequate wind fields representing the environmental conditions at alpha
ventus. Especially coherence and turbulence scale have an influence on wake meandering magnitude, which in turn affects the loads of the downstream
turbine (see also and ). Unstable and neutral atmospheric conditions imply greater turbulent
length scales and larger coherent turbulent structures than stable conditions. This leads to higher wake meandering magnitudes and higher load levels
for the downstream turbine. For the loads at the blade root, adequate capturing of wake meandering is most important, whereas for the loads at the
tower-base, both wake meandering and wake-added turbulence must be modeled. We observed that in the simulations, a direct relationship between ambient
TI conditions and wake loads exists. Consequently, higher ambient TI values lead to higher loads at the downstream turbine. In the measurements, this
relationship holds true, but it was also found that low ambient TI conditions can lead to high wake loads. This shows that there is some uncertainty in
modeling the environmental conditions and wake features that should be investigated in future.
In offshore full-scale load validation, there are many potential sources of uncertainty. Starting with the modeling of environmental conditions, we
aimed to minimize those uncertainties by making use of findings from previous research, which is available for the site. However, there are limits in
the methods used. For example, in the coherence model of the wind field, there is no directional dependency of coherence considered, and coherence is
only dependent on the velocity components u,v, and w. In the considered period of measurements and wind direction sector, alpha ventus operates in the
wake of the wind farm “Trianel Borkum I”, which is located ≈ 6.5 km east. The flow structures evolving from this farm-wake are
likely to be different from ideal freestream conditions, for which the wind field generation method was originally derived. Overall, to reduce the
input uncertainties, we followed a one-to-one simulation approach in which the measured environmental conditions are utilized directly as simulation
inputs.
Conclusions
The simulation tool FAST.Farm was validated for the prediction of power output and structural loads in single wake conditions with respect to
measurement data from the offshore wind farm alpha ventus. In addition, a wake-added turbulence model was implemented into FAST.Farm which is needed
to calculate the small-scale turbulence in the wake; this is considered of importance especially for low-ambient-turbulence intensity conditions and
for tower-base loading. It was shown that FAST.Farm predicts the mean power deficit with high accuracy compared to the measurements. Additionally, the
probability density function of generator speed was calculated in strong agreement with the measurements for the freestream and downstream
turbine. Fatigue loads were analyzed in terms of damage equivalent loads of the bending moments at the blade root in the flapwise direction and the
tower base in the fore–aft direction. Distributions over wind direction show a good match between simulations and measurements with deviations of less
than 10 % for most of the investigated wind directions.
More detailed insights into the aforementioned structural load quantities were provided by power spectral density (PSD) analyses. They show that
FAST.Farm calculates trends in the structural response in good agreement with the measurements in the frequency domain. In particular, excitation at
the tower base of the waked turbine is reproduced with FAST.Farm, which can be attributed to the wake-added turbulence feature in FAST.Farm included in
the course of this study. However, by looking at the PSD at blade root, it is indicated that not all phenomena are captured sufficiently by FAST.Farm,
leaving room for further improvements.
It was demonstrated that the proposed one-to-one simulation approach works well for the validation in offshore single wake conditions. It is concluded
that the calibration of the aeroelastic model with respect to imbalances, as well as appropriate transfer of environmental conditions to simulations, is
important. Here, a differentiation of atmospheric stability helps to refine simulation inputs such as coherence in the wind field but also indicates
that especially stable atmospheric conditions remain challenging to model for capturing the loads of a waked turbine.
Influence of the wake-added turbulence feature in FAST.Farm
Comparison of FAST.Farm simulations with activated (WAT = true) and deactivated (WAT = false) wake-added turbulence with respect to the measurements (AV5) for wind speeds of 7.5–8.5 ms-1. Statistics are shown as mean values per wind direction bin. The shaded area indicates the 15th and 85th percentiles.
Comparison of power spectral densities (PSD) between measurements (AV5) and FAST.Farm simulations with activated (WAT = true) and deactivated (WAT = false) wake-added turbulence for wind speeds of 7.5–8.5 ms-1 (wind direction 270.75–273.25∘). Thick lines show the mean of the PSDs for all considered events. The shaded area indicates the 15th and 85th percentiles.
Code and data availability
The version of FAST.Farm including the wake-added turbulence model which was used to create the results presented in this paper is different than the official version released by NREL but is available at (10.5281/zenodo.4733599). NREL has plans to incorporate an updated version of the present implementation in a future release of FAST.Farm.
The measurement data used in this study can be accessed via the “Bundesamt für Seeschifffahrt und Hydrographie” (; https://www.bsh.de/EN/TOPICS/Monitoring_systems/MARNET_monitoring_network/RAVE/rave_node.html, last access: 3 May 2021). A data usage agreement must be signed with the BSH in advance.
Author contributions
MK and JJ created with the help of VP the conceptualization of this project. MK processed, aggregated, and filtered the measurement data. MK added the wake-added turbulence feature to FAST.Farm and set up and performed all simulations. MK wrote the manuscript with support from JJ. JJ, VP, and PWC provided guidance for the research and reviewed the manuscript.
Competing interests
The contact author has declared that neither they nor their co-authors have any competing interests.
Disclaimer
The Alliance for Sustainable Energy, LLC (Alliance) is the manager and operator of NREL. NREL is a national laboratory of the US Department of Energy, Office of Energy Efficiency and Renewable Energy. This work was authored in part by the Alliance and supported by the U.S. Department of Energy under contract no. DE-AC36-08GO28308. Funding was provided by the Wind and Water Power Program. The views expressed in the article do not necessarily represent the views of the US Department of Energy or the US government. The US government retains – and the publisher, by accepting the article for publication, acknowledges that the US government retains – a nonexclusive, paid-up, irrevocable, worldwide license to publish or reproduce the published form of this work, or allow others to do so, for US government purposes.
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Acknowledgements
We would like to thank Senvion and DNV GL for providing the turbine measurements. Additional acknowledgements go to the Federal Maritime and Hydrographic Agency (BSH) for providing access to the FINO1 database and UL DEWI for providing wind measurement corrections of the FINO1 met mast.
Financial support
The research leading to these results has received funding from the German Federal Ministry for Economic Affairs and Energy (BMWi) in the framework of the national research project RAVE-OWP Control (grant no. FKZ 0324131B) and is part of the research done in the WindForS research cluster.
This open-access publication was funded by the University of Stuttgart.
Review statement
This paper was edited by Amir R. Nejad and reviewed by Valentin Chabaud and Adam Wise.
ReferencesBSH (Bundesamt für Seeschifffahrt und Hydrographie):
RAVE measurement database [data set],
available at: https://www.bsh.de/EN/TOPICS/Monitoring_systems/MARNET_monitoring_network/RAVE/rave_node.html (last access: 3 May 2021), 2021.Conti, D., Dimitrov, N., Peña, A., and Herges, T.: Calibration and validation of the Dynamic Wake Meandering model Part I: Bayesian estimation of model parameters using SpinnerLidar-derived wake characteristics, Wind Energ. Sci. Discuss. [preprint], 10.5194/wes-2020-135, in review, 2021.Davenport, A. G.:
The spectrum of horizontal gustiness near the ground in high winds,
Q. J. Roy. Meteor. Soc.,
87, 194–211, 10.1002/qj.49708737208, 1961.Doubrawa, P., Annoni, J. R., and Jonkman, J. M.:
Optimization-Based Calibration of FAST.Farm Parameters against Large-Eddy Simulations,
in: 2018 Wind Energy Symposium, American Institute of Aeronautics and Astronautics, 8–12 January 2018, Kissimmee, Florida, 10.2514/6.2018-0512, 2018.Doubrawa, P., Quon, E. W., Martinez-Tossas, L. A., Shaler, K., Debnath, M., Hamilton, N., Herges, T. G., Maniaci, D., Kelley, C. L., Hsieh, A. S., Blaylock, M. L., Laan, P., Andersen, S. J., Krueger, S., Cathelain, M., Schlez, W., Jonkman, J., Branlard, E., Steinfeld, G., Schmidt, S., Blondel, F., Lukassen, L. J., and Moriarty, P.:
Multimodel validation of single wakes in neutral and stratified atmospheric conditions,
Wind Energy,
23, 2027–2055, 10.1002/we.2543, 2020.DTU Wind Energy:
Mann 64bit turbulence generator,
available at: https://www.hawc2.dk/Download/Pre-processing-tools/Mann-64bit-turbulence-generator (last access: 13 April 2021), 2021.
IEC:
61400-1 Ed. 4: Wind energy generation systems – Part 1: Design requirements, Standard,
International Electrotechnical Commission, Geneva, Switzerland, 2019.
Jonkman, B.:
TurbSim User's Guide v1.50, Tech. Rep. NREL/TP-500-46198,
National Renewable Energy Laboratory, Golden, CO, 2009.Jonkman, B., Mudafort, R. M., Platt, A., Sprague, M., Vijayakumar, G., Buhl, M., Ananthan, S., Masciola, M., Schmidt, M. J., Branlard, E., jjonkman, Rood, J., rdamiani, ashesh2512, Bendl, K., pschuenemann, psakievich, KUSUNO, N., Martinez, T., robynnemurrayNREL, Mertz, B., Mylonas, C., IAbda, Rinker, J., Friedrich, M., Bachant, P., Hammond, R., lapadron, lssraman, mattEhall, and Kretschmer, M.:
MatthiasCK/OpenFAST: FAST.Farm with wake-added turbulence,
Zenodo [code], 10.5281/zenodo.4733599, 2021.
Jonkman, J. and Shaler, K.:
FAST.Farm User's Guide and Theory Manual, Tech. Rep. NREL/TP-5000-78485,
National Renewable Energy Laboratory, Golden, CO, 2021.Jonkman, J., Annoni, J., Hayman, G., Jonkman, B., and Purkayastha, A.:
Development of FAST.Farm: A New Multiphysics Engineering Tool for Wind Farm Design and Analysis,
in: AIAA Scitech 2017 Forum, 9–13 January 2017, Grapevine, Texas, 10.2514/6.2017-0454, 2017.
Kaufer, D. and Cheng, P. W.: Validation of an integrated simulation method with high resolution load measurements of the offshorewind turbine REpower 5M at Alpha Ventus, Journal of Ocean and Wind Energy, 1, 30–40, 2014.Keck, R.-E., de Maré, M. T., Churchfield, M. J., Lee, S., Larsen, G. C., and Aagaard Madsen, H.:
On atmospheric stability in the dynamic wake meandering model,
Wind Energy,
17, 1689–1710, 10.1002/we.1662, 2014.Keck, R. E., De Maré, M., Churchfield, M. J., Lee, S., Larsen, G., and Madsen, H. A.:
Two improvements to the dynamic wake meandering model: Including the effects of atmospheric shear on wake turbulence and incorporating turbulence build-up in a row of wind turbines,
Wind Energy,
18, 111–132, 10.1002/we.1686, 2015.Kelly, M.: From standard wind measurements to spectral characterization: turbulence length scale and distribution, Wind Energ. Sci., 3, 533–543, 10.5194/wes-3-533-2018, 2018.Larsen, G. C., Madsen, H. A., Thomsen, K., and Larsen, T. J.:
Wake meandering: a pragmatic approach,
Wind Energy,
11, 377–395, 10.1002/we.267, 2008.Larsen, T., Larsen, G., Pedersen, M., Enevoldsen, K., and Madsen, H.:
Validation of the Dynamic Wake Meander model with focus on tower loads,
J. Phys. Conf. Ser.,
854, 012027, 10.1088/1742-6596/854/1/012027, 2017.Larsen, T. J., Aagaard Madsen, H., Larsen, G. C., and Hansen, K. S.:
Validation of the dynamic wake meander model for loads and power production in the Egmond aan Zee wind farm,
Wind Energy,
16, 605–624, 10.1002/we.1563, 2013.Madsen, H. A., Larsen, G. C., Larsen, T. J., Troldborg, N., and Mikkelsen, R.:
Calibration and Validation of the Dynamic Wake Meandering Model for Implementation in an Aeroelastic Code,
J. Sol. Energy Eng.,
132, 041014, 10.1115/1.4002555, 2010.Mann, J.:
The spatial structure of neutral atmospheric surface-layer turbulence,
J. Fluid Mech.,
273, 141–168, 10.1017/S0022112094001886, 1994.Mann, J.:
Wind field simulation,
Probabilist. Eng. Mech.,
13, 269–282, 10.1016/S0266-8920(97)00036-2, 1998.Nybø, A., Nielsen, F. G., Reuder, J., Churchfield, M. J., and Godvik, M.:
Evaluation of different wind fields for the investigation of the dynamic response of offshore wind turbines,
Wind Energy,
23, 1810–1830, 10.1002/we.2518, 2020.Peña, A., Gryning, S.-E., and Mann, J.:
On the length-scale of the wind profile,
Q. J. Roy. Meteor. Soc.,
136, 2119–2131, 10.1002/qj.714, 2010.Pettas, V., Kretschmer, M., Clifton, A., and Cheng, P. W.: On the effects of inter-farm interactions at the offshore wind farm Alpha Ventus, Wind Energ. Sci. Discuss. [preprint], 10.5194/wes-2021-50, in review, 2021.
Popko, W., Robertson, A., Jonkman, J., Wendt, F., Thomas, P., Müller, K., Kretschmer, M., Hagen, T. R., Galinos, C., Le Dreff, J.-B., Gilbert, P., Auriac, B., Oh, S., Qvist, J., Sørum, S. H., Suja-Thauvin, L., Shin, H., Molins, C., Trubat, P., Bonnet, P., Bergua, R., Wang, K., Fu, P., Cai, J., Cai, Z., Alexandre, A., and Harries, R.:
Validation of Numerical Models of the Offshore Wind Turbine From the Alpha Ventus Wind Farm Against Full-Scale Measurements Within OC5 Phase III,
J. Offshore Mech. Arct.,
143, 012002, 10.1115/1.4047378, 2021.RAVE (German research initiative Research at alpha ventus RAVE):
available at: https://www.rave-offshore.de/en/ (last access: 13 January 2021), 2021.Reinwardt, I., Schilling, L., Dalhoff, P., Steudel, D., and Breuer, M.:
Dynamic wake meandering model calibration using nacelle-mounted lidar systems, Wind Energ. Sci., 5, 775–792, 10.5194/wes-5-775-2020, 2020.Reinwardt, I., Schilling, L., Steudel, D., Dimitrov, N., Dalhoff, P., and Breuer, M.: Validation of the dynamic wake meandering model with respect to loads and power production, Wind Energ. Sci., 6, 441–460, 10.5194/wes-6-441-2021, 2021.Shaler, K. and Jonkman, J.:
FAST.Farm development and validation of structural load prediction against large eddy simulations,
Wind Energy,
24, 428–449, 10.1002/we.2581, 2020.Shaler, K., Jonkman, J., Doubrawa, P., and Hamilton, N.:
FAST.Farm Response to Varying Wind Inflow Techniques: Preprint,
in: AIAA Scitech 2019 Forum, 7–11 January 2019, San Diego, CA, 10.2514/6.2019-2086, 2019a.Shaler, K., Jonkman, J., and Hamilton, N.:
Effects of Inflow Spatiotemporal Discretization on Wake Meandering and Turbine Structural Response using FAST.Farm,
J. Phys. Conf. Ser.,
1256, 012023, 10.1088/1742-6596/1256/1/012023, 2019b.Shaler, K., Debnath, M., and Jonkman, J.:
Validation of FAST.Farm Against Full-Scale Turbine SCADA Data for a Small Wind Farm,
J. Phys. Conf. Ser.,
1618, 062061, 10.1088/1742-6596/1618/6/062061, 2020.
Westerhellweg, A., Canadillas, B., Beeken, A., and Neumann, T.:
One year of LiDAR measurements at FINO1-Platform: Comparison and verification to met-mast data,
in: Proceedings of 10th German wind energy conference DEWEK, 17–18 November 2010, Bremen, 2010.
Westerhellweg, A., Riedel, V., and Neumann, T.:
Comparison of Lidar-and UAM-based offshore mast effect corrections,
in: Proceedings of EWEA, 14–17 March 2011, Brussels, 2011.Westerhellweg, A., Cañadillas, B., Kinder, F., and Neumann, T.:
Wake Measurements at alpha ventus – Dependency on Stability and Turbulence Intensity,
J. Phys. Conf. Ser.,
555, 012106, 10.1088/1742-6596/555/1/012106, 2014.Wise, A. S. and Bachynski, E. E.:
Wake meandering effects on floating wind turbines,
Wind Energy, 23, 1266–1285, 10.1002/we.2485, 2020.