Scaling in Supersonic Turbulence: Energy Spectra and Fluxes using High-Fidelity Direct Numerical Simulations

Using high-resolution GPU-accelerated direct numerical simulations, this study reveals that supersonic turbulence undergoes a fundamental shift in energy cascade mechanisms, characterized by a transition from Kolmogorov-like to Burgers-like scaling in rotational energy spectra driven by dominant cross-scale energy transfer from solenoidal to compressive modes.

The Gist

Imagine a giant, invisible ocean of gas filling the universe. Sometimes this gas flows smoothly, like a gentle river. Other times, it goes wild, churning, crashing, and forming shockwaves like sonic booms. This chaotic state is called turbulence.

When this gas moves slower than the speed of sound, we know a lot about how it behaves. But when it moves faster than sound (supersonic)—like in exploding stars or high-speed rocket engines—it becomes a mystery. Scientists have struggled to understand how energy moves through this super-fast chaos.

This paper is like a high-definition movie of that chaos, created by running a massive computer simulation. Here is what the researchers found, explained simply:

1. The Super-Computer Movie

The authors built a virtual box of gas and used a super-powerful computer (specifically, a machine with 128 advanced graphics cards) to simulate it. They didn't just guess; they solved the actual physics equations for gas moving at different speeds, from slightly below the speed of sound to three times faster.

They used a special "camera" (a mathematical method called TENO) that is sharp enough to see tiny swirls of gas and the incredibly thin, sharp lines where shockwaves crash, without blurring them out.

2. The Two Types of "Dance Moves"

In this gas, there are two main ways the particles move:

• The Spin (Rotational): Like a spinning top or a whirlpool.
• The Squeeze (Compressive): Like a piston pushing air, creating compression waves or shockwaves.

In slow (subsonic) gas, the "Spin" moves energy in a predictable, steady way, like a waterfall flowing down steps. This is the famous "Kolmogorov" pattern that scientists have known for decades.

3. The Big Surprise: The Rules Change at High Speeds

The researchers discovered that once the gas goes supersonic, the rules of the game change completely.

• The Spin gets tired: As the gas speeds up, the "Spin" energy stops flowing smoothly. Instead of a steady waterfall, it becomes a steep slide. The energy drains away faster than expected.
• The Squeeze gets weird: The "Squeeze" energy, which usually behaves like a specific type of wave (Burgers turbulence), actually gets flatter and more spread out as the speed increases.

The Analogy: Imagine a crowded dance floor.

• In slow motion, everyone spins in their own spot, and the energy stays local.
• In supersonic motion, the dancers start bumping into each other so hard that the "spinners" start transferring their energy to the "squeezers." The spinners lose their energy to the shockwaves, and the shockwaves get a weird, flatter distribution of energy.

4. The "Handoff" Between Modes

The most important discovery is a massive energy handoff.
In slow gas, the spinning motion and the squeezing motion barely talk to each other. But in supersonic gas, the spinning motion (which the researchers forced into the system) aggressively dumps its energy into the squeezing motion.

Think of it like a relay race where the runner (spin) doesn't just pass the baton to the next runner; they actually throw the baton into the air, and the other runner (squeeze) has to catch it while running through a wall. This "cross-talk" is what changes the shape of the energy patterns.

5. Shockwaves are the New Boss

As the gas gets faster, the "Squeeze" motion becomes dominated by shockwaves (sudden, violent jumps in pressure).

• The researchers found that the behavior of these shockwaves in supersonic gas follows a very old, simple mathematical rule called Burgers turbulence.
• It's as if, despite the complexity of the gas, the shockwaves simplify the chaos into a predictable pattern: the stronger the shock, the more energy it carries, following a specific "cube" relationship.

6. What This Means for the Paper's Claims

The paper concludes that you cannot use the old "slow gas" rules to understand "fast gas."

• Old View: Energy flows smoothly from big swirls to small swirls.
• New View: In supersonic gas, energy is constantly being stolen from the swirls and dumped into shockwaves and heat (pressure dilatation). This changes the entire landscape of how the gas moves.

The researchers did not claim this solves problems in medicine or specific engineering designs yet. They simply provided the "blueprint" of how energy moves in this extreme environment, showing that the interaction between spinning gas and shockwaves is the key to understanding the chaos.

In a nutshell: Supersonic turbulence isn't just "fast" turbulence; it's a different beast entirely where the spinning motion gets hijacked by shockwaves, creating a new set of rules for how energy travels through the universe.

Technical Summary

1. Problem Statement

Supersonic turbulence is critical to understanding astrophysical phenomena (e.g., star formation, supernovae) and high-speed engineering flows. However, the mechanisms governing energy transfer in these regimes remain poorly understood compared to incompressible turbulence.

• The Gap: Previous studies often relied on the Euler equations (inviscid), Large Eddy Simulations (LES), or low-order numerical schemes (like PPM) that introduce excessive numerical dissipation. This obscures the resolution of fine-scale eddies and sharp shock fronts, leading to conflicting results regarding energy scaling laws (e.g., $k^{-3/2}$ vs. $k^{-2}$).
• The Challenge: Supersonic flows involve complex interactions between shocks, acoustic modes, and nonlinear advection. The transition from subsonic to supersonic regimes fundamentally alters the flow organization, introducing significant compressive (dilatational) components and pressure-dilatation effects that are not present in incompressible flows.

2. Methodology

The authors performed Direct Numerical Simulations (DNS) of forced compressible turbulence to resolve the full range of scales from the energy injection scale down to the dissipation scale.

