Coherent Turbulence Research Articles

Demonstrating 24-hour continuous vertical monitoring of atmospheric optical turbulence.

We report what is believed to be the first example of fully continuous, 24-hour vertical monitoring of atmospheric optical turbulence. This is achieved using a novel instrument, the 24-hour Shack-Hartmann Image Motion Monitor (24hSHIMM). Optical turbulence is a fundamental limitation for applications such as free-space optical communications, where it limits the achievable bandwidth, and ground-based optical astronomy, restricting the observational precision. Knowledge of the turbulence enables us to select the best sites, design optical instrumentation and optimise the operation of ground-based optical systems. The 24hSHIMM estimates the vertical optical turbulence coherence length, time, angle and Rytov variance from the measurement of a four-layer vertical turbulence profile and a wind speed profile retrieved from meteorological forecasts. To illustrate our advance we show the values of these parameters recorded during a 36-hour, continuous demonstration of the instrument. Due to its portability and ability to work in stronger turbulence, the 24hSHIMM can also operate in urban locations, providing the field with a truly continuous, versatile turbulence monitor for all but the most demanding of applications.

A Feasibility Study of Riblet for Aeroacoustics Applications

Although the use of finlet, serration or porous surface has been shown a good level of success in reducing aerofoil trailing-edge noise, they are largely incompatible to the otherwise streamlined aerofoil body. This paper is a feasibility study to investigate the riblet, which has so far been quite successful as a drag-reducing device, for its potential to reduce the turbulent pressure sources that are important for aerofoil self-noise radiation. The completed results show that the riblet used in the current study can reduce the skin-friction coefficient, as well as the turbulence-intensity in the boundary-layer profiles. In addition, the turbulence structures in the convective field can be dissipated more rapidly when crossing the riblet surface. It is also found that (1) the riblet produces a slight reduction of the wall-pressure power spectral density level at the low and high frequency ranges, but experiences an increase at the mid frequency, (2) the riblet can reduce the lateral turbulence coherence length-scale across a large frequency range. The product of these two hydrodynamic sources represents an important mechanism for the radiation of the trailing-edge self-noise, whose low and high frequency ranges are found to be sensitive to the riblet effect where reduction can occur.

Analytical Eddy Viscosity Model for Turbulent Wave Boundary Layers: Application to Suspended Sediment Concentrations over Wave Ripples

Turbulence related to flow oscillations near the seabed, in the wave bottom boundary layer (WBBL), is the phenomenon responsible for the suspension and transport of sediments. The vertical distribution of turbulent eddy viscosity within the WBBL is a key parameter that determines the vertical distribution of suspended sediments. For practical coastal engineering applications, the most used method to parameterize turbulence consists in specifying the shape of the one-dimensional-vertical (1DV) profile of eddy viscosity. Different empirical models have been proposed for the vertical variation of eddy viscosity in the WBBL. In this study, we consider the exponential-type profile, which was validated and calibrated by direct numerical simulation (DNS) and experimental data for turbulent channel and open-channel flows, respectively. This model is generalized to the WBBL, and the period-averaged eddy viscosity is calibrated by a two-equation baseline (BSL) k-ω model for different conditions. This model, together with a β-function (where β is the inverse of the turbulent Schmidt number), is used in modeling suspended sediment concentration (SSC) profiles over wave ripples, where field and laboratory measurements of SSC show two kinds of concentration profiles depending on grain particles size. Our study shows that the convection–diffusion equation, for SSC in WBBLs over sand ripples with an upward convection term, reverts to the classical advection–diffusion equation (ADE) with an “apparent” sediment diffusivity εs*=α εs related to the sediment diffusivity εs by an additional parameter α associated with the convective sediment entrainment process over sand ripples, which is defined by two equations. In the first, α depends on the relative importance of upward convection related to coherent vortex shedding and downward settling of sediments. When the convective transfer is very small, above low-steepness ripples, α≈1. In the second, α depends on the relative importance of coherent vortex shedding and random turbulence. When random turbulence is more important than coherent vortex shedding, α≈1, and “apparent” sediment diffusivity reverts to the classical sediment diffusivity εs*≈ εs. Comparisons with experimental data show that the proposed method allows a good description of both SSC for fine and coarse sand and “apparent” sediment diffusivity εs*.

