The Energy Cascade Rate in Supersonic Magnetohydrodynamic Turbulence

Three-dimensional direct numerical simulations (DNS) are implemented to investigate the energy cascade rate in compressible isothermal magnetohydrodynamic (MHD) turbulence. Utilizing an exact law derived from the K\'arm\'an-Howarth equation, we examine the contributions of flux and non-flux terms to the cascade rate across a broad range of sonic and Alfv\'enic Mach numbers, from subsonic to supersonic regimes and varying mean magnetic fields. Cascade rates are computed using on-grid 3-D decomposition and two plasma increment approaches: signed and absolute values. Anisotropy induced by strong magnetic fields is analyzed through angular-dependent scaling of the cascade terms. Moreover, the increment calculation method significantly influences the relative contributions of flux and non-flux terms, with absolute methods tending to overestimate the latter. These findings extend current studies of compressible turbulence and offer critical insights into energy transfer mechanisms relevant to many astrophysical phenomena.