Discontinuous transition from direct to inverse energy cascade in helically reweighted three-dimensional turbulence: critical characterization, mechanism in rotating/MHD/active flows, analytical theory, and experimental fingerprints

In a one-parameter modification of the three-dimensional Navier-Stokes equations that reweights heterochiral helical triadic interactions by $\lambda\in[0,1]$ while preserving energy, helicity, and all symmetries, Sahoo, Alexakis, and Biferale (Phys. Rev. Lett. 118, 164501, 2017) discovered that the turbulent energy cascade reverses direction at a critical weight $\lambda_c\approx 0.3$ in a discontinuous manner with diverging fluctuations in the large-Reynolds-number limit — the first discontinuous transition of cascade direction ever reported. The source leaves open four linked issues: the quantitative theory of the transition (order, scaling laws, critical exponents, universality); an analytical (renormalization-group or closure) theory reproducing it; whether the underlying helical-quenching mechanism (opposite-helicity modes decorrelating faster, dynamically quenching heterochiral nonlinearities) operates in physically realized systems such as rotating turbulence, magnetohydrodynamics, and active fluids; and experimental detection of the predicted fingerprints (sign change of the third-order structure function $S_3(r)$ for $r>r_f$, inverse-range spectra, amplified forcing-scale fluctuations). Nine years of follow-up — the authors' direct continuation measuring one exponent $E(k_f)\propto(\lambda-\lambda_c)^{-2/3}$ and flux-loop states, in which large homochiral and heterochiral partial energy fluxes flow in opposite directions while their difference equals the injection rate, a Physics Reports review, an EDQNM closure at the homochiral endpoint, rotating-turbulence phase diagrams, and active-fluid and MHD helical studies — has produced only partial progress: no theory of the transition, no demonstration of the mechanism in a physical system, and no experimental detection. The problem remains open.

Background

The direction of the turbulent energy cascade is one of the foundational organizing facts of fluid turbulence. In three-dimensional (3D) flows, kinetic energy injected at an intermediate scale is transported by the nonlinear term of the Navier-Stokes equations (NSE) predominantly toward smaller scales — a direct or forward cascade — where viscosity dissipates it, with the Kolmogorov energy spectrum $E(k)\propto k^{-5/3}$ ($E(k)$ is the kinetic energy per unit wavenumber $k$). In two-dimensional (2D) turbulence, a second inviscid quadratic invariant (the enstrophy, the mean-square vorticity) blocks forward energy transfer, and energy cascades inversely to ever larger scales, where it can accumulate in coherent condensates. What fixes cascade direction — dimension, inviscid invariants, symmetries, or the structure of the nonlinear interactions — is therefore a central question of turbulence theory, with direct consequences for geophysical and astrophysical flows.

Beyond energy, ideal 3D flows also conserve the helicity $H=\langle\mathbf{u}\cdot\boldsymbol{\omega}\rangle$, the correlation of the velocity $\mathbf{u}$ with the vorticity $\boldsymbol{\omega}=\nabla\times\mathbf{u}$. Expanding the velocity in the eigenbasis of the curl operator (the helical or Craya-Herring decomposition) assigns every Fourier mode a definite helicity sign and splits every nonlinear triad (a resonant triple of wavevectors coupled by the nonlinear term) into homochiral triads, whose three modes share one helicity sign, and heterochiral triads, whose signs are mixed. Alexakis (J. Fluid Mech. 812, 752-770, 2017) introduced a one-parameter family of equations — the $\lambda$-model — in which every heterochiral triadic interaction is multiplied by a weight $\lambda\in[0,1]$: $\lambda=1$ recovers the unmodified NSE, while $\lambda=0$ keeps only homochiral interactions. The model preserves both ideal invariants (energy and helicity) and all symmetries (rotations, reflections, scale transformations) of the NSE, so any change of cascade direction in it can be traced to the relative weight of homo- versus heterochiral interactions alone. At the homochiral endpoint $\lambda=0$ the flow develops a split cascade: energy is transferred inversely to large scales while helicity cascades forward.

