In this paper it is shown that a modal detuned instability of periodic near-wall streaks originates a large-scale structure in the bulk of the turbulent channel flow. The effect of incoherent turbulent fluctuations is included in the linear operator by means of an eddy viscosity. The base flow is an array of periodic two-dimensional streaks, extracted from numerical simulations in small domains, superposed to the turbulent mean profile. The stability problem for a large number of periodic units is efficiently solved using the block-circulant matrix method proposed by Schmid et al. (Phys. Rev. Fluids, vol. 2, 2017, 113902). For friction Reynolds numbers equal or higher than $590$, it is shown that an unstable branch is present in the eigenspectra. The most unstable eigenmodes display large-scale modulations whose characteristic wavelengths are compatible with the large-scale end of the premultiplied velocity fluctuation spectra reported in previous computational studies. The wall-normal location of the large-wavelength near-wall peak in the spanwise spectrum of the eigenmode exhibits a power-law dependence on the friction Reynolds number, similarly to that found in experiments of pipes and boundary layers. Lastly, the shape of the eigenmode in the streamwise-wall-normal plane is reminiscent of the superstructures reported in the recent experiments of Deshpande et al. (J. Fluid Mech., vol. 969, 2023, A10). Therefore, there is evidence that such large-wavelength instabilities generate large-scale motions in wall-bounded turbulent flows.

An experimental investigation is conducted to study the flow patterns, spectral properties and energy fluxes in thin-layer turbulence with varying system sizes and damping rates. It is found that although a system-size vortex (an indicator of spectral condensation) occurs for small system sizes and does not for large ones, the spectra for different system sizes consistently exhibit a scaling close to $k^{-3}$ in inverse cascade (another indicator of spectral condensation). On the other hand, under a fixed system size larger than the friction-dominated length scale, the energy spectrum in the inverse cascade range changes from $k^{-3}$ to $k^{-5/3}$ as the damping rate increases, suggesting that the friction-dominated length scale may not be a suitable parameter for predicting spectral transition. At lower damping rates and large system sizes, turbulent structures grow larger via inverse cascade, manifesting as long streamers, and the small-scale vortices are suppressed. This suppression leads to a reduction of energy flux at intermediate scales and a change in the spectral shape. The dimensionless Taylor microscale is found to exhibit a monotonic dependence on the damping rate. With the reduction in the damping rate, the Taylor microscale increases to become comparable with the forcing scale, and the spectrum in inverse cascade transits to a steeper scaling, $k^{-3}$, indicating that the dimensionless Taylor microscale may be used as a diagnostic parameter for spectral transition.

Predicting the magnitude of the slip velocity of non-tracer particles with respect to the surrounding fluid is crucial to address both fundamental and practical questions involving dispersed turbulent flows. Here we derive an analytical model to predict the slip velocity of spherical particles in homogeneous isotropic turbulence. We modulate the particle equation of motion according to the inertial filtering framework, and obtain closed-form expressions for the mean slip velocity magnitude as a function of the governing parameters. These are compared against laboratory measurements and direct numerical simulations, demonstrating close agreement for both light and heavy particles, both smaller and larger than the Kolmogorov scales. The predictive value of the model and its implications are discussed, as well as the range of validity of the underlying assumptions.

