pith. sign in

arxiv: 2510.06620 · v1 · submitted 2025-10-08 · 🌊 nlin.CD · physics.flu-dyn

Ultra-chaotic property of Navier-Stokes turbulence

Shijie Qin , Kun Xu , Shijun Liao This is my paper

Pith reviewed 2026-05-18 09:41 UTC · model grok-4.3

classification 🌊 nlin.CD physics.flu-dyn
keywords ultra-chaosNavier-Stokes turbulenceKolmogorov flowinitial condition sensitivityclean numerical simulationflow symmetrystatistical dependencetwo-dimensional turbulence
0
0 comments X

The pith

Navier-Stokes turbulence statistics change sharply with tiny initial condition shifts.

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

The paper establishes that in two-dimensional Kolmogorov flow governed by the Navier-Stokes equations, minuscule changes to the initial velocity field produce not only divergent trajectories but also different flow symmetries and different long-term statistical properties. This sensitivity is demonstrated using clean numerical simulation, which keeps artificial errors far below the scale of the observed physical differences over the time windows used for averaging. If correct, the result means that the statistical behavior of Navier-Stokes turbulence cannot be regarded as uniquely determined by the governing equations; small disturbances must be taken into account even when only averages are sought. The authors note that this creates a logical tension because practical turbulence models routinely omit such disturbances.

Core claim

Applying clean numerical simulation to the two-dimensional turbulent Kolmogorov flow, the authors find that extremely small variations in the initial conditions cause large differences in the spatiotemporal evolution, in the symmetry properties of the flow, and in its time-averaged statistics. They therefore conclude that Navier-Stokes turbulence is ultra-chaotic, in the sense that its statistics themselves depend sensitively on initial data, and that small disturbances cannot be neglected when describing turbulence from a statistical viewpoint.

What carries the argument

Clean numerical simulation (CNS), a technique that suppresses artificial numerical noise to levels negligible compared with physical effects over long integration intervals used for statistics.

If this is right

  • Statistical quantities extracted from Navier-Stokes turbulence simulations can differ substantially for initial conditions that are indistinguishable at machine precision.
  • Flow symmetry is not robust under arbitrarily small initial perturbations.
  • Standard turbulence closures that start from the Navier-Stokes equations without explicit disturbance terms may yield inconsistent statistical predictions.
  • Any practical model of Navier-Stokes turbulence must incorporate the effects of small disturbances when statistics are the target quantity.

Where Pith is reading between the lines

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

  • Ensemble or stochastic representations may be required to describe the full set of possible statistical outcomes.
  • Limits on the predictability of turbulent statistics could extend well beyond the usual Lyapunov-time horizon for individual trajectories.
  • Laboratory experiments that introduce controlled micro-perturbations into grid turbulence or channel flow could provide an independent test of statistical sensitivity.

Load-bearing premise

The clean numerical simulation method can reduce all artificial numerical noise to a level negligible compared with the physical differences caused by the initial-condition variations over the time interval used for statistics.

What would settle it

Repeating the Kolmogorov-flow runs with initial perturbations an order of magnitude smaller than those reported and verifying whether the statistical divergences persist at the same relative magnitude.

Figures

Figures reproduced from arXiv: 2510.06620 by Kun Xu, Shijie Qin, Shijun Liao.

