Turbulent Coherent Structures Research Articles

Turbulent particle-gas feedback exacerbates the hazard impacts of pyroclastic density currents

Causing one-third of all volcanic fatalities, pyroclastic density currents create destruction far beyond our current scientific explanation. Opportunities to interrogate the mechanisms behind this hazard have long been desired, but pyroclastic density currents persistently defy internal observation. Here we show, through direct measurements of destruction-causing dynamic pressure in large-scale experiments, that pressure maxima exceed theoretical values used in hazard assessments by more than one order of magnitude. These distinct pressure excursions occur through the clustering of high-momentum particles at the peripheries of coherent turbulence structures. Particle loading modifies these eddies and generates repeated high-pressure loading impacts at the frequency of the turbulence structures. Collisions of particle clusters against stationary objects generate even higher dynamic pressures that account for up to 75% of the local flow energy. To prevent severe underestimation of damage intensities, these multiphase feedback processes must be considered in hazard models that aim to mitigate volcanic risk globally.

Spectral proper orthogonal decomposition of harmonically forced turbulent flows

Many turbulent flows exhibit time-periodic statistics. These include turbomachinery flows, flows with external harmonic forcing and the wakes of bluff bodies. Many existing techniques for identifying turbulent coherent structures, however, assume the statistics are statistically stationary. In this paper, we leverage cyclostationary analysis, an extension of the statistically stationary framework to processes with periodically varying statistics, to generalize the spectral proper orthogonal decomposition (SPOD) to the cyclostationary case. The resulting properties of the cyclostationary SPOD (CS-SPOD for short) are explored, a theoretical connection between CS-SPOD and the harmonic resolvent analysis is provided, simplifications for the low and high forcing frequency limits are discussed, and an efficient algorithm to compute CS-SPOD with SPOD-like cost is presented. We illustrate the utility of CS-SPOD using two example problems: a modified complex linearized Ginzburg–Landau model and a high-Reynolds-number turbulent jet.

The structure of turbulence in unsteady flow over urban canopies

The topology of turbulent coherent structures is known to regulate the transport of energy, mass and momentum in the atmospheric boundary layer (ABL). While previous research has primarily focused on characterizing the structure of turbulence in stationary ABL flows, real-world scenarios frequently deviate from stationarity, giving rise to nuanced and poorly understood changes in the turbulence geometry and associated transport mechanisms. This study sheds light on this problem by examining topological changes in ABL turbulence induced by non-stationarity and their effects on momentum transport. Results from a large-eddy simulation of pulsatile open channel flow over an array of surface-mounted cuboids are examined. The analysis reveals that the flow pulsation triggers a phase-dependent shear rate, and the ejection-sweep pattern varies with the shear rate during the pulsatile cycle. From a turbulence structure perspective, it is attributed to the changes in the geometry of hairpin vortices. An increase (decrease) in the shear rate intensifies (relaxes) these structures, leading to an increase (decrease) in the frequency of ejections and an amplification (reduction) of their percentage contribution to the total momentum flux. Furthermore, the size of the hairpin packets undergoes variations, which depend on the geometry of the constituting hairpin vortices, yet the packet inclination preserves its orientation throughout the pulsatile cycle. These observations reinforce the important role non-stationarity holds in shaping the structure of ABL turbulence and the momentum transport mechanisms it governs.

Lagrangian and Eulerian perspectives of turbulent transport mechanisms in a lateral cavity

The dynamics of turbulent flows past lateral cavities is relevant for multiple environmental applications. In rivers and coastal environments, these lateral recirculating regions constitute surface storage zones, where large-scale turbulent coherent structures control the transport and fate of contaminants. Mass transport in these flows is typically represented by one-dimensional first-order equations that predict the evolution of the spatially integrated concentration between the cavity and the main channel. These models, however, cannot represent the long-term evolution of the concentration or incorporate memory effects induced by turbulence. In this investigation, we carry out large-eddy simulations (LES) of the open-channel flow with a lateral square cavity of Mignot et al. (Phys. Fluids, vol. 28, issue 4, 2016, 045104). The model is coupled with an advection–diffusion equation and a Lagrangian particle model to investigate the transport mechanisms in the cavity and across the interface. From the simulations we provide quantitative comparisons of the physical processes from both perspectives, and investigate the effects of turbulent coherent structures on residence times and trajectories from finite-time Lyapunov exponents. From the Lagrangian results, we identify general spatial distributions of time scales in the cavity associated with the dynamics of coherent structures, providing new insights into the mechanisms that drive the global transport. We also show that an upscaled model informed by LES and based on a fractional derivative captures the evolution of concentration, and the exchange between the cavity and the main channel, providing accurate predictions of mass transport and reproducing the temporal dependence observed at larger scales.

