Powered by OpenAIRE graph
Found an issue? Give us feedback
image/svg+xml art designer at PLoS, modified by Wikipedia users Nina, Beao, JakobVoss, and AnonMoos Open Access logo, converted into svg, designed by PLoS. This version with transparent background. http://commons.wikimedia.org/wiki/File:Open_Access_logo_PLoS_white.svg art designer at PLoS, modified by Wikipedia users Nina, Beao, JakobVoss, and AnonMoos http://www.plos.org/ ZENODOarrow_drop_down
image/svg+xml art designer at PLoS, modified by Wikipedia users Nina, Beao, JakobVoss, and AnonMoos Open Access logo, converted into svg, designed by PLoS. This version with transparent background. http://commons.wikimedia.org/wiki/File:Open_Access_logo_PLoS_white.svg art designer at PLoS, modified by Wikipedia users Nina, Beao, JakobVoss, and AnonMoos http://www.plos.org/
ZENODO
Preprint
Data sources: ZENODO
addClaim

Energy Conservation, Cascade Stabilisation, and the Regularity of the 3D Navier-Stokes Equations: Computational Evidence from Galerkin Truncations

Authors: Higgins, Rod;

Energy Conservation, Cascade Stabilisation, and the Regularity of the 3D Navier-Stokes Equations: Computational Evidence from Galerkin Truncations

Abstract

We present computational evidence for the global regularity of the three-dimensional incompressible Navier–Stokes equations on the periodic torus T³ = (ℝ/2πℤ)³. Through the development of a multi-perspective scaffold array methodology — which measures the same Galerkin system from multiple truncation-level perspectives simultaneously — we discovered that the 3D spectral solver used in our investigation (and potentially in other spectral NS implementations) failed to conserve energy due to a missing imaginary factor −i in the Fourier-space trilinear coupling. This energy conservation failure caused spurious energy injection of 1–15% per unit time (completely independent of the time step Δt), producing enstrophy growth that was indistinguishable from genuine cascade blow-up. We correct this error by implementing complex Fourier coefficients with the full −i factor, achieving exact energy conservation: Σk Re(ūk · NLk) = 0 to machine precision at every truncation level. Three independent implementations (C, Python/NumPy, and scipy RK45) validate this result: initial energies agree to all digits, evolved energies agree to 9 × 10⁻⁶ relative, and the Taylor–Green vortex analytical solution is reproduced to 10⁻⁷. With the corrected solver, we observe that: The forward energy cascade stabilises at a finite wavenumber (N ≤ 14) for all tested initial conditions, with total energy monotonically decreasing and enstrophy bounded. All scaffold array contraction ratios satisfy ρ < 1 at every amplitude through A = 0.35 — no divergence is observed in any perspective. The cascade transfer rate into shell k satisfies the simple bound |Tk| ≤ 0.123 · E(t) · k3/2 (and the tighter analytical form |Tk| ≤ 0.031 · E(t) · Ω(t)1/2 · kγ−1 with γ < 2), verified at every shell and time point, with the ratio decreasing over time. We establish the following regularity mechanism: viscous diffusion at rate ν|k|² absorbs the energy cascade at every wavenumber because |k|² grows quadratically while the measured cascade transfer rate grows at most as k3/2 (and is empirically observed to decrease with k). Since total energy is finite and decreasing (dE/dt = −2νΩ ≤ 0, a mathematical identity), the cascade runs on a shrinking budget against an ever-stronger drain. The solution remains smooth because no finite-time concentration of energy is possible. We theorise that this mechanism extends to all smooth initial data and all ν > 0, and we claim that the energy conservation identity — when correctly implemented — is the structural property that prevents blow-up. Previous computational studies that did not verify energy conservation at ν = 0 may have been observing solver artefacts rather than genuine Navier–Stokes dynamics.

Powered by OpenAIRE graph
Found an issue? Give us feedback