scholarly journals Stochastic Wake Modeling Based on POD Analysis

Author(s):  
David Bastine ◽  
Lukas Vollmer ◽  
Matthias Wächter ◽  
Joachim Peinke

Abstract. In this work, large eddy simulation data is analyzed to investigate a new stochastic modeling approach for the wake of a wind turbine. The data is generated by the LES model PALM combined with an actuator disk with rotation representing the turbine. After applying a proper orthogonal decomposition (POD), three different stochastic models for the weighting coefficients of the POD modes are deduced resulting in three wake models. Their performance is investigated mainly on the basis of aeroelastic simulations of a wind turbine in the wake. Three different load cases and their statistic characteristics are compared for the original LES, truncated PODs and the stochastic wake models including different numbers of POD modes. It is shown that approximately six POD modes are enough to capture the load dynamics on large temporal scales. Modeling the weighting coefficients as independent stochastic processes leads to similar load characteristics as in the case of the truncated POD. To complete this simplified wake description, we show evidence that the small-scale dynamics can be grasped by adding to our model a homogeneous turbulent field. In this way, we present a procedure how to derive stochastic wake models from costly CFD calculations or elaborated experimental investigations. These numerically efficient models provide the added value of possible long-term studies. Depending on the aspects of interest, different minimalized models may be obtained.

2014 ◽  
Vol 71 (11) ◽  
pp. 4001-4027 ◽  
Author(s):  
Peter P. Sullivan ◽  
James C. McWilliams ◽  
Edward G. Patton

Abstract Momentum and scalar transport in the marine atmospheric boundary layer (MABL) is driven by a turbulent mix of winds, buoyancy, and surface gravity waves. To investigate the interaction between these processes, a large-eddy simulation (LES) model is developed with the capability to impose a broadband spectrum of time-varying finite-amplitude surface waves at its lower boundary. The LES model adopts a Boussinesq flow model and integrates the governing equations on a time-varying, surface-fitted, nonorthogonal mesh using cell-centered variables with special attention paid to the solution of the pressure Poisson equation near the wavy boundary. Weakly unstable MABLs are simulated with geostrophic winds increasing from 5 to 25 m s−1 and wave age varying from swell-dominated to wind-wave equilibrium. The simulations illustrate cross-scale coupling as wave-impacted near-surface turbulence transitions into shear-convective rolls with increasing distance from the water. In a regime with swell, low winds, and weak heating, wave-induced vertical velocity and pressure signals are readily observed well above the standard reference height ζa = 10 m. At wind-wave equilibrium, the small-scale wave-induced signals are detectable only near the water surface. Below ζa, a nearly-constant-flux layer is observed where the momentum flux carried by turbulence, form stress, and subgrid-scale motions shifts with varying wave age and distance above the water. The spectral content of the surface form stress is wave-age dependent, especially at low wavenumbers. The LES wind profiles deviate from Monin–Obukhov similarity theory in nonequilibrium wind-wave conditions, and entrainment is greatly enhanced by shear-induced engulfment events.


Author(s):  
Charlie Koupper ◽  
Tommaso Bacci ◽  
Bruno Facchini ◽  
Alessio Picchi ◽  
Lorenzo Tarchi ◽  
...  

To deepen the knowledge of the interaction between modern lean burn combustors and high pressure turbines, a non-reactive real scale annular trisector Combustor Simulator (CS) has been assembled at University of Florence, with the goal of investigating and characterizing the combustor aerothermal field as well as the hot streak transport towards the high pressure vanes. To generate hot streaks and simulate lean burn combustor behaviors, the rig is equipped with axial swirlers fed by a main air flow stream that is heated up to 531 K, while liners with effusion cooling holes are fed by air at ambient temperature. Detailed experimental investigations are then performed with the aim of characterizing the turbulence quantities at the exit of the combustion module, and specifically evaluating an integral scale of turbulence. To do so, an automatic traverse system is mounted at the exit of the CS and equipped to perform Hot Wire Anemometry (HWA) measurements. In this paper, two-point correlations are computed from the time signal of the axial velocity giving access to an evaluation of the turbulence timescales at each measurement point. For assessment of the advanced numerical method that is Large Eddy Simulation (LES), the same methodology is applied to a LES prediction of the CS. Although comparisons seem relevant and easily accessible, both approaches and contexts have fundamental differences: mostly in terms of duration of the signals acquired experimentally and numerically but also with potentially different acquisition frequencies. In the exercise that aims at comparing high-order statistics and diagnostics, the specificity of comparing experimental and numerical results is comprehensively discussed. Attention is given to the importance of the acquisition frequency, intrinsic bias of having a short duration signal and influence of the investigating windows. For an adequate evaluation of the turbulent time scales, it is found that comparing experiments and numerics for high Reynolds number flows inferring small-scale phenomena requires to obey a set of rules, otherwise important errors can be made. If adequately processed, LES and HWA are found to agree well indicating the potential of LES for such problems.


2012 ◽  
Vol 13 (5) ◽  
pp. 1421-1441 ◽  
Author(s):  
Charles Talbot ◽  
Elie Bou-Zeid ◽  
Jim Smith