Efficacy of the Combined Use of Bed Sill and Sacrificial Piles to Control Local Scour around Circular Bridge Piers

In this study, the results are presented and discussed from laboratory test campaigns specifically designed to investigate the behavior of sacrificial piles as a countermeasure against local scouring at a circular bridge pier and clear-water conditions with flow intensity slightly below the threshold of sediment motion. Sacrificial piles are assessed on the upstream side of the pier in two transverse and triangular arrangements. Piles can reduce scouring by deflecting the flow and creating a low-velocity wake region behind them. The efficiency of the piles against local scouring depends on the diameter of the piles, the number of piles, and the angle of the wedge. The investigation was aimed at evaluating the effectiveness of the sacrificial piles as a function of different dimensionless groups. It was found that the triangular arrangement of sacrificial piles has better results than the transverse arrangement. The results showed that the triangular sacrificial piles reduced the maximum local scour depth at the pier to 37.2% in the best configuration. Combined countermeasures were tested, which were composed of sacrificial piles and a bed sill downstream of the pier; in the best configuration, the scour depth reduction in front of the pier reached 51.1%. The increased efficiency of the combination of bed sill and sacrificial piles (BSSP) is an advantage that can reduce the risk of pier failure when the duration of the flood is short. This last result shows that a combination of BSSP may be a very effective countermeasure against local scouring at bridge piers. Finally, the coherent turbulent flow structure around the best combination of BSSP was investigated, and its effect on the bed scouring pattern was studied. A 3D analysis of the bursting process was used. Turbulence characteristics, as well as the occurrence and transition probabilities of bursting events, were calculated. The obtained results confirmed the quite effective effect of the combination of these two countermeasures in reducing the scour depth.

Direct numerical simulation of a wall source dispersion in a turbulent channel flow

A comprehensive direct numerical simulation (DNS) is conducted to analyze the turbulent mixing process of a surface source dispersion emitting from the bottom wall in a turbulent channel flow. A comparative analysis of ten test cases has been conducted to examine the turbulent mixing process and its relationship with the source strength. The study involves an examination of the dispersion process in both physical and spectral spaces, which includes an evaluation of statistical moments of the concentration field, pre-multiplied spatial spectra of the velocity and concentration fields, as well as an analysis of the coherent turbulent structures in the flow. Upon normalization by friction concentration, the DNS analysis reveals that the first- and second-order statistical moments of the concentration field remain unaffected by variations in the source strength. Conversely, the relative intensity of concentration fluctuations and the thickness of the concentration boundary layer exhibit significant susceptibility to variations in the source strength. It has been observed that the dispersion process within the channel is significantly influenced by the small-scale eddies in cases where the source strength is not particularly intense. When confronted with high source strength scenarios, analysis of the simulation results establishes that both small- and large-scale eddies are crucial players in the dispersion phenomena, as evinced by a comparison of the pre-multiplied spectra of concentration and velocity fields. Furthermore, hairpin vortices have been recognized as the primary source of turbulence in the dispersal of substances from regions of high shear to the center of the channel.

A New Pathway for Tornadogenesis Exposed by Numerical Simulations of Supercells in Turbulent Environments

