Exact formulas for arbitrary order velocity-gradient moments in isotropic turbulence
Pith reviewed 2026-05-09 20:45 UTC · model grok-4.3
The pith
Longitudinal velocity gradient moments above third order depend on both dissipation and strain self-amplification in isotropic turbulence.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The authors derive exact expressions for arbitrary-order moments of longitudinal and transverse velocity gradients in isotropic turbulence by combining isotropic tensor theory, orientational averaging, and algorithmic implementation. These expressions are written in terms of invariants of the velocity gradient tensor. In particular, longitudinal velocity gradient moments of order higher than three depend not only on tr(S²), which is proportional to the dissipation rate, but also on tr(S³), which reflects strain self-amplification.
What carries the argument
Systematic derivation via isotropic tensor theory and orientational averaging that reduces all moments to invariants of the strain-rate tensor S.
Load-bearing premise
The turbulence is isotropic, which allows averaging over all orientations to collapse the moments onto a small number of strain-rate invariants.
What would settle it
A direct numerical simulation in which the measured fourth-order longitudinal velocity gradient moment fails to match the exact linear combination of tr(S²) and tr(S³) predicted by the formulas.
read the original abstract
Statistical moments of velocity gradients provide fundamental information on the small-scale properties of turbulence. In this work, we propose a systematic method to derive exact expressions for statistical moments of arbitrary order for both longitudinal and transverse velocity gradients in isotropic turbulence. The approach is applicable to both compressible and incompressible flows and expresses the moments in terms of invariants of the velocity gradient tensor. The derivation combines isotropic tensor theory, orientational averaging, and an algorithmic implementation, enabling the computation of high-order moments in a unified framework. We show that longitudinal velocity gradient moments of order higher than three depend not only on $\mathrm{tr}(\boldsymbol{S}^2)$, which is proportional to the dissipation rate, but also on $\mathrm{tr}(\boldsymbol{S}^3)$, which reflects strain self-amplification, where $\boldsymbol{S}$ denotes the strain-rate tensor. The resulting theoretical expressions are validated through comparisons with existing theoretical results and direct numerical simulations.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes a systematic method combining isotropic tensor theory, orientational averaging over the velocity-gradient tensor, and an algorithmic implementation to derive exact expressions for arbitrary-order statistical moments of both longitudinal and transverse velocity gradients in isotropic turbulence. These expressions are reduced to functions of the strain-rate invariants tr(S²) (proportional to dissipation) and tr(S³) (reflecting strain self-amplification), with the key result that longitudinal moments of order >3 depend on both invariants. The approach is stated to apply to compressible and incompressible flows and is validated by comparison with existing theory and DNS data.
Significance. If the derivations hold, the work supplies a parameter-free, exact framework for computing high-order velocity-gradient moments under isotropy, directly implementing the representation theory of the symmetric traceless strain-rate tensor in 3D (via Cayley-Hamilton). This could strengthen theoretical understanding of small-scale intermittency and strain dynamics beyond dissipation scaling alone, with potential utility in turbulence modeling and analysis of DNS data.
minor comments (4)
- The abstract and introduction refer to an 'algorithmic implementation' for arbitrary-order moments, but no pseudocode, explicit recurrence relation, or worked example (e.g., for order 4 or 5) is described; adding this would improve reproducibility and allow readers to verify the reduction to only two invariants.
- Notation for the velocity-gradient tensor and its decomposition into strain S and rotation should be introduced with a clear equation early in the manuscript to avoid ambiguity when discussing orientational averaging.
- The validation section compares with DNS and prior theory, but the specific Reynolds numbers, grid resolutions, or forcing methods used in the DNS are not stated; these details are needed to assess the strength of the numerical support.
- A brief statement on how the method extends (or does not) to compressible flows would clarify the scope, since the invariants tr(S²) and tr(S³) are defined for the traceless strain but compressibility introduces dilatation.
Simulated Author's Rebuttal
We thank the referee for their careful reading of the manuscript and for the positive summary and significance assessment. The referee's description accurately reflects the core contributions of the work. The recommendation for minor revision is noted, and we will incorporate any such changes in the revised version.
Circularity Check
No significant circularity identified
full rationale
The paper derives exact expressions for arbitrary-order velocity-gradient moments in isotropic turbulence by combining standard isotropic tensor theory with orientational averaging over the velocity-gradient distribution. This reduces all moments to invariants of the strain-rate tensor via the Cayley-Hamilton theorem, which limits independent invariants to tr(S²) and tr(S³) in 3D. The approach is algorithmic but follows directly from representation theory of the symmetric traceless tensor without introducing fitted parameters, self-definitional loops, or load-bearing self-citations. Validation against DNS and prior results is external and does not alter the derivation chain. The central claim—that higher-order longitudinal moments depend on both dissipation and strain self-amplification—is a direct consequence of isotropy assumptions and tensor algebra, not a renaming or smuggling of inputs.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption The turbulence is isotropic
Reference graph
Works this paper leans on
-
[1]
Batchelor, G. K. 1953 The theory of homogeneous turbulence\/ . Cambridge University Press
work page 1953
-
[2]
Bell, E. T. 1934 Exponential polynomials . Ann. Math. 35 (2), 258--277
work page 1934
-
[3]
Bentkamp, L. , Lalescu, C. C. & Wilczek, M. 2019 Persistent accelerations disentangle lagrangian turbulence . Nat. Commun. 10 (1), 3550
work page 2019
-
[4]
1956 An inequality concerning the production of vorticity in isotropic turbulence
Betchov, R. 1956 An inequality concerning the production of vorticity in isotropic turbulence . J. Fluid Mech. 1 (5), 497--504
work page 1956
-
[5]
Boschung, J. 2015 Exact relations between the moments of dissipation and longitudinal velocity derivatives in turbulent flows . Phys. Rev. E 92 (4), 043013
work page 2015
-
[6]
Carbone, M. & Wilczek, M. 2022 Only two betchov homogeneity constraints exist for isotropic turbulence . J. Fluid Mech. 948 , R2
work page 2022
-
[7]
1978 The fine-scale structure of the turbulent velocity field
Champagne, F.H. 1978 The fine-scale structure of the turbulent velocity field . J. Fluid Mech. 86 , 67--108
work page 1978
-
[8]
2012 Advanced Combinatorics: The art of finite and infinite expansions\/
Comtet, L. 2012 Advanced Combinatorics: The art of finite and infinite expansions\/ . Springer
work page 2012
-
[9]
Eyink, G. L. 2006 Multi-scale gradient expansion of the turbulent stress tensor . J. Fluid Mech. 549 , 159--190
work page 2006
- [10]
- [11]
-
[12]
Fang, J. , Zheltovodov, A. A. , Yao, Y. , Moulinec, C. & Emerson, D. R. 2020 On the turbulence amplification in shock-wave/turbulent boundary layer interaction . J. Fluid Mech. 897 , A32
work page 2020
-
[13]
Hierro, J. & Dopazo, C. 2003 Fourth-order statistical moments of the velocity gradient tensor in homogeneous, isotropic turbulence . Phys. Fluids. 15 (11), 3434--3442
work page 2003
-
[14]
Ishihara, T. , Gotoh, T. & Kaneda, Y. 2009 Study of high R eynolds number isotropic turbulence by D irect N umerical S imulation . Annu. Rev. Fluid Mech. 41 , 65
work page 2009
-
[15]
Ishihara, T. , Kaneda, Y. , Yokokawa, M. , Itakura, K. & Uno., A. 2007 Small--scale statistics in high--resolution direct numerical simulation of turbulence: R eynolds number dependence of one--point velocity gradient . J. Fluid Mech. 592 , 335
work page 2007
-
[16]
Johnson, P. L. & Wilczek, M. 2024 Multiscale velocity gradients in turbulence . Annu. Rev. Fluid Mech. 56 (1), 463--490
work page 2024
- [17]
- [18]
-
[19]
Annual Review of Fluid Mechanics 43 (1), 219--245
Meneveau, Charles 2011 Lagrangian dynamics and models of the velocity gradient tensor in turbulent flows . Annual Review of Fluid Mechanics 43 (1), 219--245
work page 2011
-
[20]
Pennisi, S. & Trovato, M. 1987 On the irreducibility of professor g.f. smith's representations for isotropic functions . Int. J. Eng. Sci. 25 (8), 1059--1065
work page 1987
-
[21]
Petersen, M. R. & Livescu, D. 2010 Forcing for statistically stationary compressible isotropic turbulence . Phys. Fluids. 22 (11)
work page 2010
-
[22]
Phan-Thien, N. & Antonia, R. A. 1994 Isotropic cartesian tensors of arbitrary even orders and velocity gradient correlation functions . Phys. Fluids. 6 (12), 3818--3822
work page 1994
-
[23]
Pope, S. B. 1975 A more general effective-viscosity hypothesis . J. Fluid Mech. 72 (2), 331--340
work page 1975
-
[24]
Pope, S. B. 2000 Turbulent Flows\/ . Cambridge University Press
work page 2000
-
[25]
Robertson, H. P. 1940 The invariant theory of isotropic turbulence. In Math. proc. Camb. Philos.\/ , , vol. 36 , pp. 209--223 . Cambridge University Press
work page 1940
-
[26]
Siggia, E. D. 1981 Invariants for the one-point vorticity and strain rate correlation functions . Phys. Fluids. 24 (11), 1934--1936
work page 1981
-
[27]
Taylor, G. I. 1935 Statistical theory of turbulence . Proc. Roy. Soc. London. Ser. A, Math. Phys. Sci. 151 , 421--444
work page 1935
-
[28]
Wallace, J. M. 2009 Twenty years of experimental and direct numerical simulation access to the velocity gradient tensor: What have we learned about turbulence? Phys. Fluids. 21 (2), 021301
work page 2009
-
[29]
1946 The classical groups: Their invariants and representations\/
Weyl, H. 1946 The classical groups: Their invariants and representations\/ . Princeton University Press
work page 1946
- [30]
- [31]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.