Turbulent Structures Research Articles

Response of turbulent energy spectrum and flow structures when vortical motion of a certain scale is suppressed by artificial forcing

Response of turbulent energy spectrum and flow structures when vortical motion of a certain scale is suppressed by artificial forcing

Elucidating the boundary layer turbulence dissipation rate using high-resolution measurements from a radar wind profiler network over the Tibetan Plateau

Abstract. The planetary boundary layer (PBL) over the Tibetan Plateau (TP) exerts a significant influence on regional and global climate, while its vertical structures of turbulence and evolution features remain poorly understood, largely due to the scarcity of observations. This study examines the vertical profile of and daytime variation in the turbulence dissipation rate (ε) in the PBL and free troposphere over the TP using the high-resolution (6 min and 120 m) measurements from a radar wind profiler (RWP) network, combined with hourly data from ERA5 during the period from 1 September 2022 to 31 October 2023. Observational analyses show that the magnitude of ε below 3 km under all-sky conditions exhibits a large spatial discrepancy over the six RWP stations over the TP. Particularly, the values of ε at Minfeng and Jiuquan over the northern TP and at Dingri (alternately Tingri) over the southern TP are roughly an order of magnitude greater than those at Lijiang, Ganzi (alternately Garzê), and Hongyuan over the eastern TP. This could be partially attributed to the difference in land cover across the six RWP stations. In terms of the diurnal variation, ε rapidly intensifies from 09:00 local standard time (LST) to 14:00 LST and then gradually levels off in the late afternoon. Under clear-sky conditions, both ε and the planetary boundary layer height (zi) are greater compared with cloudy-sky conditions, which could be due to the cooling effect of clouds, which reduces the solar irradiation reaching the surface. In the lower PBL (0.3 ≤ z/zi ≤ 0.5), where z is the height above ground level, the dominant influential factor in the development of turbulence is the surface–air temperature difference (Ts−Ta). By comparison, in the upper PBL (0.6 ≤ z/zi ≤ 1.0), both Ts−Ta and vertical wind shear (VWS) affect the development of turbulence. Above the PBL (1.0 < z/zi ≤ 2.0), the shear production resulting from VWS dominates the variation in turbulence. Under cloudy-sky conditions, the reduced Ts−Ta and weakened surface sensible heat flux tend to inhibit the turbulent motion in the PBL. On the other hand, the strong VWS induced by clouds enhances the turbulence above the PBL. The findings obtained here underscore the importance of the RWP network in revealing the fine-scale structures of the PBL over the TP and gaining new insight into the PBL evolution.

Experimental study on the interaction between particles and coherent structures in dilute liquid–solid two-phase turbulent boundary layer

In order to reveal the modulation of particles on turbulent coherent structures and particle concentration variations induced by turbulence, particle image velocimetry experiments were performed in a dilute liquid–solid two-phase turbulent boundary layer. Large-scale coherent structures, including streaks and hairpin vortices, were extracted in the liquid–solid two-phase turbulent boundary layer by using multiple methods, such as proper orthogonal decomposition and finite-time Lyapunov exponent. The results showed that the addition of particles thickened the buffer layer of the average velocity profile, weakened the streamwise turbulence intensity at y+ > 10, strengthened the wall-normal turbulence in the near wall region (y+ < 52), and weakened the Reynolds stress in the logarithmic law region. A more thorough study of the turbulent structure revealed that the particles caused the mean spacing of the near-wall streak structure to increase, hairpin vortex migration to slow down, the uplift of hairpin vortex to be suppressed, and the overall size of the hairpin vortex packet to be compressed. This demonstrated that particles could modulate turbulence by inhibiting the development of hairpin vortices.

Turbulent boundary layer structures downstream of a square cylinder