• Numerical Solver: They utilized DHARA, a high-performance, GPU-accelerated Python solver.
  • Scheme: A seventh-order, low-dissipation Targeted Essentially Non-Oscillatory (TENO7) scheme was employed. This was chosen over standard WENO schemes to minimize numerical dissipation in smooth regions while maintaining robustness near shocks.
  • Hardware: Simulations were run on the Polaris supercomputer (128 NVIDIA A100 GPUs), achieving a $\sim$150x speedup compared to single-core CPU runs.
• Simulation Parameters:
  • Domain: Periodic cubic box of size $(2\pi)^3$ with a resolution of $1024^3$ grid points.
  • Mach Numbers ($M_t$): Six cases were simulated ranging from low subsonic to high supersonic: $M_t = \{0.2, 0.7, 1.0, 1.4, 1.8, 3.0\}$.
  • Forcing: Stochastic Ornstein–Uhlenbeck (OU) forcing was applied in Fourier space. The forcing ratio between solenoidal (rotational) and compressive modes was controlled via a parameter $\zeta$ (mostly solenoidal $\zeta=2/3$, some mixed $\zeta=1/3$).
  • Physics: The flow was modeled as nearly isothermal ($\gamma \approx 1$) to simulate astrophysical conditions with efficient radiative cooling.
  • Duration: Simulations ran for $>40$ eddy turnover times to ensure statistical convergence.

3. Key Contributions

1. High-Fidelity DNS: Provided the first high-resolution ($1024^3$) DNS of forced compressible turbulence across a wide $M_t$ range using a low-dissipation TENO7 scheme, avoiding the numerical artifacts of previous studies.
2. Mode-to-Mode Analysis: Applied a rigorous mode-to-mode transfer formalism to separately quantify energy fluxes for rotational (solenoidal) and compressive modes, as well as cross-mode transfers.
3. Burgers Turbulence Connection: Established a quantitative link between the compressive modes in supersonic turbulence and classical Burgers turbulence.

4. Key Results

A. Flow Structures

• As $M_t$ increases, the flow transitions from fine-scale vortical structures (subsonic) to intense, persistent shock sheets and localized compression regions.
• Vorticity generation becomes increasingly concentrated near shock fronts due to baroclinic torque, creating a coupling between compressive seeds and solenoidal cascades.

B. Energy Spectra Scaling

The study reveals a fundamental shift in scaling laws as the flow becomes supersonic:

• Rotational (Solenoidal) Spectrum ($E_R$):
  • Subsonic ($M_t < 1$): Follows the Kolmogorov $k^{-5/3}$ scaling.
  • Supersonic ($M_t > 1$): Steepens significantly, approaching a Burgers-like $k^{-2}$ scaling at $M_t=3.0$. This steepening is driven by energy draining from rotational modes to compressive modes.
• Compressive Spectrum ($E_C$):
  • Subsonic: Follows a $k^{-2}$ scaling (due to shocklets).
  • Supersonic: Becomes shallower than $k^{-2}$ (deviating from standard Burgers scaling) as $M_t$ increases.
• Density Spectrum: Flattens with increasing $M_t$, approaching $k^{-1}$ at $M_t=3.0$.
• Total Kinetic Energy: Closely follows the scaling of the rotational component.

C. Energy Fluxes and Transfers

• Cross-Scale Transfer: In supersonic regimes, there is a dominant cross-scale transfer of energy from solenoidal to compressive modes within the inertial range.
  • This cross-flux ($\tilde{\Pi}_{R \to C}$) increases with $M_t$, reaching $\sim 0.70$ of the total injection at $M_t=3.0$.
  • This energy drain reduces the forward cascade available for rotational modes, causing the rotational flux to decrease with wavenumber (breaking the constant flux assumption of Kolmogorov theory).
• Pressure Dilatation: Unlike in subsonic flows where pressure work is confined to large scales, in supersonic flows, pressure dilatation (converting kinetic to internal energy) remains non-negligible throughout the inertial range.
• Dissipation: As $M_t$ increases, normalized rotational dissipation decreases, while compressive dissipation and pressure dilatation increase.

D. Connection to Burgers Turbulence

The authors demonstrated that the compressive components in supersonic turbulence obey scaling laws remarkably similar to classical Burgers turbulence:

• RMS Compressive Velocity ($U_C$): Scales as $U_C \approx \Delta V / \sqrt{12}$.
• Compressive Energy Flux ($\Pi_C$): Scales as $\Pi_C \approx (\Delta V)^3 / (12L)$.
• Here, $\Delta V$ represents the velocity jump across shocks. This suggests that while the total flow is complex, the compressive dynamics are governed by shock-dominated physics similar to the Burgers equation.

5. Significance

• Resolution of Controversy: The results clarify previous discrepancies in the literature. The observed $k^{-2}$ scaling for rotational modes in supersonic turbulence (contrary to the $k^{-3/2}$ reported in some earlier studies) is attributed to the high resolution and low-dissipation schemes used here, which correctly capture the energy drain to compressive modes.
• Mechanism Elucidation: The paper establishes that supersonic turbulence cannot be described by a simple Kolmogorov or Burgers cascade alone. Instead, it is a hybrid system where intermodal energy conversion (solenoidal $\to$ compressive) and pressure-dilatation fundamentally reshape the energy cascade.
• Astrophysical Modeling: The derived scaling laws for compressive modes provide a robust theoretical basis for modeling energy dynamics in extreme astrophysical environments (e.g., interstellar medium, molecular clouds) where compressibility is the dominant factor.
• Computational Benchmark: The study validates the TENO7 scheme and the DHARA solver as powerful tools for high-fidelity compressible turbulence research, offering a path toward predictive theories for supersonic flows in both astrophysics and engineering.