Stefan Ivanell

Numerical Computations of Wind Turbine Wakes

Numerical simulations of the Navier-Stokes equations are performed to achieve a better understanding of the behaviour of wakes generated by wind turbines. The simulations are performed by combining ...

Simulation of wind turbine wakes using the actuator line technique

The actuator line technique was introduced as a numerical tool to be employed in combination with large eddy simulations to enable the study of wakes and wake interaction in wind farms. The technique is today largely used for studying basic features of wakes as well as for making performance predictions of wind farms. In this paper, we give a short introduction to the wake problem and the actuator line methodology and present a study in which the technique is employed to determine the near-wake properties of wind turbines. The presented results include a comparison of experimental results of the wake characteristics of the flow around a three-bladed model wind turbine, the development of a simple analytical formula for determining the near-wake length behind a wind turbine and a detailed investigation of wake structures based on proper orthogonal decomposition analysis of numerically generated snapshots of the wake.

Wake downstream of the Lillgrund wind farm - A Comparison between LES using the actuator disc method and a Wind farm Parametrization in WRF

Simulations of the Lillgrund wind farm (located between Malmo and Copenhagen) are performed using both Large Eddy Simulation (LES) and mesoscale simulations in WRF. The aim is to obtain a better un ...

Complex terrain experiments in the New European Wind Atlas

The New European Wind Atlas project will create a freely accessible wind atlas covering Europe and Turkey, develop the model chain to create the atlas and perform a series of experiments on flow in many different kinds of complex terrain to validate the models. This paper describes the experiments of which some are nearly completed while others are in the planning stage. All experiments focus on the flow properties that are relevant for wind turbines, so the main focus is the mean flow and the turbulence at heights between 40 and 300 m. Also extreme winds, wind shear and veer, and diurnal and seasonal variations of the wind are of interest. Common to all the experiments is the use of Doppler lidar systems to supplement and in some cases replace completely meteorological towers. Many of the lidars will be equipped with scan heads that will allow for arbitrary scan patterns by several synchronized systems. Two pilot experiments, one in Portugal and one in Germany, show the value of using multiple synchronized, scanning lidar, both in terms of the accuracy of the measurements and the atmospheric physical processes that can be studied. The experimental data will be used for validation of atmospheric flow models and will by the end of the project be freely available. This article is part of the themed issue ‘Wind energy in complex terrains’.

Comparison of Engineering Wake Models with CFD Simulations

The engineering wake models by Jensen [1] and Frandsen et al. [2] are assessed for different scenarios simulated using Large Eddy Simulation and the Actuator Line method implemented in the Navier-Stokes equations. The scenarios include the far wake behind a single wind turbine, a long row of turbines in an atmospheric boundary layer, idealised cases of an infinitely long row of wind turbines and infinite wind farms with three different spacings. Both models include a wake expansion factor, which is calibrated to fit the simulated wake velocities. The analysis highlights physical deficiencies in the ability of the models to universally predict the wake velocities, as the expansion factor can be fitted for a given case, but with not apparent transition between the cases.

Numerical analysis of the tip and root vortex position in the wake of a wind turbine

The stability of tip and root vortices are studied numerically in order to analyse the basic mechanism behind the break down of tip and root vortices. The simulations axe performed using the CFD pr ...

IEA wind annex XX: HAWT aerodynamic and models from wind tunnel measurements - final report

This work characterizes undocumented physical relationships that govern aerodynamic force time variations that take place in connection with rotational augmentation on rotating wind turbine blades.

Quantifying variability of Large Eddy Simulations of very large wind farms

Large Eddy Simulations are inherently dynamic as the largest scales are resolved and the smallest scales are modeled temporally. This raises challenges for simulations including very large scales such as atmospheric flows, which require very long simulation times. Simple averages fail at capturing these dynamics and potentially yield misleading interpretations concerning the capabilities of different models when tested in blind tests or in benchmarking exercises such as Wakebench, where results from different flow models are compared. This article will present results from very large wind farm simulations using Actuator Disc (AD) and Line (AL) models for two different turbine spacings with turbulent inflow. The results of each numerical flow model include a certain variability, and it will be examined if different models result in comparable probability distributions.

Numerical investigation of the wake interaction between two model wind turbines with span-wise offset

Wake interaction between two model scale wind turbines with span-wise offset is investigated numerically using Large Eddy Simulation (LES) and the results are validated against the experimental data. An actuator line technique is used for modeling the rotor. The investigated setup refers to a series of experimental measurements of two model scale turbines conducted by NTNU in low speed wind tunnel in which the two wind turbines are aligned with a span-wise offset resulting in half wake interaction. Two levels of free-stream turbulence are tested, the minimum undisturbed level of about Ti ≈ 0.23% and a high level of about Ti ≈ 10% using a passive upstream grid. The results show that the rotor characteristics for both rotors are well captured numerically even if the downstream rotor operates into stall regimes. There are however some difficulties in correct prediction of the thrust level. The interacting wake development is captured in great details in terms of wake deficit and streamwise turbulence kinetic energy. The present work is done in connection with Blind test 3 workshops organized jointly by NOWITECH and NORCOWE.

Determination of Wind Turbine Near-Wake Length Based on Stability Analysis

A numerical study on the wake behind a wind turbine is carried out focusing on determining the length of the near-wake based on the instability onset of the trailing tip vortices shed from the turbine blades. The numerical model is based on large-eddy simulations (LES) of the Navier-Stokes equations using the actuator line (ACL) method. The wake is perturbed by applying stochastic or harmonic excitations in the neighborhood of the tips of the blades. The flow field is then analyzed to obtain the stability properties of the tip vortices in the wake of the wind turbine. As a main outcome of the study it is found that the amplification of specific waves (traveling structures) along the tip vortex spirals is responsible for triggering the instability leading to wake breakdown. The presence of unstable modes in the wake is related to the mutual inductance (vortex pairing) instability where there is an out-of-phase displacement of successive helix turns. Furthermore, using the non-dimensional growth rate, it is found that the pairing instability has a universal growth rate equal to π/2. Using this relationship, and the assumption that breakdown to turbulence occurs once a vortex has experienced sufficient growth, we provide an analytical relationship between the turbulence intensity and the stable wake length. The analysis leads to a simple expression for determining the length of the near wake. This expression shows that the near wake length is inversely proportional to thrust, tip speed ratio and the logarithmic of the turbulence intensity.