Is it true that no mathematical relation exists between the Navier-Stokes equations and the multifractal model?

Dario Vincenzi; John D. Gibbon

arxiv: 2603.19125 · v2 · submitted 2026-03-19 · ⚛️ physics.flu-dyn · nlin.CD

Is it true that no mathematical relation exists between the Navier-Stokes equations and the multifractal model?

John D. Gibbon , Dario Vincenzi This is my paper

Pith reviewed 2026-05-15 08:28 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn nlin.CD

keywords Navier-Stokes equationsmultifractal modelLeray weak solutionsturbulence scalingReynolds numberinverse scalevelocity gradient norms

0 comments

The pith

A correspondence between velocity gradient norms and local scaling exponents reconciles the Navier-Stokes equations with the multifractal model.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper challenges the accepted view that the Navier-Stokes equations and the multifractal model lack any mathematical connection. It establishes a link by using L to the power of 2m norms of the velocity gradient in Leray weak solutions to identify a parameter m with the local scaling exponent h from the multifractal picture. This mapping produces the Paladin-Vulpiani inverse scale, which ties directly to the Reynolds number and serves as a bridge between the two frameworks. The approach treats m as a variable focus that selects different turbulent structures across a specific range of exponents.

Core claim

From Euler invariant scaling combined with the Navier-Stokes equations in a three-dimensional box of size L, the L^{2m}-norms of the velocity gradient establish a correspondence between m and the local multifractal exponent h. This yields the Paladin-Vulpiani inverse scale satisfying L over eta sub h,pav equals Re to the power 1 over (1 plus h), which mediates between the equations and the multifractal model for Leray weak solutions. The interval from m equals 1 to infinity maps to h from negative two-thirds to one-third.

What carries the argument

The mapping from the continuous parameter m in the L^{2m}-norms of the velocity gradient to the local scaling exponent h in the multifractal model, which functions as a sliding focus to select structures at different scales.

If this is right

The derived scale connects the box size L to Reynolds-number scaling in a way that is consistent with Euler invariants.
The range of h values overlaps the interval where thermal noise is expected to dominate and induce spontaneous stochasticity.
The parameter m provides a continuous way to probe different local structures within the same weak solution.
The reconciliation implies that multifractal predictions can be tested against properties of Leray solutions without assuming stronger smoothness.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

Simulations could vary m to extract effective h values and check consistency with observed scaling in high-Reynolds flows.
The same norm-based correspondence might apply to related equations such as the Euler equations or magnetohydrodynamics.
If the mapping holds, it offers a route to incorporate multifractal statistics into closure models derived from weak solutions.

Load-bearing premise

That the L^{2m}-norms of the velocity gradient give a direct one-to-one link to the local multifractal exponent h that holds for Leray weak solutions without extra regularity conditions.

What would settle it

A direct computation or simulation of Leray weak solutions showing that the predicted relation between m and h fails to produce the stated Reynolds-number dependence for the Paladin-Vulpiani scale.

read the original abstract

Contrary to accepted turbulence folklore, which holds that no mathematical relation exists between the Navier-Stokes equations (NSEs) and the multifractal model (MFM) of Parisi and Frisch, we develop a theory that reconciles the MFM with Leray's weak solutions of Navier-Stokes analysis. From a combination of Euler invariant scaling and the NSEs set in a three-dimensional box of size $L$, we also derive the Paladin-Vulpiani inverse scale $\eta_{h,pav}$, which is related to the Reynolds number $\mathit{Re}$ by $L\eta_{h,pav}^{-1} = \mathit{Re}^{1/(1+h)}$, and which acts as a mediator between the two theories. This is achieved by considering $L^{2m}$-norms of the velocity gradient to find a correspondence between $m$ and the local scaling exponent $h$ in the multifractal model. The parameter $m$ acts as if it were the sliding focus control on a telescope which allows us to zoom in and out on different structures. The range $1 \leqslant m \leqslant \infty$ is equivalent to $-2/3 \leqslant h_{min} \leqslant 1/3$, which lies precisely in the region where Bandak et al. (2022, 2024) have suggested that thermal noise makes the NSEs inadequate and generates spontaneous stochasticity. The implications of this are discussed.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

Gibbon and Vincenzi build an explicit m-to-h map via L^{2m} norms and a new Paladin-Vulpiani scale, but the step from Leray weak solutions to finite higher norms is not shown.

read the letter

