Ambient Turbulence Research Articles

A New Streamwise Scaling for Wind Turbine Wake Modeling in the Atmospheric Boundary Layer

In this study, we aim to investigate if there is a scaling of the streamwise distance from a wind turbine that leads to a collapse of the mean wake velocity deficit under different ambient turbulence levels. For this purpose, we perform large-eddy simulations of the wake of a wind turbine under neutral atmospheric conditions with various turbulence levels. Based on the observation that a higher atmospheric turbulence level leads to faster wake recovery and shorter near-wake length, we propose the use of the near-wake length as an appropriate normalization length scale. By normalizing the streamwise distance by the near-wake length, we obtain a collapse of the normalized wake velocity deficit profiles for different turbulence levels. We then explore the possibility of using the relationship obtained for the normalized maximum wake velocity deficit as a function of the normalized streamwise distance in the context of analytical wake modeling. Specifically, we investigate two approaches: (a) using the new relationship as a stand-alone model to calculate the maximum wake velocity deficit, and (b) using the new relationship to calculate the wake advection velocity within a physics-based wake expansion model. Large-eddy simulation of the wake of a wind turbine under neutral atmospheric conditions is used to evaluate the performance of both approaches. Overall, we observe good agreement between the simulation data and the model predictions, along with considerable savings in terms of the models’ computational costs.

Observational Analysis and Numerical Modeling of the Solar Wind Fluctuation Spectra during Intervals of Plasma Instability

We perform a statistical analysis of observed magnetic spectra in the solar wind at 1 au with localized power elevations above the level of the ambient turbulent fluctuations. We show that the elevations are seen only when the intensity of the ambient fluctuations is sufficiently low. Assuming that the spectral elevations are caused by thermal-ion instabilities, this suggests that on average the effect of the solar wind background is strong enough to suppress the instability or obscure it or both. We then carry out nonlinear numerical simulations with particle ions and an electron fluid to model a thermal-ion instability coexisting with an ambient turbulence. The parameters of the simulation are taken from a known solar wind interval where an instability was assumed to exist based on the linear theory and a bi-Maxwellian fit of the observed distribution with core and secondary-beam protons. The numerical model closely matches the position of the observed spectral elevation in the wavenumber space. This confirms that the thermal-ion instability is responsible for the elevation. At the same time, the magnitude of the elevation turns out to be smaller than in the real solar wind. When higher intensity of the turbulence is used in the simulation, which is typical of solar wind in general, the power elevation is no longer seen. This is in agreement with the reduced observability of the elevations at higher intensities. However, the simulations show that the turbulence does not simply obscure the instability but also lowers its saturation level.

On the scalar turbulent/turbulent interface of axisymmetric jets

The effect of a zero-mean-flow turbulent ambient on the geometry of the turbulent/turbulent interface (TTI) of an axisymmetric jet is investigated and compared with the traditional turbulent/non-turbulent interface (TNTI). The ambient turbulence is generated using a random jet array. Orthogonal cross-sections of the jet subjected to different levels of ambient turbulence intensities and length scales are captured using planar laser-induced fluorescence. The enhanced radial transport of concentration from the jet core towards the interface and the increased scalar fluctuations within the jet caused by the ambient turbulence result in steeper mean and root-mean-square scalar conditional jumps across the TTI layer compared with those of the TNTI, respectively. When compared with the quiescent ambient, the mean effect of background turbulence is to stretch and corrugate the interface surface area as evident from the wider probability density of the interface radial position, the lowered occurrence of zero-curvature surface elements, the greater misalignment between the radial and normal unit vectors of the interface, the increased tortuosity and the increased magnitude of the fractal exponent.

Dependence of nonlinear coupling among turbulence, geodesic acoustic modes and tearing modes on magnetic island width in the HL-2A edge plasmas

