D.J.N. Allaerts
Please Note
22 records found
1
The filtered lifting line theory presents a continuous form of the inviscid momentum equations of flow over a lifting device, such as a wing or rotor blade, using body forces without mathematical singularities. This theory is also consistent with an actuator line representation of a lifting device. In this work, we present a reformulation of the equations in terms of the local flow angle along the line, which allows solving the stand-alone equations using multivariate root-finding algorithms. This approach can be used to obtain a fast, computationally inexpensive solution of the loading distribution along a wing without the need to perform computational fluid dynamic simulations. We study the requirements in terms of resolution in the spanwise direction and establish the criteria for spacing and minimum amount of points required along the blade to obtain converged solutions. The solutions are compared to results from large-eddy simulations, and we observed excellent agreement with less than a percent difference in quantities along the blade between the methods.
Wind farm flow control has been a key research focus in recent years, driven by the idea that a collectively operating wind farm can outperform individually controlled turbines. Control strategies are predominantly applied in an open-loop manner, where the current flow conditions are used to look up precomputed steady-state set points. Closed-loop control approaches, on the other hand, take measurements from the farm into account and optimize their set points online, which makes them more flexible and resilient. This paper introduces a closed-loop model-predictive wind farm controller using the dynamic engineering model FLORIDyn to maximize the energy generated by a ten-turbine wind farm. The framework consists of an Ensemble Kalman Filter to continuously correct the flow field estimate, as well as a novel optimization strategy. To this end, the paper discusses two dynamic ways to maximize the farm energy and compares this to the current look-up table industry standard. The framework relies solely on turbine measurements without using a flow field preview. In a 3-h case study with time-varying conditions, the derived controllers achieve an overall energy gain of 3% to 4.4% with noise-free wind direction measurements. If disturbed and biased measurements are used, this performance decreases to 1.9% to 3% over the greedy control baseline with the same measurements. The comparison to look-up table controllers shows that the closed-loop framework performance is more robust to disturbed measurements but can only match the performance in noise-free conditions.
In this work, we only consider linearly stratified conditions (i.e., no inversion layer), thereby focusing on internal gravity waves in the free atmosphere and their reflections from the domain boundaries. This type of flow is governed by a single Froude number, which dictates most of the internal wave properties, such as wavelength, amplitude, and direction. This in turn will dictate the optimum domain size and Rayleigh damping layer setup. We find the effective horizontal and vertical wavelengths, (the representative wavelengths of the entire wave spectrum), to be the appropriate length scales to size the domain and damping layer thickness, and the optimal Rayleigh damping coefficient scales with the Brunt–Väisälä frequency.
Considering Froude numbers seen in wind farm applications, we propose recommendations to limit the reflections to less than 10 % of the total upwards propagating wave energy. Typically, damping is done at the top boundary, but given the non-periodic lateral boundary conditions of practical wind farm simulation domains, we find that damping the inflow-outflow boundaries is of equal importance to the top boundary. The Brunt–Väisälä frequency-normalized damping coefficient should be between 1 and 10. The damping layer thickness should be at least one effective vertical wavelength; damping layers exceeding 1.5 times the vertical wavelength are found to be unnecessary. The domain length and height should accommodate at least one effective horizontal and vertical wavelength, respectively. Moreover, Rayleigh damping does not damp the waves completely, and the non-damped energy might accumulate over the simulation time. ...
In this work, we only consider linearly stratified conditions (i.e., no inversion layer), thereby focusing on internal gravity waves in the free atmosphere and their reflections from the domain boundaries. This type of flow is governed by a single Froude number, which dictates most of the internal wave properties, such as wavelength, amplitude, and direction. This in turn will dictate the optimum domain size and Rayleigh damping layer setup. We find the effective horizontal and vertical wavelengths, (the representative wavelengths of the entire wave spectrum), to be the appropriate length scales to size the domain and damping layer thickness, and the optimal Rayleigh damping coefficient scales with the Brunt–Väisälä frequency.
Considering Froude numbers seen in wind farm applications, we propose recommendations to limit the reflections to less than 10 % of the total upwards propagating wave energy. Typically, damping is done at the top boundary, but given the non-periodic lateral boundary conditions of practical wind farm simulation domains, we find that damping the inflow-outflow boundaries is of equal importance to the top boundary. The Brunt–Väisälä frequency-normalized damping coefficient should be between 1 and 10. The damping layer thickness should be at least one effective vertical wavelength; damping layers exceeding 1.5 times the vertical wavelength are found to be unnecessary. The domain length and height should accommodate at least one effective horizontal and vertical wavelength, respectively. Moreover, Rayleigh damping does not damp the waves completely, and the non-damped energy might accumulate over the simulation time.
In recent years, the relevance of the interaction between neighboring wind farms has grown steadily. As one farm extracts energy from the wind, a downstream one can systematically experience lower wind speeds which threatens the economic viability of the farm. Significant progress has been made in understanding these farm-farm wake interactions, but we still lack methodologies to mitigate their undesired effects. In this study, we introduce Active Cluster Wake Mixing (ACWM). This novel method aims to accelerate the recovery of the cluster wake using dynamic control actions: By exciting the thrust of the individual turbines depending on their relative location, we generate non-uniform patterns of energy extraction. Phase offsets between the individual excitation signals propagate these regions through the wind farm. This results in large-scale velocity gradients inside the farm, which also affect the flow in the cluster wake region. An in-depth exploration and optimization of ACWM requires significant computational effort. Therefore, we compare three different wind farm modeling approaches in Large Eddy Simulations (LES) that differ in their computational costs regarding their suitability for further exploration of ACWM. For this purpose, we use an unoptimized ACWM scheme with two different excitation frequencies. For the first time ever we successfully show that ACWM manipulates the flow inside the wind farm with favorable effects on the wake velocity. We also demonstrate that the modeling of cluster wakes is challenging and has a significant effect on the potential gain.
The interaction of large wind farm clusters with the thermally stratified atmosphere has emerged as an important physical process that impacts the productivity of wind farms. Under stable conditions, this interaction triggers atmospheric gravity waves (AGWs) in the free atmosphere due to the vertical displacement of the atmospheric boundary layer (ABL) by the wind farm. AGWs induce horizontal pressure gradients within the ABL that alter the wind speed distribution within the farm, influencing both wind farm power generation and wake development. Additional factors, such as the growth of an internal boundary layer originating from the wind farm entrance and increased turbulence due to the wind turbines, further contribute to wake evolution. Recent studies have highlighted the considerable computational cost associated with simulating gravity wave effects within large-eddy simulations (LESs), as a significant portion of the free atmosphere must be resolved due to the large vertical spatial scales involved. Additionally, specialized boundary conditions are required to prevent wave reflections from contaminating the solution. In this study, we introduce a novel methodology to model the effects of AGWs without extending the LES computational domain into the free atmosphere. The proposed approach addresses the wave reflection problem inherently, eliminating the need for these specialized boundary conditions. We utilize the recently developed multi-scale coupled (MSC) model of Stipa et al. (2024b) to estimate the vertical ABL displacement triggered by the wind farm, and we apply the deformation to the domain of an LES that extends only to the inversion layer. The accuracy in predicting the AGW-induced pressure gradients is equivalent to the MSC model. The AGW modeling technique is verified for two distinct free-atmosphere stability conditions by comparing it to the traditional approach in which AGWs are fully resolved using a domain that extends several kilometers into the free atmosphere. The proposed approach accurately captures AGW effects on the row-averaged thrust and power distribution of wind farms while demanding 12.7 % of the computational resources needed for traditional methods. Moreover, when conventionally neutral boundary layers are studied, there is no longer a need to solve the potential temperature equation, as stability is neutral within the boundary layer. The developed approach is subsequently used to compare the global blockage and pressure disturbances obtained from the simulated cases against a solution characterized by zero boundary layer displacement, which represents the limiting case of very strong stratification above the boundary layer. This approximation, sometimes referred to as the “rigid lid”, is compared against the full AGW solution using the MSC model. This is done for different values of inversion strength and free atmosphere lapse rate, evaluating the ability of the “rigid lid” to predict blockage, wake effects, and overall wind farm performance.
As offshore wind farms grow in size, the blockage effect associated with the atmospheric gravity waves they trigger is expected to become more important. To model this, recent research has produced an Atmospheric Perturbation Model (APM), which simulates the mesoscale flow in the atmospheric boundary layer at a low computational cost compared to traditional methods. However, as a simplified reduced-order model, it can not resolve individual turbine wakes, and has to be coupled to an engineering wake model to predict farm power output. Over the years, three coupling methods have been developed, and been combined into the open-source framework WAYVE. This paper compares them, discussing both their theoretical validity and their performance. For the latter, we validate the velocities and power outputs predicted by WAYVE against 27 LES simulations. We find that the velocity matching (VM) and the pressure-based (PB) methods perform the best. Of these two, the VM method is more consistent with the APM output, while the PB method has a significantly lower computational cost.
Dynamic wind farm flow control is the art and science to maximize the energy yield of large wind farms. In this paper we will address the problem of large time delays between control actions of the different turbines in the farm and the delayed impact on the downstream turbines. We propose and show how a time-shifted cost function approach can render the receding horizon optimization problem more efficient and can mitigate the unavoidable turn-pike effect. We further show how the resulting setup can be used to break the optimization problem apart into several smaller optimization tasks to reduce the computational load. We demonstrate that the proposed changes do allow an economic model predictive control strategy to engage into collaborative wind farm control for long term gains, while a more traditional cost function approach leads to greedy turbine behavior. As a result, we take a crucial step towards a mature implementation of dynamic model based wind farm flow control.
The multi-scale coupled model
A new framework capturing wind farm–atmosphere interaction and global blockage effects
The growth in the number and size of wind energy projects in the last decade has revealed structural limitations in the current approach adopted by the wind industry to assess potential wind farm sites. These limitations are the result of neglecting the mutual interaction of large wind farms and the thermally stratified atmospheric boundary layer. While currently available analytical models are sufficiently accurate to conduct site assessments for isolated rotors or small wind turbine clusters, the wind farm's interaction with the atmosphere cannot be neglected for large-size arrays. Specifically, the wind farm displaces the boundary layer vertically, triggering atmospheric gravity waves that induce large-scale horizontal pressure gradients. These perturbations in pressure alter the velocity field at the turbine locations, ultimately affecting global wind farm power production. The implication of such dynamics can also produce an extended blockage region upstream of the first turbines and a favorable pressure gradient inside the wind farm. In this paper, we present the multi-scale coupled (MSC) model, a novel approach that allows the simultaneous prediction of micro-scale effects occurring at the wind turbine scale, such as individual wake interactions and rotor induction, and meso-scale phenomena occurring at the wind farm scale and larger, such as atmospheric gravity waves. This is achieved by evaluating wake models on a spatially heterogeneous background velocity field obtained from a reduced-order meso-scale model. Verification of the MSC model is performed against two large-eddy simulations (LESs) with similar average inflow velocity profiles and a different capping inversion strength, so that two distinct interfacial gravity wave regimes are produced, i.e. subcritical and supercritical. Interfacial waves can produce high blockage in the first case, as they are allowed to propagate upstream. On the other hand, in the supercritical regime their propagation speed is less than their advection velocity, and upstream blockage is only operated by internal waves. The MSC model not only proves to successfully capture both local induction and global blockage effects in the two analyzed regimes, but also captures the interaction between the wind farm and gravity waves, underestimating wind farm power by about only 2ĝ€¯% compared with the LES results. Conversely, wake models alone cannot distinguish between differences in thermal stratification, even if combined with a local induction model. Specifically, they are affected by a first-row overprediction bias that leads to an overestimation of the wind farm power by 13ĝ€¯% to 20ĝ€¯% for the modeled regimes.
A simple RANS closure for wind-farms under neutral atmospheric conditions
Preliminary findings
Accurately predicting wind turbine wake effects is essential for optimizing wind-farm performance and minimizing maintenance costs. This study explores the applicability of the Sparse Regression of Turbulent Stress Anisotropy (SpaRTA) framework to develop a simple yet robust Reynolds-averaged Navier-Stokes (RANS) model for wake prediction in wind energy contexts. The framework introduces two correction terms into two-equation models, with k - ϵ model being utilized in the current study. One correction term resembles the residual of the Turbulent Kinetic Energy (TKE) equation, and the other corrects the deviatoric part of the Reynolds Stress Tensor (RST). The terms are calculated from high-fidelity measurement or simulation data, and symbolic regression is used to determine the model for these terms. In this study, Large Eddy Simulation (LES) data from a single turbine is used as the training dataset, and a sample pre-selection process is employed to discover a correction model efficiently. The derived model incorporates two terms based on Pope's basis tensors and their invariants. The expression of the obtained model shows that it functions as a modification to the constant Cμ in the k - ϵ model. The model is evaluated by comparing its predicted velocity and TKE fields with the LES data used for the training. The model showed satisfactory performance in predicting both fields. Additionally, its generalizability is evaluated by testing it against a more complex six-turbine unseen case. The results indicate that the model effectively captures the velocity field and power output, but it tends to overpredict TKE, especially in the wake region.
Reproducing realistic date- and site-specific unsteady wind conditions in large-eddy simulations is becoming increasingly useful in wind energy. How to run a large-eddy simulation to match observed conditions, however, remains an open research question. One approach that has received considerable attention is mesoscale-to-microscale coupling, in which information about the mesoscale weather, most commonly acquired from a mesoscale numerical weather model, is passed on to a microscale model. In this paper, we demonstrate how the recently developed profile-assimilation technique, a form of mesoscale-to-microscale coupling, can be used to drive large-eddy simulations solely based on observed mean-flow profiles at a single location, bypassing the need for auxiliary mesoscale simulations. The new approach is evaluated for a diurnal cycle at the Scaled Wind Farm Technology site. Observed mean-flow profiles from the ground up to a height of 2 km are reconstructed by aggregating measurements from multiple instruments, and gaps in the data are infilled with natural neighbor interpolation. We perform nine simulations using various forcing approaches to deal with data limitations. The results show that it is indeed possible to drive microscale large-eddy simulation with observations using the profile-assimilation technique, notwithstanding large gaps in virtual potential temperature measurements. However, profile assimilation with vertical smoothing of the error between the desired and actual profiles is required. Without that smoothing, the microscale simulations develop unrealistically high turbulence levels under many situations. Finally, we show that simulated mesoscale data can account for missing observations, although care is needed as both data sources are not necessarily compatible.
Vertical temperature profiles influence the wind power generation of large offshore wind farms through stability-dependent effects such as blockage and gravity waves. However, numerical tools that are used to model these effects are often computationally too expensive to cover the large variety of atmospheric states occurring over time. Generally, an informed decision about which representative nonidealized situations to simulate is missing because of the lack of easily available information on representative vertical profiles, taking into account their spatiotemporal variability. Therefore, we present a novel framework that allows a smart selection of vertical temperature profiles. The framework consists of an improved analytical temperature model for the atmospheric boundary layer and lower troposphere, a subsequent clustering of these profiles to identify representatives, and last, a determination of areas with similar spatiotemporal characteristics of vertical profiles. When applying this framework on European ERA5 data, physically realistic representatives were identified for Europe, excluding the Mediterranean. Two or three profiles were found to be dominant for the open ocean, whereas more profiles prevail for land. Over the open ocean, weak temperature gradients in the boundary layer and a clear capping inversions are widespread, and stable profiles are absent except in the region of the East Icelandic Current. Interestingly, according to the ERA5 data, at its resolution, coastal areas and seas surrounded by land behave more similar to the land areas than to the open ocean, implying that a larger set of model integrations are needed for these areas to obtain representative results for offshore wind power assessments in comparison with the open ocean.
The Mesoscale to Microscale Coupling team, part of the U.S. Department of Energy Atmosphere to Electrons (A2e) initiative, has studied various important challenges related to coupling mesoscale models to microscale models for the use case of wind energy development and operation. Several coupling methods and techniques for generating turbulence at the microscale that is subgrid to the mesoscale have been evaluated for a variety of cases. Case studies included flat-terrain, complex-terrain, and offshore environments. Methods were developed to bridge the terra incognita, which scales from about 100ĝ€¯m through the depth of the boundary layer. The team used wind-relevant metrics and archived code, case information, and assessment tools and is making those widely available. Lessons learned and discerned best practices are described in the context of the cases studied for the purpose of enabling further deployment of wind energy.
This paper presents a new framework of the FLOw Redirection and Induction Dynamics (FLORIDyn) model. It is able to dynamically simulate the wake behaviour in wind farms under heterogeneous and changing environmental conditions at a low computational cost. The novelty of this work is the improved segregation of wake dynamics and wake influence: the framework creates Observation Points (OPs) at each turbine, which propagate wind field states and turbine states downstream and follow the wind direction of the free stream velocity. These observation points cover the dynamic aspects of the simulation. The OPs, along with the stored states, are now used to derive so-called Temporary Wind Farms (TWF), which approximate the effective intra-farm wind conditions at a given location. Within these TWF, the flow conditions are homogeneous and steady state. This way, arbitrary wake models can be used to calculate the farm influence on the location. The FLORIDyn framework also provides interfaces to flow field estimators, which is tested with an effective wind speed estimator. A nine turbine case is used to highlight the quality and performance of the simulation result. Compared to its predecessor, the new FLORIDyn framework decreases the computational cost by one to two orders of magnitude, which makes it a promising candidate for real-time model predictive dynamic wind farm control.
The revised FLORIDyn model
Implementation of heterogeneous flow and the Gaussian wake
Both FLORIS and FLORIDyn are parametric models which can be used to simulate wind farms, evaluate controller performance and can serve as a control-oriented model. One central element in which they differ is in their representation of flow dynamics: FLORIS neglects these and provides a computationally very cheap approximation of the mean wind farm flow. FLORIDyn defines a framework which utilizes this low computational cost of FLORIS to simulate basic wake dynamics. This is achieved by creating so-called observation points (OPs) at each time step at the rotor plane which inherit the turbine state.
In this work, we develop the initial FLORIDyn framework further considering multiple aspects. The underlying FLORIS wake model is replaced by a Gaussian wake model. The distribution and characteristics of the OPs are adapted to account for the new parametric model but also to take complex flow conditions into account. To achieve this, a mathematical approach is developed to combine the parametric model and the changing, heterogeneous world conditions and link them with each OP. We also present a computationally lightweight wind field model to allow for a simulation environment in which heterogeneous flow conditions are possible.
FLORIDyn is compared to Simulator for Offshore Wind Farm Applications (SOWFA) simulations in three- and nine-turbine cases under static and changing environmental conditions. The results show a good agreement with the timing of the impact of upstream state changes on downstream turbines. They also show a good agreement in terms of how wakes are displaced by wind direction changes and when the resulting velocity deficit is experienced by downstream turbines. A good fit of the mean generated power is ensured by the underlying FLORIS model. In the three-turbine case, FLORIDyn simulates 4 s simulation time in 24.49 ms computational time. The resulting new FLORIDyn model proves to be a computationally attractive and capable tool for model-based dynamic wind farm control. ...
Both FLORIS and FLORIDyn are parametric models which can be used to simulate wind farms, evaluate controller performance and can serve as a control-oriented model. One central element in which they differ is in their representation of flow dynamics: FLORIS neglects these and provides a computationally very cheap approximation of the mean wind farm flow. FLORIDyn defines a framework which utilizes this low computational cost of FLORIS to simulate basic wake dynamics. This is achieved by creating so-called observation points (OPs) at each time step at the rotor plane which inherit the turbine state.
In this work, we develop the initial FLORIDyn framework further considering multiple aspects. The underlying FLORIS wake model is replaced by a Gaussian wake model. The distribution and characteristics of the OPs are adapted to account for the new parametric model but also to take complex flow conditions into account. To achieve this, a mathematical approach is developed to combine the parametric model and the changing, heterogeneous world conditions and link them with each OP. We also present a computationally lightweight wind field model to allow for a simulation environment in which heterogeneous flow conditions are possible.
FLORIDyn is compared to Simulator for Offshore Wind Farm Applications (SOWFA) simulations in three- and nine-turbine cases under static and changing environmental conditions. The results show a good agreement with the timing of the impact of upstream state changes on downstream turbines. They also show a good agreement in terms of how wakes are displaced by wind direction changes and when the resulting velocity deficit is experienced by downstream turbines. A good fit of the mean generated power is ensured by the underlying FLORIS model. In the three-turbine case, FLORIDyn simulates 4 s simulation time in 24.49 ms computational time. The resulting new FLORIDyn model proves to be a computationally attractive and capable tool for model-based dynamic wind farm control.