We investigate flame–acoustic interactions in a turbulent combustor during the state of intermittency before the onset of thermoacoustic instability using complex networks. Experiments are performed in a turbulent bluff-body stabilised dump combustor where the inlet airflow rate is varied quasi-statically and continuously. We construct a natural visibility graph from the local heat release rate fluctuations ($\dot {q}'$) at each location. Comparing the average degree during epochs of high- and low-amplitude acoustic pressure oscillations ($p'$) during the state of intermittency, we detect frequency modulation in $\dot {q}'$. Through this approach, we discover unique spatial patterns of cross-variable coupling between the frequency of $\dot {q}'$ and the amplitude of $p'$. The frequency of $\dot {q}'$ increases in regions of flame anchoring owing to high-frequency excitation of the flow and flame during epochs of high-amplitude $p'$ dynamics. However, the frequency of $\dot {q}'$ decreases in regions associated with flame-front distortions by large coherent vortices. In experiments with continuously varying airflow rates, the spatial pattern of frequency modulation varies with an increase in the average amplitude of $p'$ owing to an increase in the epochs of periodic $p'$ dynamics and the size of vortices forming in the flow. Dynamic shifts in the location of flame anchoring induce low-frequency fluctuations in $\dot {q}'$ during very-high-amplitude intermittent $p'$ dynamics. Our approach using conditional natural visibility graphs thus reveals the spatial pattern of amplitude-frequency coupling between the co-evolving flame and the acoustic field dynamics in turbulent reacting flows.

Turbulent flows in three dimensions are characterized by the transport of energy from large to small scales through the energy cascade. Since the small scales are the result of the nonlinear dynamics across the scales, they are often thought of as universal and independent of the large scales. However, as famously remarked by Landau, sufficiently slow variations of the large scales should nonetheless be expected to impact small-scale statistics. Such variations, often termed large-scale intermittency, are pervasive in experiments and even in simulations, while differing from flow to flow. Here, we evaluate the impact of temporal large-scale fluctuations on velocity, vorticity and acceleration statistics by introducing controlled sinusoidal variations of the energy injection rate into direct numerical simulations of turbulence. We find that slow variations can have a strong impact on flow statistics, raising the flatness of the considered quantities. We discern three contributions to the increased flatness, which we model by superpositions of statistically stationary flows. Overall, our work demonstrates how large-scale intermittency needs to be taken into account in order to ensure comparability of statistical results in turbulence.

Predicting and perhaps mitigating against rare, extreme events in fluid flows is an important challenge. Due to the time-localised nature of these events, Fourier-based methods prove inefficient in capturing them. Instead, this paper uses wavelet-based methods to understand the underlying patterns in a forced flow over a 2-torus which has intermittent high-energy burst events interrupting an ambient low-energy ‘quiet’ flow. Two wavelet-based methods are examined to predict burst events: (i) a wavelet proper orthogonal decomposition (WPOD) based method which uncovers and utilises the key flow patterns seen in the quiet regions and the bursting episodes; and (ii) a wavelet resolvent analysis (WRA) based method that relies on the forcing structures which amplify the underlying flow patterns. These methods are compared with a straightforward energy tracking approach which acts as a benchmark. Both the wavelet-based approaches succeed in producing better predictions than a simple energy criterion, i.e. earlier prediction times and/or fewer false positives and the WRA-based technique always performs better than WPOD. However, the improvement of WRA over WPOD is not as substantial as anticipated. We conjecture that this is because the mechanism for the bursts in the flow studied is found to be largely modal, associated with the unstable eigenfunction of the Navier–Stokes operator linearised around the mean flow. The WRA approach should deliver much better improvement over the WPOD approach for generically non-modal bursting mechanisms where there is a lag between the imposed forcing and the final response pattern.

This work presents an experimental investigation of the effects of vortex shedding suppression on the properties and recovery of turbulent wakes. Four plates, properly modified so that they produce different vortex shedding strengths, are tested using high speed particle image velocimetry and hot-wire anemometry, and analysed using spectral proper orthogonal decomposition, mean-flow linear stability analysis and various turbulence statistics. When present, vortex shedding is found to exhibit a characteristic frequency that scales with the mean shear, providing a link between the mean flow and the main turbulent motion. To achieve full suppression of shedding, we combine the effects of porosity and fractal perimeter. The mean shear is then decreased to the point where the flow becomes convectively unstable and shedding vanishes. In that case, the onset of self-similarity is delayed, compared with the case with vortex shedding, and appears after another large-scale structure, the secondary vortex street, emerges. It is also found that both large- and small-scale intermittency are starkly reduced when shedding is absent. A simple theoretical representation of the wake dynamics explains the evolution of the wake properties and its connection to the coherent structures in the flow.

