The statistics and structure of dissipation in subsonic and supersonic turbulence

Using high-resolution simulations, this study reveals that kinetic energy dissipation in subsonic turbulence is vorticity-dominated, localized on small scales, and lags energy injection by approximately 1.64 turnover times, whereas supersonic dissipation is density-correlated, spans multiple scales via shocks and vorticity, and lags by only 0.48 turnover times, with distinct fractal structures identified in both regimes.

The Gist

Imagine the universe is filled with a giant, invisible ocean of gas. Sometimes this gas moves slowly and smoothly (like a gentle breeze), and other times it moves violently and chaotically (like a hurricane). This chaotic motion is called turbulence.

This paper is like a high-speed, ultra-detailed movie camera filming that gas to answer a simple but tricky question: How does the energy of this chaotic motion turn into heat?

In the real world, when you rub your hands together, friction turns motion into heat. In space, there's no air to rub against, but the gas itself has "internal friction." When the gas swirls and crashes, that motion energy disappears and turns into heat, which can trigger chemical reactions or even help form new stars.

The authors, Edward and Christoph, used supercomputers to simulate this gas at two different speeds:

1. Subsonic (Slow): Moving slower than the speed of sound (like a gentle breeze).
2. Supersonic (Fast): Moving faster than the speed of sound (like a supersonic jet).

Here is what they discovered, explained with everyday analogies:

1. The "Time Lag" Problem

Imagine you push a child on a swing (injecting energy). The child doesn't stop swinging immediately; it takes a moment for the energy to fade away.

• The Finding: The researchers found that in the slow (subsonic) gas, it takes a long time for the energy to turn into heat—about 1.6 times the time it takes for a big swirl to cross the room.
• The Finding: In the fast (supersonic) gas, the energy turns into heat much faster—less than half that time.
• The Analogy: Think of the slow gas like a spinning top. It wobbles and spins for a long time before finally falling over and stopping. The fast gas is like a crashing car; the energy is released in a violent, instant explosion.

2. What Causes the Heat? (The "Shape" of the Chaos)

The researchers looked at what structures in the gas are actually doing the work of turning motion into heat.

• In the Slow Gas (Subsonic):

  • The Culprit: Swirling vortices (like tiny tornadoes).
  • The Shape: Imagine spaghetti. The heat is generated in long, thin, twisting strands of gas that are wrapped around each other.
  • The Rule: The heat depends on how fast the gas is spinning (vorticity), not on how dense the gas is. It's like how a spinning top heats up due to its spin, regardless of how heavy it is.

• In the Fast Gas (Supersonic):

  • The Culprit: Shockwaves (like sonic booms).
  • The Shape: Imagine sheets of paper crashing into each other. When these thin sheets collide, they crumple into filaments (like crumpled paper balls).
  • The Rule: The heat depends heavily on how dense the gas is. The denser the clump of gas, the more heat is generated. It's like a heavy truck crashing into a wall creates more heat than a bicycle crash.

3. The "Resolution" Puzzle

To see these tiny details, you need a very high-resolution camera. The researchers tried different camera settings (grid sizes) from 256 pixels up to 2048 pixels.

• The Finding: For the fast gas, even a medium-resolution camera was enough to see the picture clearly. The heat was spread out everywhere.
• The Finding: For the slow gas, even their most powerful camera (2048 pixels) couldn't quite capture the tiny details perfectly. The heat is so concentrated in tiny, microscopic swirls that it's incredibly hard to simulate accurately. It's like trying to photograph a single grain of sand with a telescope meant for stars; you need a much more powerful lens to see the details.

4. The "Fractal" Dimension (How "Full" is the Space?)

The researchers asked: "If you were to paint all the places where heat is being made, how much of the room would you cover?"

• In the Slow Gas:

  • On small scales, the heat is in thin sheets (like a flat piece of paper).
  • On large scales, it fills up the whole room like a cloud.
  • Analogy: It starts as a thin sheet of smoke, but as you zoom out, it looks like a thick fog filling the entire space.

