A central problem in turbulence is understanding small-scale intermittency, which refers to the sporadic generation of intense fluctuations in velocity gradients and increments. These extreme events, strongly non-Gaussian in nature, govern dissipation, mixing and transport processes in virtually all turbulent flows. Yet, despite decades of study, a faithful and predictive characterisation of small scales remains elusive owing to the inherent mathematical intractability of the Navier–Stokes equations and the difficulty in resolving them in both simulations and experiments at high Reynolds numbers. Recent advances in high-resolution simulations and experiments have significantly reshaped this picture, particularly by providing precise data at high Reynolds numbers to probe the full tensorial structure and dynamics at small scales. In this article, we synthesise the current understanding of small-scale intermittency and universality, drawing on modern data from well-resolved simulations and experiments that resolve the full velocity-gradient tensor. The results show that, while prevailing intermittency theories capture several key trends, they fail to describe or account for observed asymmetries between longitudinal and transverse fluctuations or between strain and vorticity amplification. Evidence suggests that intermittency is closely tied to the dynamics and geometry of vorticity and strain fields, with non-locality playing an important role. We argue that a consistent picture has emerged, but a complete theory will require unifying the statistical scaling frameworks with the underlying dynamical mechanisms that govern gradient amplification. Additional implications of these findings are discussed, and several pressing open problems are identified for future work.

Turbulent pipe flow is of substantial importance in practical applications, and it remains challenging to depict the characteristic complex multiscale dynamics by a unified theoretical framework, hindered by its inherent intermittency. Inspired by a recent study of velocity circulation in turbulent channel flows from Duan, Chen & Sreenivasan (2025 J. Fluid Mech., vol. 1009, p. R4), in this study, we investigate the statistical characteristics of velocity circulation (or equally the area integral of wall-normal vorticity) over rectangular loops in concentric cylindrical shells, parallel to the pipe wall. The statistics are implemented using direct numerical simulation data at friction Reynolds numbers of $ \textit{Re}_\tau =1057$ and $2000$. Close to the pipe wall, the circulation in the inertial range resides on space-filling unifractal sets, with the Hölder exponent smaller than Kolmogorov’s $4/3$. Away from the pipe wall, the circulation displays bifractal characteristics and the Hölder exponents for high moment orders are very close to those reported in channel flows and homogeneous isotropic turbulence. The circulation statistics are only dependent on the area enclosed by the loops, and are invariant to the loop aspect ratio, once both edge lengths of the loops are in the inertial range.

Contrary to accepted turbulence folklore, which holds that no mathematical relation exists between the Navier–Stokes equations (NSEs) and the multifractal model (MFM) of Parisi and Frisch, we develop a theory that reconciles the MFM with Leray’s weak solutions of Navier–Stokes analysis. From a combination of Euler invariant scaling and the NSEs set in a three-dimensional box of side $L$, we also derive the Paladin–Vulpiani scale $\eta_{h,pav}$ which is related to the Reynolds number Re by $L\eta _{h,\textit{pa}v}^{-1} = \textit{Re}^{1/(1+h)}$, and which acts as a mediator between the two theories. This is achieved by considering $L^{2m}$-norms of the velocity gradient to find a correspondence between $m$ and the local scaling exponent $h$ in the multifractal model. The parameter $m$ acts as if it were the sliding focus control on a telescope which allows us to zoom in and out on different structures. The range $1 \leqslant m \leqslant \infty$ is equivalent to $-{{ {2}/{3}}} \leqslant h_{\textit{min}} \leqslant {{{1}/{3}}}$, which lies precisely in the region where Bandak et al. (Phys. Rev. E, 2022, vol. 105, p. 065113; Phys. Rev. Lett., 2024, vol. 132, p. 104002) have suggested that thermal noise makes the NSEs inadequate and generates spontaneous stochasticity. The implications of this are discussed.