Tom et al. (J. Fluid Mech., vol. 947, 2022, p. A7) investigated the impact of two-way coupling (2WC) on particle settling velocities in turbulence. For the limited parameter choices explored, it was found that (i) 2WC substantially enhances particle settling compared with the one-way coupled case, even at low mass loading $\varPhi _m$ and (ii) preferential sweeping remains the mechanism responsible for the particles settling faster than the Stokes settling velocity in 2WC flows. However, significant alterations to the flow structure that can occur at higher mass loadings mean that the conclusions from Tom et al. (J. Fluid Mech., vol. 947, 2022, p. A7) may not generalise. Indeed, even under very low mass loadings, the influence of 2WC on particle settling might persist, challenging the conventional assumption. We therefore explore a much broader portion of the parameter space, with simulations covering cases where the impact of 2WC on the global fluid statistics ranges from negligible to strong. We find that, even for $\varPhi _m=7.5\times 10^{-3}$, 2WC can noticeably increase the settling for some choices of the Stokes and Froude numbers. When $\varPhi _m$ is large enough for the global fluid statistics to be strongly affected, we show that preferential sweeping continues to be the mechanism that enhances particle settling rates. Finally, we compare our results with previous numerical and experimental studies. While in some cases there is reasonable agreement, discrepancies exist even between different numerical studies and between different experiments. Future studies must seek to understand this before the discrepancies between numerical and experimental results can be adequately addressed.

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 wake of two tandem square cylinders of identical width (d) is experimentally studied, with a view to understanding the dependence of the flow structure, aerodynamics forces and Strouhal number on the centre-to-centre spacing ratio L/d and Reynolds number Re, where L is the distance between the cylinder centres. Extensive measurements are carried out, using hot-wire, particle imaging velocimetry, laser-induced fluorescence flow visualization, surface-oil-flow visualization and surface pressure scanning techniques, for L/d = 1.0 ~ 5.0 and Re ≡ U∞d/ν = 2.8 × (103 ~ 104), where U∞ is the free-stream velocity and ν is the kinematic viscosity of the fluid. The flow is classified into four regimes, i.e. the extended-body (L/d ≤ 1.5–2.0), reattachment (1.5–2.0 < L/d < 2.7–3.2), co-shedding (L/d ≥ 3.0–3.4) and transition (2.7 ≤ L/d ≤ 3.3) where both reattachment and co-shedding phenomena may take place. The mean drag and fluctuating drag and lift exhibit distinct features for different flow regimes, which is fully consistent with the proposed flow classification. Comparison is made between this flow and the wake of two tandem circular cylinders, which provides valuable insight into the profound effect of the flow separation point and the presence of sharp corners on the flow development and classification.

Despite decades of research, a universal method for prediction of roughness-induced skin friction in a turbulent flow over an arbitrary rough surface is still elusive. The purpose of the present work is to examine two possibilities; first, predicting equivalent sand-grain roughness size $k_s$ based on the roughness height probability density function and power spectrum (PS) leveraging machine learning as a regression tool; and second, extracting information about relevance of different roughness scales to skin-friction drag by interpreting the output of the trained data-driven model. The model is an ensemble neural network (ENN) consisting of 50 deep neural networks. The data for the training of the model are obtained from direct numerical simulations (DNS) of turbulent flow in plane channels over 85 irregular multi-scale roughness samples at friction Reynolds number $Re_\tau =800$. The 85 roughness samples are selected from a repository of 4200 samples, covering a wide parameter space, through an active learning (AL) framework. The selection is made in several iterations, based on the informativeness of samples in the repository, quantified by the variance of ENN predictions. This AL framework aims to maximize the generalizability of the predictions with a certain amount of data. This is examined using three different testing data sets with different types of roughness, including 21 surfaces from the literature. The model yields overall mean error 5 %–10 % on different testing data sets. Subsequently, a data interpretation technique, known as layer-wise relevance propagation, is applied to measure the contributions of different roughness wavelengths to the predicted $k_s$. High-pass filtering is then applied to the roughness PS to exclude the wavenumbers identified as drag-irrelevant. The filtered rough surfaces are investigated using DNS, and it is demonstrated that despite significant impact of filtering on the roughness topographical appearance and statistics, the skin-friction coefficient of the original roughness is preserved successfully.

Proper orthogonal decomposition (POD) is applied to three-dimensional (3-D) velocity fields collected from large-eddy simulations (LES) of a baffled stirred tank. In the LES, the tank operates with a Rushton-type impeller under turbulent conditions (at least in the near-impeller region) and the working fluid exhibits either Newtonian or shear-thinning rheology. The most energetic POD modes are analysed, and a POD reconstruction based on the higher modes is proposed to approximate the fluctuating component of the velocity field. Subsequently, the POD reconstruction is used to identify vortical structures and characterise them in terms of their shape. The structures are identified by considering a frame-invariant formulation of a popular, Eulerian, local-region-type method: the $Q$-criterion. Statistics of shape-related parameters are then investigated to address the morphology of the structures. It is found that: (i) regardless of the working fluid rheology, it seems feasible to decompose the 3-D field into its mean, most energetic periodic and fluctuating components using POD, allowing, for instance, reduced-order modelling of the energetic periodic motions for mixing enhancement purposes, and (ii) vortical structures related to turbulence are mostly tubular. Finding (ii) implies that, as starting point, phenomenological models for the interaction between fluid particles (drops and bubbles) and vortices should consider the latter as cylindrical structures rather than of spherical shape, as classically assumed in these models.