Particle image velocimetry is used to investigate the coherent structures within a turbulent boundary layer (TBL) perturbed by a slender square cylinder. A moderate Reynolds number of Reτ=3691 is considered. The cylinder is immersed in different vertical regions of the TBL, from the viscous sublayer to the logarithmic layer, depending on the value of the cylinder–wall gap ratio within the range of G/D=0.0−2.0. Here, G is the distance between the wall and the lower surface of the cylinder, and D is the cross-section size of the square cylinder. The results reveal a substantial influence of the TBL statistics and structures. Utilizing velocity correlation and proper orthogonal decomposition (POD), we highlight the intriguing influence of large-/very large-scale turbulent structures induced by the cylinder. When the cylinder is positioned on or in close proximity to the wall, periodic shedding of vortices is not apparent. In the presence of a cylinder, the large-/very large-scale turbulent structures are reduced in size compared with those in the undisturbed TBL, owing to the interaction between the flow around the cylinder and the boundary layer turbulence. The lower modes of the POD in the disturbed TBL are similar to those of the undisturbed flow field. However, the large-/very large-scale turbulent structures are increasingly affected by the strengthening of the shedding vortices with increasing G/D when the cylinder is located in the logarithmic region. As a result, the Kármán vortex street becomes the dominant structure associated with the higher POD modes at sufficiently high cylinder–wall gap ratios.

Embedded Coherent Structures from Magnetohydrodynamics to Sub-ion Scales in Turbulent Solar Wind at 0.17 au

We study intermittent coherent structures in solar wind turbulence from MHD to kinetic plasma scales using Parker Solar Probe data during its first perihelion (at 0.17 au) in the highly Alfvénic slow solar wind. We detect coherent structures using Morlet wavelets. For the first time, we apply a multiscale analysis in physical space. At MHD scales within the inertial range, times scales τ ∈ (1, 102) s, we find (i) current sheets including switchback boundaries and (ii) Alfvén vortices. Within these events are embedded structures at smaller scales: typically Alfvén vortices at ion scales, τ ∈ (0.08, 1) s, and compressible vortices at sub-ion scales, τ ∈ 8(10−3, 10−2) s. The number of coherent structures grows toward smaller scales: we observe ∼200 events during a 5 hr time interval at MHD scales, ∼103 at ion scales, and ∼104 at sub-ion scales. In general, there are multiple structures of ion and sub-ion scales embedded within one MHD structure. There are also examples of ion and sub-ion scale structures outside MHD structures. To quantify the relative importance of different types of structures, we do a statistical comparison of the observed structures with the expectations of models of the current sheets and vortices. The results show the dominance of Alfvén vortices at all scales in contrast to the widespread view of the dominance of current sheets. This means that Alfvén vortices are important building blocks of Alfvénic solar wind turbulence.

Turbulent flow structures and Reynolds stress anisotropy in an asymmetric sinuous mobile channel

Turbulent flow structures and Reynolds stress anisotropy in an asymmetric sinuous mobile channel

Development of a multiple laser-sheet imaging technique for the analysis of three-dimensional turbulent explosion flame structures

The present work describes the development of a multiple laser-sheet imaging technique for the study of turbulent, premixed flame surface structures at high Karlovitz stretch factor values. Experiments were conducted using CH4 air mixtures at 365 K and 0.5 MPa and up to a root mean square turbulence velocity of 1.5 m/s. A high-speed Nd:YAG laser capable of pulsing up to 60 kHz in conjunction with a high-speed camera and a rotating mirror was used to reconstruct time-resolved three-dimensional turbulent flames. This has, for the first time, enabled the direct measurement of wrinkled flame surface areas, along with reaction progress variable and flame brush thickness. These are important parameters for the characterization of turbulent burning rates and provide more insight into the dynamic nature of the flames and their structures. In addition, the current data aid toward direct comparison with results from combustion simulation studies.

Turbulence Embedded Into the Ionosphere by Electromagnetic Waves

AbstractWhen charged particles are accelerated from Earth's magnetosphere and precipitate into the atmosphere, their impact with neutral gas creates the aurora. Structured electric fields drive the acceleration processes but they are also passed down to the ionosphere, meaning that turbulence can in part be embedded into the ionosphere rather than emerge through instability processes locally. Applying a point‐cloud analysis technique adapted from observational cosmology, we show how observed turbulence in the ionosphere matches electrical current signatures in the pulsating aurora in a series of conjunctions between space‐ and ground‐based instruments. We propose that the temporal spectrum of pulsations in the pulsating aurora is the driver of a clearly observed energy injection into the ionosphere's unstable bottomside. Precipitating electrons produce electric fields through charge deposition, and we observe wave characteristics that are present in this pattern. Next, the relative electron‐ion drifts excite the Farley‐Buneman instability, the distribution of whose waves are organized according to the local electric field. It is the temporal characteristics of chorus wave interactions in the magnetosphere that is imparted, via precipitating electrons, to the pulsating aurora, and so we propose that chorus wave interactions are capable of embedding turbulent structure into the ionosphere. This structure (now pressure gradients) dissipate energy in the E‐region through turbulent processes, observed by the icebear coherent scatter radar.