The anomalous scaling of passive scalar fluctuations is experimentally investigated in turbulent pipe flow with a Taylor-microscale Péclet number of $\mathcal{O}(10^5)$, where the turbulence is known to deviate from the homogeneous isotropic turbulence. The scalar structure functions and intermittency in the mixing are examined. The experimental results consolidate that the scaling exponents of scalar structure functions saturate at high-order even moments, evidenced previously in homogeneous isotropic turbulence with a Taylor-microscale Péclet number of $\mathcal{O}(10^3)$. The saturation scaling exponent decreases to approach unity as the Taylor-microscale Péclet number increases. This saturation scaling exponent is further corroborated by the fractal codimension of sharp scalar fronts.

Experiments were conducted to investigate the characteristics of turbulent spots formed in transitional boundary layers developed over a flat plate and an axisymmetric cone placed in similar hypersonic free-stream environment of Mach number 5.85. The free-stream unit Reynolds number in the present work varied in the range of $(3.0{-}6.0)\times 10^6$ m−1. Heat transfer measurement along the surface of both the test models was used to ascertain the state of boundary layer and to calculate the intermittency associated with the transitional boundary layer. Turbulent spots generated in the transitional boundary layer were characterised in terms of their leading–trailing-edge velocities, their streamwise length scales and their generation rates on both the test models. The leading edge of the turbulent spots developed over both the test models were found to be convecting at a speed equivalent to 90 % of the boundary layer edge speed. The trailing edge of the spots developed on a planar boundary layer traversed at a lower speed than its axisymmetric counterpart. Streamwise length scales of a turbulent spot developed in a planar boundary layer grew at a higher rate when compared with the axisymmetric boundary layer. Turbulent spot generation rates for both planar and axisymmetric boundary layers was found to be in the range of $10\,00\,000{-}30\,00\,000$ spots m−1/s−1.

This paper presents a theoretical and computational investigation into how a propagating three-dimensional vortex modifies ambient turbulence. Using rapid distortion theory and numerical simulations, the study explores both local and non-local changes in the external vorticity field resulting from fluid displacement and stretching. Cases involving structured and unstructured turbulence reveal that the vortex introduces permanent distortions along its path, and alters the far field turbulence through reflux effects. The findings extend classical models by quantifying the impact of vortex-induced strain and displacement on turbulence, offering new insights into turbulent–turbulent interfaces and the role of coherent structures in modulating external turbulent fields.

Slender fibres, including textile-derived microplastics, are abundant in aquatic environments and often extend beyond the Kolmogorov length scale. While breakup at dissipative scales has been characterised by velocity-gradient statistics, no closure existed for inertial-range spans where eddy turnover sets the clock. Here we develop a turbulence-informed kinetic theory of fibre fragmentation bridging turbulence forcing and slender-beam mechanics. First, we derive a load-to-curvature mapping showing that spanwise forcing generates peak bending moments scaling as $\sim U_L L^2$, with $U_L$ the velocity increment across fibre length $L$. Second, we construct a breakup hazard $h(L)$ from curvature-threshold exceedances over eddy-time blocks, which identifies a turbulence-defined critical span $\ell _c$. For $L\gt \ell _c$, breakup is eddy-time-limited, $h(L)=O(\bar \varepsilon ^{1/3}L^{-2/3})$ with $\bar \varepsilon$ the mean turbulent energy dissipation rate, whereas for $L\lt \ell _c$, it is a rare-event process with $h(L)\propto L^{5/3+\alpha }$, $\alpha$ denoting the small correction from intermittency. Embedding this hazard in a self-similar binary kernel yields a closed population-balance equation for the fragment distribution $n(L,t)$ with sources and sinks. The framework produces explicit predictions: intermittency-corrected curvature scalings, critical spans set by material and flow parameters, start-up and halving times linked to surf-zone conditions and scaling profiles in the cascade. The steady-state bulk distribution on the subcritical branch, with vertical removal induced by horizontal convergence, follows $n(L)\propto L^{-8/3-\alpha }\simeq L^{-2.7}$, in striking agreement with the mean slope $\simeq -2.68$ observed for environmental microfibres in recent surveys. The reported variability of slopes is naturally explained in our framework by the coexistence of supercritical and subcritical branches together with $L$-dependent removal-driven sinks.