• In the Fast Gas:

  • The heat is mostly in lines and filaments (like crumpled wires).
  • Analogy: Imagine a room filled with thousands of crumpled pieces of wire. They don't fill the room like a fog; they form a messy, tangled web of lines.

Why Does This Matter?

Understanding exactly how and where this gas turns into heat is crucial for astrophysicists.

• Star Formation: If the gas gets too hot, it can't collapse to form stars. If it stays cool, stars are born.
• Chemistry: Heat triggers chemical reactions that create the building blocks of life.
• Weather & Engineering: The math used here also helps us understand turbulence in jet engines and weather patterns on Earth.

In Summary:
This paper is a deep dive into the "friction" of the universe. It tells us that slow, gentle turbulence is a game of spinning vortices that takes a long time to cool down, while fast, violent turbulence is a game of crashing shockwaves that turns energy into heat almost instantly. The slow version is so complex and tiny that even our best supercomputers struggle to map it perfectly, while the fast version is easier to see but happens in a very different, "crumpled" way.

Technical Summary

Here is a detailed technical summary of the paper "The statistics and structure of dissipation in subsonic and supersonic turbulence" by Troccoli and Federrath (2026).

1. Problem Statement

Turbulence is a fundamental physical process governing energy transfer from large scales (driving) to small scales (dissipation). While turbulence drives star formation and interstellar medium (ISM) thermodynamics, the specific structure and statistics of kinetic energy dissipation ($\varepsilon_{kin}$) remain poorly understood.

• The Gap: It is unclear how dissipation is spatially organized (morphology), how it correlates with fluid properties (density vs. vorticity), and whether numerical simulations can accurately capture these statistics across different regimes (subsonic vs. supersonic).
• The Challenge: Dissipation is highly intermittent and occurs on small scales. Achieving numerical convergence for dissipation rates, particularly in subsonic turbulence, is computationally demanding and often plagued by numerical artifacts.

2. Methodology

The authors employed high-resolution 3D hydrodynamic simulations to study fully developed turbulence.

• Numerical Setup:

  • Code: Modified version of the FLASH code using the HLL5R 5-wave Riemann solver.
  • Grid: Uniform, triply-periodic 3D grids with resolutions ranging from $256^3$to$2048^3$ cells.
  • Physics: Compressible hydrodynamics with an isothermal equation of state ($p = c_s^2 \rho$).
  • Viscosity: Explicit kinematic viscosity ($\nu$) was included to control the Reynolds number ($Re = 2500$) and minimize numerical dissipation.
  • Driving: An Ornstein-Uhlenbeck (OU) process drove turbulence at large scales ($k_{driv}=2$) to maintain constant injection rates.

• Novel Dissipation Measurement Technique:

  • Instead of relying solely on numerical diffusion, the authors evolved a separate kinetic energy equation (Eq. 4) alongside the standard continuity and momentum equations.
  • The "perfect" kinetic energy ($E_{kin}$) was evolved without dissipation terms.
  • The actual kinetic energy was calculated from the standard momentum solution ($\frac{1}{2}\rho|v|^2$).
  • Dissipation Rate ($\varepsilon_{kin}$): Defined locally as the difference between the evolved $E_{kin}$ and the actual kinetic energy per unit time step. This method captures both explicit viscous and numerical dissipation.

• Analysis Tools:

  • Time Evolution: Tracking injection vs. dissipation rates to determine time lags.
  • Correlations: 2D Probability Distribution Functions (PDFs) of $\varepsilon_{kin}$ vs. density ($\rho$) and vorticity ($|\nabla \times v|$).
  • Spectral Analysis: Power spectra of velocity and dissipation rate to identify scaling laws and convergence.
  • Fractal Analysis: Mass-radius method to determine the fractal dimension ($D$) of dissipation structures at different scales.

3. Key Contributions

