Richard Stevens

Physics of Fluids

University of Twente

Research Highlights

Selected papers illustrating wind-farm turbulence, atmospheric-boundary-layer coupling, high-performance simulation, and canonical turbulence. For the complete chronological record, see the publications page.

Large-eddy simulation of turbine wakes in an extended wind farm

Our research uses large-eddy simulation, direct numerical simulation, physical modeling, and high-performance computing to study turbulent flows in wind energy, atmospheric boundary layers, and canonical fluid systems. The highlights below show how we connect fundamental turbulence physics to predictive models for wind-farm performance, atmospheric exchange, flow variability, noise, and heat transport.

Large-eddy simulation of turbulent wakes in an extended wind farm. Such simulations reveal how turbine-scale wakes, farm-scale blockage, and atmospheric-boundary-layer dynamics interact.

Research map

Featured highlights

Recent highlights

Impact of atmospheric turbulence on performance and loads of wind turbines: knowledge gaps and research challenges

B. Kosović, S. Basu, J. Berg, L.K. Berg, S.E. Haupt, X.G. Larsén, J. Peinke, R.J.A.M. Stevens, P. Veers & S. Watson, Wind Energy Science, 2026

This community review identifies how atmospheric-boundary-layer turbulence affects wind-turbine power production, loads, design, and wind-farm operation, and outlines key knowledge gaps.

Mean turbulent momentum fluxes and wind deficits in nocturnal stable atmospheric boundary layers

Z. Shen, L. Liu, X.-Y. Lu & R.J.A.M. Stevens, Journal of Fluid Mechanics, 2025

This paper connects stable stratification, turbulent momentum fluxes, and wind deficits in nocturnal boundary layers, with direct relevance for offshore wind-farm inflow and wake behavior.

The global properties of nocturnal stable atmospheric boundary layers

Z. Shen, L. Liu, X. Lu & R.J.A.M. Stevens, Journal of Fluid Mechanics, 2024

This study analyzes the global structure of nocturnal stable atmospheric boundary layers and clarifies how stratification modifies turbulent transport and boundary-layer depth.

Simulation and modeling of wind farms in baroclinic atmospheric boundary layers

J.H. Kasper, A. Stieren & R.J.A.M. Stevens, Journal of Renewable and Sustainable Energy, 2024

This work extends wind-farm simulation and modeling to baroclinic boundary layers, where wind direction and shear vary with height and alter wake recovery and power production.

Modeling wind farm noise emission and propagation: Effects of flow and layout

J. Colas, A. Emmanuelli, D. Dragna & R.J.A.M. Stevens, Renewable Energy, 2026

This paper models how wind-farm layout and flow conditions jointly determine turbine noise emission and propagation.

Low-frequency wind speed variations and their impact on wind farm performance

Y. Liu & R.J.A.M. Stevens, Journal of Renewable and Sustainable Energy, in press

This study examines how slow atmospheric wind-speed variations affect wind-farm power production and farm-scale performance variability. Forthcoming - a DOI will be added once formally published; see the publications page for the current record.

Wind-farm turbulence and wake physics

Large wind farms operate in a multiscale turbulent flow. Turbine wakes interact with neighboring turbines, merge into farm-scale wakes, and alter the exchange of momentum between the surface layer and the atmosphere above. We use large-eddy simulation and reduced-order modeling to understand these processes and improve predictive wind-farm models. See also the wind-farm LES and analytical wind-farm modeling pages.

Effects of turbine spacing on the power output of extended wind-farms

R.J.A.M. Stevens, D.F. Gayme & C. Meneveau, Wind Energy, 2016

This study uses large-eddy simulation to quantify how streamwise and spanwise turbine spacing influence the power output of very large wind farms. It links layout effects to the vertical flux of kinetic energy into the turbine region.

Key idea: Wind-farm layout controls not only local wake losses but also farm-scale replenishment of kinetic energy from above.

Comparison of large eddy simulations using actuator disk or actuator line models with wind tunnel experiments

R.J.A.M. Stevens, L.A. Martínez Tossas & C. Meneveau, Renewable Energy, 2018

This paper compares actuator-disk and actuator-line turbine models in LES against wind-tunnel measurements. It clarifies when simplified turbine representations are sufficient and when more detailed actuator-line modeling is needed.

Key idea: High-fidelity wind-farm simulations must be validated against experiments to ensure that wake physics and turbine interactions are captured correctly.

Atmospheric coupling and wind-energy limits

