High-Reynolds-number turbulence generates extreme velocity gradients whose statistical asymmetry — vorticity (enstrophy) is markedly more intermittent than strain (dissipation) despite exactly equal means — remains unexplained. Direct numerical simulations show that strain and vorticity are coupled nonlocally through the Biot–Savart relation: intense vorticity is amplified predominantly by nonlocal strain, and the conditional response of strain to intense vorticity grows sub-proportionately, ⟨Σ∣Ω⟩∝Ωγ\langle\Sigma|\Omega\rangle\propto\Omega^\gamma with 0<γ<10<\gamma<1, an empirical input that underlies the most successful current predictions of velocity-gradient PDF tails. The open problem is the dynamical origin of this asymmetric nonlocal coupling: a quantitative theory (or a controlled refutation) of the exponent γ\gamma and its Reynolds-number dependence, an explanation of why vorticity is more intermittent than strain and of how nearby but non-colocated vortex tubes and strain sheets cooperate to produce extreme gradients, and validation through incorporation of the resulting conditional statistics into statistical models and establishment of the high-Reynolds-number asymptotics.

Background

Turbulent flows at high Reynolds number spontaneously generate intense, rare fluctuations of velocity gradients — small-scale intermittency — which break Kolmogorov's 1941 mean-field picture and underlie the non-Gaussian statistics of dissipation, mixing, and transport. The velocity gradient tensor Aij=∂ui/∂xjA_{ij}=\partial u_i/\partial x_j decomposes into a symmetric part, the strain-rate tensor SijS_{ij} (the local rate of deformation), and an antisymmetric part characterized by the vorticity ω=∇×u\boldsymbol{\omega}=\nabla\times\mathbf{u} (twice the local angular velocity of a fluid element). Two scalar intensities follow: the enstrophy density Ω=ωiωi\Omega=\omega_i\omega_i, measuring rotation, and the pseudo-dissipation Σ=2SijSij\Sigma=2S_{ij}S_{ij}, proportional to the viscous dissipation rate of kinetic energy (ε=νΣ\varepsilon=\nu\Sigma, with ν\nu the kinematic viscosity). In homogeneous isotropic turbulence (HIT) the two have exactly equal means, ⟨Ω⟩=⟨Σ⟩=1/τK2\langle\Omega\rangle=\langle\Sigma\rangle=1/\tau_K^2, where τK\tau_K is the Kolmogorov microscale time; yet their probability density functions (PDFs) differ drastically in the tails — extreme vorticity events are far more probable than extreme strain events. That vorticity is more intermittent than strain is a long-standing observation, quantified using massive direct numerical simulations (DNS — exact numerical solution of the Navier–Stokes equations on grids resolving all scales, here up to 12288312288^3 grid points) by Buaria, Pumir, Bodenschatz and Yeung (2019).

The gradients evolve through three nonlinear mechanisms. Vortex stretching amplifies vorticity (ωiSijωj>0\omega_iS_{ij}\omega_j>0); strain self-amplification amplifies strain through the cubic invariant SijSjkSki>0S_{ij}S_{jk}S_{ki}>0; and the pressure Hessian (the tensor of second pressure derivatives, which via the pressure Poisson equation is a nonlocal functional of the whole velocity field) redistributes gradient fluctuations toward the mean. Buaria, Pumir and Bodenschatz (2022), analyzing DNS at Taylor-scale Reynolds number RλR_\lambda (the standard turbulence-intensity parameter built on the Taylor microscale) from 140 to 1300, established the statistical budget of intense dissipation: intense strain is self-amplified, depleted by vortex stretching, and depleted by the pressure Hessian, with the net amplification in the strain eigenbasis produced by the most compressive eigenvalue. Intense strain organizes into sheet-like structures adjacent to, but never co-located with, tube-like regions of intense vorticity.

The two fields are moreover nonlocally coupled: since incompressibility determines u\mathbf{u} from ω\boldsymbol{\omega} through the Biot–Savart integral (the inverse-curl relation), the strain at any point is a functional of the entire vorticity field. Decomposing strain into a part induced by vorticity inside a sphere of radius RR (local) and the rest (nonlocal) — a device introduced by Hamlington, Schumacher and Dahm (2008) — Buaria and Pumir (2021) showed that intense vorticity is amplified predominantly by nonlocal strain, while the locally induced strain attenuates it (the inviscid self-attenuation mechanism of Buaria, Pumir and Bodenschatz (2020)); the characteristic scale of the amplifying strain shrinks as a power law of the vorticity magnitude. All of this is measured phenomenology; none of it is yet a dynamical theory.

The sharpest quantitative expression of the missing theory is an asymmetry of conditional statistics (the conditional expectation ⟨X∣Y⟩\langle X|Y\rangle is the mean of XX given a fixed value of YY). In DNS of HIT, given intense strain, the accompanying enstrophy is proportionately intense, ⟨Ω∣Σ⟩∝Σ\langle\Omega|\Sigma\rangle\propto\Sigma; but given intense vorticity, the accompanying strain grows only sub-proportionately,