Abstract A simulation of a supercell storm produced for a prior study on tornado predictability is reanalyzed for the purpose of examining the fine-scale details of tornadogenesis. It is found that the formation of a tornado-like vortex in the simulation differs from how such vortices have been understood to form in previous numerical simulations. The main difference between the present simulation and past ones is the inclusion of a turbulent boundary layer in the storm’s environment in the present case, whereas prior simulations have used a laminar boundary layer. The turbulent environment contains significant near-surface vertical vorticity (ζ > 0.03 s−1 at z = 7.5 m), organized in the form of longitudinal streaks aligned with the southerly ground-relative winds. The ζ streaks are associated with corrugations in the vertical plane in the predominantly horizontal, westward-pointing environmental vortex lines; the vortex-line corrugations are produced by the vertical drafts associated with coherent turbulent structures aligned with the aforementioned southerly ground-relative winds (longitudinal coherent structures in the surface layer such as these are well known to the boundary layer and turbulence communities). The ζ streaks serve as focal points for tornadogenesis, and may actually facilitate tornadogenesis, given how near-surface ζ in the environment can rapidly amplify when subjected to the strong, persistent convergence beneath a supercell updraft. Significance Statement In high-resolution computer simulations of supercell storms that include a more realistic, turbulent environment, the means by which tornado-like vortices form differs from the mechanism identified in prior simulations using a less realistic, laminar environment. One possibility is that prior simulations develop intense vortices for the wrong reasons. Another possibility could be that tornadoes form in a wide range of ways in the real atmosphere, even within supercell storms that appear to be similar, and increasingly realistic computer simulations are finally now capturing that diversity.

Buoyancy effects on very-large-scale motions and amplitude modulation in convective atmospheric surface layers

A large amount of high-Reynolds-number experimental data was obtained from long-term observations of the atmospheric surface layer in the Qingtu Lake Observation Array under different thermal stability conditions. To reveal the thermal effects, the variations of the scales of turbulent coherent structures, energy distribution and amplitude modulation as a function of thermal stability are investigated. Results show that in the convective atmospheric surface layer, the superposition of buoyancy and shear drives leads to a large enhancement of turbulent kinetic energy, an increase in the scales of coherent structures, and in particular, a more significant enhancement of the very-large-scale motions kinetic energy. This results in an increased energy proportion of very-large-scale motions, and thus enhancing the amplitude modulation effect on small scales. Moreover, the variations in the coherent structure scale, energy proportion of the very-large-scale motions, amplitude modulation coefficient and the corresponding modulating and carrier signal scales with the thermal stability parameter can all be characterized by the same parametric equation form. This study contributes to a better understanding of turbulence dynamics driven by both buoyancy and shear.

Impact force of non-Boussinesq gravity-current head against a circular cylinder above the seabed

We study the impact force of the gravity-current head on a circular cylinder above the seabed for the non-Boussinesq effect using a large-eddy simulation model well validated against laboratory data. Most existing studies of the gravity current were for the slight density difference between the current and the ambient fluid, assuming the density is constant using the Boussinesq approximation. The aim of the present study, however, is about gravity current of significant density difference. We detail the coherent turbulence structure of the gravity current and the transient oscillation of the impact force on the cylinder over a range of density differences varying from (ρs−ρa)/ρa=0.02 to 2 for the elevation relative to the diameter of the cylinder E/D=0.05, 0.4 and 1.0. The velocity and vorticity field in the current show a gravity-current head separating from the primary current driven by excess pressure at the back, with a structural similarity almost independent of the density difference. The gravity-current head’s frontal and wake velocities are determined from the numerical simulations and related to the buoyancy speed. The impact force on the cylinder rises to a peak immediately after the arrival of the gravity-current head. We propose a set of new drag and lift coefficients using a gravity-driven pressure as a reference for the impact force of the non-Boussinesq current. The peak impact force, made dimensionless by the gravity-driven pressure, is found to be relatively independent of the density difference, regardless of how large the current density may be.

Dynamic Analysis of Tip Leakage Phenomena in Axial Flow Pumps Using a Square-Cavity Jet Model

In the field of pump impeller studies, tip leakage flow (TLF) and the resultant tip leakage vortex (TLV) significantly influence hydraulic efficiency, cavitation, and noise generation. This paper builds a novel square-cavity jet model combined with Large Eddy Simulation (LES) technology to obtain precise the dynamic properties of the TLV, significantly simplifying the computational resources required for numerical simulations. The novel square-cavity jet model simplifies a single blade channel to a square-cavity, and then adds a longitudinal slit on the top wall of the square-cavity. The analysis of both instantaneous and time-averaged flow fields indicates that the interaction between the main flow and the jet is the primary source of TLV generation. This study successfully captures the formation process of the TLV and accurately reveals its turbulent coherent structures. The evolution of the TLV is divided into three main parts: the first part is the jet slot, predominantly characterized by negative vorticity flow. The second part is the TLV formation, which is mainly composed of significant negative streamwise vortices. The third part is the development of the TLV, where positive and negative vorticities begin to interact, resulting in a more complex overall structure. The entire evolution of the TLV phenomenon starts with a concentrated negative vortex, which, after breakdown, develops at a certain angle to the slot and continuously advances towards the sidewall, ultimately resulting in the formation of a large-scale intermingled group of small-scale positive and negative vortices. This research not only provides a new physical model for investigating the tip leakage phenomenon in axial flow pumps but also offers a powerful tool and methodology for future studies in similar complex flow domains.

