arxiv: 2603.28801 · v2 · submitted 2026-03-26 · ⚛️ physics.flu-dyn · nlin.CD· physics.comp-ph

Recognition: 2 theorem links

· Lean Theorem

Sparse M\"untz--Sz\'asz Recovery for Boundary-Anchored Velocity Profiles: A Short-Record Roughness Diagnostic in Turbulence

D Yang Eng

Authors on Pith no claims yet

Pith reviewed 2026-05-15 00:23 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn nlin.CDphysics.comp-ph

keywords sparse recoveryMuntz-Szasz dictionaryturbulence roughnessscaling exponentsvelocity incrementsvorticity alignmentshort-record diagnostics

0 comments

The pith

Sparse regression recovers local scaling exponents from short boundary-anchored velocity profiles in turbulence as a roughness diagnostic.

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

The paper introduces a sparse convex-relaxation method that estimates effective local scaling exponents from short velocity-increment profiles anchored at boundaries with only about 40 points. It solves an l1-regularized regression inside a mixed Muntz-Szasz and Jacobi dictionary and treats the output as a finite-scale directional roughness diagnostic rather than a pointwise exponent. On isotropic data from the Johns Hopkins Turbulence Database the approach shows strong internal self-consistency, with F1 scores near 0.93 on subsampled benchmarks and balanced accuracy of 0.928 on synthetic controls. The recovered exponents remain only weakly tied to dissipation but consistently smaller in regions of higher vorticity, and directional tests reveal a positive contrast aligned with the vorticity axis. The method is positioned as a geometric complement to energetic observables that can resolve low-order anisotropic organization without long records.

Core claim

A sparse l1-regularized regression in a mixed Muntz-Szasz/Jacobi dictionary applied to short boundary-anchored velocity-increment profiles yields effective local scaling exponents that serve as a finite-scale directional roughness diagnostic. On isotropic turbulence data this detector achieves F1 approximately 0.93 under internal subsampling against N=200 labels and balanced accuracy 0.928 at N=40 on synthetic controls. Across Reynolds numbers 433 to 1300 the fixed-window sharp fraction stays between 30 and 50 percent, higher vorticity correlates with smaller detected exponents, and vorticity-aligned directional contrast Pi_alpha averages 0.093 with a statistically detectable quadrupolar leg

What carries the argument

l1-regularized regression inside a mixed Muntz-Szasz/Jacobi dictionary applied to boundary-anchored profiles, functioning as a finite-scale directional roughness diagnostic

If this is right

The fixed-window fraction of sharp profiles remains between 30 and 50 percent across the tested Reynolds-number range.
Higher vorticity is associated with systematically smaller recovered roughness exponents in conditioned samples.
A positive vorticity-aligned contrast Pi_alpha appears on true axes but not on fake axes, together with a low-order quadrupolar component.
Positive Pi_alpha persists at both the smallest and largest tested radii in a seeded scale-transfer scan.

Where Pith is reading between the lines

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

The same short-record detector could be applied to non-isotropic or wall-bounded flows to test whether the vorticity-roughness link survives changes in geometry.
The observed low-order quadrupolar pattern suggests the diagnostic may capture early signatures of anisotropic organization that energetic spectra miss.
Because the method works with N approximately 40, it opens the possibility of real-time roughness tracking in experimental time series where full inertial-range records are unavailable.

Load-bearing premise

Internal subsampling benchmarks and synthetic controls are taken to validate short-record performance and directional contrast without external calibration or independent checks on whether the dictionary isolates true scaling exponents.

What would settle it

An independent long-record calculation of local Holder exponents on the same high-vorticity centers that yields a directional contrast Pi_alpha statistically indistinguishable from zero would falsify the claim that the detector isolates genuine roughness structure.

Figures

Figures reproduced from arXiv: 2603.28801 by D Yang Eng.