Particle Acceleration in Relativistic Alfvénic Turbulence

Strong magnetically dominated Alfvénic turbulence is an efficient engine of nonthermal particle acceleration in a relativistic collisionless plasma. We argue that in the limit of strong magnetization, the type of energy distribution attained by accelerated particles depends on the relative strengths of turbulent fluctuations δ B 0 and the guide field B 0. If δ B 0 ≪ B 0, the particle magnetic moments are conserved, and the acceleration is provided by magnetic curvature drifts. Curvature acceleration energizes particles in the direction parallel to the magnetic field lines, resulting in log-normal tails of particle energy distribution functions. Conversely, if δ B 0 ≳ B 0, interactions of energetic particles with intense turbulent structures can scatter particles, creating a population with large pitch angles. In this case, magnetic mirror effects become important, and turbulent acceleration leads to power-law tails of the energy distribution functions.

Effects of turbulence on the blade root load of the tandem wind turbine in neutral atmospheric boundary layer

Abstract A model characterizing a turbulent wind field, consistent with neutral atmospheric boundary layer attributes, was developed for a field-based experimental wind turbine using stochastic processes combined with shear flow profiles. By incorporating interactions of vertical thermal exchange, terrain-induced surface roughness, and Coriolis forces, the model simulates environmental impacts. A sequence of wind turbines situated within the atmospheric wind field underwent computational analysis leveraging Large Eddy Simulation (LES) alongside an Actuator Line Model (ALM) to examine atmospheric turbulence alterations both preceding and following the turbines. Analysis of the power spectra concerning wind velocity, turbulence intensity, and components of Reynolds stress upwind and downwind demonstrated the turbines’ blade effects on surrounding turbulence patterns. The research demonstrated that turbulence kinetic energy experienced enhancement in downwind turbines as a result of the upwind turbine’s impact, particularly evident in the lower-frequency spectrum. Augmenting the distance between turbines initially aids in the wake’s turbulent structure restoration, with a marked effect on tip vortices that hold greater amounts of high-frequency energy.

Generation of inflow turbulence using an improved synthetic eddy method

The generation of turbulent inflow conditions is a key issue in large eddy simulation (LES) or direct numerical simulation (DNS). In this paper, an improved synthetic eddy method (SEM) is proposed to generate inflow turbulence for LES and DNS. The improvements about SEM focus on the eddy radius and eddy distributions. First, the eddy radius is improved to reduce the nonphysical vortex structure on the wall caused by overestimation of the radius. Second, a sampling method using Gaussian distribution is proposed to improve the distribution of eddies, which accurately captures the randomness of turbulent structure size and is close to the actual flow field. The improved method is applied to the direct numerical simulation of the supersonic turbulent boundary layer at Mach 2.7 and the 24° compression ramp. Results indicate that the predictions yielded by the improved method are in good agreement with both DNS and experimental data. Compared to the original method, the improved method exhibits a more rapid recovery of the friction coefficient and effectively shortens the development distances. The improved SEM has enhanced the efficiency and accuracy of generating inlet turbulence, which can provide inlet turbulence boundary conditions for LES and DNS.

Exact coherent structures in two-dimensional turbulence identified with convolutional autoencoders

