pith. sign in

arxiv: 2604.26290 · v1 · submitted 2026-04-29 · ⚛️ physics.flu-dyn · physics.comp-ph

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

Pith reviewed 2026-05-07 12:43 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn physics.comp-ph
keywords supersonic turbulenceenergy spectracompressible turbulencedirect numerical simulationMach numberenergy cascaderotational compressive modesturbulent dissipation
0
0 comments X

The pith

In supersonic turbulence, rotational energy spectra steepen toward minus-two scaling with rising Mach number while compressive spectra become shallower.

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

This paper presents high-resolution computer simulations of turbulent flows moving faster than sound to track how energy moves from large scales to small ones. The central finding is that higher flow speeds relative to sound speed change the power-law shape of energy spectra differently for two kinds of motion: swirling motions develop a steeper spectrum, while squeezing and expanding motions develop a shallower one. These shifts occur because energy moves across scales from swirling components to compressing ones and because pressure forces perform significant work on the fluid. The results matter for predicting behavior in contexts such as star-forming regions and high-speed vehicle flows, where accurate cascade descriptions affect mixing, heating, and shock formation.

Core claim

High-fidelity direct numerical simulations at 1024 cubed resolution of forced compressible turbulence at turbulent Mach numbers from 0.2 to 3.0 demonstrate that the rotational kinetic energy spectrum steepens from a minus five-thirds power law toward a minus two power law as the Mach number increases. The compressive energy spectrum becomes shallower, deviating from expected minus two scaling. These modifications result from a dominant cross-scale transfer of energy from swirling to compressing flow components in the inertial range, together with significant pressure work contributions. Scaling laws for the root-mean-square compressive velocity and the compressive energy flux follow those of

What carries the argument

The dominant cross-scale transfer of energy from swirling to compressing flow components within the inertial range, together with pressure forces performing work on the fluid.

If this is right

  • Rotational dissipation rates decrease as the turbulent Mach number rises.
  • Compressive dissipation rates and pressure work increase with rising Mach number.
  • Overall energy injection rates depend on the type of forcing applied rather than on Mach number.
  • Root-mean-square compressive velocity and compressive energy flux obey the same scaling relations found in classical shock-dominated turbulence.

Where Pith is reading between the lines

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

  • These mode-specific transfers suggest that large-scale models of sustained supersonic turbulence should track rotational and compressive contributions separately to capture the correct dissipation balance.
  • The Mach-dependent shift may alter mixing and heating rates in other high-speed compressible systems beyond the regimes simulated here.
  • Similar inter-component energy movement could appear in decaying turbulence or in flows with different forcing, providing a testable extension of the cascade picture.

Load-bearing premise

The computer simulations must resolve the smallest eddies and sharp discontinuities without the solution method introducing artificial changes to the measured energy distributions.

What would settle it

An independent higher-resolution simulation or laboratory measurement at high Mach number that finds the rotational energy spectrum staying near the minus five-thirds law instead of approaching minus two would disprove the reported change in scaling.

Figures

Figures reproduced from arXiv: 2604.26290 by Dhananjay Singh, Harshit Tiwari, Mahendra K. Verma, Rajesh Ranjan.

Figure 1
Figure 1. Figure 1: FIG. 1. Time evolution of (a) total kinetic energy view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Flow structures in subsonic and transonic turbulence at view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Flow structures in supersonic turbulence at view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. For view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. For view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. For forced compressible turbulence with (a) view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Plots of normalized energy transfer terms with turbulent Mach number view at source ↗
Figure 8
Figure 8. Figure 8: (b) presents the rotational and compressive RMS velocities (UR, UC ) and the maximum energy fluxes (ΠR, ΠC ). The compressive components align re￾markably well with Burgers equation predictions (solid and dashed black lines). A best-fit analysis yields UC = 0.54(∆V ) 1.08 ≈ ∆V √ 12 , (44) ΠC = 0.038 (∆V ) 3.10 ≈ (∆V ) 3 12L . (45) Interestingly, the rotational component exhibits similar scaling trends with… view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. Scalability of view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10. Case study using supersonic Taylor–Green vortex simulation. (a) Time evolution of the volume-averaged kinetic view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11. Normalized pressure error along the view at source ↗
Figure 12
Figure 12. Figure 12: FIG. 12. Kelvin-Helmholtz instability diagnostics. (a), (b): Time evolution of volume-integrated kinetic energy view at source ↗
read the original abstract

