Geometric Solution of Turbulence as Diffusion in Loop Space

TL;DR

This work presents loop-space calculus as a unified framework that recasts nonlinear turbulence and YM gradient flow into a linear diffusion problem over loop space, revealing the Euler ensemble as a universal attractor dual to a solvable string theory. It derives two intertwined spectra of intermittency and decay exponents tied to the nontrivial zeros of the Riemann zeta function and predicts log-periodic oscillations in correlation functions, supported by high-precision DNS and experimental hints. The framework extends to MHD turbulence, predicts a Prandtl-number–driven phase transition, and yields an exact Hodge-dual matrix surface solution to YM loop equations, establishing a deep duality between turbulence and gauge theory confinement. This geometric, number-theoretic perspective connects disparate nonlinear systems and offers concrete experimental and numerical tests for validating a universal attractor in decaying turbulence and a confinement-like structure in YM theory.

Abstract

Strongly nonlinear dynamics, from fluid turbulence to quantum chromodynamics, have long constituted some of the most challenging problems in theoretical physics. This review describes a unified theoretical framework, the loop space calculus, which offers an analytical approach to these problems. The central idea is a shift in perspective from pointwise fields to integrated loop observables, a transformation that recasts the governing nonlinear equations into a universal linear diffusion equation in the space of loops. This framework, supported by recent mathematical analysis, is analytically solvable and yields an exact, parameter-free solution for decaying hydrodynamic turbulence--the Euler ensemble--which is shown to be dual to a solvable string theory. The theory's predictions include: (i) the unification of spatial and temporal scaling laws, which are shown to be governed by two related, infinite spectra of intermittency and decay exponents, respectively, both derived from the nontrivial zeros of the Riemann zeta function; (ii) a first-order phase transition in magnetohydrodynamic (MHD) turbulence; and (iii) the formation of quantized, concentric shells in passive scalar mixing. The theory also predicts log-periodic oscillations in correlation functions, effects not described by standard phenomenology, for which there is now emerging evidence from high-precision turbulence experiments. The appearance of identical mathematical structures as solutions to the turbulent regime of Yang-Mills gradient flow points towards the broad applicability of this approach. The framework also yields a new type of minimal surface that solves the QCD loop equations in the limit of large loops.

Figures (7)

• Figure 1: Examples of the regular star polygons, $\{q/p\}$, that form the discrete target space of the dual string theory. The turbulent state, or Euler ensemble, is a statistical average over random walks on all such polygons.
• Figure 2: Verification of the Euler-ensemble theory against new high-resolution 4K DNS data SreeniAkash2025. (Top Left): The predicted parabolic relation between time $t$ and the integral length scale $L_{M}(t)$. The simulation data is in excellent agreement with a parabolic fit. (Top Right): The parameter-free prediction for the effective index $f(x)=\langle r\partial_{r}\log(\Delta v^{2})\rangle$. The theory captures the universal curvature seen in the DNS data. (Bottom): Verification of the full energy decay law at $Re_{\lambda}=145$. The DNS data (blue line) deviates from the leading-order term of the theory ($E \propto L_M^{-5/2}$, gray dash-dotted line). However, it is in excellent agreement with the complete theoretical solution (orange line), which includes all sub-leading exponents from the Mellin-Barnes integral. Only the two unknown scales (energy and length) were fitted to match the DNS data with the theory; no dimensionless parameters were changed. Such a perfect fit with only two parameters for about a thousand data points provides a high level of confidence.
• Figure 3: Log--log plot of the universal function $\xi^2\Phi\left(\lfloor*\rfloor{\frac{1}{ \xi}}\right)$ where $\xi = \frac{2 \pi r}{\left(\sqrt{2 \tilde{\nu} (t+ t_0)} -\sqrt{2 \tilde{\nu} t_0}\right)}$. Here $\Phi(n) = \sum_{0<p < n} \varphi(p)$ is the Euler totient summatory function.
• Figure 4: Three phases of the wavevector scale $|f(\Pr)|$ in MHD decaying turbulence. The red-dashed line represents a metastable phase.
• Figure 5: The 'ramp-cliff' structure in the time trace of temperature fluctuations ($\Delta\theta$) in a heated turbulent jet. These asymmetric patterns, characterized by a gradual rise (the 'ramp') followed by a sharp drop (the 'cliff'), are a key signature of large-scale coherent structures imprinting on small-scale scalar fields. This empirical observation is analogous to the sawtooth profile predicted for the quantized scalar shells in the Euler ensemble (Figure \ref{['fig::PhiXiPlot']}) if you reflect this picture across the right-hand side.( Figure adapted from sreenivasan2018turbulent

Theorems & Definitions (3)

• Theorem 1: Operator-Holonomy Identity
• proof
• proof