1. Direct Dissipation Measurement: Implementation of a robust method to isolate and measure local kinetic energy dissipation rates by comparing evolved energy equations with standard hydrodynamic solutions.
2. Resolution Convergence Study: A systematic study up to $2048^3$resolution, highlighting that subsonic dissipation is extremely difficult to converge numerically, whereas supersonic dissipation converges at lower resolutions ($\sim 1024^3$).
3. Regime-Specific Statistics: Quantitative differentiation of dissipation mechanisms in subsonic (vorticity-dominated) vs. supersonic (shock-dominated) regimes.

4. Key Results

A. Time Evolution and Lag

• Convergence: Subsonic dissipation rates ($\varepsilon_{kin}$) showed strong resolution dependence, only converging at $N \ge 1024$. Supersonic rates converged faster.
• Injection-Dissipation Lag: There is a measurable time delay between energy injection and dissipation:
  • Subsonic ($M=0.2$): Lag $\approx 1.64 \pm 0.21$ turbulent turnover times ($t_{turb}$).
  • Supersonic ($M=5$): Lag $\approx 0.48 \pm 0.07$ $t_{turb}$.
  • Implication: Energy dissipates $\sim 3\times$ faster in supersonic turbulence due to shock "jumping" the cascade and non-local energy transfer.

B. Correlations (Morphology)

• Subsonic Regime:
  • No correlation with density ($\rho$).
  • Strong correlation with vorticity: $\varepsilon_{kin} \propto |\nabla \times v|^2$.
  • Dissipation occurs in shear layers wrapped around vortex tubes.
• Supersonic Regime:
  • Strong correlation with density: $\varepsilon_{kin} \propto \rho^{3/2}$.
  • Dissipation is dominated by shocks (thin, dense sheets).
  • Vorticity plays a secondary role, generated at shock intersections.

C. Spectral Analysis

• Velocity Spectra: Subsonic follows Kolmogorov scaling ($k^{-5/3}$); Supersonic follows Burgers scaling ($k^{-2}$).
• Dissipation Spectra:
  • Subsonic: Does not converge even at $2048^3$. Dissipation is highly concentrated at small scales ($k \sim 30-50$), but numerical noise persists at grid scales.
  • Supersonic: Converges at $N \ge 1024$. Dissipation is distributed across a broad range of scales with weak peaks at $k \sim 10$ (shock length) and $k \sim 100$ (shock width).

D. Fractal Dimension ($D$)

• Subsonic:
  • Near dissipation scale ($\ell_\nu$): $D \approx 1.99$ (Sheet-like structures).
  • Large scales: $D \approx 2.56$ (Volume-filling).
  • Interpretation: Dissipation forms flattened vortex filaments embedded in shear layers.
• Supersonic:
  • Near dissipation scale: $D \approx 1.60$.
  • Large scales: $D \approx 1.35$.
  • Interpretation: A mix of 2D shock sheets and 1D filamentary intersections. The lower dimension at large scales indicates a dominance of filamentary structures formed by intersecting shocks.

5. Significance and Future Directions

• Astrophysical Impact: The results provide critical constraints for modeling the thermodynamics and chemistry of the Interstellar Medium (ISM). Since dissipation drives heating and chemical reactions, the distinct timescales and spatial structures in subsonic vs. supersonic regimes must be accounted for in star formation models.
• Numerical Challenges: The study highlights that achieving numerical convergence for subsonic dissipation is exceptionally difficult, requiring resolutions beyond current standard capabilities ($>2048^3$) for full accuracy.
• Future Work: The authors suggest extending this methodology to Magnetohydrodynamics (MHD) to include magnetic fields and further investigating the "bottleneck effect" in subsonic turbulence. They also propose quantifying the fraction of total dissipation arising from the most dissipative volume fractions to better model intermittent heating.

In summary, this paper establishes that while supersonic turbulence dissipates energy rapidly via shocks across many scales, subsonic turbulence dissipates energy slowly via vorticity-dominated shear layers on small scales, with fundamentally different statistical and fractal properties.