The Research Position
The framework is the classical one: Type I limits are ancient Navier–Stokes solutions with bounded scale-invariant energy (Seregin; Koch–Nadirashvili–Seregin–Šverák), and gauged Type II limits are ancient Euler solutions. On the critical-concentration branch (vorticity ~ |x|−2 at the core), Grujić's evasion and transfer theorems (arXiv:2607.08866, 2609.05720) reduce regularity to a strain-aligned, inflow-controlled core with a logarithmic modulus on the vorticity direction; the alignment condition carries a smallness clause that is a condition on the core Reynolds number κ/ν.
Strategy classes closed by explicit fields
- Liouville theorems based on vorticity coherence. The Burgers vortex has exactly unidirectional vorticity and sustained stretching. Coherence does not stop a cascade.
- Necessary conditions based on a one-point direction statistic. A six-patch divergence-free field has an exactly isotropic direction second moment and enstrophy-weighted stretching 0.999.
- Regularity theorems hypothesising a pointwise-in-K spectral closure or uniform intra-shell dispersion. Both hypotheses fail for explicit smooth data: a spectral-gap datum, localised vortex blobs, the dissipation range, and the exact shear flow.
Where the difficulty sits
In the covariance between strain and vorticity direction at the flux scale, and in concentration. On the critical-concentration branch, in a core-Reynolds-number smallness for non-axisymmetric cores: the aligned homogeneous critical profiles are the swirl-free axisymmetric family, so the residual is a non-symmetric core whose direction is driven at no less than about 1% of κ/|x|2 per unit time.
The residual, stated exactly
- Liouville for ancient Type I Navier–Stokes solutions.
- Exclusion of a rotating |x|−2 core as an ancient Type I solution.
- Liouville for gauged Type II ancient Euler limits. Smooth compactly supported stationary Euler flows exist (Gavrilov 2019; Constantin–La–Vicol 2019), so any Liouville statement must use how the limit arose.
Computational Infrastructure and Measurements
The 3D Fourier–Galerkin triad kernel with complex coefficients and the −i factor conserves energy. The earlier real-coefficient implementation did not: it injected up to 15% at N = 7, and the defect was found and corrected on 2026-03-24. An FFT pseudo-spectral solver and the triad kernel agree to 2.8×10−17 at N = 4 over T = 0.01.
Instruments and what they measure
| Instrument | Measures | Validation |
|---|---|---|
| Triad kernel v3 (C) | 3D Galerkin truncation, exact Fourier trilinear coupling, Leray projection | Energy conserved; cross-checked against FFT solver to 2.8×10−17 |
| Spectral3D solver (Python) | 2/3-dealiased RK4 pseudo-spectral DNS | ν = 0 energy drift 10−7 |
A1 solver (shear_profile.py) | Angular shear of homogeneous critical profiles ω = |x|−2Ψ(x̂) | Validated to 10−14 |
E9 (pincer_probe.py) | Covariance term of the vortex stretching on DNS | Betchov's identity in explicit form |
verify_review_claims.py | Spectral-closure, intra-shell dispersion, amplitude scaling, and Burgers-vortex checks | Self-contained; NumPy only; seconds |
| Thermalisation check | Ω/E against the lattice mean 〈k2〉 in inviscid truncations | Within 1% at every N from 3 to 8; unchanged under time-step halving |
Research Programme: Seven Zenodo Records
Each record is linked at its current version. The one-line descriptions state what the record establishes.
Global Regularity of 3D Navier–Stokes: An Energy Argument
A conditional bound ΩA ≤ C E(0)5/3/ν4/3 on the active-shell enstrophy under three cascade hypotheses, together with the explicit smooth data for which those hypotheses fail.
Angular Relaxation on the Integer Lattice
An effective PDE for angular-resolved shell energy fitted to DNS in 1D, 2D and 3D, reproducing the 2D data within 2%; the 3D closure is a fitted description of the data, not a derived equation.
Energy Conservation and the 3D Galerkin Solver
The corrected energy-conserving 3D Galerkin solver, the history of its −i defect, and two lattice lemmas on triad geometry; the rigorous content is conditional and bounds nothing about a particular solution.
Unified Spectral Framework for Arithmetic Oscillations
Truncated Odlyzko–te Riele sums with computed ζ′(ρn) over 2,001,052 zeros: the Mertens sum crosses 1 at K = 58 and 2 near K = 1,040; zeta zeros recovered by Fourier analysis of M(x)/√x.
Scaffold Array Diagnostics: Euler, SQG & MHD
Galerkin truncations of Euler, MHD and SQG relax to the Lee–Kraichnan absolute equilibrium, which accounts for the observed enstrophy growth; a finite Galerkin system says nothing about blow-up of the three PDEs.
Shell Model, Diagnostic H, First 3D Kernel
A hand-constructed cubic shell model with a multi-level diagnostic H and the first 3D kernel whose mode-enumeration, triad-precomputation and Leray-projection structure became the basis of the corrected solver; no conclusion about 3D Navier–Stokes follows from it.
Unified Adaptation Theorem
Sufficient conditions for convergence of sequentially composed contractive updates on a shared parameter: block contraction with a non-negative coupling matrix plus LaSalle, with trajectory diagnostics B = 1 + 2I and the convergence-score lemma.
Beyond Navier–Stokes: Adaptation and Arithmetic
| Record | What it establishes | Diagnostics |
|---|---|---|
| Unified Adaptation Theorem v3 | Sequential composition of contractive updates converges under block contraction with a non-negative coupling matrix (Bertsekas–Tsitsiklis) plus LaSalle. PCGrad-type projection preserves first-order descent of each objective but admits non-unique limit points, stalls in one dimension, and period-two orbits at moderate step sizes. | B = 1 + 2I for the balance residual; S ≥ 1 − 2σT/2/(1−σ)2; balanced-matrix eigenvalues; local Hessian of B. |
| Arithmetic Oscillations v2 | Truncated Odlyzko–te Riele sums Σn≤K 2/|ρnζ′(ρn)| with computed ζ′(ρn) over Odlyzko's 2,001,052 zeros cross 1 at K = 58 and 2 near K = 1,040; growth is consistent with the conjectural (log T)5/4 law. | Recovery of the first zeta zeros by Fourier analysis of M(x)/√x. |
Key Results
Each finding is paired with what it clarifies about the regularity problem.
Stretching Drops Out of the Direction Equation
The tangential strain forcing on the vorticity direction, Pξ⊥Sξ, contains no stretching component: α = ξ·Sξ drops out and only the off-diagonal shear rotates the direction. The forcing splits into an alignment part and a Calderón–Zygmund commutator part bounded by the bmo norm of the direction.
The Open Step Is a Bounded-Constant Comparison
The Dirichlet energy of the vorticity direction on a ball controls its mean-square oscillation with the first nonzero Neumann eigenvalue, μN ≈ 4.333: Dr ≥ (μN/r2) Wr. This kinematic Poincaré inequality holds for any field. The direction-oscillation bootstrap closes if and only if C⋆(κ/ν)‖ξ‖bmo ≤ μN/2, an r-independent condition on the core Reynolds number.
Coherence Does Not Stop a Cascade
The Burgers vortex has exactly unidirectional vorticity and sustained stretching α = γ > 0. Aligned vorticity constrains only the self-induced velocity; strain induced by vorticity elsewhere stretches an aligned tube.
No One-Point Direction Statistic Is Necessary for Stretching
A divergence-free field of six Burgers-type patches has an exactly isotropic direction second moment (‖Q‖ = 0, resultant 0) and enstrophy-weighted stretching 0.999. On DNS at the enstrophy peak the mean-strain term is at most 6% of the stretching at the flux scale; the strain–direction covariance carries it (Betchov's identity in explicit form).
A Pointwise Kolmogorov Closure Is Not a Property of Smooth Flows
The hypothesis |ΠK| ≤ αK EK3/2 for every shell fails at t = 0 for a datum with a spectral gap (E10 = 0, Π10 = −1.04), for localised vortex blobs of size ℓ at their own scale by a factor of order ℓ−3 (13.6× to 55× for ℓ = 0.25 to 0.15), and in the exponential dissipation range of generic data (7.5×, 30×, 98× at K = 40, 48, 52).
Uniform Intra-Shell Dispersion Fails for Exact Solutions
For the shear flow u = f(x1)e2 every shell has two modes, so the dispersion constant is nK/2: 9, 31, 175, 625, 2517 at K = 1, 2, 5, 10, 20.
Conditional Active-Shell Enstrophy Bound
Under uniform dispersion, the pointwise closure and cascade contiguity, the enstrophy carried by the active shells satisfies ΩA ≤ C E(0)5/3/ν4/3, and any uniform-in-time enstrophy bound gives regularity by Prodi–Serrin.
Zenodo →Aligned Critical Profiles Are the Axisymmetric Family
For homogeneous critical profiles ω = |x|−2Ψ(x̂), the angular shear G = Pξ⊥Σξ vanishes identically (to 10−15) for swirl-free axisymmetric profiles and for nothing else found: the linearisation at four such profiles has null space equal to the symmetry count (l ≤ 4), next singular value 0.7–1.9% of the largest, and six nonlinear minimisations from random profiles converge to that family. A computed result within band-limited profiles.
Inviscid Truncations Thermalise
Ω/E reaches the lattice mean 〈k2〉 to within 1% at every N from 3 to 8, unchanged under time-step halving. This is the Lee–Kraichnan absolute equilibrium.
Zenodo →Conjectures
Exploratory conjectures from the adaptation-theorem experiments. The status on each tile records the outcome of the experiment reported on that page.
Reproducible Experiments
Every instrument can be run independently. Navier–Stokes instruments use Python (NumPy/SciPy) and C. Shell-model and adaptation experiments use Simplex.
Navier–Stokes Instruments10 experiments
| Experiment | Measures | Result |
|---|---|---|
| verify_review_claims.py · check 1 | Pointwise Kolmogorov closure at t = 0, spectral-gap datum | E10 = 0, Π10 = −1.04: closure fails |
| verify_review_claims.py · check 2 | Intra-shell dispersion for the shear flow u = f(x1)e2 | nK/2 = 9, 31, 175, 625, 2517 at K = 1, 2, 5, 10, 20 |
| verify_review_claims.py · check 4 | Burgers vortex: direction and stretching | Unidirectional vorticity, sustained α = γ > 0 |
| N = 192 closure script | Closure for localised blobs and the dissipation range | 13.6×–55× (ℓ = 0.25–0.15); 7.5×, 30×, 98× at K = 40, 48, 52 |
| Six-patch field | Direction second moment vs enstrophy-weighted stretching | ‖Q‖ = 0, resultant 0, stretching 0.999 |
| E9 · pincer_probe.py | Mean-strain vs covariance share of the stretching on DNS | Mean-strain term ≤ 6% at the flux scale |
| A1 · shear_profile.py | Angular shear G of homogeneous critical profiles | G = 0 to 10−15 on the swirl-free axisymmetric family only |
| pseudospectral_compare.py | FFT solver vs triad kernel v3 | Agree to 2.8×10−17 at N = 4, T = 0.01 |
| experiment_energy_audit.c | Energy conservation of the triad kernel at ν = 0 | Real-coefficient kernel injects up to 15% at N = 7; v3 conserves |
| experiment_leray_angle_variance.c | Angular variance Var(cos2αT) of triads on Z3 | Positive at every shell computed; tends to 4/45 |
Shell-Model Experiments (Simplex)8 experiments
| Experiment | System | Record |
|---|---|---|
| exp_ns_6mode_solved.sx | Cubic shell model, 3 wavenumbers | Finite-window enstrophy trajectories under the stated protocol |
| exp_ns_8mode_solve.sx | Cubic shell model, 4 wavenumbers | Finite-window enstrophy trajectories under the stated protocol |
| exp_ns_12mode_solve.sx | Cubic shell model, 6 wavenumbers | Finite-window enstrophy trajectories under the stated protocol |
| exp_ns_16mode_solve.sx | Cubic shell model, 8 wavenumbers | Finite-window enstrophy trajectories under the stated protocol |
| Mode scaling summary | Shell models, 3–12 wavenumbers | Window-length, threshold and solver dependence of the excursion amplitude |
| exp_ns_scaffold.sx | Multi-level diagnostic H on shell-model trajectories | In-sample classifier; weights learned from reported outcomes |
| exp_ns_viscosity_spectrum.sx | Shell model across viscosity | Trajectory record across ν |
| exp_ns_3d_vortex.sx | First real-coefficient 3D kernel | Mode enumeration and triad structure later carried into the v3 kernel |
Adaptation Theorem9 experiments
| Experiment | Examines | Result |
|---|---|---|
| exp_contraction.sx | Contraction of individual updates (C1) | 5/5 subsystems contract |
| exp_gradient_interference.sx | PCGrad-type projection | Projection identity on 500 conflicting pairs |
| exp_lyapunov.sx | Lyapunov decrease along composed updates | 0 increases observed |
| exp_invariants.sx | Invariant set under strong gradients | 0 violations / 20K steps |
| exp_timescale.sx | Two-timescale separation | Separation maintained |
| exp_interaction_matrix.sx | Block-coupling matrix estimated from trajectories | Converges in 5 cycles |
| exp_convergence_order.sx | Convergence score S | S → 0.000264 |
| exp_iratio_proof.sx | The identity B = 1 + 2I | Holds to 2.2×10−16 for K = 2…20 |
| exp_balance_residual.sx | Gradient flow on the balance residual B | Reaches machine epsilon; local Hessian 2JTJ/D |
Exploratory: Beliefs, Dynamics, Compilers10 experiments
| Experiment | Examines | Result |
|---|---|---|
| exp_anima_deep.sx | Belief interaction and consolidation | Single-run observation |
| exp_anima_correlated.sx | Correlated beliefs with misaligned targets | Exploratory |
| exp_belief_cascade.sx | Chain, circular and delayed topologies | Chain topology recovered |
| exp_skeptical_annealing.sx | Skeptic vs believer across horizons | Conjecture 6.3 refuted |
| exp_chaos_boundary.sx | Convergence score S on the logistic map | Exploratory; periodic and chaotic orbits not separated |
| exp_s_vs_lyapunov.sx | S against the Lyapunov exponent | Complementary diagnostics |
| exp_nash_equilibrium.sx | Iterated Prisoner's Dilemma with an added term | Experiment outcome; not a theorem |
| exp_code_gates.sx | Code-structure convergence via gating | S → 0 at step 50 |
| exp_compiler_passes.sx | Per-program compiler pass adaptation | Per-program adaptation |
| exp_sensitivity.sx | Learning-rate sensitivity | Stable across 3 orders of magnitude |
Run Any Experiment
Experiments use Python, C, or Simplex depending on the domain. Each is standalone and independently reproducible.
Python (Navier–Stokes instruments)
# Requirements
pip install numpy scipy
# The closure, dispersion, scaling and Burgers-vortex checks (seconds)
python3 numerical/verify_review_claims.py
# Angular shear of homogeneous critical profiles (A1)
python3 numerical/shear_profile.py
# Covariance term of the stretching on DNS (E9)
python3 numerical/pincer_probe.py
Code: lab-code
C (Triad kernel v3 and lattice geometry)
# Any C99 compiler
clang -O3 ns-proof-experiments/experiment_leray_angle_variance.c \
solvers/v3/triad_kernel_v3.c -o build/leray_variance -lm
./build/leray_variance
Code: lab-code
Simplex (Shell models and adaptation theorem)
# 1. Compile the Simplex source to LLVM IR
./sxc experiment.sx -o experiment.ll
# 2. Link with the runtime into a native binary
clang -O2 experiment.ll standalone_runtime.c \
-o experiment -lm -lssl -lcrypto
# 3. Run
./experiment
Need the compiler? Download pre-built binaries or build from source.
Full Setup Guide →Citation
@misc{higgins2026energy,
title = {Global Regularity of 3D Navier--Stokes:
An Energy Argument (v5)},
author = {Higgins, Rod},
year = {2026},
doi = {10.5281/zenodo.22710241},
url = {https://lab.senuamedia.com/papers/energy-regularity.html}
}
@misc{higgins2026unified,
title = {Unified Adaptation Theorem: Convergence of Composed
Adaptive Systems via Interaction Matrices and
Higher-Order Convergence Diagnostics (v3)},
author = {Higgins, Rod},
year = {2026},
doi = {10.5281/zenodo.22714409},
url = {https://lab.senuamedia.com/papers/unified-adaptation-theorem.html}
}