Renormalized Lambert–W Cascade and Finite-Time Amplification and Blowup for moving coterminal reconstructed b Dynamics on T 3 for the 3D Navier Stokes Equations

Terry Moschandreou
Terry Moschandreou * § Doctor of Philosophy Applied Mathematics
§ Intermediate Science and Mathematics

Send Message

To: Author

Renormalized Lambert–W Cascade and Finite-Time Amplification and Blowup for moving coterminal reconstructed b Dynamics on T 3 for the 3D Navier Stokes Equations

Article Fingerprint

ReserarchID

V89U5

Renormalized Lambert–W Cascade and Finite-Time Amplification and Blowup for moving coterminal reconstructed b Dynamics on T 3 for the 3D Navier Stokes Equations Banner

Key Research Insights

Synthesized scholarly intelligence & interactive research assistant
  • English
  • Afrikaans
  • Albanian
  • Amharic
  • Arabic
  • Armenian
  • Azerbaijani
  • Basque
  • Belarusian
  • Bengali
  • Bosnian
  • Bulgarian
  • Catalan
  • Cebuano
  • Chichewa
  • Chinese (Simplified)
  • Chinese (Traditional)
  • Corsican
  • Croatian
  • Czech
  • Danish
  • Dutch
  • Esperanto
  • Estonian
  • Filipino
  • Finnish
  • French
  • Frisian
  • Galician
  • Georgian
  • German
  • Greek
  • Gujarati
  • Haitian Creole
  • Hausa
  • Hawaiian
  • Hebrew
  • Hindi
  • Hmong
  • Hungarian
  • Icelandic
  • Igbo
  • Indonesian
  • Irish
  • Italian
  • Japanese
  • Javanese
  • Kannada
  • Kazakh
  • Khmer
  • Korean
  • Kurdish (Kurmanji)
  • Kyrgyz
  • Lao
  • Latin
  • Latvian
  • Lithuanian
  • Luxembourgish
  • Macedonian
  • Malagasy
  • Malay
  • Malayalam
  • Maltese
  • Maori
  • Marathi
  • Mongolian
  • Myanmar (Burmese)
  • Nepali
  • Norwegian
  • Pashto
  • Persian
  • Polish
  • Portuguese
  • Punjabi
  • Romanian
  • Russian
  • Samoan
  • Scots Gaelic
  • Serbian
  • Sesotho
  • Shona
  • Sindhi
  • Sinhala
  • Slovak
  • Slovenian
  • Somali
  • Spanish
  • Sundanese
  • Swahili
  • Swedish
  • Tajik
  • Tamil
  • Telugu
  • Thai
  • Turkish
  • Ukrainian
  • Urdu
  • Uzbek
  • Vietnamese
  • Welsh
  • Xhosa
  • Yiddish
  • Yoruba
  • Zulu
Reading Preferences
Font Size
Line Spacing
Background
This converted HTML version may contain rendering inconsistencies. Please refer to the PDF for the authoritative version, or click here to provide feedback.

Abstract

This article extracts and consolidates the renormalized Lambert–W branch-point cascade, its distinguished phase reduction, the exact characteristic invariant and finite-time amplification mechanism, and the extended reconstructed bi equation on T3 . For a moving coterminal approach the principal proofs and terminal reconstruction estimates are demonstrated. The presentation separates exact finite-depth statements from coupled-depth asymptotics and records the hypotheses required for the extended PDE reconstruction. This paper further supports a recent paper [1] published by the corresponding author which claims that the Navier Stokes equations lose smoothness in finite time from initial smooth data.

1. Introduction and organization

Papers explicitly stating and showing that there are finite time singularities of the 3D incompressible Navier Stokes equations were not found in my literature search at the time of the published work in [1] and before it’s inception. There is the famous result of Terence Tao in [3] for a finite time blowup of the fluid equations but for an averaged NS equation. Going back to 2021 Moschandreou [4] showed that there are no-finite time singularities but left open the possibility of finite time blowup on sets of measure zero in T3(see also [5]). Singular solution behaviour was determined for the Euler and related partial differential equations in [6], [7] and [8]. Recently the existence of a nonsmooth periodic attractor is shown in [9]. Critical and important findings among papers in the theory of Navier Stokes equations span the years (1934-1984) and are [10],[11], [12] and [13]. Important references in differential equations theory are found in [14],[15],[16],[17] and a strong PDE source textbook is found in [18]. Next a solid coverage of vorticity and Incompressible Navier Stokes theory is found in [19]. Excellent sources on the Navier Stokes problem can be found in [20], [21],[22], [23] and [24]. A compendium on the LambertW function and it’s theory is found in [25]. And finally a very good reference in functional analysis and in particular on the Fredholm alternative and Fredholm compatability theory as used can be located here[26] The analysis is organized around four linked structures: the distinguished phase θ=t−z; the real renormalized Lambert map Δ↦1+W0(−e−1−Δ); the characteristic invariant that transfers branch-distance collapse to amplification of b3; and the corrected extended b3 equation with periodic velocity reconstruction. The finite-depth cascade and finite-time characteristic proofs are presented first. The later sections retain the essential full-b3 PDE reconstruction, comparison, trapping, and terminal closure arguments while suppressing duplicate derivations and exploratory restatements.

2. The Renormalized Lambert Cascade, Geometric Reconstruction, and Finite-Time Amplification Mechanism

The purpose of this paper is to present in condensed form the principal mathematical architecture developed first in a longer construction which I used to simplify with full import to the present paper. For example some steps were implicit but have been worked out thoroughly filling in completely any gaps in the flow of the reading of this particular work. The central mechanism begins with a renormalized Lambert–W cascade, passes through a distinguished scalar product variable

b3:=uxuy,

and culminates in a characteristic regime in which the Lambert-generated logarithmic coefficient produces a positive quadratic amplification term. As stated in the conclusion of [1] this leads to:
”In particular, the analysis shows that singular behavior is associated with the branch structure of the Lambert W function rather than unbounded growth of the solution itself at time t1<t2. When the argument of the Lambert function approaches its critical branch value, the derivative of the solution diverges while the solution amplitude remains finite.” [1]
Further in the conclusion there it states explicitly: ”However there is a ‖b‖2 term which leads to blowup in derivatives as shown at t=t1 and for t=t2(t2>t1) it was proven analytically and numerically that there is finite time amplitude blowup in the solution b3”.

The analysis is carried out on the periodic domain

T3=(R/2πZ)3,

with velocity

u=(ux,uy,uz).

A central point of the construction is that b3 is not introduced merely as an auxiliary scalar. It remains the exact physical product

b3=uxuy

in the reconstructed velocity field. The scalar characteristic analysis and the geometric velocity reconstruction are therefore linked through an identity rather than through an asymptotic identification.

The reduced paper retains five principal structures:

renormalized Lambert--W cascade⇓branch variable QN and logarithmic coefficient ΓN⇓full extended b3-PDE⇓terminal characteristic trapping and finite-time amplification⇓periodic square-root velocity reconstruction on T3.

The purpose of this introduction is to display this architecture before the proofs are developed in detail.

2.1. Historical symmetric precursor and the active reconstruction

An earlier symmetric stage of the construction used a common transverse profile

uy=uz=Q.

At that level,

b1=uyuz=Q2,b2=uxuz=uxQ,b3=uxuy=uxQ,

so that

b2=b3.

This symmetric representation remains useful for understanding the origin of several scalar formulas and preterminal periodic constructions. However, it is not the active terminal reconstruction used in the final analysis.

The active variables are instead

b=b3=uxuy,A=uy,C=uz,ux=bA.

Thus

uxuy=b

holds identically, whereas uy and uz are no longer assumed to coincide.

The two transverse components acquire different roles:

uy=Ais determined by the square-root geometry,uz=Cis reconstructed from exact incompressibility.

This distinction is fundamental to the presentation.

2.2. Renormalized Lambert–W cascade

The scalar branch variable is generated by the principal real branch W0. Let

ΔN,0>0

denote the initial branch distance and define

ΔN,k+1=1+W0(−e−1−ΔN,k).

The terminal Lambert regime corresponds to

ΔN,k↓0.

The local branch-point expansion is

1+W0(−e−1−d)=2d−23d+O(d3/2),d↓0.

Consequently, at every fixed finite depth k,

ΔN,k≍ΔN,02−k.

The distinguished terminal quantity will be denoted

QN:=ΔN,k.

Its logarithmic derivative is

ΓN:=∂ξlog⁡QN=QN′QN.

When the branch distance is represented by a terminal coordinate s for which

QN(s)≍s2−k,

one obtains

ΓN(s)∼2−ks.

Thus the renormalized cascade produces two simultaneous effects:

QN⟶0,ΓN⟶+∞.

At the frequently used finite depth k=3,

QN≍s1/8,ΓN∼18s.

It is essential to distinguish the branch variable QN from the physical transverse velocity components. In the active reconstruction,

QN≠uyandQN≠uzin general.

Rather,

QN: Lambert branch variable,A,C: reconstructed physical velocity components.

The role of QN is to generate ΓN, which enters the b3-equation. The physical velocity behavior is determined independently by the geometric reconstruction.

2.3. Distinguished phase θ=t−z

The terminal analysis is organized around

θ=t−z.

After locating the distinguished branch position θ∗,N, define

s=t−z−θ∗,N.

Along the reduced characteristic for which

z˙=b,

one obtains

s˙=1−b.

This elementary identity becomes the trapping law.

If the terminal region is entered with

b(se)≥1+δ,δ>0,

then

s˙=1−b≤−δ<0.

Hence

s(t)↓0.

Thus the characteristic is driven toward the Lambert endpoint once the amplitude has crossed the threshold b>1.

2.4. The corrected extended b3-equation

With

A=uy,C=uz,ux=bA,

define the horizontal differential operator

X=A∂x+bA∂y

and the transport operator

D=∂t+b∂z+C2X.

The corrected extended amplitude equation can be written

Db=ΓNb2+2CG−RN,

where

RN=2bXC+bCt+CzC.

Equivalently, before the final compression of the reconstruction terms, one may write

DNbN+EN=QN′QNbN2+2CNG.

Along a characteristic,

DNbN=b˙N,

and therefore

b˙N=QN′QNbN2+2CNG−EN.

The additional terms are not discarded. Their role is measured relative to the positive source.

Define

SN=QN′QNbN2+2CNG.

If

|EN|≤ϑSN,0<ϑ<1,

then

b˙N=SN−EN≥SN−|EN|≥(1−ϑ)SN.

Hence

b˙N≥(1−ϑ)[QN′QNbN2+2CNG].

On the positive branch,

CN>0,G≥0,

so

2CNG≥0

and consequently

b˙N≥(1−ϑ)QN′QNbN2.

This is the Riccati-type inequality that drives the terminal growth.

The logical reduction is therefore

full extended b3-PDE⇓b˙N=QN′QNbN2+2CNG−EN⇓|EN|≤ϑ(positive source)b˙N≥(1−ϑ)[QN′QNbN2+2CNG]⇓CNG≥0b˙N≥(1−ϑ)QN′QNbN2.

Thus the reconstruction error is used explicitly in obtaining the quadratic lower bound.

2.5. The explicit reconstruction error

At an intermediate stage of the periodic reconstruction, the additional terms can be collected as

EN=bNCN,t+CN,zCN+bNANCN,x+bN2ANCN,y+2CN∇bN⋅∇CN.

The terminal objective is

EN=o(QN′QNbN2).

Equivalently,

|EN|(QN′/QN)bN2⟶0.

A normalized estimate of the form

|EN|LNbN2≤K[|CN,t+CN,z|LNbN+|CN,x|LNbN+|CN,y|LN+|∇bN||∇CN|LNbN2],

with

LN=QN′QN,

reduces the problem to weighted reconstruction bounds.

When

|CN,t|+|∇CN|≤MC

and

QN′QN→∞,bN→∞,

the first three normalized terms vanish.

For the remaining term, the terminal derivative estimate

|∇bN|≤CbpNQN′QNbN

gives

|∇bN||∇CN|(QN′/QN)bN2≤CbMCpNbN.

The comparison argument provides

bNpN→∞,

and hence

pNbN→0.

Therefore

|EN|(QN′/QN)bN2→0.

Consequently, for every fixed

0<ϑ<1,

one eventually has

|EN|≤ϑQN′QNbN2,

and hence

b˙N≥(1−ϑ)QN′QNbN2>0.

This estimate is the bridge from the full extended PDE to the scalar terminal comparison.

2.6. Release-layer and square-root-core estimates

The terminal reconstruction naturally separates into a release layer and a deeper square-root core.

In the release layer,

|Rrel|ΓNb2=O(s7/4).

In the square-root core,

|Rcore|ΓNb2=O(s).

Both estimates vanish as

s↓0.

Thus

|RN|ΓNb2→0,

and the geometric transition into the terminal square-root regime does not destroy the positive Lambert-generated quadratic term.

Rather, the error becomes progressively smaller relative to that term.

2.7. Finite-depth trapping and amplification

The two fundamental terminal inequalities are

s˙=1−b

and

b˙≥(1−ϑ)ΓNb2.

For finite cascade depth k,

ΓN∼2−ks.

Since

−s˙=b−1,

one obtains

−dbds=b˙b−1≥(1−ϑ)2−ksb2b−1.

For b>1,

b2b−1>b,

and therefore

−dbds>(1−ϑ)2−kbs.

Integration gives

b(s)≥be(ses)(1−ϑ)2−k.

Consequently,

s↓0⟹b(s)→+∞.

At k=3,

b(s)≥be(ses)(1−ϑ)/8.

The remaining physical characteristic time is finite because

dt=−dsb−1.

The lower power-law bound gives

1b−1≲sp,p>0,

and hence

∫0sedsb−1<∞.

Thus

T∗,N<∞,

with

s(t)→0,QN(t)→0,b3(t)→+∞as t↑T∗,N.

The same result can be expressed directly in terms of QN. In the large-amplitude terminal regime,

−dbNdQN≥pNbNQN,pN>0.

Integration from a terminal-entry point

(Qe,N,be,N)

gives

bN(QN)≥be,N(Qe,NQN)pN.

Hence

QN↓0⟹bN→+∞.

The feedback mechanism is therefore

QN↓0⇒ΓN↑⇒bN↑⇒s˙=1−bN<0⇒s↓0⇒QN↓0.

2.8. Square-root velocity geometry

The terminal geometric specialization begins with

2xux+uy=0,y=x2.

Using

ux=bA,uy=A,

gives

2xbA+A=0,

and therefore

A2=−2xb.

Thus

uy=A=−2xb,ux=bA

on the selected real branch.

Moreover,

uxuy=bA2=−12x.

Since

y=x2,

one obtains

|uxuy|=12|x|=12y.

Thus the square-root behavior follows directly from the geometric constraint and the exact product identity b=uxuy.

2.9. Periodic realization through χ(x)

On T3, the local coordinate x is replaced near a selected periodic zero by a smooth periodic function χ(x). The active relation becomes

2χ(x)ux+uy=0

or equivalently

A2=−2χ(x)b.

Therefore

uy2=−2χb,ux2=−b2χ.

The exact ratio is

uxuy=−12χ(x).

Suppose χ has a simple zero and

c1s≤χ(x)2≤c2s.

Then

|χ|≍s1/2.

If

b≍s−α,

then

uy≍s1/4−α/2,ux≍s−1/4−α/2.

At k=3,

α=18,

and hence

b≍s−1/8,uy≍s3/16,ux≍s−5/16.

Consequently,

uy→0,|ux|→∞,uxuy=b→+∞.

The singular behavior is therefore anisotropic.

2.10. Ordinary spatial distance and the J3 scaling

A convenient periodic choice is

χ(x)=sin⁡(x−x∗).

Set

sgeo(x)=χ(x)2=sin2⁡(x−x∗).

Near a distinguished zero, let

r=x∗−x→0+.

Because χ has a simple zero,

χ(x)=−χ′(x∗)r+O(r2),

so

|χ|≍r,sgeo≍r2.

For the k=3 scaling

b≍sgeo−1/8,

one obtains

|b|≍r−1/4,
|uy|≍r3/8,

and

|ux|≍r−5/8.

The exponent identities

1−58=38

and

−58+38=−14

imply simultaneously

|χux|≍r3/8→0

and

|uxuy|≍r−1/4→∞.

In particular,

ux=o(1|χ|).

Thus the divergence of ux is weaker than the reciprocal vanishing scale 1/|χ|, and the identity

uy=−2χux

remains completely consistent with

uy→0.

For

b>0,

reality of

uy2=−2χb

requires

χ≤0.

Near a simple zero of sin⁡(x−x∗), one may therefore select a one-sided interval

J3=(x∗−δ,x∗),

on which

χ(x)<0.

Either real square-root branch

uy=±−2χb

may be selected, with the corresponding sign of

ux=buy.
Exact sine geometry/cancellation:valid for general x,Terminal power-law asymptotics:local near x=x∗+mπ.

2.11. Exact continuity and reconstruction of uz

The square-root law determines A=uy, but it does not prescribe C=uz.

Exact incompressibility requires

∂x(bA)+Ay+Cz=0.

Thus

Cz=−∂x(bA)−Ay.

Differentiating

A2=−2χb

gives

AyA=12byb

and

AxA=12χ′χ+12bxb.

Hence

∂x(bA)=b2A(bxb−χ′χ).

Therefore

Cz=−b2A(bxb−χ′χ)−A2byb.

In the distinguished terminal core,

b=b(s),s=t−z−θ∗,N,

one has

bx=by=0,Ay=0.

Since

∂z=−∂s,

continuity becomes

Cs=∂x(bA).

Using

A2=−2χb,

this simplifies to

Cs=Aχ′4χ2.

If

b=BNs−αFN(s)

and

IN(s)=∫0sr−α/2FN(r)1/2dr,

then

C(x,s)=C0(x)+χ′(x)−2BNχ(x)4χ(x)2IN(s).

For k=3 in the square-root tube,

C=c0+O(s3/16),c0>0.

Hence

uz=C→c0>0

even though the Lambert branch variable satisfies

QN→0.

This is another reason QN and uz must be kept conceptually distinct.

2.12. Distinguished-phase cancellation in the reconstruction

Every function depending on t,z only through

s=t−z−θ∗,N

satisfies

(∂t+∂z)f=0.

Consequently,

Ct+Cz=0

in the local distinguished-phase reconstruction.

The remainder

RN=2bXC+bCt+CzC

therefore reduces to

RN=2bXC.

The explicit square-root-core reconstruction yields

|XC|ΓNb=O(s)+O(s3/2),

and hence

|RN|ΓNb2=O(s)→0.

Thus exact incompressibility and the distinguished phase combine to produce precisely the small normalized error required by the terminal amplification argument.

2.13. General geometric reparametrization G(x)

The local parabolic geometry can be generalized by writing

X(x)=G(x),y=G(x)2.

Then

dydx=2G(x)G′(x).

The exactly matched horizontal velocity condition is

2G(x)G′(x)ux+uy=0.

For the full periodic sine choice

G(x)=sin⁡x,

one has

G′(x)=cos⁡x

and therefore

2G(x)G′(x)=2sin⁡xcos⁡x=sin⁡(2x).

Thus the exact periodic coefficient replacing the local approximation 2x is

sin⁡(2x).

Near

x=0,
sin⁡(2x)=2x+O(x3),

so the original local relation is recovered to leading order.

2.14. The derivative term E and its exact cancellation

Consider

E=uxuz,zx+uyuz,zy.

The curve

y=sin2⁡x

has tangent derivative

Dτ=∂x+sin⁡(2x)∂y.

For a general smooth scalar field F,

ddxF(x,sin2⁡x,z,t)=Fx+sin⁡(2x)Fy.

Thus the curve itself does not imply

Fx=sin⁡(2x)Fy.

The exact cancellation instead follows from a global characteristic structure.

Define

η=y+sin2⁡x.

Suppose

uz,z(x,y,z,t)=H(η,z,t)=H(y+sin2⁡x,z,t).

Then

uz,zx=Hηsin⁡(2x),

and

uz,zy=Hη.

Therefore

uz,zx=sin⁡(2x)uz,zy.

Define

D−=∂x−sin⁡(2x)∂y.

Then

D−uz,z=0.

If the horizontal velocity satisfies

sin⁡(2x)ux+uy=0,

then

(ux,uy)=ux(1,−sin⁡2x),

and hence

ux∂x+uy∂y=uxD−.

Consequently,

E=uxuz,zx+uyuz,zy=uxD−(uz,z)=0.

Thus

E=0

is an exact structural identity.

For a general G(x), take

η=y+G(x)2.

Then

ηx=2GG′,ηy=1.

If

uz,z=H(y+G(x)2,z,t)

and

2GG′ux+uy=0,

then

E=(2GG′ux+uy)Hη=0.

2.15. Geometry of Dτ and D−

The two operators are

Dτ=∂x+sin⁡(2x)∂y

and

D−=∂x−sin⁡(2x)∂y.

Their direction vectors are

τ=(1,sin⁡2x),v−=(1,−sin⁡2x).

Their dot product is

τ⋅v−=1−sin2⁡(2x)=cos2⁡(2x).

Thus the directions are not generally perpendicular. They are perpendicular only when

cos⁡(2x)=0,

i.e.

x=π4+kπ2,k∈Z.

Near the distinguished zero

x=0,

one has

sin⁡(2x)→0,

so

Dτ→∂x,D−→∂x.

They become parallel rather than perpendicular.

Their sum and difference are

Dτ+D−=2∂x

and

Dτ−D−=2sin⁡(2x)∂y.

Therefore

∂x=12(Dτ+D−)

and, whenever

sin⁡(2x)≠0,
∂y=Dτ−D−2sin⁡(2x).

The Dτ-characteristics satisfy

dydx=sin⁡(2x),

and therefore

y−sin2⁡x=C.

The D−-characteristics satisfy

dydx=−sin⁡(2x),

and therefore

y+sin2⁡x=C.

Thus the two families are

y=C+sin2⁡x

and

y=C−sin2⁡x.

They are reflected across the horizontal line y=C.

The horizontal velocity associated with

sin⁡(2x)ux+uy=0

is parallel to D−. Therefore

E=uxD−(uz,z),

so the condition

D−(uz,z)=0

has the simple geometric meaning that uz,z is constant in the horizontal velocity direction.

2.16. Global definition of the sine structure on T3

Let

η(x,y)=y+sin2⁡x.

Although the real representative of y depends on the chosen lift from T3 to R3, the composition

H(y+sin2⁡x,z,t)

defines a single-valued function on the torus provided H is periodic in its first spatial argument.

Assume

H(η+2π,z,t)=H(η,z,t)

and

H(η,z+2π,t)=H(η,z,t).

Define

h(x,y,z,t)=H(y+sin2⁡x,z,t).

Then

h(x+2π,y,z,t)=h(x,y,z,t),

because

sin2⁡(x+2π)=sin2⁡x.

Also,

h(x,y+2π,z,t)=h(x,y,z,t)

by periodicity of H in η, and

h(x,y,z+2π,t)=h(x,y,z,t)

by periodicity in z.

Hence

h(x,y,z,t)=H(y+sin2⁡x,z,t)

is defined globally on all of T3.

The structural prescription is

uz,z=H(y+sin2⁡x,z,t).

Differentiating gives everywhere on the torus

uz,zx=Hηsin⁡(2x),

and

uz,zy=Hη.

Therefore

uz,zx=sin⁡(2x)uz,zy

throughout

T3.

Equivalently,

(∂x−sin⁡(2x)∂y)uz,z=0on T3.

If

sin⁡(2x)ux+uy=0

globally, then

E=uxuz,zx+uyuz,zy=(sin⁡(2x)ux+uy)uz,zy=0.

Hence

E=0on T3.

2.17. Periodic primitive for uz

If

uz,z=H(y+sin2⁡x,z,t)

is to be the z-derivative of a periodic velocity component uz, rather than merely a periodic scalar field, one additional compatibility condition is required.

A derivative of a periodic function has zero mean over one period. Therefore it is necessary that

∫02πH(η,z,t)dz=0for every (η,t).

Under this condition, define

uz(x,y,z,t)=U0(x,y,t)+∫0zH(y+sin2⁡x,ζ,t)dζ,

where

U0(x,y,t)

is periodic in x and y.

Then

uz,z=H(y+sin2⁡x,z,t),

and

uz(x,y,z+2π,t)−uz(x,y,z,t)=∫zz+2πH(y+sin2⁡x,ζ,t)dζ=0.

Thus

uz(x,y,z+2π,t)=uz(x,y,z,t).

The sine-characteristic prescription therefore defines uz,z globally on T3, and, under the zero-z-mean condition, it admits a globally periodic primitive uz.

3. X=0 in [2](page 261 there)

Here we prove that X=0 in [2](page 261 there) With the derivative convention specified on the constrained geometry

y=x2,

the cancellation does occur directly, without introducing any H. Along this constraint,

dydx=2x.

Thus, when an x-derivative is converted to a y-derivative by the chain rule as intended,

∂∂x=dydx∂∂y=2x∂∂y.

Therefore, for uz,

uz,x=2xuz,y,

and, applying the same rule to uz,z,

uz,zx=2xuz,zy.

At the same time,

uy=−2xux.

Since x is independent of t and z,

uy,t=−2xux,t,

and

uy,zt=−2xux,zt.

Now substitute all four identities explicitly into X. Starting from

X=((δ−1)uxuz,t+2ρuzux,t)uz,zx+((δ−1)uyuz,t+2ρuzuy,t)uz,zy−uz,t(uzuz,xux,zt+uzuz,yuy,zt−uz,xuz,zux,t−uz,yuz,zuy,t)+uzuz,xux,tuz,tz+uzuz,yuy,tuz,tz.

Substitute

uy=−2xux,uy,t=−2xux,t,uy,zt=−2xux,zt,

and

uz,x=2xuz,y,uz,zx=2xuz,zy.

We obtain

(3.1)X=((δ−1)uxuz,t+2ρuzux,t)(2xuz,zy)+((δ−1)(−2xux)uz,t+2ρuz(−2xux,t))uz,zy−uz,t(uz(2xuz,y)ux,zt+uzuz,y(−2xux,zt)−(2xuz,y)uz,zux,t−uz,yuz,z(−2xux,t))+uz(2xuz,y)ux,tuz,tz+uzuz,y(−2xux,t)uz,tz.

Now expand every term. The first line becomes

2x(δ−1)uxuz,tuz,zy+4xρuzux,tuz,zy.

The second line becomes

−2x(δ−1)uxuz,tuz,zy−4xρuzux,tuz,zy.

Thus these cancel exactly:

2x(δ−1)uxuz,tuz,zy−2x(δ−1)uxuz,tuz,zy=0,

and

4xρuzux,tuz,zy−4xρuzux,tuz,zy=0.

Now consider the four terms inside the large parentheses:

uz(2xuz,y)ux,zt+uzuz,y(−2xux,zt)−(2xuz,y)uz,zux,t−uz,yuz,z(−2xux,t).

Expand:

=2xuzuz,yux,zt−2xuzuz,yux,zt−2xuz,yuz,zux,t+2xuz,yuz,zux,t.

Both pairs cancel:

2xuzuz,yux,zt−2xuzuz,yux,zt=0,

and

−2xuz,yuz,zux,t+2xuz,yuz,zux,t=0.

Therefore the entire uz,t-term is

0.

Finally, the last two terms are

2xuzuz,yux,tuz,tz

and

−2xuzuz,yux,tuz,tz.

Hence

2xuzuz,yux,tuz,tz−2xuzuz,yux,tuz,tz=0.

Putting everything together,

X=2x(δ−1)uxuz,tuz,zy−2x(δ−1)uxuz,tuz,zy⏟0+4xρuzux,tuz,zy−4xρuzux,tuz,zy⏟0−uz,t[2xuzuz,yux,zt−2xuzuz,yux,zt⏟0+−2xuz,yuz,zux,t+2xuz,yuz,zux,t⏟0]+2xuzuz,yux,tuz,tz−2xuzuz,yux,tuz,tz⏟0.

Consequently,

X=0.

So under the simultaneous constrained differentiation rules

y=x2,dydx=2x,uy=−2xux,

we have

uz,x=2xuz,y,uz,zx=2xuz,zy,

while

uy,t=−2xux,t,uy,zt=−2xux,zt,

and every term in X cancels pairwise, giving

X≡0.

This is the direct calculation referred to; no H(y+x2,z,t) representation is needed.

4. Generalized J3 Geometry and Exact Cancellation of X

We begin with the local definition of the distinguished geometric set J3 used in the original construction. In the local coordinates it is defined by

J3(x,y,z,t)={(x,t)∈(R3,R+):2xux+uy=0andy=x2}.

Thus, on J3,

uy=−2xux,y=x2.

The second relation gives

dydx=2x.

Consequently, under the constrained chain-rule differentiation used on J3,

∂∂x=dydx∂∂y=2x∂∂y.

In particular,

uz,x=2xuz,y,uz,zx=2xuz,zy.

This is the local form of the cancellation geometry. To periodicize the local geometry, replace the local coordinate x by a smooth periodic function

χ=χ(x).

The geometric curve is then taken to be

y=χ(x)2.

Differentiating gives

dydx=2χ(x)χ′(x).

Therefore the exact periodic analogue of the local relation

2xux+uy=0

which is compatible with the chain rule is

(4.1)2χ(x)χ′(x)ux+uy=0.

Equivalently,

(4.2)uy=−2χ(x)χ′(x)ux.

Because χ=χ(x) is independent of t and z, differentiation of (4.2) gives

(4.3)uy,t=−2χχ′ux,t,

and

(4.4)uy,zt=−2χχ′ux,zt.

Introduce

η=y+χ(x)2

and let

(4.5)uz(x,y,z,t)=H(η,z,t)=H(y+χ(x)2,z,t).

Then

ηx=2χχ′,ηy=1.

Hence

uz,x=Hηηx=2χχ′Hη,

whereas

uz,y=Hη.

Therefore

(4.6)uz,x=2χχ′uz,y,

or equivalently

(4.7)uz,x−2χχ′uz,y=0.

Similarly,

uz,z=Hz,

and therefore

uz,zx=2χχ′Hηz,uz,zy=Hηz.

Thus

(4.8)uz,zx=2χχ′uz,zy,

or

(4.9)uz,zx−2χχ′uz,zy=0.

4.1. Exact factorization of X

Consider

(4.10)X=((δ−1)uxuz,t+2ρuzux,t)uz,zx+((δ−1)uyuz,t+2ρuzuy,t)uz,zy−uz,t(uzuz,xux,zt+uzuz,yuy,zt−uz,xuz,zux,t−uz,yuz,zuy,t)+uzuz,xux,tuz,tz+uzuz,yuy,tuz,tz.

In [2] work was completed showing X vanishes on a space which can be extended to the whole Torus and is re-emphasized in this paper. The reader should confirm this to themselves by substituting the equations defining the space Ji into the X operator.

Using

uy=−2χχ′ux,
uy,t=−2χχ′ux,t,

and

uy,zt=−2χχ′ux,zt,

equation (4.10) becomes

(4.11)X=((δ−1)uxuz,t+2ρuzux,t)(uz,zx−2χχ′uz,zy)+[−uz,t(uzux,zt−uz,zux,t)+uzux,tuz,tz](uz,x−2χχ′uz,y).

But from the invariant representation (4.5),

uz,x−2χχ′uz,y=0

and

uz,zx−2χχ′uz,zy=0.

Consequently,

X=((δ−1)uxuz,t+2ρuzux,t)(0)+[−uz,t(uzux,zt−uz,zux,t)+uzux,tuz,tz](0),

and therefore

(4.12)X=0.

Thus the exact general cancellation condition is