In this paper, we study the aeroacoustic instability which occurs in a deep axisymmetric cavity in a turbulent pipe flow. This phenomenon is the axisymmetric counterpart of the classical whistling of a rectangular deep cavity subject to a grazing flow. The whistling of such axisymmetric cavity originates from the interaction of the coherent fluctuations of the vorticity at the cavity's opening with one of its trapped azimuthal or radial acoustic modes. We focus here on the situation involving the first pure azimuthal mode, which is trapped in the cavity. As a consequence of the rotational symmetry of the configuration, azimuthal modes are actually pairs of degenerate eigenmodes, or almost degenerate in the presence of small asymmetries. Therefore, the aeroacoustic instabilities exhibit more complex mechanisms than in the case of a rectangular deep cavity. In particular, we show that self-sustained spinning modes induce a symmetry breaking of the mean flow and we will elucidate the details of this phenomenon. To that end, simultaneous acoustic and time-resolved stereoscopic particle image velocimetry (PIV) measurements are performed. They reveal that when large-amplitude aeroacoustic waves spin around the cavity, a quasi-steady mean flow starts whirling slowly in the opposite direction to the wave propagation. A linear perturbation analysis around an axisymmetric mean flow confirms the experimental observations: although the incoming pipe flow is not swirling, the hydrodynamic component of the aeroacoustic wave induces such whirling motion of the mean flow because of the forcing from the steady part of the coherent Reynolds stress tensor.

We propose an experimental study on the gravitational settling velocity of dense, sub-Kolmogorov inertial particles under different background turbulent flows. We report phase Doppler particle analyser measurements in a low-speed wind tunnel uniformly seeded with micrometre scale water droplets. Turbulence is generated with three different grids (two consisting of different active-grid protocols while the third is a regular static grid), allowing us to cover a very wide range of turbulence conditions in terms of Taylor-scale-based Reynolds numbers ($Re_\lambda \in [30\unicode{x2013}520]$), Rouse numbers ($Ro \in [0\unicode{x2013}5]$) and volume fractions ($\phi _v \in [0.5\times 10^{-5}\unicode{x2013}2.0\times 10^{-5}]$). We find, in agreement with previous works, that enhancement of the settling velocity occurs at low Rouse number, while hindering of the settling occurs at higher Rouse number for decreasing turbulence energy levels. The wide range of flow parameters explored allowed us to observe that enhancement decreases significantly with the Taylor–Reynolds number and is significantly affected by the volume fraction $\phi _v$. We also studied the effect of large-scale forcing on settling velocity modification. The possibility of changing the inflow conditions by using different grids allowed us to test cases with fixed $Re_\lambda$ and turbulent intensity but with different integral length scale. Finally, we assess the existence of secondary flows in the wind tunnel and their role on particle settling. This is achieved by characterising the settling velocity at two different positions, the centreline and close to the wall, with the same streamwise coordinate.

We use experiments to explore the effect of surfactants on bubble-induced turbulence (BIT) at different scales, considering how the bubbles affect the flow kinetic energy, anisotropy and extreme events. To this end, high-resolution particle shadow velocimetry measurements are carried out in a bubble column in which the flow is generated by bubble swarms rising in water for two different bubble diameters (3 and 4 mm) and moderate gas volume fractions (0.5 %–1.3 %). We use tap water as the base liquid and add 1-Pentanol as an additional surfactant with varying bulk concentration, leading to different bubble shapes and surface boundary conditions. The results reveal that with increasing surfactant concentration, the BIT generated increases in strength, even though bubbles of a given size rise more slowly with surfactants. We also find that the level of anisotropy in the flow is enhanced with increasing surfactant concentration for bubbles of the same size, and that for the same surfactant concentration, smaller bubbles generate stronger anisotropy in the flow. Concerning the intermittency quantified by the normalized probability density functions of the fluid velocity increments, our results indicate that extreme values in the velocity increments become more probable with decreasing surfactant concentration for cases with smaller bubbles and low gas void fraction, while the effect of the surfactant is much weaker for cases with larger bubble and higher void fractions.