Statistical spectral features of the dependence of geodesic acoustic modes (GAMs) and their nonlinear couplings with ambient turbulence on the magnetic island (MI) width (W) in the edge region of HL-2A tokamak plasmas are analyzed. Experimental observations have indicated that the modulation influence as well as the strength of nonlinear interactions between GAMs and turbulence generally shows a gradual decay while the couplings between MIs and the latter are increased simultaneously as the MI becomes larger. The MIs mainly reduce the couplings between GAMs and potential fluctuations, whereas the changes in the nonlinear interactions between density fluctuations and MIs are more evident. Moreover, it is found that there exists a nonmonotonic relationship between the turbulence correlation length and island width, in which it exhibits a minimum around W ∼ 3.7 cm, suggesting that the MI around such a scale would have a significant suppression effect on turbulent transport. These findings promote the understanding of the nonlinear interactions between MIs and turbulence in the edge of fusion plasmas.

A low-order wake interaction modeling framework for the performance of ocean current turbines under turbulent conditions

Understanding the effects of ambient turbulence (expressed often in terms of the turbulence intensity It) is critical to the development of predictive models for the performance of Ocean Current Turbine (OCTs). This paper describes a new, wake interaction modeling framework capable of capturing the detailed effects of turbulence on various performance parameters associated with OCTs that may be arranged in any arbitrary configuration. The model accounts for the effects of turbulence on the structure of the turbine wakes, specifically the extents of near- and far-wake regions, and the dependence of the transition point between the two regions on It. The analytical description for turbine wake is combined with an existing wake interaction model, the Unrestricted Wind Farm Layout Optimization (UWFLO) model to predict the global power output from an array of OCTs. The resulting modelling framework accurately captures the effect of inlet turbulence on the OCT farm power and efficiency, and can be applied to any array configuration. Results from the model are validated against Large Eddy Simulations (LES) in which the OCTs are modeled using the Blade Element Momentum (BEM) model, while the inlet flow is superposed with a synthetic turbulence field designed to approximate turbulence properties obtained from observational measurements of the Gulf Stream. The simulations show that OCT wakes recover faster at higher levels of inlet turbulence due to the enhanced entrainment and mixing between ambient flow and the wake, an effect that is captured by the modified UWFLO model.

An empirical model accounting for added turbulence in the wake of a full-scale turbine in realistic tidal stream conditions

The increased turbulence behind the turbine induces a fatigue load on the downstream turbine blades. An accurate estimate of turbulent intensity in turbine wake is paramount to optimizing the placement of turbines in tidal turbine parks. A simple empirical model is developed using data fitting from numerical simulation of a non-rotational actuator disc model to estimate the added turbulence in realistic tidal stream conditions similar to the Alderney Race, with a rotor diameter to depth ratio of 40%. The study shows a self-similar Gaussian shape streamwise turbulent intensity in a lateral direction similar to the velocity deficit profile. The turbulent wake radius expands according to a power-law depending on the ambient turbulent intensity. The added turbulence model justifies that the major turbulence source in the near wake is attributed to the rotor as it is weakly dependent on the ambient turbulence in the flow. As no known existing empirical model for tidal turbine wake added turbulence, the model is compared with existing models of unbounded wind turbines. The proposed model estimates the average turbulence in the far wake can assist turbine placement in a tidal farm.

Offshore wind resource assessment off the coast of Daejeong, Jeju Island using 30-year wind estimates