Sahoo, Alexakis, and Biferale (Phys. Rev. Lett. 118, 164501, 2017) simulated this model at increasing Reynolds number $Re$ (the ratio of inertial to viscous forces; $Re\to\infty$ is the fully turbulent limit) and box size. They found that for $\lambda>\lambda_c\approx 0.3$ — close to $1/3$, where the cumulative weight of homochiral triads equals that of heterochiral ones — the cascade is purely forward with $E(k)\propto k^{-5/3}$, while for $\lambda<\lambda_c$ energy cascades inversely (with forward helicity transfer and an inverse-range spectrum $E(k)\propto k^{-7/3}$). At any finite $Re$ a bidirectional window coexists near $\lambda_c$, but it shrinks with growing $Re$ and box size, and the authors concluded that the transition is discontinuous in the infinite-$Re$ limit, with fluctuations of the energy flux and of the energy at the forcing scale growing sharply near criticality — the first discontinuous transition of cascade direction ever reported, in contrast with quasi-2D systems (for example the 2D magnetohydrodynamic analogue of Seshasayanan and Alexakis, Phys. Rev. E 93, 013104, 2016) where the analogous transition is continuous with measurable critical exponents. Because $\lambda$ introduces no new length scale, both inertial ranges (above and below the forcing scale) remain scale-invariant; the authors argue that this scale similarity forces the two ranges to share one cascade direction, which is the mechanism behind the discontinuity. A practical diagnostic is the third-order longitudinal structure function $S_3(r)=\langle(\delta u_\parallel(r))^3\rangle$ (the mean cube of the velocity increment over separation $r$), which the von Kármán-Howarth relation ties to cascade direction (negative for forward, positive for inverse transfer); in the model its sign changes across $\lambda_c$ for separations $r>r_f=2\pi/k_f$ ($k_f$ the forcing wavenumber), and it is easily measured in experiments.

The source's own conclusion defines what remains open. First, whether the proposed physical mechanism — opposite-helicity modes decorrelating faster than same-helicity ones, dynamically quenching heterochiral nonlinearities — operates in real helical systems: rotating turbulence (whose linear eigenmodes, the inertial waves, are precisely the helical modes, opposite signs propagating in opposite directions), magnetohydrodynamic (MHD) turbulence, and active fluids. Second, the theoretical characterization of the transition (order, scaling laws, critical behavior), since the simulations demonstrate but do not explain it. Third, the use of the $\lambda$-model as a benchmark for analytical turbulence theory — renormalization-group (RG) approaches or closures (phenomenological approximations that close the hierarchy of spectral correlations, such as EDQNM, the eddy-damped quasi-normal Markovian approximation). Fourth, experimental detection of the inverse-transfer fingerprints. Later work sharpened the theoretical stakes: the review by Alexakis and Biferale (Phys. Rep. 767-769, 1-101, 2018) lists the absence of theories that "support the existence of criticality, explain quantitatively how these out-of-equilibrium (phase) transitions take place and predict their exponents and their universality class (if any)" among the field's leading open problems.

Problem Statement

The problem, as posed by the concluding discussion of Sahoo, Alexakis, and Biferale (Phys. Rev. Lett. 118, 164501, 2017) and sharpened by the follow-up literature, is to turn the empirically discovered discontinuous cascade-reversal transition of the $\lambda$-reweighted 3D Navier-Stokes equations into established physics. Its four requested results, raised as separate open directions by the source and its follow-up literature, share one objective: establishing the transition as a quantitatively characterized out-of-equilibrium phase transition whose mechanism and fingerprints reach beyond the decimated model. The four results are separately answerable sub-goals: each may be pursued and judged on its own, and settling any one of them is a valid partial resolution of the problem, while a complete solution to the problem as posed settles all four. In particular, for item 3, resolving the mechanism question for at least one of the three physical systems is itself a valid resolution of that item, while a complete account covers each system.

1. Critical characterization. Establish independently, by derivation, simulation, or both, the order of the transition in the combined infinite-Reynolds-number, infinite-box limit, and its quantitative critical behavior: the scaling laws of the energy at the forcing scale $E(k_f)$ and of the flux fluctuations as $\lambda\to\lambda_c$, finite-size and finite-Reynolds scaling functions, critical exponents, and whether a universality class exists (and, if so, which). The order and the exponents must be obtained in this work; restating the values already reported in the source and its follow-up literature does not discharge this item.

