Major Jeffrey Camlin, USA (Ret.)
Red Dawn Academic Press | Science Policy Research Unit, University of Sussex
Contributor Role: Author
DOI: 10.63968/post-bio-ai-epistemics.v2n1.013 | Preprint: arXiv:2510.09805
Abstract:
The incompressible Navier–Stokes equations on \(\mathbb{T}^3\) exhibit a structural asymmetry: the spatial domain inherits the compact geometry of \(\mathbb{R}^3/\mathbb{Z}^3\), while the temporal axis remains unbounded and analytically unconstrained. Classical approaches treat time as a neutral parameter, a clock labeling solution states without participating in the analytic structure. On periodic domains, this separation forfeits geometric constraints that the lattice structure naturally provides. We construct a coupled system \((U, \varphi)\) via a bounded vorticity-response functional and prove that the Beale–Kato–Majda and Prodi–Serrin regularity criteria are invariant under this bounded temporal lifting. Numerical validation on grids up to \(256^3\) confirms invariance to machine precision across Reynolds numbers spanning \(10^3\) to \(10^8\).
Keywords: Navier–Stokes equations · BKM criterion · Prodi–Serrin criterion · temporal lifting · vorticity-response functional · regularity invariance · \(\mathbb{T}^3\) torus · bounded diffeomorphism
MSC 2020: 35Q30 · 76D05 · 65M70 · 68T07 · 68T27 · 03D45
ACM: I.2.0
Article Info:
Journal: The Scholarly Journal of Post-Biological Epistemics
Volume: 2 · Issue: 1 · Pages: 1–
ISSN: 3069-499X
Preprint: arXiv:2510.09805
Final manuscript: March 15, 2026
License: CC BY 4.0
Cite this article (APA 7):
Camlin, J. (2026). Invariance of BKM and Prodi–Serrin integrals under bounded temporal lifting. The Scholarly Journal of Post-Biological Epistemics, 2(1). https://doi.org/10.63968/post-bio-ai-epistemics.v2n1.013
Related Works:
· Global Regularity for Navier–Stokes on T³ via Bounded Vorticity–Response Functionals (subsequent work)
· iDNS: Numerical Validation of 3D Navier–Stokes on T³ (numerical confirmation)
Contents
1. Introduction
1.1 The Construction
1.2 Bounded Temporal Lifting
2. Preliminaries
2.1 Function Spaces and Navier–Stokes Equations
2.2 Bounded Temporal Lifting
3. Main Theorem
4. Numerical Validation
5. Example Algorithm
References
The incompressible Navier–Stokes equations on \(\mathbb{T}^3 = \mathbb{R}^3/\mathbb{Z}^3\) inherit compact geometry from the spatial lattice, yet classical formulations treat time as an unbounded external parameter. This asymmetry is artificial: the lattice structure \(\mathbb{Z}^3\) induces geometric constraints — symmetry axes, fundamental domains, spectral localization — that the temporal direction does not share.
We embed time into this geometry via bounded temporal lifting: a coupled construction \((U, \varphi)\) where the velocity field and temporal diffeomorphism are solved simultaneously through a bounded vorticity-response functional \(\Phi\). Physical time is recovered as \(t = \varphi(\tau)\).
The main results establish that the Beale–Kato–Majda integral and Prodi–Serrin criteria are invariant under this lifting.
Let \(\Phi : \mathbb{R}_{\geq 0} \to [\varphi_{\min}, \varphi_{\max}]\) be a bounded vorticity-response functional with \(0 < \varphi_{\min} \leq \varphi_{\max} < \infty\). Bounded Temporal Lifting constructs a coupled system \((U, \varphi)\) satisfying
\[ \frac{dt}{d\tau} = \varphi'(\tau) = \Phi\bigl(\|\nabla \times U(\tau)\|_{L^\infty}\bigr), \qquad \varphi(0) = 0. \]The velocity field \(U\) and diffeomorphism \(\varphi\) evolve simultaneously in the auxiliary parameter \(\tau\). Physical time is recovered as \(t = \varphi(\tau)\).
The bounds on \(\Phi\) ensure \(\varphi : [0,\infty) \to [0,\infty)\) is a \(C^1\) diffeomorphism. Near high-vorticity regions, \(\varphi'\) remains bounded below, preventing the auxiliary coordinate from stalling. The construction is motivated by the Path Lifting Lemma: a loop on \(S^1\) lifts to a smooth path on the universal cover \(\mathbb{R}\). Here, the periodic structure of \(\mathbb{T}^3\) provides the covering geometry.
The construction proceeds in three steps:
A bounded vorticity-response functional \(\Phi : \mathbb{R}_{\geq 0} \to [\varphi_{\min}, \varphi_{\max}]\) couples the velocity field \(U\) to the temporal diffeomorphism \(\varphi\):
\[ U(x,\tau) = u(x,\varphi(\tau)), \qquad \varphi'(\tau) = \Phi\bigl(\|\nabla \times U(\tau)\|_{L^\infty}\bigr), \qquad \varphi(0) = 0. \]Physical time is recovered as \(t = \varphi(\tau)\). The bounds \(0 < \varphi_{\min} \leq \varphi' \leq \varphi_{\max}\) ensure \(\varphi\) is a \(C^1\) diffeomorphism of \([0,\infty)\) onto itself.
Figure 1. Geometry of bounded temporal lifting. Blue: Fundamental cell \([0,1)^3\). Red: Axis \(\mathcal{A}\) at centroid \(x_c = 1/2\). Green: Lifted trajectory \(U(\mathbf{x}, \tau)\). Magenta: Spectral projection via \(\varphi(\tau)\). The fibers \((z, \tau)\) through \(\mathcal{A}\) embed time into the torus.
This section establishes the theoretical foundation for the bounded temporal lifting framework.
Let \(\mathbb{T}^3 := \mathbb{R}^3 / \mathbb{Z}^3\) denote the three-torus. We use standard Lebesgue spaces \(L^p(\mathbb{T}^3)\) and Sobolev spaces \(H^s(\mathbb{T}^3)\) for \(s \ge 0\). The divergence-free subspace is defined by
\[ H^s_{\mathrm{div}}(\mathbb{T}^3) := \{\, u \in H^s(\mathbb{T}^3)^3 : \nabla \cdot u = 0 \,\}. \]The incompressible Navier–Stokes equations on \(\mathbb{T}^3\) are:
\[ \partial_t u + (u \cdot \nabla)u + \nabla p - \nu \Delta u = 0, \qquad \nabla \cdot u = 0, \]for velocity \(u(x,t) \in \mathbb{R}^3\), pressure \(p(x,t) \in \mathbb{R}\), viscosity \(\nu > 0\), and initial data \(u(x,0) = u_0(x) \in H^s_{\mathrm{div}}(\mathbb{T}^3)\) with \(s\) sufficiently large. We follow the classical framework of Leray (1934) and Hopf (1951).
Let \(\varphi \in C^\infty([0,\infty))\) with \(\varphi' > 0\). Define the lifted trajectory by
\[ U(x,\tau) := u(x,\varphi(\tau)), \qquad t = \varphi(\tau). \]Unlike classical time reparametrization — a neutral coordinate change — bounded temporal lifting is chosen adaptively to smooth derivative discontinuities at singular times and restore global \(C^\infty\) regularity. On the periodic domain \(\mathbb{T}^3 = \mathbb{R}^3/\mathbb{Z}^3\), the lifting geometry is not arbitrary: natural choices align the computational trajectory with the centroid of the fundamental domain, exploiting lattice symmetries for enhanced stability.
Theorem 3.1 (Temporal Lift Equivalence Theorem).
Let \(u(x,t)\) be a Leray–Hopf (resp. classical) solution of the incompressible Navier–Stokes equations on \(\mathbb{T}^3 = \mathbb{R}^3/\mathbb{Z}^3\) with initial data \(u_0(x) \in H^s_{\mathrm{div}}(\mathbb{T}^3)\).
Let \(\varphi \in C^\infty(\mathbb{R})\) be strictly increasing with \(0 < c \leq \varphi'(\tau) \leq C < \infty\). Define the lifted solution by
\[ U(x,\tau) := u(x,\varphi(\tau)), \qquad P(x,\tau) := p(x,\varphi(\tau)). \]Then \(U\) is a Leray–Hopf (resp. classical) solution of the lifted Navier–Stokes system
\[ \varphi'(\tau)\,\partial_\tau U + (U\cdot\nabla)U + \nabla P - \nu \Delta U = 0, \qquad \nabla \cdot U = 0, \]which preserves the Leray–Hopf energy structure and all regularity criteria up to constants depending only on \(c\) and \(C\).
In particular, the Prodi–Serrin and Beale–Kato–Majda blowup criteria remain invariant under such lifts. If \(\varphi'\) is allowed to vanish or blow up, singularities may be shifted to infinite lifted time \(\tau\), but the system then leaves the class of uniformly parabolic Navier–Stokes equations.
Proof.
The proof proceeds by a change of variables in the weak formulation. Let \(\psi \in C_c^{\infty}(\mathbb{T}^3 \times [0,T))^3\) satisfy \(\nabla \cdot \psi = 0\).
For \(u(x,t)\) a Leray–Hopf solution, the weak form is
\[ \int_0^T \!\!\int_{\mathbb{T}^3} \Bigl( u \cdot \partial_t \psi + (u\cdot\nabla)u \cdot \psi + \nu \nabla u : \nabla \psi \Bigr) \,dx\,dt = 0. \]Substitute \(t = \varphi(\tau)\) and define \(\tilde{\psi}(x,\tau) = \psi(x,\varphi(\tau))\). Since \(dt = \varphi'(\tau)\, d\tau\) and \(\partial_t \psi = \varphi'(\tau)\, \partial_\tau \tilde{\psi}\) by the chain rule, integration yields
\[ \int_0^{\tilde{T}} \!\!\int_{\mathbb{T}^3} \Bigl( U \cdot (\varphi'(\tau)\,\partial_\tau \tilde{\psi}) + (U\cdot\nabla)U \cdot \tilde{\psi} + \nu \nabla U : \nabla \tilde{\psi} \Bigr) \,dx\,d\tau = 0, \]which is precisely the weak form of the lifted system.
For the energy inequality, the same substitution gives
\[ \frac{1}{2}\|U(\tau)\|_{L^2}^2 + \nu \int_0^\tau \|\nabla U(s)\|_{L^2}^2 \,\varphi'(s)\,ds \leq \frac{1}{2}\|U(0)\|_{L^2}^2, \]preserving the Leray–Hopf structure with \(\varphi'(s)\) entering as a time weight.
Regularity criteria depending on \(L^p_t L^q_x\) norms are preserved:
\[ \int_0^{\tilde{T}} \|U\|_{L^q}^p \,\varphi'(\tau)\, d\tau = \int_0^T \|u\|_{L^q}^p\, dt. \]Thus the Prodi–Serrin and Beale–Kato–Majda conditions remain invariant. \(\square\)
Remark 3.2 (Geometric Interpretation of the Regularity Threshold).
The Sobolev embedding theorem establishes that classical solutions to Navier–Stokes on \(\mathbb{T}^3\) require \(u \in H^s\) with \(s > \tfrac{5}{2}\). This critical exponent admits a geometric decomposition:
\[ \frac{5}{2} = \underbrace{2}_{\dim(\mathbb{T}^2)} + \underbrace{\frac{1}{2}}_{\text{centroid coordinate}}, \]where \(\dim(\mathbb{T}^2) = 2\) is the dimension of spatial slices on which Navier–Stokes is globally regular, and \(\tfrac{1}{2}\) is the coordinate of the geometric centroid in the fundamental domain \([0,1)^3\).
This decomposition suggests that the "missing half-derivative" separating weak (Leray–Hopf) solutions from classical solutions may be related to the symmetry structure of the periodic domain. Bounded temporal lifting, by aligning the computational trajectory with this centroid, exploits lattice symmetries that merit further investigation in the context of spectral methods.
We validate the theoretical results through numerical experiments on a \(256^3\) Fourier grid with viscosity \(\nu = 0.01\) and Taylor–Green initial data. Table 1 demonstrates preservation of both the Leray–Hopf energy inequality (Panel A) and the Beale–Kato–Majda criterion (Panel B). Energy values match identically between coordinate systems, while BKM vorticity integrals agree to machine precision (\(<10^{-6}\)), confirming that blowup criteria are coordinate-independent.
In extended validation experiments spanning Reynolds numbers from \(10^3\) to \(10^8\), the BKM integral converges to approximately \(36.9\) across all regimes, suggesting the presence of a geometric invariant arising from the torus structure.
| Panel A: Energy Conservation | |||||
| Physical time | Lifted time | ||||
| \(t\) | \(\|u\|_{L^2}^2\) | \(\int\|\nabla u\|^2\) | \(\tau\) | \(\|U\|_{L^2}^2\) | \(\int\|\nabla U\|^2\varphi'\) |
| 5 | 1.229 | 0.243 | 10 | 1.229 | 0.243 |
| 10 | 1.205 | 0.491 | 20 | 1.205 | 0.491 |
| 15 | 1.178 | 0.734 | 30 | 1.178 | 0.734 |
| 20 | 1.149 | 0.972 | 40 | 1.149 | 0.972 |
| 25 | 1.122 | 1.206 | 50 | 1.122 | 1.206 |
| Panel B: Beale–Kato–Majda Criterion | |||||
| \(t\) | \(\int\|\omega\|_{L^\infty}\) | \(\tau\) | \(\int\|\Omega\|_{L^\infty}\varphi'\) | \(|\text{Diff}|\) | |
| 5.0 | 2.76 | 10.2 | 2.76 | \(8.3\times10^{-7}\) | |
| 10.0 | 5.63 | 18.7 | 5.63 | \(1.2\times10^{-7}\) | |
| 15.0 | 8.54 | 25.3 | 8.54 | \(2.9\times10^{-7}\) | |
| 20.0 | 11.49 | 31.1 | 11.49 | \(4.7\times10^{-7}\) | |
| 25.0 | 14.47 | 36.4 | 14.47 | \(6.1\times10^{-7}\) | |
Table 1. Numerical validation of theorem preservation properties. Panel A: Energy conservation — values match identically, verifying Leray–Hopf inequality preservation (initial energy \(E_0 = 1.250\)). Panel B: BKM criterion — vorticity integrals agree to machine precision, confirming blowup condition invariance. Note: the lifted integral uses weight \(\varphi'\) (corrected from \(\varphi'^{-1}\) in v1).
Algorithm 1: Adaptive Bounded Temporal Lifting Procedure on \(\mathbb{T}^3\)
Require: Initial velocity field \(u_0(x)\) on \(\mathbb{T}^3\), viscosity \(\nu\), time step \(\Delta t\), total time \(T\)
Ensure: Lifted trajectory \(U(x,\tau)\) and temporal map \(\varphi(\tau)\)
1: Initialize \(\tau \gets 0\), \(u(x,0) \gets u_0(x)\)
2: for \(t \gets 0\) to \(T\) step \(\Delta t\) do
3: Compute \(\varphi'(t) \gets f(\|\nabla u(x,t)\|)\)
4: Update lifted time: \(\tau \gets \tau + \varphi'(t)\Delta t\)
5: Set \(U(x,\tau) \gets u(x,t)\)
6: Integrate lifted system: \(\varphi'(\tau)\,\partial_\tau U + (U\!\cdot\!\nabla)U + \nabla P - \nu\Delta U = 0\)
7: end for
8: return \(U(x,\tau)\) and \(\varphi(\tau)\)
Funding. This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.
Declaration of AI Instrumentality. During the preparation of this manuscript, the author utilized Artificial Intelligence (Red Dawn MLAR Architecture, 2026) as a research instrument for mathematical verification and exposition refinement. The author reviewed and verified all content and takes full responsibility for the published text.
[1] R. A. Adams, J. J. F. Fournier. Sobolev Spaces, 2nd Edition. Academic Press, 2003.
[2] J. T. Beale, T. Kato, A. Majda. Remarks on the breakdown of smooth solutions for the 3D Euler equations. Commun. Math. Phys., 94:61–66, 1984. doi:10.1007/BF01212349
[3] L. C. Evans. Partial Differential Equations, 2nd Edition. AMS, 2010.
[4] A. Hatcher. Algebraic Topology. Cambridge Univ. Press, 2002. URL
[5] E. Hopf. Über die Anfangswertaufgabe für die hydrodynamischen Grundgleichungen. Math. Nachr., 4:213–231, 1951.
[6] J. Leray. Sur le mouvement d'un liquide visqueux emplissant l'espace. Acta Math., 63:193–248, 1934.
[7] G. Prodi. Un teorema di unicità per le equazioni di Navier–Stokes. Ann. Sc. Norm. Super. Pisa, 13:429–435, 1959.
[8] J. Serrin. On the interior regularity of weak solutions of the Navier–Stokes equations. Arch. Ration. Mech. Anal., 9:187–195, 1962. doi:10.1007/BF00253344
[9] R. Temam. Navier–Stokes Equations: Theory and Numerical Analysis. North-Holland, 1977 (2001 reprint).
— End of Article —