Heavy particles suspended in turbulent flow possess inertia and are ejected from violent vortical structures by centrifugal forces. Once piled up along particle paths, this small-scale mechanism leads to an effective large-scale drift. This phenomenon, known as ‘turbophoresis’, causes particles to leave highly turbulent regions and migrate towards calmer regions, explaining why particles transported by non-homogeneous flows tend to concentrate near the minima of turbulent kinetic energy. It is demonstrated here that turbophoretic effects are just as crucial in statistically homogeneous flows. Although the average turbulent activity is uniform, instantaneous spatial fluctuations are responsible for inertial-range inhomogeneities in the particle distribution. Direct numerical simulations are used to probe particle accelerations, specifically how they correlate to local turbulent activity, yielding an effective coarse-grained dynamics that accounts for particle detachment from the fluid and ejection from excited regions through a space- and time-dependent non-Fickian diffusion. This leads to cast fluctuations in particle distributions in terms of a scale-dependent Péclet number ${\textit {Pe}}_\ell$, which measures the importance of turbulent advection compared with inertial turbophoresis at a given scale $\ell$. Multifractal statistics of energy dissipation indicate that $ {\textit {Pe}}_\ell \sim \ell ^\delta /\tau _{p}$ with $\delta \approx 0.84$. Numerical simulations support this behaviour and emphasise the relevance of the turbophoretic Péclet number in characterising how particle distributions, including their radial distribution function, depends on $\ell$. This approach also explains the presence of voids with inertial-range sizes, and the fact that their volumes have a non-trivial distribution with a power-law tail $p(\mathcal {V}) \propto \mathcal {V}^{-\alpha }$, with an exponent $\alpha$ that tends to 2 as ${\textit {Pe}}_\ell \to 0$.

This study investigates the influence of surface wave characteristics, specifically wave steepness and directional spreading, on intermittency in deep-water gravity wave turbulence through long-term numerical simulations of three-dimensional potential fully nonlinear periodic gravity waves. We conducted this investigation by estimating the scaling exponent of the surface elevation under different sea state conditions. With our numerical methods, we were able to evaluate the scaling exponents of the structure-function up to 12th order. The observed increased intermittency in directionally narrower sea states and in higher steepness conditions aligns with known effects of quasi-resonant wave–wave interactions and wave breaking. Comparative analyses reveal that both the conventional She–Leveque model and the multifractal models, also used to represent intermittency in wave turbulence of a different nature, exhibit a strong correlation in this study. This observation underscores the universality of intermittency phenomena within wave turbulence.

Turbulent entrainment at the turbulent/non-turbulent interface (TNTI) plays an important role in understanding the turbulent diffusion. While entrainment in fully developed canonical turbulent flows has been extensively studied, the evolution of entrainment in spatially developing flows remains poorly understood. In this work, characteristics of entrainment and the effect of vortices on entrainment of the shear layer separated from a wall-mounted fence are studied by the experiment in a water channel. The shedding vortex experiences a series of stages, including generation, growth, deformation and breakdown into smaller vortices. With the development of the flow, entrainment varies correspondingly. The prograde vortex near the TNTI is found to suppress entrainment but have little effect on the detrainment process, while the retrograde vortex promotes entrainment and suppresses detrainment as well. Consequently, the local entrainment velocity is decreased by the prograde vortex and increased more significantly by the retrograde vortex. Along the streamwise direction, the time-mean entrainment velocity is smallest where the prograde vortex is strongest in the vortex deformation stage. However, the largest time-mean entrainment velocity is located where the enstrophy gradient near the TNTI is greatest after reattachment, rather than where the retrograde vortex is strongest shortly after the breakdown of the shedding vortex, because the scarcity of retrograde vortices in the vicinity of the TNTI makes their long-time cumulative contribution not as significant as their local enhancement. The present study reveals how entrainment evolves in the separated and reattaching flow, and improves our understanding of the effect of vortices on entrainment.

Velocity gradient tensor, $A_{ij}\equiv \partial u_i/\partial x_j$, in a turbulence flow field is modelled by separating the treatment of intermittent magnitude ($A = \sqrt {A_{ij}A_{ij}}$) from that of the more universal normalised velocity gradient tensor, $b_{ij} \equiv A_{ij}/A$. The boundedness and compactness of the $b_{ij}$-space along with its universal dynamics allow for the development of models that are reasonably insensitive to Reynolds number. The near-lognormality of the magnitude $A$ is then exploited to derive a model based on a modified Ornstein–Uhlenbeck process. These models are developed using data-driven strategies employing high-fidelity forced isotropic turbulence data sets. A posteriori model results agree well with direct numerical simulation data over a wide range of velocity-gradient features and Reynolds numbers.