Convolutional autoencoders are used to deconstruct the changing dynamics of two-dimensional Kolmogorov flow as $Re$ is increased from weakly chaotic flow at $Re=40$ to a chaotic state dominated by a domain-filling vortex pair at $Re=400$ . ‘Latent Fourier analysis’ (Page et al., Phys. Rev. Fluids6, 2021, p. 034402) reveals a detached class of bursting dynamics at $Re=40$ which merge with the low-dissipation dynamics as $Re$ is increased to $100$ and provides an efficient representation within which to find unstable periodic orbits (UPOs) using recurrent flow analysis. Focusing on initial guesses with energy in higher latent Fourier wavenumbers allows a significant number of high-dissipation-rate UPOs associated with the bursting events to be found for the first time. At $Re=400$ , the UPOs discovered at lower $Re$ move away from the attractor, and an entirely different embedding structure is formed within the network devoid of small-scale vortices. Here latent Fourier projections identify an associated ‘large-scale’ UPO which we believe to be a finite- $Re$ continuation of a solution to the Euler equations.

Three-dimensional generative adversarial networks for turbulent flow estimation from wall measurements

Different types of neural networks have been used to solve the flow sensing problem in turbulent flows, namely to estimate velocity in wall-parallel planes from wall measurements. Generative adversarial networks (GANs) are among the most promising methodologies, due to their more accurate estimations and better perceptual quality. This work tackles this flow sensing problem in the vicinity of the wall, addressing for the first time the reconstruction of the entire three-dimensional (3-D) field with a single network, i.e. a 3-D GAN. With this methodology, a single training and prediction process overcomes the limitation presented by the former approaches based on the independent estimation of wall-parallel planes. The network is capable of estimating the 3-D flow field with a level of error at each wall-normal distance comparable to that reported from wall-parallel plane estimations and at a lower training cost in terms of computational resources. The direct full 3-D reconstruction also unveils a direct interpretation in terms of coherent structures. It is shown that the accuracy of the network depends directly on the wall footprint of each individual turbulent structure. It is observed that wall-attached structures are predicted more accurately than wall-detached ones, especially at larger distances from the wall. Among wall-attached structures, smaller sweeps are reconstructed better than small ejections, while large ejections are reconstructed better than large sweeps as a consequence of their more intense footprint.

Influence of Time-Varying Freestream Conditions on the Dynamics of Unsteady Boundary-Layer Separation

Unsteady flow separation of a turbulent boundary layer under dynamic pressure gradients is investigated using the Large-Eddy Simulation technique. The unsteadiness is introduced by prescribing an oscillating freestream vertical-velocity profile at the top boundary of the domain. Although previous studies, including Ambrogi et al. (Journal of Fluid Mechanics, Vol. 945, Aug. 2022, p. A10) and Ambrogi et al. (Journal of Fluid Mechanics, Vol. 972, Oct. 2023, p. A36), focused on the kinematics of the flow and the effects of the oscillation frequency on flow separation, the goal of this paper is to analyze the effects of three time-varying freestream-forcing profiles while the oscillation frequency is kept the same for all Cases. Whereas in Case A the freestream-velocity profile changes from suction–blowing to blowing–suction in a complete cycle, Cases B and C are both suction–blowing only and the strength of the adverse pressure gradient is modulated in time. Moreover, the boundary layer in Case B never approaches a zero-pressure gradient condition. A closed separation bubble is formed for all Cases; however, its dimensions change depending on the far-field forcing. The time evolution of turbulent kinetic energy (TKE) reveals an advection mechanism of turbulent structures out of the domain for all Cases. Whereas in Case A and C the high-TKE region, generated in the separated shear layer, is washed out of the domain as a rigid body, in Case B the separation bubble remains present and the advection mechanism of TKE is characterized by a breathing pattern.

Aerodynamic characteristics of the retro propulsion landing burn of vertically landing launchers

In the frame of the European funded H2020 project RETALT (retro propulsion-assisted landing technologies), the unsteady aerodynamics of vertically descending and landing launchers have been investigated. In this paper, experimental data of the landing burn tested in the Vertical Free-Jet Facility Cologne at DLR in Cologne are presented. The landing burn was simulated with a cold gas jet of pressurized air opposing the wind tunnel free stream. Tests with several jet conditions were compared to results without active jet. Proper orthogonal decomposition of schlieren recordings and spectral analyses of their time histories are performed and are compared to frequencies in pressure measurements. Dominant frequencies were found, which are strongest at Mach 0.8. Especially, a Strouhal number of 0.2 was found to be most dominant. The intensity of the dominant frequencies can be lowered if the engine is active. The normalized root mean square pressure fluctuations are between 0.1 and 0.3 during the landing maneuver. Additionally, the steady flow features scale well with the ambient pressure ratio and the momentum flux ratio. The unsteady flow field dynamics of the subsonic retro propulsion flow field can likely be linked to large-scale turbulent structures in the supersonic jet, triggering large-scale pressure fluctuations and altering the overall flow field.