We present three-dimensional velocity gradient statistics from turbulent Rayleigh–Bénard convection experiments in a horizontally extended cell of aspect ratio 25, a paradigm for mesoscale convection with its organisation into large-scale patterns. The Rayleigh number ${\textit{Ra}}$ ranges from $3.7 \times 10^5$ to $4.8 \times 10^6$, the Prandtl number ${\textit{Pr}}$ from 5 to 7.1. Spatio-temporally resolved volumetric data are reconstructed from moderately dense Lagrangian particle tracking measurements. All nine components of the velocity gradient tensor from the experiments show good agreement with those from direct numerical simulations, both conducted at ${\textit{Ra}} = 1 \times 10^6$ and ${\textit{Pr}} = 6.6$. As expected, with increasing ${\textit{Ra}}$, the flow in the bulk approaches isotropic conditions in the horizontal plane. The focus of our analysis is on non-Gaussian velocity gradient statistics. We demonstrate that statistical convergence of derivative moments up to the sixth order is achieved. Specifically, we examine the probability density functions (PDFs) of components of the velocity gradient tensor, vorticity components, kinetic energy dissipation and local enstrophy at different heights in the bottom half of the cell. The probability of high-amplitude derivatives increases from the bulk to the bottom plate. A similar trend is observed with increasing ${\textit{Ra}}$ at fixed height. Both indicate enhanced small-scale intermittency of the velocity field. We also determine derivative skewness and flatness. The PDFs of the derivatives with respect to the horizontal coordinates are found to be more symmetric than the ones with respect to the vertical coordinate. The conditional statistical analysis of the velocity derivatives with respect to up-/down-welling regions and the rest did not display significant difference, most probably due to the moderate Rayleigh numbers. Furthermore, doubly logarithmic plots of the PDFs of normalised energy dissipation and local enstrophy at all heights show that the left tails follow slopes of 3 / 2 and 1 / 2, respectively, in agreement with numerical results. In general, the left tails of the dissipation and local enstrophy distributions show higher probability values with increasing proximity towards the plate, in comparison with those in the bulk.

Wall pressure fluctuations (WPFs) over aerodynamic surfaces contribute to the physical origin of noise generation and vibrational loading. Understanding the generation mechanism of WPFs, especially those exhibiting extremely high amplitudes, is important for advancing design and control in practical applications. In this work, we systematically investigate extreme events of WPFs in turbulent boundary layers and the compressibility effects thereon. The compressibility effects, encompassing extrinsic and intrinsic ones, ranging from weak to strong, are achieved by varying Mach numbers and wall temperatures. A series of datasets at moderate Reynolds numbers obtained from direct numerical simulation are analysed. It is found that the intermittency of WPFs depends weakly on extrinsic compressibility effects, whereas intrinsic compressibility effects significantly enhance intermittency at small scales. Coherent structures related to extreme events are identified using volumetric conditional average. Under extrinsic compressibility effects, extreme events are associated with the weak dilatation structures induced by interactions of high- and low-speed motions. When intrinsic compressibility effects dominate, these events are associated with the strong alternating positive and negative dilatation structures embedded in low-speed streaks. Furthermore, Poisson-equation-based pressure decomposition is performed to partition pressure fluctuations into components governed by distinct physical mechanisms. By analysing the proportion of each pressure component in extreme events, it is found that the contributions of the slow pressure and viscous pressure exhibit weak dependence on the compressibility effects, especially the extrinsic ones, and the varying trend of contributions of the rapid pressure with compressibility effects is opposite to that of the compressible pressure component.