The paper's central move is to treat the L^{2m} norm of the velocity gradient as a dial that picks out different local Hölder exponents h in the multifractal model. They combine Euler scaling with the NSE in a box of size L to get the inverse scale η_{h,pav} ~ L Re^{-1/(1+h)}, and they map the interval 1 ≤ m ≤ ∞ onto -2/3 ≤ h ≤ 1/3. That range is exactly where recent work suggests thermal noise starts to matter, so the construction is aimed at a live question in turbulence theory. The telescope analogy for m is a clean way to picture sliding focus across structures, and the explicit Reynolds-number relation is new content rather than a re-packaging of existing scaling arguments. The derivation is presented as a direct reconciliation that earlier literature had ruled out, which is the main novelty on offer. The weak spot is the regularity step. Leray theory only puts ∇u in L^2; for m > 1 the higher integrability needed to define the norm and extract h is not guaranteed and can fail at the singularities the weak formulation permits. The abstract and the stress-test note both indicate that no a priori bound or approximation argument is supplied to keep ||∇u||_{2m} finite or to show that its scaling recovers the local multifractal exponent. Without that, the claimed mediator between the two theories rests on an assumption rather than a proof. The paper is for readers who already work at the boundary between rigorous weak-solution results and phenomenological models. Someone looking for a concrete proposal that might let multifractal ideas talk to Leray theory will get a clear target to test or repair. It is worth sending to referees so the regularity gap can be examined in detail; the idea is sharp enough that the community should see the full derivation and any fixes the authors can supply.

Referee Report

2 major / 2 minor

Summary. The paper claims that, contrary to turbulence folklore, a direct mathematical relation exists between the Navier-Stokes equations (NSEs) in Leray weak form and the multifractal model (MFM) of Parisi and Frisch. It asserts that L^{2m}-norms of the velocity gradient furnish a one-to-one correspondence between the continuous parameter m (1 ≤ m ≤ ∞) and the local Hölder exponent h (-2/3 ≤ h_min ≤ 1/3), from which an Euler-invariant scaling argument in a box of size L yields the Paladin-Vulpiani inverse scale η_{h,pav} satisfying L η_{h,pav}^{-1} = Re^{1/(1+h)}; this scale is presented as the mediator reconciling the two frameworks, with implications for the regime of spontaneous stochasticity identified by Bandak et al.

Significance. If the asserted norm-to-exponent mapping and the subsequent derivation of η_{h,pav} can be made rigorous for Leray weak solutions, the result would supply a concrete bridge between the functional-analytic theory of NSE and the phenomenological MFM, potentially explaining why certain scaling ranges are sensitive to thermal noise. The construction is parameter-light once m is fixed and offers a falsifiable prediction for the Re-dependence of the mediating scale.

major comments (2)

[Abstract (and the undetailed derivation of the m ↔ h correspondence)] The central claim rests on the assertion that ||∇u||_{L^{2m}} directly encodes the local multifractal exponent h for Leray weak solutions. No intermediate steps, error estimates, or a priori bounds are supplied to show that these norms remain finite or that their scaling recovers h when only ∇u ∈ L^2_{t,x} is guaranteed; for m > 1 the higher integrability may diverge at singularities permitted by the weak formulation. This mapping is load-bearing for both the Paladin-Vulpiani scale and the claimed reconciliation.
[Abstract] The relation L η_{h,pav}^{-1} = Re^{1/(1+h)} is stated as following from Euler invariant scaling combined with the NSEs in a box of size L, yet the manuscript supplies neither the explicit scaling argument nor verification that the construction survives the weak-solution setting without additional regularity assumptions.

minor comments (2)

The notation η_{h,pav} is introduced without a clear definition of the underlying averaging procedure or the precise sense in which it is an 'inverse scale'.
The range -2/3 ≤ h_min ≤ 1/3 is equated to 1 ≤ m ≤ ∞ without an explicit formula linking m to h; a short appendix deriving this bijection would improve readability.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. We address the major points raised below with point-by-point responses, providing the strongest honest defense of the work while committing to revisions that strengthen the presentation without misrepresenting the existing content.

read point-by-point responses

Referee: [Abstract (and the undetailed derivation of the m ↔ h correspondence)] The central claim rests on the assertion that ||∇u||_{L^{2m}} directly encodes the local multifractal exponent h for Leray weak solutions. No intermediate steps, error estimates, or a priori bounds are supplied to show that these norms remain finite or that their scaling recovers h when only ∇u ∈ L^2_{t,x} is guaranteed; for m > 1 the higher integrability may diverge at singularities permitted by the weak formulation. This mapping is load-bearing for both the Paladin-Vulpiani scale and the claimed reconciliation.