Turbulence modulation by charged inertial particles in channel flow

Large amounts of small inertial particles embedded in a turbulent flow are known to modify the turbulent statistics and structures, a phenomenon referred to as turbulence modulation. While particle electrification is ubiquitous in particle-laden turbulence and significantly alters particle behaviour, the effects of inter-particle electrostatic forces on turbulence modulation and the underlying physical mechanisms remain unclear. To fill this gap, we perform a series of point-particle direct numerical simulations of turbulent channel flows at a friction Reynolds number of approximately 540, laden with uncharged and charged bidisperse particles. The results demonstrate that, compared to flows laden with uncharged particles, the presence of inter-particle electrostatic forces leads to substantial changes in both turbulent intensities and structures. In particular, the inner-scaled mean streamwise fluid velocity is found to shift towards lower values, indicating a noticeable increase in fluid friction velocity. Turbulent intensities appear to be further suppressed through facilitating the particles to extract momentum from the fluid phase and increasing extra turbulent kinetic dissipation by particles. Importantly, the overall drag is enhanced by indirectly strengthening the contribution of particle stress, even though the contribution of the total fluid stress is decreased. On the other hand, the magnitude of the large-scale motions is weakened by simultaneously reducing turbulent production and increasing particle feedback around the scales of the large-scale motions. Meanwhile, the average streaky fluid structures in the streamwise–spanwise planes and inclined fluid structures in the streamwise–wall-normal planes become expanded and flattened, respectively.

Kármán Street-Boundary Layer Interactions In Turbulent Flow

Understanding the intricate interplay between turbulent boundary layers and wake structures caused by obstacles in the flow is key to assessing the local shear stress fields as well as the drag. These types of flows can be relevant to natural phenomena including sediment and nutrient transport. While existing literature predominantly explores the flow downstream of cylinders mounted perpendicular to bounding surfaces or downstream of cylinders positioned outside the boundary layer region, this work experimentally studies the turbulent boundary layer passing over a cylindrical obstacle positioned near and parallel to the surface of an underlying solid wall. Planar particle image velocimetry (PIV) was performed in an open recirculating water channel for three test cases at a friction Reynolds number Reτ = 1320. Key parameters investigated include the cylinder diameter relative to the boundary layer thickness (D/δ = 0.09 and 0.01) and the gap between the cylinder and the wall (G/δ = 0.09), where the range of D/δ used is smaller than most examples found in the literature. Instantaneous and averaged velocity and vorticity fields are examined based on the PIV measurements. The results in Fig. 1 show how von Kármán vortex streets and turbulent boundary layer structures interact to form unique flow features. The larger cylinder generates a clear separation downstream of its axis as well as a wake that remains coherent over a distance δ downstream. The smaller cylinder with D/δ = 0.01 generates a much narrower and weaker wake, but still shows a clear influence on time-averaged velocity profile. Future experiments will evaluate flow conditions further downstream as well as additional values of G/δ for the same D/δ values.

Volumetric Reconstruction Of Stagnated Wake Structures In Adverse Pressure Gradient

To gather insights about turbulent structures of a wake in adverse pressure gradient, a wind tunnel model is investigated with PIV methods. A flat plate of L = 1.058 m at zero angle of attack is placed in the MUB wind tunnel at TU Braunschweig. The wake convects downstream into a carefully designed diffuser, where it undergoes a pressure increase of ∆c p = +0.6. Wakes, with the common example of a cylinder shedding van-Karman vortices, are generally prone to instability. The wake instability of the sharp-edged plate is attenuated by the adverse pressure gradient, manifesting in large coherent structures that carry up to 37 % of the turbulent kinetic energy. While the primary shedding motion is effectively a quasi-2D motion, a secondary mode of is also observed: Over the span S = 1.3 m of the wake, it organizes in three distinct cells. Their motion is similar to "owl-eyes", that are known to appear on separated flow over smooth 2D bodies. The presence of these modes reduces the average flow velocity in the wake centre up to stagnation in the mean sense. Because stagnation and shedding does occur far away from any solid surfaces, this test case offers excellent observability of the phenomenon with PIV methods. RANS results of this setup cannot predict the flow stagnation accurately. Because the shedding frequency is low compared to freestream velocity, scale and time resolving IDDES results need to run many convective time units until the observed instability develops. While the original goal of observing small scale turbulent wake structures is compromised by the large shedding modes, the observation gives an insight on wake-bursting phenomena that occur on multi-element high-lift devices. The mode shedding and consequent wake stagnation is sensitive to the Reynolds-number (Re = [1.6 ... 3.2] · 10 6 ) as well as to the initial boundary layer thickness of the flat plate.