2. Analytical theory. Construct an analytical theory of the $\lambda$-Navier-Stokes equations — the source suggests renormalization-group approaches or closures as examples — that derives, rather than presupposes, the existence and location of $\lambda_c$, the discontinuous character of the transition, and the growth of fluctuations near criticality; accounting for the quantitative behavior requested in item 1 is a desirable target but is not required for this item, which may be judged on whether the theory predicts the transition and its character. The source poses this model explicitly as a benchmark against which analytical theories of turbulence should be validated.

3. Mechanism reach in physical systems. Determine whether the helical-quenching mechanism — opposite-helicity modes decorrelating faster than same-helicity modes, dynamically quenching heterochiral nonlinearities — actually operates in physically realized helical systems (rotating turbulence, where the linear eigenmodes are helical inertial waves; MHD turbulence; active fluids), or whether the inverse cascades observed in those systems arise from other mechanisms. A natural diagnostic is the mode-resolved comparison of same- versus opposite-helicity sector correlations and of homochiral versus heterochiral partial energy fluxes, which would help distinguish helical quenching from competing mechanisms (quasi-2D resonant transfer, magnetic-helicity dynamics, odd viscosity — a non-dissipative chiral transport coefficient of chiral active fluids — and confinement).

4. Experimental fingerprints. Determine whether the predicted inverse-transfer fingerprints — the sign change of $S_3(r)$ for separations $r>r_f$, the coexistence of inverse-range ($k^{-7/3}$-type) and forward-range spectral slopes, and amplified fluctuations at the forcing scale — can be detected in laboratory or geophysical turbulent flows. Demonstrating the detectability of at least one of the three fingerprints in at least one laboratory or geophysical flow is itself a valid resolution of this item, while a complete account assesses each fingerprint.

References