Long-term offshore wind resource assessment was performed for the offshore site off the Daejeong coast on Jeju Island of South Korea. A jacket type offshore met mast was installed at the site to measure wind and meteorological conditions for 1 year from 1 July 2015 to 30 June 2016. The wind conditions were measured at 32.5 m, 72.5 m and 100.5 m above mean sea level (AMSL). Fundamental wind characteristics such as average wind speed, wind direction frequency, turbulence intensity and so on were analyzed utilizing 1-year wind measurements. In order to correct the limited period of 1-year measurements to long-term conditions, the Measure-Correlate-Predict (MCP) method was applied to estimate 30-year wind conditions using nearby Modern-Era Retrospective Analysis for Research and Applications, Version 2 (MERRA-2) reanalysis data. Also, wind turbine classes suitable for the offshore site were determined based on extreme wind speed estimations and the measured ambient turbulence intensity. Five commercial wind turbines were chosen for estimation of the 30-year annual energy productions (AEPs) and capacity factors (CFs). As a result, it was estimated that the 30-year average wind speed was 7.86 m/s at 100.5 m AMSL. The wind turbine class II rating was estimated to be suitable for enduring extreme wind speed for a return period of 50 years. For the 30-year results, the minimum CF of 30.48% was estimated, while the maximum nearly reached 45%.

Isotope effects on turbulence and zonal flows in HL-2A edge plasmas

The effect of isotope mass on the interaction between turbulence and geodesic acoustic mode (GAM) zonal flows has been investigated in HL-2A ohmically heated deuterium (D) and hydrogen (H) plasmas using a double-step Langmuir probe array. The experimental results indicate that the level of GAM zonal flows and the turbulence eddy size together with the eddy tilting angle are all increased in the edge region in D plasmas compared to those in H plasmas under similar discharge parameters involving plasma current, magnetic field and line-averaged density. Evidence shows that in D plasmas, the nonlinear energy transfer is the main cause of the stronger excitation of GAM zonal flows, which extract more energy from ambient turbulence and, consequently, lead to lower turbulent transport and better confinement in D plasmas. The experimental findings may contribute to the understanding of the isotopic physics and associated turbulent transport in tokamak plasmas.

Assessment of offshore wind power including effects of sea surface roughness using observed wind statistics

The paper addresses how the offshore wind power can be assessed based on observed wind statistics at an offshore location. Results are exemplified for a Pierson-Moskowitz wave amplitude spectrum together with mean wind speed statistics from 11 sites in the North Sea and the North Atlantic. The ambient turbulence level, for example relevant for design of the wind turbine blades, is also estimated by adopting a wind gust spectrum. The sea surface roughness is specified using a formula given in terms of significant wave height and spectral wave steepness. The present method should be useful for making initial assessments of available wind power and ambient turbulence level at offshore sites based on wind statistics.

Impact of the mass isotope on plasma confinement and transport properties in the HL-2A tokamak

The impact of the mass isotope on plasma confinement and transport properties has been investigated in Ohmically-heated hydrogen and deuterium plasmas in the HL-2A tokamak. Experimental results show that under similar discharge parameters the deuterium plasma has better confinement and lower turbulent transport than the hydrogen one, and concomitantly, it is found that the magnitude of geodesic acoustic mode zonal flows, the tilting angle of the Reynolds stress tensor and the turbulence correlation lengths are all larger in the edge region of the deuterium plasma. The results provide direct experimental evidence on the importance of the nonlinear energy coupling between ambient turbulence and zonal flows for governing the isotope effects in fusion plasmas.

Combined three-dimensional flow field measurements and motion tracking of freely moving spheres in a turbulent boundary layer

A combination of time-resolved tomographic particle image velocimetry, refractive index matching technique and machine vision algorithms was used to measure the translational and rotational motion of freely moving, nearly neutrally buoyant spheres in a fully developed turbulent boundary layer (TBL). Located in the buffer and logarithmic layers, the hydrogel spheres ( $\sim$ 70 inner wall units in diameter) were refractive index matched with the water and tagged by ‘spokes’. Besides translational motion, the spheres exhibited significant rotation. The spheres were surrounded by typical coherent structures observed in TBLs, among them hairpin packets and transverse and longitudinal vortices that induced ejections and sweeps. While the majority of instantaneous sphere Reynolds numbers did not exceed 100, and vortex shedding was not observed, the results showed that the spheres may affect the evolution of hairpin packets in TBLs due to their finite size. The instantaneous rotation-, wall- and shear-induced lift forces, as well as the drag forces, acting on the spheres were estimated using available correlations for the lift and drag coefficients. Results hinted at negative shear-induced lift due to flow separation at a smaller critical Reynolds number than incorporated in the correlations that do not include the effect of ambient turbulence. The results indicated further that the drag force aided by the rotation-induced lift force was instrumental in keeping one of the spheres aloft. For the wall-ward moving spheres, lift forces opposed sphere motion. As a result, the spheres approached the wall with velocities lower than their quiescent settling velocity.

