Global Regularity of 3D Navier–Stokes: An Energy Argument

Rod Higgins · Senuamedia · v5, September 2026 · DOI: 10.5281/zenodo.22710241

Abstract

Current version: v5 · DOI 10.5281/zenodo.22710241

Contents

1. The Equation
2. Hypotheses and the Active-Shell Bound
3. Scope of the Hypotheses
4. Relation to Prior Work
5. Appendix: Decomposition of the Stretching
6. References

1. The Equation

The 3D incompressible Navier–Stokes equations on \(\mathbb{T}^3\):

with constant viscosity \(\nu > 0\). Kinetic energy and enstrophy:

Taking the \(L^2\) inner product of the equation with \(u\):

The nonlinear and pressure terms vanish by incompressibility. Integrating:

The enstrophy is integrable in time, \(\Omega \in L^1(0,\infty)\). This is compatible with \(\Omega(t)\) becoming unbounded on a short interval, which is exactly what a singularity would require.

The conditional statement. If \(\sup_{t \geq 0}\Omega(t) < \infty\), then \(u \in L^\infty(H^1) \hookrightarrow L^\infty(L^6)\) by Sobolev embedding in \(d = 3\), placing the solution in the Prodi–Serrin class \((p = \infty, q = 6)\), and the solution is smooth. Everything below concerns what would be needed to bound \(\Omega(t)\); no such bound is established for the total enstrophy.

2. Hypotheses and the Active-Shell Bound

In Fourier space, \(u = \sum_{\mathbf{k} \in \mathbb{Z}^3} \hat{u}_{\mathbf{k}}\,e^{i\mathbf{k}\cdot x}\). The shell energy at wavenumber \(K\) is

where \(B_K = \{K - \tfrac12 \le |\mathbf{k}| < K + \tfrac12\}\). The nonlinearity transfers energy between shells:

where \(\Pi_K\) is the energy flux from shells \(\leq K\) to shells \(> K\). The flux conserves total energy, \(\sum_K (\Pi_{K-1} - \Pi_K) = 0\); the viscous term removes \(2\nu K^2 E_K\) from each shell.

2.1 The three hypotheses

All three are quantified over the solution for all \(t \ge 0\). They are hypotheses on the cascade, not consequences of the equations; Section 3 gives explicit smooth flows on which (H1) and (H2) fail.

2.2 The active range under (H2)

Under (H2), a shell can sustain net forward transfer only if the flux can exceed the viscous term:

Call a shell satisfying (5) active, and write \(\mathcal A(t)\) for the active set.

2.3 The conditional active-shell bound

This is a statement about \(\Omega_{\mathcal A}\) only. Inactive shells (\(K > K_d\), or shells below the threshold (5)) can carry energy and enstrophy, and the argument gives no control over them; in particular it gives no bound on the total enstrophy \(\Omega(t)\). The Prodi–Serrin remark of Section 1 would apply to any uniform-in-time bound on \(\Omega\), not to (7).

3. Scope of the Hypotheses

Each of the following is an explicit computation on a smooth divergence-free field (64³ and 192³ pseudo-spectral evaluation, 2/3-dealiased; numerical/verify_review_claims.py, checks 1 and 2, and the \(N = 192\) script). They delimit where (H1) and (H2) can hold.

(H2) fails at \(t = 0\) for a datum with a spectral gap

For a smooth divergence-free field with one empty shell, \(E_{10} = 0\) while the flux through that shell is \(\Pi_{10} = -1.04\). The right-hand side of (4) vanishes and the left does not. A pointwise closure of the form (4) cannot hold for such data at any \(\alpha\).

(H2) fails for localised vortices at their own scale

For a vortex blob of size \(\ell\), the flux at wavenumber \(K \sim 1/\ell\) exceeds \(\alpha K E_K^{3/2}\) by a factor of order \(\ell^{-3}\): measured ratios are 13.6× to 55× for \(\ell = 0.25\) to \(0.15\). A closure that holds uniformly in \(K\) therefore excludes spatial concentration a priori, which is precisely the configuration a singularity would require.

(H2) fails in the dissipation range of generic data

In the exponential dissipation range of a generic smooth datum the ratio \(|\Pi_K|/(K E_K^{3/2})\) grows without bound as \(E_K\) decays: 7.5×, 30× and 98× at \(K = 40, 48, 52\).

(H1) fails for exact shear-flow solutions