Supersonic turbulence is vital to astrophysical and high-speed engineering flows, yet its energy transfer mechanisms remain poorly understood. We present high-resolution ($1024^3$) direct numerical simulations (DNS) of forced compressible turbulence across a range of turbulent Mach numbers ($M_t = 0.2$ to $3.0$). Using the GPU-accelerated solver \texttt{DHARA} with a seventh-order, low-dissipation Targeted Essentially Non-Oscillatory (TENO) scheme, we resolve both fine-scale eddies and sharp shock fronts. Our results reveal a fundamental shift in the energy cascade in the supersonic regime. As $M_t$ increases, the rotational kinetic energy spectrum steepens from a Kolmogorov-like $k^{-5/3}$ scaling toward a Burgers-like $k^{-2}$ scaling. Conversely, the compressive energy spectrum becomes shallower, deviating from Burgers scaling. We show that these spectral modifications are driven by a dominant cross-scale transfer of energy from solenoidal to compressive modes within the inertial range, alongside significant contributions from pressure dilatation. Scaling laws for the root-mean-square compressive velocity ($U_C$) and compressive energy flux ($\Pi_C$) are found to mirror classical Burgers turbulence. Finally, we show that while energy injection rates depend on forcing type rather than Mach number, increased $M_t$ leads to decreased rotational dissipation and increased compressive dissipation and pressure dilatation. These findings elucidate intermodal energy cascade mechanisms, advancing our understanding of energy transfers in supersonic turbulence.

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

2 major / 1 minor

Summary. The manuscript reports 1024^3 DNS of forced compressible turbulence at Mt = 0.2–3.0 using the GPU-accelerated DHARA solver with a seventh-order low-dissipation TENO scheme. It claims that rotational kinetic energy spectra steepen from Kolmogorov-like k^{-5/3} toward Burgers-like k^{-2} as Mt increases, while compressive spectra become shallower; these changes are attributed to dominant cross-scale solenoidal-to-compressive energy transfer and pressure dilatation within the inertial range. Additional results include Burgers-like scaling for rms compressive velocity U_C and compressive flux Π_C, Mach-dependent shifts in dissipation (decreased rotational, increased compressive and pressure dilatation), and forcing-type dependence of injection rates.

Significance. If the reported spectral modifications and intermodal transfers are shown to be free of numerical artifacts, the work would provide valuable quantitative insight into energy cascade mechanisms in supersonic turbulence, relevant to astrophysical and high-speed engineering flows. The high-resolution DNS approach with explicit decomposition into rotational/compressive modes and direct computation of cross-scale fluxes represents a clear strength.

major comments (2)
  1. [Abstract and Methods] Abstract and Methods: The claim that the 1024^3 TENO7 simulations 'fully resolve both fine-scale eddies and sharp shock fronts' is not supported by any grid-convergence tests, resolution-sensitivity studies, or explicit checks that inertial-range spectral exponents and flux plateaus remain unchanged under increased resolution or higher scheme order. At Mt=3 the compressive discontinuities are thin; without such evidence the reported steepening of rotational spectra toward k^{-2} cannot be unambiguously attributed to physics rather than residual numerical dissipation.
  2. [Results (spectral analysis)] Results (spectral analysis): No error bars, uncertainty estimates, or details on the selection of inertial-range fitting windows are provided for the reported spectral exponents. This omission makes it impossible to assess the statistical significance of the claimed transition from k^{-5/3} to k^{-2} for rotational modes or the shallowing of compressive spectra.