The modulation of coherent structures by the near-wall motions of particles

Particle–wall interaction generates strong particle near-wall motion, including collision bounce and impact splashing. To distinguish the effect of particles and particle near-wall motions on the turbulent coherent structure, this study carried out three different cases of sand-laden two-phase flow measurements: a uniform sand release at the top, local-laying sand bed and global-laying sand bed (Liu et al., J. Fluid Mech., vol. 943, 2022, A8). Based on large field of view particle image velocimetry/particle tracking velocimetry measurements, we obtained the velocity field of a two-dimensional gas–solid two-phase dilute faction flow $(\varPhi _{v} \sim O(10^{-4}))$ with a friction Reynolds number $R e_{\tau }$ of 3950. Results indicate that particles weaken the high- and low-velocity iso-momentum zones and hairpin vortices, resulting in the increased length scale of the coherent structure. However, the collision bounce and impact splashing break up the inner iso-momentum zone and hairpin vortices while enhancing them in the outer region, thus reducing the structure scale. In addition, the upward-moving particles increase the large-scale structure inclination angle, while the downward-moving particles decrease it. The linear coherence spectrum analysis suggests that the particles themselves do not change the structural self-similarity, but their saltation motions disrupt the similarity of the near-wall structure, making the inclination angle decrease with the scale, and the generated ascending particles reduce the aspect ratio of the streamwise to wall-normal direction in the outer region.

The Chimera Revisited: Wall- and Magnetically-Bounded Turbulent Flows

This review is a first attempt at bringing together various concepts from research on wall- and magnetically-bounded turbulent flows. Brief reviews of both fields are provided: The main similarities identified are coherent (turbulent) structures, flow generation, and transport barriers. Examples are provided and discussed.

Experimental study of coherent structures and bed shear stress on a smooth bed under a solitary wave

The coherent structures induced by a solitary wave on a smooth bed were studied in a wave flume. The solitary wave was near the point of incipient breaking and the wave Reynolds number was in the intermittent turbulent regime. The velocity field inside the bottom boundary layer was measured using a Particle Image Velocimetry (PIV) system. The bed shear stress was obtained from the measured velocities in the viscous sublayer by employing the Newton's law of viscosity. It was found that the generation and maintenance of turbulent coherent structures in a solitary wave boundary layer was like that observed in oscillatory flow boundary layers at similar wave Reynolds numbers, a process that was first described as the bursting phenomenon in the wall boundary layer of steady flows. The effect of coherent structures on the instantaneous velocity profile and bed shear stress was examined. Ejections and sweeps were found to play a crucial role in producing large bed shear stress fluctuations. Quadrant analysis and probability distribution plots of the bed shear stress showed that sweep (Q4) events were responsible for producing instantaneous bed shear stress values as high as four to five times the standard deviation above the mean, while ejection (Q2) events produced bed shear stress values mostly below the mean. Both types of events were more numerous than outward interaction (Q1) and inward interaction (Q3). It was also found that a high proportion of the measured velocity profiles that did not follow the law of the wall were associated with large vortices in the boundary layer. It was shown that the instantaneous velocities did not conform to a law-of-the-wall velocity profile unless the measured velocity field was averaged over an area large compared to the vortices.

Influence of turbulent coherent structures on the performance and wake of a wind turbine