Wind-farm performance is controlled by the atmospheric boundary layer. Stability, low-level jets, baroclinicity, geostrophic forcing, and turbulent momentum transport determine how much kinetic energy is available to the farm and how quickly wakes recover. See also the turbulent boundary layer page.

  • Understanding wind farm power densities - also featured above.
  • The global properties of nocturnal stable atmospheric boundary layers, and Mean turbulent momentum fluxes and wind deficits in nocturnal stable atmospheric boundary layers - also listed under Recent highlights above.

Universal Wind Profile for Conventionally Neutral Atmospheric Boundary Layers

L. Liu, S.N. Gadde & R.J.A.M. Stevens, Physical Review Letters, 2021

This study develops and validates a universal velocity profile for conventionally neutral atmospheric boundary layers. It provides a theoretical basis for wind profiles relevant to wind-energy applications and atmospheric modeling.

Key idea: Accurate wind-farm modeling requires physically consistent descriptions of the atmospheric boundary layer, not only turbine-wake models.

Impact of Negative Geostrophic Wind Shear on Wind Farm Performance

A. Stieren, J.H. Kasper, S.N. Gadde & R.J.A.M. Stevens, PRX Energy, 2022

This paper shows how the direction of geostrophic wind shear in the atmosphere changes wind-farm power production, highlighting the importance of atmospheric forcing beyond idealized neutral or stable boundary layers.

Key idea: Large wind farms are sensitive to the vertical structure of geostrophic forcing, not only to the near-surface wind speed.

Effect of low-level jet height on wind farm performance

S.N. Gadde & R.J.A.M. Stevens, Journal of Renewable and Sustainable Energy, 2021

This work shows how the height of a nocturnal low-level jet relative to the turbine rotor changes wake recovery and wind-farm performance.

Key idea: Atmospheric structure can directly determine whether wake recovery improves or degrades.

Multiscale wind-farm variability

Wind-farm power varies over seconds, minutes, hours, and longer atmospheric time scales. We study how coherent atmospheric motions, turbulent structures, and turbine-array interactions shape aggregate power fluctuations.

Temporal structure of aggregate power fluctuations in large-eddy simulations of extended wind-farms

R.J.A.M. Stevens & C. Meneveau, Journal of Renewable and Sustainable Energy, 2014

This paper characterizes the temporal spectrum of aggregate power output fluctuations from large-eddy simulations of extended wind farms.

Key idea: Aggregate wind-farm power fluctuations have a distinct spectral signature that reflects farm-scale turbulence.

A wavenumber-frequency spectral model for atmospheric boundary layers

M. Wilczek, R.J.A.M. Stevens, Y. Narita & C. Meneveau, Journal of Physics: Conference Series, 2014

This paper develops a wavenumber-frequency spectral model for turbulence in atmospheric boundary layers, relevant to wind-power forecasting.

Key idea: A compact spectral model can connect atmospheric turbulence statistics to wind-power fluctuation forecasting.

Wind-farm noise and environmental impact

Wind-farm flow physics affects more than power production. Wakes, turbine layout, atmospheric turbulence, and rotor operating conditions influence noise emission, propagation, and amplitude modulation.

  • Modeling wind farm noise emission and propagation: Effects of flow and layout - also listed under Recent highlights above.

Wake-induced variations in noise levels and amplitude modulation for two interacting wind turbines

J. Colas, A. Emmanuelli, D. Dragna & R.J.A.M. Stevens, Journal of the Acoustical Society of America, 2026

This study examines how the wake of an upstream turbine changes the noise level and amplitude modulation of a downstream turbine. It provides a mechanistic link between wake interaction and perceived noise variability.

Key idea: Turbine wakes can modify downstream noise characteristics, creating flow-dependent acoustic variability.

Three-dimensional effects of the wake on wind turbine sound propagation using parabolic equation

H. Bommidala, J. Colas, A. Emmanuelli, D. Dragna, C. Khodr, B. Cotté & R.J.A.M. Stevens, Journal of Sound and Vibration, 2025

This study models the three-dimensional effects of a turbine wake on downstream sound propagation using a parabolic-equation method.

Key idea: Wake-induced flow structures shape how turbine noise propagates in three dimensions, not just along the ground.

Impact of a Two-Dimensional Steep Hill on Wind Turbine Noise Propagation

J. Colas, A. Emmanuelli, D. Dragna, P. Blanc-Benon, B. Cotté & R.J.A.M. Stevens, Wind Energy Science, 2024