Characteristics of the turbulent/non-turbulent interface (TNTI) and entrainment in separated and reattaching flows induced by an oscillating fence are investigated using time-resolved particle image velocimetry. Disturbed flows are classified into subcritical, transitional, critical and supercritical cases based on the ratio of the oscillation frequency to the natural vortex shedding frequency. In the recirculation zone, distinct vortices across different cases lead to significant variations in TNTI characteristics. In the subcritical case, the TNTI evolution resembles that in the stationary fence case but with intensified height fluctuations due to the undulation of separated shear layer. For other cases, the mean TNTI height increases with the oscillation frequency, while height fluctuation diminishes. The TNTI thickness varies with nearby vortices, scaling with the Taylor microscale. After the reattachment, TNTI height distributions converge into two groups: subcritical and transitional cases exhibit larger fluctuations and positively skewed probability density functions (PDFs), while critical and supercritical cases show smaller fluctuations and basically symmetric PDFs. The TNTI thickness becomes consistent across various cases, matching the adjacent small-scale vortex size. Besides, the nibbling mechanism of entrainment aligns well with the flow development. The minimum mean entrainment velocity coincides with the strongest prograde vortex while the maximum occurs at $x\approx 1.2x_{{r}}$ (where $x$ denotes the streamwise coordinate and $x_{{r}}$ is the mean reattachment position) in all cases. Engulfment is enhanced near the reattachment location by oscillations in the transitional and critical cases, but is suppressed in the supercritical cases due to the weakness of vortex structures at higher oscillation frequencies.

In the standard picture of fully developed turbulence, highly intermittent hydrodynamic fields are nonlinearly coupled across scales, where local energy cascades from large scales into dissipative vortices and large density gradients. Microscopically, however, constituent fluid molecules are in constant thermal (Brownian) motion, but the role of molecular fluctuations in large-scale turbulence is largely unknown, and with rare exceptions, it has historically been considered irrelevant at scales larger than the molecular mean free path. Recent theoretical and computational investigations have shown that molecular fluctuations can impact energy cascade at Kolmogorov length scales. Here, we show that molecular fluctuations not only modify energy spectrum at wavelengths larger than the Kolmogorov length in compressible turbulence, but also significantly inhibit spatio-temporal intermittency across the entire dissipation range. Using large-scale direct numerical simulations of computational fluctuating hydrodynamics, we demonstrate that the extreme intermittency characteristic of turbulence models is replaced by nearly Gaussian statistics in the dissipation range. These results demonstrate that the compressible Navier–Stokes equations should be augmented with molecular fluctuations to accurately predict turbulence statistics across the dissipation range. Our findings have significant consequences for turbulence modelling in applications such as astrophysics, reactive flows and hypersonic aerodynamics, where dissipation-range turbulence is approximated by closure models.

Surface quasi-geostrophic (SQG) theory describes the two-dimensional active transport of a scalar field, such as temperature, which – when properly rescaled – shares the same physical dimension of length/time as the advecting velocity field. This duality has motivated analogies with fully developed three-dimensional turbulence. In particular, the Kraichnan – Leith – Batchelor similarity theory predicts a Kolmogorov-type inertial range scaling for both scalar and velocity fields, and the presence of intermittency through multifractal scaling was pointed out by Sukhatme & Pierrehumbert (2002 Chaos 12, 439–450), in unforced settings. In this work, we refine the discussion of these statistical analogies, using numerical simulations with up to $16\,384^2$ collocation points in a steady-state regime dominated by the direct cascade of scalar variance. We show that mixed structure functions, coupling velocity increments with scalar differences, develop well-defined scaling ranges, highlighting the role of anomalous fluxes of all the scalar moments. However, the clean multiscaling properties of SQG transport are blurred when considering velocity and scalar fields separately. In particular, the usual (unmixed) structure functions do no follow any power-law scaling in any range of scales, neither for the velocity nor for the scalar increments. This specific form of the intermittency phenomenon reflects the specific kinematic properties of SQG turbulence, involving the interplay between long-range interactions, structures and geometry. Revealing the multiscaling in single-field statistics requires us to resort to generalised notions of scale invariance, such as extended self-similarity and a specific form of refined self-similarity. Our findings emphasise the fundamental entanglement of scalar and velocity fields in SQG turbulence: they evolve hand in hand and any attempt to isolate them destroys scaling in its usual sense. This perspective sheds new lights on the discrepancies in spectra and structure functions that have been repeatedly observed in SQG numerics for the past 20 years.

