Wave turbulence in a rotating channel: numerical implementation and results

The analysis of Scott ( J. Fluid Mech. , vol. 741, 2014, pp. 316–349) is implemented numerically. Decaying turbulence is confined to a channel between two infinite, parallel, rotating walls. The Rossby and Ekman numbers are supposed small, the former condition making nonlinearity small, while the latter allows the turbulence to persist for the many rotational periods needed for the small nonlinearity to be effective. The flow is expressed as a combination of inertial waveguide modes, indexed by a two-dimensional wave vector $\boldsymbol{k}$ and an integer n . The $n = 0$ modes form a two-dimensional component of the flow, whereas the remainder is the wave component, on which attention is focused in this article. Assuming statistical axisymmetry and homogeneity in directions parallel to the walls, the second-order moments of the mode amplitudes yield a spectral matrix ${A_{nm}}(k,t)$ (where $k = |\boldsymbol{k} |$ ), of which the diagonal elements describe the distribution of energy over different modes. Wave-turbulence analysis provides an equation governing the time evolution of ${A_{nn}}$ , $n \ne 0$ , the wave spectra, which forms the basis for the present work. The initial distribution of energy is Gaussian and depends on a parameter $\varXi$ , the initial spectral width. The problem has two other parameters, ${\beta _w}$ and ${\beta _v}$ , which correspond to two distinct viscous dissipative mechanisms: wall damping due to boundary layers and volumetric damping by viscous effects throughout the flow. Results obtained by numerical solution include the time evolution of the total wave energy, E , and the detailed description of its distribution over k and n provided by ${A_{nn}}(k)$ .