⟨Σ∣Ω⟩∝Ωγ,0<γ<1,\langle\Sigma|\Omega\rangle\propto\Omega^\gamma,\qquad 0<\gamma<1,

with γ≈0.6\gamma\approx0.6 at Rλ=140R_\lambda=140 rising slowly to ≈0.75\approx0.75 at Rλ=1300R_\lambda=1300 (Buaria and Pumir 2022). Because the strain that dissipates an intense vortex grows more slowly than the vorticity itself, vorticity fluctuations reach heavier distributional tails than strain; this asymmetry, combined with a Burgers-vortex-type balance (a steady stretched-vortex solution balancing viscosity against the strain acting on the filament), is the empirical input of the most successful current prediction of velocity-gradient PDF tails and of the smallest turbulence scales. The 2022 source paper posed the corresponding tasks as open: a quantitative analysis of the nonlocal strain–vorticity coupling; an explanation of why vorticity is more intermittent than strain; incorporation of the measured conditional power laws and prefactors (its Table II: ⟨X∣Σ⟩≈cXΣp\langle X|\Sigma\rangle\approx c_X\Sigma^p for the amplification terms XX) into statistical models such as PDF closures (probability-density-function transport methods used in turbulence modeling); and establishment of the asymptotic high-RλR_\lambda behavior of the conditional statistics and alignment. A 2026 synthesis review by the same authors confirms that the chief surviving issue is precisely the physical origin of the exponent γ\gamma and of its Reynolds-number dependence — how the vorticity field, through the Biot–Savart coupling, generates a strain field with these specific scaling properties.

Problem Statement

The problem is to convert the empirically measured nonlocal strain–vorticity coupling of high-Reynolds-number turbulence into a quantitative dynamical theory — grounded in the incompressible Navier–Stokes equations and the Biot–Savart relation — with one overarching objective: to explain, rather than measure, how nonlocal coupling between strain and vorticity generates extreme velocity gradients. A complete answer comprises three linked components.

(a) Origin of the asymmetric conditional coupling. Explain quantitatively why conditioning on intense strain and on intense vorticity is asymmetric — ⟨Ω∣Σ⟩∝Σ\langle\Omega|\Sigma\rangle\propto\Sigma but ⟨Σ∣Ω⟩∝Ωγ\langle\Sigma|\Omega\rangle\propto\Omega^\gamma with 0<γ<10<\gamma<1 in the far tail (ΩτK2≫1\Omega\tau_K^2\gg1) — by deriving, bounding, or otherwise predicting from the dynamics the exponent γ\gamma and its dependence on the Taylor-scale Reynolds number (measured: γ≈0.6\gamma\approx0.6 at Rλ=140R_\lambda=140 rising to ≈0.75\approx0.75 at Rλ=1300R_\lambda=1300), or by a controlled refutation of this scaling. The scaling carries three independently refutable sub-claims, and a refutation must state which it overturns and which it leaves standing: (i) that the far-tail conditional response is sublinear and takes the measured exponents; (ii) that γ\gamma depends on RλR_\lambda as measured; (iii) that the reported velocity-gradient PDF-tail predictions rely on this conditional scaling. A refutation that removes the coupling's dynamical significance — showing the sublinear scaling to be an artifact of the fitted range carrying no dynamical constraint — resolves this component on its own and dispenses with components (b) and (c), whose premises it removes. A refutation confined to sub-claim (i), (ii), or (iii) resolves only the demand it targets, and components (b) and (c) must still be addressed for the scaling that survives.

(b) Mechanism of extreme-gradient generation and of the intermittency asymmetry. Identify which configurations of the vorticity field, through the Biot–Savart coupling, generate the observed conditional scaling of strain; explain how nearby but non-colocated tube-like vorticity structures and sheet-like strain structures cooperate to produce extreme gradients; and thereby resolve, as a consequence of the mechanism in (a), why vorticity (enstrophy) is more intermittent than strain (dissipation) despite their exactly equal means.

(c) Validation and completion. This component has two demands that must be satisfied jointly. First, show that the conditional power laws and prefactors implied by the theory, when incorporated into a statistical model of the dissipation and enstrophy fields (for example a PDF closure), reproduce the observed intermittency and non-Gaussian tail statistics of velocity gradients; this demand is discharged only by a stated, falsifiable comparison against the statistics being reproduced, identifying the quantity and range over which agreement is claimed — a qualitative similarity of shape is not sufficient. Second, establish the asymptotic (infinite-Reynolds-number) behavior of the conditional exponents, prefactors, and alignment statistics, including whether γ\gamma converges as RλR_\lambda increases beyond currently available values; this demand may be met under an explicit, declared closure or extrapolation assumption, provided that assumption is stated and its consequences for the exponents, prefactors, and alignment statistics are derived, and a purely fitted or undeclared extrapolation does not discharge it.