1. G. Sahoo, A. Alexakis, and L. Biferale, "Discontinuous transition from direct to inverse cascade in three-dimensional turbulence", Physical Review Letters 118, 164501 (2017). DOI: 10.1103/PhysRevLett.118.164501; arXiv:1701.07699; https://doi.org/10.1103/PhysRevLett.118.164501
2. A. Alexakis, "Helically decomposed turbulence", Journal of Fluid Mechanics 812, 752-770 (2017). DOI: 10.1017/jfm.2016.831; https://doi.org/10.1017/jfm.2016.831
3. A. Alexakis and L. Biferale, "λ-Navier–Stokes turbulence", Philosophical Transactions of the Royal Society A 380, 20210243 (2022). DOI: 10.1098/rsta.2021.0243; arXiv:2109.01073; https://doi.org/10.1098/rsta.2021.0243
4. A. Alexakis and L. Biferale, "Cascades and transitions in turbulent flows", Physics Reports 767-769, 1-101 (2018). DOI: 10.1016/j.physrep.2018.08.001; arXiv:1808.06186; https://doi.org/10.1016/j.physrep.2018.08.001
5. S. J. Benavides and M. D. Bustamante, "Phase dynamics and their role determining energy flux in hydrodynamic shell models", Physical Review Fluids 11, 044601 (2026). DOI: 10.1103/y9cr-4698; arXiv:2507.03397; https://doi.org/10.1103/y9cr-4698
6. A. Briard, L. Biferale, and T. Gomez, "Closure theory for the split energy-helicity cascades in homogeneous isotropic homochiral turbulence", Physical Review Fluids 2, 102602 (2017). DOI: 10.1103/PhysRevFluids.2.102602; arXiv:1707.06191; https://doi.org/10.1103/PhysRevFluids.2.102602
7. A. van Kan and A. Alexakis, "Critical transition in fast-rotating turbulence within highly elongated domains", Journal of Fluid Mechanics 899, A33 (2020). DOI: 10.1017/jfm.2020.443; arXiv:1912.05394; https://doi.org/10.1017/jfm.2020.443
8. P. Clark Di Leoni, A. Alexakis, L. Biferale, and M. Buzzicotti, "Phase transitions and flux-loop metastable states in rotating turbulence", Physical Review Fluids 5, 104603 (2020). DOI: 10.1103/PhysRevFluids.5.104603; https://doi.org/10.1103/PhysRevFluids.5.104603
9. A. Alexakis, P. D. Mininni, and R. Marino, "Energy cascades in rotating and stratified turbulence in anisotropic domains", Journal of Fluid Mechanics 1031, A7 (2026). DOI: 10.1017/jfm.2026.11301; arXiv:2508.00340; https://doi.org/10.1017/jfm.2026.11301
10. J. Słomka and J. Dunkel, "Spontaneous mirror-symmetry breaking induces inverse energy cascade in 3D active fluids", Proceedings of the National Academy of Sciences USA 114, 2119-2124 (2017). DOI: 10.1073/pnas.1614721114; https://doi.org/10.1073/pnas.1614721114
11. J. Słomka, P. Suwara, and J. Dunkel, "The nature of triad interactions in active turbulence", Journal of Fluid Mechanics 841, 702-731 (2018). DOI: 10.1017/jfm.2018.108; arXiv:1708.01394; https://doi.org/10.1017/jfm.2018.108
12. X. M. de Wit, M. Fruchart, T. Khain, F. Toschi, and V. Vitelli, "Pattern formation by turbulent cascades", Nature 627, 515-521 (2024). DOI: 10.1038/s41586-024-07074-z; arXiv:2304.10444; https://doi.org/10.1038/s41586-024-07074-z
13. M. Linkmann, G. Sahoo, M. McKay, A. Berera, and L. Biferale, "Effects of magnetic and kinetic helicities on the growth of magnetic fields in laminar and turbulent flows by helical Fourier decomposition", The Astrophysical Journal 836, 26 (2017). DOI: 10.3847/1538-4357/836/1/26; https://doi.org/10.3847/1538-4357/836/1/26
14. F. Plunian, A. Teimurazov, R. Stepanov, and M. K. Verma, "Inverse cascade of energy in helical turbulence", Journal of Fluid Mechanics 895, A13 (2020). DOI: 10.1017/jfm.2020.307; arXiv:2005.04701; https://doi.org/10.1017/jfm.2020.307
15. P. Z. Stasiak, A. Baggaley, C. F. Barenghi, G. Krstulovic, and L. Galantucci, "Inverse energy transfer in three-dimensional quantum vortex flows" (2025). DOI: 10.48550/arXiv.2503.04450; arXiv:2503.04450; https://arxiv.org/abs/2503.04450
16. C. Zhou, H.-S. Dou, L. Niu, and W. Xu, "Inverse energy cascade in turbulent Taylor-Couette flows", Physics of Fluids 37, 014110 (2025). DOI: 10.1063/5.0250908; https://doi.org/10.1063/5.0250908
17. S. Angriman, P. J. Cobelli, M. Bourgoin, S. G. Huisman, R. Volk, and P. D. Mininni, "Broken mirror symmetry of tracer's trajectories in turbulence", Physical Review Letters 127, 254502 (2021). DOI: 10.1103/PhysRevLett.127.254502; https://doi.org/10.1103/PhysRevLett.127.254502
18. V. B. Jha, K. Seshasayanan, and V. Dallas, "Cascades transition in generalised two-dimensional turbulence", Journal of Fluid Mechanics 1008, A23 (2025). DOI: 10.1017/jfm.2025.115; https://doi.org/10.1017/jfm.2025.115
19. W. J. T. Bos, "Cascade of mesostrophy in turbulence with reduced vortex stretching", Physical Review Fluids 11, 074605 (2026). DOI: 10.1103/qhck-4dky; https://doi.org/10.1103/qhck-4dky
20. K. Seshasayanan and A. Alexakis, "Critical behavior in the inverse to forward energy transition in two-dimensional magnetohydrodynamic flow", Physical Review E 93, 013104 (2016). DOI: 10.1103/PhysRevE.93.013104; https://doi.org/10.1103/PhysRevE.93.013104