Wavelet applications for turbulence characterisation of real tidal flows measured with an ADCP

Understanding turbulence is crucial to the design of tidal energy converters as it influences loads, fatigue life and power production. Acoustic Doppler Current Profilers (ADCPs) are normally employed in measurement campaigns to analyse turbulence at a particular location. Coherent turbulent structures are known to lead to the highest turbulence-induced stresses and are typically analysed by Fourier transform. However, this method can only yield parameters which are averaged over the signal length or the analysing window, and hence may not fully represent real flow turbulence, which is often non-stationary. This study applies wavelet time–frequency analysis to ADCP data from an active tidal test site to examine the turbulence spectrum and coherence characteristics of the flow, comparing these to Fourier methods. Results show that the average quantities from wavelet analysis agree favourably with Fourier analysis. Moreover, wavelet analysis yields instantaneous spectra and coherence, which show that at this site, high energy bursts occur at frequencies known to be significant for TEC loading and have turbulence intensities up to 80% higher than the average. We show that these intermittent, coherent bursts are obscured by the averages associated with Fourier analysis, highlighting the need for methods such as wavelets to understand interactions between fluid flow and TEC.

Dynamic Mechanisms Associated with the Structure and Evolution of Roll Vortices and Coherent Turbulence in the Hurricane Boundary Layer: A Large Eddy Simulation During the Landfall of Hurricane Harvey