For \(u = f(x_1)\,e_2\), an exact solution of the equations, every shell contains exactly two modes, so the dispersion constant is \(\rho_K = n_K/2\): 9, 31, 175, 625, 2517 at \(K = 1, 2, 5, 10, 20\). No uniform bound \(\rho_K \le \rho\) exists for this family.

The upshot: (H2) is a no-concentration closure, and (H1)–(H3) are quantified over a solution assumed smooth for all time. The bound (7) is a statement about flows whose cascade satisfies the hypotheses; it is not a route to global regularity.

4. Relation to Prior Work

The energy identity \(2\nu\int\Omega\,dt \leq E(0)\) has been known since Leray (1934). The enstrophy growth rate \(d\Omega/dt \leq C_\nu\Omega^3\) (cubic, from Gagliardo–Nirenberg) is sharp for arbitrary divergence-free vector fields (Lu and Doering 2008): no improvement is possible within the standard functional-analytic framework, and any improvement must come from structure specific to solutions.

Cheskidov and Shvydkoy (2014) formalised the Kolmogorov dissipation wavenumber as a regularity criterion: solutions are regular provided their dissipation wavenumber is \(L^{5/2}\)-integrable in time. That result is conditional on an a priori bound on the dissipation wavenumber; the active-range estimate (6) is a statement of the same type, conditional on (H1)–(H3).

Tao (2016) constructed an averaged Navier–Stokes equation that obeys the identical energy identity, has the same scaling and the same local existence theory, and blows up in finite time. Any argument for regularity must therefore use structure of the actual nonlinearity that the averaged equation does not share; hypotheses of the type (H1)–(H3) are one way of naming such structure, and Section 3 shows they are not available in general.

Constantin and Fefferman (1993) proved regularity under Lipschitz control of the vorticity direction \(\omega/|\omega|\); Miller (2020) used the tracelessness of the strain to improve quantitative enstrophy bounds. The companion notes on the vorticity direction (the direction identities and the coercivity inequality) continue that line.

Appendix: Decomposition of the Stretching

The vortex-stretching integral \(S = \int \omega\cdot\nabla u\cdot\omega\,dx\) can be split into local and non-local parts; the non-local part admits an estimate in which the tracelessness of the strain removes the isotropic component of \(\omega\otimes\omega\). This appendix records that decomposition and the estimate for the non-local term. It is a lemma about one term, not a bound on \(S\) or on \(d\Omega/dt\).

A.1 Decomposition

Decompose the velocity into low and high frequencies at a cutoff \(K_0\):

where \(u^<\) contains modes \(|\mathbf{k}| \leq K_0\) and \(u^>\) contains \(|\mathbf{k}| > K_0\). Then

\(S_{HH}\) is the high-frequency self-interaction; \(S_{HL}\) is low-frequency strain acting on high-frequency vorticity. The cross terms involve \(\omega^<\) and are bounded in terms of \(K_0^2 E(0)\).

A.2 The non-local term

Traceless strain. By incompressibility the strain \(S_{ij}^< = (\partial_i u_j^< + \partial_j u_i^<)/2\) is traceless, so

the isotropic part of \(\omega^>\otimes\omega^>\) contributes nothing to \(S_{HL}\).

Anisotropy. With \(a_{ij}^> = \omega_i^>\omega_j^>/|\omega^>|^2 - \delta_{ij}/3\) (traceless, \(\|a^>\|_{L^\infty} \leq 2/3\)),

Bounds. \(\|S^<\|_{L^2} \leq \|\nabla u^<\|_{L^2} \leq K_0\sqrt{2E(0)}\), and by Gagliardo–Nirenberg in 3D, \(\|\omega^>\|_{L^4}^2 \leq C\,(\Omega^>)^{1/4}(P^>)^{3/4}\). Combining gives (9). The exponent \(1/4\) on \(\Omega^>\) (rather than the \(3/4\) of the classical bound \(|S| \le c\,\Omega^{3/4}P^{3/4}\)) comes from replacing \(\|S\|_{L^2} \le \sqrt{2\Omega}\) by the low-frequency bound \(\|S^<\|_{L^2} \le K_0\sqrt{2E(0)}\). The local term \(S_{HH}\) is not improved by this argument, and the palinstrophy \(P^>\) is not controlled by the energy identity; the lemma therefore does not by itself yield a bound on \(d\Omega/dt\). ◼

References