A physics-based model for wind turbine wake expansion in the atmospheric boundary layer

Analytical wind turbine wake models are widely used to predict the wake velocity deficit. In these models, the wake growth rate is a key parameter specified mainly with empirical formulations. In this study, a new physics-based model is proposed and validated to predict the wake expansion downstream of a turbine based on the incoming ambient turbulence and turbine operating conditions. The new model utilises Taylor diffusion theory, the Gaussian wake model, turbulent mixing layer theory and the analogy between wind turbine wake expansion and scalar diffusion. These components ensure that the model conserves mass and momentum in the far wake and accounts for the ambient turbulence and turbine-induced turbulence effects on the wake expansion. To account for the turbulence relevant scales that contribute to the wake expansion, the model uses the root-mean-square of the low-pass filtered radial velocity component. A simplified version that only requires the unfiltered velocity standard deviation and turbulence integral scale is also proposed. In addition, a new relation for the near-wake length is derived. The model performance is validated using large-eddy simulation data of a wind turbine wake under neutral atmospheric conditions with a wide range of incoming turbulence levels. The results show that the proposed model yields reasonable predictions of the wake width, maximum velocity deficit and near-wake length. In the case with a relatively low incoming streamwise turbulence intensity of 0.05, the ambient and turbine-induced terms in the model contribute almost equally to the wake width, rendering them both crucial for reasonable wake predictions.

Turbulence Optimized Park model with Gaussian wake profile

A new analytical wake model for wind farm design is presented. The new model, TurbOPark, is based on the Park model but features non-linear wake expansion and a Gaussian deficit profile. The modelled wake expansion depends on the intensity of the local turbulence, which is a combination of ambient atmospheric turbulence and wake-generated turbulence. Moving downstream from the rotor, the intensity of the wake-generated turbulence decreases, and the wake expansion slows down. This leads to more persistent wakes in wind farms and clusters of wind farms with large distances between turbines. In the present study, we focus on the London Array wind farm and show examples of how results from Park and TurbOPark compare to SCADA data. We look at the array efficiency of the wind farm, the production of individual turbines, and the impact of neighbouring wind farms. We find that predictions from TurbOPark generally are closer to observations than predictions from the Park model.

Tailoring wind turbine wake models to incoming free-stream turbulence

Analytical wake models describe the main wake statistics using a reduced number of parameters: for instance, the Gaussian wake model relies on the single parameter k ∗ to predict the mean wake velocity. However, the robustness of these models to different inflows has been seldom tested and the effects of ambient turbulence on these coefficients are reported only in a limited number of studies.With this objective, we measured the wake generated by a model-scale wind turbine operating in a wind tunnel at a ReD ⋍ 105 under different ambient conditions, with turbulence intensities between 3% and 12% and integral time scales between 1/10 and 10 times the convective time-scale. The wake generated by the turbine has been measured by means of planar PIV in a region spanning between 1.25 and 8.75 rotor diameters downstream of the turbine.The results show that a description of the turbine wake with the single parameter k ∗ is not sufficient to obtain accurate predictions, while the quality of the wake estimations can be dramatically improved by introducing a virtual origin x 0. Analysis of the higher-order wake statistics shows how the former can be related to the intensity of the Reynolds shear stress in the wake, and the latter to the erratic motion of the tip-vortices.