Abstract This paper assesses the performance of the Weather Research and Forecasting Model (WRF) as a tool for multiscale atmospheric simulations. Tests are performed in real and idealized cases with multiple configurations and with resolutions ranging from the mesoscale (gridcell size ~10 km) for the real cases to local scales (gridcell size ~50 m) for both real and idealized cases. All idealized simulations and the finest real-case simulations use the turbulence-resolving large-eddy simulation mode of WRF (WRF-LES). Tests in neutral conditions and with idealized forcing are first performed to assess the model’s sensitivity to grid resolutions and subgrid-scale parameterizations and to optimize the setup of the real cases. An increase in horizontal model resolution is found to be more beneficial than an increase in vertical resolution. WRF-LES is then tested, using extensive observational data, in real-world cases over complex terrain through nested simulations in which the mesoscale domains drive the LES domains. Analysis of the mesoscale simulations indicates that the data needed to force the largest simulated domain and to initialize surface parameters have the strongest influence on the results. Similarly, LES model fields are primarily influenced by their mesoscale meteorological forcing. As a result, the nesting of LES models down to a 50-m resolution does not improve all aspects of hydrometeorological predictions. Advantages of using fine-resolution LES are noted at nighttime (under stable conditions) and over heterogeneous surfaces when local properties are required or when resolving small-scale surface features is desirable.


2020 ◽  
Vol 148 (3) ◽  
pp. 1029-1048 ◽  
Author(s):  
Mikhail Varentsov ◽  
Igor Esau ◽  
Tobias Wolf

Abstract Detailed temperature maps are required in various applications. Any temperature interpolation over complex terrain must account for differences in land cover and elevation. Local circulations and other small-scale factors can also perturb the temperature. This study considers the surface air temperature T mapping with geostatistical kriging. The kriging methods are implemented for both T and temperature anomalies ΔT, defined as difference between T at a given location and T at the same elevation in the free atmosphere. The study explores the parallelized atmospheric large-eddy simulation (LES) model (PALM) as a source for variogram and external drift in the kriging methods. Ten kriging methods for the temperature mapping have been considered: ordinary kriging (OK) of T and ΔT with variogram derived from the observations (methods 1 and 2, correspondingly); OK of T and ΔT with variogram derived from LES data (3 and 4); universal kriging with external drift (KED) that utilizes the LES data (5 and 6); a weighted combination of KED of T and ΔT (method 7); and methods 5, 6, and 7 enhanced with additional “virtual” points in remote areas (methods 8, 9, and 10). These 10 methods are evaluated for eight typical weather situations observed in Bergen, Norway. Our results demonstrate considerable added value of the LES information for the detailed meteorological temperature mapping. The LES data improve both the variogram and the resulting temperature maps, especially in the remote mountain parts of the domain and along the coast.


2011 ◽  
Vol 2011 ◽  
pp. 1-7 ◽  
Author(s):  
Takanori Uchida ◽  
Yuji Ohya ◽  
Kenichiro Sugitani

The wake of a wind turbine generator (WTG) operated at the optimal tip speed ratio is compared to the wake of a WTG with its rotor replaced by a stationary disk. Numerical simulations are conducted with a large eddy simulation (LES) model using a nonuniform staggered Cartesian grid. The results from the numerical simulations are compared to those from wind-tunnel experiments. The characteristics of the wake of the stationary disk are significantly different from those of the WTG. The velocity deficit at a downstream distance of (: rotor diameter) behind the WTG is approximately 30 to 40% of the inflow velocity. In contrast, flow separation is observed immediately behind the stationary disk (), and the velocity deficit in the far wake () of the stationary disk is smaller than that of the WTG.


Energies ◽  
2019 ◽  
Vol 12 (7) ◽  
pp. 1306
Author(s):  
Camilo Sedano ◽  
Frederik Berger ◽  
Hamid Rahimi ◽  
Omar Lopez Mejia ◽  
Martin Kühn ◽  
...  

With the ongoing increase in the size of wind turbines, experimental investigations have become more complicated and expensive. Therefore, computational models have proven to be a viable solution for design purposes. This article aims to validate CFD simulations of an experimental model wind turbine (MoWiTO 1.8) using Delayed Detached Eddy Simulation (DDES) and Improved DDES (IDDES) turbulence modelling approaches. For the purpose of validation, integral quantities (such as power, thrust, torque and blade-root bending moment in the flapwise direction) measured in the wind tunnel are compared with numerical results obtained with OpenFOAM. In general, the computational results show a very good agreement with the measurements for most of the monitored quantities. In particular, the blade-root bending moment presents the largest difference, taking into account that the simulation assumes the turbine blades are rigid. Nevertheless, the simulation does achieve in recreating the turbulent behavior as can be evidenced by the Power Spectral Density graphs, and the wake’s velocity measurements. In general, the IDDES turbulent model achieves a better agreement to the experimental results, while maintaining a very similar computational time as the DDES model.


2007 ◽  
Vol 64 (12) ◽  
pp. 4445-4457 ◽  
Author(s):  
M. Antonelli ◽  
R. Rotunno

