Senuamedia Lab
Navier–Stokes Regularity Research
The open step, stated exactly
C⋆(κ/ν) · ‖ξ‖bmo ≤ μN/2
μN ≈ 4.333 is the first Neumann eigenvalue of the ball. The direction-oscillation bootstrap closes if and only if this bounded-constant comparison on the core Reynolds number holds. It is not a scaling gap.

Navier–Stokes Regularity: Where the Difficulty Sits

A first singularity of 3D Navier–Stokes is Type I (blow-up at the scaling rate) or Type II (faster). Zooming in on a Type I singularity gives an ancient Navier–Stokes solution with bounded scale-invariant energy; on a gauged Type II singularity, Seregin's modulated zoom gives an ancient Euler solution. Global regularity is therefore equivalent to Liouville-type statements for these two limit objects. This programme locates, with explicit fields and computed identities, which strategies toward those statements are closed and where the remaining difficulty sits.

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

  1. Liouville for ancient Type I Navier–Stokes solutions.
  2. Exclusion of a rotating |x|−2 core as an ancient Type I solution.
  3. 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

InstrumentMeasuresValidation
Triad kernel v3 (C)3D Galerkin truncation, exact Fourier trilinear coupling, Leray projectionEnergy 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 DNSBetchov's identity in explicit form
verify_review_claims.pySpectral-closure, intra-shell dispersion, amplitude scaling, and Burgers-vortex checksSelf-contained; NumPy only; seconds
Thermalisation checkΩ/E against the lattice mean ⟨k2⟩ in inviscid truncationsWithin 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.

Energy Argument · v5

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 / Effective PDE · v2

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, Cascade Stabilisation · v19

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.

Arithmetic Oscillations · v2

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.

Three Fluid Equations · v2

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.

Holistic Scaffold Framework · v3

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 · v3

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

RecordWhat it establishesDiagnostics
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.

Unified Adaptation Theorem →   Arithmetic Oscillations →

Key Results

Each finding is paired with what it clarifies about the regularity problem.

Direction identities · proved

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.

What it clarifies: direction dynamics and amplitude dynamics separate; the direction is driven by shear alone.

Coercivity · proved

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.

What it clarifies: the remaining step is a constant comparison, not a scaling gap.

Explicit field · verified

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.

What it clarifies: Liouville theorems based on vorticity coherence are closed as a strategy class.

Explicit field · computed

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).

What it clarifies: the difficulty sits in the strain–direction covariance at the flux scale.

Spectral closure · computed

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).

What it clarifies: such a closure excludes concentration a priori and cannot be a hypothesis of a regularity theorem.

Exact solution · verified

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.

What it clarifies: uniform dispersion cannot be assumed for all smooth solutions; it delimits where the conditional bound applies.

Energy Argument v5 · conditional

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.

What it clarifies: exactly which cascade hypotheses would suffice, and, with the two results above, where they hold.

Zenodo →
Critical profiles · computed

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.

What it clarifies: Grujić's alignment condition holds for the leading profile only on the axisymmetric family; the residual is a non-symmetric core whose direction is driven at no less than about 1% of κ/|x|2 per unit time.

Truncated Euler · measured

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.

What it clarifies: the enstrophy growth observed in truncated Euler is equilibration, not evidence about the PDE.

Zenodo →

Conjectures

Exploratory conjectures from the adaptation-theorem experiments. The status on each tile records the outcome of the experiment reported on that page.