Turbulence above offshore wind farms measured by aircraft

The flow inside and around large offshore wind farms may range from smaller and shorter structures present due to the mechanical turbulence generated by wind turbines to larger, and longer structures indicative of the mesoscale flow. In this study, we explore the variation in turbulence structures above large offshore wind farms located in the North Sea using the data obtained from a research aircraft. The aircraft was flown upstream, downstream, and above the wind farm clusters. We evaluate the normalized autocorrelation of the wind velocity component in the mean wind direction for transverse separations and the rate of energy dissipation in the inertial subrange for neutral and stable stratification. During the neutral stratification, there is high ambient turbulence in the atmosphere and large values of the energy dissipation rate. The intensity of small-scale turbulence structures is increased above and downstream of the wind farm and it dominates the mesoscale flow structures. While in stable stratification, mesoscale flow structures are not only dominant upstream of the wind farm but also in the downstream flow. Apart from the atmospheric stratification, the wind farm layout also affects the turbulent structures, the rate of energy dissipation, and time for wake recovery.

The influence of rotor downwash on spray distribution under a quadrotor unmanned aerial system

• Computational Fluid Dynamics model of sprays interaction with propeller airflow. • CFD model for multi nozzle/propeller crop sprayer, crosswind and forward flight. • Experimental validation of CFD model using a spray pattonator. • CFD model can be used to model the distribution of the spray’s impact with the crop. • Modelled distributions can be used for spray mission planners. • Spray distribution heavily effected by propeller downwash. • Higher propeller thrust leads to a peakier and less uniform distribution. This paper investigates how sprays are influenced by rotor downwash produced by a quadrotor Unmanned Aerial System when crop spraying. Computational Fluid Dynamics simulations, using RANS modelling for the flow field and Lagrangian particle tracking for the spray, are conducted for several single-rotor and multi-rotor cases with the spray injected underneath the centre of the rotor. The rotor is modelled using a Blade Element actuator disc model. The accuracy of the computational approach is demonstrated by good agreement with experiment for both thrust and deposited spray pattern for a single rotor in an indoor environment. Both the experimental and computational results show that the peak in the spray distribution under a single rotor increases as rotor thrust increases, whilst increasing the rotor height above the ground causes this peak to decrease. The validated Computational Fluid Dynamics method is then used to simulate flight conditions for single-rotor and multi-rotor cases. These show the existence of a critical flight speed, above which the spray impinging on the ground no longer contains a notable peak. This behaviour is seen to be due to the streamtube detaching from the ground and no longer carrying the spray directly to the ground plane. Above this critical speed the spray is seen to become suspended in the air behind the Unmanned Aerial System. This behaviour makes realistic simulation more difficult, as details of the ambient turbulence conditions would be needed to model the subsequent spray transport. The reliance on turbulence to transport the spray is undesirable from a practical point of view due to the increased likelihood of significant spray drift.

Wake interactions of two horizontal axis tidal turbines in tandem

The wake interactions between two horizontal axis model tidal turbines for two different spacings were investigated by flow measurements using acoustic Doppler velocimetry (ADV). The analysis of the velocity time-series data showed that a velocity sampling time period of at least 1000 integral time scales is required for the convergence of turbulence statistics in the wake of the turbines. The velocity measurement showed that a distinct three-layer (core, inner, and outer) structure was observed in the wake of a single turbine, whereas a simple wake pattern without such a structure was observed downstream of the two turbines. The comparison between the measurements for the different turbine spacings showed that the mean velocity fields, turbulent stress fields, velocity spectra plots, and integral length scales behind the second turbine were nearly independent of the turbine spacing. It was concluded that the increased turbulent momentum transfer across the wake boundary, mainly due to the wake meandering motion by the first turbine, overshadows the effect of the first turbine. This suggests that the ambient turbulence intensity induced by the upstream turbine could be a key factor for the wake characteristics within a tidal turbine array.