This work aims at studying the mechanisms behind the occurrence of extreme dissipation events in a channel flow, identifying nonlinear optimal perturbations as potential precursors of these events. Nonlinear optimal perturbations with respect to a generic turbulent instantaneous snapshot are computed for the first time using a direct-adjoint algorithm in the channel flow at $Re_{\tau }\approx 180$. The resulting initial perturbation displays the upstream tilting characteristic of Orr's mechanism and is positioned along the interfaces between two opposite-sign velocity streaks of the pre-existing turbulent field. Such a perturbation induces a sudden breakdown of the pre-existing structures and a heavier tail in the dissipation probability density function distribution. Different mechanisms are at play during this process: the high shear present at the interface between coherent low- and high-momentum regions is exploited to break down the larger structures and drive energy to small scales. This energy cascade is fed by an enhanced lift-up effect that produces intense streaks near the wall. It is found that the optimal perturbation grows exponentially during the first phase of its evolution reflecting the existence of a secondary modal instability of the streaks. To corroborate the results, the conditional spatiotemporal proper orthogonal decomposition (POD) analysis of Hack & Schimdt (J. Fluid Mech., vol. 907, 2021, A9) is performed both in the perturbed and in the unperturbed flow, showing a clear agreement between the two cases and with the reference study. Thus, the optimal perturbation at initial time can be considered as a precursor of extreme events.

The equation for the fluid velocity gradient along a Lagrangian trajectory immediately follows from the Navier–Stokes equation. However, such an equation involves two terms that cannot be determined from the velocity gradient along the chosen Lagrangian path: the pressure Hessian and the viscous Laplacian. A recent model handles these unclosed terms using a multi-level version of the recent deformation of Gaussian fields (RDGF) closure (Johnson & Meneveau, Phys. Rev. Fluids, vol. 2 (7), 2017, 072601). This model is in remarkable agreement with direct numerical simulations (DNS) data and works for arbitrary Taylor Reynolds numbers $\textit {Re}_\lambda$. Inspired by this, we develop a Lagrangian model for passive scalar gradients in isotropic turbulence. The equation for passive scalar gradients also involves an unclosed term in the Lagrangian frame, namely the scalar gradient diffusion term, which we model using the RDGF approach. However, comparisons of the statistics obtained from this model with DNS data reveal substantial errors due to erroneously large fluctuations generated by the model. We address this defect by incorporating into the closure approximation information regarding the scalar gradient production along the local trajectory history of the particle. This modified model makes predictions for the scalar gradients, their production rates, and alignments with the strain-rate eigenvectors that are in very good agreement with DNS data. However, while the model yields valid predictions up to $\textit {Re}_\lambda \approx 500$, beyond this, the model breaks down.

Turbulence in the restricted nonlinear (RNL) dynamics is analysed and compared with direct numerical simulations (DNS) of Poiseuille turbulence at Reynolds number $R=1650$. The structures are obtained by proper orthogonal decomposition (POD) analysis of the two components of the flow partition used in RNL dynamics: the streamwise mean flow and fluctuations. POD analysis of the streamwise mean flow indicates that the dominant POD modes, in both DNS and RNL dynamics, are roll-streaks harmonic in the spanwise direction. However, we conclude that these POD modes do not occur in isolation but rather are Fourier components of a coherent roll-streak structure. POD analysis of the fluctuations in DNS and RNL dynamics reveals similar complex structures consisting in part of oblique waves collocated with the streak. The origin of these structures is identified by their correspondence to POD modes predicted using a stochastic turbulence model (STM). These predicted POD modes are dominated by the optimally growing structures on the streak, which the STM predicts correctly to be of sinuous oblique wave structure. This close correspondence between the roll-streak structure and the associated fluctuations in DNS, RNL dynamics and the STM implies that the self-sustaining mechanism operating in DNS is essentially the same as that in RNL dynamics, which has been associated previously with optimal perturbation growth on the streak.