To study the physics of small-scale properties of homogeneous isotropic turbulence at increasingly high Reynolds numbers, direct numerical simulation results have been obtained for forced isotropic turbulence at Taylor-scale Reynolds number $R_\lambda =2500$ on a $32\,768^3$ three-dimensional periodic domain using a GPU pseudo-spectral code on a 1.1 exaflop GPU supercomputer (Frontier). These simulations employ the multi-resolution independent simulation (MRIS) technique (Yeung & Ravikumar 2020, Phys. Rev. Fluids, vol. 5, 110517) where ensemble averaging is performed over multiple short segments initiated from velocity fields at modest resolution, and subsequently taken to higher resolution in both space and time. Reynolds numbers are increased by reducing the viscosity with the large-scale forcing parameters unchanged. Although MRIS segments at the highest resolution for each Reynolds number last for only a few Kolmogorov time scales, small-scale physics in the dissipation range is well captured – for instance, in the probability density functions and higher moments of the dissipation rate and enstrophy density, which appear to show monotonic trends persisting well beyond the Reynolds number range in prior works in the literature. Attainment of range of length and time scales consistent with classical scaling also reinforces the potential utility of the present high-resolution data for studies of short-time-scale turbulence physics at high Reynolds numbers where full-length simulations spanning many large-eddy time scales are still not accessible. A single snapshot of the $32\,768^3$ data is publicly available for further analyses via the Johns Hopkins Turbulence Database.

This work experimentally explores the alignment of the vorticity vector and the strain-rate tensor eigenvectors at locations of extreme upscale and downscale energy transfer. We show that the turbulent von Kármán flow displays vorticity–strain alignment behaviour across a large range of Reynolds numbers, which is very similar to previous studies on homogeneous, isotropic turbulence. We observe that this behaviour is amplified for the largest downscale energy transfer events, which tend to be associated with sheet-like geometries. These events are also shown to have characteristics previously associated with high flow field nonlinearity and singularities. In contrast, the largest upscale energy transfer events display very different structures which showcase a strong preference for vortex compression. Notably, in both cases we find that these trends are strengthened as the probed scales approach the Kolmogorov scale. We then show further evidence for the argument that strain self-amplification is the most salient feature in characterising the cascade direction. Finally, we identify possible invariant behaviour for the largest energy transfer events, even at scales near the Kolmogorov scale.

Turbulent flows exhibit large intermittent fluctuations from inertial to dissipative scales, characterised by multifractal statistics and breaking the statistical self-similarity. It has recently been proposed that the Navier–Stokes turbulence restores a hidden form of scale invariance in the inertial interval when formulated for a dynamically (nonlinearly) rescaled quasi-Lagrangian velocity field. Here we show that such hidden self-similarity extends to the large-eddy-simulation (LES) approach in computational fluid dynamics (CFD). In particular, we show that classical subgrid-scale models, such as implicit or explicit Smagorinsky closures, respect the hidden scale invariance at all scales – both resolved and subgrid. In the inertial range, they reproduce the hidden scale invariance of Navier–Stokes statistics. These properties are verified very accurately by numerical simulations and, beyond CFD, turn LES into a valuable tool for fundamental turbulence research.

The wake merging of two side-by-side porous discs with varying disc spacing is investigated experimentally in a wind tunnel. Two disc designs used in the literature are employed: a non-uniform disc and a mesh disc. Hot-wire anemometry is utilised to acquire two spanwise profiles at 8 and 30 disc diameters downstream and along the centreline between the dual-disc configuration up to 40 diameters downstream. The spanwise Castaing parameter profiles confirm the appearance of rings of internal intermittency at the outermost parts of the wakes. These rings are the first feature to interact between the discs. After this point, the turbulence develops to a state whereby an inertial range is observable in the spectra. Farther downstream, the internal intermittency shows the classical features of homogeneous, isotropic turbulence. These events are repeatable and occur in the same order for both types of porous discs. This robustness allows us to develop a general map of the merging of the two wakes.