minor comments (1)
  1. [Abstract] The abstract states that 'scaling laws for U_C and Π_C mirror classical Burgers turbulence' but does not quote the explicit functional forms or the range of Mt over which they hold; these should be stated explicitly in the abstract or early results section.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their constructive and detailed review of our manuscript. We address each major comment below and have revised the manuscript accordingly to strengthen the presentation of our results.

read point-by-point responses
  1. Referee: [Abstract and Methods] Abstract and Methods: The claim that the 1024^3 TENO7 simulations 'fully resolve both fine-scale eddies and sharp shock fronts' is not supported by any grid-convergence tests, resolution-sensitivity studies, or explicit checks that inertial-range spectral exponents and flux plateaus remain unchanged under increased resolution or higher scheme order. At Mt=3 the compressive discontinuities are thin; without such evidence the reported steepening of rotational spectra toward k^{-2} cannot be unambiguously attributed to physics rather than residual numerical dissipation.

    Authors: We agree that explicit resolution-sensitivity analysis would strengthen the attribution of the observed spectral changes to physical mechanisms. The 1024^3 resolution and TENO7 scheme were chosen following standard practices for DNS of compressible turbulence at these Mach numbers to capture both eddies and shocks, but the submitted manuscript does not include direct comparisons. In the revised version we will add a dedicated subsection (in Methods or an appendix) presenting spectra and flux plateaus computed at 512^3 and 1024^3 resolutions. This will demonstrate that the inertial-range exponents and the reported steepening of rotational spectra remain unchanged, thereby supporting that the trends are not dominated by residual numerical dissipation. revision: yes

  2. Referee: [Results (spectral analysis)] Results (spectral analysis): No error bars, uncertainty estimates, or details on the selection of inertial-range fitting windows are provided for the reported spectral exponents. This omission makes it impossible to assess the statistical significance of the claimed transition from k^{-5/3} to k^{-2} for rotational modes or the shallowing of compressive spectra.

    Authors: We concur that quantitative uncertainty measures and transparent fitting criteria are necessary for assessing the robustness of the Mach-dependent spectral transitions. In the revised manuscript we will report error bars on all fitted exponents, obtained from the standard deviation across multiple independent time windows in the statistically stationary regime. We will also explicitly state the wavenumber intervals chosen for each fit (identified where the energy flux is approximately constant) and include a brief sensitivity test showing how the exponents vary with modest changes in the fitting window. These additions will allow readers to evaluate the statistical significance of the transition from Kolmogorov-like to Burgers-like scaling. revision: yes

Circularity Check

0 steps flagged

No significant circularity; results emerge from DNS and post-processing

full rationale

The paper's central claims on spectral steepening, intermodal energy transfer, and Burgers-like scalings for compressive quantities are obtained by integrating the compressible Navier-Stokes equations at 1024^3 resolution with the DHARA TENO7 solver, then computing rotational/compressive decompositions, spectra, and fluxes directly from the velocity and pressure fields. No step defines a quantity in terms of a result it is later used to predict, no fitted parameters are renamed as predictions, and no load-bearing uniqueness theorems or ansatzes are imported via self-citation. The observed modifications to the cascade are therefore independent outputs of the numerical experiment rather than tautological restatements of its inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claims rest on the compressible Navier-Stokes equations and the numerical properties of the TENO scheme; no new physical entities are introduced and no parameters are fitted to produce the target spectra.

axioms (2)
  • standard math Compressible Navier-Stokes equations govern the evolution of density, velocity, and pressure in the simulated flows
    Invoked implicitly as the basis for all DNS of compressible turbulence
  • domain assumption The seventh-order TENO scheme provides sufficient low-dissipation resolution of both turbulence and shocks at the employed grid resolution
    Required for the claim that observed spectral modifications are physical rather than numerical