This paper quantifies how terrain, specifically a steep hill, affects atmospheric propagation of wind-turbine noise.

Key idea: Local terrain features can meaningfully change how far and in what direction wind-turbine noise propagates.

High-performance simulation and open-source tools

High-fidelity turbulence simulations require scalable numerical methods and efficient use of modern supercomputers. We develop and use simulation tools for wind-farm LES and canonical DNS, including the open-source AFiD framework.

AFiD direct numerical simulation of Rayleigh-Bénard convection

AFiD-GPU: a versatile Navier-Stokes Solver for Wall-Bounded Turbulent Flows on GPU Clusters

X. Zhu et al., Computer Physics Communications, 2018

AFiD-GPU is a high-performance implementation of the AFiD incompressible Navier-Stokes solver for GPU clusters. It enables large-scale simulations of canonical turbulent flows such as Rayleigh-Bénard convection, Taylor-Couette flow, channel flow, and plane Couette flow.

Key idea: Open, scalable simulation tools are essential for connecting fundamental turbulence physics with high-resolution numerical experiments.

Comparison of computational codes for direct numerical simulations of turbulent Rayleigh-Bénard convection

G.L. Kooij, M.A. Botchev, E.M.A. Frederix, B.J. Geurts, S. Horn, D. Lohse, E.P. van der Poel, O. Shishkina, R.J.A.M. Stevens & R. Verzicco, Computers & Fluids, 2018

This paper cross-validates independent direct numerical simulation codes for wall-bounded thermal turbulence, including the AFiD code.

Key idea: Cross-code validation is essential for trusting high-fidelity DNS results across research groups.

Canonical turbulence and thermal convection

Canonical turbulent flows provide controlled systems for studying transport, coherent structures, and scaling behavior. These studies support the physical understanding and numerical methods used across our work on wind-energy and environmental flows. See also the thermal convection page.

Temperature field showing turbulent thermal superstructures in Rayleigh-Bénard convection

Turbulent thermal superstructures in Rayleigh-Bénard convection

R.J.A.M. Stevens, A. Blass, X. Zhu, R. Verzicco & D. Lohse, Physical Review Fluids, 2018

This study identifies large-scale, long-lived thermal superstructures in turbulent Rayleigh-Bénard convection. It shows that sufficiently large domains are required to capture the organization and transport role of these coherent structures.

Key idea: Turbulent convection contains persistent large-scale organization that can strongly influence heat transport.

How wide must Rayleigh–Bénard cells be to prevent finite aspect ratio effects in turbulent flow?

R.J.A.M. Stevens, R. Hartmann, R. Verzicco & D. Lohse, Journal of Fluid Mechanics, 2024

This paper examines how domain width affects turbulent Rayleigh-Bénard convection. It clarifies when finite-aspect-ratio effects contaminate heat transport and flow organization.

Key idea: Reliable simulations of turbulent convection require domains large enough to contain the relevant large-scale structures.

Direct numerical simulation of sheared thermal convection flow field, winner of the SURFsara 2017 Visualization Competition

Scaling relations for heat and momentum transport in sheared Rayleigh-Bénard convection

G.S. Yerragolam, C.J. Howland, R.J.A.M. Stevens, R. Verzicco, O. Shishkina & D. Lohse, Journal of Fluid Mechanics, 2024

This study derives how an imposed mean shear changes the classical scaling of heat and momentum transport in convection.

Key idea: Shear changes the balance between heat transport, momentum transport, and coherent turbulent structures.

Related work also covers rotating Rayleigh-Bénard convection and Taylor-Couette turbulence; see those research pages and the publications list for the corresponding papers.

Public-facing results

Some results provide especially clear examples of why wind-farm flow physics matters beyond specialist turbulence research.

  • Effect of low-level jet height on wind farm performance, and Modeling wind farm noise emission and propagation - also listed above.

Enhanced wind-farm performance using windbreaks

L. Liu & R.J.A.M. Stevens, Physical Review Fluids, 2021

Large-eddy simulations showed that carefully designed low windbreaks can increase power production in large wind farms. The important point is not simply that barriers speed up the flow above them, but that the optimal windbreak for a wind farm is much lower than the optimal windbreak for a single isolated turbine.

Key idea: Small changes to the flow entering a wind farm can have farm-scale consequences for power production.

For the complete chronological list of peer-reviewed articles, see the publications page.