(4.13)y=χ(x)2,2χ(x)χ′(x)ux+uy=0,uz=H(y+χ(x)2,z,t),

and these relations imply

(4.14)X≡0.

4.2. Specialization to χ(x)=sin⁡x

Now choose

(4.15)χ(x)=sin⁡x.

Then

χ′(x)=cos⁡x,

and hence

2χ(x)χ′(x)=2sin⁡xcos⁡x=sin⁡(2x).

The generalized J3 geometry therefore becomes

(4.16)y=sin2⁡x,sin⁡(2x)ux+uy=0.

Equivalently,

(4.17)uy=−sin⁡(2x)ux.

The invariant representation is

(4.18)uz=H(y+sin2⁡x,z,t).

Indeed,

uz,x=sin⁡(2x)Hη,uz,y=Hη,

and therefore

(4.19)uz,x−sin⁡(2x)uz,y=0.

Likewise,

uz,zx=sin⁡(2x)Hηz,uz,zy=Hηz,

so

(4.20)uz,zx−sin⁡(2x)uz,zy=0.

Substituting (4.19) and (4.20) into the factorized expression for X yields

(4.21)X≡0.

Therefore the exact periodic realization is

y=sin2⁡x,sin⁡(2x)ux+uy=0,uz=H(y+sin2⁡x,z,t),

and it gives the exact cancellation

X=0.

4.3. Recovery of the original local J3 geometry

The original local geometry is recovered near a simple zero of χ(x)=sin⁡x, for example near x=0. Since

sin⁡x=x+O(x3),

we have

sin2⁡x=x2+O(x4),

while

sin⁡(2x)=2x+O(x3).

Thus

y=sin2⁡x

reduces locally to

y=x2,

and

sin⁡(2x)ux+uy=0

reduces locally to

2xux+uy=0.

Hence

J3:y=x2,2xux+uy=0

is precisely the local square-root geometry recovered from the corrected periodic formulation

y=χ(x)2,2χ(x)χ′(x)ux+uy=0.

For the sine choice this reads

y=sin2⁡x,sin⁡(2x)ux+uy=0.

The factor χ′(x) is therefore essential in the exact global chain-rule formulation, while the original coefficient 2x is recovered automatically in the local limit near a simple zero of χ.

5. The space ΣC

For the local J3,

y=x2,2xux+uy=0.

The velocity relation gives

uy=−2xux.

Therefore

ux∂x+uy∂y=ux(∂x−2x∂y).

Define

D−loc=∂x−2x∂y.

Now extend this vector field into the surrounding (x,y)-space and calculate its characteristics:

dydx=−2x.

Integrating,

y=−x2+C,

so

y+x2=C.

Thus the characteristic leaves generated by the J3-derived velocity direction are

ΣCloc={(x,y,z):y+x2=C}.

The periodic version works the same way. Start with

y=χ(x)2

and the corrected exact matched relation

2χχ′ux+uy=0.

Then

uy=−2χχ′ux

and therefore

ux∂x+uy∂y=ux(∂x−2χχ′∂y)⏟D−.

Hence

D−=∂x−2χχ′∂y.

Its characteristics satisfy

dydx=−2χ(x)χ′(x).

But

ddxχ(x)2=2χ(x)χ′(x),

so

ddx(y+χ(x)2)=0.

Therefore

y+χ(x)2=C.

For

χ(x)=sin⁡x,

this becomes

ΣC={(x,y,z):y+sin2⁡x=C(mod2π)}.

So there is a precise derivation:

J3:y=χ2,2χχ′ux+uy=0⇓uy=−2χχ′ux⇓ux∂x+uy∂y=uxD−⇓D−=∂x−2χχ′∂y⇓dydx=−2χχ′⇓y+χ2=C⇓ΣC.

There is an interesting geometric feature here: the original J3 curves in one direction,

y=χ2,

whereas the D−-characteristics generated by its velocity relation curve in the opposite direction,

y=C−χ2.

Thus ΣC is not obtained by simply “shifting J3.” It is obtained by taking the velocity direction specified on J3, extending that direction into the ambient torus, and integrating its characteristic curves. So I would phrase the conclusion carefully:

J3 determines D− on J3; an ambient extension of D− then generates the family ΣC.

The extension step is an additional construction; ΣC does not follow from the set equation

y=χ2

alone. On each leaf

ΣC:y+sin2⁡x=C,

take

uz=H(C,z,t)=H(y+sin2⁡x,z,t).

Then, since

C=y+sin2⁡x,

we get

uz,x=sin⁡(2x)HC,uz,y=HC,

and hence

uz,x−sin⁡(2x)uz,y=0.

Differentiating in z,

uz,zx=sin⁡(2x)HCz,uz,zy=HCz,

so

uz,zx−sin⁡(2x)uz,zy=0.

At the same time, from the corrected J3-derived velocity relation extended to the ambient torus,

uy=−sin⁡(2x)ux,

and because sin⁡(2x) is independent of t and z,

uy,t=−sin⁡(2x)ux,t,
uy,zt=−sin⁡(2x)ux,zt.

Now recall the factorization of the algebraic quantity X:

X=((δ−1)uxuz,t+2ρuzux,t)(uz,zx−sin⁡(2x)uz,zy)+[−uz,t(uzux,zt−uz,zux,t)+uzux,tuz,tz](uz,x−sin⁡(2x)uz,y).

But on every ΣC,

uz,zx−sin⁡(2x)uz,zy=0

and

uz,x−sin⁡(2x)uz,y=0.

Therefore both factors multiplying the two large brackets vanish, and so

X=0on every ΣC.

Because

T3=⨆C∈S1ΣC,

every point of T3 lies on one of these leaves. Thus, provided the relations

uy=−sin⁡(2x)ux

and

uz=H(y+sin2⁡x,z,t)

are imposed throughout the ambient torus, the cancellation is not confined to one distinguished J3 surface. It becomes

X(x,y,z,t)≡0throughout T3.

The key mechanism is therefore

ΣC:y+sin2⁡x=C,uz=H(C,z,t),D−=∂x−sin⁡(2x)∂y,D−C=0,D−uz=0,D−uz,z=0,⟹X=0.

So: the family ΣC is precisely what allows the J3-derived cancellation mechanism to be propagated through the ambient space, once that extension is made consistently.

5.1. Preterminal periodic reconstruction and terminal release

The organization is

affine or tuned preterminal reconstruction⇓periodic continuation and compatibility⇓phase conversion to θ=t−z⇓release of the historical transverse tuning⇓A2=−2χb⇓uy=A,ux=b/A⇓uz=C from exact continuity⇓terminal characteristic comparison.

Thus the preterminal periodic machinery and the terminal square-root geometry are not imposed simultaneously as competing ansatzes. They belong to successive regimes of the construction.

5.2. Why the reconstruction does not create the amplification

The origin of the amplification is

ΓNb2.

The velocity reconstruction does not manufacture this source. It contributes the remainder

RN=2bXC+bCt+CzC.

The purpose of the reconstruction estimates is to prove

|RN|≪ΓNb2.

In the release layer,

|Rrel|ΓNb2=O(s7/4),

and in the square-root core,

|Rcore|ΓNb2=O(s).

Therefore

RN=o(ΓNb2).

The reconstruction embeds the scalar amplification mechanism into a divergence-free periodic velocity field while preserving

b=uxuy.

5.3. Complete introduction-level mechanism

The entire reduced architecture can now be summarized without ambiguity.

principal real Lambert branch W0⇓ΔN,k+1=1+W0(−e−1−ΔN,k)⇓QN≍s2−k⇓ΓN=∂ξlog⁡QN∼2−ks−1⇓b˙=ΓNb2+2CG−RN⇓RN=o(ΓNb2)⇓b˙≥(1−ϑ)ΓNb2⇓b>1⟹s˙=1−b<0⇓s↓0,QN↓0,b↑⇓b(s)≳s−(1−ϑ)2−k⇓T∗,N<∞,b(t)→+∞⇓b=uxuy⇓A2=−2χb⇓uy=A,ux=b/A⇓∇⋅u=0 reconstructs uz=C.

At k=3, the corresponding terminal hierarchy is

QN≍s1/8,ΓN≍s−1,b≍s−1/8,

together with

uy≍s3/16,ux≍s−5/16,uxuy≍s−1/8.

The exponent identity

−516+316=−18

verifies directly that

uxuy=b.

Likewise,

|χ|≍s1/2

gives

|χux|≍s1/2s−5/16=s3/16→0,

so

uy=−2χux→0

is consistent with the divergence of ux.

The physical interpretation is therefore not one of isotropic growth of all velocity components. The terminal reconstruction is anisotropic:

uy→0,|ux|→∞,uz→c0>0

in the local k=3 square-root scaling considered above, while

uxuy=b3→+∞.

Meanwhile,

QN→0

continues to control the singular logarithmic coefficient

ΓN=∂ξlog⁡QN,

without being identified with either uy or uz.

The resulting structure is therefore a self-strengthening terminal mechanism:

QN↓0⟹ΓN↑⟹b3↑⟹s↓0⟹QN↓0.

Within the hypotheses developed in the finite-depth analysis, this feedback reaches the terminal branch in finite characteristic time. The scalar amplification is then transferred to the velocity reconstruction through the exact identity

b3=uxuy.

This is the central mathematical architecture developed in the remainder of the paper.

5.4. The apparent division-by-zero problem

Write

b=b3=uxuy,A=uy,C=uz.

Then

ux=bA.

The horizontal operator appearing in the reduced b-equation is

X=A∂x+bA∂y.

The remainder is

(5.1)RN=2bXC+bCt+CzC.

Thus the concern is genuine: if

A=uy⟶0

on the zero set I3, division by A cannot simply be declared harmless.

6. The zero set I3

I3 should be defined as the zero set of the geometric profile χ used in the square-root reconstruction.

For

χ:T→R,

define

I3:=Z(χ)={x∈T:χ(x)=0}.

Although A=uy→0 and

ux=bA

is pointwise singular at the terminal set, the parabolic wedge geometry is sufficient to retain local L2 integrability of the velocity. Thus the nonmoving coterminal construction does not fail at the level of kinetic energy.

The distinction appears at the derivative level. For the exact nonmoving k=3 terminal profile, the same wedge extending to ρ=0 produces a nonintegrable spatial contribution from ∂zux, so that the corresponding local H1 estimate required by the Leray dissipation integral is not obtained.

The moving-coterminal construction is introduced precisely to resolve this derivative-level obstruction. It preserves the local L2 velocity control already present in the nonmoving wedge while replacing the fixed spatial approach to I3 by a time-dependent terminal geometry for which the spatial gradient norm is finite at every preterminal time and its growth is integrable in time. The required Leray condition is therefore

∫0T∗‖∇u(t)‖L2(T3)2dt<∞.

Accordingly, the role of the moving coterminal construction is not to restore finite kinetic energy—that property is already compatible with the nonmoving I3 wedge—but to obtain the stronger spacetime derivative estimate required by the Leray energy inequality. Thus the distance variable used in the L2 analysis is

ρ(x):=distT⁡(x,I3).

Near a simple zero x∗∈I3,

χ(x∗)=0,χ′(x∗)≠0,

and Taylor expansion gives

χ(x)=χ′(x∗)(x−x∗)+O((x−x∗)2).

Consequently,

|χ(x)|≍ρ(x)

locally near I3. This is precisely what permits us to replace powers of |χ| by powers of the distance ρ in the shrinking-wedge L2 calculations.

For the periodic choice used elsewhere in this paper,

χ(x)=sin⁡x,

the zero set on T=R/(2πZ) is

I3={[0],[π]}.

Equivalently, on the representative interval [0,2π),

I3={0,π}.

Both zeros are simple because

χ′(0)=1,χ′(π)=−1.

In the full three-dimensional torus, the corresponding geometric zero set is the union of the two periodic yz-tori

I3=I3×T2={(x,y,z)∈T3:sin⁡x=0}.

This also clarifies an important notational point: I3 is not the Lambert branch zero s=0. It is the spatial zero set of χ(x). The terminal parabolic geometry couples the two through

χ(x)2≍s,

or equivalently near I3,

ρ(x)2≍s.

One must establish either sufficient cancellation in the numerator or sufficient integrability of the quotient.

The construction does not give an exact matching zero of b. In fact, in the terminal core b grows while A vanishes. Hence the relevant mechanism is not pointwise cancellation of b/A.

It must instead be an integrability mechanism.

6.1. Exact square-root reconstruction

In the terminal region the reconstruction imposes

(6.1)A2=−2χb.

On the selected branch put

q=−χ>0.

Then

A2=2qb,

and therefore

(6.2)A=2qb,ux=bA=b2q.

This immediately shows why ux need not vanish with A.

For k=3,

b(s)∼Bs−1/8.

Consequently, with q and s still treated as independent variables,

(6.3)A∼2Bq1/2s−1/16,

and

(6.4)ux=bA∼B2q−1/2s−1/16.

This is the correct expression to use in a functional estimate.

It is important not to substitute q∼s1/2 before taking partial derivatives.

6.2. What happens along the square-root trajectory?

If one evaluates on

q2≍s,

then

q≍s1/2,

and (6.3)–(6.4) become

A≍s1/4s−1/16=s3/16,

while

ux≍s−1/4s−1/16=s−5/16.

Thus

(6.5)uy=A≍s3/16→0,ux=bA≍s−5/16→∞.

Therefore the proposed construction does not remove the quotient singularity pointwise.

The question is instead whether

ux∈Lloc2

and, more importantly,

RN∈Lloc2.

6.3. Ordinary spatial distance

Let

ρ=dist⁡(x,I3)

locally.

If χ has a simple zero,

χ(x∗)=0,χ′(x∗)≠0,

then

|χ(x)|≍ρ.

Hence

q≍ρ.

In the square-root core

s≍ρ2.

Therefore

b≍ρ−1/4,
A≍ρ3/8,

and

(6.6)ux=bA≍ρ−5/8.

So there is indeed a singular quotient at I3.

But

|ux|2≍ρ−5/4

is not the complete integrability test, because the singular terminal region is not a product neighborhood of fixed s-width.

6.4. The terminal region is a shrinking wedge

Introduce

(6.7)Ωϵ={0<ρ<ϵ,c1ρ2≤s≤c2ρ2,y∈T}.

The condition y∈T means that no shrinking or localization is being imposed in the y-direction. The entire periodic y-circle is included. Since

s=t−z−θ∗,N,

at fixed t,

|ds|=|dz|.

Also a simple zero of χ gives

dx≍dρ.

Thus

(6.8)dV≍dρdsdy.

6.5. Direct L2 estimate for the dangerous quotient

Before imposing the tube relation,

ux≍ρ−1/2s−1/16.

Therefore

|ux|2≲ρ−1s−1/8.

Consequently

∫Ωϵ|ux|2dV≲∫0ϵρ−1∫c1ρ2c2ρ2s−1/8dsdρ.

Now

∫c1ρ2c2ρ2s−1/8ds=87(c27/8−c17/8)ρ7/4.

Hence

∫Ωϵ|ux|2dV≲C∫0ϵρ3/4dρ.

Thus

(6.9)∫Ωϵ|ux|2dV≲Cϵ7/4<∞.

Therefore

(6.10)bA=ux∈Lloc2(Ωϵ).

This directly answers the first part of the objection: division by A creates a pointwise singularity, but not an L2 singularity on the proposed shrinking core.

6.6. Relation to weighted Sobolev/Hardy spaces

There is also a natural abstract formulation.

Near I3,

A≍ρ3/8,

so

1A≍ρ−3/8.

Thus

bA

belongs naturally to a distance-weighted space.

A useful notation is

Lβ2(Ω)={f:∫Ωρ2β|f|2dV<∞}.

The classical Hardy principle controls division by distance for functions having appropriate vanishing or trace properties; weighted Hardy–Sobolev inequalities extend this to distance weights and more general singular sets.

But here we should not claim that Hardy’s inequality by itself solves the problem. In particular, b does not vanish on I3.

Instead, the explicit reconstruction supplies

A2=2ρb

up to uniformly comparable factors, and therefore

bA=b2ρ.

The shrinking-wedge calculation (6.9) then directly establishes the needed integrability. Weighted Hardy theory provides an appropriate surrounding functional framework, but the explicit estimate is the substantive proof.

6.7. Continuity must also survive the quotient

Exact incompressibility is

∂x(bA)+Ay+Cz=0.

In the local terminal construction

Ay=0.

Since

s=t−z−θ∗,∂z=−∂s,

we obtain

(6.11)Cs=∂x(bA).

Thus the singular quotient is not hidden: its x-derivative explicitly determines C.

With

b=Bs−1/8

and

A=2Bqs−1/16,

we have

bA=B2qs−1/16.

Since

qx=−χ′,
∂x(bA)=B22χ′q−3/2s−1/16.

Hence

Cs=P(x)s−1/16,

where

(6.12)P(x)=2B4χ′(x)q−3/2.

Integrating in s,

C=C0(x)+P(x)I(s),

with

I(s)=∫0sτ−1/16dτ=1615s15/16.

Thus

(6.13)C=C0(x)+42B15χ′(x)q−3/2s15/16.

This is the exact local continuity reconstruction at leading order.

6.8. The reconstructed C remains finite in the core

On

q≍s1/2,

we have

q−3/2≍s−3/4.

Therefore the correction in (6.13) behaves as

s−3/4s15/16=s3/16.

Hence

(6.14)C=c0+O(s3/16).

In particular, if

c0>0,

then C remains bounded and nonzero sufficiently close to the terminal set.

So exact continuity does not force C to inherit the divergence of b/A.

6.9. The remainder simplifies in the distinguished phase

Recall

RN=2bXC+bCt+CzC.

Every t,z-dependent reconstructed quantity depends through

s=t−z−θ∗,N.

Therefore

(∂t+∂z)f(s)=0.

In particular,

(6.15)Ct+Cz=0.

Hence

(6.16)RN=2bXC.

Moreover, in the local reconstruction

Cy=0.

Thus

XC=ACx+bACy=ACx.

This is particularly important for:

(6.17)bACy=0

exactly in the local core.

Thus the most visibly dangerous occurrence of 1/A in XC vanishes structurally rather than merely being estimated.

6.10. Computation of Cx

From

C=C0+PI,
P=2B4χ′q−3/2,

and q(x)=−χ(x)>0, we obtain

Px=2B4[χ″q−3/2+32(χ′)2q−5/2].

Hence

Cx=C0′+PxI.

The leading singular term is

Cx∼2B432(χ′)2q−5/21615s15/16.

Therefore

(6.18)Cx∼22B5(χ′)2q−5/2s15/16.

6.11. Direct estimate of RN

Use

b∼Bs−1/8,
A∼CAq1/2s−1/16,

and

Cx∼CCq−5/2s15/16.

Then

RN=2bACx

has scaling

RN≍s−1/8q1/2s−1/16q−5/2s15/16.

The spatial powers give

q−2,

and the s-powers give

−18−116+1516=34.

Therefore

(6.19)RN=O(q−2s3/4).

Since q≍ρ,

(6.20)RN=O(ρ−2s3/4).

6.12. Why a trajectory calculation looks critical

On the single curve

s=cρ2,

equation (6.20) becomes

RN≍ρ−2(ρ2)3/4=ρ−1/2.

Thus

|RN|2≍ρ−1.

If one integrates only along ρ,

∫0ϵρ−1dρ=∞.

This is exactly why the trajectory calculation appeared logarithmically critical.

But a trajectory is not the three-dimensional spatial measure.

6.13. Full shrinking-wedge L2 calculation

Using (6.20),

|RN|2≲ρ−4s3/2.

Hence

∫Ωϵ|RN|2dV≲∫0ϵ∫c1ρ2c2ρ2ρ−4s3/2dsdρ.

But

∫c1ρ2c2ρ2s3/2ds=25(c25/2−c15/2)ρ5.

Therefore

∫Ωϵ|RN|2dV≲C∫0ϵρdρ.

Consequently

(6.21)∫Ωϵ|RN|2dV≤Cϵ2.

Equivalently,

(6.22)‖RN‖L2(Ωϵ)≤Cϵ.

Here this is substantially stronger than saying that

RN/(ΓNbN2)→0

along a trajectory.

6.14. Why a singular C0 is neither necessary nor desirable

One possibility considered earlier was to choose C0(x) singular enough to cancel the leading term in Cx.

Along the tube,

ACx≍ρ−1/4.

Since

A≍ρ3/8,

cancellation would require

C0′≍−ρ−5/8.

Then

C0−C∗≍ρ3/8,

so C0 itself remains continuous.

But

|C0′|2≍ρ−5/4,

and hence

C0∉Hloc1

in the ordinary transverse variable.

Therefore this attempted cancellation would trade one problem for another.

The wedge estimate (6.21) shows that no such cancellation is needed.

6.15. Higher Lambert terms do not provide the missing cancellation either

Write more generally

bN(s)=BNs−1/8FN(s),

with

FN(s)=1+aNsμ+O(sμ+δ).

Then

FN1/2=1+aN2sμ+⋯,

and

IN(s)=∫0sτ−1/16FN(τ)1/2dτ.

Direct integration gives

(6.23)IN(s)=1615s15/16[1+15aN2(15+16μ)sμ+O(sμ+δ∗)].

These terms improve only the higher-order correction.

They do not remove the leading

ρ−2s3/4

behavior of RN.

Thus the L2 closure comes from the shrinking geometry, not from a hidden Lambert cancellation.

6.16. The release layer must also be controlled

The terminal core cannot simply be glued directly to the regular preterminal reconstruction.

This paper introduces

(6.24)Aζ2=(1−ζ)A∗2+ζ(−2χb).

If

b≍s−α,

matching a bounded nonzero A∗ requires

|χ|≍sα.

For k=3,

α=18.

Thus the release occurs at

(6.25)|χ|≍s1/8.

The final core instead satisfies

|χ|≍s1/2.

Since

s1/2≪s1/8,

there is substantial scale separation between release and terminal core.

6.17. Uniform release cutoff

Set again

q=−χ>0

and define

η=qsα.

Choose a fixed smooth cutoff

ζ=Z(η)

with

Z=1

on the inner side and

Z=0

on the outer side.

Then

ζx=−Z′(η)χ′sα,

so

(6.26)ζx=O(s−α).

Similarly,

(6.27)ζxx=O(s−2α).

The constants are independent of N provided the cutoff and geometric bounds for χ are fixed uniformly.

6.18. Release nondegeneracy

Write

PN=(1−ζ)A∗,N2+2ζqbN,AN=PN1/2.

Assume

0<a−≤A∗,N≤a+

and

0<B−≤sαbN(s)≤B+.

On the active release,

q≍sα,

so

qbN≍1.

Consequently

(6.28)0<A−≤AN≤A+<∞

with constants independent of N.

The quotient b/A is not encountering A=0 inside the active release layer.

The vanishing of A occurs deeper inside the terminal core, where the shrinking-wedge estimate applies.

6.19. Exact cutoff commutator

Differentiate (6.24):

2AAx=ζx(2qb−A∗2)+2ζqxb.

Therefore

(6.29)Ax=ζx(2qb−A∗2)+2ζqxb2A.

The first term is the cutoff commutator:

CA=ζx(2qb−A∗2)2A.

Because

ζx=O(s−α),qb=O(1),A≍1,

we obtain

CA=O(s−α).

The geometric term has exactly the same order:

ζqxbA=O(s−α).

Thus

(6.30)Ax=O(s−α).

At k=3,

Ax=O(s−1/8).

Likewise,

(6.31)Axx=O(s−1/4).

6.20. Exact continuity in the release layer

Instead of independently cutting off C, reconstruct it from

Cs=∂x(bA).

Since

bx=0,
(6.32)Cs=−bAxA2.

In the release,

b=O(s−α),Ax=O(s−α),A−1=O(1),

so

(6.33)Cs=O(s−2α).

For

α=18,
Cs=O(s−1/4).

Integration yields

(6.34)C−Crel=O(s3/4).

An important consequence is that there is no independent C-cutoff commutator. The interpolation enters C through exact incompressibility.

6.21. Release remainder

Because

Ct+Cz=0,

the release remainder is again

Rrel=2bXC.

The release calculation gives

XC=O(s1−3α),

so

Rrel=O(s1−4α).

For

α=18,

this gives

(6.35)Rrel=O(s1/2).

Therefore

|Rrel|2=O(s),

which is plainly locally integrable.

Thus the release cutoff does not destroy the L2 estimate.

6.22. Matching release to core

Choose the cutoff flat at its inner endpoint:

ζ=1,∂xjζ=0.

Then

Arel2=−2χb=Acore2.

Choosing the same positive branch gives

(6.36)Arel=Acore.

Moreover,

Axrel=−χ′bA.

Differentiating the core identity

Acore2=−2χb

gives exactly the same result:

Axcore=−χ′bAcore.

Hence

(6.37)Axrel=Axcore.

Because both sides reconstruct C from

Cs=∂x(b/A),

choosing the same integration datum gives

(6.38)Crel=Ccore,(Cs)rel=(Cs)core.

And b itself is never interpolated, so

(6.39)brel=bcore.

6.23. Why the core cutoff does not create another commutator

The release takes place at

q≍s1/8.

The terminal core is at

q≍s1/2.

Since

s1/2≪s1/8,

the release cutoff has already reached

ζ≡1

before the terminal core is entered.

Thus throughout the intermediate region and the final core,

A2=−2χb

holds exactly.

There is consequently no second shrinking cutoff at

q≍s1/2.

That removes a potentially dangerous source of derivative commutators.

6.24. The resulting functional framework

The appropriate local framework can therefore be summarized as follows.

Define

ρ(x)=|χ(x)|

near the simple zero set I3, and use the anisotropic/parabolic terminal region

Ωϵ={0<ρ<ϵ, c1ρ2<s<c2ρ2}.

The reconstructed velocity satisfies

A=uy≍ρ1/2s−1/16,
ux=bA≍ρ−1/2s−1/16,

and

C=uz=c0+O(ρ−3/2s15/16).

The natural anisotropic weighted control is therefore of the form

ρ1/2s1/16ux=O(1),

together with the corresponding weighted bounds for Cx and RN.

Rather than requiring

ux∈L∞,

which is false, one proves

ux∈L2(Ωϵ)

and

RN∈L2(Ωϵ).

This is compatible with the general philosophy of distance-weighted Hardy–Sobolev spaces, where distance to a singular or boundary set supplies the weight.

But the present result is not merely an appeal to that theory: the relevant integrals have been computed directly.

6.25. Uniform-N version

Proposition 5.1 (Uniform geometric normalization of the terminal core). To turn the local estimate into a uniform theorem, impose the following genuinely quantitative hypotheses:

(6.40)0<B−≤BNFN(s)≤B+,
(6.41)|∂sjFN(s)|≤Cjs−j,

for the required derivative orders,

(6.42)0<a−≤A∗,N≤a+,

and simple-zero geometry

(6.43)0<κ−≤|χ′(x∗)|≤κ+,

together with uniform bounds for the higher derivatives of χ.

Finally, require the parabolic core constants

(6.44)0<c1<c2<∞

to be independent of N.

Under these assumptions all constants in (6.9), (6.21), and the release estimates can be chosen independently of N.

6.26. Precise release-to-core theorem

Theorem 5.1 (Uniform local L2 control across the release and square-root core). Let

bN(s)=BNs−1/8FN(s)

satisfy uniform positive upper and lower bounds and the required uniform derivative estimates. Let χ have a uniformly nondegenerate simple zero I3, and let

ρ=|χ|.

Assume that the terminal reconstruction is

AN2=−2χbN

inside the square-root region and that the release is performed through

AN,ζ2=(1−ζ)A∗,N2+ζ(−2χbN)

with a fixed smooth cutoff active at

|χ|≍s1/8,

flat at its endpoints.

Reconstruct CN from exact incompressibility,

CN,s=∂x(bNAN),

and assume the distinguished phase

s=t−z−θ∗,N,

so that

CN,t+CN,z=0.

Suppose finally that the terminal core is the parabolic region

c1ρ2≤s≤c2ρ2

with c1,c2 independent of N.

Then

bNAN=ux,N∈Lloc2,

despite

AN→0

on I3.

The release interpolation satisfies

AN,x=O(s−1/8),CN,s=O(s−1/4),

uniformly in N.

The release remainder satisfies

Rrel,N=O(s1/2)

and is uniformly locally L2.

In the square-root core,

Rcore,N=O(ρ−2s3/4),

and

‖Rcore,N‖L2(Ωϵ)≤Cϵ

with C independent of N.

Consequently,

(6.45)supN‖RN‖L2(Ωrel∪Ωcore,ϵ)≤Crel+Ccoreϵ<∞.

6.27. What the integrated functional calculation resolves

The shrinking-wedge calculation addresses the division-by-A objection in the full spatial geometry rather than only along a trajectory. Although A=uy vanishes and b/A=ux diverges pointwise, the explicit square-root reconstruction gives

bA∈Lloc2

on the parabolic core. Moreover, in the local core

Cy=0,bACy=0,

and the remaining compact core remainder satisfies

‖Rcore,N‖L2(Ωϵ)≤Cϵ.

No exact matching zero of ux and no singular choice of C0 is required.

The uniform-N version is separated from the fixed-N functional statement and is tied to the uniform coefficient bounds and parabolic-core normalization stated later. The viscous derivative calculation is retained as a completed part of the terminal reconstruction. In particular,

Cxx=O(s−13/16),Cyy=0,

so

Cxx+Cyy=O(s−13/16)=o(s−2),

and therefore

b(Cxx+Cyy)=O(s−15/16)=o(s−17/8).

Together with

ν(Cbss+bCss)∼9νBc064s−17/8>0,

this yields

νCΔb+νbΔC=9νBc064s−17/8(1+o(1))>0.

Thus the functional L2 framework, release-to-core matching, and terminal viscous derivative estimates are mutually consistent parts of the reconstruction analysis.

6.28. Why the reconstruction does not create the amplification

The origin of the amplification is

ΓNb2.

The velocity reconstruction does not manufacture this source. It contributes the remainder

RN=2bXhC+bCt+CzC.

The purpose of the reconstruction estimates is to prove

|RN|≪ΓNb2.

In the release layer,

|Rrel|ΓNb2=O(s7/4),

and in the square-root core,

|Rcore|ΓNb2=O(s).

Therefore, in the normalized terminal asymptotic sense,

RN=o(ΓNb2).

This statement must be supplemented near I3 by the functional estimate of Section 5.24, because A=uy→0 and b/A=ux is pointwise singular there. Theorem 5.1 proves the local L2 control and the release-to-core matching on the explicitly defined parabolic region. Uniformity in N is then reduced to the separate geometric normalization criterion isolated in Proposition 5.1. The identity

b=uxuy

is preserved. The later viscous terminal sections supply the corresponding second-derivative and multidimensional estimates used by the reconstruction.

6.29. Complete introduction-level mechanism

The entire reduced architecture can now be summarized without ambiguity.

principal real Lambert branch W0⇓ΔN,k+1=1+W0(−e−1−ΔN,k)⇓QN≍s2−k⇓ΓN=∂ξlog⁡QN∼2−ks−1⇓b˙=ΓNb2+2CG−RN⇓RN=o(ΓNb2)⇓b˙≥(1−ϑ)ΓNb2⇓b>1⟹s˙=1−b<0⇓s↓0,QN↓0,b↑⇓b(s)≳s−(1−ϑ)2−k⇓T∗,N<∞,b(t)→+∞⇓b=uxuy⇓A2=−2χb⇓uy=A,ux=b/A⇓∇⋅u=0 reconstructs uz=C.

At k=3, the corresponding terminal hierarchy is

QN≍s1/8,ΓN≍s−1,b≍s−1/8,

together with

uy≍s3/16,ux≍s−5/16,uxuy≍s−1/8.