Wind energy in the atmospheric boundary layer serves as the primary source for energy absorption and structural load on wind turbines. However, the impact of turbulent coherent structures on the aerodynamic performance and wake characteristics of wind turbines has not been comprehensively evaluated. In this study, the proper orthogonal decomposition (POD) method is employed to assess the influence of turbulent coherent structures of varying scales on the aerodynamic performance and wake characteristics of wind turbines in the neutral atmospheric boundary layer. The results show that turbulent coherent structures are the main factor that determines the wind velocity fluctuation, aerodynamic performance and wake characteristics of wind turbine in the atmospheric boundary layer. When considering the 13th or lower order POD mode, the wind velocity fluctuation increases with the increase of energy content (more POD modes) of the turbulent coherent structures. When considering the first 19 POD modes, the dynamic loads and power of wind turbine fluctuate with high frequencies, the thrust fluctuates in an amplitude range between 2.4 % and 13.9 % around the mean value, and the power fluctuates from 4.5 % to 28.6 % of the mean value. When considering the first 40 POD modes, the average power generation of the wind turbine increases by 26 % compared to the case with no turbulent structures considered. The study of turbine wake shows that turbulent coherent structures can expand the wind turbine wake approximately to a width of 2.5D and a height of 3D (D is the diameter of the wind turbine), offset the wake approximately to 2D, and move forward the position of the wake vortex beginning to dissipation approximately to 7D behind the wind turbine. In addition, turbulent coherent structures can accelerate the wake velocity recovery by increasing the momentum exchange between the atmospheric boundary layer and wind turbine wake.

Experimental investigation about the aerodynamic noise caused by hypersonic turbulence

In this study, experimental investigations pertaining to aerodynamic noise are conducted in a hypersonic quiet wind tunnel, which possesses the unique capability of operating in both quiet and conventional modes, achieved by controlling the bleed slot to regulate the boundary layer flow. The measurement of aerodynamic noise and analysis of coherent structures of hypersonic turbulence are accomplished employing pressure probes and the nano-tracer-based planar laser scattering technique. The characterization of aerodynamic noise within the freestream unveils that noise generated by turbulent structures on the nozzle wall is the primary factor affecting the turbulence intensity of the freestream flow. Regarding the acoustic radiation noise generated by the turbulence on the nozzle wall, a peak frequency ranging from 54 to 57 kHz is observed in its spectrum within a unit Reynolds number range of 0.3×107–1.1×107 m−1. Under different operation modes, the freestream pitot-pressure data feature the similar spectral slope of f−2.2 and f−12 before and after the aforementioned peak. Through wavelet analysis, scales of the coherent structure in the turbulent boundary layer are examined, revealing that the maximum energy streamwise scale of the coherent structure is approximately 2.47 mm under the unit Reynolds number of 0.7×107 m−1. The product of this scale with the characteristic frequency of 56 kHz obtained from fluctuating pressure measurements is close to the local sound speed of about 140 m/s.

Non-hydrostatic simulation of two-layer density stratified flows over a sill

The topography between the ocean and marginal sea basins is often complex with narrow channels and shallow sills. Two-layer exchange flows are usually observed in the channels with internal waves and small-scale instabilities at the interface. A high-resolution, non-hydrostatic model (NHWAVE) is used to simulate the interfacial mixing processes. Excellent agreement with the experiment data is obtained, especially for the simulation of vertical velocity and internal wave dispersion. In the mixing processes, turbulent coherent structures vary from spanwise uniform vortices to full 3D vortices with random shapes of spanwise and streamwise vortices. These vortices correspond to the surface velocity divergence with wave-like perturbations and the crests of the vortices associate with the down zero crossings of the divergence. The non-hydrostatic effect is also analyzed by an aspect ratio. The results show that the most obvious non-hydrostatic areas are located at the lee side of the sill, ranging from x = 0–1 m. And the effect of non-hydrostatic is more significant with larger values of the vertical velocity gradient.

PROPELLER INFLUENCE ON SURFACE PRESSURE FLUCTUATION UNDER BOUNDARY LAYER INGESTION CONDITION

Research on propeller noise under boundary layer ingestion condition is essential for designing new aircraft concepts. Measurements of pressure spectra on a flat plate below a propeller give important information about the noise coming from the blades. These also are important for designing surface noise mitigation structures, such as acoustic black holes, meta-materials, and liners. This paper investigates the effects of a propeller on the flat plate surface pressure under boundary layer ingestion, using bottom of the wind tunnel channel as flat plate. Tripping devices on the wind tunnel bottom wall control the boundary layer thicknesses. The surface pressure is measured using an instrumented plate positioned underneath the propeller. The plate has 87 pressure taps to measure the unsteady and steady pressure using the remote microphone probe (RMP) technique. The results show the influence of the propeller inflow, the advance ratio, and the number of blades on the different regions underneath the propeller, and on coherent turbulent structures of the incoming flow. The boundary layer thicknesses affects the static pressure upstream of the propeller - creating high pressure regions. The propeller affects the turbulent structures downstream of the blades, increasing its correlation length.