**Figure 3.** Figure 3: FIG. 3. Subsampling study on JHTDB profiles. The panel [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗

**Figure 4.** Figure 4: FIG. 4. Baseline comparison at [PITH_FULL_IMAGE:figures/full_fig_p008_4.png] view at source ↗

**Figure 5.** Figure 5: FIG. 5. Targeted control experiments: (a) magnitude vs. [PITH_FULL_IMAGE:figures/full_fig_p009_5.png] view at source ↗

**Figure 6.** Figure 6: FIG. 6. Dissipation–singularity correlation ( [PITH_FULL_IMAGE:figures/full_fig_p010_6.png] view at source ↗

**Figure 7.** Figure 7: FIG. 7. Lagrangian singularity evolution ( [PITH_FULL_IMAGE:figures/full_fig_p011_7.png] view at source ↗

**Figure 8.** Figure 8: FIG. 8. Reynolds-number sweep at fixed physical [PITH_FULL_IMAGE:figures/full_fig_p012_8.png] view at source ↗

read the original abstract

We present a sparse convex-relaxation framework for estimating effective local scaling exponents from short boundary-anchored velocity-increment profiles ($N\approx40$). The detector solves an $\ell_1$-regularized regression in a mixed M\"untz--Sz\'asz/Jacobi dictionary and is interpreted throughout as a finite-scale, directional roughness diagnostic rather than a pointwise H\"older exponent. On isotropic datasets from the Johns Hopkins Turbulence Database, an internal subsampling benchmark against $N=200$ detector labels gives $F_1\approx0.93$ across nine unweighted reruns, and a balanced synthetic control gives balanced accuracy $0.928$ at $N=40$, indicating useful short-record self-consistency without constituting an external calibration. Across $Re_\lambda\approx433$--$1300$, the fixed-window sharp fraction remains of order $30$--$50\%$, but a scale-normalized control does not isolate a clean Reynolds-number trend. The recovered $\hat{\alpha}$ is only weakly associated with dissipation, whereas higher vorticity is consistently associated with smaller detected roughness exponents in conditioned samples. Directional controls on 60 high-vorticity centers further show a positive vorticity-aligned contrast $\Pi_\alpha$ (mean $0.093$, bootstrap 95\% CI $[0.028,0.158]$), stronger on the true vorticity axis than on fake axes, together with a statistically detectable low-order quadrupolar component in a joint Legendre fit. A seeded scale-transfer scan shows positive $\Pi_\alpha$ at both the smallest and largest tested radii, supporting finite-range persistence without a strong theorem-level nonlocal claim. The method is therefore best viewed as a finite-scale geometric diagnostic complementary to energetic observables, capable of resolving directional structure and low-order anisotropic organization in high-vorticity regions.

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

The paper introduces a sparse l1 recovery in a mixed Müntz-Szász/Jacobi dictionary for short boundary-anchored velocity profiles as a finite-scale roughness diagnostic, with solid internal F1 and accuracy numbers but only self-consistency checks.

read the letter

The paper's main contribution is a sparse convex-relaxation approach that fits a mixed Müntz-Szász and Jacobi dictionary to short boundary-anchored velocity-increment profiles of length around 40 points. It treats the recovered exponents as a directional roughness diagnostic for turbulence rather than literal Hölder values. This setup is new in applying those basis functions to this specific turbulence context, and the results include an F1 score of roughly 0.93 on subsampled Johns Hopkins isotropic data against longer labels, plus 0.928 balanced accuracy on synthetic controls. They also find a vorticity-aligned positive contrast of 0.093 with a bootstrap interval from 0.028 to 0.158, plus some quadrupolar structure, all while staying careful about finite-scale limits and self-consistency. The main limitation is that the benchmarks are all internal, using subsampling and synthetics without testing on independent datasets or against direct roughness measurements. The l1 regularization strength is a free parameter, and without full details on how it was chosen or if sensitivity checks were done, that could affect reproducibility. The Reynolds number dependence doesn't yield a clean trend, which they note but don't push further. This kind of work would appeal to experimental fluid dynamicists who often work with limited time series and need tools to extract geometric features from high-vorticity regions. It doesn't tackle fundamental turbulence problems but provides a complementary diagnostic. I would send it for peer review. The internal numbers are presented with intervals and the framing is proportionate, so referees could usefully check the implementation and suggest external tests.

Referee Report

1 major / 2 minor

Summary. The paper introduces a sparse convex-relaxation method using ℓ1-regularized regression over a mixed Müntz–Szász/Jacobi dictionary to recover effective local scaling exponents from short (N≈40) boundary-anchored velocity-increment profiles. It positions the detector as a finite-scale directional roughness diagnostic and reports internal self-consistency benchmarks (F1≈0.93 on Johns Hopkins subsamples, balanced accuracy 0.928 on synthetic controls) together with a statistically detectable positive vorticity-aligned contrast Π_α (mean 0.093, 95% CI [0.028,0.158]) in high-vorticity regions.

Significance. If the internal benchmarks hold, the work supplies a practical geometric complement to energetic diagnostics that can resolve low-order directional structure and quadrupolar anisotropy inside intense vorticity regions at finite scales. The explicit disclaimer of external calibration and the demonstration of persistence across tested radii strengthen its utility as a short-record tool.

major comments (1)

[Abstract] Abstract and validation description: the reported F1 and accuracy figures rest on subsampling and synthetic controls, yet the manuscript supplies no explicit protocol for data exclusion, hyperparameter selection for the ℓ1 strength, or checks against post-hoc tuning; these choices are load-bearing for the self-consistency claim.

minor comments (2)

[Methods] Notation: the mixed dictionary is introduced without an explicit basis-function list or orthogonality statement; adding a short table of the retained Müntz–Szász and Jacobi terms would improve reproducibility.
[Results] Figure clarity: the directional contrast plots would benefit from explicit indication of the fake-axis controls on the same panel to allow immediate visual comparison of the reported Π_α difference.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the positive assessment and the recommendation for minor revision. We address the single major comment below by committing to explicit documentation of the validation protocols.

read point-by-point responses

Referee: [Abstract] Abstract and validation description: the reported F1 and accuracy figures rest on subsampling and synthetic controls, yet the manuscript supplies no explicit protocol for data exclusion, hyperparameter selection for the ℓ1 strength, or checks against post-hoc tuning; these choices are load-bearing for the self-consistency claim.

Authors: We agree that the validation protocols require more explicit documentation. In the revised manuscript we will add a dedicated subsection (Methods, new §3.4) that specifies: (i) the exact data-exclusion criteria applied to the Johns Hopkins subsamples (velocity-increment magnitude thresholds, boundary-anchoring tolerance, and any outlier rejection rules); (ii) the procedure used to select the ℓ1 regularization strength, including the grid range, selection criterion (e.g., cross-validation on held-out synthetic profiles or information criterion), and the final value retained; and (iii) additional robustness checks consisting of a sensitivity sweep over a factor-of-ten range in the regularization parameter together with the resulting variation in F1 and balanced accuracy. These additions will be referenced concisely in the abstract. The core numerical results remain unchanged. revision: yes

Circularity Check

0 steps flagged

No significant circularity identified

full rationale

The paper's central results derive from applying an ℓ1-regularized regression in a Müntz–Szász/Jacobi dictionary to short velocity-increment profiles drawn from the external Johns Hopkins Turbulence Database. Reported metrics (F1≈0.93 from N=200 subsampling, balanced accuracy 0.928 on synthetic controls, and directional contrast Π_α with bootstrap CI) are obtained via independent data splits and controls rather than by reducing any prediction to a fitted parameter or self-referential equation. No self-definitional steps, fitted-input predictions, load-bearing self-citations, or ansatz smuggling appear; the work explicitly disclaims external calibration and positions outputs as finite-scale self-consistency diagnostics.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The central claim rests on the assumption that velocity increments admit a useful sparse representation in the chosen dictionary and that internal validation metrics reflect true diagnostic utility; limited abstract information prevents exhaustive enumeration of parameters.

free parameters (1)

l1 regularization strength
Controls sparsity in the regression; value not specified in abstract but required for the convex relaxation to function.

axioms (1)

domain assumption Velocity-increment profiles admit a sparse representation in the mixed Muntz-Szasz/Jacobi dictionary that isolates effective scaling exponents.
Invoked to justify the regression recovering meaningful roughness measures from short records.

pith-pipeline@v0.9.0 · 5653 in / 1382 out tokens · 50136 ms · 2026-05-15T00:23:14.110259+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/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

sparse convex-relaxation framework ... mixed Müntz–Szász/Jacobi dictionary ... finite-scale, directional roughness diagnostic
IndisputableMonolith/Foundation/AlexanderDuality.lean alexander_duality_circle_linking unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

Singular Separation Condition (SSC) and population gap κ_pop(α)

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.

Reference graph

Works this paper leans on

26 extracted references · 26 canonical work pages

[1]

M¨ untz–Sz´ asz Geometric Separation and Computational Complexity of Sparse Recovery in Ultra- Coherent Dictionaries,

D Y. Eng, “M¨ untz–Sz´ asz Geometric Separation and Computational Complexity of Sparse Recovery in Ultra- Coherent Dictionaries,” TechRxiv doi:https://doi.or g/10.36227/techrxiv.177223034.43626170/v1(2026)

work page doi:10.36227/techrxiv.177223034.43626170/v1(2026 2026
[2]

Geometric Dynamics of Turbulence: A Minimal Oscillator Structure from Non-local Closure,

A. Sevilla, “Geometric Dynamics of Turbulence: A Minimal Oscillator Structure from Non-local Closure,” arXiv:2603.18913 (2026)

work page arXiv 2026
[3]

A public turbulence database cluster and applications to study Lagrangian evolution of velocity increments in turbulence,

Y. Li, E. Perlman, M. Wan,et al., “A public turbulence database cluster and applications to study Lagrangian evolution of velocity increments in turbulence,”J. Tur- bulence9, N31 (2008)

work page 2008
[4]

On the singularity structure of fully developed turbulence,

G. Parisi and U. Frisch, “On the singularity structure of fully developed turbulence,” inTurbulence and Pre- dictability in Geophysical Fluid Dynamics and Climate Dynamics(Proceedings of the International School of Physics “Enrico Fermi”, Course LXXXVIII, Varenna, Italy, 1983), edited by M. Ghil, R. Benzi, and G. Parisi, pp. 84–87, North-Holland, Amsterdam (1985)

work page 1983
[5]

Frisch,Turbulence: The Legacy of A

U. Frisch,Turbulence: The Legacy of A. N. Kolmogorov (Cambridge University Press, Cambridge, 1995)

work page 1995
[6]

The phenomenol- ogy of small-scale turbulence,

K. R. Sreenivasan and R. A. Antonia, “The phenomenol- ogy of small-scale turbulence,”Annu. Rev. Fluid Mech. 29, 435–472 (1997)

work page 1997
[7]

Compressed sensing,

D. L. Donoho, “Compressed sensing,”IEEE Trans. Inf. Theory52, 1289–1306 (2006)

work page 2006
[8]

SciPy 1.0: Fundamental Algorithms for Scientific Computing in Python,

P. Virtanenet al., “SciPy 1.0: Fundamental Algorithms for Scientific Computing in Python,”Nature Methods17, 261–272 (2020)

work page 2020
[9]

A refinement of previous hypotheses concerning the local structure of turbulence in a viscous incompressible fluid at high Reynolds number,

A. N. Kolmogorov, “A refinement of previous hypotheses concerning the local structure of turbulence in a viscous incompressible fluid at high Reynolds number,”J. Fluid Mech.13, 82–85 (1962)

work page 1962
[10]

Universal scaling laws in fully developed turbulence,

Z.-S. She and E. L´ evˆ eque, “Universal scaling laws in fully developed turbulence,”Phys. Rev. Lett.72, 336– 339 (1994)

work page 1994
[11]

The multifractal nature of turbulent energy dissipation,

C. Meneveau and K. R. Sreenivasan, “The multifractal nature of turbulent energy dissipation,”J. Fluid Mech. 224, 429–484 (1991)

work page 1991
[12]

Wavelets and multifractal formalism for singular signals: Application to turbulence data,

J. F. Muzy, E. Bacry, and A. Arneodo, “Wavelets and multifractal formalism for singular signals: Application to turbulence data,”Phys. Rev. Lett.67, 3515–3518 15 (1991)

work page 1991
[13]

Wavelet based multifractal formalism: appli- cations to DNA sequences, satellite images of the cloud structure, and stock market data,

A. Arneodo, B. Audit, N. Decoster, J.-F. Muzy, and C. Vaillant, “Wavelet based multifractal formalism: appli- cations to DNA sequences, satellite images of the cloud structure, and stock market data,” inThe Science of Dis- asters(Springer, 2002), pp. 26–102

work page 2002
[14]

p-exponent andp-leaders, Part I: Negative pointwise regularity,

S. Jaffard, C. Melot, R. Leonarduzzi, H. Wendt, P. Abry, S. G. Roux, and M. E. Torres, “p-exponent andp-leaders, Part I: Negative pointwise regularity,”Physica A448, 300–318 (2016)

work page 2016
[15]

Examination of hypotheses in the Kolmogorov refined turbulence theory through high resolution simu- lations. Part 1. Velocity field,

L.-P. Wang, S. Chen, J. G. Brasseur, and J. C. Wyn- gaard, “Examination of hypotheses in the Kolmogorov refined turbulence theory through high resolution simu- lations. Part 1. Velocity field,”J. Fluid Mech.309, 113– 156 (1996)

work page 1996
[16]

Study of high– Reynolds number isotropic turbulence by direct numer- ical simulation,

T. Ishihara, T. Gotoh, and Y. Kaneda, “Study of high– Reynolds number isotropic turbulence by direct numer- ical simulation,”Annu. Rev. Fluid Mech.41, 165–180 (2009)

work page 2009
[17]

Energy dissipation rate and energy spectrum in high resolution direct numerical simulations of tur- bulence in a periodic box,

Y. Kaneda, T. Ishihara, M. Yokokawa, K. Itakura, and A. Uno, “Energy dissipation rate and energy spectrum in high resolution direct numerical simulations of tur- bulence in a periodic box,”Phys. Fluids15, L21–L24 (2003)

work page 2003
[18]

A probabilistic and RIPless theory of compressed sensing,

E. J. Cand` es and Y. Plan, “A probabilistic and RIPless theory of compressed sensing,”IEEE Trans. Inf. Theory 57, 7235–7254 (2011)

work page 2011
[19]

Simultane- ous analysis of Lasso and Dantzig selector,

P. J. Bickel, Y. Ritov, and A. B. Tsybakov, “Simultane- ous analysis of Lasso and Dantzig selector,”Ann. Stat. 37, 1705–1732 (2009)

work page 2009
[20]

Fully developed turbulence and the multifractal conjecture,

R. Benzi and L. Biferale, “Fully developed turbulence and the multifractal conjecture,”J. Stat. Phys.135, 977– 990 (2009)

work page 2009
[21]

Extended self-similarity in turbulent flows,

R. Benzi, S. Ciliberto, R. Tripiccione, C. Baudet, F. Mas- saioli, and S. Succi, “Extended self-similarity in turbulent flows,”Phys. Rev. E48, R29–R32 (1993)

work page 1993
[22]

Beyond Kolmogorov cascades,

B. Dubrulle, “Beyond Kolmogorov cascades,”J. Fluid Mech.867, P1 (2019)

work page 2019
[23]

Multi-fractality, universality and singular- ity in turbulence,

B. Dubrulle, “Multi-fractality, universality and singular- ity in turbulence,”Fractal Fract.6, 613 (2022)

work page 2022
[24]

Local estimates of H¨ older exponents in turbulent vector fields,

F. H. Nguyen, J.-P. Laval, P. Kestener, A. Cheskidov, R. Shvydkoy, and B. Dubrulle, “Local estimates of H¨ older exponents in turbulent vector fields,”Phys. Rev. E99, 053114 (2019)

work page 2019
[25]

A proof of Onsager’s conjecture,

P. Isett, “A proof of Onsager’s conjecture,”Ann. Math. 188, 871–963 (2018)

work page 2018
[26]

Potentially singular solutions of the 3D axisymmetric Euler equations,

G. Luo and T. Y. Hou, “Potentially singular solutions of the 3D axisymmetric Euler equations,”Proc. Natl. Acad. Sci. U.S.A.111, 12968–12973 (2014)

work page 2014