A spatially developing flat-plate boundary layer free from and two-way coupled with inertial solid particles is simulated to investigate the interaction between particles and the turbulent/non-turbulent interface. Particle Stokes numbers based on the outer scale are $St=2$ (low), 11 (moderate) and 53 (high). The Eulerian–Lagrangian point-particle approach is deployed for the simulation of particle-laden flow. The outer edge of the turbulent/non-turbulent interface layer is detected as an iso-surface of vorticity magnitude. Results show that the particles tend to accumulate below the interface due to the centrifugal effect of large-scale vortices in the outer region of wall turbulence and the combined barrier effect of potential flow. Consequently, the conditionally averaged fluid velocity and vorticity vary more significantly across the interface through momentum exchange and the feedback of force in the enstrophy transport. The large-scale structures in the outer layer of turbulence become smoother and less inclined in particle-laden flow due to the modulation of turbulence by the inertial particles. As a result, the geometric features of the interface layer are changed, namely, the spatial undulation increases, the fractal dimension decreases and the thickness becomes thinner in particle-laden flow as compared with unladen case. These effects become more pronounced as particle inertia increases.

The progressive destabilisation of the incompressible flow in a cubical cavity driven by a constant shear stress is investigated numerically. To that end, one of the square faces of the cube is subjected to a constant shear stress parallel to two opposite edges of that face. The three-dimensional steady basic flow loses its mirror symmetry through a supercritical pitchfork bifurcation leading to a pair of steady stable non-symmetric flow states that are antisymmetric to each other. Upon increase of the strength of the driving, these non-symmetric equilibria become unstable via a Hopf bifurcation resulting in two limit cycles. The bifurcations are investigated using classical linear stability analyses as well as nonlinear simulations. Upon a further increase of the driving shear stress, the limit cycles destabilise through bursts triggering a complex interplay between the unstable equilibria. The transition to turbulence resembles the Pomeau–Manneville scenario.

Direct numerical simulations up to Reλ = 1445 show that the scaling exponents for the enstrophy and the dissipation rate extrema are different and depend on the Reynolds number. A similar Reynolds number dependence of the scaling exponents is observed for the moments of the dissipation rate, but not for the moments of the enstrophy. Significant changes in the exponents occur at approximately Reλ ≈ 250, where Reλ is the Taylor based Reynolds number, which coincides with structural changes in the flow, in particular the development of large-scale shear layers. A model for the probability density functions (PDFs) of the enstrophy and dissipate rate is presented, which is an extension of our existing model (Proc. R. Soc. A, vol. 476, 2020, p. 20200591) and is based on the mentioned development of large-scale layer regions within the flow. This model is able to capture the observed Reynolds number dependencies of the scaling exponents, in contrast to the existing theories which yield constant exponents. Moreover, the model reconciles the scaling at finite Reynolds number with the theoretical limit, where the enstrophy and dissipation rate scale identically at infinite Reynolds number. It suggests that the large-scale shear layers are vital for understanding the scaling of the extrema. Furthermore, to reach the theoretical limit, the scaling exponents must remain Reynolds number dependent beyond the present Reλ range.

We investigate the role of pressure, via its Hessian tensor ${\boldsymbol {H}}$, on amplification of vorticity and strain-rate and contrast it with other inviscid nonlinear mechanisms. Results are obtained from direct numerical simulations of isotropic turbulence with Taylor-scale Reynolds number in the range 140–1300. Decomposing ${\boldsymbol {H}}$ into local isotropic (${\boldsymbol {H}}^{I}$) and non-local deviatoric (${\boldsymbol {H}}^{D}$) components reveals that ${\boldsymbol {H}}^{I}$ depletes vortex stretching, whereas ${\boldsymbol {H}}^{D}$ enables it, with the former slightly stronger. The resulting inhibition is significantly weaker than the nonlinear mechanism which always enables vortex stretching. However, in regions of intense vorticity, identified using conditional statistics, contribution from ${\boldsymbol {H}}$ prevails over nonlinearity, leading to overall depletion of vortex stretching. We also observe near-perfect alignment between vorticity and the eigenvector of ${\boldsymbol {H}}$ corresponding to the smallest eigenvalue, which conforms with well-known vortex-tubes. We discuss the connection between this depletion, essentially due to (local) ${\boldsymbol {H}}^{I}$, and recently identified self-attenuation mechanism (Buaria et al., Nat. Commun., vol. 11, 2020, p. 5852), whereby intense vorticity is locally attenuated through inviscid effects. In contrast, the influence of ${\boldsymbol {H}}$ on strain-amplification is weak. It opposes strain self-amplification, together with vortex stretching, but its effect is much weaker than vortex stretching. Correspondingly, the eigenvectors of strain and ${\boldsymbol {H}}$ do not exhibit any strong alignments. For all results, the dependence on Reynolds number is very weak. In addition to the fundamental insights, our work provides useful data and validation benchmarks for future modelling endeavours, for instance in Lagrangian modelling of velocity gradient dynamics, where conditional ${\boldsymbol {H}}$ is explicitly modelled.