Numerical investigation of communicating turbulent boundary layers through porous media

The reduction of trailing edge noise can be achieved through the use of porous trailing edges. It has been found that the unsteady interaction between the suction and pressure side of the airfoil is of relevance for noise reduction. Therefore, this study aims at investigating such ow physics. The proposed numerical setup consists of two vertically stacked temporally developing channel ows (4 ˇ (2 + t) ˇ 4, where and t are the half-charnnel height and the porous medium thickness respectively). The ow eld is solved using the lattice Boltzmann method through PowerFLOW. The two channel ows communicate through a fully resolved porous medium. In this work, the porous medium is a 75% porous Schwarz' P triply periodic minimal surface. Results show span-wise coherent turbulent structures developing on both sides of the porous medium. This is consistent with previous work on turbulent channels over porous media which enable ow in the stream-wise direction. However, the communication between both turbulent boundary layers results in these span-wise structures to be out-of-phase on either side of the permeable material. This phenomenon occurs at a frequency range where the noise reduction of porous trailing edges is the highest when compared to previous studies. A shift in turbulent kinetic energy carrying scales, with more energy being present in larger scales, is also observed. It can be concluded that this would increase the likelihood of interference occurring in the various locations of hydrodynamic pressure uctuation scattering and as such increase the noise reduction of porous trailing edges.

Effects of different working parameters on the turbulent drag reduction regulated by annular plasma synthetic jet actuators

Reduction of turbulent drag force is one of important works in the design of airplanes and hypersonic aircrafts. The annular plasma synthetic jet (APSJ) has become an interesting and popular flow control method in reducing the drag of turbulent boundary layers. In this paper, a comprehensive experimental study is carried out on the turbulent drag reduction regulated by an array of annular plasma synthetic jet actuators. The effects of the operating parameters such as the actuation voltage, the pulse frequency and the incoming wind speed on the drag reduction rate are studied and discussed in detail. The performances of the plasma actuator array are evaluated and summarized at multiple working conditions. Finally, the coherent structures of turbulence and the hairpin-like vortices are depicted and discussed. The results show that the optimal drag reduction rate is achieved, as the actuation voltage, pulse frequency and incoming wind speed are Vpp = 7 kV, fp= 50 Hz, and U∞ = 7 m/s, respectively. A resonant coupling phenomenon is observed when the pulse frequency of the actuators approaches the characteristic frequency of the coherent structure. The upward sweep flow induced by APSJ actuators may lead to a reduction of the turbulent drag force, but the downward wash flow leads to an increase in the drag. Present study could provide solid experimental data and a helpful guidance for the drag reduction of an airplane.

Effect of neutral interactions on parallel transport and blob dynamics in gyrokinetic scrape-off layer simulations

The effect of neutral interactions on scrape-off layer (SOL) turbulence is investigated in a continuum gyrokinetic code that has been coupled to a continuum kinetic model of neutral transport. This extends the work of a previous paper [T. N. Bernard et al., Phys. Plasmas 9, 052501 (2022)], which compared two NSTX SOL simulations in simple helical geometry, one with neutrals and one without. The former included electron-impact ionization, charge exchange, and wall recycling. Here, the case with neutrals is compared to a gyrokinetic-only simulation that includes an effective ionization source to separate the effect of sourcing from charge exchange collisions. It is observed that sourcing accounts for many features of the simulated SOL with neutrals, including density and temperature magnitudes and reduced normalized density fluctuations, but differences persist. In particular, a flatter density profile results due to changes in parallel transport when neutral collisions are included, illustrating the importance of neutral drag on global plasma properties. An analysis of coherent turbulent structures, or blobs, in these simulations demonstrates the case with neutrals has slower and larger blobs. A series of seeded blob simulations corroborates the blob velocity observation. In general, the blob motion does not contribute significantly to radial transport in these simulations.