pith-pipeline@v0.9.0 · 5587 in / 1528 out tokens · 61485 ms · 2026-05-07T12:43:56.112114+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

83 extracted references · 83 canonical work pages

  1. [1]

    This extension makes it possible to calculate fluxes sep- arately for rotational and compressive motions and to directly quantify cross-mode energy transfers

    and Verma [16] developed the mode-to-mode trans- fer formalism for incompressible turbulence, which was later extended to compressible flows by Singhet al.[43]. This extension makes it possible to calculate fluxes sep- arately for rotational and compressive motions and to directly quantify cross-mode energy transfers. The above discussion highlights that ...

  2. [2]

    Non-dimensionalization The equations have been nondimensionalized using ref- erence quantities: densityρ 0, temperatureT 0, velocity u0, and lengthl0. The dimensionless numbers, thus, gov- erning the system are: Reynolds number Re 0 = ρ0u0l0 µ ,(7) Mach numberM 0 = u0 c = u0√γR∗T0 ,(8) Prandtl number Pr = µCp Kc (9) wherecis the speed of sound,R ∗ is the ...

  3. [3]

    External forcing We apply an external stochastic forcing at large scales to sustain turbulence, following a generalised Orn- stein–Uhlenbeck (OU) process in Fourier space similar to Schmidtet al.[10]. The evolution of the Fourier- transformed force ˆF(k, t)is governed by the stochastic differential equation: dˆF(k, t) =g ζ " − ˆF(k, t)dt t0 +P ζ(k)·[f 0(k...

  4. [4]

    We note, however, that those studies defined the kinetic energy asw2/2, where w= √ρuis the density-weighted velocity

    and Federrathet al.[35]. We note, however, that those studies defined the kinetic energy asw2/2, where w= √ρuis the density-weighted velocity. The density spectrum also flattens with increasingMt, approaching k−1 atM t = 3.0, which agrees with high Mach number results in the literature [9, 35]. In all cases, the total kinetic energy remains close to the s...

  5. [5]

    The domain is[−π, π] 3 with periodic boundary conditions in all directions

    Supersonic Taylor-Green vortex The 3D Taylor-Green vortex (TGV) problem [69] is considered in the supersonic regime [68, 70] without ex- ternal forcing. The domain is[−π, π] 3 with periodic boundary conditions in all directions. We considerM0 = 1.25,γ= 1.4,Re 0 = 1600,Pr = 0.71and grid size of 2563 for all runs. We track the time evolution of the volume-a...

  6. [6]

    We per- form the isentropic vortex test forγ= 1.4on a square domainx, y∈[−6,6], with periodic boundary conditions in both directions

    Isentropic Vortex We consider the well-known isentropic vortex problem in two dimensions [71, 72], a standard benchmark for assessing the accuracy and dissipation characteristics of numerical schemes solving the Euler equations. We per- form the isentropic vortex test forγ= 1.4on a square domainx, y∈[−6,6], with periodic boundary conditions in both direct...

  7. [7]

    The initial condition consists of a shear layer in the horizon- 17 FIG

    Kelvin-Helmholtz Instability We solve the two-dimensional Euler equations to study the classical Kelvin-Helmholtz (KH) instability [73]. The initial condition consists of a shear layer in the horizon- 17 FIG. 10. Case study using supersonic Taylor–Green vortex simulation. (a) Time evolution of the volume-averaged kinetic energyE u(t)for different reconstr...

  8. [8]

    B. W. Carroll,An Introduction to Modern Astrophysics, 2nd ed. (Pearson, 2006)

  9. [9]

    B. G. Elmegreen and J. Scalo, Interstellar turbulence I: Observations and processes, Annu. Rev. Astron. Astro- phys.42, 211 (2004)

  10. [10]

    Mac Low and R

    M.-M. Mac Low and R. S. Klessen, Control of star forma- tionbysupersonicturbulence,Reviewsofmodernphysics 76, 125 (2004)

  11. [11]

    von Kármán, Compressibility effects in aerodynamics, Journal of the Aeronautical Sciences8, 337 (1941)

    T. von Kármán, Compressibility effects in aerodynamics, Journal of the Aeronautical Sciences8, 337 (1941)

  12. [12]

    M. R. Krumholz and C. F. McKee, A general theory of turbulence-regulated star formation, from spirals to ultraluminous infrared galaxies, Astrophys. J.630, 250 (2005). 18 FIG. 12. Kelvin-Helmholtz instability diagnostics. (a), (b): Time evolution of volume-integrated kinetic energyE(t)for resolutions512 2 and1024 2, respectively. (c), (d): Velocity power ...

  13. [13]

    Padoan and Å

    P. Padoan and Å. Nordlund, The star formation rate of supersonic magnetohydrodynamic turbulence, Astro- phys. J.730, 40 (2011)

  14. [14]

    Fujisawa, Review of plasma turbulence experiments, Proc

    A. Fujisawa, Review of plasma turbulence experiments, Proc. Jpn. Acad., Ser. B97, 103 (2021)

  15. [15]

    Ingenito and C

    A. Ingenito and C. Bruno, Physics and regimes of super- sonic combustion, AIAA journal48, 515 (2010)

  16. [16]

    A. G. Kritsuk, M. L. Norman, P. Padoan, and R. Wag- ner, The statistics of supersonic isothermal turbulence, Astrophys. J.665, 416 (2007)

  17. [17]

    Schmidt, C

    W. Schmidt, C. Federrath, M. Hupp, S. Kern, and J. C. Niemeyer, Numerical simulations of compressively driven interstellar turbulence-i. isothermal gas, Astron. Astro- phys.494, 127 (2009)

  18. [18]

    Federrath, On the universality of supersonic turbu- lence, Mon

    C. Federrath, On the universality of supersonic turbu- lence, Mon. Not. R. Astron. Soc.436, 1245 (2013)

  19. [19]

    A. N. Kolmogorov, Dissipation of Energy in Locally Isotropic Turbulence, Dokl Acad Nauk SSSR32, 16 (1941)

  20. [20]

    A. N. Kolmogorov, The local structure of turbulence in incompressible viscous fluid for very large Reynolds num- bers, Dokl Acad Nauk SSSR30, 301 (1941)

  21. [21]

    R. H. Kraichnan, The structure of isotropic turbulence at very high Reynolds numbers, J. Fluid Mech.5, 497 (1959)

  22. [22]

    G.Dar, M.K.Verma,andV.Eswaran,Energytransferin two-dimensional magnetohydrodynamic turbulence: for- malismandnumericalresults,PhysicaD157,207(2001)

  23. [23]

    M. K. Verma, Statistical theory of magnetohydrody- namic turbulence: recent results, Phys. Rep.401, 229 (2004)

  24. [24]

    Lesieur,Turbulence in Fluids(Springer-Verlag, Dor- drecht, 2008)

    M. Lesieur,Turbulence in Fluids(Springer-Verlag, Dor- drecht, 2008)

  25. [25]

    Fluid Mech.47, 525 (1971)

    R.H.Kraichnan,Inertial-rangetransferintwo-andthree- dimensional turbulence, J. Fluid Mech.47, 525 (1971)

  26. [26]

    J. A. Domaradzki and R. S. Rogallo, Local Energy Trans- fer and Nonlocal Interactions in Homogeneous, Isotropic Turbulence, Phys. Fluids A2, 414 (1990)

  27. [27]

    Zhou, Degrees of locality of energy transfer in the inertial range, Phys

    Y. Zhou, Degrees of locality of energy transfer in the inertial range, Phys. Fluids5, 1092 (1993)

  28. [28]

    Zhou, Interacting scales and energy transfer in isotropic turbulence, Phys

    Y. Zhou, Interacting scales and energy transfer in isotropic turbulence, Phys. Fluids A5, 2511 (1993)

  29. [29]

    M. K. Verma, A. Ayyer, O. Debliquy, S. Kumar, and A. V. Chandra, Local shell-to-shell energy transfer via nonlocal interactions in fluid turbulence, Pramana-J. Phys.65, 297 (2005)

  30. [30]

    Kida and S

    S. Kida and S. A. Orszag, Energy and spectral dynamics in forced compressible turbulence, J. Sci. Comput.5, 85 (1990)

  31. [31]

    S. K. Lele, Compressibility effects on turbulence, Annu. Rev. Fluid Mech.26, 211 (1994)

  32. [32]

    Erlebacher, M

    G. Erlebacher, M. Y. Hussaini, H.-O. Kreiss, and S. Sarkar, The analysis and simulation of compressible turbulence, Theoret. Comput. Fluid Dyn.2, 73 (1990)

  33. [33]

    J. P. John, D. A. Donzis, and K. R. Sreenivasan, Does dissipative anomaly hold for compressible turbulence?, Journal of Fluid Mechanics920, A20 (2021)

  34. [34]

    Sarkar, The pressure–dilatation correlation in com- pressible flows, Physics of Fluids A: Fluid Dynamics4, 2674 (1992)

    S. Sarkar, The pressure–dilatation correlation in com- pressible flows, Physics of Fluids A: Fluid Dynamics4, 2674 (1992)

  35. [35]

    Miura and S

    H. Miura and S. Kida, Acoustic energy exchange in com- pressible turbulence, Phys. Fluids7, 1732 (1995). 19

  36. [36]

    D. S. Praturi and S. S. Girimaji, Effect of pressure- dilatation on energy spectrum evolution in compressible turbulence, Phys. Fluids31(2019)

  37. [37]

    Jagannathan and D

    S. Jagannathan and D. A. Donzis, Reynolds and Mach numberscalinginsolenoidally-forcedcompressibleturbu- lence using high-resolution direct numerical simulations, J. Fluid Mech.789, 669 (2016)

  38. [38]

    I - hydrodynamics, turbulence, and waves, Phys

    G.P.ZankandW.H.Matthaeus,Theequationsofnearly incompressible fluids. I - hydrodynamics, turbulence, and waves, Phys. Fluids A3, 69 (1991)

  39. [39]

    J.Wang, Y.Yang, Y.Shi, Z.Xiao, X.T.He,andS.Chen, Cascade of Kinetic Energy in Three-Dimensional Com- pressible Turbulence, Phys. Rev. Lett.110, 214505 (2013)

  40. [40]

    Schmidt and P

    W. Schmidt and P. Grete, Kinetic and internal energy transfer in implicit large-eddy simulations of forced com- pressible turbulence, Phys. Rev. E100, 043116 (2019)

  41. [41]

    Sakurai and T

    Y. Sakurai and T. Ishihara, Direct numerical simulations of compressible turbulence in a periodic box: Effect of isothermal assumptions on turbulence statistics, Phys. Fluids36, 085152 (2024)

  42. [42]

    Federrath, J

    C. Federrath, J. Roman-Duval, R. S. Klessen, W. Schmidt, and M.-M. Mac Low, Comparing the statistics of interstellar turbulence in simulations and observations-Solenoidal versus compressive turbulence forcing, Astron. & Astrophys.512, A81 (2010)

  43. [43]

    Galtier and S

    S. Galtier and S. Banerjee, Exact relation for correlation functions in compressible isothermal turbulence, Phys. Rev. Lett.107, 134501 (2011)

  44. [44]

    Aluie, Compressible turbulence: the cascade and its locality, Phys

    H. Aluie, Compressible turbulence: the cascade and its locality, Phys. Rev. Lett.106, 174502 (2011)

  45. [45]

    Aluie, Scale decomposition in compressible turbu- lence, Phys

    H. Aluie, Scale decomposition in compressible turbu- lence, Phys. D: Nonlinear Phenom.247, 54 (2013)

  46. [46]

    Aluie, S

    H. Aluie, S. Li, and H. Li, Conservative cascade of ki- netic energy in compressible turbulence, Astrophys. J. 751, L29 (2012)

  47. [47]

    Zhao and H

    D. Zhao and H. Aluie, Inviscid criterion for decomposing scales, Physical Review Fluids3, 054603 (2018)

  48. [48]

    A. G. Kritsuk, R. Wagner, and M. L. Norman, Energy cascade and scaling in supersonic isothermal turbulence, J. Fluid Mech.729, R1 (2013)

  49. [49]

    J. P. Graham, R. Cameron, and M. Schüssler, Turbulent small-scale dynamo action in solar surface simulations, Astrophys. J.714, 1606 (2010)

  50. [50]

    Singh, H

    D. Singh, H. Tiwari, L. Sharma, and M. K. Verma, Math- ematical formulation of mode-to-mode energy transfers and energy fluxes in compressible turbulence, Phys. Rev. Fluids10, 114609 (2025)

  51. [51]

    Colella and P

    P. Colella and P. R. Woodward, The piecewise parabolic method (ppm) for gas-dynamical simulations, Journal of Computational Physics54, 174 (1984)

  52. [52]

    S.Pirozzoli,Numericalmethodsforhigh-speedflows,An- nual Review of Fluid Mechanics43, 163 (2011)

  53. [53]

    L. Fu, X. Y. Hu, and N. A. Adams, A family of high-order targeted eno schemes for compressible-fluid simulations, J. Comput. Phys.305, 333 (2016)

  54. [54]

    Tiwari, L

    H. Tiwari, L. Sharma, and M. K. Verma, On the absence of the ultimate regime in turbulent thermal convection, Proceedings of the National Academy of Sciences122, e2513474122 (2025)

  55. [55]

    J. D. Anderson,Computational Fluid Dynamics: The Basics With Applications(McGraw-Hill, New York, 1995)

  56. [56]

    A.KurganovandE.Tadmor,Newhigh-resolutioncentral schemes for nonlinear conservation laws and convection– diffusion equations, J. Comput. Phys.160, 241 (2000)

  57. [57]

    Kurganov and D

    A. Kurganov and D. Levy, A third-order semidiscrete central scheme for conservation laws and convection- diffusion equations, SIAM J. Sci. Comput.22, 1461 (2000)

  58. [58]

    Gottlieb, C.-W

    S. Gottlieb, C.-W. Shu, and E. Tadmor, Strong stability- preserving high-order time discretization methods, SIAM review43, 89 (2001)

  59. [59]

    Jiang and C.-W

    G.-S. Jiang and C.-W. Shu, Efficient implementation of weighted eno schemes, J. Comput. Phys.126, 202 (1996)

  60. [60]

    Borges, M

    R. Borges, M. Carmona, B. Costa, and W. S. Don, An improved weighted essentially non-oscillatory scheme for hyperbolic conservation laws, J. Comput. Phys.227, 3191 (2008)

  61. [61]

    D. Levy, G. Puppo, and G. Russo, Compact central weno schemes for multidimensional conservation laws, SIAM J. Sci. Comput.22, 656 (2000)

  62. [62]

    Okuta, Y

    R. Okuta, Y. Unno, D. Nishino, S. Hido, and C. Loomis, CuPy: A NumPy-Compatible Library for NVIDIA GPU Calculations, inProc. Workshop on Machine Learning Systems (LearningSys) in The Thirty-first Annual Con- ference, NeurIPS(2017). [56]The movies are available here, YouTube (2025)

  63. [63]

    J.Wang, T.Gotoh,andT.Watanabe,Spectraandstatis- tics in compressible isotropic turbulence, Phys. Rev. Flu- ids2, 013403 (2017)

  64. [64]

    J. Wang, M. Wan, S. Chen, C. Xie, and S. Chen, Effect of shock waves on the statistics and scaling in compressible isotropic turbulence, Phys. Rev. E97, 043108 (2018)

  65. [65]

    J. Wang, M. Wan, S. Chen, and S. Chen, Kinetic energy transfer in compressible isotropic turbulence, J. Fluid Mech.841, 581 (2018)

  66. [66]

    M. K. Verma, Variable energy flux in turbulence, J. Phys. A: Math. Theor.55, 013002 (2022)

  67. [67]

    S. A. Orszag, Lectures on the statistical theory of turbu- lenceinfluiddynamics,inLes Houches Summer School of Theoretical Physics, edited by R. Balian and J. L. Peube (Gordon Breach, New York, 1973) p. 235

  68. [68]

    M.K.Verma,Energy transfers in Fluid Flows: Multiscale and Spectral Perspectives(Cambridge University Press, Cambridge, 2019)

  69. [69]

    J. M. Burgers, A Mathematical Model Illustrating the Theory of Turbulence, Advances in Applied Mechanics 1, 171 (1948)

  70. [70]

    P. G. Saffman, inTopics in Nonlinear Physics, edited by N. J. Zabusky (Springer-Verlag, Berlin, 1968) pp. 485– 614

  71. [71]

    S. Alam, P. Sahu, and M. Verma, Universal functions for Burgers turbulence, Phys. Rev. Fluids7, 074605 (2022)

  72. [72]

    C. R. Harris, K. J. Millman, S. J. Van Der Walt, R. Gom- mers, P. Virtanen, D. Cournapeau, E. Wieser, J. Taylor, S. Berg, N. J. Smith,et al., Array programming with NumPy, Nature585, 357 (2020)

  73. [73]

    Dalcín, R

    L. Dalcín, R. Paz, M. Storti, and J. D’Elía, MPI for Python: Performance improvements and MPI-2 exten- sions, J. Parallel Distrib. Comput.68, 655 (2008)

  74. [74]

    Chapelier, D

    J.-B. Chapelier, D. J. Lusher, W. Van Noordt, C. Wen- zel, T. Gibis, P. Mossier, A. Beck, G. Lodato, C. Brehm, M. Ruggeri,et al., Comparison of high-order numeri- cal methodologies for the simulation of the supersonic taylor–green vortex flow, Phys. Fluids.36(2024). 20

  75. [75]

    M. E. Brachet, D. I. Meiron, S. A. Orszag, B. G. Nickel, R. H. Morf, and U. Frisch, Small-scale structure of the taylor–green vortex, J. Fluid Mech.130, 411 (1983)

  76. [76]

    D. J. Lusher and N. D. Sandham, Assessment of low- dissipative shock-capturing schemes for the compressible taylor–green vortex, AIAA Journal59, 533 (2021)

  77. [77]

    Nonomura, N

    T. Nonomura, N. Iizuka, and K. Fujii, Freestream and vortex preservation properties of high-order weno and wcns on curvilinear grids, Comput. Fluids.39, 197 (2010)

  78. [78]

    Song, A.S.Ghate, K.V

    H. Song, A.S.Ghate, K.V. Matsuno, J. R.West, A.Sub- ramaniam, and S. K. Lele, A robust compact finite differ- ence framework for simulations of compressible turbulent flows, J. Comput. Phys.519, 113419 (2024)

  79. [79]

    Zhou,Hydrodynamic Instabilities and Turbulence: Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin– Helmholtz Mixing(Cambridge University Press, Cam- bridge, 2024)

    Y. Zhou,Hydrodynamic Instabilities and Turbulence: Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin– Helmholtz Mixing(Cambridge University Press, Cam- bridge, 2024)

  80. [80]

    San and K

    O. San and K. Kara, Evaluation of riemann flux solvers for weno reconstruction schemes: Kelvin–helmholtz in- stability, Comput. Fluids.117, 24 (2015)

Showing first 80 references.