The exponent identity

−516+316=−18

verifies directly that

uxuy=b.

Likewise,

|χ|≍s1/2

gives

|χux|≍s1/2s−5/16=s3/16→0,

so

uy=−2χux→0

is consistent with the divergence of ux.

The physical interpretation is therefore not one of isotropic growth of all velocity components. The terminal reconstruction is anisotropic:

uy→0,|ux|→∞,uz→c0>0

in the local k=3 square-root scaling considered above, while

uxuy=b3→+∞.

Meanwhile,

QN→0

continues to control the singular logarithmic coefficient

ΓN=∂ξlog⁡QN,

without being identified with either uy or uz.

The resulting structure is therefore a self-strengthening terminal mechanism:

QN↓0⟹ΓN↑⟹b3↑⟹s↓0⟹QN↓0.

Within the hypotheses developed in the finite-depth analysis, this feedback reaches the terminal branch in finite characteristic time. The scalar amplification is then transferred to the velocity reconstruction through the exact identity

b3=uxuy.

This is the central mathematical architecture developed in the remainder of the paper.

7. Leray Energy Verification and Moving Coterminal Reconstruction

7.1. Editorial roadmap and logical status

The preceding J3 geometry supplies the local square-root reconstruction used below. The purpose of the present inserted development is to test that reconstruction against the Leray–Hopf energy class and then to record, in chronological mathematical order, the redesign forced by that test.

There are two logically distinct conclusions. First, the literal nonmoving coterminal wedge extending to ρ=0 has finite local kinetic energy but fails the required local H1 spatial estimate: the most singular fixed-coordinate derivative produces a divergent spatial integral. Second, a moving preterminal/coterminal front can change the lower radial endpoint and makes the terminal Lt2Hx1 power count integrable, provided the moving collar is constructed with the matching, derivative, incompressibility, and periodic-compatibility properties audited below.

Accordingly, the moving-front calculation resolves the local singular-exponent obstruction to the Leray dissipation integral at the level of the explicit collar power count. This does not yet close the global construction: the final periodic divergence correction must be localized so that it does not destroy the exact outer Lambert comparison field. Earlier unsuccessful or conditional constructions are retained because they explain why the final one-front, exact-overlap, explicit-collar architecture is required.

The controlling terminal time scale used in the successful power count is

sf(t)≍(T∗−t)8/9,r(t)≍sf(t)1/8≍(T∗−t)1/9.

The final explicit collar is simultaneously restricted by s≍r2, so its local volume is O(r3) with this paper’s exact measure dρdsdy. This distinction is essential and supersedes any earlier intermediate estimate that informally identified the local phase s in the collar with the temporal front parameter sf.

Reading convention for the audit.

Each construction below is treated as a response to the obstruction established immediately before it. When a proposed construction does not close that obstruction, the text records the precise residual defect before passing to the next redesign. Thus an intermediate ansatz is not to be read as an established solution merely because a later section replaces it. The chain of implications terminates in the explicit moving-collar estimate and then in the still-separate periodic compatibility obligation stated in the consolidated conclusion.

8. Leray Energy as a Global Consistency Test

The Leray energy relation now becomes a global consistency test for the favorable-viscosity b3 mechanism. The key is that the favorable sign in the scalar b3-equation and the dissipative sign in the velocity energy equation are not contradictory: they concern different quantities.

For the unforced incompressible Navier–Stokes system on T3,

ut+(u⋅∇)u+∇p=νΔu,∇⋅u=0,

a smooth solution satisfies

(8.1)12‖u(t)‖L2(T3)2+ν∫0t‖∇u(τ)‖L2(T3)2dτ=12‖u0‖L2(T3)2.

For a Leray–Hopf weak solution this becomes the corresponding energy inequality.

The reconstruction is

b=b3=uxuy,A=uy,ux=bA,C=uz,

and in the k=3 terminal core we have,

b≍s−1/8,uy≍s3/16,ux≍s−5/16,

with C=uz→c0>0.

So the first Leray question is:

Does ‖u(t)‖2 remain finite as t↑T∗?

The shrinking geometry is crucial. Near the distinguished zero,

s≍ρ2.

Therefore

ux2≍s−5/8≍ρ−5/4.

But the terminal wedge has s-width of order ρ2. Consequently its contribution to kinetic energy behaves schematically like

(8.2)∫0ϵρ−5/4ρ2dρ=∫0ϵρ3/4dρ<∞.

Thus

ux∈Lloc2

despite ux→∞ pointwise. This is precisely the mechanism we emphasize: the singular quotient is locally L2 on the full shrinking wedge.

So finite kinetic energy is compatible with this pointwise amplification.

The more stringent Leray test is the dissipation:

(8.3)∫0T∗‖∇u(t)‖22dt<∞.

This is where we now need to concentrate.

For example, if

ux≍s−5/16,

then the most singular distinguished derivative is formally

∂zux≍s−21/16,

and hence

(8.4)|∂zux|2≍s−21/8.

If we simply substitute s≍ρ2, this becomes

|∂zux|2≍ρ−21/4.

The spatial wedge factor ds≍ρ2 would give

(8.5)∫0ϵρ−21/4ρ2dρ=∫0ϵρ−13/4dρ,

which diverges.

But (8.5) is not yet the Leray dissipation integral, because Leray requires spacetime integration:

∫0T∗∫T3|∇u|2dxdt.

The trapping law supplies the missing temporal geometry. We have:

−s˙=HN+γNb

and in the terminal regime

b≳Bs−λ.

Hence

(8.6)dt≲sλds.

We use exactly this relation to prove that the characteristic reaches s=0 in finite time.

Therefore the next calculation must be a genuine four-dimensional terminal integral, not merely a spatial path estimate:

(8.7)∫terminal spacetime tube(|∇ux|2+|∇uy|2+|∇uz|2)dxdydzdt.

There is also an important conceptual point concerning our viscosity discussion. In the b-equation we found that the older transformed viscosity can have the favorable asymptotic sign

(8.8)V[b,C]∼9νB64s−17/8>0.

That does not mean viscosity injects kinetic energy. When the original vector NS equation is multiplied by u and integrated over the torus,

(8.9)ν∫T3u⋅Δudx=−ν∫T3|∇u|2dx.

So viscosity remains globally dissipative. The positive sign in (8.8) occurs only after the nonlinear product transformation b=uxuy. Diffusion of ux,uy can increase their product locally while still decreasing total kinetic energy.

That resolves the apparent sign paradox.

The lower-barrier theorem says essentially

b≥Bs−λ⟹b→∞

in finite characteristic time, with viscosity reinforcing the comparison at first contact.

Leray now imposes the independent requirement

(8.10)supt<T∗‖u(t)‖22+2ν∫0T∗‖∇u(t)‖22dt≤‖u0‖22

for an unforced Leray–Hopf solution.

Therefore we now have three separate statements to prove:

b3→∞

from the lower-barrier mechanism;

supt<T∗‖u(t)‖2<∞

from the shrinking terminal geometry; and most importantly,

(8.11)∫0T∗‖∇u(t)‖22dt<∞.

The first two are compatible with what we have so far. The third is the decisive Leray audit.

And this is stronger than simply checking whether the pointwise viscosity has a favorable sign. The standard energy theory gives uniform Lt∞Lx2 and Lt2Hx1 control for Leray–Hopf solutions; such energy control by itself does not imply smoothness.

So the calculation I would do next is explicit: insert

ux=bA,uy=A,uz=C,

with

b∼Bs−1/8,A2=−2χb,χ2≍s,

compute all nine derivatives

ui,j,i,j=x,y,z,

using the exact fixed-coordinate differentiation rule, and then evaluate

ν∫T∗−δT∗∫Ω(t)|∇u|2dVdt.

That calculation will tell us whether the proposed terminal blowup is actually compatible with the Leray energy inequality, rather than merely having finite instantaneous L2 energy.

9. Fixed Coterminal Wedge: Bounded Kinetic Energy but Borderline Dissipation

Transition.

The global energy identity by itself does not decide whether the reconstructed terminal field belongs to the Leray class. We therefore pass first to the literal fixed coterminal wedge and separate the kinetic-energy question from the spatial-derivative question. The conclusion of this section is intentionally diagnostic: finite local L2 energy survives, whereas the derivative estimate remains unresolved until the exact terminal measure is inserted in the next section.

Carrying the calculation through reveals a very sharp issue: with the current k=3 terminal reconstruction and the stated parabolic tube geometry, the instantaneous kinetic energy can remain finite, but the Leray dissipation appears borderline divergent in time. This is the point we need to resolve.

For reference, the Leray–Hopf requirement is

u∈L∞(0,T;L2)∩L2(0,T;H1),

with

12‖u(t)‖22+ν∫0t‖∇u(τ)‖22dτ≤12‖u0‖22

in the unforced case.

The k=3 terminal reconstruction is

b=b3≍s−1/8,uy≍s3/16,ux≍s−5/16,

together with the parabolic geometry

χ2≍s.

this paper also explicitly uses

−s˙=HN+γNb

and obtains finite-time arrival at s=0.

The important point is that we must differentiate the fixed-coordinate formula, not simply differentiate s−5/16 after imposing χ2≍s.

Write

q=−χ>0.

From

A2=2qb,b=Bs−1/8,

we have

uy=A=2Bqs−1/16,

and

(9.1)ux=bA=B2q−1/2s−1/16.

Only after differentiation do we impose

(9.2)q2≍λs.

Because s=t−z−θ∗,

∂z=−∂s.

Therefore

(9.3)ux,z=116B2q−1/2s−17/16.

Consequently,

(9.4)|ux,z|2=B512q−1s−17/8.

Using

q≍s1/2,

gives

(9.5)|ux,z|2≍s−21/8.

Now introduce the transverse distance

ρ≍|q|.

Since

s≍ρ2,

we obtain

(9.6)|ux,z|2≍ρ−21/4.

Suppose the terminal wedge is

c1ρ2≤s≤c2ρ2.

For each ρ, its s-width is therefore

Δs≍ρ2.

Hence the spatial contribution from this single derivative alone is

∫|ux,z|2dV≳C∫ρmin(t)ϵρ−21/4ρ2dρ,

assuming the remaining periodic direction has a nonvanishing width. Thus

∫|ux,z|2dV≳C∫ρmin(t)ϵρ−13/4dρ.

Integration gives

(9.7)‖∇u(t)‖22≳Cρmin(t)−9/4.

This is why the finite-N cutoff matters enormously.

Now use the trapping dynamics. In the terminal regime,

−s˙=HN+γNb∼γNBs−1/8.

Thus

dt∼−1γNBs1/8ds.

Integrating from s to zero,

T∗−t≍s9/8.

Therefore

(9.8)smin(t)≍(T∗−t)8/9.

Since

ρmin2≍smin,

we obtain

(9.9)ρmin(t)≍(T∗−t)4/9.

Substitute this into (9.7):

‖∇u(t)‖22≳C[(T∗−t)4/9]−9/4.

The exponents cancel exactly:

4994=1.

Hence

(9.10)‖∇u(t)‖22≳CT∗−t.

And therefore

∫T∗−δT∗‖∇u(t)‖22dt≳C∫T∗−δT∗dtT∗−t.

The latter is logarithmically divergent:

(9.11)∫T∗−δT∗‖∇u(t)‖22dt=+∞

under these geometric assumptions.

This is a substantially stronger conclusion than our previous pointwise viscosity comparison.

Notice also that there is no contradiction with the finite kinetic energy calculation. For ux,

|ux|2≍q−1s−1/8.

On q≍ρ, s≍ρ2,

|ux|2≍ρ−5/4.

Consequently,

∫wedge|ux|2∼∫0ϵρ−5/4ρ2dρ=∫0ϵρ3/4dρ<∞.

Thus we get the interesting pair

supt<T∗‖u(t)‖2<∞,

but apparently

(9.12)∫0T∗‖∇u(t)‖22dt=∞.

That would put the proposed terminal field outside the Leray–Hopf energy class if the estimate truly applies to a spacetime region of the assumed size. Leray–Hopf solutions require finite integrated H1 dissipation.

9.1. There is one important escape route we must audit

Equation (9.11) used a terminal tube having a nonvanishing measure in the remaining y-direction. If the full multidimensional reconstruction shrinks in another transverse direction as s→0, there is an extra measure factor.

For example, suppose that transverse width behaves like

Ly(s)≍sβ.

Since

smin≍(T∗−t)8/9,

the previous borderline estimate acquires

Ly≍(T∗−t)8β/9.

Then schematically,

‖∇u(t)‖22≲(T∗−t)−1+8β/9,

which is time-integrable for every

β>0.

So the borderline logarithm means that the exact three-dimensional geometry is decisive. Even a genuine shrinking width in the second transverse direction can change the Leray conclusion.

There is another possibility: the fixed-N reconstruction may modify the derivative estimate before the formal ρ=0 core is reached. That also has to be included rather than extrapolating the limiting asymptotic all the way to zero.

Therefore I would not yet write in this paper that the construction violates the Leray inequality. What we have proved conditionally is more precise:

(9.13)If the k=3 terminal law persists down to ρmin(t)≍(T∗−t)4/9,and the second transverse width remains bounded below,then ‖∇u(t)‖22≳(T∗−t)−1,so the Leray dissipation diverges logarithmically.

This pinpoints the next calculation very cleanly: we need to extract the exact three-dimensional terminal tube/reconstruction and determine its y-width (or corresponding second transverse Jacobian). That will decide whether the logarithmic divergence in (9.11) survives or is canceled by the shrinking geometry.

10. Exact Terminal Volume Element and the Nonmoving Coterminal Obstruction

Transition.

The preceding wedge calculation identifies the derivative as the possible obstruction but does not yet use the exact (ρ,s,y) volume geometry. We now recompute the dissipation with that measure. This removes the ambiguity: the nonmoving terminal wedge retains local L2 velocity control but fails the required local H1 estimate for the literal terminal profile. The next section therefore asks whether the finite-N release geometry already supplies a positive lower cutoff that could remove this divergence.

I extracted the terminal region actually used in this paper. It answers the y-width question quite clearly:

ΩN,ε={0<ρ<ε,c1ρ2≤s≤c2ρ2,y∈T}.

At fixed t,

s=t−z−θ∗,N,|ds|=|dz|,

and near the simple zero of χ,

dx≍dρ.

Therefore this paper explicitly uses

dV≍dρdsdy.

Most importantly, y∈T is the whole periodic y-direction. There is no shrinking y-width in the stated terminal core.

So the possible rescue I mentioned—an additional factor Ly(s)≍sβ, β>0—is not yet present in the current paper geometry.

10.1. Recompute the Leray dissipation using the exact volume element

Before imposing s≍ρ2,

ux=b2q≍ρ−1/2s−1/16,

because q=|χ|≍ρ and b≍s−1/8. This is exactly the fixed-coordinate scaling used by this manuscript.

Since

∂z=−∂s,

we obtain

ux,z≍ρ−1/2s−17/16,

and therefore

(10.1)|ux,z|2≍ρ−1s−17/8.

Now integrate over the actual core:

(10.2)∫ΩN,ε|ux,z|2dV≍∫Tdy∫0ερ−1∫c1ρ2c2ρ2s−17/8dsdρ.

The s-integral is

∫c1ρ2c2ρ2s−17/8ds=89(c1−9/8−c2−9/8)ρ−9/4.

Thus

(10.3)∫ΩN,ε|ux,z|2dV≍C∫0ερ−13/4dρ.

But

∫0ερ−13/4dρ=∞.

Therefore something stronger than our previous time-integrated calculation follows:

(10.4)ux,z∉Lloc2(ΩN,ε)

for the literal terminal power law extended all the way to ρ=0.

So under the stated terminal geometry,

u∉Hloc1

in that terminal core once the singular profile has actually reached the zero set.

This should be contrasted with the valid L2 calculation for ux. It has

|ux|2≲ρ−1s−1/8,

and hence

∫0ερ−1∫c1ρ2c2ρ2s−1/8dsdρ≲ε7/4.

That is precisely why we obtain

ux∈Lloc2,

despite the pointwise divergence.

But differentiating changes the integrability completely:

ux∈Lloc2⇏∇ux∈Lloc2.

10.2. What this means for Leray

There is an important temporal distinction.

For every preterminal time t<T∗, if the physical solution has not yet populated the arbitrarily small-ρ portion of this asymptotic core, there can still be an effective cutoff

ρ≥ρmin(t)>0.

Then

‖∂zux(t)‖22≍Cρmin(t)−9/4.

Using the k=3 trapping relation

b≍s−1/8,−s˙≍s−1/8,

gives

T∗−t≍smin9/8,

and hence

ρmin(t)≍smin1/2≍(T∗−t)4/9.

Therefore

(10.5)‖∂zux(t)‖22≍(T∗−t)−1.

Consequently,

(10.6)∫T∗‖∂zux(t)‖22dt∼∫T∗dtT∗−t=∞.

So the full y-width does not cancel the logarithmic divergence. It confirms it.

10.3. But this exposes an even deeper issue in the geometry

There is something we need to be careful about before declaring the construction incompatible with Leray.

We defined the core at fixed time as

c1ρ2≤s≤c2ρ2,y∈T.

Thus, mathematically, for every arbitrarily small ρ>0, the region contains corresponding arbitrarily small

s∼ρ2.

There is no ρmin(t) in that definition.

Therefore, if the asymptotic field

b(s)≍s−1/8

is literally imposed throughout that entire core at any fixed preterminal time, then (10.3) says

‖∇u(t)‖2=∞

already at that time—not merely in the limit t↑T∗.

That would be much more serious than logarithmic divergence at T∗.

Accordingly, this paper now needs to specify exactly how the moving terminal front intersects

ΩN,ε(t).

There must be a distinction between the geometric test region

0<ρ<ε,c1ρ2≤s≤c2ρ2

and the portion of that region actually occupied by the singular asymptotic profile at a given t<T∗.

Without such a moving cutoff, the current k=3 reconstruction proves local L2 velocity but does not give the H1 regularity required for a Leray–Hopf solution.

There is one other feature worth noting. It explicitly says that in the local core

Cy=0,

and its three-dimensional integration nevertheless retains the entire y∈T direction. Thus there is no hidden y-Jacobian or y-localization in the stated reconstruction that supplies the missing factor.

So we have narrowed the problem substantially:

current terminal core: y∈T⇓no second transverse shrinking factor⇓|∂zux|2≍ρ−1s−17/8⇓∫Ωε|∂zux|2dV≍∫0ερ−13/4dρ=∞.

The next mathematical question is therefore no longer the y-width. this paper itself answers that: it is full width. The decisive question is whether its finite-N, preterminal/release construction supplies a genuine moving lower cutoff ρmin,N(t)>0 and, if so, whether its approach to zero is slow enough that

∫0T∗,Nρmin,N(t)−9/4dt<∞.

That is the quantity we should extract and calculate next.

11. Finite-N Release Geometry and the Missing Positive Moving Cutoff

Transition.

The nonmoving obstruction can only be avoided if the actual preterminal geometry prevents the singular wedge from reaching ρ=0 at every t<T∗. We therefore inspect the finite-N release theorem rather than adding a cutoff by assumption. The audit below shows that the existing release interpolation does not itself provide the needed positive moving lower endpoint. This negative conclusion motivates, rather than already proves, the moving-preterminal construction of the following section.

I checked the finite-N release/core construction specifically for the missing moving cutoff. The result is important:

We did not previously have a positive moving lower cutoff ρmin,N(t)>0.

The release layer does something different.

The paper defines the release by

AN,ζ2=(1−ζ)A∗,N2+ζ(−2χbN),

with

η=qs1/8,ζ=Z(η).

Thus the active release occurs at

q≍s1/8,

whereas the terminal parabolic core occurs at

q≍s1/2.

this paper explicitly observes that

s1/2≪s1/8

and therefore the release cutoff has already become identically 1 before the parabolic core is reached. In fact it emphasizes that there is no second cutoff at q≍s1/2.

That means the release interpolation cannot provide the ρmin(t) we were looking for.

12. What the finite-N theorem actually assumes

At every fixed finite N, Theorem 106.1 takes

bN(s)=BNs−1/8FN(s)

and defines

(12.1)ΩN,ε={0<ρN<ε,c1,NρN2≤s≤c2,NρN2,y∈T}.

There is no lower endpoint

ρN≥ρmin,N(t).

The interval is literally

0<ρN<ε.

Moreover, the uniform-N argument only seeks constants c1,c2 independent of N:

c1χN2≤s≤c2χN2.

It does not introduce a finite-N exclusion of χN=0.

So the cutoff cannot presently be obtained from either:

finite N,release interpolation,oruniform-N normalization.

13. Now redo the H1 calculation exactly

this paper gives, before imposing the tube relation,

(13.1)ux,N=bNAN≍ρ−1/2s−1/16.

At fixed x,

∂z=−∂s,

so

ux,z≍ρ−1/2s−17/16.

Therefore

(13.2)|ux,z|2≍ρ−1s−17/8.

Using this paper’s exact volume element

dV≍dρdsdy,

we get

∫ΩN,ε|ux,z|2dV≍|T|∫0ερ−1∫c1ρ2c2ρ2s−17/8dsdρ=C∫0ερ−13/4dρ.

Hence

(13.3)∫ΩN,ε|ux,z|2dV=∞.

This occurs already at the spatial level.

So the time integration cannot repair it:

∫0T∗∫ΩN,ε|∇u|2dVdt

contains a nonnegative contribution that is already infinite whenever the literal terminal profile occupies all of (12.1).

14. Compare with what we originally proved

There is no contradiction in its L2 velocity calculation.

For velocity itself,

|ux|2≍ρ−1s−1/8,

and therefore

∫Ωε|ux|2dV≍∫0ερ−1∫c1ρ2c2ρ2s−1/8dsdρ≍∫0ερ3/4dρ<∞.

That is exactly the local L2 theorem actually proved originally.

The issue is specifically one derivative higher:

(14.1)ux∈Lloc2butux,z∉Lloc2

for the stated terminal profile.

14.1. A useful general criterion

We can make the obstruction more transparent.

Suppose generally

b≍s−α,

while retaining

A2=2ρb.

Then

ux=bA≍ρ−1/2s−α/2.

Thus

ux,z≍ρ−1/2s−1−α/2,

and

|ux,z|2≍ρ−1s−2−α.

Integrating over

s≍ρ2

gives

∫c1ρ2c2ρ2s−2−αds≍ρ−2−2α.

Hence

(14.2)∫|∂zux|2dV≍∫0ερ−3−2αdρ.

For every amplification exponent

α>0,

this diverges. In fact even α=0 gives ρ−3.

So merely changing 1/8 to another positive finite-depth Lambert exponent will not solve this particular H1 issue while retaining the same square-root reconstruction and the same parabolic core.

14.2. Does the favorable viscosity help?

This is where two facts must remain separate.

We proved that its scalar viscous contribution has favorable sign. In its terminal calculation it obtains

νCΔb+νbΔC=9νBc064s−17/8(1+o(1))>0.

That can reinforce the b3 lower barrier.

But it does not imply

∇u∈L2.

These are different statements:

favorable sign in transformed b3 equation

versus

ν∫|∇u|2<∞ in the original velocity energy relation.

The first does not supply the second.

14.2.1. What would actually be needed

Our earlier proposed moving cutoff would need to alter the region to something like

(14.3)Ωε(t)={ρmin(t)<ρ<ε,c1ρ2<s<c2ρ2,y∈T}.

Then

‖∂zux(t)‖22≍Cρmin(t)−9/4

for α=1/8.

If

ρmin(t)≍(T∗−t)β,

then

∫T∗‖∂zux(t)‖22dt∼∫T∗(T∗−t)−9β/4dt.

This is finite exactly when

(14.4)β<49.

The characteristic estimate we previously obtained,

ρmin(t)≍(T∗−t)4/9,

is precisely the borderline value and produces

∫T∗dtT∗−t,

a logarithmic divergence.

14.2.2. Current mathematical status

We originally supported the following chain:

ΩN,ε={0<ρ<ε, c1ρ2≤s≤c2ρ2, y∈T}⇓ux∈Lloc2but∂zux∉Lloc2for the literal s−1/8 core.

The original release layer cannot fix this because it explicitly terminates before the q≍s1/2 core and introduces no second core cutoff.

The promising route is to revisit what the Leray requirement actually demands before versus at the singular time and determine whether a classical solution on every 0<T<T∗ can have finite dissipation while the limiting terminal profile itself is outside H1. For that, we need to distinguish carefully between the terminal asymptotic profile and the actual t<T∗ velocity field. That distinction is now the decisive point.

15. Moving Preterminal Core and Leray-Compatible Dissipation Audit

**Transition.**Because the finite-N release layer does not supply the required lower endpoint, we now introduce a genuinely moving preterminal scale and determine exactly what rate would make the Leray dissipation time-integrable. At this stage the calculation is a criterion and power-counting audit, not yet a complete reconstruction theorem. Its output is the admissible motion law and the list of geometric properties that the subsequent construction must realize.

The preceding fixed terminal wedge

ΩN,ε={0<ρN<ε,c1,NρN2≤s≤c2,NρN2,y∈T}

is sufficient for the local L2-integrability of the velocity component ux,N=bN/AN, but it is not sufficient by itself for the H1-integrability required by the Leray energy estimate.

Indeed, at cascade depth k=3,

bN(s)≍s−1/8,

and the square-root reconstruction

AN2=−2χNbN

gives, before imposing the parabolic relation,

ux,N=bNAN≍ρN−1/2s−1/16,ρN=|χN|.

Since

s=t−z−θ∗,N,∂z=−∂s,

one obtains

∂zux,N≍ρN−1/2s−17/16,

and hence

|∂zux,N|2≍ρN−1s−17/8.

On the fixed wedge,

dV≍dρNdsdy,

and therefore

∫ΩN,ε|∂zux,N|2dV≍∫0ερN−1∫c1,NρN2c2,NρN2s−17/8dsdρN≍∫0ερN−13/4dρN=+∞.

Thus the fixed terminal wedge proves

ux,N∈Lloc2,

but it does not prove

∇uN∈Lloc2.

This distinction is essential for comparison with the Leray energy inequality.

15.1. Moving preterminal truncation

For every preterminal time t<T∗,N, introduce a dynamically determined inner radius

ρmin,N(t)>0

and define the moving terminal core by

ΩN,ε(t)={ρmin,N(t)<ρN<ε,c1,NρN2<s<c2,NρN2,y∈T}.

The quantity ρmin,N(t) is not introduced as an arbitrary regularizing cutoff. It must be determined from the finite-N preterminal geometry and must satisfy

ρmin,N(t)>0(t<T∗,N),

while

ρmin,N(t)⟶0ast↑T∗,N.

Thus the singular terminal set is approached only in the limit t↑T∗,N.

For each fixed t<T∗,N, the moving core remains separated from ρN=0, and therefore the reconstructed velocity remains H1 on the truncated terminal region.

15.2. Dissipation generated by the moving core

Using

|∂zux,N|2≍ρN−1s−17/8,

we obtain

∫ΩN,ε(t)|∂zux,N|2dV≍∫ρmin,N(t)ερN−1∫c1,NρN2c2,NρN2s−17/8dsdρN≍CN∫ρmin,N(t)ερN−13/4dρN.

Since

∫ρ−13/4dρ=−49ρ−9/4,

we obtain

‖∂zux,N(t)‖L2(ΩN,ε(t))2≍CNρmin,N(t)−9/4

to leading order as

ρmin,N(t)↓0.

Consequently a necessary terminal integrability test for the Leray dissipation is

∫T∗,Nρmin,N(t)−9/4dt<∞.

This condition must be derived from the actual finite-N preterminal dynamics rather than imposed independently.

15.3. Power-law criterion

Suppose that the moving inner radius satisfies

ρmin,N(t)≍(T∗,N−t)βN

near T∗,N.

Then

‖∂zux,N(t)‖22≍(T∗,N−t)−9βN/4.

Therefore

∫T∗,N‖∂zux,N(t)‖22dt<∞

precisely at the level of this contribution when

9βN4<1,

or equivalently

βN<49.

The value

βN=49

is critical. In that case

‖∂zux,N(t)‖22≍1T∗,N−t,

and consequently

∫T∗,N‖∂zux,N(t)‖22dt

diverges logarithmically.

Thus the Leray-compatible regime requires the dynamically generated inner scale to approach zero more slowly than

(T∗,N−t)4/9.

15.4. Relation with the characteristic trapping law

The moving spatial scale must now be compared with the distinguished characteristic dynamics.

For the k=3 terminal law,

bN(s)≍BNs−1/8,

and the characteristic equation has the form

−s˙=HN+γNbN.

If the bN-term dominates asymptotically, then

−s˙≍γNBNs−1/8.

Hence

dt≍−s1/8ds,

and integration gives

T∗,N−t≍s9/8.

Equivalently,

s(t)≍(T∗,N−t)8/9.

If one were to identify the moving inner radius directly with the parabolic characteristic scale

ρmin,N(t)2≍s(t),

then

ρmin,N(t)≍(T∗,N−t)4/9.

This is exactly the critical exponent:

βN=49.

Consequently,

‖∂zux,N(t)‖22≍(T∗,N−t)−1,

and the corresponding dissipation is logarithmically divergent.

Therefore the simple identification

ρmin,N2≍s

does not by itself close the Leray dissipation estimate.

15.5. Required geometric refinement

The preceding calculation isolates the exact additional condition required of the moving preterminal geometry.

If

ρmin,N(t)≍s(t)γN,

then, since

s(t)≍(T∗,N−t)8/9,

we have

ρmin,N(t)≍(T∗,N−t)8γN/9.

The dissipation condition

βN<49

therefore becomes

8γN9<49,

or

γN<12.

Hence a sufficient geometric condition at the level of the ∂zux contribution is

ρmin,N(t)≳s(t)γ,0<γ<12.

Because

sγ>s1/2(s↓0, γ<1/2),

this means that the preterminal field must remain farther from the zero set than the limiting parabolic scale until the singular time.

The limiting parabolic geometry

ρ2≍s

may then emerge only at

t=T∗,N,

rather than being occupied all the way down to ρ=0 at every preterminal time.

15.6. Compatibility with the favorable viscous sign

This moving-core refinement does not alter the favorable-sign calculation for the transformed b3-viscosity.

The lower-barrier comparison remains a pointwise statement in the active terminal region. At a hypothetical first downward contact, the viscous operator retains the maximum-principle ordering derived previously.

Likewise, in the terminal core the leading longitudinal contribution has the favorable sign

νCΔb+νbΔC∼9νBc064s−17/8>0.

The role of the moving inner scale is different. It controls the global-in-space dissipation integral

ν∫T3|∇u|2dx

for t<T∗,N.

Thus there is no contradiction between

positive transformed viscosity in the b3-comparison

and

finite dissipative Leray integral for the velocity field.

They impose distinct conditions on the construction.

15.7. Revised status of the Leray closure

The moving-core formulation converts the previous H1 obstruction into a precise geometric proof obligation.

The construction must establish, from the finite-N equations and not by assumption, a function

ρmin,N(t)>0,t<T∗,N,

such that

ρmin,N(t)→0ast↑T∗,N,

and

∫T∗,Nρmin,N(t)−9/4dt<∞.

A sufficient power-law condition is

ρmin,N(t)≳(T∗,N−t)β,0<β<49.

Equivalently, relative to the characteristic branch variable, a sufficient condition is

ρmin,N(t)≳s(t)γ,0<γ<12.

The critical parabolic scaling

γ=12

produces logarithmic divergence and therefore does not close the Leray dissipation estimate.

Accordingly, the moving cutoff is not a technical regularization. Its existence and rate are a new geometric theorem that must be derived from the finite-N preterminal reconstruction.