We investigate the statistics of turbulence in emulsions of two immiscible fluids of the same density. We compute velocity increments between points conditioned to be located in the same phase or in different phases, and examine their probability density functions (PDFs) and the associated structure functions (SFs). This enables us to demonstrate that the presence of the interface reduces the skewness of the PDF at small scales and therefore the magnitude of the energy flux towards the dissipative scales, which is quantified by the third-order SF. The analysis of the higher-order SFs shows that multiphase turbulence is more intermittent than single-phase turbulence. In particular, the local scaling exponents of the SFs display a saturation below the Kolmogorov–Hinze scale, which indicates the presence of large velocity gradients across the interface. Interestingly, the statistics of the velocity differences in the carrier phase recovers that of single-phase turbulence when the viscosity of the dispersed phase is high.

The whistling induced by a low-Mach turbulent flow through a deep axisymmetric cavity in a duct is investigated theoretically and experimentally. The experiments include acoustic measurements and stereoscopic particle image velocimetry (PIV). The paper focuses on the effect of a mean swirl on the dynamics of the azimuthal aeroacoustic modes. The mean swirl in the cavity has two origins: one component is imposed by a controlled tangential air injection upstream of the cavity, and the other component spontaneously arises under the action of the self-sustained azimuthal aeroacoustic mode, as explained in the companion paper, Part 1 (Faure-Beaulieu, Xiong, Pedergnana & Noiray, J. Fluid Mech., vol. 971, 2023, A21). Experiments show that the dynamics of the aeroacoustic wave is influenced by the imposed swirl. In particular, the spinning wave propagating against the swirl is promoted. To explain this, a linear perturbation analysis is performed around an incompressible mean swirling flow obtained from RANS simulations. It reveals that the dominant shear layer modes of azimuthal order 1 and −1 involved in the whistling phenomenon are helical modes winding respectively with and against the swirl, and spinning respectively in counterswirl and co-swirl directions. The counterswirl hydrodynamic mode is the least damped of the two, which is in agreement with the experimental observations. Finally, a low-order model based on the wave equation is derived. With only a few parameters, it fully reproduces the experimental observations for a wide range of imposed swirl intensity in the duct flow, and it allows us to disentangle the mechanisms responsible for this complex aeroacoustic instability.

We study fluctuations of all co-existing energy exchange/transfer/transport processes in stationary periodic turbulence including those that average to zero and are not present in average cascade theories. We use a Helmholtz decomposition of accelerations that leads to a decomposition of all terms in the Kármán–Howarth–Monin–Hill (KHMH) equation (scale-by-scale two-point energy balance) causing it to break into two energy balances, one resulting from the integrated two-point vorticity equation and the other from the integrated two-point pressure equation. The various two-point acceleration terms in the Navier–Stokes difference (NSD) equation for the dynamics of two-point velocity differences have similar alignment tendencies with the two-point velocity difference, implying similar characteristics for the NSD and KHMH equations. We introduce the two-point sweeping concept and show how it articulates with the fluctuating interscale energy transfer as the solenoidal part of the interscale transfer rate does not fluctuate with turbulence dissipation at any scale above the Taylor length but with the sum of the time derivative and the solenoidal interspace transport rate terms. The pressure fluctuations play an important role in the interscale and interspace turbulence transfer/transport dynamics as the irrotational part of the interscale transfer rate is equal to the irrotational part of the interspace transfer rate and is balanced by two-point fluctuating pressure work. We also study the homogeneous/inhomogeneous decomposition of interscale transfer. The statistics of the latter are skewed towards forward cascade events whereas the statistics of the former are not. We also report statistics conditioned on intense forward/backward interscale transfer events.