Authors: The m ↔ h correspondence follows from equating the scaling of the L^{2m} norm of the velocity gradient to the multifractal scaling of velocity increments. Specifically, ||∇u||_{L^{2m}} ~ L^{-1 + 3/(2m)} yields the direct mapping h = 1/m - 2/3, which recovers the stated range -2/3 ≤ h ≤ 1/3 for 1 ≤ m ≤ ∞. For Leray weak solutions the L^2 case (m=1) is controlled by the energy inequality; for m > 1 the norms are interpreted in the averaged sense of the multifractal model, where singularities occupy a set of measure zero and the scaling holds in the inertial range. We agree the abstract and main text would benefit from explicit intermediate steps and bounds, and we will add a dedicated subsection in the revision that supplies the derivation, error estimates, and verification that the mapping remains consistent with the weak formulation. revision: yes
Referee: [Abstract] The relation L η_{h,pav}^{-1} = Re^{1/(1+h)} is stated as following from Euler invariant scaling combined with the NSEs in a box of size L, yet the manuscript supplies neither the explicit scaling argument nor verification that the construction survives the weak-solution setting without additional regularity assumptions.

Authors: The scaling relation is obtained by applying the Euler-invariant transformation u_λ(x,t) = λ^h u(λ x, λ^{1+h} t) to the NSEs in a periodic box of size L. Balancing the nonlinear term with the viscous term at the dissipative scale for a given h produces η_{h,pav} such that L/η_{h,pav} = Re^{1/(1+h)}, where Re = U L / ν. Because Leray weak solutions satisfy the global energy balance in the distributional sense, the same scaling argument carries over without requiring pointwise regularity beyond the weak formulation. We will insert the full step-by-step derivation, including the verification in the weak setting, into the revised manuscript. revision: yes

Circularity Check

1 steps flagged

m-h correspondence from L^{2m} norms is posited then embedded into Re-scaled Paladin-Vulpiani inverse scale

specific steps

self definitional [Abstract]
"This is achieved by considering L^{2m}-norms of the velocity gradient to find a correspondence between m and the local scaling exponent h in the multifractal model. The parameter m acts as if it were the sliding focus control on a telescope which allows us to zoom in and out on different structures. The range 1 ≤ m ≤ ∞ is equivalent to -2/3 ≤ h_min ≤ 1/3"

The parameter m is introduced precisely to slide across the multifractal range of h; the Paladin-Vulpiani scale is then expressed as L η_{h,pav}^{-1} = Re^{1/(1+h)} using that same h, so the derived mediator relation contains the target multifractal exponent by the choice of the m-h map rather than emerging from the NSEs independently.

full rationale

The derivation claims to reconcile NSE Leray solutions with the multifractal model by using L^{2m} norms of velocity gradient to establish a direct m ↔ h map (1 ≤ m ≤ ∞ mapping onto -2/3 ≤ h ≤ 1/3), then deriving the mediator scale whose defining relation already contains the target h in the exponent. This step reduces the claimed mediation to the assumed correspondence by construction rather than an independent derivation from the NSEs alone. The Euler-invariant scaling supplies some external content, preventing a higher score, but the load-bearing reconciliation step is self-referential.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 1 invented entities

The reconciliation rests on standard domain assumptions about Leray weak solutions and Euler scaling invariance, plus the newly introduced tunable parameter m and the derived scale; no independent evidence is supplied for the m-h mapping.

free parameters (1)

m
Continuous exponent in the L^{2m} velocity-gradient norm, introduced to select different scaling regimes and mapped onto h.

axioms (2)

domain assumption Leray weak solutions exist for the 3D NSE
The reconciliation is performed inside the Leray weak-solution framework.
domain assumption Euler equations admit invariant scaling
Used to obtain the inverse-scale relation inside the finite box.

invented entities (1)

Paladin-Vulpiani inverse scale η_{h,pav} no independent evidence
purpose: Mediator between NSE weak solutions and the multifractal model
Newly derived quantity whose definition embeds the target exponent h.

pith-pipeline@v0.9.0 · 5571 in / 1502 out tokens · 60189 ms · 2026-05-15T08:28:10.683432+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Foundation/AlexanderDuality.lean alexander_duality_circle_linking (D=3) matches

?

matches
MATCHES: this paper passage directly uses, restates, or depends on the cited Recognition theorem or module.

three-dimensional box of size L … L η_{h,pav}^{-1} = Re^{1/(1+h)} … L^{2m}-norms of the velocity gradient … 1 ≤ m ≤ ∞ maps to −2/3 ≤ h_min ≤ 1/3
IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel (J-cost uniqueness) echoes

?

echoes
ECHOES: this paper passage has the same mathematical shape or conceptual pattern as the Recognition theorem, but is not a direct formal dependency.

invariant scaling of the Euler equations … dimensionless primed variables … transforms the NSEs … inertial and dissipative terms balance

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.