Figure 1
Figure 1. Figure 1: Vorticity fields ω at t = 100 of the 2D turbulent Kolmogorov flow governed by (1) and (2) in the case of nK = 16 and Re = 2000, given by CNS subject to the initial conditions (3) (left, marked by Flow CNS-1), (4) (middle, marked by Flow CNS-2), and (5) (right, marked by Flow CNS-3), respectively. solution and thus is negligible in a long enough interval of time, we use the CNS [14–20] to solve Eq. (1) subj… view at source ↗
Figure 2
Figure 2. Figure 2: Comparison of time histories of the spatially averaged enstrophy dissipation rate [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Comparison of probability density functions (PDFs) of (top-left) kinetic energy [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Comparison of the spatiotemporally averaged (left) kinetic energy dissipation rate [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Comparison of the spatiotemporally averaged (top-left) horizontal velocity [PITH_FULL_IMAGE:figures/full_fig_p007_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Comparison of the temporal averaged kinetic energy spectra [PITH_FULL_IMAGE:figures/full_fig_p007_6.png] view at source ↗
Figure 6
Figure 6. Figure 6: figure 6: the three turbulent flows have different kinetic energy spectra, say, different distributions [PITH_FULL_IMAGE:figures/full_fig_p008_6.png] view at source ↗
read the original abstract

A chaotic system is called ultra-chaos when its statistics have sensitivity dependence on initial condition and/or other small disturbances. In this paper, using two-dimensional turbulent Kolmogorov flow as an example, we illustrate that tiny variation of initial condition of Navier-Stokes equations can lead to huge differences not only in spatiotemporal trajectory but also in flow symmetry and its statistics. Here, in order to avoid the influence of artificial numerical noise, we apply ``clean numerical simulation'' (CNS) which can guarantee that the numerical noise can be reduced to such a desired low level that they are negligible in a time interval long enough for calculating statistics. This discovery highly suggests that the Navier-Stokes turbulence (i.e. turbulence governed by the Navier-Stokes equations) might be an ultra-chaos, say, small disturbances must be considered even from viewpoint of statistics. This however leads to a paradox in logic, since small disturbances, which are unavoidable in practice, are unfortunately neglected by the Navier-Stokes turbulence. Some fundamental characteristics of turbulence model are discussed and suggested in general meanings.

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.

Referee Report

3 major / 2 minor

Summary. The paper claims that Navier-Stokes turbulence exhibits an 'ultra-chaotic' property, in which long-time statistics (including flow symmetry and quantities such as energy spectra) are sensitive to tiny variations in initial conditions. This is illustrated numerically for two-dimensional turbulent Kolmogorov flow by comparing two clean numerical simulations (CNS) started from slightly different initial data; the authors argue that CNS keeps artificial numerical noise negligible over the averaging interval, implying that the observed statistical differences are physical and that standard turbulence models neglect unavoidable small disturbances, leading to a logical paradox.

Significance. If the central claim is substantiated with quantitative error control, the result would challenge the conventional assumption that turbulence statistics are robust and insensitive to initial conditions, with potential consequences for predictability, ensemble modeling, and the foundations of statistical turbulence theory. The methodological choice to employ CNS is a clear strength for controlling truncation and round-off effects, but the absence of explicit bounds on residual numerical error relative to the reported statistical differences currently limits the strength of the evidence.

major comments (3)
  1. [Numerical method and results sections] The description of the CNS procedure and the results section do not contain a quantitative comparison showing that the CNS error estimator (or the magnitude of the artificial noise term) remains smaller than the observed differences in long-time statistics (e.g., symmetry measures or energy spectra) over the chosen averaging window. Without this comparison the central claim that the statistical divergence is physical rather than numerical cannot be verified.
  2. [Results and figures] No error bars, standard deviations, or statistical significance measures are reported for the differences between the two CNS runs in any of the presented statistics; the abstract and results therefore provide only a qualitative illustration rather than a quantified demonstration of sensitivity.
  3. [Discussion and conclusions] The generalization from the 2-D Kolmogorov flow example to the broader statement that 'Navier-Stokes turbulence might be an ultra-chaos' is not supported by any test or discussion of three-dimensional cases or alternative forcings, which is required for the claim to be load-bearing.
minor comments (2)
  1. [Introduction] The term 'ultra-chaos' is introduced in the abstract and introduction without a precise mathematical definition or citation to related concepts in the literature on chaotic systems and statistical sensitivity.
  2. [Figures] Figure captions and axis labels should explicitly identify the two initial conditions being compared and indicate the time interval over which statistics are accumulated.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the constructive and detailed comments. We address each major comment below, indicating where revisions will be made to strengthen the manuscript while maintaining the integrity of our original claims.

read point-by-point responses

  1. Referee: [Numerical method and results sections] The description of the CNS procedure and the results section do not contain a quantitative comparison showing that the CNS error estimator (or the magnitude of the artificial noise term) remains smaller than the observed differences in long-time statistics (e.g., symmetry measures or energy spectra) over the chosen averaging window. Without this comparison the central claim that the statistical divergence is physical rather than numerical cannot be verified.

    Authors: We agree that an explicit quantitative comparison would provide stronger verification. In the revised manuscript we will add a new subsection in the numerical method section that reports the CNS error estimator bounds (including truncation and round-off contributions) and directly compares these bounds to the observed differences in symmetry measures and energy spectra over the full averaging window. This comparison will show that residual numerical noise remains at least an order of magnitude below the reported statistical divergences, confirming the physical nature of the sensitivity. revision: yes

  2. Referee: [Results and figures] No error bars, standard deviations, or statistical significance measures are reported for the differences between the two CNS runs in any of the presented statistics; the abstract and results therefore provide only a qualitative illustration rather than a quantified demonstration of sensitivity.

    Authors: We accept this observation. The revised results section will include error bars or standard deviations for all key statistical quantities, obtained by dividing the long-time averaging interval into multiple sub-intervals and computing variability across them. We will also add a brief discussion of the statistical significance of the differences between the two CNS runs to move beyond qualitative illustration. revision: yes

  3. Referee: [Discussion and conclusions] The generalization from the 2-D Kolmogorov flow example to the broader statement that 'Navier-Stokes turbulence might be an ultra-chaos' is not supported by any test or discussion of three-dimensional cases or alternative forcings, which is required for the claim to be load-bearing.

    Authors: The manuscript presents the two-dimensional Kolmogorov flow as a specific, well-controlled example to illustrate the ultra-chaotic property and uses cautious language ('might be', 'highly suggests') rather than asserting universality. We will expand the discussion section to provide additional reasoning based on the structure of the Navier-Stokes equations and the generic presence of small disturbances, while explicitly noting that three-dimensional tests and other forcings remain important future work. We maintain that the current illustrative evidence is sufficient for the paper's stated scope and does not require immediate three-dimensional computations to support the suggestion. revision: partial

Circularity Check

0 steps flagged

No significant circularity; result is empirical observation from controlled simulations

full rationale

The paper presents numerical evidence from clean numerical simulations (CNS) of 2D Kolmogorov flow, comparing long-time statistics (energy spectra, symmetry measures) across runs started from slightly perturbed initial conditions. The central claim—that Navier-Stokes turbulence exhibits sensitivity in its statistics—is an output of these experiments rather than a closed mathematical derivation. No equation reduces a claimed prediction to a parameter fitted from the same data, no self-definitional loop exists, and the CNS noise-control guarantee is invoked as a methodological precondition whose validity is external to the present statistics (prior method papers). The argument chain is therefore self-contained against external benchmarks and does not collapse to its inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 1 invented entities

The central claim depends on the assumption that CNS truly isolates physical sensitivity from numerical artifacts and on the interpretation that observed differences in symmetry and statistics are physically meaningful rather than transient.

axioms (1)
  • domain assumption Clean numerical simulation can keep numerical noise negligible over the time interval required to compute statistics.
    This premise is required to attribute differences solely to initial conditions.
invented entities (1)
  • ultra-chaos no independent evidence
    purpose: Label for chaotic systems whose statistics are sensitive to initial conditions or small disturbances.
    New term introduced to distinguish statistical sensitivity from ordinary trajectory chaos.

pith-pipeline@v0.9.0 · 5707 in / 1242 out tokens · 42729 ms · 2026-05-18T09:41:17.482997+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

27 extracted references · 27 canonical work pages

  1. [1]

    Clay Mathematics Institute of Cambridge, Massachusetts, The Millennium Prize Problems, https://www.claymath.org/millennium-problems/(2000)

    work page 2000

  2. [2]

    R. G. Deissler, Is Navier-Stokes turbulence chaotic?, Phys. Fluids 29 (1986) 1453 – 1457

    work page 1986

  3. [3]

    Boffetta, S

    G. Boffetta, S. Musacchio, Chaos and predictability of homogeneous-isotropic turbulence, Phys. Rev. Lett. 119 (5) (2017) 054102

    work page 2017

  4. [4]

    Berera, R

    A. Berera, R. D. J. G. Ho, Chaotic properties of a turbulent isotropic fluid, Phys. Rev. Lett. 120 (2) (2018) 024101

    work page 2018

  5. [5]

    J. Ge, J. Rolland, J. C. Vassilicos, The production of uncertainty in three-dimensional Navier-Stokes turbulence, J. Fluid Mech. 977 (2023) A17

    work page 2023

  6. [6]

    E. N. Lorenz, Deterministic nonperiodic flow, J. Atmos. Sci. 20 (2) (1963) 130–141

    work page 1963

  7. [7]

    S. Liao, S. Qin, Ultra-chaos: an insurmountable objective obstacle of reproducibility and replication, Adv. Appl. Math. Mech. 14 (4) (2022) 799–815

    work page 2022

  8. [8]

    Y. Yang, S. Qin, S. Liao, Ultra-chaos of a mobile robot: A higher disorder than normal-chaos, Chaos Solitons Fractals 167 (2023) 113037

    work page 2023

  9. [9]

    S. Qin, S. Liao, A kind of Lagrangian chaotic property of the Arnold-Beltrami-Childress flow, J. Fluid Mech. 960 (2023) A15

    work page 2023

  10. [10]

    Zhang, S

    B. Zhang, S. Liao, Ultra-chaos in a meandering jet flow, Physica D 455 (2023) 133886

    work page 2023

  11. [11]

    A. M. Obukhov, Kolmogorov flow and laboratory simulation of it, Russian Math. Surveys 38 (4) (1983) 113–126

    work page 1983

  12. [12]

    G. J. Chandler, R. R. Kerswell, Invariant recurrent solutions embedded in a turbulent two- dimensional Kolmogorov flow, J. Fluid Mech. 722 (2013) 554–595

    work page 2013

  13. [13]

    W. Wu, F. G. Schmitt, E. Calzavarini, L. Wang, A quadratic Reynolds stress development for the turbulent Kolmogorov flow, Phys. Fluids 33 (2021) 125129

    work page 2021

  14. [14]

    Liao, On the reliability of computed chaotic solutions of non-linear differential equations, Tellus Ser

    S. Liao, On the reliability of computed chaotic solutions of non-linear differential equations, Tellus Ser. A-Dyn. Meteorol. Oceanol. 61 (4) (2009) 550–564

    work page 2009

  15. [15]

    Liao, Clean Numerical Simulation, Chapman and Hall/CRC, 2023

    S. Liao, Clean Numerical Simulation, Chapman and Hall/CRC, 2023

    work page 2023

  16. [16]

    T. Hu, S. Liao, On the risks of using double precision in numerical simulations of spatio- temporal chaos, J. Comput. Phys. 418 (2020) 109629

    work page 2020

  17. [17]

    S. Qin, S. Liao, Influence of numerical noises on computer-generated simulation of spatio- temporal chaos, Chaos Solitons Fractals 136 (2020) 109790

    work page 2020

  18. [18]

    S. Qin, S. Liao, Large-scale influence of numerical noises as artificial stochastic disturbances on a sustained turbulence, J. Fluid Mech. 948 (2022) A7. 12

    work page 2022

  19. [19]

    S. Liao, S. Qin, Noise-expansion cascade: an origin of randomness of turbulence, J. Fluid Mech. 1009 (2025) A2

    work page 2025

  20. [20]

    S. Liao, S. Qin, Physical significance of artificial numerical noise in direct numerical simu- lation of turbulence, J. Fluid Mech. 1008 (2025) R2

    work page 2025

  21. [21]

    S. Qin, Y. Yang, Y. Huang, X. Mei, L. Wang, S. Liao, Is a direct numerical simulation (DNS) of Navier-Stokes equations with small enough grid spacing and time-step definitely reliable/correct?, Journal of Ocean Engineering and Science 9 (2024) 293 – 310

    work page 2024

  22. [22]

    S. B. Pope, Turbulent Flows, IOP Publishing, 2001

    work page 2001

  23. [23]

    Boffetta, R

    G. Boffetta, R. E. Ecke, Two-dimensional turbulence, Annu. Rev. Fluid Mech. 44 (2012) 427–51

    work page 2012

  24. [24]

    S. Qin, S. Liao, A self-adaptive algorithm of the clean numerical simulation (CNS) for chaos, Adv. Appl. Math. Mech. 15 (5) (2023) 1191–1215

    work page 2023

  25. [25]

    L. D. Landau, E. M. Lifshitz, Course of Theoretical Physics: Fluid Mechanics (Vol. 6), Addision-Wesley, Reading, 1959

    work page 1959

  26. [26]

    X. Yang, K. Xu, Wave-particle based multiscale modeling and simulation of non-equilibrium turbulent flows, arXiv preprint arXiv:2503.07207

    work page arXiv

  27. [27]

    X. Yang, K. Xu, Wave-particle turbulence simulation of spatially developing round jets: Turbulent flow modeling and method validation, arXiv preprint arXiv:2507.14524

    work page arXiv