Large-scale turbulence in fluid layers and other quasi-two-dimensional compressible systems consists of planar vortices and waves. Separately, wave turbulence usually produces a direct energy cascade, while solenoidal planar turbulence transports energy to large scales by an inverse cascade. Here, we consider turbulence at finite Mach numbers when the interaction between acoustic waves and vortices is substantial. We employ solenoidal pumping at intermediate scales and show how both direct and inverse energy cascades are formed starting from the pumping scale. We show that there is an inverse cascade of kinetic energy up to a scale $\ell$, where a typical velocity reaches the speed of sound; this creates shock waves, which provide for a compensating direct cascade. When the system size is less than $\ell$, the steady state contains a system-size pair of long-living condensate vortices connected by a system of shocks. Thus turbulence in fluid layers processes energy via a loop: Most energy first goes to large scales via vortices and is then transported by waves to small-scale dissipation.

An inverse cascade - energy transfer to progressively larger scales - is a salient feature of two-dimensional turbulence. If the cascade reaches the system scale, it creates a coherent flow expected to have the largest available scale and conform with the symmetries of the domain. In a doubly periodic rectangle, the mean flow with zero total momentum was therefore believed to be unidirectional, with two jets along the short side; while for an aspect ratio close to unity, a vortex dipole was expected. Using direct numerical simulations, we show that in fact neither the box symmetry is respected nor the largest scale is realized: the flow is never purely unidirectional since the inverse cascade produces coherent vortices, whose number and relative motion are determined by the aspect ratio. This spontaneous symmetry breaking is closely related to the hierarchy of averaging times. Long-time averaging restores translational invariance due to vortex wandering along one direction, and gives jets whose profile, however, can be deduced neither from the largest-available-scale argument, nor from the often employed maximum-entropy principle or quasi-linear approximation.

Using data from a large-scale three-dimensional simulation of supersonic isothermal turbulence, we have tested the validity of an exact flux relation derived analytically from the Navier--Stokes equation by Falkovich, Fouxon and Oz [2010 New relations for correlation functions in Navier--Stokes turbulence. J. Fluid Mech. 644, 465]. That relation, for compressible barotropic fluids, was derived assuming turbulence generated by a large-scale force. However, compressible turbulence in simulations is usually initialized and maintained by a large-scale acceleration, as in gravity-driven astrophysical flows. We present a new approximate flux relation for isothermal turbulence driven by a large-scale acceleration, and find it in reasonable agreement with the simulation results.