Roll vortices are a series of large-scale turbulent eddies that nearly align with the mean wind direction and prevail in the hurricane boundary layer. In this study, the one-way nested WRF-LES model simulation results from Li et al. (J Atmos Sci 78(6):1847–1867, https://doi.org/10.1175/JAS-D-20-0270.1, 2021) are used to examine the structure and generation mechanism of roll vortices and associated coherent turbulence in the hurricane boundary layer during the landfall of Hurricane Harvey from 00 UTC 25 to 18 UTC 27 August 2017. Results indicate that roll vortices prevail in the hurricane boundary layer. The intense roll vortices and associated large turbulent eddies above them (at a height of ~ 200 to 3000 m) accumulate within a hurricane radius of 20–40 km. Their intensity is proportional to hurricane intensity during the simulation period. Before and during hurricane landfall, strong inflow convergence leads to horizontal advection of roll vortices throughout the entire hurricane boundary layer. Combined with the strong wind shear, the strongest roll vortices and associated large turbulent eddies are generated near the eyewall with suitable thermodynamic (Richardson number at around − 0.2 to 0.2) and dynamic conditions (strong negative inflow wind shear). After landfall, the decayed inflow weakens the inflow convergence and quickly reduces the strong roll vortices and associated large turbulent eddies. Diagnosis of vertical turbulent kinetic energy indicates that atmospheric pressure perturbation, caused by horizontal convergence, transfers the horizontal component of turbulence to the vertical component with a mean wavelength of about 1 km. The buoyancy term is weak and negative, and the large turbulent eddies are suppressed.

Coherent structures, turbulence intermittency, and anisotropy of gravity currents propagating on a rough and porous bed

This study experimentally investigated the impacts of rough and porous (RP) bed and sedimentation processes on the coherent structures, turbulence intermittency, and anisotropy of saline and turbidity currents. The results reveal that the local current concentration responds immediately (saline current) or languidly (turbidity current) to turbulence bursting events. Inside the dense current, the turbulent momentum fluxes in the streamwise (u′u′¯) and vertical (w′w′¯) directions transfer downstream and downward, which favor the sweep events. Inside the ambient water, u′u′¯ and w′w′¯ transfer upstream and upward, contributing to the formation of ejection events. At the current–ambient water interface, u′u′¯ and w′w′¯ do not tend to transfer in particular directions resulting in almost equal quantities of sweep and ejection events. The Gram–Charlier series expansion is strictly applicable to probability density functions (PDFs) of the sweep and ejection events but not suitable ideally for PDFs of the outward and inward interaction events. The primary anisotropy invariant map (AIM) of gravity currents starts from the two-component plain strain limit (near the bed). It is followed by the three-dimensional isotropy (inside the dense current and ambient water) and the axisymmetric contraction limit (current–ambient water interface). Finally, it ends in two-dimensional isotropy (near the free surface). This AIM is sensitive to the RP boundary and the sedimentation processes. Along the streamwise direction, the RP boundary causes alternations between the anisotropic and isotropic turbulence, but the arranged pattern of the rough units determines the period of this alternation.

Large-Scale Coherent Turbulence Structures in the Atmospheric Boundary Layer over Flat Terrain

Abstract We investigate characteristics of large-scale coherent motions in the atmospheric boundary layer using field measurements made with two long-range scanning wind lidars. The joint scans provide quasi-instantaneous wind fields over a domain of ∼50 km2, at two heights above flat but partially forested terrain. Along with the two-dimensional wind fields, two-point statistics and spectra are used to identify and characterize the scales, shape, and anisotropy of coherent structures—as well as their influence on wind field homogeneity. For moderate to high wind speeds in near-neutral conditions, most of the observed structures correspond to narrow streaks of low streamwise momentum near the surface, extending several hundred meters in the streamwise direction; these are associated with positive vertical velocity ejections. For unstable conditions and moderate winds, these structures become large-scale rolls, with longitudinal extent exceeding the measuring domain (>∼5 km); they dominate the conventional surface-layer structures in terms of both physical scale and relative size of velocity-component variances, appearing as quasi-two-dimensional structures throughout the entire boundary layer. The observations shown here are consistent with numerical simulations of atmospheric flows, field observations, and laboratory experiments under similar conditions. Significance Statement Coherent structures have attracted the interest of researchers for decades, being viewed as the closest to “order” that we can find within the chaos of turbulence. In the turbulent atmospheric boundary layer, micro- and mesoscale coherent structures come in many shapes and sizes, such as convective cells, rolls, or streaks. In this study we used dual lidars (remote sensing measurements), developing analysis of their tandem usage to characterize in detail some of the large-scale coherent structures generated over flat terrain. This allowed us to better understand the mechanisms that generate such structures and describe their influence on the morphology of the turbulent atmospheric boundary layer across a good deal of its depth.

Unfrozen Skewed Turbulence for Wind Loading on Structures

The paper introduces an algorithm to generate a three-variate four-dimensional wind turbulence field suited for yawed wind dynamic load simulation. At large yaw angles, a relaxation of Taylor’s hypothesis of frozen turbulence becomes relevant as well as the flow phase lag in the along-wind direction, which modulates the real and imaginary parts of the coherence. To capture such a general wind action on a structure, a modified spectral representation method is used where the coherence of turbulence is described as a complex-valued function. The one-point and two-point co-spectra are implemented in the simulation setup using a square-root-free Cholesky decomposition of the spectral matrix. The numerical procedure is illustrated based on turbulence characteristics derived from data collected during storm Aina (2017) on the Norwegian coast by three-dimensional sonic anemometers. During this event, a remarkable 3-hour stationary time series with a mean wind speed of 24 m s−1 at a height of 49 m above ground was recorded. Since no computational grid is needed, the velocity fluctuations with representative spatio-temporal characteristics can be directly simulated on structural elements of slender structures. Such an algorithm may be essential for the design of super-long span bridges in coastal areas.

Mitigation of turbulent noise sources by riblets

A feasibility study is presented on the application of riblets to suppress turbulent pressure sources leading to a reduction in the generation of aerofoil self-noise. It is shown that riblets can reduce skin friction as well as the turbulence intensity inside the boundary layer. In addition, near-wall turbulence structures were found to be dissipated quite rapidly when crossing the riblets surface. It is found that riblets (1) slightly reduce the wall fluctuating pressure power spectral density level at the low and high frequency ranges, but cause an increase at the mid frequency range, and (2) reduce the lateral turbulence coherence length scale across a large frequency range. As a result, riblets have the potential to reduce trailing edge noise at the low and high frequency regions. The fundamental mechanism by which riblets reduce the turbulent intensity level inside the boundary layer is investigated by studying the spatio-temporal evolution of turbulent spots, which are commonly regarded as the building blocks of a turbulent boundary layer. In this paper two mechanisms are proposed. First, the enhanced momentum achieved at the becalmed region of each turbulent spot will lead to a reduction in the overall turbulence intensity obtained when turbulent spots merge downstream. Second, the internal turbulence level at the rear part of a turbulent spot can be reduced directly by an enhanced re-laminarisation effect.

Effects of atmospheric stability on the structural response of a 12 MW semisubmersible floating wind turbine

AbstractGlobal dynamic response models used for the design of wind turbines are largely based on neutral stability, which is not representative of real atmospheric conditions. Offshore wind farms, for example, have been seen to experience predominantly unstable conditions, especially at lower wind speeds. In the current work, we use four wind generation models under stable, neutral and unstable atmospheric stability conditions to study the low‐frequency content of the global responses of a semisubmersible floating wind turbine (FWT). To represent the wind fields, we use the Kaimal Spectrum and Exponential Coherence Model (Kaimal), the Mann Spectral Tensor Model (Mann), a point measurement based model (TIMESR) and large‐eddy simulation (LES). At the low‐frequency range, both atmospheric stability and the turbulence wind model significantly affect the response of the FWT. In all the cases studied throughout this paper, the structural response under unstable conditions is higher than under stable or neutral conditions. The TIMESR and the Kaimal models fitted to the FINO‐1 offshore meteorological mast measurements show more similar responses than the Mann model; surge and pitch are higher for the TIMESR and Kaimal models, and yaw is lower. When fitted to LES, TIMESR and Kaimal predict surge and pitch responses closer to LES, but they underestimate the responses related to yaw, opposite to what the Mann model does. The responses are directly related to turbulence intensity and coherence, which are affected by atmospheric stability. Therefore, based on the analyses carried out through this study, the structural analysis of FWT should account for the effect of atmospheric stability.

Large-eddy simulation of supercritical free-surface flow in an open-channel contraction

Large-eddy simulations (LES) of supercritical flow in a straight-wall, open-channel contraction of 6° and contraction ratio of 2:1 are performed. The LES code solves the filtered Navier-Stokes equations for two-phase flows (water-air) and employs the level-set method. The simulation was validated by replicating a previously reported experiment. Contours of the time-averaged velocities indicate that the flow loses energy and momentum in the contracting channel. Further, secondary currents in the contraction are redistributing momentum and are responsible for local up-and down-flows. The turbulent kinetic energy reaches very high values at the entrance of the contraction, mainly contributed by the streamwise normal stress. The flow contains coherent turbulence structures which are responsible for carrying low-momentum from the bed and the water surface towards the channel centre. Flow deceleration results in significant turbulence anisotropy in the contracted section. It is shown that mainly pressure drag contributes to the energy loss in the contraction.

Intermittency of Kolmogorov and Coherent Turbulence in the Mountain Atmospheric Boundary Layer (Review)

Intermittency of Kolmogorov and Coherent Turbulence in the Mountain Atmospheric Boundary Layer (Review)

Orbital angular momentum evolution of twisted multi-Gaussian Schell model beams in anisotropic turbulence

Twist phase is an unconventional phase that allows a partially coherent beam to carry orbital angular momentum (OAM). In this work, based on the Mercer’s expansion and the nonnegativity of the cross-spectral density function, we obtain a twistable condition for the multi-Gaussian Schell model (MGSM) correlation structure. The expression of an elliptical twisted multi-Gaussian Schell model (TMGSM) beam through non-Kolmogorov anisotropic atmospheric turbulence is rigorously derived, and its spectral density and OAM variation in anisotropic turbulence are studied in detail. The results show that the spectral density and OAM flux density distributions rotate with increasing transmission distance, and increasing the multi-Gaussian modulus and twist factor can significantly suppress the negative effects of atmospheric turbulence. The difference of turbulence anisotropy, beam modulus, and coherence length can affect the distribution of OAM. Furthermore, the twist phase affects both the OAM distribution and magnitude. Our work is helpful for exploring new forms of OAM sources, and promote the application of free-space optical communications.

Effect of secondary currents on the flow and turbulence in partially filled pipes

Large-eddy simulations of turbulent flow in partially filled pipes are conducted to investigate the effect of secondary currents on the friction factor, first- and second-order statistics and large-scale turbulent motion. The method is validated first and simulated profiles of the mean streamwise velocity, normal stresses and turbulent kinetic energy (TKE) are shown to be in good agreement with experimental data. The secondary flow is stronger in half- and three-quarters full pipes compared with quarter full or fully filled pipe flows, respectively. The origin of the secondary flow is examined by both the TKE budget and the steamwise vorticity equation, providing evidence that secondary currents originate from the corner between the free surface and the pipe walls, which is where turbulence production is larger than the sum of the remaining terms of the TKE budget. An extra source of streamwise vorticity production is found at the free surface near the centreline bisector, due to the two-component asymmetric turbulence there. The occurrence of dispersive stresses (due to secondary currents) reduces the contribution of the turbulent shear stress to the friction factor, which results in a reduction of the total friction factor of flows in half and three-quarters full pipes in comparison to a fully filled pipe flow. Furthermore, the presence of significant secondary currents inhibits very-large-scale motion (VLSM), which in turn reduces the strength and scales of near-wall streaks. Subsequently, near-wall coherent structures generated by streak instability and transient growth are significantly suppressed. The absence of VLSM and less coherent near-wall turbulence structures is supposedly responsible for the drag reduction in partially filled pipe flows relative to a fully filled pipe flow at an equivalent Reynolds number.

Study on flow characteristics and diversity index of diamond-type boulder cluster with different spacing ratios.

The placement of boulder or boulder cluster in rivers can increase or repair the complexity of river structure and the diversity of hydraulic conditions, which is very important for the habitat of many aquatic organisms. In this study, the diamond-type boulder cluster was modeled as four hemispheres exposed to a fully developed turbulent open channel flow. Numerical simulation was conducted to investigate the time-averaged flow characteristics, three-dimensional coherent structures, turbulence characteristics, and flow diversity index at different spacing ratios L/D (the ratio of the distance L to the diameter D, 1.0 ≤ L/D ≤ 3.5, where L is the center-to-center distance between two adjacent hemispheres and D is the diameter of the hemisphere). The results show that with the increase of the spacing ratio, the shear layer on the side of the gap flow gradually strengthens, and the single Karman vortex street in the wake region of the hemisphere array is suppressed. The time-averaged peak velocity in the gap flow gradually decreases with the increase of the spacing ratio, and the single of the recirculation zone behind the hemisphere array transforms into the recirculation zone behind each hemisphere, and the length of the each recirculation zone increases to the same. The turbulence intensity of the array first increases with the increase of the spacing ratio and then gradually decreases to a constant, reaching the peak intensity at L/D = 2. Based on the Shannon entropy concept, the flow diversity index in the zone of influence (ZOI) is calculated by considering the velocity and turbulence kinetic energy. The flow diversity index is the largest in the ZOI at the spacing ratio of 1.5.

Optimization of incoherent beam combining in turbulent atmosphere

Incoherent beam combining possesses many advantages not only in the design and control of system structure, but also in the propagation under some condition. In present paper we investigate analytically the propagation of coherent and incoherent combined beam. Results show that when the maximum on-axis intensity is selected as the criterion, the optimization of the number of the channels for incoherent beam combining is obtained according to the transmitter aperture size and the coherence length of turbulence. The number of the channels for incoherent beam combining has an optimal value when the transmitter and turbulence coherence length is unchanged. The propagation of incoherent beam combining is better, at least not worse than that of coherent beam combining under some conditions. The investigation will help to select and optimize the parameters of the incoherent beam combining in some applications.

Direct numerical simulation of the flow around a sphere immersed in a flat-plate turbulent boundary layer

This study investigates flow past a sphere immersed in a flat-plate turbulent boundary layer by using direct numerical simulations, with the objective of clarifying the effects of a wall-proximity sphere on turbulent coherent structures and turbulence statistics. Three cases are evaluated with gap ratios (G/D) of 0.1, 0.25, and 0.5 at a Reynolds number of ReD=2500, where D is the diameter of the sphere and G is the gap width between the bottom of the sphere and the flat plate. The results show that the wake of the sphere plays an important role in the streamwise region 0<x/D<10. The near-wall streaks break into small-scale point-like or patch-like structures in the near-wake region, with the most significant effect at G/D=0.1. This can be attributed to the interactions between the shedding vortex behind the sphere and vortical structures within the flat-plate boundary layer. Detail analysis of turbulence statistics indicates that the flat-plate boundary layer thickness is increased at x/D<−0.5 owing to the blockage effect of the sphere, whereas it is decreased at x/D>0.5 because of the high-speed jet flow around it. In addition, the presence of a wall-proximity sphere significantly affects the skin friction coefficient. The budgets of the turbulent kinetic energy show that turbulence production and viscous dissipation are augmented due to the formation of small-scale vortices and interactions among them in the near-wake region, especially at G/D=0.1.

Data-driven enhancement of coherent structure-based models for predicting instantaneous wall turbulence

Predictions of the spatial representation of instantaneous wall-bounded flows, via coherent structure-based models, are highly sensitive to the geometry of the representative structures employed by them. In this study, we propose a methodology to extract the three-dimensional (3-D) geometry of the statistically significant eddies from multi-point wall-turbulence datasets, for direct implementation into these models to improve their predictions. The methodology is employed here for reconstructing a 3-D statistical picture of the inertial wall coherent turbulence for all canonical wall-bounded flows, across a decade of friction Reynolds number (Reτ). These structures are responsible for the Reτ-dependence of the skin-friction drag and also facilitate the inner-outer interactions, making them key targets of structure-based models. The empirical analysis brings out the geometric self-similarity of the large-scale wall-coherent motions and also suggests the hairpin packet as the representative flow structure for all wall-bounded flows, thereby aligning with the framework on which the attached eddy model (AEM) is based. The same framework is extended here to also model the very-large-scaled motions, with a consideration of their differences in internal versus external flows. Implementation of the empirically-obtained geometric scalings for these large structures into the AEM is shown to enhance the instantaneous flow predictions for all three velocity components. Finally, an active flow control system driven by the same geometric scalings is conceptualized, towards favourably altering the influence of the wall coherent motions on the skin-friction drag.

Acoustics of rotor ingesting an axisymmetric turbulent boundary layer

The noise generation by a five-bladed rotor ingesting an axisymmetric turbulent boundary layer at the tail end of a body of revolution (BOR) is investigated numerically at Reynolds number of 1.9 × 106 based on the BOR length and free-stream Mach number of 0.059. The turbulent boundary layer on the nose and midsection of the BOR is computed using wall-modeled LES whereas that in the acoustically important tail-cone section under a strong adverse pressure gradient is wall resolved. The radiated acoustic field is calculated based on the Ffowcs Williams-Hawkings equation. The computed flow statistics and sound pressure spectra agree well with the experimental measurements at Virginia Tech. In addition to broadband turbulence-ingestion noise, haystacking spectral peaks near the blade passing frequency and its harmonics are well captured and related to correlated blade interactions with coherent turbulence structures. The accuracy of the classical Sears theory for rotor-noise prediction is evaluated using the upwash velocity field from the LES. The theory is shown to predict well the sound pressure levels at low to mid-frequencies including the haystacking peaks but less accurately at higher frequencies. The spatial and frequency characteristics of blade unsteady-loading dipoles and their correlations are investigated, and the appropriate Mach number scaling for the ingestion noise is identified.

The COTUR project: remote sensing of offshore turbulence for wind energy application

Abstract. The paper presents the measurement strategy and data set collected during the COTUR (COherence of TURbulence with lidars) campaign. This field experiment took place from February 2019 to April 2020 on the southwestern coast of Norway. The coherence quantifies the spatial correlation of eddies and is little known in the marine atmospheric boundary layer. The study was motivated by the need to better characterize the lateral coherence, which partly governs the dynamic wind load on multi-megawatt offshore wind turbines. During the COTUR campaign, the coherence was studied using land-based remote sensing technology. The instrument setup consisted of three long-range scanning Doppler wind lidars, one Doppler wind lidar profiler and one passive microwave radiometer. Both the WindScanner software and LidarPlanner software were used jointly to simultaneously orient the three scanner heads into the mean wind direction, which was provided by the lidar wind profiler. The radiometer instrument complemented these measurements by providing temperature and humidity profiles in the atmospheric boundary layer. The scanning beams were pointed slightly upwards to record turbulence characteristics both within and above the surface layer, providing further insight on the applicability of surface-layer scaling to model the turbulent wind load on offshore wind turbines. The preliminary results show limited variations of the lateral coherence with the scanning distance. A slight increase in the identified Davenport decay coefficient with the height is partly due to the limited pointing accuracy of the instruments. These results underline the importance of achieving pointing errors under 0.1∘ to study properly the lateral coherence of turbulence at scanning distances of several kilometres.