Multi-Scale Flow Estimation Within Stereoscopic PIV Processing Based On Deep Learning

Classical methods of particle-image velocimetry (PIV) based on state-of-the-art cross-correlation algorithms constitute a spatial filtering operation that leads to a loss of spatial resolution of the estimated velocity fields. This is associated with numerous limitations, e.g., when it comes to aligning high-fidelity numerical simulations, since the detection of small turbulent scales is often precluded. An optical flow neural network adapted for PIV processing, called RAFT-PIV, can overcome the drawback of reduced spatial resolution while maintaining the reliability and robustness of the classical approach. The present work aims to extend the application of the aforementioned deep learning based method by integrating it within a stereoscopic PIV processing workflow. Experimental data from an optically accessible single-cylinder combustion engine the flow field of which is characterized by complex, three-dimensional flow structures and a broad range turbulent scales are processed using RAFT-PIV. The 2D-3C vector fields obtained from stereoscopic reconstruction are compared against a classical cross-correlation based PIV evaluation, providing insight into instantaneous quantities and turbulent structures. A qualitative comparison of the instantaneous in-plane and out-of-plane velocity components reveals that the reconstructed flow obtained from the deep learning based estimation includes a wide range of resolved turbulent length scales. It outperforms the classical method. Furthermore, the findings of the present study demonstrate that the pixel accuracy has an effect on both physical quantities and turbulence characteristics. In particular, an increase in turbulent kinetic energy is observed across the engine cycles of the investigated engine flow states. Moreover, a higher fraction of three-dimensional isotropic turbulent structures is derived from RAFT-PIV compared to the results from the cross-correlation based technique.

High-Resolution Two-Phase Velocimetry Of Aspherical Particles Using Wavelet-Based Optical Flow Velocimetry (WOFV) Benchmarked With Fully Resolved Direct Numerical Simulations

Several techniques exist to measure the carrier and dispersed phases in multi-phase flows. Particle image velocimetry (PIV) analyzes the carrier phase by cross-correlating particle images, but its resolution is limited by the size of the interrogation window, making it difficult to resolve fine-scale turbulent structures. Particle tracking algorithms capture translational motion of dispersed particles well, but struggle with rotational motion, especially for irregular and aspherical particles. Currently, there is no method that universally measures both the carrier phase velocity field and the translational and rotational motions of dispersed particles. This study evaluates wavelet-based optical flow velocimetry (wOFV) for motion estimation in multi-phase flows, focusing on dispersed ellipsoidal particles and their surrounding turbulent carrier flow using synthetic image data. The research proceeds in two phases: first, the rigid motion of ellipses, including translational and rotational components, is analyzed using wOFV-generated dense motion fields. The results highlight the critical role of the regularization weighting parameter λ. Higher λ values improve translational motion estimation, while an optimal λ avoids under-regularization and non-physical structures in rotational motion. wOFV maintains accuracy in combined motion scenarios with optimal λ values. In the second phase, wOFV captures the turbulent carrier flow around an aspherical particle, benchmarked against DNS data and compared to PIV. wOFV outperforms PIV in resolving finer structures near the particle surface and in accurately representing the wake region. Error analyses confirm wOFV’s superior performance, with optimal results within a specific λ range. In conclusion, wOFV is a highly effective tool for the analysis of multi-phase flow dynamics, providing greater resolution and accuracy than PIV, especially in complex scenarios involving aspherical particles.