Wall turbulence consists of various sizes of vortical structures that induce flow circulation around a wide range of closed Eulerian loops. Here we investigate the multiscale properties of circulation around such loops in statistically homogeneous planes parallel to the wall. Using a high-resolution direct numerical simulation database of turbulent channels at Reynolds numbers of $Re_\tau =180$, 550, 1000 and 5200, circulation statistics are obtained in planes at different wall-normal heights. Intermittency of circulation in the planes of the outer flow ($y^+ \gtrsim 0.1Re_\tau$) takes the form of universal bifractality as in homogeneous and isotropic turbulence. The bifractal character simplifies to space-filling character close to the wall, with scaling exponents that are linear in the moment order, and lower than those given by the Kolmogorov paradigm. The probability density functions of circulation are long-tailed in the outer bifractal region, with evidence showing their invariance with respect to the loop aspect ratio, while those in the inner region are closely Gaussian. The unifractality near the wall implies that the circulation there is not intermittent in character.

Intermittency as it occurs in fast dynamos in the magnetohydrodynamics (MHD) framework is evaluated through the examination of relations between normalized moments at third order (skewness $S$) and fourth order (kurtosis $K$) for both the velocity and magnetic field, and for their local dissipations. As investigated by several authors in various physical contexts such as fusion plasmas (Krommes 2008 Phys. Plasmas 15, 030703), climate evolution (Sura & Sardeshmukh 2008 J. Phys. Oceano. 38, 639-647), fluid turbulence or rotating stratified flows (Pouquet et al. 2023 Atmosphere 14, 01375), approximate parabolic $K(S)\sim S^\alpha$ laws emerge whose origin may be related to the applicability of intermittency models to their dynamics. The results analyzed herein are obtained through direct numerical simulations of MHD flows for both Taylor–Green and Arnold–Beltrami–Childress forcing at moderate Reynolds numbers, and for up to $3.14 \times 10^5$ turn-over times. We observe for the dissipation $0.2 \lesssim \alpha \lesssim 3.0$, an evaluation that varies with the field, the forcing and when filtering for high-skewness intermittent structures. When using the She & Lévêque (1994) Phys. Rev. Lett. 72, 336-339 intermittency model, one can compute $\alpha$ analytically; we then find $\alpha \approx 2.5$, clearly differing from a (strict) parabolic scaling, a result consistent with the numerical data.

Dispersion in spatio-temporal random flows is dominated by the competition between spatial and temporal velocity resets along particle paths. This competition admits a range of normal and anomalous dispersion behaviours characterised by the Kubo number, which compares the relative strength of spatial and temporal velocity resets. To shed light on these behaviours, we develop a Lagrangian stochastic approach for particle motion in spatio-temporally fluctuating flow fields. For space–time separable flows, particle motion is mapped onto a continuous time random walk (CTRW) for steady flow in warped time, which enables the upscaling and prediction of the large-scale dispersion behaviour. For non-separable flows, we measure Lagrangian velocities in terms of a new sampling variable, the average number of velocity transitions (both temporal and spatial) along pathlines, which renders the velocity series Markovian. Based on this, we derive a Lagrangian stochastic model that represents particle motion as a coupled space–time random walk, that is, a CTRW for which the space and time increments are intrinsically coupled. This approach sheds light on the fundamental mechanisms of particle motion in space–time variable flows, and allows for its systematic quantification. Furthermore, these results indicate that alternative strategies for the analysis of Lagrangian velocity data using new sampling variables may facilitate the identification of (hidden) Markov models, and enable the development of reduced-order models for otherwise complex particle dynamics.