Despite elegant theoretical descriptions and numerous potential applications in nature, the various features of wave turbulence and its associated transfers across scales are still debated. One reason is the difficulty, both experimentally and numerically, of studying a chaotic nonlinear wavy system with sufficient precision, over a sufficient duration, and with sufficient separation between the injection scale, the system scale and the dissipation scale, to dwell in detail on the description of a statistically converged state. The numerical study of Thomas & Ding (J. Fluid Mech., 2023, this issue) sheds new light on one key aspect of this global problem: the upscale transfer of waves in rotating shallow water equations, which is of relevance to atmosphere and ocean dynamics. Their results first confirm the theoretical predictions, with a robust inverse wave cascade, a predominant role of quartic resonances and a slope value of the wave spectrum in agreement with the expected range. But their deeper analysis of the results also highlights some weakness in the theoretical foundations, exhibiting strong intermittency and departure from ideal Gaussian statistics. This work thus calls for improved modelling, acknowledging that important large-scale climatic features could actually result from the piling up of small-scale waves, unresolved by typical Earth system models. Beyond rotating shallow water systems, this study more generically resonates with open questions in the field of wave turbulence and its geophysical applications, including recent research in deep rotating and stably stratified systems.

Exact budget equations are derived for the coherent and stochastic contributions to the second-order structure function tensor. They extend the anisotropic generalised Kolmogorov equations (AGKE) by considering the coherent and stochastic parts of the Reynolds stress tensor, and are useful for the statistical description of turbulent flows with periodic or quasi-periodic features, like, for example, the alternate shedding after a bluff body. While the original AGKE describe production, transport, inter-component redistribution and dissipation of the Reynolds stresses in the combined space of scales and positions, the new equations, called $\varphi$AGKE, contain the phase $\varphi$ as an additional independent variable, and describe the interplay among the mean, coherent and stochastic fields at the various phases. The newly derived $\varphi$AGKE are then applied to a case where an exactly periodic external forcing drives the flow: a turbulent plane channel flow modified by harmonic spanwise oscillations of the wall to reduce drag. The phase-by-phase action of the oscillating transversal Stokes layer generated by the forcing on the near-wall turbulent structures is observed, and a detailed description of the scale-space interaction among mean, coherent and stochastic fields is provided thanks to the $\varphi$AGKE.

We introduce maps of Cauchy–Green strain tensor eigenvalues to barycentric coordinates to quantify and visualize the full geometry of three-dimensional deformation in stationary and non-stationary fluid flows. As a natural extension of Lagrangian coherent structure diagnostics, which provide separate scalar fields and a one-dimensional quantification of fluid deformation, our barycentric mapping visualizes the role of all three Cauchy–Green eigenvalues (or rates of stretching) in a single plot through a novel stretching coordinate system. The coordinate system is based on the distance from three distinct limiting states of deformation that correspond with the dimension of the underlying invariant manifolds. One-dimensional axisymmetric deformation (sphere to rod deformation) corresponds to one-dimensional unstable manifolds, two-dimensional axisymmetric deformation (sphere to disk deformation) corresponds to two-dimensional unstable manifolds and the rare three-dimensional isometric case (sphere to sphere translation and rotation) corresponds to shear-free elliptic Lagrangian coherent structures (LCSs). We provide methods to visualize the degree to which fluid deformation approximates these limiting states, and tools to quantify differences between flows based on the compositional geometry of invariant manifolds in the flow. We also develop a simple analogue for bilinearly representing and plotting both rates of stretching and rotation as a single vector. As with other LCS techniques, these diagnostics define frame-indifferent material features in the flow. We provide multiple computed examples of LCS and momentum transport barriers, and show advantages over other coherent structure diagnostics.

