Skip to content

Turbulence: a pattern that repeats over many decades of scale, and the two sizes where it stops #19

Closed

Description

@ahrazzle

Question

A turbulent flow looks statistically the same at every size between the scale where energy is stirred in and the scale where it turns into heat. Can a small model predict both ends of that range from first principles, and reproduce the measured slopes in between, including where they deviate?

  • Lower end in air and water: the Kolmogorov scale η = (ν³/ε)^¼, about a millimetre in the atmosphere and upper ocean.
  • Lower end in the solar wind: the ion scale. The proton inertial length is about 100 km at 1 AU, which by Taylor's hypothesis puts the spectral break near 0.6 Hz.
  • Lower end in molecular clouds: the sonic scale, where the velocity spread falls to the sound speed, about 0.07 pc from σ = 0.72 S^0.5 km/s at 10 K.

Why it matters

Turbulence is the best-measured case of real self-similarity with both cutoffs physically understood. The interstellar electron-density spectrum follows one Kolmogorov-like power law over many decades (the "big power law in the sky", Armstrong, Rickett and Spangler, ApJ 443, 209 (1995)). The repetition is not perfect, and the model can show where:

  • intermittency (She and Leveque, PRL 72, 336 (1994));
  • a change of cascade in molecular clouds (Larson's σ ∝ L^0.38, MNRAS 194, 809 (1981), against the supersonic 0.5 of Solomon et al. 1987);
  • a 2D-like k⁻³ range in the atmosphere above ~500 km (Nastrom and Gage, J. Atmos. Sci. 42, 950 (1985)).

It also meets the quantum thread: in superfluid helium the classical cascade hands over to quantised vortices, with circulation h/m.

First steps

  • turbulence/model.js: η and L/η ~ Re^¾; the solar-wind break from the ion inertial length; the molecular-cloud sonic scale.
  • Checks against published numbers: the exact 4/5 law; the Kolmogorov constant (~1.5); She–Leveque exponents against measured ones; the solar-wind break at 0.1–1 Hz with a kinetic-range slope near −2.8 (Bruno and Carbone, Living Rev. Solar Phys. 10, 2 (2013)); the sonic scale against the ~0.1 pc coherent cores and filament widths. The numbers in this issue are from recall in a session whose literature search was exhausted, so verify them first.
  • Larson's 0.38 against 0.5: which cascade the cloud data favour, and where the analogy with lab turbulence breaks (compressible, magnetised, self-gravitating).
  • Cloud edges: their perimeter-area dimension (~1.36, Falgarone et al. 1991) does not follow simply from a Kolmogorov spectrum (a fractional-Brownian estimate gives 7/6). Test whether the agreement people cite is real or a coincidence.

Activity

  1. ahrazzle commented on Oct 18, 2026

    @ahrazzle
    ContributorAuthor

    I'm taking this one. The plan follows the issue's first steps, plus a cascade model so that the "where it deviates" part is computed rather than quoted:

    • turbulence/model.js, with no dependencies:
      • the Kolmogorov scale, and L/η ~ Re^¾;
      • the solar-wind ion break, from the proton inertial length and gyroradius, using Taylor's hypothesis;
      • the molecular-cloud sonic scale, from linewidth–size relations;
      • a seeded Sabra shell model (L'vov et al. 1998) for the cascade itself. It is deterministic and integer-free except for IEEE arithmetic.
    • tests/turbulence.mjs:
      • the shell model's constant energy flux, which is its exact analogue of the 4/5 law;
      • its structure-function exponents ζ_p, compared with K41 (p/3), She–Leveque, and the published Sabra and experimental values;
      • η for measured ε in air and the upper ocean;
      • the solar-wind break, against 0.1 to 1 Hz;
      • the sonic scale, against core and filament widths.
    • A README section on where the repetition is real and where it breaks:
      • intermittency;
      • Larson's 0.38 against 0.5;
      • the 2D-like k⁻³ range;
      • the cloud-edge perimeter dimension against the fractional-Brownian 7/6.

    The issue's numbers are from recall, so the first step is to check each one against its primary source. Runtime will stay small enough for CI. If anyone has already started on this, say so here and I'll stop.


    Generated by Claude Code

  2. ahrazzle commented on Oct 18, 2026

    @ahrazzle
    ContributorAuthor

    #34 is up for this issue.

    Several of this issue's numbers, which I wrote from recall, did not hold up. CORRECTIONS 19–20 record them:

    • Solomon et al.'s 0.72 goes with the effective radius, not with S.
    • The 0.6 Hz solar-wind break needs sin θ_BV. With it the estimate is 0.46 Hz, inside the observed range.
    • The 7/6 for cloud edges came from the formula for all contours at a level. A perimeter–area fit measures single loops, (3 − H)/2, which a computation on Gaussian surfaces in the PR confirms.

    The issue's question about cloud edges has an answer in #34. Projected Kolmogorov density gives 1.15, not 1.36. The laboratory 1.36 comes from interface sections, a different measurement, so the match does not show shared physics.

    Every literature value in the PR was checked only against search text. The PR body lists the ones that need a check against the papers.


    Generated by Claude Code

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Metadata

Metadata

Assignees

No one assigned

    Labels

    enhancementNew feature or request

    Type

    No type

    Projects

    No projects

      Milestone

      No milestone

      Relationships

      None yet

      Development

      No branches or pull requests

      Issue actions