A new algorithm for simulating homogeneous decaying turbulence in an incompressible fluid subjected to uniform system rotation has been proposed. The complex helical waves decomposition is introduced to diagonalize the Coriolis term of the Navier Stokes equations in the Fourier space. The integral factor technique has been extended to the Coriolis term. Since the Coriolis term is integrated analytically in the algorithm, time increment is not restricted by the system rotation. It is proved that the new algorithm is effective at high rotation rate. The effect of system rotation is observed clearly in high rotation cases. It is also indicated that the energy spectrum depends on the modulus of wave vector κ and cos θ, where θ is the angle between the rotation vector Ω and the wave vector κ.
Read full abstract