A new time-dependent analysis of the global and local fluctuating velocity signals in grid turbulence is conducted to assess the scaling laws for non-equilibrium turbulence. Experimental datasets of static- and active-grid turbulence with different Rossby numbers $R_o({=}U/\varOmega M$: $U$ is the mean velocity, $\varOmega$ is the mean rotation rate and $M$ is the grid mesh size) are considered. Although the global (long-time-averaged) non-dimensional dissipation rate $C_\varepsilon$ is independent of the Reynolds number $Re_\lambda$ based on the global Taylor microscale, the local (short-time-averaged) non-dimensional dissipation rate $\left \langle C_\varepsilon (t_i) \right \rangle$ ($t_i$ is the local time) both in the static- and active-grid turbulence clearly show the non-equilibrium scaling $\left \langle C_\varepsilon (t_i)\right \rangle / \sqrt {Re_0} \propto \left \langle Re_\lambda (t_i) \right \rangle ^{-1}$ ($\left \langle Re_\lambda (t_i) \right \rangle$ and $Re_0$ are the Reynolds numbers based on the local Taylor microscale $\lambda (t_i)$ and the global integral length scale, respectively), which has only been confirmed for global statistics in the near field of grid turbulence. The local value of $\left \langle L(t_i) / \lambda (t_i) \right \rangle$ ($L(t_i)$ is the local integral length scale) shifts from the equilibrium to non-equilibrium scaling as $\left \langle Re_\lambda (t_i) \right \rangle$ increases, further confirming that the non-equilibrium scalings are recovered for local statistics both in the static- and active-grid turbulence. The local values of $\left \langle C_\varepsilon (t_i) \right \rangle$ and $\left \langle L(t_i) / \lambda (t_i) \right \rangle$ follow the theoretical predictions for global statistics (Bos & Rubinstein, Phys. Rev. Fluids, vol. 2, 2017, 022601).

We take advantage of the extremely small kinematic viscosity of superfluid $^4$He to investigate the propagation of macroscopic vortex rings at Reynolds numbers between $2 \times 10^4$ and $4 \times 10^6$. These inhomogeneous flow structures are thermally generated by releasing short power pulses into a small volume of liquid, open to the surrounding bath through a vertical tube $2$ mm in diameter. We study specifically the ring behaviour between $1.30$ and $1.80$ K using the flow visualization and second sound attenuation techniques. From the obtained data sets, containing more than $2600$ realizations, we find that the rings remain well-defined in space and time for distances up to at least $40$ tube diameters, and that their circulation depends significantly on the travelled distance, in a way similar to that observed for turbulent vortex rings propagating in Newtonian fluids. Additionally, the ring velocity and circulation appear to be influenced solely by a single, experimentally accessible parameter, combining the liquid temperature with the magnitude and duration of the power pulse. Overall, our results support the view that macroscopic vortex rings moving in superfluid $^4$He closely resemble their Newtonian analogues, at least in the absence of significant thermal effects and at sufficiently large flow scales.

Despite decades of investigations, there is still no consensus on whether inertial particles augment or dampen turbulence. Here, we perform the first experimental study in which the particle concentration is varied systematically across a broad range of volume fractions $\varPhi _{v}$, from nominally one-way coupled to heavily two-way coupled regimes, keeping all other parameters constant. We utilize a zero-mean flow chamber where steady, homogeneous and approximately isotropic air turbulence is realized, with a Taylor-microscale Reynolds number $Re_{\lambda } = 150\unicode{x2013}300$. We consider spherical solid particles of two sizes, both much smaller than the Kolmogorov length, and yielding Stokes numbers $St_{\eta } = 0.3$ and 2.6 based on the Kolmogorov time scale. By adjusting the turbulent intensity, the settling velocity parameter is kept constant for both cases, $Sv_{\eta } = V_{t}/u_{\eta }\approx 3$ (where $V_{t}$ is the still-air terminal velocity, and $u_{\eta }$ is the Kolmogorov velocity scale). Unlike previous studies focused on massively inertial particles, we find that the turbulent kinetic energy increases with particle loading, being more than doubled at $\varPhi _{v} =5\times 10^{-5}$. This is attributed to the energy input associated with gravitational settling: the particles release their potential energy into the fluid and increase its dissipation rate, while the time scale associated with the inter-scale energy transfer is not strongly changed. Two-point statistics indicate that the energy-containing eddies become vertically elongated in the presence of falling particles, and that the latter redistribute the energy more homogeneously across the scales compared to unladen turbulence. This is rooted in an enhanced cascade, as shown by the nonlinear inter-scale energy transfer rate.