Until that theorem is proved, this paper establishes the L2-integrability of the terminal velocity on the shrinking parabolic wedge, but the full Leray

L2(0,T∗,N;H1(T3))

closure remains conditional on the moving-core estimate above.

16. Status of the Moving-Core Approach and the Theorem Still Required

**Transition.**The moving-core power count removes the fixed-wedge exponent obstruction only conditionally. Before choosing a concrete front, we therefore isolate what has actually been established and what remains to be constructed. This status section prevents the preceding criterion from being mistaken for a completed solution and provides the hypotheses that the release-scale front must satisfy.

This is the mathematically safer rewrite because it does not claim that introducing the moving region has solved the Leray problem. It identifies the new theorem that would have to be proved.

There is also a useful connection to the existing release geometry. We currently have the release at

q≍s1/8,

while the final parabolic core is

q≍s1/2.

Since 1/8<1/2, the release scale already has exactly the qualitative property

s1/8≫s1/2.

So rather than inventing an unrelated ρmin, a potentially natural revision would be to investigate whether the finite-N release/front itself can define

ρmin,N(t)≍s(t)1/8

until sufficiently close to T∗,N. If that were genuinely derived from the reconstruction, then

β=8918=19<49,

and the contribution above would behave as

‖∂zux(t)‖22≍(T∗−t)−1/4,

which is time-integrable:

∫T∗(T∗−t)−1/4dt<∞.

That is promising, but the current paper says the release is already completed before entering the q≍s1/2 core. Therefore we cannot simply identify the existing release boundary with ρmin; we would have to redesign the release-to-core transition and then recheck continuity, the remainder RN, viscosity, and all nine velocity derivatives.

That redesigned transition is the next calculation to do before rewriting the rest of the paper around the moving cutoff.

17. Release Scale as a Moving Coterminal Front

**Transition.**The previous status audit leaves one constructive question: can a scale already present in the Lambert/release geometry serve as the required moving lower endpoint? We test the release scale for that role. The favorable dissipation exponent obtained here solves the leading local power-counting problem, but it immediately creates an interface-matching problem; that new problem is carried explicitly into the next section rather than being regarded as already resolved.

Proceeding with the redesign, the natural candidate is to keep the release scale active as the inner preterminal boundary rather than allowing the reconstructed square-root core to extend immediately to ρ=0.

We already have defined the release variable

η=qs1/8,q=|χ|,

so the release occurs at q≍s1/8. It also states that the existing cutoff becomes identically 1 before the parabolic core q≍s1/2 is reached.

17.1. Define the moving front

Take

(17.1)ρmin,N(t)=κNsf(t)1/8,κN>0,

where sf(t) is the distinguished phase distance of the moving front.

The reconstructed singular core is then used only on

(17.2)ΩN,ε(t)={ρmin,N(t)<ρN<ε,c1,NρN2<s<c2,NρN2,y∈T}.

Inside

0≤ρN≤ρmin,N(t)

we retain a regular preterminal reconstruction instead of imposing the limiting square-root field all the way to ρ=0.

We define the parabolic core with 0<ρN<ε.

See figure 1 below which is used throughout the analyses:

17.2. Leray dissipation becomes integrable at this scale

We already established

|∂zux|2≍ρ−1s−17/8.

Therefore

(17.3)‖∂zux(t)‖22≍Cρmin(t)−9/4.

With

ρmin≍sf1/8,

this becomes

(17.4)‖∂zux(t)‖22≲Csf−9/32.

The characteristic law

−s˙f≍sf−1/8

gives

Regions with different layers used in this paper.
Figure 1. Regions with different layers used in this paper.
dt≍sf1/8dsf.

Consequently

∫T∗‖∂zux(t)‖22dt≲C∫0s0sf−9/32sf1/8dsf=C∫0s0sf−5/32dsf.

Since

−532>−1,

we obtain

(17.5)∫T∗‖∂zux(t)‖22dt<∞.

Equivalently,

sf(t)≍(T∗−t)8/9

implies

ρmin(t)≍(T∗−t)1/9,

and hence

(17.6)‖∂zux(t)‖22≲(T∗−t)−1/4.

That contribution is Leray-integrable.

But this only checks one derivative. We need the whole gradient.

17.3. Audit uy=A

In the square-root region,

(17.7)A=2ρb≍ρ1/2s−1/16.

The z-derivative satisfies

Az≍ρ1/2s−17/16,

so

|Az|2≍ρs−17/8.

Integration over the parabolic slice gives

‖Az(t)‖22≲∫ρminερ∫c1ρ2c2ρ2s−17/8dsdρ≍∫ρminερ−5/4dρ≍ρmin−1/4.

Thus

(17.8)‖uy,z(t)‖22≲ρmin−1/4.

For ρmin≍sf1/8,

‖uy,z‖22≲sf−1/32.

Hence

(17.9)∫‖uy,z‖22dt≲∫0sf−1/32+1/8dsf=∫0sf3/32dsf<∞.

So uy,z is substantially less restrictive than ux,z.

17.4. The transverse derivative ux,x

This one must be treated carefully because differentiating ρ−1/2 produces another singular factor.

Locally, because χ has a simple zero,

|ρx|≍1.

From

ux≍ρ−1/2s−1/16,

we obtain

ux,x≍ρ−3/2s−1/16.

Therefore

(17.10)|ux,x|2≍ρ−3s−1/8.

Integrating over s≍ρ2,

‖ux,x(t)‖22≍∫ρminερ−3∫c1ρ2c2ρ2s−1/8dsdρ≍∫ρminερ−3ρ7/4dρ=∫ρminερ−5/4dρ.

Hence

(17.11)‖ux,x(t)‖22≍ρmin−1/4.

Again, with ρmin≍sf1/8,

‖ux,x(t)‖22≲sf−1/32,

and this is spacetime integrable.

So far the worst derivative remains

ux,z.

17.5. uy,x

From

uy=A≍ρ1/2s−1/16,

we get

Ax≍ρ−1/2s−1/16.

Thus

|Ax|2≍ρ−1s−1/8.

This is precisely the same integrand that appeared in the L2 estimate for ux:

∫ρ−1∫c1ρ2c2ρ2s−1/8dsdρ≍∫ρ3/4dρ.

Therefore

(17.12)uy,x∈L2

even without relying on the moving lower cutoff.

The y-derivatives vanish in the local y-independent core used by this manuscript.

17.6. What about uz=C?

this paper reconstructs

CN(x,s)=C0,N(x)+42BN15χN′(x)qN−3/2s15/16,

and on q2≍s,

CN=c0,N+O(s3/16).

It also takes the local datum independent of y, so

Cy=0.

Its derivative estimate is

(17.13)Cx=C0′(x)+O(ρ−5/2s15/16).

Squaring the singular part gives

|Cx|2≲ρ−5s15/8.

Integrating through the parabolic width,

∫|Cx|2dV≲∫ρminερ−5∫c1ρ2c2ρ2s15/8dsdρ≍∫ρminερ−5ρ23/4dρ=∫ρminερ3/4dρ.

Thus

(17.14)Cx∈L2

without difficulty.

For Cz=−Cs,

Cs≍ρ−3/2s−1/16,

and therefore

|Cz|2≍ρ−3s−1/8.

Its spatial integral is

∫ρminερ−5/4dρ≍ρmin−1/4,

which is again spacetime-integrable for the proposed sf1/8 front.

We therefore obtain, for the outer moving square-root region, the schematic derivative hierarchy

(17.15)derivativeleading dependence on ρminux,zρmin−9/4ux,xρmin−1/4uy,zρmin−1/4uy,xO(1)uz,zρmin−1/4uz,xO(1)

with the local y-derivatives zero.

Thus the single decisive term really is

ux,z.

And the release-scale front

ρmin≍sf1/8

is more than sufficient for its spacetime integrability.

17.7. But there is a new interface issue

We cannot simply truncate the field at

ρ=ρmin(t).

A sharp cutoff would produce distributional derivatives and destroy precisely the H1 estimate we are trying to establish.

We therefore need a smooth moving interpolation. A natural choice is

ηN(x,t)=ρN(x)ρmin,N(t)

and a fixed smooth function