Abstract This paper describes results from a large-eddy simulation (LES) model used in an idealized setting to simulate the onset of the sea breeze. As the LES is capable of simulating boundary layer–scale, three-dimensional turbulence along with the mesoscale sea-breeze circulation, a parameterization of the planetary boundary layer was unnecessary. The basic experimental design considers a rotating, uniformly stratified, resting atmosphere that is suddenly heated at the surface over the “land” half of the domain. To focus on the simplest nontrivial problem, the diurnal cycle, effects of moisture, interactions with large-scale winds, and coastline curvature were all neglected in this study. The assumption of a straight coastline allows the use of a rectangular computational domain that extends to 50 km on either side of the coast, but only 5 km along the coast, with 100-m grid intervals so that the small-scale turbulent convective eddies together with the mesoscale sea breeze may be accurately computed. Through dimensional analysis of the simulation results, the length and velocity scales characterizing the simulated sea breeze as functions of the externally specified parameters are identified.


Author(s):  
Charlie Koupper ◽  
Laurent Gicquel ◽  
Florent Duchaine ◽  
Tommaso Bacci ◽  
Bruno Facchini ◽  
...  

To deepen the knowledge of the interaction between modern lean burn combustors and high pressure (HP) turbines, a nonreactive real scale annular trisector combustor simulator (CS) has been assembled at University of Florence (UNIFI), with the goal of investigating and characterizing the combustor aerothermal field as well as the hot streak transport toward the HP vanes. To generate hot streaks and simulate lean burn combustor behaviors, the rig is equipped with axial swirlers fed by a main air flow stream that is heated up to 531 K, while liners with effusion cooling holes are fed by air at ambient temperature. Detailed experimental investigations are then performed with the aim of characterizing the turbulence quantities at the exit of the combustion module, and specifically evaluating an integral scale of turbulence. To do so, an automatic traverse system is mounted at the exit of the CS and equipped to perform hot wire anemometry (HWA) measurements. In this paper, two-point correlations are computed from the time signal of the axial velocity giving access to an evaluation of the turbulence timescales at each measurement point. For assessment of the advanced numerical method that is large Eddy simulation (LES), the same methodology is applied to a LES prediction of the CS. Although comparisons seem relevant and easily accessible, both approaches and contexts have fundamental differences: mostly in terms of duration of the signals acquired experimentally and numerically but also with potentially different acquisition frequencies. In the exercise that aims at comparing high-order statistics and diagnostics, the specificity of comparing experimental and numerical results is comprehensively discussed. Attention is given to the importance of the acquisition frequency, intrinsic bias of having a short duration signal and influence of the investigating windows. For an adequate evaluation of the turbulent time scales, it is found that comparing experiments and numerics for high Reynolds number flows inferring small-scale phenomena requires to obey a set of rules, otherwise important errors can be made. If adequately processed, LES and HWA are found to agree well indicating the potential of LES for such problems.


Energies ◽  
2021 ◽  
Vol 14 (13) ◽  
pp. 3745
Author(s):  
Tristan Revaz ◽  
Fernando Porté-Agel

Large-eddy simulation (LES) with actuator models has become the state-of-the-art numerical tool to study the complex interaction between the atmospheric boundary layer (ABL) and wind turbines. In this paper, a new evaluation of actuator disk models (ADMs) for LES of wind turbine flows is presented. Several details of the implementation of such models are evaluated based on a test case studied experimentally. In contrast to other test cases used in previous similar studies, the present test case consists of a wind turbine immersed in a realistic turbulent boundary-layer flow, for which accurate data for the turbine, the flow, the thrust and the power are available. It is found that the projection of the forces generated by the turbine into the flow solver grid is crucial for rotor predictions, especially for the power, and less important for the wake flow prediction. In this context, the projection of the forces into the flow solver grid should be as accurate as possible, in order to conserve the consistency between the computed axial velocity and the projected axial force. Also, the projection of the force is found to be much more important in the rotor plane directions than in the streamwise direction. It is found that for the case of a wind turbine immersed in a realistic turbulent boundary-layer flow, the potential spurious numerical oscillations originating from sharp force projections are not harmful to the results. By comparing an advanced model which computes the non-uniform distribution of the turbine forces over the rotor with a simple model which assumes uniform effects of the turbine forces, it is found that both can lead to accurate results for the far wake flow and the thrust and power predictions. However, the comparison shows that the advanced model leads to better results for the near wake flow. In addition, it is found that the simple model overestimates the rotor velocity prediction in comparison to the advanced model. These elements are explained by the lack of local feedback between the axial velocity and the axial force in the simple model. By comparing simulations with and without including the effects of the nacelle and tower, it is found that the consideration of the nacelle and tower is relatively important both for the near wake and the power prediction, due to the shadow effects. The grid resolution is not found to be critical once a reasonable resolution is used, i.e. in the order of 10 grid points along each direction across the rotor. The comparison with the experimental data shows that an accurate prediction of the flow, thrust, and power is possible with a very reasonable computational cost. Overall, the results give important guidelines for the implementation of ADMs for LES.


Sign in / Sign up

Export Citation Format

Share Document