Numerical simulations presented in this work were performed on a single-processor workstation (2.5GHz, 1Gb RAM). We performed a direct numerical simulation, integrating the dynamical equations of motion (4) and (5) using pseudo-spectral method with resolution of wavenumbers. Numerical integrator used for advancing in time was RK7(8) presented in . Time step was where is the period of the shortest wave on the axis. Approximate processor time for this work was 4.5 weeks.
In our numerical experiment, we force the system in the -space ring with and . This ring is located at the low wavenumber part of the -space in order to generate energy cascade toward large 's, but we deliberately avoid forcing even longer waves () because our experience shows that this would lead to undesirable strong anisotropic effects. In the ring, we fix the shape to coincide with the ZF spectrum, , and hence we set , . These fixed amplitudes were then multiplied by random phase factors. Thus, surface and velocity potential were set to , where were uniformly distributed in and
Nonlinearity parameter was set to , which is a sufficient value to produce a resonance broadening for supporting energy cascade.