The dynamics of an axisymmetric turbulent jet in ambient turbulence interpreted from the passive scalar field statistics

The effect of approximately homogeneous isotropic turbulence on the dynamics of an axisymmetric turbulent jet (Re = 10 600 and 5800) in an ambient with a negligible mean flow is interpreted from the statistics of the passive scalar field. The ambient turbulence is generated by a random jet array and scalar concentrations are measured in orthogonal cross sections of the jet using planar laser-induced fluorescence. Statistics of the scalar field of the jet in a turbulent ambient are compared to those in a quiescent ambient, using classical Eulerian averages and those conditioned on the jet centroid. A two-region model for the jet structure in ambient turbulence is proposed based on the centroidal statistics. Following the developing region of the jet, the ambient turbulence disrupts the jet structure, due to modulation of the jet interface, meandering of the jet by large eddies, and entrainment of the turbulent ambient fluid, resulting in a faster concentration decay and reduced entrainment compared to the quiescent ambient. Further downstream, once the ambient turbulence has destroyed the jet, only molecular and turbulent diffusion modify the scalar concentrations. The regions' relative lengths depend primarily on the relative turbulence intensity (ξ) between the ambient and the jet, as assessed using centroidal analysis, which removes the effect of the relative length scale (l) on the jet behavior in the turbulent ambient. The centroidal scalar statistics reveal self-similarity and self-preservation in mean scalar properties before jet breakup, which occurs abruptly once ξ> 0.5. The smaller scales of ambient turbulence modulate the jet boundary and, when entrained, lead to a wider range of centerline concentrations and root mean square concentrations, which are hypothesized to increase local concentration gradients within the jet and reduce jet mixing.

Atmospheric turbulent structures and fire sweeps during shrub fires and implications for flaming zone behaviour

Background Wildfires propagate through vegetation exhibiting complex spread patterns modulated by ambient atmospheric wind turbulence. Wind gusts at the fire-front extend and intensify flames causing direct convective heating towards unburnt fuels resulting in rapid acceleration of spread. Aims To characterise ambient and fire turbulence over gorse shrub and explore how this contributes to fire behaviour. Methods Six experimental burns were carried out in Rakaia, New Zealand under varying meteorological conditions. The ignition process ensured a fire-line propagating through dense gorse bush (1 m high). Two 30-m sonic anemometer towers measured turbulent wind velocity at six different levels above the ground. Visible imagery was captured by cameras mounted on uncrewed aerial vehicles at 200 m AGL. Key results Using wavelet decomposition, we identified different turbulent time scales that varied between 1 and 128 s relative to height above vegetation. Quadrant analysis identified statistical distributions of atmospheric sweeps (downbursts of turbulence towards vegetation) with sustained events emanating from above the vegetation canopy and impinging at the surface with time scales up to 10 s. Conclusions Image velocimetry enabled tracking of ‘fire sweeps’ and characterised for the first time their lifetime and dynamics in comparison with overlying atmospheric turbulent structures. Implications This methodology can provide a comprehensive toolkit when investigating coupled atmosphere–fire interactions.

The Effect of Solar Wind Turbulence on Parallel and Oblique Firehose Instabilities

Abstract We consider the firehose instability coexisting with the omnipresent ambient solar wind turbulence. The characteristic temporal and spatial scales of the turbulence are comparable to those of the instability. Therefore, turbulence may violate the common assumption of a uniform and stationary background used to describe instabilities and make the properties of the instabilities different. To investigate this effect, we perform three-dimensional hybrid simulations with particle-in-cell ions and a quasi-neutralizing electron fluid. We find that the turbulence significantly reduces the growth rates and saturation levels of both instabilities. Comparing the cases with and without turbulence, the former results in a higher temperature anisotropy in the asymptotic marginally stable state at large times. In the former case, the distribution function averaged over the simulation box is also closer to the initial one.