This is the second paper in a cycle investigating the exact solution of loop equations in decaying turbulence. We perform numerical simulations of the Euler ensemble, suggested in the previous work, as a solution to the loop equations. We designed novel algorithms for simulation, which take a small amount of computer RAM so that only the CPU time grows linearly with the size of the system and its statistics. This algorithm allows us to simulate systems with billions of discontinuity points on our fractal curve, dual to the decaying turbulence. For the vorticity correlation function, we obtain \textbf{quantum fractal laws} with regime changes between universal discrete values of indexes. The traditional description of turbulence with fractal or multifractal scaling laws does not apply here. The measured conditional probabilities of fluctuating variables are smooth functions of the logarithm of scale, with statistical errors being negligible at our large number of random samples $T = 167,772,160$ with $N = 200,000,000$ points on a fractal curve. The quantum jumps arise from the analytical distribution of scaling variable $X = \frac{\cot^2(\pi p/q)}{q^2}$ where $\frac{p}{q}$ is a random fraction. This distribution is related to the totient summatory function and is a discontinuous function of $X$. In particular, the energy spectrum is a devil's staircase with tilted uneven steps with slopes between $-1$ and $-2.03$. The logarithmic derivative of energy decay $n(t)$ as a function of time is jumping down the stairs of universal levels between $1.50$ and $1.53$. The quantitative verification of this quantization would require more precise experimental data.

This paper brings to fruition recently found \cite{migdal2023exact} exact reduction of decaying turbulence in the Navier-Stokes equation in $3 + 1$ dimensions to a Number Theory problem of finding the statistical limit of the Euler ensemble. We reformulate the Euler ensemble as a Markov chain and show the equivalence of this formulation to the \textbf{quantum} statistical theory of $N$ fermions on a ring, with an external field related to the random fractions of $\pi$. We find the solution of this system in the turbulent limit $N\to \infty, \nu \to 0$ in terms of a complex trajectory (instanton) providing a saddle point to the path integral over the density of these fermions. This results in an analytic formula for the observable correlation function of vorticity in wavevector space. This is a full solution of decaying turbulence from the first principle without approximations or fitted parameters. The energy spectrum decays as $k^{-\frac{7}{2}}$. We compute resulting integrals in \Mathematica{} and present effective index $n(t) = -t\partial_t \log E $ for the energy decay as a function of time Fig.\ref{fig::NPlot}. The asymptotic value of the effective index in energy decay $n(\infty) = \frac{5}{4}$, but the universal function $n(t)$ is neither constant nor linear.

We have found an infinite dimensional manifold of exact solutions of the Navier-Stokes loop equation for the Wilson loop in decaying Turbulence in arbitrary dimension $d >2$. This solution family is equivalent to a fractal curve in complex space $\mathbb C^d$ with random steps parametrized by $N$ Ising variables $\sigma_i=\pm 1$, in addition to a rational number $\frac{p}{q}$ and an integer winding number $r$, related by $\sum \sigma_i = q r$. This equivalence provides a \textbf{dual} theory describing a strong turbulent phase of the Navier-Stokes flow in $\mathbb R_d$ space as a random geometry in a different space, like ADS/CFT correspondence in gauge theory. From a mathematical point of view, this theory implements a stochastic solution of the unforced Navier-Stokes equations. For a theoretical physicist, this is a \textbf{quantum} statistical system with integer-valued parameters, satisfying some number theory constraints. Its long-range interaction leads to critical phenomena when its size $N \rightarrow \infty$ or its chemical potential $\mu \rightarrow 0$. The system with fixed $N$ has different asymptotics at odd and even $N\rightarrow \infty$, but the limit $\mu \rightarrow 0$ is well defined. The energy dissipation rate is analytically calculated as a function of $\mu$ using methods of number theory. It grows as $\nu/\mu^2$ in the continuum limit $\mu \rightarrow 0$, leading to anomalous dissipation at $\mu \propto \sqrt{\nu} \to 0$. The same method is used to compute all the local vorticity distribution, which has no continuum limit but is renormalizable in the sense that infinities can be absorbed into the redefinition of the parameters. The small perturbation of the fixed manifold satisfies the linear equation we solved in a general form. This perturbation decays as $t^{-\lambda}$, where the anomalous dimensions $\lambda$ satisfy the spectral equation \eqref{spectralEq}. The spectrum becomes a continuum in the statistical limit $N \to \infty$, presumably leading to multifractal phenomena.

We develop a quantitative microscopic theory of decaying Turbulence by studying the dimensional reduction of the Navier-Stokes loop equation for the velocity circulation. We have found an infinite dimensional manifold of solutions of the Navier-Stokes loop equation\cite{M93, M23PR} for the Wilson loop in decaying Turbulence in arbitrary dimension $d >2$. This family of solutions corresponds to a fractal curve in complex space $\mathbb C^d$, described by an algebraic equation between consecutive positions plus a nonlinear periodicity condition. We derive the constrained SDE for the evolution of the fractal curve at a fixed moment of physical time as a function of an auxiliary stochastic time. We expect this stochastic process to cover our fixed manifold of the solutions of the decaying Turbulence. The energy density of the fluid decays as $\mathcal E_0/t$, where $\mathcal E_0$ is an initial dissipation rate. Presumably, we have found a new phase of extreme Turbulence yet to be observed in real or numerical experiments.

We study the Kelvinons: monopole ring solutions to the Euler equations, regularized as the Burgers vortex in the viscous core. There is finite anomalous dissipation in the inviscid limit. However, in the anomalous Hamiltonian, some terms are growing as logarithms of Reynolds number; these terms come from the core of the Burgers vortex. In our theory, the turbulent multifractal phenomenon is similar to asymptotic freedom in QCD, with these logarithmic terms summed up by an RG equation. The small effective coupling does not imply small velocity; on the contrary, velocity is large compared to its fluctuations, which opens the way for a quantitative theory. In the leading order in the perturbation theory in this effective coupling constant, we compute running multifractal dimensions for high moments of velocity circulation, in good agreement with the data for quantum Turbulence and available data for classical Turbulence. The logarithmic dependence of fractal dimensions on the loop size comes from the running coupling in anomalous dimensions. This slow logarithmic drift of fractal dimensions would be barely observable at Reynolds numbers achievable at modern DNS.