1. J. Leray, Sur le mouvement d'un liquide visqueux emplissant l'espace, Acta Math. 63 (1934), 193–248.
2. E. Hopf, Über die Anfangswertaufgabe für die hydrodynamischen Grundgleichungen, Math. Nachr. 4 (1951), 213–231.
3. G. Prodi, Un teorema di unicità per le equazioni di Navier–Stokes, Ann. Mat. Pura Appl. 48 (1959), 173–182.
4. J. Serrin, On the interior regularity of weak solutions of the Navier–Stokes equations, Arch. Rational Mech. Anal. 9 (1962), 187–195.
5. O. A. Ladyzhenskaya, The Mathematical Theory of Viscous Incompressible Flow, 2nd English ed., Gordon and Breach, 1969.
6. A. N. Kolmogorov, The local structure of turbulence in incompressible viscous fluid for very large Reynolds numbers, Dokl. Akad. Nauk SSSR 30 (1941), 301–305.
7. R. H. Kraichnan, Inertial-range transfer in two- and three-dimensional turbulence, J. Fluid Mech. 47 (1971), 525–535.
8. J. T. Beale, T. Kato, and A. Majda, Remarks on the breakdown of smooth solutions for the 3-D Euler equations, Comm. Math. Phys. 94 (1984), 61–66.
9. F. Waleffe, The nature of triad interactions in homogeneous turbulence, Phys. Fluids A 4 (1992), 350–363.
10. P. Constantin and C. Fefferman, Direction of vorticity and the problem of global regularity for the Navier–Stokes equations, Indiana Univ. Math. J. 42 (1993), 775–789.
11. P. Constantin, C. Fefferman, and A. J. Majda, Geometric constraints on potentially singular solutions for the 3-D Euler equations, Comm. Partial Differential Equations 21 (1996), 559–571.
12. L. Escauriaza, G. Seregin, and V. Šverák, \(L^{3,\infty}\)-solutions of the Navier–Stokes equations and backward uniqueness, Russian Math. Surveys 58 (2003), 211–250.
13. G. L. Eyink, Locality of turbulent cascades, Physica D 207 (2005), 91–116.
14. L. Lu and C. R. Doering, Limits on enstrophy growth for solutions of the three-dimensional Navier–Stokes equations, Indiana Univ. Math. J. 57 (2008), 2693–2727.
15. L. Biferale, S. Musacchio, and F. Toschi, Inverse energy cascade in three-dimensional isotropic turbulence, Phys. Rev. Lett. 108 (2012), 164501.
16. L. Biferale and E. S. Titi, On the global regularity of a helical-decimated version of the 3D Navier–Stokes equations, J. Stat. Phys. 151 (2013), 1089–1098.
17. A. Cheskidov and R. Shvydkoy, A unified approach to regularity problems for the 3D Navier–Stokes and Euler equations: the use of Kolmogorov's dissipation range, J. Math. Fluid Mech. 16 (2014), 263–273.
18. T. Tao, Finite time blowup for an averaged three-dimensional Navier–Stokes equation, J. Amer. Math. Soc. 29 (2016), 601–674.
19. J. G. Ballouz and N. T. Ouellette, Tensor geometry in the turbulent cascade, J. Fluid Mech. 835 (2018), 1048–1064.
20. T. Buckmaster and V. Vicol, Nonuniqueness of weak solutions to the Navier–Stokes equation, Ann. of Math. 189 (2019), 101–144.
21. E. Miller, A regularity criterion for the Navier–Stokes equation involving only the middle eigenvalue of the strain tensor, Arch. Rational Mech. Anal. 235 (2020), 99–139.
22. J. G. Ballouz and N. T. Ouellette, Geometric constraints on energy transfer in the turbulent cascade, Phys. Rev. Fluids 5 (2020), 034603.
23. M. Carbone and A. D. Bragg, Is vortex stretching the main cause of the turbulent energy cascade?, J. Fluid Mech. 883 (2020), R2.
24. R. Higgins, Energy Conservation, Cascade Stabilisation, and the Global Regularity of the 3D Navier–Stokes Equations, 2026.
25. R. Higgins, Angular Relaxation on the Integer Lattice and Global Regularity of 3D Navier–Stokes, 2026.

Citation

@article{higgins2026energy,
  title={Global Regularity of 3D Navier--Stokes: An Energy Argument},
  author={Higgins, Rod},
  year={2026},
  note={v5},
  doi={10.5281/zenodo.22710241},
  url={https://lab.senuamedia.com/papers/energy-regularity.html}
}