References

  1. Dhawal Buaria, Alain Pumir, and Eberhard Bodenschatz, 'Generation of intense dissipation in high Reynolds number turbulence', Philosophical Transactions of the Royal Society A 380, 20210088 (2022), DOI 10.1098/rsta.2021.0088, https://doi.org/10.1098/rsta.2021.0088 (arXiv:2107.01719)
  2. Dhawal Buaria, Alain Pumir, Eberhard Bodenschatz, and P. K. Yeung, 'Extreme velocity gradients in turbulent flows', New Journal of Physics 21, 043004 (2019), DOI 10.1088/1367-2630/ab0756, https://doi.org/10.1088/1367-2630/ab0756
  3. Peter E. Hamlington, Jörg Schumacher, and Werner J. A. Dahm, 'Direct assessment of vorticity alignment with local and nonlocal strain rates in turbulent flows', Physics of Fluids 20, 111703 (2008), DOI 10.1063/1.3021055, https://doi.org/10.1063/1.3021055
  4. Dhawal Buaria, Eberhard Bodenschatz, and Alain Pumir, 'Vortex stretching and enstrophy production in high Reynolds number turbulence', Physical Review Fluids 5, 104602 (2020), DOI 10.1103/PhysRevFluids.5.104602, https://doi.org/10.1103/PhysRevFluids.5.104602
  5. Dhawal Buaria, Alain Pumir, and Eberhard Bodenschatz, 'Self-attenuation of extreme events in Navier–Stokes turbulence', Nature Communications 11, 5852 (2020), DOI 10.1038/s41467-020-19530-1, https://doi.org/10.1038/s41467-020-19530-1
  6. Dhawal Buaria and Alain Pumir, 'Nonlocal amplification of intense vorticity in turbulent flows', Physical Review Research 3, L042020 (2021), DOI 10.1103/PhysRevResearch.3.L042020, https://doi.org/10.1103/PhysRevResearch.3.L042020
  7. Dhawal Buaria and Alain Pumir, 'Vorticity-Strain Rate Dynamics and the Smallest Scales of Turbulence', Physical Review Letters 128, 094501 (2022), DOI 10.1103/PhysRevLett.128.094501, https://doi.org/10.1103/PhysRevLett.128.094501
  8. Dhawal Buaria and Alain Pumir, 'Role of pressure in the dynamics of intense velocity gradients in turbulent flows', Journal of Fluid Mechanics 973, A23 (2023), DOI 10.1017/jfm.2023.786, https://doi.org/10.1017/jfm.2023.786
  9. P.-F. Yang, H. Xu, A. Pumir, and G. W. He, 'Structure and role of the pressure Hessian in regions of strong vorticity in turbulence', Journal of Fluid Mechanics 983, R2 (2024), DOI 10.1017/jfm.2024.143, https://doi.org/10.1017/jfm.2024.143
  10. Dhawal Buaria, John M. Lawson, and Michael Wilczek, 'Twisting vortex lines regularize Navier-Stokes turbulence', Science Advances 10, eado1969 (2024), DOI 10.1126/sciadv.ado1969, https://doi.org/10.1126/sciadv.ado1969
  11. Dhawal Buaria and Alain Pumir, 'Universality of extreme events in turbulent flows', Physical Review Fluids 10, L042601 (2025), DOI 10.1103/PhysRevFluids.10.L042601, https://doi.org/10.1103/PhysRevFluids.10.L042601
  12. Dhawal Buaria and Alain Pumir, 'Turbulence intermittency and velocity gradients', Journal of Fluid Mechanics 1034, P1 (2026), DOI 10.1017/jfm.2026.11494, https://doi.org/10.1017/jfm.2026.11494
  13. G. E. Elsinga, T. Ishihara, and J. C. R. Hunt, 'Intermittency across Reynolds numbers – the influence of large-scale shear layers on the scaling of the enstrophy and dissipation in homogenous isotropic turbulence', Journal of Fluid Mechanics 974, A17 (2023), DOI 10.1017/jfm.2023.799, https://doi.org/10.1017/jfm.2023.799 (published title reads 'homogenous', sic)
  14. Miguel P. Encinar, 'A length scale for non-local multi-scale gradient interactions in isotropic turbulence', Journal of Fluid Mechanics 971, A40 (2023), DOI 10.1017/jfm.2023.706, https://doi.org/10.1017/jfm.2023.706
  15. Benjamin Musci, Bérengère Dubrulle, Jean LeBris, Damien Geneste, Pierre Bragança, Jean-Marc Foucaut, Christophe Cuvier, and Adam Cheminet, 'Experimental measurement of the vorticity–strain alignment around extreme energy transfer events', Journal of Fluid Mechanics 1016, A53 (2025), DOI 10.1017/jfm.2025.10447, https://doi.org/10.1017/jfm.2025.10447
  16. Perry L. Johnson and Michael Wilczek, 'Multiscale Velocity Gradients in Turbulence', Annual Review of Fluid Mechanics 56, 463–490 (2024), DOI 10.1146/annurev-fluid-121021-031431, https://doi.org/10.1146/annurev-fluid-121021-031431
  17. Rishita Das and Sharath S. Girimaji, 'Data-driven model for Lagrangian evolution of velocity gradients in incompressible turbulent flows', Journal of Fluid Mechanics 984, A39 (2024), DOI 10.1017/jfm.2024.235, https://doi.org/10.1017/jfm.2024.235