Z(η)={0,η≤1,1,η≥2,

with smooth transition for 1<η<2.

Define

(17.16)AN2=(1−Z(ηN))Areg,N2+Z(ηN)(−2χNbN).

Then

AN=Areg,N

in the regular inner region and

AN2=−2χNbN

in the outer square-root region.

The transition layer has thickness

Δρ≍ρmin.

Its derivatives satisfy

(17.17)Zx=Z′(η)ρxρmin=O(ρmin−1),

and

Zt=−Z′(η)ρρ˙minρmin2.

On the transition support ρ≍ρmin,

(17.18)Zt=O(ρ˙minρmin).

For

ρmin≍sf1/8,
ρ˙minρmin=18s˙fsf.

Since

|s˙f|≍sf−1/8,

we get

(17.19)|Zt|=O(sf−9/8).

This is highly singular pointwise.

That means the moving front solves the spatial H1 problem but introduces a potentially serious time-dependent commutator into the reconstructed equations.

This is now the next obstruction to audit.

17.8. Why the interface may nevertheless be manageable

The transition volume also shrinks.

Because

ρ≍ρmin≍sf1/8,

the transition thickness in x is

O(sf1/8).

If its relevant s-width is still parabolic,

Δs≍ρmin2≍sf1/4,

then the spatial transition volume scales like

(17.20)|Ωtr(t)|≍sf1/8sf1/4=sf3/8

up to the fixed y-length.

So every new commutator has to be tested as

amplitude2×sf3/8×dt,dt≍sf1/8dsf.

For example, Zt alone would give

|Zt|2|Ωtr|dt≍sf−9/4sf3/8sf1/8dsf=sf−7/4dsf,

which diverges.

But Zt never appears alone: it multiplies the difference

Acore2−Areg2.

Therefore the crucial matching condition is the order to which the two reconstructions agree at the moving interface.

If

Acore2−Areg2=O(sfμ),

then the squared moving-cutoff contribution scales like

sf−7/4+2μdsf.

It is integrable precisely if

−74+2μ>−1,

i.e.

(17.21)μ>38.

This gives us a concrete design criterion.

17.8.1. Revised construction target

We should therefore require the inner regular reconstruction to match the outer square-root reconstruction at least to order

(17.22)Acore2−Areg2=O(sf3/8+δ),δ>0.

Preferably we should impose flat matching of sufficiently high order, just as used in the flat endpoint for its original release. The present release construction explicitly uses a cutoff flat at the inner endpoint to match A, Ax, C, and Cs.

So there is an existing mechanism we can adapt rather than invent from scratch.

The resulting revised architecture is therefore

regular inner field⟶flat moving transition⟶square-root outer core

with

ρmin,N(t)≍sf(t)1/8.

This gives the desired Leray scaling for the outer core, but we have not yet closed the construction: the next calculation has to prove that the moving transition’s Zt, Zx, and Zxx commutators are integrable and that exact incompressibility can still reconstruct C without introducing a worse term.

That is the decisive calculation before replacing the corresponding release/core sections of the paper.

18. Moving-Interface Matching and Cutoff-Commutator Audit

**Transition.**The release-scale front gives the desired time-integrable lower cutoff, but a moving cutoff introduces new spatial and temporal commutator terms. We now quantify those terms. Where the naive interpolation fails, the failure is used to determine the matching order required of the amplitude. The resulting high-order requirement is then implemented, rather than merely assumed, in the next section.

Continuing from the moving-front construction, the key point is that we can sharpen the previous estimate considerably. The dangerous quantity is not Zt by itself; it is Zt multiplied by the mismatch between the inner regular reconstruction and the outer square-root reconstruction.

The existing work already uses exactly this philosophy in its stationary release: it interpolates A2, chooses the cutoff flat at an endpoint, and reconstructs C from incompressibility rather than imposing a separate C-cutoff. We can adapt that architecture to a moving interface.

18.1. Moving interpolation

Put

rN(t):=ρmin,N(t)≍sf(t)1/8,ηN=ρNrN(t).

Choose Z∈C∞ with

Z(η)=0(η≤1),Z(η)=1(η≥2),

and define

(18.1)AN2=Areg,N2+Z(ηN)DN,

where

(18.2)DN:=Acore,N2−Areg,N2,Acore,N2=−2χNbN.

The transition region is

TN(t)={rN(t)<ρN<2rN(t)}.

We require Areg,N>0, so that the regular inner reconstruction does not divide by zero.

The current paper’s release already interpolates A2 rather than A itself. That remains the preferable choice.

18.2. Exact derivatives of the moving cutoff

Since

η=ρr(t),

we have

ηx=ρxr,

and

ηt=−ρr˙r2=−ηr˙r.

Hence

(18.3)Zx=Z′(η)ρxr=O(r−1),

while

(18.4)Zt=−Z′(η)ηr˙r=O(|r˙|r).

Similarly,

Zxx=Z″(η)ρx2r2+Z′(η)ρxxr,

and therefore

(18.5)Zxx=O(r−2).

Now

r=sf1/8

gives

r˙=18sf−7/8s˙f.

With

|s˙f|≍sf−1/8,

we obtain

|r˙|≍sf−1.

Consequently,

(18.6)|r˙|r≍sf−9/8.

Thus

(18.7)Zt=O(sf−9/8),Zx=O(sf−1/8),Zxx=O(sf−1/4).

The time derivative is indeed the most singular cutoff derivative.

18.3. Matching order

Suppose on the moving transition

(18.8)DN=O(sfμ).

Then the new cutoff contribution to ∂t(A2) is

(18.9)ZtDN=O(sfμ−9/8).

The transition width in ρ is

Δρ≍r≍sf1/8.

For the part of the construction following the parabolic relation, its corresponding s-width scales as

Δs≍r2≍sf1/4.

Therefore

(18.10)|TN(t)|≍sf3/8

up to the fixed y-length.

Squaring (18.9),

|ZtDN|2=O(sf2μ−9/4).

Spatial integration gives

(18.11)∫TN(t)|ZtDN|2dV=O(sf2μ−15/8).

Finally,

dt≍sf1/8dsf,

so

(18.12)∫∫TN(t)|ZtDN|2dVdt≲∫0sf2μ−7/4dsf.

Therefore this term is integrable exactly when

2μ−74>−1,

or

(18.13)μ>38.

So our previous 3/8 threshold is confirmed.

18.4. Spatial cutoff derivative is easier

The first spatial commutator satisfies

ZxDN=O(sf−1/8sfμ)=O(sfμ−1/8).

Squaring and integrating in spacetime gives

∫|ZxDN|2dVdt≲∫0sf2μ−1/4sf3/8sf1/8dsf=∫0sf2μ+1/4dsf.

This is finite for

μ>−58.

In particular every positive matching order is sufficient.

Thus

ZtDN is much more restrictive than ZxDN.

18.5. Second spatial derivative

For viscosity we also encounter

ZxxDN=O(sfμ−1/4).

Its squared spacetime integral behaves like

∫0sf2μ−1/2sf3/8sf1/8dsf=∫0sf2μdsf.

Therefore

ZxxDN∈Lt,x2

for every

μ>−12.

Again, μ>3/8 automatically handles it.

But differentiating (18.1) twice also produces

2ZxDx,

so we need a derivative matching condition, not just a zeroth-order one.

If

Dx=O(sfμx),

then

ZxDx=O(sfμx−1/8),

which is spacetime L2 provided

μx>−58.

Flat matching can make this substantially better.

18.6. Choose a stronger matching order

Instead of using the minimal condition

μ>38,

we can design the interpolation with

(18.14)DN=O(sf1/2+δ),δ≥0.

The simplest useful target is

(18.15)DN=O(sf1/2).

Then

ZtDN=O(sf−5/8),

and its squared spacetime integral behaves like

∫0sf−3/4dsf<∞.

Thus 1/2-order matching gives a genuine margin beyond the critical 3/8.

18.7. But we need to check whether Acore2 naturally has that scale

At the moving front,

ρ≍r≍sf1/8.

The outer reconstruction is

Acore2=2ρb.

With

b≍sf−1/8,

we get

(18.16)Acore2≍sf1/8sf−1/8≍1.

This is extremely useful.

The moving interface is therefore exactly where the square-root reconstruction has returned to an O(1) transverse amplitude. That agrees with the rationale for this paper’s original release: matching a bounded nonzero A∗ requires q≍s1/8.

So there is no leading-order amplitude mismatch that must diverge.

We may choose the regular inner field so that

(18.17)Areg,N2=Acore,N2+O(sf1/2)

at the moving interface.

This is a matching condition on the regular reconstruction, not an arbitrary cancellation of a singular field.

18.8. Exact incompressibility

Now reconstruct C=uz from

∇⋅u=0.

Since

ux=bA,uy=A,

we impose

(18.18)∂x(bA)+Ay+Cz=0.

If the moving local construction remains y-independent,

Ay=0,

so

(18.19)Cz=−∂x(bA).

Since

s=t−z−θ∗,Cz=−Cs,

we get

(18.20)Cs=∂x(bA).

Thus we should not independently interpolate C.

We interpolate only A2, then solve (18.20). This retains the strongest feature of the current paper’s release construction: incompressibility is exact by construction. The existing paper explicitly follows this procedure.

18.9. Moving-cutoff contribution to Cs

Because

bA=b(A2)−1/2,

and bx=0 in the distinguished local core,

∂x(bA)=−b2A3∂x(A2).

From

A2=Areg2+ZD,

we have

∂x(A2)=∂xAreg2+ZxD+ZDx.

The new cutoff term is therefore

(18.21)Cscut=−b2A3ZxD.

In the moving transition,

A≍1,b≍sf−1/8.

If

D=O(sfμ),

then

(18.22)Cscut=O(sf−1/8sf−1/8sfμ)=O(sfμ−1/4).

With our proposed

μ=12,

this gives

(18.23)Cscut=O(sf1/4).

So the moving interpolation does not introduce a singular continuity correction. It actually vanishes.

That is encouraging.

18.10. Horizontal transport remainder

Recall this paper’s compact remainder

RN=2bNXhCN+bNCN,t+CN,zCN,

with

Xh=A∂x+bA∂y.

In the local y-independent reconstruction,

Cy=0,

so

XhC=ACx.

this paper’s current core estimate uses precisely this simplification.

The subtlety is now

Ct+Cz.

In the old stationary reconstruction, if C=C(x,s) with

s=t−z−θ∗,

then

Ct+Cz=0.

But our new A depends additionally on

r(t)=ρmin(t).

Consequently the reconstructed C will generally have the form

C=C(x,s,r(t)).

Then

Ct=Cs+Crr˙,Cz=−Cs,

so

(18.24)Ct+Cz=Crr˙.

This is the genuinely new remainder created by the moving front.

It cannot be set to zero.

18.11. Estimate the new moving-front remainder

We therefore acquire

(18.25)Rmov=bCrr˙C.

To preserve the scalar Lambert mechanism we need

(18.26)|Rmov|Γb2=|Crr˙|CΓb⟶0.

At k=3,

Γ≍sf−1,b≍sf−1/8,

so

Γb≍sf−9/8.

Also

|r˙|≍sf−1.

Hence

(18.27)|Rmov|Γb2≍|Cr|sf1/8.

Thus the requirement is surprisingly mild:

(18.28)Cr=o(sf−1/8).

In particular, if the flat matching yields

Cr=O(1),

then

(18.29)|Rmov|Γb2=O(sf1/8)→0.

So the moving-front time dependence does not automatically destroy the Lambert quadratic source.

18.12. An important point

We impose first-order flat matching of the inner and outer reconstructions across the moving collar. In particular, at the matching interface,

Cin=Cout,∂ρCin=∂ρCout,

and the interpolation cutoff H is chosen flat at its endpoints,

H′(0)=H′(1)=0.

More generally, when higher derivatives enter the reconstruction we impose

H(j)(0)=H(j)(1)=0,1≤j≤m.

The matching amplitudes are chosen so that the collar interpolation satisfies the uniform estimate

|Cr|≤M

with M independent of the terminal parameter sf.

18.13. Revised list of conditions

We now have a fairly precise moving-front theorem to aim for. It would require:

rN(t)≍sf(t)1/8,
DN=Acore,N2−Areg,N2=O(sf1/2),

together with suitable first-derivative matching,

AN≥a−>0throughout the moving transition,

and

CN,r=O(1)

or, more generally,

CN,r=o(sf−1/8).

Under those assumptions we obtain simultaneously:

∫T∗‖∂zux(t)‖22dt<∞,

and the other first derivatives considered above are no worse;

∇⋅u=0

exactly, because C is reconstructed from continuity;

and

|Rmov|Γb2→0.
18.13.1. The remaining problem is viscosity across the moving transition

This is the next calculation we should do before declaring the redesign successful. The viscosity operator is

V[b,C]=ν[Δb+bCΔC−2∇C⋅∇b].

The moving interpolation changes C, hence ΔC. We therefore have to calculate the new Cxx, Cxz, and Czz terms generated by Zx, Zxx, Zt, determine their powers of sf, and compare them with

νΔb∼9νB64s−17/8.

If all moving-interface viscosity terms are lower order while the first-contact ordering remains intact, then we would have a coherent candidate for replacing the current fixed-core/release section. If not, the proposed moving cutoff would need another modification.

19. High-Order Extension of the Square-Root Amplitude

**Transition.**The commutator audit shows that low-order matching is insufficient. We therefore replace the naive moving interpolation by a flat, high-order extension of A2. This resolves the amplitude-side cutoff defect, but exact continuity and periodicity expose a separate solvability condition. In addition, keeping the singular outer b unchanged inside the front leaves a derivative obstruction; these two residual issues determine the next redesign.

The clean way is to make Areg,N an explicit high-order extension of the outer square-root amplitude across the moving interface, rather than choosing an unrelated bounded function.

Let

FN(x,s):=Acore,N2(x,s)=−2χN(x)bN(s),

and let the moving interface be

ρN=rN(t),rN(t)=κNsf(t)1/8.

On this interface, FN=O(1), because bN≍sf−1/8 and ρN≍sf1/8. This is exactly the scale at which this paper’s existing release matches a bounded nonzero transverse amplitude.

The useful observation is that we do not need to approximate FN throughout the inner region. We need it to agree sufficiently accurately only in the moving transition where derivatives of the cutoff are supported.

19.1. Use a flat matching collar

Choose fixed smooth functions Φ, Ψ∈C∞(R) such that

Φ(η)={0,η≤12,1,η≥1,Ψ=1−Φ,

with Φ flat at η=1:

Φ(j)(1)=0,j≥1.

Set

ηN=ρNrN(t).

Choose a positive smooth periodic background amplitude

aN(x,t)≥a−>0.

Now define

(19.1)Areg,N2=Φ(ηN)FN+Ψ(ηN)aN2.

This has two immediate properties:

(19.2)Areg,N2=FNfor ρN≥rN(t),

and

(19.3)Areg,N2=aN2for ρN≤12rN(t).

Thus it is regular and strictly positive at the zero set while agreeing exactly with the square-root reconstruction at the outer edge.

But there is an even cleaner construction for our purpose.

19.2. Put the moving interpolation outside the matching collar

Recall that we had written

AN2=Areg,N2+Z(ηN)DN,DN=FN−Areg,N2.

Choose Z′ to be supported only where

ηN≥1.

But by (19.2),

Areg,N2=FN

throughout this support. Therefore

DN=0

wherever

Zx,Zt,Zxx

are nonzero.

Consequently,

(19.4)ZtDN=ZxDN=ZxxDN=0.

Even better,

DN,x=DN,xx=0

on the active moving-cutoff support.

Hence the three estimates we previously sought,

DN=O(sf1/2),
DN,x=O(sf3/8),
DN,xx=O(sf1/4),

hold there trivially—in fact with exact zero.

This removes the moving commutator rather than merely estimating it.

19.3. Simplify further: we no longer need two cutoffs

The preceding construction shows that the most natural definition is simply

(19.5)AN2(x,t)=Φ(ρN(x)rN(t))[−2χN(x)bN(s)]+[1−Φ(ρN(x)rN(t))]aN2(x,t).

Then:

ρN≥rN(t)⟹AN2=−2χNbN,

whereas

ρN≤12rN(t)⟹AN2=aN2.

The transition occurs in

12rN(t)<ρN<rN(t).

So we can dispense with the artificial decomposition

Areg2+ZD.

Equation (19.5) itself is the regularized moving reconstruction.

19.4. Positivity

This is immediate if we stay on the selected branch where

−χNbN>0

and choose

aN2≥a−2>0.

Since

0≤Φ≤1,

(19.5) is a convex combination of positive quantities:

AN2≥min{2|χN|bN,aN2}.

In the transition,

ρN≍rN≍sf1/8

and

bN≍sf−1/8,

so

2|χN|bN≍1.

Hence

(19.6)AN≥A−>0

uniformly through the moving collar, provided the constants are uniformly controlled.

19.5. A note on the moving collar

Definition 14.1 (Moving matching collar). Let

ρ(x)=distT⁡(x,I3),

and let rN(t)>0 denote the moving coterminal front. For a positive collar width ℓN(t), define

CN(t)={(x,y,z)∈T3:rN(t)≤ρ(x)≤rN(t)+ℓN(t)}.

The normalized collar coordinate is

ηN(x,t)=ρ(x)−rN(t)ℓN(t)∈[0,1].

The inner and outer reconstructions are joined in CN(t) by a flat cutoff H∈C∞(R), with

H=0for η≤0,H=1for η≥1,

and

H(j)(0)=H(j)(1)=0,j≥1.

Thus, for example,

CN=(1−H(ηN))Cin,N+H(ηN)Cout,N.

The matching is said to be flat to order m if

∂ρj(Cout,N−Cin,N)=O(ℓNm+1−j),0≤j≤m,

uniformly in the collar as t↑T∗,N.

Inside the regular region,

AN=aN≥a−.

Thus division by AN is removed there.

19.6. Periodicity

Take

χN(x)

periodic, as done earlier, and define the distance variable locally through the periodic profile itself:

ρN=|χN(x)|

rather than using a nonperiodic Euclidean coordinate globally.

Then

Φ(ρNrN)

is periodic in x.

If aN(x,t) is periodic, equation (19.5) is periodic.

Therefore

AN(x+1,t)=AN(x,t)

in the normalized torus convention, or 2π-periodically if that is the convention used in the relevant section.

19.7. Smoothness and the absolute-value issue

There is one subtle point. The function

|χN|

is not smooth through a simple zero.

But in the innermost region

ρN≤12rN

we have

Φ=0.

Therefore the reconstructed amplitude is simply

AN2=aN2

in an open neighborhood of χN=0.

All derivatives of Φ vanish there.

Consequently the nonsmoothness of |χN| at its zero never enters a differentiated active cutoff.

Alternatively, work separately on the two signed charts

qN±=∓χN>0,

which is even cleaner analytically.

Thus the construction can be made C∞.

19.8. Exact continuity

Now define

uy,N=AN,ux,N=bNAN.

Do not prescribe CN=uz,N independently.

Instead solve

(19.7)CN,z=−∂x(bNAN)−∂yAN.

In the local y-independent construction,

AN,y=0,

so

(19.8)CN,z=−∂x(bNAN).

Hence

∇⋅uN=∂x(bNAN)+∂yAN+CN,z=0.

Therefore

(19.9)∇⋅uN=0

identically.

This preserves the existing strategy of obtaining C from exact incompressibility.

19.9. But periodicity of C creates a genuine solvability condition

This point should not be skipped.

Solving

Cz=−∂x(b/A)−Ay

locally is straightforward. Solving it periodically in z requires

(19.10)∫Tz[∂x(bA)+Ay]dz=0.

Otherwise integrating Cz around the z-circle does not return to the same C.

Thus the exact periodic reconstruction is

(19.11)CN(x,y,z,t)=C0,N(x,y,t)−∫z0z[∂x(bNAN)+AN,y]dζ,

provided (19.10) holds.

This is a real global compatibility condition. The moving cutoff does not automatically prove it.

So the construction is explicit locally, but the torus-periodic C-closure still requires verification.

19.10. Initial data

Choose the moving modification to be inactive initially.

For example, introduce a smooth activation function

ΛN(t)∈[0,1]

with

ΛN(t)=0

on an initial interval and

ΛN(t)=1

only after the solution enters the terminal regime.

Then replace (19.5) by

(19.12)AN2=(1−ΛN)Apre,N2+ΛN[Φ(ρNrN)(−2χNbN)+(1−Φ)aN2].

At t=0,

ΛN(0)=0,

hence

(19.13)AN(x,0)=Apre,N(x,0).

So the finite-N initial data are preserved exactly.

However, this introduces another time cutoff ΛN′(t). It should be supported strictly away from T∗, where all quantities are regular. Then its derivatives are bounded and do not affect terminal asymptotics.

19.11. The important correction to our previous D-strategy

There is a conceptual improvement here.

Previously we were trying to prove

D=O(sf1/2),Dx=O(sf3/8),Dxx=O(sf1/4).

But an explicit flat extension lets us do better:

(19.14)D=Dx=Dxx=0

at the interface where the moving transition joins the square-root field.

This is preferable to merely achieving the power-law estimates.

It means the moving transition can be designed so that there is no extra release/core commutator at its outer flat endpoint, paralleling the flat-endpoint philosophy already used by this paper.

19.12. What this does to the H1 estimate

The outer square-root field starts at

ρ≳r(t)≍sf1/8.

Thus the previous estimate remains

‖∂zux(t)‖22≲Cr(t)−9/4.

Since

r(t)≍sf1/8,
‖∂zux(t)‖22≲Csf−9/32.

With

dt≍sf1/8dsf,

we get

(19.15)∫T∗‖∂zux(t)‖22dt≲∫0sf−5/32dsf<∞.

Inside

ρ<12r(t),

A=aN is regular and bounded away from zero, so

ux=baN.

But here we encounter the next decisive issue: if the same singular

bN(s)≍s−1/8

is retained throughout this inner regular region, then regularizing A alone does not necessarily regularize ux’s z-derivative. In particular,

∂z(baN)

still contains bz∼s−9/8.

So we have successfully constructed an explicit positive, smooth, periodic Areg,N, but that by itself is not enough to establish the full Leray bound.

The inner reconstruction must also specify what happens to bN.

That leaves us with the sharper structural choice:

either the moving front truncates/modifies the singular bN inside ρ<rN(t),

or one must prove that the spacetime geometry of the inner region makes the retained bN-derivatives integrable.

So next I will do this. It is more fundamental than further tuning Areg,N, because simply keeping

bN(s)=BNs−1/8

everywhere while replacing A by a bounded aN may move, rather than eliminate, the dissipation obstruction.

20. Why the Inner b-Field Must Also Be Regularized

Transition. Regularizing only the square-root amplitude does not regularize ux=b/A if the singular b-profile is retained in the inner region. This section therefore regularizes b as well and audits direct matching. The calculation reveals that a separate b-matching scale or a moving cutoffwith nontrivial time derivative produces a new mismatch. Rather than accumulating independent cutoffs, the next section replaces them by one geometric front with an exact-overlap collar.

Proceeding, the inner-region calculation shows something stronger: regularizing only A is insufficient if the singular law b∼s−1/8 is retained unchanged inside the moving front. We therefore need a moving regularization of b as well, while matching it exactly to the Lambert-driven outer solution.

20.1. Why b must also change inside the front

In the regular inner region we proposed

AN=aN,aN≥a−>0.

If we nevertheless retain

bN(s)=BNs−1/8,

then

ux,N=bNaN.

Even taking aN locally constant,

∂zux,N=BN8aNs−9/8,

and therefore

|∂zux,N|2≍s−9/4.

So replacing the vanishing A by a positive amplitude does not remove the longitudinal singular derivative.

The clean redesign is therefore

(breg,Areg)⟶moving flat transition⟶(bLam,Acore).

This is a substantive modification of the current paper, not something already proved.

20.2. A natural regularized b

Let

sf(t)>0

denote the moving terminal scale and define

(20.1)breg,N(s,t)=BN(s+sf(t))−1/8.

For every t<T∗,N,

breg,N(0,t)=BNsf(t)−1/8<∞.

As

t↑T∗,N,sf(t)↓0,

we recover formally

breg,N(s,t)⟶BNs−1/8

for every fixed s>0.

Thus the singular profile is obtained only in the terminal limit.

Its derivatives are

breg,s=−BN8(s+sf)−9/8,

and

(20.2)breg,ss=9BN64(s+sf)−17/8.

In particular,

(20.3)breg,ss>0.

So the favorable convexity needed by the viscous term is retained.

20.3. But direct matching at s∼sf has an O(1) mismatch

The outer Lambert profile is

bout=BNs−1/8.

At s=sf,

bout(sf)=BNsf−1/8,

whereas

breg(sf)=BN(2sf)−1/8=2−1/8BNsf−1/8.

The relative mismatch is

1−2−1/8,

which does not vanish.

So (20.1) is an excellent inner regularization, but we should not glue it directly to the outer Lambert profile.

We need a matching collar.

20.4. Introduce a separate b-matching scale

Let

(20.4)σf(t)=sf(t)γ,0<γ<1.

Since sf↓0,

σf≫sf.

At

s≍σf,

we have

sfs≍sf1−γ→0.

Now expand:

(s+sf)−1/8=s−1/8(1+sfs)−1/8.

Hence

breg=bout[1−18sfs+O(sf2s2)].

Therefore

(20.5)bout−breg=O(sfs−9/8).

At s≍sfγ,

(20.6)bout−breg=O(sf1−9γ8).

This mismatch tends to zero provided

1−9γ8>0,

i.e.

(20.7)γ<89.

So there is a large admissible range.

20.5. Match derivatives too

Differentiate (20.5):

bout,s−breg,s=O(sfs−17/8),

so at s=sfγ,

(20.8)Δbs=O(sf1−17γ8).

For this absolute difference to vanish we require

(20.9)γ<817.

For the second derivative,

bout,ss−breg,ss=O(sfs−25/8),

and hence

(20.10)Δbss=O(sf1−25γ8).

Absolute C2 matching would require

(20.11)γ<825.

This suggests choosing a comfortably smaller value, for example

(20.12)γ=14.

Then

σf=sf1/4,

and

Δb=O(sf23/32),
Δbs=O(sf15/32),
(20.13)Δbss=O(sf7/32).

All three vanish.

That is a particularly useful matching scale.

20.6. Explicit smooth b-interpolation

Choose Y∈C∞ with

Y(ξ)=0(ξ≤1),Y(ξ)=1(ξ≥2),

flat at both endpoints, and put

ξ=sσf(t)=ssf(t)1/4.

Define

(20.14)b~N=(1−Y(ξ))breg,N+Y(ξ)bout,N.

Then

s≤σf⟹b~N=breg,N,

while

s≥2σf⟹b~N=bout,N.

Thus the Lambert profile remains exactly unchanged in the outer region.

20.7. Spatial cutoff derivatives

Because

ξ=sσf

and sz=−1,

Yz=−Y′(ξ)σf=O(sf−1/4),

and

(20.15)Yzz=O(sf−1/2).

But these multiply the small mismatch (20.13).

Thus

YzΔb=O(sf−1/4sf23/32)=O(sf15/32),

and

YzzΔb=O(sf−1/2sf23/32)=O(sf7/32).

Similarly,

YzΔbs=O(sf−1/4sf15/32)=O(sf7/32).

All vanish.

Therefore the second spatial derivative of the interpolated b differs from the corresponding regular/outer derivatives by vanishing commutator terms:

(20.16)b~zz=(1−Y)breg,zz+Ybout,zz+O(sf7/32).

Since both principal second derivatives are positive,

breg,zz>0,bout,zz>0,

we obtain, sufficiently near the terminal regime,

(20.17)b~zz>0.

So the favorable longitudinal viscous sign can survive this interpolation.

20.8. Important time derivative

Now

σf=sf1/4,

so

σ˙fσf=14s˙fsf.

Under

|s˙f|≍sf−1/8,

we have

(20.18)|σ˙f|σf≍sf−9/8.

On the transition support ξ=O(1),

Yt=O(sf−9/8).

Multiplying by the mismatch,

(20.19)YtΔb=O(sf−9/8sf23/32)=O(sf−13/32).

This diverges pointwise, but much more slowly than the Lambert quadratic source.

At s≍σf=sf1/4,

b2≍s−1/4≍sf−1/16.

Also

Γ≍1s≍sf−1/4.

Hence

(20.20)Γb2≍sf−5/16=sf−10/32.

Here we discover something important:

sf−13/32

is actually more singular than

sf−10/32.

Thus the simple time-dependent interpolation (20.14) is not automatically a lower-order perturbation of the Lambert source.

So γ=1/4, although excellent for spatial C2 matching, does not automatically close the evolution equation.

20.9. General γ reveals the obstruction

For general

σf=sfγ,

we found

Δb=O(sf1−9γ/8).

Meanwhile

Yt=O(sf−9/8),

so

(20.21)YtΔb=O(sf−1/8−9γ/8).

At s≍sfγ,

(20.22)Γb2≍s−5/4=sf−5γ/4.

The ratio is therefore

(20.23)|YtΔb|Γb2=O(sf−1/8−9γ/8+10γ/8)=O(sf(γ−1)/8).

But for every

0<γ<1,

we have

γ−18<0.

Hence

(20.24)|YtΔb|Γb2→∞.

This is a decisive negative result for this particular moving s-cutoff design.

20.10. What this teaches us

We should not add an independently time-moving interpolation of b through

Y(ssf(t)γ).

Its time derivative produces a commutator that can dominate the Lambert source.

That does not invalidate the moving-ρ idea. It tells us that the scalar amplitude b should not be regularized by a second independently moving s-cutoff.

The better possibility is to regularize b using the same geometric moving variable as A:

η=ρr(t),r(t)≍sf1/8,

and arrange the inner breg to match bout to sufficiently high order on the support of ηt.

That returns us to the useful condition

mismatch×r˙r,

but now we can force the mismatch to vanish exactly at a flat geometric endpoint.

So the calculation has narrowed the redesign substantially:

one moving geometric interface in ρ, not independent moving cutoffs in both ρ and s.

The next construction should therefore define the pair

(Areg,N,breg,N)

jointly across the same flat η=ρ/rN(t) collar, with exact jet matching to the Lambert pair at its outer endpoint. That is the route that avoids the YtΔb obstruction we just found.

21. One Geometric Front and an Exact-Overlap Collar

Transition. The preceding two-scale b interpolation solves one defect by creating another. We therefore consolidate the geometry: one front separates a regular inner field from the exact Lambert outer field, while an overlap collar is chosen so that no moving cutoff differentiates a mismatch of leading order. This removes the earlier cutoff-time obstruction. What remains is to construct a scalar regularized phase whose derivatives and viscous terms are controlled, which is the purpose of the next section.

The previous calculation tells us exactly what to avoid: we should not introduce a second independently moving s-cutoff for b. Instead, use one geometric front

ηN=ρNrN(t),rN(t)≍sf(t)1/8,

and construct AN and bN jointly across the same collar.

There is, however, an important structural constraint: merely making a cutoff flat at its endpoint does not make the mismatch vanish throughout the support of its derivatives. The cleanest construction is therefore to build an overlap collar in which the regular and Lambert fields agree identically, so every moving-cutoff commutator vanishes there.

21.1. Three-region construction

Let

rN(t)=κNsf(t)1/8,ηN=ρNrN(t).

Use three regions:

(21.1)IN(t)={ηN<1},regular inner region,ON(t)={1<ηN<2},overlap/matching collar,EN(t)={ηN>2},Lambert outer region.

The essential design requirement is

(21.2)(breg,N,Areg,N)=(bLam,N,Acore,N)on ON(t).

Not merely asymptotically. Exactly.

Then a cutoff supported inside ON creates no mismatch commutator at all.

21.2. Joint cutoff

Choose

Z∈C∞(R),Z(η)=0(η≤1),Z(η)=1(η≥2).

Define

(21.3)bN=(1−Z)breg,N+ZbLam,N,

and

(21.4)AN2=(1−Z)Areg,N2+ZAcore,N2.

Ordinarily, differentiating (21.3) produces

Zt(bLam−breg),Zx(bLam−breg),

and similarly for A2.

But (21.2) gives on supp⁡Z′,

bLam−breg=0,

and

Acore2−Areg2=0.

Therefore

(21.5)Zt(bLam−breg)=Zx(bLam−breg)=Zxx(bLam−breg)=0,

and likewise

(21.6)Zt(Acore2−Areg2)=Zx(Acore2−Areg2)=Zxx(Acore2−Areg2)=0.

This removes the obstruction found in the previous calculation.

21.3. But can the regular field agree with the singular field in the overlap?

Yes, because for every fixed preterminal time

t<T∗,N,

the overlap is separated from the singular set:

ρN≥rN(t)>0.

Thus the Lambert reconstruction is smooth there.

The singularity occurs only as

rN(t)↓0.

So there is no contradiction between:

breg=bLamfor ρ≥rN(t)

and having a genuinely regular extension for

0≤ρ<rN(t).

This converts the problem into an extension problem.

21.4. Explicit inner extension by Taylor jets

Fix the inner interface

ρ=r(t).

For brevity suppress N.

Define the boundary jets

Bj(s,t)=∂ρjbLam(ρ,s,t)|ρ=r(t),

and

Fj(s,t)=∂ρjAcore2(ρ,s,t)|ρ=r(t).

For the distinguished scalar construction,

bLam=bLam(s),

so actually

(21.7)Bj=0,j≥1.

That simplifies the b-extension dramatically.

Choose a smooth shape function

H(ξ),ξ=ρr(t),

with

H(ξ)=0near ξ=0,

and

H(ξ)=1near ξ=1.

Then define

(21.8)breg(ρ,s,t)=bin(s,t)+H(ρr(t))[bLam(s)−bin(s,t)].

If H≡1 on an interval

1−δ<ρr<1,

then

(21.9)breg=bLam

on an open collar immediately inside the interface.

All jets consequently agree there.

The outer cutoff Z can be supported entirely inside this exact-agreement collar.

21.5. What should bin be?

Here we must not repeat the previous mistake. A choice such as

bin=B(s+sf)−1/8

has strong time dependence and creates a new evolution error.

A cleaner possibility is to define the inner value by the front value of the outer field:

(21.10)bin(t)=bLam(sf(t)).

Then bin is spatially constant in the innermost region.

Consequently

(21.11)∇bin=0,Δbin=0.

This immediately removes the inner spatial H1 singularity coming from bz.

At every t<T∗,

bin(t)<∞.

But as

t↑T∗,
(21.12)bin(t)∼Bsf(t)−1/8→∞.

Thus we obtain finite preterminal spatial gradients while retaining terminal amplitude growth.

This is a much better candidate than retaining s−1/8 inside the front.

21.6. Corresponding Ain

In the innermost region choose

(21.13)Ain=aN>0

constant locally, or a smooth positive periodic function with uniformly bounded derivatives.

Then

uxin=bin(t)aN,uyin=aN.

Spatially,

∇uxin=0

if aN is constant in the local chart.

Hence the inner region contributes no singular spatial dissipation.

21.7. Extend A2 to the outer field

Set

Fcore=Acore2=−2χbLam.

Choose another smooth collar function K with

K=0near ρ=0,
K=1on an open neighborhood of ρ=r(t).

Define

(21.14)Areg2=aN2+K(ρr(t))[Fcore−aN2].

Near the interface,

K=1,

so

(21.15)Areg2=Fcore

exactly.

Near the center,

K=0,

so

(21.16)Areg2=aN2>0.

Thus positivity and exact outer matching are both achieved.

21.8. Collapse the redundant outer interpolation

At this point something useful happens.

Because both regular extensions already equal the outer fields on an open collar, we do not actually need another physical interpolation there.

We can define directly

(21.17)bN={breg,N,ρN<rN(t),bLam,N,ρN≥rN(t),

and

(21.18)AN2={Areg,N2,ρN<rN(t),−2χNbLam,N,ρN≥rN(t).

Because the two definitions agree on an open neighborhood of the interface, these piecewise formulas define a single smooth function.

There is no derivative jump.

There is no distributional interface term.

There is no ZtD commutator.

This is cleaner than retaining an explicit Z in the final formulas.

21.9. The crucial H1 decomposition

Now split

T3=IN(t)∪CN(t)∪EN(t),

where:

  • IN: constant/smooth inner region;

  • CN: extension collar;

  • EN: outer Lambert region.

For the inner region,

(21.19)‖∇uN‖L2(IN)2=O(1)

under the locally constant choice.

For the outer region we already found

(21.20)‖∂zux‖L2(EN)2≲rN(t)−9/4.

With

rN≍sf1/8,
rN−9/4=sf−9/32.

Since

dt≍sf1/8dsf,
(21.21)∫T∗‖∂zux‖22dt≲∫0sf−5/32dsf<∞.

Thus the remaining issue is entirely concentrated in the extension collar.

21.10. Collar derivative scaling

The collar has width

Δρ≍r.

Therefore a generic interpolation between two O(1)-different amplitudes produces

∂ρA=O(r−1),

and hence

|∂ρA|2=O(r−2).

Multiplying by the collar width r,

(21.22)∫collar|∂ρA|2dρ=O(r−1).

The relevant s-width associated with the moving front is of order

r2.

Therefore

(21.23)‖∂ρA‖L2(C(t))2=O(r).

This actually tends to zero.

For b, however,

bin≍sf−1/8.

At

r=sf1/8,

this is

bin≍r−1.

If the collar interpolated an O(r−1) mismatch over width r, we would obtain

bρ=O(r−2),

which would be dangerous.

But our construction deliberately avoids this: bin must be chosen so that it matches the outer b at the collar to high order.

Since the outer b=b(s) has no ρ-dependence,

bρout=0.

Therefore the optimal choice is simply

breg(ρ,s,t)=bLam(s)

through the entire matching collar, with the flattening to bin(t) occurring farther inside.

That separates the singular outer field from the moving-interface derivatives.

21.11. We have therefore obtained a nested geometry

The revised construction naturally has four regions:

(21.24)I:0≤ρ≤θ1r,b=bin(t),A=aN,II:θ1r<ρ<θ2r,inner flattening layer,III:θ2r<ρ<r,b=bLam(s),A2=−2χb,IV:ρ≥r,outer Lambert reconstruction,

with

0<θ1<θ2<1.

Region III is an exact overlap.

Therefore the moving boundary at ρ=r(t) produces no commutator whatsoever.

The only nontrivial transition is Region II, which is strictly inside the moving front.

21.12. A new issue appears—and it is the correct one

Region II must interpolate between

bin(t)

and

bLam(s).

For that transition to have small gradients, these two values must be close on Region II.

That requires

(21.25)s≈sf(t)

inside the flattening layer.

In other words, the inner regularization cannot be defined purely by ρ/r(t). It must use the geometry of the moving front in the (ρ,s)-plane.

This is actually more natural.

Define a front function

ΣN(ρ,t)

such that the singular outer field is replaced by its value at

s=ΣN(ρ,t)

whenever

s<ΣN(ρ,t).

A particularly simple regularization is

s^N=s+ΣN(ρ,t)

or, better for exact outer matching, a smooth maximum

(21.26)s^N=smax⁡{s,ΣN(ρ,t)}.

Then define

(21.27)bN=BNs^N−1/8.

In the outer region,

s>ΣN⟹s^N=s,

so

bN=BNs−1/8

exactly.

Inside,

s<ΣN⟹s^N≈ΣN>0,

so bN remains finite.

This is superior to independently interpolating the amplitude itself.

21.13. The problem has therefore reduced to constructing one scalar regularized phase

Instead of independently regularizing A and b, define

s^N(ρ,s,t)>0

such that

s^N=s

on the outer Lambert region,

s^N≥sf(t)>0

inside the moving front, and its first and second spatial derivatives have controlled spacetime norms.

Then set

bN=BNs^N−1/8,

followed by

AN2=−2χNbN

only where this quantity stays safely positive/nondegenerate, and use the positive inner A-extension closer to χN=0.

This is a much more unified redesign.

The next decisive calculation is therefore to choose an explicit smooth-max phase s^N, compute

s^z,s^zz,s^x,s^xx,s^t,

and then determine whether

bN=BNs^N−1/8

simultaneously preserves the favorable viscosity, finite spacetime H1, and the Lambert source asymptotically. That is now a single, well-defined calculation rather than several independent cutoff problems.

22. Unified Regularized Phase and the Convexity Constraint

Transition. The one-front architecture reduces the reconstruction problem to a single regularized phase. We now test smooth-max and equivalent regularizations, including their derivative and viscosity signs. The convexity calculation shows which pointwise sign requirements are too strong and which estimates are actually needed. After this structural reduction, the remaining question is quantitative: whether the inner interpolation itself has finite Lt2Hx1 cost.

Proceeding with the unified regularized-phase construction, there is a useful candidate, but the derivative audit exposes an important constraint on how the phase floor may depend on ρ.

Take

α=18,bN=BNs^N−α.

We want s^N=s exactly in the outer Lambert region, while s^N>0 in the regularized inner region.

22.1. A smooth maximum with exact outer recovery

Choose h∈C∞(R) satisfying

h(ξ)=0(ξ≥1),h(ξ)=−ξ(ξ≤−1),

with

0≤1+h′(ξ)≤1.

For a transition width δN>0, define

(22.1)s^N=s+δNh(s−ΣN(ρ,t)δN).

Write

ξ=s−Σδ.

Then:

  • if s≥Σ+δ, h=0, so

    s^=s;
  • if s≤Σ−δ, h=−ξ, so

    s^=Σ.

Thus the outer Lambert profile is recovered exactly, rather than asymptotically.

For the moment take δ proportional to Σ:

(22.2)δ=εΣ,0<ε<1.

Then throughout the transition

s^≍Σ.

22.2. Derivatives of the regularized phase

It is convenient to write the smooth maximum abstractly as

s^=Mδ(s,Σ).

Define

m:=∂ss^.

By construction,

0≤m≤1,

with

m=1

in the outer region and

m=0

in the inner region.

For a fixed-time spatial derivative,

sz=−1,sx=0.

If Σ=Σ(ρ,t), then

(22.3)s^z=−m

provided Σ has no z-dependence.

Hence

|s^z|≤1.

A second derivative costs one inverse transition width:

(22.4)s^zz=O(δ−1)=O(Σ−1)

inside the smoothing layer.

For x,

(22.5)s^x=(1−m)Σx,

and schematically

(22.6)s^xx=(1−m)Σxx+O(Σx2Σ).

These formulas immediately show that the geometry of Σ(ρ,t) is decisive.

22.3. Derivatives of b=Bs^−1/8

For general α>0,

b=Bs^−α.

Therefore

(22.7)bj=−αBs^−α−1s^j,

and

(22.8)bjj=α(α+1)Bs^−α−2s^j2−αBs^−α−1s^jj.

At α=1/8,

(22.9)bjj=9B64s^−17/8s^j2−B8s^−9/8s^jj.

In the exact outer region,

s^=s,s^z=−1,s^zz=0,

so

(22.10)bzz=9B64s−17/8,

exactly reproducing the favorable term used earlier.

22.4. Viscosity inside the smooth-max layer

In the smoothing layer,

s^≍Σ,

and

s^zz=O(Σ−1).

Hence both pieces of (22.9) have the same possible scale:

s^−17/8s^z2=O(Σ−17/8),

and

s^−9/8s^zz=O(Σ−17/8).

Therefore

(22.11)bzz=O(Σ−17/8).

But there is a crucial sign issue: unlike the exact outer region, the second term

−B8s^−9/8s^zz

need not be positive.

So an arbitrary smooth maximum does not automatically preserve

bzz>0.

We need a structural condition on the smoothing function.

22.5. Exact convexity criterion

From (22.9),

bzz≥0

iff

98s^z2s^≥s^zz.

Thus the precise condition is

(22.12)s^s^zz≤98s^z2.

This is much better than estimating the two terms separately.

Equivalently,

(22.13)∂zz(s^−1/8)≥0.

So we should design the regularized phase so that b=Bs^−1/8 itself is convex in the distinguished z-direction.

22.6. An even cleaner construction: regularize b, then define s^

This observation suggests reversing the construction.

Rather than picking an arbitrary smooth maximum for s^ and hoping its composition is convex, choose a smooth convex regularization of

s−1/8

directly.

Let

βf(t)>0

be the finite inner amplitude, and choose a smooth convex function

Bf(s)

such that

Bf(s)=Bs−1/8

for s≥smatch(t),

while

Bf(s)=βf(t)

in the innermost regular region.

Then set

(22.14)bN=Bf(s),s^N=(BNbN)8.

Automatically,

bzz=bss≥0.

However, a globally convex function cannot be exactly constant on one open interval and then join a decreasing convex s−1/8 profile with the required derivative orientation without checking monotonicity carefully. So this reversal is useful conceptually, but it does not by itself solve the matching problem.

We should instead preserve convexity through the transition explicitly.

22.7. Hermite matching of bs

Let the inner regular region end at

s=s0(t),

and the exact Lambert region begin at

s=s1(t)>s0(t).

We want

bs(s0)=0

and at the outer endpoint

bs(s1)=−B8s1−9/8<0.

But if

bss≥0,

then bs must be nondecreasing as s increases.

It cannot move from

0

to a negative value while remaining nondecreasing.

Therefore:

(22.15)a spatially constant inner b cannot be joined to Bs−1/8 while maintaining bss≥0 everywhere.

This is an important obstruction.

Our previous proposal of a completely flat inner b is therefore incompatible with preserving the favorable pointwise convexity through the entire transition.

22.8. We do not actually require pointwise positive viscosity everywhere

Fortunately, the lower-barrier argument requires something weaker.

At a first downward contact,

w=b−b―

satisfies

w=0,∇w=0,Δw≥0.

For the older transformed operator,

V[b,C]−V[b―,C]=νΔw+νΔCCw−2ν∇C⋅∇w.

Thus at contact,

(22.16)V[b,C]−V[b―,C]=νΔw≥0.

Therefore the transition does not need

Δb≥0

pointwise everywhere.

What matters for the lower-barrier proof is:

  1. the operator remains smooth and parabolic;

  2. the same C is used in the comparison;

  3. the barrier itself has

    V[b―,C]>0

    in the region where the barrier argument is applied.

This means we can keep the smoothing layer outside the actual comparison tube.

That is the key geometric simplification.

22.9. Separate the regularization layer from the blowup comparison tube

Define two moving scales:

rreg(t)<rcmp(t).

Use

0≤ρ≲rreg

for regularization and interpolation.

Use

ρ≥rcmp(t)

for the exact Lambert/barrier region.

Require a nonempty exact-overlap collar

rreg(t)<ρ<rcmp(t)

where

(22.17)b=Bs−1/8,A2=−2χb

exactly.

Then:

  • the smoothing layer need not preserve the Lambert Riccati mechanism;

  • the barrier proof is performed only where the original equations and favorable asymptotics are exact;

  • smoothness and Leray estimates are handled in the inner region separately.

This is mathematically cleaner than demanding one interpolation satisfy every property simultaneously.

22.10. Choose the two scales

Take

(22.18)rreg(t)=κ0sf1/8,rcmp(t)=κ1sf1/8,0<κ0<κ1.

Then both scales have the same favorable Leray exponent, but there is a fixed relative-width exact overlap.

For the exact outer region,

‖∂zux(t)‖22≲rcmp−9/4≍sf−9/32.

Hence

(22.19)∫T∗‖∂zux(t)‖22dt<∞.

The comparison argument can begin at

ρ=rcmp(t)

without entering the smoothing layer.

22.11. Time derivative of the regularized phase

We still need the evolution equation to remain meaningful in the inner layer.

For

s^=M(s,Σ(ρ,t)),
st=1,

so

(22.20)s^t=m+(1−m)Σt+width-variation term.

If the smoothing width is proportional to Σ, the last term is controlled by

O(|Σt|).

Thus

|s^t|≲1+|Σt|.

For a floor tied to the moving front, we should avoid choosing

Σ=sf

if its derivative creates an uncontrollable scalar residual. Instead, because the regularized region is no longer used to prove the Lambert differential inequality, we only require its resulting forcing/residual to belong to the function class allowed by the full NS formulation.

This is a weaker and more appropriate requirement.

22.12. What is now established structurally

The audit has changed the architecture in an important way:

(22.21)smooth regular inner field⇓moving interpolation layer⇓exact overlap collar⇓exact Lambert/barrier region.

The exact Lambert region begins at

ρ≍sf1/8,

which passes the outer Lt2Hx1 scaling test.

The smoothing layer is no longer required to preserve

bzz>0

or the exact Lambert source. It only has to remain smooth and contribute finite dissipation.

The first-contact viscosity ordering remains exact in the comparison region:

(22.22)V[b,C]≥V[b―,C].

And the transformed-viscosity asymptotic there remains

(22.23)V[b―,C]=νBλ(λ+1)s−λ−2(1+o(1))>0.

22.13. The remaining quantitative problem is now the inner interpolation’s H1 cost

We need to construct one explicit inner interpolation and compute

∫0T∗∫IN(t)∪TN(t)|∇uN|2dxdt.

If that is finite while the exact outer Lambert/barrier region remains unchanged, then the moving-front modification passes the full local Leray dissipation audit.

23. Inner Interpolation and the Lt2Hx1 Audit

Transition. The unified phase removes the need to solve the moving-cutoff problem repeatedly, but it does not by itself establish the Leray norm. We therefore audit the inner and collar derivatives component by component. This section identifies the derivative orders that are integrable and the reconstruction of C=uz that is still needed. The following section then fixes a concrete cutoff and performs the same calculation with exact volume scaling and torus compatibility visible.

We can now perform the inner interpolation Lt2Hx1 audit explicitly. There is a particularly useful choice: interpolate the velocity variables b and A across a collar whose spatial width is comparable to the moving radius r(t), and determine exactly how closely the inner values must track the outer boundary values.

Let

r(t)=κsf(t)1/8,sf(t)≍(T∗−t)8/9,

so

(23.1)r(t)≍(T∗−t)1/9.

Take constants

0<θ0<θ1<1,

and define the inner transition collar

(23.2)T(t)={θ0r(t)<ρ<θ1r(t),c1ρ2<s<c2ρ2,y∈T}.

For ρ≥θ1r(t), retain the exact Lambert/square-root reconstruction.

23.1. Outer values at the collar

At

ρ≍r,s≍r2,

the k=3 profile gives

(23.3)bout≍s−1/8≍r−1/4,

not r−1. This corrects an overestimate in the previous exploratory calculation.

Similarly,

(23.4)Aout≍ρ1/2s−1/16≍r1/2r−1/8=r3/8,

and

(23.5)ux,out=bA≍r−1/4−3/8=r−5/8.

These are exactly the earlier parabolic-core exponents.

23.2. Why a bounded positive Ain is actually problematic

Suppose we choose

Ain=a0>0.

At the outer side of the collar,

Aout≍r3/8→0.

Thus the amplitude jump is

ΔA=a0−O(r3/8)=O(1).

Across a collar of width O(r),

Aρ=O(r−1).

The spatial volume of the parabolic collar is

(23.6)|T(t)|≍r⋅r2=r3

up to the fixed y-length.

Therefore

∫T(t)|Aρ|2dV≲r−2r3=r.

Hence

(23.7)‖∂ρA‖L2(T(t))2=O(r),

which is harmless.

So a bounded inner A is not ruled out by H1 alone.

The more restrictive quantity is again

ux=bA.

23.3. Construct an explicit collar

Choose

H∈C∞(R),

with

H(η)=0(η≤θ0),H(η)=1(η≥θ1),

and define

η=ρr(t).

Set

(23.8)b=(1−H)bin+Hbout,

and

(23.9)A=(1−H)Ain+HAout.

Take

(23.10)bin(t)≍r(t)−1/4,

so that its magnitude matches the outer value.

For example,

(23.11)bin(t)=βr(t)−1/4,β>0.

Then throughout the collar

b=O(r−1/4).

23.4. bρ

Since

Hρ=O(r−1),

we obtain

(23.12)bρ=Hρ(bout−bin)+H(bout)ρ.

If bout=b(s) and s is independent of x in the local coordinates, then

(bout)ρ=0.

Without fine matching,

bout−bin=O(r−1/4),

so

(23.13)bρ=O(r−5/4).

Therefore

(23.14)∫T|bρ|2dV≲r−5/2r3=O(r1/2).

This tends to zero.

Its time integral is certainly finite because

r(t)≍(T∗−t)1/9.

Thus the b-interpolation itself passes the spatial-gradient test.

23.5. Aρ

We have

Aout=O(r3/8),Ain=O(1),

so

Aρ=O(r−1).

As above,

(23.15)∫T|Aρ|2dV=O(r).

Again harmless.

23.6. The decisive quantity ux,ρ

Now

ux=bA,

so

(23.16)ux,ρ=bρA−bAρA2.

Because Aout→0, the worst part occurs near the outer edge of the collar.

There,

A≍r3/8,b≍r−1/4.

If

Aρ=O(r−1),

then

(23.17)bAρA2=O(r−1/4r−1r−3/4)=O(r−2).

Consequently

|ux,ρ|2=O(r−4).

Multiplying by the collar volume r3,

(23.18)∫T|ux,ρ|2dV=O(r−1).

Now

r(t)≍(T∗−t)1/9,

so

r(t)−1≍(T∗−t)−1/9.

Therefore

(23.19)∫T∗r(t)−1dt<∞.

So even this apparently severe transition derivative remains spacetime integrable.

This is an important positive result.

23.7. ux,z in the collar

Here the cutoff H(ρ/r) is independent of z at fixed t. Hence

Hz=0.

The z-dependence enters through the outer s-profile:

bout=Bs−1/8,Aout≍ρ1/2s−1/16.

At

s≍r2,

we have

(23.20)(bout)z=O(s−9/8)=O(r−9/4),

and

(23.21)(Aout)z=O(ρ1/2s−17/16)=O(r1/2r−17/8)=O(r−13/8).

The outer ux derivative has scale

(23.22)(ux,out)z≍ρ−1/2s−17/16≍r−1/2r−17/8=r−21/8.

Squaring and multiplying by r3,

(23.23)r−21/4r3=r−9/4.

Thus

(23.24)‖ux,z(t)‖L2(T)2≲Cr(t)−9/4.

This is exactly the same critical outer scaling we already obtained.

Using

r(t)≍(T∗−t)1/9,

we get

r−9/4≍(T∗−t)−1/4.

Therefore

(23.25)∫T∗‖ux,z(t)‖L2(T)2dt<∞.

23.8. uy=A

The largest transverse interpolation derivative gave

‖Aρ(t)‖22=O(r).

For the longitudinal derivative,

Az=O(r−13/8),

so

(23.26)∫T|Az|2dV≲r−13/4r3=O(r−1/4).

Since

r−1/4≍(T∗−t)−1/36,

this is easily time-integrable.

23.9. Reconstruct C=uz

Exact continuity gives

Cz=−∂x(bA)−Ay.

For the local y-independent construction,

Ay=0,

so

(23.27)Cz=−ux,x.

Thus the estimate already obtained for the transverse derivative of ux gives

(23.28)‖Cz(t)‖22≲Cr−1

for the crude interpolation.

Therefore

(23.29)∫T∗‖Cz(t)‖22dt≲∫T∗r(t)−1dt<∞.

So exact incompressibility does not, at the first-derivative level, create a worse Lt2Lx2 singularity.

23.10. But Cx requires one more derivative

This is where we must be careful.

From

Cz=−ux,x,

we can control Cz using one derivative of ux.

But obtaining Cx by differentiating the integral formula

C(x,z,t)=C0(x,t)−∫z0zux,x(x,ζ,t)dζ

gives

(23.30)Cx=C0,x−∫z0zux,xxdζ.

Thus a direct estimate of Cx from this representation involves

ux,xx,

a second derivative.

We cannot infer

Cx∈Lt2Lx2

from the first-derivative estimates alone.

This is the first remaining gap in the full nine-derivative audit.

23.11. Use a periodic divergence solver instead

There is a cleaner global approach.

Let

ux=bA,uy=A.

Define

g=−∂xux−∂yuy.

We need

∂zC=g.

On the torus this has a periodic solution iff the zero z-Fourier mode vanishes:

(23.31)∫Tzg(x,y,z,t)dz=0.

When this compatibility holds, define

(23.32)C=∂z−1g+C0(x,y,t),

where ∂z−1 acts on nonzero z-Fourier modes.

Then

Cz=g

exactly.

However,

Cx=∂z−1gx+C0,x,

still requires control of gx, hence generally a second horizontal derivative of ux.

So the inverse derivative in z does not magically give full H1 regularity from ux,A∈H1.

This is a genuine functional-analytic issue.

23.12. A better way: reconstruct the divergence-free field by a potential

To avoid losing a derivative, define the x,z components through a stream function ψN:

(23.33)ux=∂zψN,uz=−∂xψN,

and take

uy=AN

with

AN,y=0.

Then automatically

∂xux+∂zuz=0,

and hence

(23.34)∇⋅u=0.

Now impose the product relation

b=uxuy

by requiring

(23.35)ψz=bA.

Thus

(23.36)ψ(x,z,t)=ψ0(x,t)+∫z0zb(x,ζ,t)A(x,ζ,t)dζ,

and

uz=−ψx.

This is equivalent locally to the continuity reconstruction, but it makes the derivative requirement transparent:

∇uz

involves second derivatives of ψ.

Hence full

u∈H1

requires roughly

ψ∈H2.

So we must audit the second derivatives of b/A sufficiently to establish this.

23.13. The next exponent: ux,xx

The crude transition estimate was

ux,x=O(r−2).

Differentiating once more across a collar of width r gives conservatively

(23.37)ux,xx=O(r−3).

Its squared collar integral scales like

(23.38)r−6r3=r−3.

Since

r(t)≍(T∗−t)1/9,
r−3≍(T∗−t)−1/3,

which is still time-integrable:

(23.39)∫T∗r(t)−3dt<∞.

This is encouraging.

A third transverse derivative would give

r−4

pointwise and roughly

r−5

after squaring/integrating over the collar, corresponding to

(T∗−t)−5/9,

which is also integrable.

The basic moving scale r∼(T∗−t)1/9 therefore gives substantial derivative room.

23.14. General collar derivative criterion

Suppose a quantity of amplitude

F(t)≍r−a

is interpolated over width r.

Its m-th transverse derivative scales like

∂ρmF=O(r−a−m).

The collar volume is O(r3), so

(23.40)‖∂ρmF‖22=O(r3−2a−2m).

Since

r(t)≍(T∗−t)1/9,

the time integral converges provided

2a+2m−39<1.

Equivalently,

(23.41)a+m<6.

This is a useful general criterion for the redesigned collar.

For

ux≍r−5/8,

we have

a=58.

Thus transverse derivatives up through

m=5

still satisfy

58+5<6.

So the moving r(t)∼(T∗−t)1/9 collar has ample formal spacetime integrability for the finite number of derivatives needed to reconstruct C.

23.15. Result of the local Leray scaling audit

Subject to smooth interpolation profiles with the assumed collar-width scaling, we now obtain:

∫0T∗∫I(t)∪T(t)|∇ux|2dVdt<∞,
∫0T∗∫I(t)∪T(t)|∇uy|2dVdt<∞,

and the derivative-loss needed to reconstruct uz is not ruled out by scaling: even the conservative second-horizontal-derivative estimate

‖ux,xx(t)‖22=O(r−3)

is time-integrable.

Together with the outer estimate

‖∂zux(t)‖22≲(T∗−t)−1/4,

the redesigned moving front therefore passes the local power-counting test for finite Leray dissipation.

What we have not yet proved is the corresponding exact norm estimate for uz, because that requires choosing an explicit interpolation H, performing the full second-derivative calculation rather than the conservative O(r−3) bound, and imposing the periodic z-compatibility condition (23.31).

That is now narrow enough to finish directly: choose a concrete C∞ flat cutoff H, compute ux,x and ux,xx exactly in scaled coordinates η=ρ/r(t), reconstruct C, and prove its three first derivatives are in Lt2Lx2. If that succeeds, the local Lt2Hx1 reconstruction gap is closed for the redesigned moving collar.

24. Explicit Collar Calculation, Exact Volume Scaling, and Periodic Compatibility

**Transition.**The preceding scaling audit is now converted into an explicit collar construction. We interpolate ux directly, compute the relevant derivatives with the exact collar volume, and reconstruct the remaining component subject to incompressibility. This closes the local singular-exponent part of the Leray audit. The final subsection of this calculation also isolates the one issue that is not closed locally: a periodic divergence correction must preserve the exact outer Lambert comparison field.

We can now make the collar calculation explicit. There is also one correction to the preceding power count: because this paper’s local volume element is dρdsdy, the collar-volume estimate O(r3) is correct only when the collar is simultaneously restricted to s≍r2. We will impose that explicitly.

24.1. Fix a concrete smooth flat cutoff

A convenient C∞ step is obtained from

ϑ(q)={0,q≤0,e−1/q,q>0,

and

S(q)=ϑ(q)ϑ(q)+ϑ(1−q).

Then

S(q)=0(q≤0),S(q)=1(q≥1),

and every derivative vanishes at q=0,1.

Choose

0<θ0<θ1<1,η=ρr(t),

and define

(24.1)H(η)=S(η−θ0θ1−θ0).

Therefore

(24.2)Hρ=1rHη,Hρρ=1r2Hηη,

with

‖Hη‖∞+‖Hηη‖∞<∞.

All derivatives vanish at the two collar boundaries. Hence the interpolated field glues smoothly to both neighboring pieces.

24.2. It is better to interpolate ux directly

Previously we wrote separate interpolations of b and A. For the H1 audit this unnecessarily complicates the quotient

ux=bA.

Instead define

U:=ux

first.

On the exact outer square-root region,

(24.3)Uout=bA=B2ρ−1/2s−1/16.

Choose an inner profile Uin(t) of the same boundary magnitude,

(24.4)Uin(t)=μr(t)−5/8,

where μ>0 is fixed.

Then put

(24.5)U(ρ,s,t)=(1−H(η))Uin(t)+H(η)Uout(ρ,s).

Afterward define

A=bU.

This guarantees

b=UA

identically instead of hoping two independent interpolations preserve the desired quotient.

We must still check positivity:

U>0,b>0

implies

A=bU>0.

24.3. Exact first transverse derivative

Because Uin is independent of ρ,

(24.6)Uρ=Hρ(Uout−Uin)+H(Uout)ρ.

From (24.3),

(24.7)(Uout)ρ=−12B2ρ−3/2s−1/16.

Inside the collar,

ρ≍r,s≍r2.

Therefore

Uout=O(r−5/8),

and

(Uout)ρ=O(r−13/8).

Also

Hρ=O(r−1),

so

Hρ(Uout−Uin)=O(r−13/8).

Consequently

(24.8)Uρ=O(r−13/8).

This improves the deliberately crude O(r−2) estimate used earlier.

Squaring,

|Uρ|2=O(r−13/4).

The collar is

ρ≍r,s≍r2,

so

|T(t)|≍r3.

Hence

(24.9)‖Uρ(t)‖L2(T)2=O(r−1/4).

With

r(t)≍(T∗−t)1/9,

this becomes

O((T∗−t)−1/36),

which is integrable.

24.4. Exact second transverse derivative

Differentiate (24.6):

(24.10)Uρρ=Hρρ(Uout−Uin)+2Hρ(Uout)ρ+H(Uout)ρρ.

Now

(24.11)(Uout)ρρ=34B2ρ−5/2s−1/16.

Thus on the collar,

(Uout)ρρ=O(r−21/8).

The other terms satisfy

Hρρ(Uout−Uin)=O(r−2r−5/8)=O(r−21/8),

and

Hρ(Uout)ρ=O(r−1r−13/8)=O(r−21/8).

Therefore

(24.12)Uρρ=O(r−21/8).

Consequently

(24.13)‖Uρρ(t)‖22=O(r−21/4r3)=O(r−9/4).

Since

r−9/4≍(T∗−t)−1/4,

we obtain

(24.14)∫T∗‖Uρρ(t)‖22dt<∞.

This is precisely the second-horizontal-derivative estimate needed for the simplest continuity reconstruction of C.

24.5. Longitudinal derivative

Since H(ρ/r(t)) has no z-dependence at fixed t,

Uz=H(Uout)z

if the inner field is z-independent.

Using

∂z=−∂s,

we have

(24.15)(Uout)z=116B2ρ−1/2s−17/16.

Thus

(24.16)(Uout)z=O(r−1/2r−17/8)=O(r−21/8).

Therefore

(24.17)‖Uz(t)‖L2(T)2=O(r−9/4).

Again,

∫T∗r−9/4dt<∞.

So both

UρρandUz

have the same favorable spacetime exponent.

24.6. Mixed derivative Uρz

We will need this for Cz and Cρ.

From

Uz=H(Uout)z,
Uρz=Hρ(Uout)z+H(Uout)ρz.

Now

(Uout)ρz=−132B2ρ−3/2s−17/16.

Hence

(Uout)ρz=O(r−3/2r−17/8)=O(r−29/8).

Likewise

Hρ(Uout)z=O(r−1r−21/8)=O(r−29/8).

Therefore

(24.18)Uρz=O(r−29/8).

Its squared collar norm is

(24.19)O(r−29/4r3)=O(r−17/4).

Since

r−17/4≍(T∗−t)−17/36,

and

1736<1,

we still have

(24.20)Uρz∈Lt2Lx2

over the moving collar.

This is an important strengthening.

24.7. Reconstruct C explicitly

In the local chart, let

ρ=ρ(x)

with

0<cρ≤|ρx|≤Cρ

on the chosen signed side of the simple zero.

Assume local y-independence:

Ay=0.

Exact continuity is

Cz=−Ux.

Since

Ux=Uρρx,

define

(24.21)C(x,z,t)=C0(x,t)−∫z0zUρ(x,ζ,t)ρx(x)dζ.

Then

Cz=−Ux

exactly and hence

(24.22)∇⋅u=0.

24.8. Estimate Cz

Immediately,

|Cz|≲|Uρ|.

Thus

(24.23)‖Cz(t)‖22=O(r−1/4)

for the pure collar contribution, and therefore

Cz∈Lt2Lx2.

24.9. Estimate Cx

Differentiate (24.21):

(24.24)Cx=C0,x−∫z0z[Uρρρx2+Uρρxx]dζ.

Assume the local coordinate satisfies

|ρx|+|ρxx|≤K.

Then, on a z-interval of uniformly bounded length, Cauchy–Schwarz gives

|Cx|2≲|C0,x|2+∫(|Uρρ|2+|Uρ|2)dζ.

After integrating in x,y,

(24.25)‖Cx(t)‖22≲‖C0,x(t)‖22+C(r−9/4+r−1/4).

Hence the dominant contribution is

(24.26)‖Cx(t)‖22=O(r−9/4),

provided the integration datum is uniformly controlled.

Therefore

(24.27)Cx∈Lt2Lx2.

24.10. Cy

Under the local y-independent construction and a y-independent integration datum,

(24.28)Cy=0.

Thus all three first derivatives of C pass the local scaling audit:

(24.29)∇C∈Lt2Lx2.

24.11. What about A=b/U?

We now have to make sure our convenient interpolation of U=ux does not merely move the problem into A.

Take

b=O(r−1/4)

through the collar and

U=O(r−5/8).

Then

(24.30)A=bU=O(r3/8),

which matches the square-root outer scale.

This is better than interpolating to an O(1) inner A: the entire collar now respects the natural vanishing scale.

If

bρ=O(r−5/4),Uρ=O(r−13/8),

then

Aρ=bρU−bUρU2.

The two terms have scale

r−5/4r5/8=r−5/8,

and

r−1/4r−13/8r10/8=r−5/8.

Hence

(24.31)Aρ=O(r−5/8).

Its squared collar norm is

(24.32)r−5/4r3=r7/4,

which actually vanishes.

Similarly, the outer z-scaling gives

Az=O(r−13/8),

so

(24.33)‖Az(t)‖22=O(r−1/4),

again integrable in time.

Thus

(24.34)∇A∈Lt2Lx2

on the collar.

24.12. The resulting local H1 estimate

For

u=(U,A,C),

the estimates above yield schematically

‖∇U(t)‖22≲Cr−9/4,
‖∇A(t)‖22≲Cr−1/4,

and

‖∇C(t)‖22≲C(1+r−9/4).

Therefore

(24.35)‖∇u(t)‖L2(T(t))2≲C(1+r(t)−9/4).

With

r(t)≍(T∗−t)1/9,

we obtain

(24.36)‖∇u(t)‖22≲C[1+(T∗−t)−1/4].

Consequently

(24.37)∫t0T∗‖∇u(t)‖L2(T(t))2dt<∞.

This closes the power-counting part of the moving-collar local Leray dissipation problem.

24.13. But there are two compatibility conditions we cannot skip

First, equation (24.21) defines a periodic C only if

(24.38)∫TzUxdz=0.

Equivalently,

(24.39)∂x∫TzUdz=0.

Thus the z-mean of U must be independent of x.

A purely local collar interpolation does not automatically satisfy this.

Second, if C is reconstructed from continuity, it is not free to simultaneously satisfy an independently prescribed C-formula from the older construction. We must use the reconstructed C consistently in

V[b,C]=ν[Δb+bCΔC−2∇C⋅∇b]

and in the full momentum equations.

24.14. The clean way to enforce the torus condition

Let Uraw denote the interpolated field (24.5), and define its z-mean

U―raw(x,t)=1|T|∫TzUraw(x,z,t)dz.

Choose an x-independent target mean mN(t), and set

(24.40)UN=Uraw−U―raw(x,t)+mN(t).

Then

∫TzUNdz=|T|mN(t),

so

(24.41)∂x∫TzUNdz=0.

Consequently a periodic CN solving

CN,z=−UN,x

exists.

Moreover,

UN,x=Uraw,x−∂xU―raw.

The mean-subtraction operator is bounded on Sobolev spaces, so it does not worsen the derivative order or the power-counting exponents derived above.

This gives a genuine route to exact periodic incompressibility.

But it introduces the next nontrivial question: the correction

−U―raw+mN(t)

is global in z. It may alter the exact Lambert relation

Uout=bA

outside the moving collar.

The next decisive step is to construct the periodic compatibility correction so that it is supported in the regular inner region—or equivalently use a divergence correction/Bogovskiı̆-type construction localized away from the exact Lambert comparison region. If such a localized correction has the same Lt2Hx1 bounds, then we can preserve the exact outer Lambert field and obtain global periodic incompressibility simultaneously.

That localized divergence correction is now the precise remaining problem; importantly, the singular-exponent obstruction itself has passed the moving-collar audit.

25. Consolidated Leray conclusion of the moving-coterminal audit

**Logical conclusion of the sequence.**The preceding sections should now be read as one chain rather than as competing constructions. The fixed wedge establishes the derivative obstruction; the finite-N release audit shows that the original release interpolation does not remove it; the moving-front criterion identifies an integrable scale; the interface and high-order audits eliminate naive cutoff choices; the one-front/regularized-phase construction removes those cutoff mismatches; and the explicit collar calculation verifies the resulting local Leray power count. What remains is not another local exponent problem but the global periodic compatibility correction described below.

The developmental calculations above establish the following conditional but precise conclusion for the terminal reconstruction.

For the nonmoving coterminal region

ΩN,ε={0<ρ<ε, c1ρ2≤s≤c2ρ2, y∈T},

the velocity itself can remain locally square integrable, while the fixed-time H1 contribution of the most singular derivative diverges as the lower endpoint tends to zero. Thus bounded kinetic energy alone does not close the Leray audit.

For the redesigned moving collar with

r(t)≍(T∗−t)1/9,θ0r(t)<ρ<θ1r(t),c1ρ2<s<c2ρ2,

the explicit interpolation calculation yields the controlling estimate

‖∇u(t)‖L2(T(t))2≲C(1+r(t)−9/4)≲C(1+(T∗−t)−1/4).

Consequently,

∫T∗−δT∗‖∇u(t)‖L2(T(t))2dt<∞,

because 1/4<1. This is the precise sense in which the moving coterminal construction repairs the terminal spacetime power-counting obstruction that occurs in the nonmoving wedge.

The same calculation leaves a global compatibility obligation. Exact incompressibility may be reconstructed locally through Cz=−Ux (in the y-independent collar), while periodicity requires the corresponding torus mean condition. The mean-subtraction construction preserves the Sobolev power count but is global in z and may alter the exact outer Lambert relation. Therefore the remaining closure problem is to construct a localized periodic divergence correction, supported away from the exact Lambert comparison region, with the same Lt2Hx1 bounds. Thus the local terminal singular-exponent audit is successful for the explicit moving collar, while the global Leray closure remains conditional on the localized periodic divergence correction and on preserving the exact outer comparison field.

26. Distinguished phase θ=t−z and continuity-corrected cascade

26.1. Revised specialization rule

In the revised construction the phase specialization θ=t−z is retained, but the old conclusion uy=uz=ΔN,k is only used as stepping stone in it’s analyses leading to the parabolic construction and towards the terminal branch point. The Lambert cascade supplies the phase coefficient and the b-equation; the velocity components are reconstructed from A2=−2χb, ux=b/A, and exact continuity for C=uz.

27. Specialization to θ=t−z in the Revised Geometry

The Lambert branch variable is still specialized to

(27.1)ΔN,k=ΔN,k(θ),θ=t−z.

Therefore

(27.2)ΔN,k,x=0,ΔN,k,y=0,ΔN,k,z=−ΔN,k,θ.

What changes is the velocity reconstruction. The revised field is not obtained by identifying either transverse component with ΔN,k. Instead set

(27.3)b=b3=uxuy,A=uy,C=uz,ux=bA.

On the periodic terminal branch impose

(27.4)A2=−2χ(x)b.

Equivalently,

2χ(x)ux+uy=0.

The exact incompressibility equation is

(27.5)∂x(bA)+Ay+Cz=0.

Let

s=θ−θ∗,N=t−z−θ∗,N

and suppose in the local terminal core that

b=b(s).

Then

bx=by=0,A=A(x,s),Ay=0.

Since

∂z=−∂s,

equation (27.5) gives

Cs=∂x(bA).

From

A2=−2χb

and bx=0,

2AAx=−2χ′b,

so

Ax=−χ′bA.

Also

bA=−A2χ.

Hence

∂x(bA)=−12∂x(Aχ)=−Ax2χ+Aχ′2χ2=χ′b2Aχ+Aχ′2χ2.

Because

b=−A22χ,

the first term is

χ′b2Aχ=−Aχ′4χ2.

Therefore

(27.6)Cs=Aχ′4χ2.

This directly verifies incompressibility:

ux,x+uy,y+uz,z=∂x(b/A)+0−Cs=0.

Thus

(27.7)∇⋅u=0

exactly in the revised distinguished-phase reconstruction.

If

b=BNs−αFN(s),

then

A=−2BNχ(x)s−α/2FN(s)1/2.

Define

IN(s)=∫0sr−α/2FN(r)1/2dr.

Integrating (27.6) gives

(27.8)C(x,s)=C0(x)+χ′(x)−2BNχ(x)4χ(x)2IN(s).

If C0=c0>0 is constant on the local core chart, then

Ct+Cz=0

exactly, because all t,z-dependence occurs through s=t−z−θ∗,N. Thus

(27.9)Ct+Cz=0.

28. The Logarithmic Derivative in the Modified Construction

The logarithmic derivative is

(28.1)∂zlog⁡ΔN,k=−ΔN,k,θΔN,k.

Near the branch endpoint,

(28.2)ΔN,k∼CN,ksαk

and

(28.3)ΔN,k,θ∼CN,kαksαk−1.

Consequently,

(28.4)∂zlog⁡ΔN,k∼−αks.

This is precisely the singular coefficient entering the reduced characteristic amplitude equation.

Thus the same branch derivative controls both

(28.5)b˙3=−b32∂zlog⁡ΔN,k

and the surviving continuity-corrected velocity term

(28.6)ux−HN,k=(x−x0)ΔN,k,θ.

29. Common-Time Behavior of the Continuity-Corrected Velocity

We now evaluate the same coupled sequence used for the amplitude:

(29.1)tN=T∗,N−N−N.

Write

(29.2)τN=N−N.

The local time law gives

(29.3)τN∼DN,ksN1+αN,

where

(29.4)αN=2−k(N).

Therefore

(29.5)sN∼(τNDN,k)1/(1+αN).

The derivative satisfies

(29.6)ΔN,k,θ∼CN,kαNsNαN−1.

Taking logarithms,

(29.7)log⁡|ΔN,k,θ|=log⁡CN,k+log⁡αN+αN−11+αN(log⁡τN−log⁡DN,k)+o(1).

Since

(29.8)log⁡τN=−Nlog⁡N

and

(29.9)αN→0,

the dominant term is

(29.10)Nlog⁡N.

Under the bounded-prefactor assumptions verified numerically in the reduced characteristic model,

(29.11)log⁡|ΔN,k(N),θ(tN)|=(1+o(1))Nlog⁡N.

Equivalently,

(29.12)|ΔN,k(N),θ(tN)|=NN(1+o(1)).

Thus

(29.13)|ΔN,k(N),θ(tN)|→+∞.

The numerical ratios

(29.14)log⁡|ΔN,k,θ(tN)|Nlog⁡N

were found to be

(29.15)0.97547, 0.98761, 0.99377, 0.99687, 0.99843, 0.99921, 0.99961, 0.99980, 0.99990,

for increasing values of N, strongly confirming

(29.16)log⁡|ΔN,k,θ(tN)|Nlog⁡N→1.

30. Consequence for the Revised Physical Velocity Components

The old continuity formula

ux=HN,k+(x−x0)ΔN,k,θ

belongs to the historical reconstruction in which the transverse Lambert profile itself was inserted directly as a velocity component. It is not the active formula in the square-root reconstruction.

The active identities are

uy=A,ux=bA,A2=−2χb.

Hence

(30.1)ux2=−b2χ

on the branch χ<0, b>0. Equivalently,

(30.2)|uxuy|=12|χ|.

In the terminal tube

c1s≤χ2≤c2s,

we have

|χ|≍s1/2.

If

b≍s−α,

then

uy2=2|χ|b≍s1/2−α,

and therefore

(30.3)uy≍s1/4−α/2.

Similarly,

ux2=b2|χ|≍s−α−1/2,

so

(30.4)ux≍s−1/4−α/2.

For k=3,

α=18,

and hence

(30.5)uy≍s3/16,ux≍s−5/16,ux2≍s−5/8.

The continuity reconstruction gives

(30.6)uz=C=c0+O(s3/16),c0>0.

Thus the revised physical velocity amplification is carried directly by the product amplitude b together with the shrinking geometric factor |χ|−1/2. It is not transferred from the derivative ΔN,k,θ through an integration function HN,k.

31. Behavior of the Two Reconstructed Transverse Components

The two transverse components have different roles:

uy=A=−2χb,

while

uz=C

is determined by exact continuity.

At k=3,

uy→0

as s→0+, with

uy≍s3/16,

whereas

uz→c0>0.

Thus the revised construction does not have a common transverse terminal value.

The Lambert cascade still controls the singular coefficient

ΓN=∂θlog⁡QN

appearing in the b-equation. It should not be reinterpreted as saying that QN itself equals uy or uz.

This separation is essential:

QN is the Lambert branch variable,A,C are reconstructed velocity components.

32. Summary Theorem for the Revised Reduced Construction

Conditional theorem. Let the renormalized Lambert cascade define QN and

ΓN=∂θlog⁡QN.

Assume that the preterminal periodic reconstruction and the compact Floquet/Fredholm phase homotopy can be continued to the distinguished phase

s=t−z−θ∗,N,

with entry data

b(se)≥1+δ,C(se)≥c−>0.

Assume further that the tuning ratio is released before the square-root geometry is imposed, and that the geometric release connects smoothly to a terminal core satisfying

A2=−2χb,c1s≤χ2≤c2s.

Then:

  1. b=uxuy

    is preserved exactly.

  2. The local terminal reconstruction satisfies

    ∇⋅u=0

    exactly.

  3. The release layer satisfies

    |Rrel|ΓNb2=O(s7/4)

    for k=3.

  4. The square-root core satisfies

    |Rcore|ΓNb2=O(s).
  5. If

    G≥0,ΓN(s)=18s(1+o(1)),

    then for sufficiently small s,

    b˙≥(1−ϑ)ΓNb2,s˙=1−b.
  6. Consequently,

    b(s)≥be(ses)p

    for every fixed p<1/8 sufficiently close to 1/8, after shrinking the terminal interval if necessary.

  7. At the formal k=3 square-root scaling,

    b3≍s−1/8,uy≍s3/16,ux≍s−5/16,uz=c0+O(s3/16).

This is a theorem for the revised reduced reconstruction conditional on the stated continuation hypotheses. It is not yet a theorem for a full three-dimensional Navier–Stokes solution.

33. Analysis of the Renormalized Lambert Branch-Point Cascade and Its Characteristic Blow-Up Mechanism

We give a rigorous analysis of the renormalized Lambert recursion

(33.1)Δk+1=1+W0(−e−1−Δk),Δ0≥0.

The purpose of the renormalization is to map the Lambert branch point back to itself. Unlike the literal composition W0∘W0, which leaves the real branch domain after the first iteration near −1/e, the renormalized map preserves the nonnegative real half-line.

We prove that the map is real-valued and invariant on [0,∞), that

(33.2)Δk+1=2Δk−23Δk+O(Δk3/2),

(see Appendix A for the proof)and consequently that, for every fixed finite k,

(33.3)Δk=21−2−kΔ02−k(1+O(Δ02−k)).

(See Appendix B for the proof) The exponents are therefore

(33.4)1,12,14,18,…,

and the leading constants satisfy

(33.5)Ck+1=2Ck,C0=1,

giving

(33.6)Ck=21−2−k.

We then derive the corresponding spatial derivative and logarithmic derivative asymptotics. Finally, for the reduced characteristic equation

(33.7)b˙3=b32∂θlog⁡Δk,θ˙=1−b3,

we derive the exact characteristic invariant

(33.8)1b3−log⁡1b3+log⁡Δk=C.

This invariant implies

(33.9)Δk=Kb3−1e−1/b3,

and therefore

(33.10)b3∼KΔk.

For every finite k, the characteristic reaches the singular layer in finite time and

(33.11)b3(t)≍(T∗−t)−1/(2k+1).

We also explain carefully why the limit k→∞ cannot simply be taken term-by-term and why a coupled scaling such as k∼log2⁡n is required if the cascade depth and the small initial branch distance are allowed to vary simultaneously.

34. Statement of the Renormalized Problem

The central object is the map

(34.1)F(Δ)=1+W0(−e−1−Δ).

The renormalized recursion is

(34.2)Δk+1=F(Δk),Δ0≥0.

The subscript n can be retained when the initial quantity depends on the geometric parameter n. For the local analysis it is convenient to suppress n temporarily and write

(34.3)Δ0=Δn,0,Δk=Δn,k.

The essential claim to be proved is

(34.4)Δk=21−2−kΔ02−k(1+O(Δ02−k))

for every fixed finite k as

(34.5)Δ0↓0.

The proof consists of four logically separate parts:

  1. (i) prove that F is real and maps [0,∞) into [0,∞);

  2. (ii) derive the precise local Puiseux expansion of F;

  3. (iii) iterate the expansion with controlled error;

  4. (iv) transfer the resulting spatial singularity to the characteristic equation.

35. Why the Renormalization Is Necessary

Before proving the new recursion, we explain why it differs fundamentally from the ordinary Lambert composition.

Consider

(35.1)X(Δ)=−1ee−Δ.

At

(35.2)Δ=0,

we have

(35.3)X(0)=−1e.

The principal Lambert function satisfies

(35.4)W0(−1/e)=−1.

Thus the ordinary first iterate produces

(35.5)W0(X(Δ))⟶−1.

But the branch point of the next Lambert function is

(35.6)−1e,

not −1.

Since

(35.7)−1<−1e,

the second ordinary real Lambert composition is not defined near Δ=0.

The renormalized transformation instead defines

(35.8)Δ1=1+W0(X(Δ0)).

Consequently,

(35.9)Δ0=0⟹Δ1=0.

Hence the branch point is mapped back to itself:

(35.10)0⟼0.

This is the structural reason the square-root singularity can now be iterated.

36. Exact Parametrization of the Renormalized Map

Let

(36.1)q=F(Δ).

Then

(36.2)q=1+W0(−e−1−Δ).

Therefore

(36.3)W0(−e−1−Δ)=q−1.

By the defining relation

(36.4)W0(x)eW0(x)=x,

we obtain

(36.5)(q−1)eq−1=−e−1−Δ.

Multiplying by e gives

(36.6)(q−1)eq=−e−Δ.

Equivalently,

(36.7)(1−q)eq=e−Δ.

Taking logarithms,

(36.8)log⁡(1−q)+q=−Δ.

Hence

(36.9)Δ=−log⁡(1−q)−q.

Define

(36.10)H(q)=−log⁡(1−q)−q.

Then the renormalized map is equivalently characterized by

(36.11)H(F(Δ))=Δ.

This identity is exact.

It is also extremely useful because it removes Lambert W from the local inversion problem.

37. Positivity and the Real Invariant Domain

We now prove that the recursion is real for every finite number of iterations.

Lemma 32.1. For every Δ≥0,

(37.1)−e−1−Δ∈[−1/e,0).

Consequently,

(37.2)W0(−e−1−Δ)∈[−1,0),

and therefore

(37.3)F(Δ)∈[0,1).

Proof. If Δ≥0, then

(37.4)−1−Δ≤−1.

Hence

(37.5)0<e−1−Δ≤e−1.

Multiplication by −1 reverses the inequality:

(37.6)−1e≤−e−1−Δ<0.

The principal real Lambert branch W0 is real and increasing on [−1/e,∞). Therefore

(37.7)W0(−e−1−Δ)∈[−1,0).

Adding 1 gives

(37.8)0≤1+W0(−e−1−Δ)<1.

Thus

(37.9)F(Δ)∈[0,1).

◻

Corollary 32.1. If

(37.10)Δ0≥0,

then

(37.11)Δk∈[0,1)

for every integer k≥1.

Proof. The result follows by induction from the previous lemma.◻

Thus

(37.12)[0,∞)⟶[0,1)⟶[0,1)⟶⋯

is an invariant real recursion.

This is precisely the property that the unrenormalized Lambert composition does not possess.

38. Branch-Point Fixed Point, Cascade Derivatives, and Characteristic Invariants

38.1. The Fixed Point at the Branch Point

The point Δ=0 is a fixed point.

Indeed,

(38.1)F(0)=1+W0(−e−1)
(38.2)=1−1
(38.3)=0.

Therefore

(38.4)F(0)=0.

However, this fixed point is not an ordinary differentiable fixed point.

The derivative is singular.

This is the source of the square-root behavior.

39. Exact Equation for the Puiseux Expansion

From the previous exact identity,

(39.1)Δ=−log⁡(1−q)−q,q=F(Δ).

For |q|<1,

(39.2)−log⁡(1−q)=q+q22+q33+q44+⋯.

Therefore

(39.3)Δ=q22+q33+q44+⋯.

In particular,

(39.4)Δ=q22+O(q3).

Thus

(39.5)q2=2Δ(1+O(q)).

Since q→0 as Δ→0,

(39.6)q=2Δ(1+O(Δ)).

Consequently,

(39.7)F(Δ)=2Δ+O(Δ).

This already proves the square-root law.

40. Higher-Order Puiseux Expansion

We now derive the first several coefficients.

Write

(40.1)s=2Δ.

Seek an expansion of the form

(40.2)q=s+as2+bs3+cs4+O(s5).

Since

(40.3)Δ=s22,

we substitute into

(40.4)Δ=−log⁡(1−q)−q.

Using

(40.5)−log⁡(1−q)−q=q22+q33+q44+q55+O(q6),

and comparing powers of s, one obtains

(40.6)a=−13,
(40.7)b=136,
(40.8)c=1270.

Therefore

(40.9)F(Δ)=2Δ−23Δ+218Δ3/2+2135Δ2+O(Δ5/2).

The first term is the only one needed for the exponent cascade, but the higher-order terms are useful for error control.

41. A Two-Sided Square-Root Estimate

The asymptotic formula can be strengthened to a two-sided estimate.

Lemma 37.1. There exist constants δ∗>0 and 0<c1<c2<∞ such that for all

(41.1)0<Δ<δ∗,
(41.2)c1Δ≤F(Δ)≤c2Δ.

Proof. From

(41.3)F(Δ)=2Δ(1+O(Δ)),

there exists δ∗>0 such that

(41.4)|O(Δ)|≤12

whenever

(41.5)0<Δ<δ∗.

Therefore

(41.6)122Δ≤F(Δ)≤322Δ.

Thus we may take, for example,

(41.7)c1=12,c2=32.

◻

This estimate proves that the map has genuine square-root scaling and not merely a formal asymptotic resemblance to a square root.

42. Derivative of the Renormalized Lambert Map

Differentiate

(42.1)Δ=−log⁡(1−q)−q.

We obtain

(42.2)1=q1−qdqdΔ.

Therefore

(42.3)dqdΔ=1−qq.

Since

(42.4)q=F(Δ),

we have

(42.5)F′(Δ)=1−F(Δ)F(Δ).

As

(42.6)F(Δ)∼2Δ,

we obtain

(42.7)F′(Δ)∼12Δ.

More precisely,

(42.8)F′(Δ)=12Δ−23+O(Δ).

Thus the map is continuous at zero but has an unbounded derivative there.

43. Monotonicity

Since

(43.1)0≤F(Δ)<1

for Δ≥0, the exact derivative formula gives

(43.2)F′(Δ)>0

for Δ>0.

Hence

(43.3)F is strictly increasing on (0,∞).

This fact will be useful when transferring monotonicity from Δ0(θ) to Δk(θ).

44. The Main Puiseux Cascade Theorem

We now prove the central formula.

Theorem 40.1. Let

(44.1)Δk+1=F(Δk),Δ0>0,

where

(44.2)F(Δ)=1+W0(−e−1−Δ).

For every fixed integer k≥0,

(44.3)Δk=CkΔ02−k(1+O(Δ02−k))

as Δ0↓0, where

(44.4)C0=1,Ck+1=2Ck.

Consequently,

(44.5)Ck=21−2−k,

and therefore

(44.6)Δk=21−2−kΔ02−k(1+O(Δ02−k)).

Proof. We proceed by induction.

For k=0,

(44.7)Δ0=1⋅Δ0,

so

(44.8)C0=1.

Assume that

(44.9)Δk=CkΔ02−k(1+O(Δ02−k)).

Set

(44.10)x=Δ02−k.

Then

(44.11)Δk=Ckx(1+O(x)).

Using the local expansion

(44.12)F(y)=2y(1+O(y)),

we obtain

(44.13)Δk+1=F(Δk)
(44.14)=2Δk(1+O(Δk)).

Since

(44.15)Δk=Ckx(1+O(x)),

we have

(44.16)Δk=Ckx1/2(1+O(x)).

Therefore

(44.17)Δk+1=2Ckx1/2(1+O(x1/2)).

But

(44.18)x1/2=Δ02−(k+1).

Hence

(44.19)Δk+1=Ck+1Δ02−(k+1)(1+O(Δ02−(k+1))),

where

(44.20)Ck+1=2Ck.

This completes the induction.

It remains to solve the recurrence for Ck.

Taking logarithms,

(44.21)log2⁡Ck+1=12(1+log2⁡Ck).

Set

(44.22)ak=log2⁡Ck.

Then

(44.23)ak+1=12(1+ak),a0=0.

The solution is

(44.24)ak=1−2−k.

Therefore

(44.25)Ck=21−2−k.

The theorem follows.◻

45. Explicit First Four Iterates

The theorem gives

(45.1)Δ1∼21/2Δ01/2,

that is,

(45.2)Δ1∼2Δ0.

For the second iterate,

(45.3)C2=23/4,

and

(45.4)2−2=14.

Therefore

(45.5)Δ2∼23/4Δ01/4.

For the third iterate,

(45.6)C3=27/8,

and

(45.7)2−3=18.

Therefore

(45.8)Δ3∼27/8Δ01/8.

For the fourth iterate,

(45.9)C4=215/16,

and

(45.10)2−4=116.

Therefore

(45.11)Δ4∼215/16Δ01/16.

Thus the exponent sequence is

(45.12)1⟶12⟶14⟶18⟶116⟶⋯.

46. Derivative of the k-th Renormalized Iterate

Suppose

(46.1)Δk=CkΔ0αk(1+o(1)),αk=2−k.

Differentiating with respect to θ gives

(46.2)∂θΔk=CkαkΔ0αk−1∂θΔ0(1+o(1)).

Thus

(46.3)∂θΔk∼21−2−k2−kΔ02−k−1∂θΔ0.

Therefore

(46.4)∂θΔk

has singular exponent

(46.5)2−k−1.

Since

(46.6)2−k−1>−1

for every finite k, the singularity is locally integrable.

47. The Logarithmic Derivative

Divide the preceding formula by

(47.1)Δk∼CkΔ02−k.

The constant Ck cancels:

(47.2)∂θΔkΔk∼2−k∂θΔ0Δ0.

Hence

(47.3)∂θlog⁡Δk∼2−k∂θlog⁡Δ0.

This identity is particularly important for the characteristic equation.

The logarithmic derivative becomes more concentrated near the branch point even though the explicit prefactor 2−k decreases.

48. Integrability of the Finite-k Spatial Singularity

Suppose

(48.1)Δ0(θ)∼A(θ−θ∗)

with

(48.2)A>0.

Then

(48.3)Δk(θ)∼CkA2−k(θ−θ∗)2−k.

Consequently,

(48.4)∂θΔk∼Ck2−kA2−k(θ−θ∗)2−k−1.

The singularity is locally integrable because

(48.5)2−k−1>−1.

Indeed,

(48.6)∫0δs2−k−1ds=2kδ2−k.

Thus

(48.7)∫0δs2−k−1ds<∞

for every finite k.

However,

(48.8)2−k−1→−1.

Thus the infinite cascade approaches the borderline exponent.

49. Geometric Initial Branch Distance

We now restore the geometric variable θ and write

(49.1)Δn,0(θ)=Rn(θ)−Ln,

where

(49.2)Rn(θ)=11−sin⁡θ+εn

and

(49.3)LN=log⁡(1+2N).

Let θ∗,n satisfy

(49.4)Rn(θ∗,n)=Ln.

Assume

(49.5)Rn′(θ∗,n)=An>0.

Then Taylor’s theorem gives

(49.6)Δn,0(θ)=An(θ−θ∗,n)+O((θ−θ∗,n)2).

Hence

(49.7)Δn,0(θ)∼An(θ−θ∗,n)

as

(49.8)θ↓θ∗,n

from the right.

50. The Spatial Lambert Cascade

Combining the geometric expansion with the Puiseux theorem gives

(50.1)Δn,k(θ)∼21−2−k[An(θ−θ∗,n)]2−k.

Thus

(50.2)Δn,k(θ)∼21−2−kAn2−k(θ−θ∗,n)2−k.

Differentiating,

(50.3)∂θΔn,k∼21−2−k2−kAn2−k(θ−θ∗,n)2−k−1.

The logarithmic derivative is

(50.4)∂θlog⁡Δn,k∼2−kθ−θ∗,n.

Notice that the leading slope An cancels.

51. The Reduced Characteristic System

Consider the reduced transport equation

(51.1)bt+bbz=−b2∂zΔkΔk.

Define a characteristic Z(t) by

(51.2)Z˙(t)=b3(t),b3(t)=b(t,Z(t)).

Introduce

(51.3)θ(t)=t−Z(t).

Then

(51.4)θ˙=1−Z˙
(51.5)=1−b3.

Therefore

(51.6)θ˙=1−b3.

52. Evolution of the Characteristic Amplitude

Along the characteristic,

(52.1)b˙3=bt+bbz.

The PDE gives

(52.2)b˙3=−b32∂zΔkΔk.

Because

(52.3)θ=t−z,

we have

(52.4)∂z=−∂θ.

Therefore

(52.5)b˙3=b32∂θΔkΔk.

Equivalently,

(52.6)b˙3=b32∂θlog⁡Δk.

53. The Reciprocal Characteristic Variable

Set

(53.1)q(t)=1b3(t).

Then

(53.2)q˙=−1b32b˙3
(53.3)=−∂θlog⁡Δk.

Hence

(53.4)q˙=−∂θlog⁡Δk.

At the same time,

(53.5)θ˙=1−b3=1−1q=q−1q.

Thus

(53.6)θ˙=q−1q.

54. Derivation of the Exact Characteristic Invariant

We calculate

(54.1)dqdθ=q˙θ˙.

Using

(54.2)q˙=−∂θlog⁡Δk

and

(54.3)θ˙=q−1q,

we obtain

(54.4)dqdθ=−qq−1∂θlog⁡Δk.

Multiply by

(54.5)q−1q.

Then

(54.6)(1−1q)dq=−dlog⁡Δk.

Since

(54.7)(1−1q)dq=dq−dqq=dq−dlog⁡q,

we obtain

(54.8)dq−dlog⁡q=−dlog⁡Δk.

Integrating,

(54.9)q−log⁡q+log⁡Δk=C.

Since

(54.10)q=1b3,

we obtain the exact invariant

(54.11)1b3−log⁡1b3+log⁡Δk=C.

55. Characteristic Invariants, Critical Scaling, and Logarithmic Singularities

55.1. Explicit Form of the Characteristic Invariant

Exponentiating

(55.1)q−log⁡q+log⁡Δk=C

gives

(55.2)eqq−1Δk=eC.

Hence

(55.3)Δk=eCqe−q.

Define

(55.4)K=eC>0.

Then

(55.5)Δk=Kqe−q.

Since

(55.6)q=1b3,

we obtain

(55.7)Δk=Kb3e−1/b3.

56. The Blow-Up Branch

Suppose

(56.1)b3,0>1.

Then

(56.2)0<q0<1.

For

(56.3)0<q<1,

we have

(56.4)q−1q<0.

Therefore

(56.5)θ˙<0.

Thus the characteristic moves toward decreasing θ.

Suppose the branch point lies at

(56.6)θ∗<θ0.

Then the characteristic moves toward θ∗.

On the same branch,

(56.7)q↓0

corresponds to

(56.8)b3=1q⟶+∞.

From

(56.9)Δk=Kqe−q,

we have

(56.10)Δk⟶0⟺q⟶0.

Therefore

(56.11)Δk⟶0⟹b3⟶+∞.

57. Asymptotic Relation Between b3 and Δk

As

(57.1)q⟶0,

we have

(57.2)e−q=1−q+O(q2).

Therefore

(57.3)Δk=Kq(1−q+O(q2)).

Hence

(57.4)Δk=Kq+O(q2).

Thus

(57.5)q=ΔkK+O(Δk2).

Since

(57.6)b3=1q,

we obtain

(57.7)b3=KΔk+O(1).

In particular,

(57.8)b3∼KΔk.

58. Transfer of the Spatial Singularity to b3

Suppose

(58.1)Δk(θ)∼Ck(θ−θ∗)αk,αk=2−k.

Then

(58.2)b3∼KCk(θ−θ∗)−αk.

Therefore

(58.3)b3(θ)∼KCk(θ−θ∗)−2−k.

Thus the spatial Lambert singularity is transferred directly into the characteristic amplitude.

59. Finite-Time Arrival at the Branch Point

Let

(59.1)s(t)=θ(t)−θ∗.

Then

(59.2)s(t)>0

before blow-up and

(59.3)s(t)⟶0

at the singular time.

The characteristic equation gives

(59.4)s˙=1−b3.

Since

(59.5)b3⟶∞,

we have

(59.6)s˙∼−b3.

Using

(59.7)b3∼KCks−αk,

we obtain

(59.8)s˙∼−KCks−αk.

Therefore

(59.9)dt∼−CkKsαkds.

Integrating from s(t) to 0,

(59.10)T∗−t∼CkK∫0s(t)rαkdr
(59.11)=CkK(1+αk)s(t)1+αk.

Since

(59.12)αk=2−k>0,

the right-hand side is finite.

Thus

(59.13)T∗<∞.

60. Temporal Blow-Up Exponent

From

(60.1)T∗−t∼CkK(1+αk)s1+αk,

we obtain

(60.2)s∼[K(1+αk)Ck(T∗−t)]1/(1+αk).

Since

(60.3)b3∼KCks−αk,

we obtain

(60.4)b3(t)∼KCk[K(1+αk)Ck(T∗−t)]−αk/(1+αk).

Therefore

(60.5)b3(t)≍(T∗−t)−αk/(1+αk).

Since

(60.6)αk=2−k,

we obtain

(60.7)b3(t)≍(T∗−t)−1/(2k+1).

61. Finite-k Blow-Up Theorem

Theorem 58.1. Assume that for some finite integer k≥0:

  1. (a) Δk(θ)>0 for θ>θ∗ sufficiently close to θ∗;

  2. (b)

    (61.1)Δk(θ)=Ck(θ−θ∗)2−k(1+o(1));
  3. (c) b3,0>1;

  4. (d) the characteristic satisfies

    (61.2)θ˙=1−b3,b˙3=b32∂θlog⁡Δk.

Then the characteristic reaches θ∗ in a finite time T∗, and

(61.3)b3(t)⟶+∞

as t↑T∗. More precisely,

(61.4)b3(t)≍(T∗−t)−1/(2k+1).

Proof. The characteristic invariant gives

(61.5)Δk=Kqe−q,q=1b3.

Since b3,0>1, we have 0<q0<1. The characteristic branch satisfying q↓0 corresponds to θ↓θ∗.

As θ↓θ∗,

(61.6)Δk⟶0,

hence

(61.7)q⟶0,

and therefore

(61.8)b3⟶∞.

Using

(61.9)Δk∼Ck(θ−θ∗)2−k,

we obtain

(61.10)b3∼KCk(θ−θ∗)−2−k.

The equation

(61.11)θ˙=1−b3

therefore gives

(61.12)ddt(θ−θ∗)∼−KCk(θ−θ∗)−2−k.

The resulting time integral is

(61.13)∫0s0s2−kds,

which is finite.

Hence the characteristic reaches θ∗ in finite time.

The asymptotic integration gives

(61.14)b3(t)≍(T∗−t)−1/(2k+1).

◻

62. The Infinite Cascade and Its Limitation

For every fixed k,

(62.1)αk=2−k>0.

But

(62.2)αk⟶0.

Therefore

(62.3)12k+1⟶0.

Thus the finite-k temporal blow-up exponent does not converge to a nonzero Riccati exponent.

This is an important mathematical point.

The spatial singularity approaches

(62.4)(θ−θ∗)−1

in the sense that

(62.5)2−k−1⟶−1.

However, the corresponding characteristic exponent behaves differently because the relation between b3 and Δk must also be included.

Therefore one cannot simply substitute k=∞ into the finite-k temporal formula.

63. The Correct Double-Scaling Parameter

The cascade formula is

(63.1)Δk∼21−2−kΔ02−k.

Write

(63.2)Δ02−k=exp⁡(2−klog⁡Δ0).

Thus the relevant parameter is

(63.3)Λk=2−klog⁡1Δ0.

If

(63.4)Λk⟶0,

then

(63.5)Δ02−k⟶1.

If

(63.6)Λk⟶Λ∈(0,∞),

then

(63.7)Δ02−k⟶e−Λ.

If

(63.8)Λk⟶∞,

then

(63.9)Δ02−k⟶0.

Thus the infinite cascade has three distinct asymptotic regimes.

64. The Coupled Scaling k∼log2⁡n

Suppose now that

(64.1)Δn,0∼2−2nξ,ξ>0.

Then

(64.2)log⁡1Δn,0=2nlog⁡2−log⁡ξ+o(1).

Therefore

(64.3)Λn,k=2−k(2nlog⁡2−log⁡ξ+o(1)).

To keep Λn,k of order one, we require

(64.4)n2−k=O(1).

Thus

(64.5)2k≍n.

Equivalently,

(64.6)k=log2⁡n+O(1).

This is the critical cascade depth.

65. Uniformity Warning

The fixed-k theorem proves

(65.1)Δk=CkΔ02−k(1+O(Δ02−k))

for every fixed k.

It does not automatically prove that the error remains uniform as

(65.2)k⟶∞.

This distinction is essential.

For a full infinite-cascade theorem one would need estimates of the form

(65.3)|ΔkCkΔ02−k−1|≤E(Δ0,k),

with an explicit bound on E that remains controlled in the double-scaling regime

(65.4)k∼log2⁡n.

Such a bound requires iterating the error equation rather than merely iterating the leading asymptotic term.

66. A More Precise Logarithmic Form

The recurrence can be analyzed logarithmically.

Write

(66.1)Δk+1=2Δk(1+r(Δk)),

where

(66.2)r(Δ)=O(Δ).

Taking logarithms,

(66.3)log⁡Δk+1=12log⁡Δk+12log⁡2+log⁡(1+r(Δk)).

Therefore

(66.4)log⁡Δk+1=12log⁡Δk+12log⁡2+O(Δk).

Iterating this relation gives

(66.5)log⁡Δk=2−klog⁡Δ0+(1−2−k)log⁡2+Ek,

where

(66.6)Ek=∑j=0k−12−(k−1−j)O(Δj).

Hence

(66.7)log⁡Δk=2−klog⁡Δ0+(1−2−k)log⁡2+Ek.

Exponentiating,

(66.8)Δk=21−2−kΔ02−keEk.

Thus the leading formula is exact up to the multiplicative error eEk.

67. Control of the Iterated Error for Fixed k

For fixed k, we have

(67.1)Δj=O(Δ02−j).

Therefore

(67.2)Δj=O(Δ02−(j+1)).

Hence

(67.3)Ek=O(∑j=0k−12−(k−1−j)Δ02−(j+1)).

For fixed k, the dominant contribution comes from the largest exponent of Δ0, and in particular

(67.4)Ek=O(Δ02−k).

Therefore

(67.5)eEk=1+O(Δ02−k).

This recovers

(67.6)Δk=21−2−kΔ02−k(1+O(Δ02−k)).

68. The Geometric Profile and the Branch Location

Return to the corrected profile

(68.1)Rn(θ)=11−sin⁡θ+εn.

On the left of π/2 write

(68.2)θ=π2−x,x>0.

Then

(68.3)Rn=11−cos⁡x+εn=1εn+x2/2+O(x4).

The branch equation Rn=Ln therefore gives

(68.4)x∗,n=2(1Ln−εn)[1+o(1)],

with the condition

(68.5)Lnεn<1.

If Lnεn⟶λ<1, then

(68.6)θ∗,n−=π2−2(1−λ)Ln+o(Ln−1/2).

69. The Branch Slope

Differentiation gives

(69.1)Rn′(θ)=cos⁡θ(1−sin⁡θ+εn)2.

At the left branch,

(69.2)An:=Rn′(θ∗,n−)∼2Ln3/21−Lnεn>0.

This is the corrected branch slope used in the Taylor expansion of Δn,0.

70. The Complete Local Spatial Formula

Taylor expansion gives

(70.1)Δn,0=An(θ−θ∗,n−)+O((θ−θ∗,n−)2).

Therefore

(70.2)Δn,k=21−2−kAn2−k(θ−θ∗,n−)2−k(1+o(1)).

Thus

(70.3)Δn,k=Cn,k(θ−θ∗,n−)2−k(1+o(1)),

where

(70.4)Cn,k=21−2−kAn2−k.

Since the corrected construction takes uz=Δn,k, this formula is also the local physical z-velocity profile near the branch.

71. The Corresponding Logarithmic Singularity

Taking logarithms and differentiating yields

(71.1)∂θlog⁡Δn,k=2−kθ−θ∗,n−+o(1θ−θ∗,n−).

Thus the new reciprocal-sine input changes the branch location and slope but not the universal leading simple-pole coefficient generated by the renormalized Lambert cascade.

72. Trapping, Corrected Full PDE, and Continuity Reconstruction

72.1. Characteristic Trapping Lemma

Lemma 70.1. Assume

(72.1)b3,0>1,θ0>θ∗,

and suppose the renormalized b3-characteristic satisfies

(72.2)θ˙=1−b3,b˙3=b32∂θlog⁡Δk,

while remaining on the branch b3>1.

Then θ(t) is strictly decreasing. If it is bounded below by θ∗ and

(72.3)Δk(θ(t))⟶0,

then

(72.4)θ(t)⟶θ∗,b3(t)⟶+∞.

Proof. Since b3>1,

(72.5)θ˙=1−b3<0.

Thus θ(t) decreases. The exact invariant

(72.6)Δk=Kqe−q,q=1b3,

shows that on the branch 0<q<1,

(72.7)Δk⟶0⟹q⟶0.

Therefore

(72.8)b3=1q⟶+∞.

Because Δk(θ)>0 for θ>θ∗ and Δk(θ∗)=0, the limiting phase is θ∗.◻

73. Finite-Time Blow-Up Is a Consequence of the Invariant

Let

(73.1)s=θ−θ∗.

If

(73.2)Δk∼Cksαk,αk=2−k,

then the invariant gives

(73.3)b3∼KCk−1s−αk.

Hence

(73.4)s˙=1−b3=−KCk−1s−αk(1+o(1)).

Therefore

(73.5)dtds=−CkKsαk(1+o(1)).

Since

(73.6)∫0s0sαkds<∞,

the characteristic reaches s=0 at finite time. Thus

(73.7)b3(t)⟶+∞atT∗<∞

for every fixed finite cascade depth.

74. Interpretation of the Exponent Sequence

The recursion acts on the singular exponent according to

(74.1)αk+1=12αk.

Since

(74.2)α0=1,

we obtain

(74.3)αk=2−k.

Thus the map is

(74.4)α⟼α2.

The corresponding derivative exponent is

(74.5)αk−1.

Therefore

(74.6)0,−12,−34,−78,−1516,…

and

(74.7)αk−1⟶−1.

Thus the derivative singularity approaches the critical non-integrable exponent from above.

75. Why the Renormalized Cascade Is Different from Ordinary Lambert Composition

The distinction can now be stated precisely.

For ordinary composition,

(75.1)Yk+1=W0(Yk),

the branch point of W0 is fixed at

(75.2)−1e.

But

(75.3)W0(−1/e)=−1,

so the first iteration sends the branch point outside the domain of the next real composition.

For the renormalized map,

(75.4)Δk+1=1+W0(−e−1−Δk),

the branch point is translated:

(75.5)−1⟼0.

The next argument is reconstructed as

(75.6)−e−1−Δk+1.

Therefore the next Lambert operation again encounters the branch point at

(75.7)Δk+1=0.

This is why the square-root singularity is genuinely repeated.

76. The Renormalized Map as a Branch-Point Dynamical System

It is useful to regard

(76.1)Δk+1=F(Δk)

as a dynamical system.

The fixed point is

(76.2)Δ=0.

The map is continuous at zero:

(76.3)F(0)=0.

But

(76.4)F′(Δ)∼12Δ

diverges.

Thus zero is not a conventional hyperbolic fixed point.

It is a singular fixed point.

The local dynamics are

(76.5)Δk+1∼2Δk.

Hence sufficiently small Δk is transformed into a larger quantity.

For example, if

(76.6)0<Δ<1,

then

(76.7)Δ>Δ.

Thus the forward iteration moves away from zero.

This observation is important:

(76.8)The cascade creates singularity amplification, but the renormalized orbit does not approach zero under forward iteration.

The spatial singularity occurs because Δ0 approaches zero as a function of θ, while each finite cascade transforms the small-distance dependence through a fractional power.

77. This Distinction Is Essential

The statement

(77.1)Δk∼Δ02−k

does not mean that

(77.2)Δk⟶0

as k⟶∞ for a fixed Δ0∈(0,1).

Indeed,

(77.3)Δ02−k⟶1.

Thus

(77.4)Δk⟶2

at the level of the leading asymptotic coefficient.

The cascade exponent becomes singular because the dependence on Δ0 becomes flatter:

(77.5)Δ01/2k.

The corresponding derivative becomes more singular with respect to Δ0.

This is a crucial conceptual distinction between amplitude and sensitivity.

78. The Correct Interpretation of the Infinite Limit

For fixed Δ0>0,

(78.1)limk→∞Δ02−k=1.

Thus

(78.2)limk→∞21−2−kΔ02−k=2.

The singular behavior is instead seen in

(78.3)∂Δ0Δ02−k=2−kΔ02−k−1.

For fixed Δ0>0, this tends to zero.

But if

(78.4)Δ0⟶0

simultaneously with

(78.5)k⟶∞,

the result can be nontrivial.

This is exactly why the double-scaling variable

(78.6)2−klog⁡(1/Δ0)

is essential.

79. The Double-Scaling Limit

Suppose

(79.1)Δ0=Δ0(k)

and define

(79.2)Λk=2−klog⁡1Δ0(k).

Then

(79.3)Δ0(k)2−k=e−Λk.

Therefore

(79.4)Δk∼21−2−ke−Λk.

If

(79.5)Λk⟶Λ,

then

(79.6)Δk⟶2e−Λ.

Thus the infinite cascade can retain a nontrivial dependence only when the initial branch distance is exponentially small relative to the cascade depth.

80. Application to the n-Dependent Construction

Suppose

(80.1)Δn,0∼2−2nξn,

with

(80.2)ξn≍1.

Then

(80.3)Λn,k=2−klog⁡1Δn,0∼2−k(2nlog⁡2−log⁡ξn).

Hence

(80.4)Λn,k≍n2−k.

The nontrivial regime is

(80.5)n2−k≍1.

Therefore

(80.6)k=log2⁡n+O(1).

This is the correct interpretation of the proposed n-dependent amplification.

81. What Has Now Been Proved

The analysis above establishes the following rigorous chain.

  1. The renormalized map

    (81.1)F(Δ)=1+W0(−e−1−Δ)

    is real for every Δ≥0.

  2. It satisfies

    (81.2)F(0)=0.
  3. It maps

    (81.3)[0,∞)

    into

    (81.4)[0,1).
  4. Its exact inverse relation is

    (81.5)Δ=−log⁡(1−F(Δ))−F(Δ).
  5. Its local behavior is

    (81.6)F(Δ)=2Δ−23Δ+O(Δ3/2).
  6. For every fixed finite k,

    (81.7)Δk=21−2−kΔ02−k(1+O(Δ02−k)).
  7. The singular exponents therefore satisfy

    (81.8)αk=2−k.
  8. If

    (81.9)Δ0(θ)∼A(θ−θ∗),

    then

    (81.10)Δk(θ)∼21−2−kA2−k(θ−θ∗)2−k.
  9. The corresponding derivative satisfies

    (81.11)∂θΔk∼Ck2−kA2−k(θ−θ∗)2−k−1.
  10. The reduced characteristic system has the exact invariant

    (81.12)1b3−log⁡1b3+log⁡Δk=C.
  11. Consequently,

    (81.13)Δk=Kb3−1e−1/b3.
  12. Hence

    (81.14)b3∼K/Δk.
  13. For every fixed finite k, a characteristic with b3,0>1 reaches the branch point in finite time.

  14. The corresponding temporal blow-up rate is

    (81.15)b3(t)≍(T∗−t)−1/(2k+1).
  15. The infinite-cascade limit requires a coupled scaling and cannot be obtained by simply replacing k by ∞ in the finite-k formula.

  16. The natural n-dependent critical regime is

    (81.16)k∼log2⁡n

    when

    (81.17)Δn,0∼2−2nξn.

82. Final Theorem for the Renormalized Lambert Mechanism

We can summarize the rigorous reduced result as follows.

Theorem 80.1 (Renormalized Lambert–Characteristic Blow-Up Theorem). Let

(82.1)F(Δ)=1+W0(−e−1−Δ),Δ≥0,

and define

(82.2)Δk+1=F(Δk),Δ0>0.

Then:

  1. (i) The recursion is real for every finite k and

    (82.3)Δk∈[0,1)(k≥1).
  2. (ii) The branch point is a fixed point:

    (82.4)F(0)=0.
  3. (iii) The local expansion is

    (82.5)F(Δ)=2Δ−23Δ+O(Δ3/2).
  4. (iv) For every fixed finite k,

    (82.6)Δk=21−2−kΔ02−k(1+O(Δ02−k))

    as Δ0↓0.

  5. (v) Consequently, if

    (82.7)Δ0(θ)=A(θ−θ∗)(1+o(1)),A>0,

    then

    (82.8)Δk(θ)=21−2−kA2−k(θ−θ∗)2−k(1+o(1)).
  6. (vi) Hence

    (82.9)∂θlog⁡Δk=2−kθ−θ∗(1+o(1)).
  7. (vii) Consider the characteristic system

    (82.10)θ˙=1−b3,b˙3=b32∂θlog⁡Δk.

    Then

    (82.11)1b3−log⁡1b3+log⁡Δk=C.
  8. (viii) If b3,0>1 and the characteristic lies on the branch approaching θ∗, then

    (82.12)θ(t)↓θ∗,b3(t)⟶+∞.
  9. (ix) The singularity is reached in finite time T∗ and

    (82.13)b3(t)≍(T∗−t)−1/(2k+1).

Proof. Items (i)–(v) follow from the invariant-domain lemma, the Puiseux expansion, and the induction theorem.

Item (vi) follows by differentiating the spatial asymptotic formula.

Item (vii) follows from the reciprocal transformation q=1/b3 and direct integration of the characteristic equations.

Item (viii) follows from the branch choice b3,0>1, which gives θ˙<0, together with the invariant

(82.14)Δk=Kqe−q.

Item (ix) follows from

(82.15)Δk∼Ck(θ−θ∗)2−k,

which implies

(82.16)b3∼KCk(θ−θ∗)−2−k,

and therefore

(82.17)θ˙∼−KCk(θ−θ∗)−2−k.

The resulting time integral is finite and gives

(82.18)b3(t)≍(T∗−t)−1/(2k+1).

◻

83. The Full PDE analyses for b3 PDE according to the correct Equation (4.19) in Theoretical and Computational Fluid Mechanics Volume II 2026

This section collects the terminal b3 analysis in a compact form. The purpose is to retain the corrected PDE, the exact reconstruction identities, the characteristic mechanism, the error estimate needed for closure, and the finite-time terminal conclusion, while removing repeated intermediate reformulations.

83.1. Corrected extended b3 equation

Set

(83.1)b:=b3=uxuy,A:=uy,C:=uz,ux=bA.

The corrected forcing inherited from Eq. (4.19) is

(83.2)Fc=2C[A0+(sin⁡xsin⁡ysin⁡z)2].

Thus its sign is controlled by the sign of C=uz. In the terminal branch used below we assume

(83.3)C≥c∗>0.

Let

(83.4)Q=ΔN,k(ξ),ΓN,k:=∂ξlog⁡Q=Q′Q,ξ=C4t−C3z−C2y−C1x.

The principal nonlinear source is therefore ΓN,kb2. With the notation

(83.5)Xh:=A∂x+bA∂y,D:=∂t+b∂z+C2Xh,

the extended equation can be organized as

(83.6)Db=ΓN,kb2+Fc−RN,

where RN contains the reconstruction terms produced by C and the transverse geometry. In the notation used in the terminal reconstruction,

(83.7)RN=2bXhC+bCt+CzC,

up to the already absorbed terms in the chosen corrected formulation. The analysis below requires only the normalized estimate

(83.8)|RN|ΓN,kb2⟶0(s↓0).

This is the single perturbative statement needed for the terminal comparison.

83.2. Exact incompressible reconstruction

The two-correction representation may be written

(83.9)A=Q+S−D,C=Q+S+D,

with ux=b/A. It preserves b=uxuy identically. More importantly, incompressibility is imposed exactly:

(83.10)∂x(bA)+Ay+Cz=0.

Thus C is not an independent terminal ansatz; once b and A are fixed, it is reconstructed from

(83.11)Cz=−∂x(bA)−Ay.

Periodic reconstruction additionally requires the corresponding zero-mean compatibility in z.

For the square-root specialization used in the terminal tube,

(83.12)A2=−2χ(x)b,χ<0,

so that

(83.13)uy=A,ux=bA,uy2=−2χb,ux2=−b2χ.

Here and below νx,νy denote the corresponding velocity components ux,uy; this notation is used only in this displayed identity to avoid ambiguity with the viscosity coefficient ν.

If the distinguished phase is

(83.14)s=t−z−θ∗,N,

then functions depending on (t,z) only through s satisfy

(83.15)(∂t+∂z)f(s)=0.

This cancellation is the main reason that the Ct+Cz part of (83.7) is small in the terminal branch.

83.3. Lambert branch scale and tuned terminal regime

The renormalized Lambert cascade is

(83.16)ΔN,j+1=1+W0(−e−1−ΔN,j).

Near the branch point,

(83.17)1+W0(−e−1−d)=2d−23d+O(d3/2),d↓0.

At fixed finite depth k this gives

(83.18)QN(s)≍s2−k,ΓN,k(s)∼2−ks.

For k=3,

(83.19)QN≍s1/8,ΓN,3∼18s.

The logarithmic singularity in ΓN,k is therefore generated directly by the Lambert branch geometry.

The tuning of the reconstructed branch is chosen so that the principal phase dynamics and the quadratic source have the same sign in the terminal region. The corrected forcing (83.2) is then a positive seed when C>0, but it is not the leading terminal term. Indeed, once b is large,

(83.20)Fc=O(C),ΓN,kb2≍b2s,

so the forcing is lower order under the terminal estimates below. Its principal role is to assist entry into the positive-amplitude regime.

83.4. Characteristic reduction and terminal trapping

Let (X(t),Y(t),Z(t)) be a characteristic of D:

(83.21)X˙=C2A,Y˙=C2bA,Z˙=b.

Along it, (83.6) becomes

(83.22)b˙=ΓN,kb2+Fc−RN.

For the distinguished phase (83.14),

(83.23)s˙=1−b.

Hence b>1 implies s˙<0.

Assume that after an entry time te,

(83.24)b(te)>1,Fc≥0,|RN|≤ϑΓN,kb2,0<ϑ<1.

Then

(83.25)b˙≥(1−ϑ)ΓN,kb2>0,s˙=1−b<0.

The terminal region is therefore self-trapping: b increases while s decreases.

Using (83.18),

(83.26)−dbds=b˙b−1≥(1−ϑ)2−ksb2b−1.

Since b/(b−1)>1 for b>1,

(83.27)−ddslog⁡b≥(1−ϑ)2−ks.

Integration from (se,be) yields the useful lower bound

(83.28)b(s)≥be(ses)(1−ϑ)2−k.

Thus b(s)⟶+∞ as s↓0 within the stated hypotheses.

The physical time is finite. Since

(83.29)dt=−dsb−1,

and (83.28) gives a positive power of s−1,

(83.30)T∗,N−te=∫0sedsb(s)−1<∞.

Therefore

(83.31)T∗,N<∞,s(t)↓0,QN(t)↓0,b(t)↑+∞

for every finite configuration satisfying the terminal assumptions.

83.5. The single reconstruction estimate required for closure

The lengthy reconstruction analysis reduces to proving (83.8). In the distinguished phase, the first part of (83.7) vanishes when C depends on (t,z) only through s:

(83.32)Ct+Cz=0.

The remaining term is

(83.33)RN=2bXhC.

Consequently,

(83.34)|RN|ΓN,kb2≤2|XhC|ΓN,kb.

The normalized reconstruction estimate obtained in the terminal tube is

(83.35)|XhC|ΓN,kb=O(s)+O(s3/2),

and hence

(83.36)|RN|ΓN,kb2=O(s)⟶0.

This estimate is self-strengthening: once the characteristic enters smaller s, the relative reconstruction error becomes smaller.

More generally, if the exact cancellation (83.32) is replaced by a perturbed reconstruction, it is sufficient to verify

(83.37)|Ct+Cz|CΓN,kbo0,|XhC|ΓN,kbo0.

These two dimensionless estimates are the compact replacement for the repeated term-by-term remainder calculations.

83.6. Entry, logarithmic branch distance, and bootstrap

A convenient branch variable is the logarithmic distance

(83.38)q:=−log⁡Q.

Since Q↓0 corresponds to q↑∞, the terminal problem can be separated into an entry regime and an amplitude-dominated regime. The corrected forcing (83.2) supplies a positive seed on the C>0 branch. Once b reaches a fixed threshold be>1, the quadratic term dominates and the bootstrap (83.24)–(83.25) applies.

A minimal bootstrap formulation is therefore:

(83.39)C≥c∗>0,b≥be>1,ΓN,k>0,|RN|≤ϑΓN,kb2,0<ϑ<1.

Under these assumptions the inequalities in (83.25) preserve the bootstrap until s=0. No separate large-amplitude reconstruction is required beyond the uniform validity of the normalized error bound.

83.7. Two-sided terminal scale and finite-depth exponent

If, in addition to the lower comparison, one has a matching upper estimate

(83.40)(1−ϑ)ΓN,kb2≤b˙≤(1+ϑ)ΓN,kb2

in the terminal region, then b is trapped between nearby powers of s. Formally neglecting the vanishing relative error gives

(83.41)−dbds∼2−ksb,b(s)≍s−2−k.

Substitution into dt=−ds/(b−1) gives

(83.42)T∗−t≍s1+2−k,

and therefore

(83.43)b(t)≍(T∗−t)−1/(2k+1).

For k=3 this gives the finite-depth formal terminal law

(83.44)b(s)≍s−1/8,b(t)≍(T∗−t)−1/9.

The rigorous lower comparison (83.28) should be distinguished from the sharper asymptotic (83.41), which requires the corresponding two-sided control.

83.8. Terminal velocity geometry

In the parabolic terminal tube,

(83.45)χ2≍s.

Combining (83.12) with a general amplitude law b≍s−α gives

(83.46)A=uy≍s1/4−α/2,ux=ux=bA≍s−1/4−α/2.

For the k=3 scale α=1/8,

(83.47)b≍s−1/8,uy≍s3/16,νx≍s−5/16.

Hence

(83.48)νy⟶0,|νx|⟶∞,νxνy=b⟶+∞.

The exponents are consistent because

(83.49)−516+316=−18.

Moreover, |χ|≍s1/2 implies

(83.50)|χνx|≍s3/16⟶0,

so the square-root relation remains compatible with the divergence of ux.

The component C=uz is reconstructed from (83.11). In the local square-root core the resulting form is

(83.51)C(x,s)=C0(x)+χ′(x)−2Bχ(x)4χ(x)2∫0sr−α/2F(r)1/2dr,

where F contains the bounded correction to the leading b profile. At k=3 this gives, along the terminal tube,

(83.52)C=uz=c0+O(s3/16),c0>0,

which is consistent with the positivity assumption used in the corrected forcing.

83.9. Critical coupled depth

The preceding terminal argument is a finite-depth statement. The cascade itself has a distinct coupled-depth transition. If ΔN,0∼λN−2, then the characteristic depth is of order log2⁡log⁡N. In particular,

(83.53)kN−log2⁡log⁡N⟶−∞⟹ΔN,kN⟶0,

whereas

(83.54)kN−log2⁡log⁡N⟶+∞⟹ΔN,kN⟶Δ∗>0.

The positive fixed point satisfies

(83.55)(1−Δ∗)e2Δ∗=1,Δ∗=0.796812130020020….

Thus the finite-k terminal law and a coupled N,k⟶∞ limit are different questions. The terminal theorem below applies to each finite configuration satisfying its hypotheses; it does not by itself identify every coupled infinite-depth limit.

83.10. Compact terminal theorem

Theorem 81.1 (Finite-depth terminal amplification). Fix finite N and finite cascade depth k. Assume that on a terminal characteristic interval:

  1. QN>0 and ΓN,k∼2−k/s with s=t−z−θ∗,N>0;

  2. C=uz≥c∗>0 and the corrected forcing (83.2) is nonnegative;

  3. b(te)=be>1 at some entry time te;

  4. the reconstruction satisfies

    |RN|≤ϑΓN,kb2,0<ϑ<1.

Then the characteristic remains trapped in the terminal regime and

(83.56)b˙≥(1−ϑ)ΓN,kb2>0,s˙=1−b<0.

Moreover,

(83.57)b(s)≥be(ses)(1−ϑ)2−k,

and the characteristic reaches s=0 in finite physical time. Consequently,

(83.58)QN(t)⟶0,b(t)⟶+∞as t↑T∗,N<∞.

If the normalized remainder tends to zero two-sidedly so that (83.40) holds with arbitrarily small ϑ, then the sharper finite-depth scaling is

(83.59)b(s)≍s−2−k,b(t)≍(T∗,N−t)−1/(2k+1).

Proof. The first inequality follows immediately from (83.22), nonnegativity of Fc, and the remainder hypothesis. Since be>1, monotonicity gives b>1 thereafter, and (83.23) yields s˙<0. Dividing the two characteristic equations gives (83.26); integration gives the stated lower bound. The time integral (83.30) converges. The sharper scaling follows from the two-sided comparison and integration of the asymptotic logarithmic equation.◻

83.11. What remains distinct from the scalar terminal mechanism

The preceding theorem isolates the scalar characteristic mechanism and the reconstruction estimate required for it. For the full three-dimensional Navier–Stokes conclusion, three further issues remain logically separate and should not be hidden inside the scalar comparison:

(83.60)(i)uniform periodic reconstruction and derivative bounds,(ii)the Leray--Hopf energy/dissipation requirements,(iii)one common periodic pressure satisfying all three momentum equations.

The moving-coterminal analysis developed earlier in this paper addresses the local terminal Lt2Hx1 obstruction by replacing the nonmoving wedge extending to ρ=0 with a moving preterminal core. That construction is complementary to the present b3 characteristic theorem: the latter supplies the scalar amplification mechanism, whereas the former controls the local energy geometry of the reconstructed velocity. Neither step should be used as a substitute for the remaining global pressure and periodic reconstruction checks.

84. Lower-Barrier Viscosity Analysis and the Favorable Laplacian Sign

The viscous term requires a separate treatment in the lower-barrier argument. In particular, it is not necessary to prove that viscosity is asymptotically smaller than the Lambert production term in order to preserve the lower barrier. What is required is the sign of the viscous contribution at a hypothetical first downward contact between the solution and the barrier.

This distinction is important. A comparison of the absolute sizes of

νΔbandΓN,kb2

does not by itself determine whether viscosity destroys the lower-barrier argument. The relevant quantity at first contact is instead

Δ(b−b―),

and this quantity has a favorable sign.

We first state the argument for the ordinary Laplacian term and then record the corresponding calculation for the earlier transformed viscosity operator.

84.1. The lower barrier

Let

s=s(x,z,t)>0

denote the terminal phase variable and suppose that, in the terminal region, the characteristic geometry satisfies

(84.1)−Ds=HN+γNb,γN>0.

Assume also that the Lambert coefficient obeys

(84.2)ΓN,k≥G−s,G−>0,

throughout the comparison region.

After the reconstruction errors have been estimated, suppose that the scalar equation gives the differential inequality

(84.3)Db≥(1−θR)ΓN,kb2+νΔb,0<θR<1.

We introduce the lower barrier

(84.4)b―(s)=Bs−λ,B>0,λ>0.

Along the characteristic,

Db―=−λBs−λ−1Ds.

Using (84.1),

(84.5)Db―=λBs−λ−1(HN+γNb).

At a contact point at which

b=b―=Bs−λ,

the dominant transport contribution is therefore

(84.6)Db―=λγNB2s−2λ−1+λBHNs−λ−1.

On the other hand, by (84.2),

(84.7)(1−θR)ΓN,kb―2≥(1−θR)G−B2s−2λ−1.

Consequently, the Lambert production dominates the leading characteristic motion of the barrier whenever

(84.8)λγN<(1−θR)G−.

Equivalently, it is sufficient to choose

(84.9)0<λ<(1−θR)G−γN.

The remaining issue is the sign of viscosity at a first downward contact.

84.2. First-contact geometry

Define

(84.10)w=b−b―.

Suppose initially

w>0

and assume, for contradiction, that the solution crosses below the barrier. Let (x∗,t∗) be the first point of downward contact. Then

(84.11)w(x∗,t∗)=0,∇w(x∗,t∗)=0,Δw(x∗,t∗)≥0.

Hence

(84.12)b=b―,∇b=∇b―,

and, crucially,

(84.13)Δb≥Δb―.

Thus the viscous term does not have to be discarded or treated as a small error. At the first downward contact it contributes in the favorable direction provided the barrier itself has positive Laplacian.

84.3. Positive Laplacian of the power barrier

For the distinguished terminal phase

(84.14)s=t−z−θ∗,

we have

(84.15)sz=−1,szz=0.

More generally, if the phase is affine in the active spatial direction and

(84.16)|∇s|2=qs2>0,Δs=0,

then

Δb―=b―ss|∇s|2.

Since

b―=Bs−λ,

we obtain

(84.17)b―s=−λBs−λ−1

and

(84.18)b―ss=Bλ(λ+1)s−λ−2>0.

Therefore

(84.19)Δb―=qs2Bλ(λ+1)s−λ−2>0.

For the distinguished normalization qs2=1,

(84.20)νΔb―=νBλ(λ+1)s−λ−2>0.

Combining (84.13) and (84.20) gives

(84.21)νΔb≥νΔb―=νBλ(λ+1)s−λ−2>0

at a first downward contact.

This is the favorable viscosity sign needed in the comparison argument.

Notice that (84.21) is a first-contact statement. It does not assert that

Δb≥0

at every point of the terminal region. Such a global assertion is neither needed nor implied by the comparison argument.

84.4. The k=3 exponent

For the finite-depth k=3 terminal scaling,

(84.22)ΓN,3∼18s.

The natural power appearing in the scalar calculation is

(84.23)λ=18.

For this exponent,

λ(λ+1)=1898=964.

Hence

(84.24)Δb―=9B64s−17/8,

and therefore

(84.25)νΔb―=9νB64s−17/8>0.

Thus the fact that the Laplacian has the stronger singular order

s−17/8

than the Lambert source

(84.26)ΓN,3b―2∼B28s−5/4

does not by itself invalidate the lower-barrier argument.

Indeed, the comparison mechanism is not based upon

|νΔb―|≪|ΓN,3b―2|.

It is based upon the fact that, at first downward contact,

(84.27)νΔb≥νΔb―>0.

Hence the viscous contribution acts in the same direction as the desired lower-barrier inequality.

84.5. The earlier transformed viscosity

The same point can be formulated using the earlier two-stage product-variable viscosity structure.

Let

(84.28)b=b3=uxuy,C=uz.

For the earlier transformed operator, the third component is

(84.29)V3=νΔ(Cb)−2ν(G2)3.

With

(84.30)(G2)3=(1+C)∇C⋅∇b,

we obtain

(84.31)V3=νΔ(Cb)−2ν(1+C)∇C⋅∇b.

After division by C>0, define

(84.32)V[b,C]=V3C.

Using

(84.33)Δ(Cb)=CΔb+bΔC+2∇C⋅∇b,

we find

(84.34)V[b,C]=ν[Δb+bCΔC−2∇C⋅∇b].

The importance of this representation is that the lower-order terms are linear in b and ∇b. Therefore the first-contact subtraction is particularly simple.

At a first downward contact set

w=b−b―.

Then

(84.35)w=0,∇w=0,Δw≥0.

Subtracting the transformed viscosity evaluated on the barrier gives

(84.36)V[b,C]−V[b―,C]=νΔw+νΔCCw−2ν∇C⋅∇w.

At the contact point, the last two terms vanish exactly. Hence

(84.37)V[b,C]−V[b―,C]=νΔw≥0.

Consequently,

(84.38)V[b,C]≥V[b―,C]

at first downward contact.

Thus the transformed viscosity preserves the same comparison structure as the ordinary Laplacian. It remains only to determine the sign of V[b―,C].

84.6. Evaluation of the transformed viscosity on the barrier

For

b―=Bs−λ,

equation (84.34) gives

(84.39)V[b―,C]=νΔb―+νb―CΔC−2ν∇C⋅∇b―.

For the affine phase s=t−z−θ∗,

(84.40)b―z=λBs−λ−1,

and therefore

(84.41)V[b―,C]=νBλ(λ+1)s−λ−2+νBs−λΔCC−2νBλCzs−λ−1.

The first term is strictly positive:

(84.42)νBλ(λ+1)s−λ−2>0.

The reconstructed vertical component in the square-root terminal tube has the form

(84.43)C=c0+O(s3/16),c0>0.

Along the parabolic terminal scaling, the previously derived derivative order is

(84.44)Cz=O(s−13/16).

Therefore the mixed term satisfies

(84.45)Czs−λ−1=O(s−λ−29/16).

By contrast, the positive Laplacian contribution has order

(84.46)s−λ−2=s−λ−32/16.

Hence

(84.47)s−λ−32/16≫s−λ−29/16(s↓0).

Thus the gradient correction is lower order than the positive barrier Laplacian.

If, in addition, the reconstructed C-field satisfies the uniform terminal estimate

(84.48)ΔCC=o(s−2),

then

(84.49)νBs−λΔCC=o(s−λ−2).

Equations (84.41)–(84.49) therefore imply

(84.50)V[b―,C]=νBλ(λ+1)s−λ−2[1+o(1)]>0.

At k=3, with λ=1/8,

(84.51)V[b―,C]∼9νB64s−17/8>0.

Combining this with (84.38),

(84.52)V[b,C]≥V[b―,C]>0

at a sufficiently terminal first downward contact.

84.7. Closure of the first-contact argument

We may now combine the production and viscosity calculations.

At a hypothetical first downward contact,

(84.53)b=b―,∇(b−b―)=0,Δ(b−b―)≥0.

The scalar inequality gives

(84.54)Db≥(1−θR)ΓN,kb―2+V[b,C].

By (84.38),

(84.55)Db≥(1−θR)ΓN,kb―2+V[b―,C].

For sufficiently small s, (84.50) gives

(84.56)V[b―,C]>0.

Consequently,

(84.57)Db>(1−θR)ΓN,kb―2.

On the other hand, the barrier derivative is given by (84.5). Therefore, after controlling the lower-order HN-term, condition

(84.58)λ<(1−θR)G−γN

gives

(84.59)Db>Db―

at the hypothetical first downward contact.

But a first crossing from

b>b―

to

b<b―

requires

(84.60)D(b−b―)≤0.

Equations (84.59) and (84.60) contradict one another. Hence the downward crossing cannot occur within the region in which the stated hypotheses hold.

We therefore obtain the conditional lower-barrier conclusion

(84.61)b(s)≥Bs−λ,0<λ<(1−θR)G−γN,

provided the reconstruction estimates, positivity of C, and the uniform derivative condition (84.48) hold throughout the comparison region.

84.8. Interpretation and remaining uniform estimate

The preceding calculation clarifies the role of viscosity in the terminal comparison argument.

It is incorrect to infer

“viscosity is more singular than the Lambert source”

and conclude from this alone that

“viscosity destroys the lower-barrier mechanism.”

The relevant comparison is instead

(84.62)V[b,C]−V[b―,C]=νΔ(b−b―)≥0

at first downward contact, together with

(84.63)V[b―,C]>0.

Thus the singular positive Laplacian of the power barrier is favorable to the lower-bound comparison.

The remaining issue is not the sign of the principal barrier Laplacian. That sign is explicit:

(84.64)Δ(Bs−λ)=Bλ(λ+1)s−λ−2>0.

Rather, for the full multidimensional transformed operator, the remaining uniform reconstruction task is to justify estimates of the form

(84.65)C≥c−>0,ΔCC=o(s−2),Cz=O(s−13/16)

throughout the complete terminal comparison tube, and not merely along a single distinguished parabolic path.

Accordingly, the rigorous logical structure is

(84.66)Lambert production lower bound+first-contact minimum principle+Δb―>0+uniform control of the reconstructed C-terms⟹preservation of the lower barrier.

This is the form of the viscosity argument required in the terminal lower-barrier analysis. It preserves the favorable sign calculation without making the stronger and unnecessary assertion that viscosity is a perturbatively small term relative to the Lambert production.

85. Conclusion

The corrected terminal architecture can be summarized without the repeated developmental variants:

(85.1)Lambert branch geometry⟹QN↓0,ΓN,k∼2−ks,corrected b3-PDE⟹b˙≥(1−ϑ)ΓN,kb2,b>1⟹s˙=1−b<0,s↓0⟹|RN|ΓN,kb2↓0,henceb↑∞in finite characteristic time.

For k=3, the associated formal finite-depth terminal hierarchy is

(85.2)QN≍s1/8,ΓN,3≍s−1,b3≍s−1/8,

with the square-root reconstruction

(85.3)uy≍s3/16,νx≍s−5/16,C=uz=c0+O(s3/16).

The essential point is the separation of roles: the cascade determines the branch scale, ΓN,k drives the scalar amplitude, the phase equation provides trapping, exact incompressibility reconstructs uz, and the square-root relation determines the anisotropic ux–uy geometry. This condensed formulation retains the mathematical content needed for the terminal argument while removing the repeated derivations that led to it.

86. Appendix A: Lambert–W Branch Expansion

Theorem 83.1 (Square-root expansion of the renormalized Lambert map). Define

F(Δ)=1+W0(−e−1−Δ),Δ>0.

As Δ↓0,

F(Δ)=2Δ−23Δ+192Δ3/2+O(Δ2).

In particular,

F(Δ)∼2Δ,

and there exist constants c1,c2,Δ∗>0 such that

c1Δ≤F(Δ)≤c2Δ,0<Δ<Δ∗.

Thus one application of the renormalized principal-branch Lambert map halves the small-Δ exponent. This is the branch-point estimate used in the finite-depth cascade analysis.

87. Appendix B: Fixed-Depth Renormalized Cascade

Theorem 84.1 (Finite-depth cascade asymptotic). Let

Δj+1=1+W0(−e−1−Δj),Δ0↓0.

For every fixed finite integer k≥1,

Δk=21−2−kΔ02−k(1+O(Δ02−k)).

Consequently, the exponent cascade is

1⟶12⟶14⟶18⟶⋯⟶2−k,

while the corresponding leading prefactors are

1,21/2,23/4,27/8,…,21−2−k.

The statement is a fixed-finite-depth asymptotic: the constants implicit in the error term may depend on k. It therefore does not by itself establish a uniform asymptotic in a coupled regime k=k(Δ0)⟶∞.

88. Appendix C: Characteristics of the Reduced b3-Equation

Theorem 85.1 (Characteristic system and distinguished phase). Let

A=uy,C=uz,b=b3=uxuy,

and

b1=uyuz=AC,b2=uxuz=bCA.

For the transport operator

∂t+uzb1∂x+uzb2∂y+b∂z,

the characteristic curves satisfy

t˙=1,x˙=AC2,y˙=bC2A,z˙=b.

Accordingly, along each characteristic,

dbdτ=bt+uzb1bx+uzb2by+bbz,

so the characteristic derivative reproduces exactly the transport part of the reduced b3-equation. For the distinguished phase

θ=t−z,

one has independently of the transverse characteristic motion

θ˙=t˙−z˙=1−b.

This phase identity is the characteristic relation used in the terminal trapping analysis.

References

26 Cites in Article
  1. 1. Terry Moschandreou (2026). Exploration of Finite Time Singularities of the 3D Navier Stokes Equations over a Periodic Domain $mathbb{T}^3$. Global Journal of Science Frontier Research
  2. 2. T. Moschandreou,K. Afas,K. Nguyen,K. Kyritsis (2025). Theoretical and Computational Fluid Mechanics: Existence, Blowup and Discrete Exterior Calculus Problems, Volume II, Chapman & Hall/CRC, Numerical Analysis and Scientific Computing Series, 2025, ISBN 978-1-032-98931-0. Theoretical and Computational Fluid Mechanics: Existence, Blowup and Discrete Exterior Calculus Problems, 978-1.
  3. 3. Terence Tao (2016). Finite time blowup for an averaged three-dimensional Navier-Stokes equation. Journal of the American Mathematical Society, 29(3), 601-674.
  4. 4. Terry Moschandreou (2021). No Finite Time Blowup for 3D Incompressible Navier Stokes Equations via Scaling Invariance. Mathematics and Statistics, 9(3), 386-393.
  5. 5. L. Caffarelli,R. Kohn,L. Nirenberg (1982). Partial regularity of suitable weak solutions of the navier‐stokes equations. Communications on Pure and Applied Mathematics, 35(6), 771-831.
  6. 6. Guo Luo,Thomas Hou (2014). Potentially singular solutions of the 3D axisymmetric Euler equations. Proceedings of the National Academy of Sciences, 111(36), 12968-12973.
  7. 7. T. Hou,R. Li (2006). Dynamic strain-induced singularity formation and its opposition to blow-up for the 3D incompressible Euler equations. Journal of Nonlinear Science, 16(3), 223-267.
  8. 8. J. Chen,T. Hou (2022). Stable Nearly Self-Similar Blowup of the 2D Boussinesq and 3D Euler Equations with Smooth Data II: Rigorous Numerics. Communications on Pure and Applied Mathematics, 23(4), 2343-2426.
  9. 9. T. Moschandreou (2026). Existence of a Periodic Attractor for the 3D Navier–Stokes Equations on T3 , Advances in Pure Mathematics, 2026. Advances in Pure Mathematics
  10. 10. Jean Leray (1934). Sur le mouvement d'un liquide visqueux emplissant l'espace. Acta Mathematica, 63, 193-248.
  11. 11. E. Hopf (1951). Über die Anfangswertaufgabe für die hydrodynamischen Gleichungen. Mathematische Nachrichten, 4(1-6), 213-231.
  12. 12. James Serrin (1962). On the interior regularity of weak solutions of the Navier-Stokes equations. Archive for Rational Mechanics and Analysis, 9(1), 187-195.
  13. 13. J. Beale,T. Kato,A. Majda (1984). Remarks on the breakdown of smooth solutions for the 3-D Euler equations. Communications in Mathematical Physics, 94(1), 61-66.
  14. 14. W. Reid (1972). Riccati Differential Equations.
  15. 15. A. Polyanin,V. Zaitsev (2003). Handbook of Exact Solutions for Ordinary Differential Equations, Chapman and Hall/CRC, Boca Raton, 2nd edition, 2003.
  16. 16. E. Kamke (1959). Differentialgleichungen: Lösungsmethoden und Lösungen, Akademische Verlagsgesellschaft, Leipzig, 6th edition, 1959. Differentialgleichungen: Lösungsmethoden und Lösungen
  17. 17. E. Coddington,N. Levinson (1955). Theory of Ordinary Differential Equations, McGraw-Hill, New York, 1955.
  18. 18. Lawrence Evans (2010). Four important linear partial differential equations. Graduate Studies in Mathematics, 19, 17-90.
  19. 19. A. Majda,A. Bertozzi (2002). Vorticity and Incompressible Flow.
  20. 20. R. Temam (2001). Navier–Stokes Equations: Theory and Numerical Analysis.
  21. 21. G. Prodi (1959). Un teorema di unicità per le equazioni di Navier-Stokes. Rendiconte del Seminario Matematico della Università di Padova, 29, 269-274.
  22. 22. O. Ladyzhenskaya (1969). The Mathematical Theory of Viscous Incompressible Flow, Gordon and Breach, New York, 2nd edition, 1969.
  23. 23. L. Escauriaza,G. Seregin,V. Šverák (2003). L 3,∞ -solutions of the Navier-Stokes equations and backward uniqueness. Russian Mathematical Surveys, 58(2), 211-250.
  24. 24. Peter Constantin,Charles Fefferman (1993). Direction of vorticity and the problem of global regularity for the 3D Navier–Stokes equations. Indiana University Mathematics Journal, 42(3), 775.
  25. 25. R. Corless,G. Gonnet,D. Hare,D. Jeffrey,D. Knuth (1996). On the LambertW function. Advances in Computational Mathematics, 5(1), 329-359.
  26. 26. Haim Brezis (2011). Functional Analysis, Sobolev Spaces and Partial Differential Equations. Functional Analysis, Sobolev Spaces and Partial Differential Equations

Funding

No external funding was declared for this work.

Conflict of Interest

The authors declare no conflict of interest.

Ethical Approval

No ethics committee approval was required for this article type.

Data Availability

Not applicable for this article.

How to Cite This Article

Terry Moschandreou. 2026. "Renormalized Lambert–W Cascade and Finite-Time Amplification and Blowup for moving coterminal reconstructed b Dynamics on T 3 for the 3D Navier Stokes Equations". Global Journal of Science Frontier Research - F: Mathematics & Decision GJSFR-F Volume 26 (N/A).

More Citation Formats

Download Citation

Journal Specifications

Crossref Journal DOI 10.17406/GJSFR

Print ISSN 0975-5896

e-ISSN 2249-4626

Keywords
Classification
MSC 35Q30
MSC 35B44
Version of record

v1.2

Language
English
Order Article Reprint
Experience in AR

Explore published articles in an immersive Augmented Reality environment. Our platform converts research papers into interactive 3D books, allowing readers to view and interact with content using AR and VR compatible devices.

Read in 3D

Your published article is automatically converted into a realistic 3D book. Flip through pages and read research papers in a more engaging and interactive format.

Article Matrices
Total Views: 35
Total Downloads: 0
All Trends

Request Access

Please fill out the form below to request access to this research paper. Your request will be reviewed by the editorial or author team.
X

This is the heading

Lorem ipsum dolor sit amet, consectetur adipiscing elit. Ut elit tellus, luctus nec ullamcorper mattis, pulvinar dapibus leo.

High-quality academic research articles on global topics and journals.

Renormalized Lambert–W Cascade and Finite-Time Amplification and Blowup for moving coterminal reconstructed b Dynamics on T 3 for the 3D Navier Stokes Equations

Terry Moschandreou
Terry Moschandreou <p>Intermediate Science and Mathematics</p>