Details →

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
ExperimentMeasuresResult
verify_review_claims.py · check 1Pointwise Kolmogorov closure at t = 0, spectral-gap datumE10 = 0, Π10 = −1.04: closure fails
verify_review_claims.py · check 2Intra-shell dispersion for the shear flow u = f(x1)e2nK/2 = 9, 31, 175, 625, 2517 at K = 1, 2, 5, 10, 20
verify_review_claims.py · check 4Burgers vortex: direction and stretchingUnidirectional vorticity, sustained α = γ > 0
N = 192 closure scriptClosure for localised blobs and the dissipation range13.6×–55× (ℓ = 0.25–0.15); 7.5×, 30×, 98× at K = 40, 48, 52
Six-patch fieldDirection second moment vs enstrophy-weighted stretching‖Q‖ = 0, resultant 0, stretching 0.999
E9 · pincer_probe.pyMean-strain vs covariance share of the stretching on DNSMean-strain term ≤ 6% at the flux scale
A1 · shear_profile.pyAngular shear G of homogeneous critical profilesG = 0 to 10−15 on the swirl-free axisymmetric family only
pseudospectral_compare.pyFFT solver vs triad kernel v3Agree to 2.8×10−17 at N = 4, T = 0.01
experiment_energy_audit.cEnergy conservation of the triad kernel at ν = 0Real-coefficient kernel injects up to 15% at N = 7; v3 conserves
experiment_leray_angle_variance.cAngular variance Var(cos2αT) of triads on Z3Positive at every shell computed; tends to 4/45
▶Shell-Model Experiments (Simplex)8 experiments
ExperimentSystemRecord
exp_ns_6mode_solved.sxCubic shell model, 3 wavenumbersFinite-window enstrophy trajectories under the stated protocol
exp_ns_8mode_solve.sxCubic shell model, 4 wavenumbersFinite-window enstrophy trajectories under the stated protocol
exp_ns_12mode_solve.sxCubic shell model, 6 wavenumbersFinite-window enstrophy trajectories under the stated protocol
exp_ns_16mode_solve.sxCubic shell model, 8 wavenumbersFinite-window enstrophy trajectories under the stated protocol
Mode scaling summaryShell models, 3–12 wavenumbersWindow-length, threshold and solver dependence of the excursion amplitude
exp_ns_scaffold.sxMulti-level diagnostic H on shell-model trajectoriesIn-sample classifier; weights learned from reported outcomes
exp_ns_viscosity_spectrum.sxShell model across viscosityTrajectory record across ν
exp_ns_3d_vortex.sxFirst real-coefficient 3D kernelMode enumeration and triad structure later carried into the v3 kernel
▶Adaptation Theorem9 experiments
ExperimentExaminesResult
exp_contraction.sxContraction of individual updates (C1)5/5 subsystems contract
exp_gradient_interference.sxPCGrad-type projectionProjection identity on 500 conflicting pairs
exp_lyapunov.sxLyapunov decrease along composed updates0 increases observed
exp_invariants.sxInvariant set under strong gradients0 violations / 20K steps
exp_timescale.sxTwo-timescale separationSeparation maintained
exp_interaction_matrix.sxBlock-coupling matrix estimated from trajectoriesConverges in 5 cycles
exp_convergence_order.sxConvergence score SS → 0.000264
exp_iratio_proof.sxThe identity B = 1 + 2IHolds to 2.2×10−16 for K = 2…20
exp_balance_residual.sxGradient flow on the balance residual BReaches machine epsilon; local Hessian 2JTJ/D
▶Exploratory: Beliefs, Dynamics, Compilers10 experiments
ExperimentExaminesResult
exp_anima_deep.sxBelief interaction and consolidationSingle-run observation
exp_anima_correlated.sxCorrelated beliefs with misaligned targetsExploratory
exp_belief_cascade.sxChain, circular and delayed topologiesChain topology recovered
exp_skeptical_annealing.sxSkeptic vs believer across horizonsConjecture 6.3 refuted
exp_chaos_boundary.sxConvergence score S on the logistic mapExploratory; periodic and chaotic orbits not separated
exp_s_vs_lyapunov.sxS against the Lyapunov exponentComplementary diagnostics
exp_nash_equilibrium.sxIterated Prisoner's Dilemma with an added termExperiment outcome; not a theorem
exp_code_gates.sxCode-structure convergence via gatingS → 0 at step 50
exp_compiler_passes.sxPer-program compiler pass adaptationPer-program adaptation
exp_sensitivity.sxLearning-rate sensitivityStable across 3 orders of magnitude
All Experiments →

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}
}