REVIEW 3 major objections 5 minor 1 cited by
Circulation Fluctuations of Elementary Turbulent Vortices
T0 review · 3 major / 5 minor · reviewed 2026-08-05 · deepseek-v4-flash
Pith's one-line read Fast dissipation modes explain why vortex circulations are fat-tailed.
desk verdict The paper's real contribution is a closed-form fat-tailed PDF for elementary vortex circulation from a GMC split of the dissipation field; it deserves review, but the headline 'no-fit' agreement rests on the untested ℓ≈λ identification. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the decomposition $\phi = \phi_< + \phi_>$ of the lognormal GMC scalar field into slow (low-wavenumber) and fast (high-wavenumber) modes, with the crossover $\ell$ taken to be the Taylor microscale. The key identity is $\bar{\Gamma} = \xi_>^2 \tilde{\Gamma}$ (Eq. 3.10 with $\beta=2$), which makes the redefined circulation a product of a lognormal fast-mode variable and a Gaussian variable. The $\beta=2$ choice follows from requiring the conditioned vortex density $\sigma \eta^2$ to equal the GMC density $\xi$ in law, using the $Z_2$ symmetry of $\phi_>$. This product structure carries the argument: the quadratic log-moments (3.25), the parameter-free closed-form PDF (3.32), and the Gauss
What would settle it
Compute the $q^2$ coefficient in $\ln \left( \frac{\langle |\Gamma|^q \rangle}{A_q} \right)$ for moments beyond $q=6$ or at a Reynolds number beyond 2556 and compare it with $\frac{\mu}{2} \ln\left(\frac{\lambda}{\eta}\right)$. The paper's prediction fixes that coefficient exactly; any statistically significant deviation would falsify the $\beta=2$ product structure and the $\ell \approx \lambda$ identification. A second check is to measure the two-point correlation of standardized elementary circulations at separations below the Taylor microscale, where the fast-mode decomposition predicts a decay governed by the $(r/\ell)^{-3/2}$ behavior of the $\phi_>$ correlator rather than a flat G
Extended reading notes
Core claim
The paper's central claim is that elementary vortex circulation should be redefined as $\bar{\Gamma}(x) = \xi_>^2(x) \tilde{\Gamma}(x)$, where $\tilde{\Gamma}$ is the original Gaussian circulation field and $\xi_>$ is the fast-mode factor of the vortex density obtained by splitting the GMC lognormal field at a crossover scale $\ell \approx \lambda$ (the Taylor microscale). With the exponent $\beta = 2$, the conditioned vortex number density equals, in law, the full GMC density, so the original vortex gas model is recovered unchanged at inertial scales. At a single point the new circulation is a lognormal variable times a Gaussian, hence non-Gaussian and fat-tailed, while for separations larger than \e
Load-bearing premise
The prediction stands or falls on the assumption that the crossover between slow and fast dissipation modes sits at the Taylor microscale ($\ell \approx \lambda$); if the true crossover sits elsewhere, the constant $c$ and the whole predicted PDF shift.
Editorial extensions
If this is right
- The closed-form PDF (3.32) becomes a parameter-free prediction once \ell = \lambda; the paper's DNS comparison shows it reproduces the measured distributions for five datasets spanning R_\lambda = 433 to 2556.
- The vortex circulation scale \tilde{\Gamma}_0 scales as u'L / R_\lambda^2, i.e. proportionally to the kinematic viscosity, so the characteristic circulation of elementary vortices decreases sharply with Reynolds number.
- The vortex number density inherits GMC scaling: the density autocorrelation decays as r^{-\mu/4} and coarse-grained density moments follow exponents (\mu/8)q(1-q), providing alternative measurements of the intermittency exponent.
- The original inertial-range predictions of the vortex gas model are preserved; the fast-mode absorption modifies only the sub-Taylor-scale single-point statistics, keeping the model's earlier successes intact.
- For moment orders beyond about q=6 the free-field approximation used in the construction is expected to break down, so the model identifies where deviations from the predicted moment parabola should appear.
Reading between the lines
- Extending the same construction to two-point statistics would give a concrete prediction for how elementary circulations become correlated below the Taylor scale: the fast-mode factor contributes a decay set by the (r/\ell)^{-3/2} behavior of the \phi_> correlator, which could be measured directly in DNS.
- If the crossover scale is only approximately the Taylor microscale, the same moment relation provides a way to measure the effective \ell from the q^2 coefficient, allowing a test of whether \ell = \lambda survives at Reynolds numbers beyond 2556.
- A coarse-grained or large-eddy formulation of the model should see a renormalization between vortex density and circulation fluctuations as the filter scale changes; the \beta=2 absorption suggests that the effective elementary circulation depends on resolution even though total circulation statistics remain fixed.
- The same algebraic structure may apply to other intermittent small-scale turbulent observables: any field built as a lognormal multiplier times a Gaussian component will show fat tails at one point while factorizing at large separations.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper extends the vortex gas model (VGM) of homogeneous isotropic turbulence by replacing the Gaussian elementary circulation field \tilde{\Gamma}(x) with \bar{\Gamma}(x)=\xi_>^2(x)\tilde{\Gamma}(x), where \xi_> is the fast-mode GMC factor obtained by splitting the dissipation scalar field at a crossover scale \ell. A symmetry argument and the empirical law (3.15) select \beta=2. The resulting single-point circulation PDF, Eq. (3.32), is a lognormal-Gaussian product with variance fixed by c=(\ell/\eta)^{\mu/2}, while multi-point correlations factorize at separations large compared with \ell. The model is validated against JHU DNS data at R_\lambda=433, 610, 1278, and 2556, using moments (3.25), the PDF (3.32) with no fitted parameters, density correlations (3.33), and coarse-grained density scalings (3.35)-(3.38). The key physical input is the working hypothesis \ell\approx\lambda (Taylor microscale), stated in Sec. III(v).
Significance. If correct, this is a conceptually important result: it resolves an apparent paradox in the VGM—non-Gaussian circulation of individual vortices coexisting with Gaussian-like multi-point correlations—by tracing the non-Gaussianity to the GMC structure of the dissipation field. The paper also provides a closed-form, parameter-free prediction for the elementary circulation PDF that is tested extensively at four Reynolds numbers using public DNS data. The derivation is internally consistent given the stated assumptions, and the numerical validations cover a broad parameter range. The main limitation is that the headline 'no fitting' PDF depends on the unverified identification \ell\approx\lambda; this is a load-bearing working hypothesis rather than an independently measured or tested scale.
major comments (3)
- [Sec. III(v), Eqs. (3.20) and (3.31)] The entire shape of the predicted PDF (3.32) is fixed by the constant c=(\ell/\eta)^{\mu/2}, and the paper sets \ell\approx\lambda based only on a citation to Ishihara et al. [35] that the Taylor microscale is the typical thickness of dissipation layers. This is explicitly labeled a 'working hypothesis.' No independent measurement or sensitivity analysis is provided. Since \sigma_z^2=\mu\ln(\lambda/\eta) determines the fatness of the predicted distribution, a modest error in \ell/\lambda changes the PDF in a way that could degrade the apparent agreement in Fig. 3. This is not a mere technicality: it is the load-bearing scale choice for the central parameter-free claim.
- [Sec. IV, Fig. 2(a) and Eq. (3.25)] The validation of the moment relation (3.25) fixes the quadratic coefficient a priori to \ln c, i.e., to the value implied by \ell\approx\lambda. This prevents the DNS data from independently testing the most distinctive prediction of the model. The authors should report unconstrained quadratic fits of ln(<|\bar{\Gamma}|^q>/A_q) versus q for each R_\lambda, with confidence intervals for the q^2 coefficient, and compare these with \mu/2 \ln(\lambda/\eta). Such a test would either confirm the working hypothesis or reveal its limitations. Without it, the 'no fitting' agreement in Fig. 3 remains conditional on an untested input.
- [Sec. III, Eqs. (3.15)-(3.16)] The selection \beta=2 uses the empirical equality \sigma(x|\epsilon)\eta^2 \stackrel{d}{=} \xi(x), which is taken from the authors' earlier DNS study [20]. The present paper then validates the model on DNS data that include the same type of data (and very likely overlapping R_\lambda=433 dataset). This is not a fit of Eq. (3.32), but it is a model-selection step informed by the validation dataset. The authors should clarify the degree of overlap and, if possible, provide an independent check of \beta=2—for example, comparing measured conditional densities with \xi_0\xi_<\xi_>^{1-\beta} for \beta=2 against \beta=0 or other values.
minor comments (5)
- [Eq. (3.23)] The notation 'c q2' should be typeset as c^{q^2}; similarly, the exponent in Eq. (3.26) is c^{q^2/2}. Please correct the typography to avoid confusion.
- [Sec. IV, Fig. 3] The excellent visual agreement in Fig. 3 is stated qualitatively. Please provide a quantitative measure (e.g., Kolmogorov-Smirnov statistic, mean squared log error, or chi-square per degree of freedom) and state the number of vortex samples used at each Reynolds number.
- [Sec. IV, Fig. 4] The predicted exponent -\mu/4 is about -0.0425, which is very small. The 'reasonable scaling range' over intermediate scales should be supported by a fit with confidence intervals and by stating the range of r/\eta used; otherwise the agreement is hard to assess visually.
- [Sec. IV, vortex detection] The analysis depends on the swirling-strength threshold (|Im(\lambda)|>\sigma_\lambda/8) and, for densities, on the KDE bandwidth (8\eta). Please state whether the circulation PDF and the reported scalings are robust to reasonable variations of these choices, and reference the earlier papers where the threshold was calibrated.
- [Eq. (3.22)] The asymptotic evaluation leading to the power-law tail of the \phi_> correlator is sketched; a few more steps or a reference would help the reader verify the sign and the prefactors, especially since Eq. (3.22) is used to justify the Gaussian factorization property.
Circularity Check
No significant circularity: the vortex-circulation PDF is derived and validated without fitting; the ℓ≈λ working hypothesis is an explicit, testable assumption rather than a circular input.
full rationale
The paper's central new result is Eq. (3.32), a closed-form PDF for the normalized elementary vortex circulation. The derivation is self-contained given the GMC setup: Γ⋆ = X·Y with X lognormal (σ_z² = μ ln(ℓ/η)) and Y Gaussian, producing the integral in (3.32). No parameter of the PDF is fitted to the circulation data used in Fig. 3: μ is taken from Ref. [40], ℓ is set to the Taylor microscale λ as an explicitly stated working hypothesis (Sec. III(v)), and λ/η is computed from DNS. The step selecting β=2 uses the empirical identity (3.15) from the authors' prior DNS study [20] and the Z₂ symmetry; although this is a self-citation, it is an externally checkable DNS observation about a different observable (vortex density vs. dissipation), not an unverified theorem, and the subsequent PDF validation is a genuinely new, independent test. The only notable weakness is the ℓ≈λ assumption, which fixes the q² coefficient in (3.25) and the width σ_z in (3.31); if the true GMC mode-split scale differed from λ, the predicted PDF would shift. However, choosing a scale by physical argument and then testing the resulting parameter-free curve is not circular: the scale is not inferred from the PDF data. Figures 2 and 4–6 similarly test derived scalings; where a prefactor is fitted (Fig. 5), the scaling exponent is still predicted. Hence no step reduces by construction to its inputs.
Assumptions & free parameters
free parameters (3)
- Gamma_tilde_0 (typical elementary circulation scale) =
Fitted per dataset via the linear coefficient in Eq. (3.25); values inferred from Fig. 2(a).
- KDE bandwidth =
8 eta
- Vortex detection threshold =
|Im(lambda)| > sigma_lambda / 8
assumptions (7)
- domain assumption The vortex gas model assumptions: circulation is approximated by the flux of thin vortex tubes crossing the contour, and the probability of finding a vortex tube depends on the local dissipation rate.
- domain assumption The fields xi(x) and Gamma_tilde(x) are independent.
- domain assumption Elementary circulation Gamma_tilde is a Gaussian random field with power-law spectrum |k|^{alpha-2} and alpha = 2 - mu/4 - zeta_2.
- domain assumption The free scalar field approximation for phi(x) is adequate for circulation moments up to about q = 6.
- domain assumption The empirical law sigma(x|epsilon) eta^2 has the same probability law as xi(x), taken from the authors' prior DNS study.
- ad hoc to paper The crossover scale ell between fast and slow modes is approximately the Taylor microscale lambda.
- standard math Standard Gaussian multiplicative chaos theory for log-correlated Gaussian fields, including source functional formulas.
Cite this review
Pith. "Pith review of Circulation Fluctuations of Elementary Turbulent Vortices." pith.science (2026). https://pith.science/paper/TTVRUYKC
@misc{pith2026250813958,
author = {Pith},
title = {Pith review of: Circulation Fluctuations of Elementary Turbulent Vortices},
year = {2026},
howpublished = {\url{https://pith.science/paper/TTVRUYKC}},
note = {Machine review of arXiv:2508.13958}
}
read the original abstract
Thin vortex tubes, with core sizes within the dissipation range, profuse in a homogeneous and isotropic turbulent flow. Their intersections with an arbitrary plane define, as a mathematical construct, a dilute gas of localized, intermittently distributed, two-dimensional vortex spots. While their planar density fluctuations are described by a field-theoretical extension of log-normal single-point statistics, known as Gaussian multiplicative chaos (GMC), they carry circulations which are Gaussian-correlated throughout the inertial range. It is puzzling, then, to find that the circulations of individual vortices are fat-tailed distributed, an apparent paradox that we fix within the GMC framework. The solution, validated through the examination of direct numerical simulation data for a broad range of Reynolds numbers, unveils, as a surprising phenomenological result, an existing coupling between the circulation of vortex structures and the short-distance properties of their spatial distribution fluctuations at sub-Taylor microscales.
Figures
Figures from the paper (3 more)
Forward citations
Cited by 1 Pith paper
-
Optimal Surfaces for Turbulent Circulation Statistics
A vortex gas model with log-correlated vortex densities predicts that turbulent circulation statistics are governed by the minimal surface bounded by the contour, matching DNS for nonplanar loops.
Reference graph
Works this paper leans on
- [20]
-
[35]
T. Ishihara, Y. Kaneda, and J.C.R. Hunt, Flow Turbul. Combust.91, 895 (2013)
work page 2013
- [1]
- [2]
- [3]
- [4]
- [5]
- [6]
Show all 47 references
-
[7]
L´ evy,Th´ eorie de l’Addition des Variables Al´ eatoires, Gauthier-Villars (1937)
P. L´ evy,Th´ eorie de l’Addition des Variables Al´ eatoires, Gauthier-Villars (1937)
1937
-
[8]
Mandelbrot,The Fractal Geometry of Nature, W.H
B. Mandelbrot,The Fractal Geometry of Nature, W.H. Freeman (1982)
1982
-
[9]
Wiener, Proc
N. Wiener, Proc. Natl. Acad. Sci. U.S.A.7, 253 (1921)
1921
-
[10]
Kac, Bull
M. Kac, Bull. Amer. Math. Soc.72, 52 (1966). 19
1966
-
[11]
Kolmogorov, Dokl
A.N. Kolmogorov, Dokl. Akad. Nauk SSSR30, 301 (1941); Proc. R. Soc. Lond. A434, 9 (1991)
1941
-
[12]
Obukhov, J
A.M. Obukhov, J. Fluid Mech.13, 77 (1962)
1962
-
[13]
Kolmogorov, J
A.N. Kolmogorov, J. Fluid Mech.13, 82 (1962)
1962
-
[14]
Smagorinsky, Mon
J. Smagorinsky, Mon. Weather Rev.91, 99 (1963)
1963
-
[15]
Kadanoff, Phys
L.P. Kadanoff, Phys. Phys. Fiz.2, 263 (1966)
1966
-
[16]
Sreenivasan and J
K.R. Sreenivasan and J. Schumacher, Annu. Rev. Condens. Matter Phys.16, 121 (2025)
2025
-
[17]
Apolin´ ario, L
G.B. Apolin´ ario, L. Moriconi, R.M. Pereira e V.J. Valad˜ ao, Phys. Rev. E102, 041102(R) (2020)
2020
-
[18]
Moriconi, Phys
L. Moriconi, Phys. Rev. E103, 062137 (2021)
2021
-
[19]
Moriconi, R.M
L. Moriconi, R.M. Pereira e V.J. Valad˜ ao, Phys. Rev. E106, L023101 (2022)
2022
-
[21]
Moriconi, R.M
L. Moriconi, R.M. Pereira e V.J. Valad˜ ao, Phys. Rev. E109, 045106 (2024)
2024
-
[22]
Iyer and L
K.P. Iyer and L. Moriconi, Phys. Fluids36, 085139 (2024)
2024
-
[23]
Z.-S. She, E. Jackson e S.A. Orszag, Nature344, 226 (1990)
1990
-
[24]
Farge, G
M. Farge, G. Pellegrino e K. Schneider, Phys. Rev. Lett.87, 054501 (2001)
2001
-
[25]
Yokokawa, K
M. Yokokawa, K. Itakura, A. Uno, T. Ishihara e Y. Kaneda, Phys. Fluids15, L21 (2003)
2003
-
[26]
Ishihara, T
T. Ishihara, T. Gotoh e Y. Kaneda, Annu. Rev. Fluid Mech.41, 165 (2009)
2009
-
[27]
B¨ urger, M
K. B¨ urger, M. Treib, R. Westermann, S. Werner, C.C. Lalescu, A. Szalay, C. Meneveau, G.L. Eyink,Vortices within vortices: hierarchical nature of vortex tubes in turbulence,arxiv.org/ abs/1210.3325(2012)
2012 arXiv
-
[28]
Ghira, G.E
A.A. Ghira, G.E. Elsinga, and C.B. da Silva, Phys. Rev. Fluids7, 104605 (2022)
2022
-
[29]
Migdal, Phys
A.A. Migdal, Phys. Rep.102, 199 (1983)
1983
-
[30]
Migdal, Int
A.A. Migdal, Int. J. Mod. Phys. A9, 1197 (1994)
1994
-
[31]
N. Cao, S. Chen, and K.R. Sreenivasan, Phys. Rev. Lett.76, 616 (1996)
1996
-
[32]
Iyer, K.R
K.P. Iyer, K.R. Sreenivasan e P.K. Yeung, Phys. Rev. X9, 041006 (2019)
2019
-
[33]
Iyer, S.S
K.P. Iyer, S.S. Bharadwaj e K.R. Sreenivasan, Proc. Nat. Acad. Sci. U.S.A.118, e2114679118 (2021)
2021
-
[34]
Moriconi, Proc
L. Moriconi, Proc. Nat. Acad. Sci. U.S.A.118, e2117968118 (2021)
2021
-
[36]
Rhodes and V
R. Rhodes and V. Vargas, Probab. Surv.11, 315 (2014). 20
2014
-
[37]
Kholmyansky, L
M. Kholmyansky, L. Moriconi, and A. Tsinober, Phys. Rev. E76, 026307 (2007)
2007
-
[38]
Y. Li, E. Perlman, M. Wan, Y. Yang, R. Burns, C. Meneveau, R. Burns, S. Chen, A. Szalay e G. Eyink, J. Turbul.9N31 (2008)
2008
-
[39]
Zhou, R.J
J. Zhou, R.J. Adrian, S. Balachandar, and T.M. Kendall, J. Fluid Mech.387, 353 (1999)
1999
-
[40]
Tang, R.A
S.L. Tang, R.A. Antonia, L. Djenidi e Y. Zhou, J. Fluid Mech.891, A26 (2020)
2020
-
[41]
Frisch,Turbulence, Cambridge University Press (1995)
U. Frisch,Turbulence, Cambridge University Press (1995)
1995
-
[42]
Yeung, D.A
P.K. Yeung, D.A. Donzis, and K.R. Sreenivasan, J. Fluid Mech.700, 5 (2012)
2012
-
[43]
Yeung, X.M
P.K. Yeung, X.M. Zhai, and K.R. Sreenivasan, Proc. Nat. Acad. Sci. U.S.A.112, 12633 (2015)
2015
-
[44]
Yeung, K.R
P.K. Yeung, K.R. Sreenivasan, and S.B. Pope, Phys. Rev. Fluids3, 064603 (2018)
2018
-
[45]
Yeung, K
P.K. Yeung, K. Ravikumar, S. Nichols, and R. Uma-Vaideswaran, Comput. Phys. Commun. 306, 109364 (2025). [46]https://turbulence.idies.jhu.edu/home
2025
-
[47]
Okabe, B
A. Okabe, B. Boots, K. Sugihara, and S.N. Chiu,Spatial Tessellations: Concepts and Appli- cations of Voronoi Diagrams(2 ed.), Wiley, (2000). 21
2000
-
[2556]
Gaussian puzzle
ForR λ = 610, two different simulations are analyzed: one in which the truncation wavenumberk max obeysk maxη≈2.67, and a higher resolution one, in whichk maxη≈5.34 (hereafter identified as “h.r.”). Following previous studies on the VGM [19–21], individual vortices are detecte...
2023
Reviewed August 5, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.