Freely decaying turbulence in a finite domain at finite Reynolds number

We perform direct numerical simulations to study the effects of the finite Reynolds number and domain size on the decay law of Saffman turbulence. We observe that the invariant for Saffman turbulence, u(2)l(3), and non-dimensional dissipation coefficient, C-E = E/(u(3)/l), are sensitive to finite domain size; here, u is the rms velocity, l is the integral length scale, and E is the energy dissipation rate. Consequently, the exponent n in the decay law u(2) similar to t(-n) for Saffman turbulence deviates from 6/5. Due to the finite Reynolds number and the domain size, Saffman turbulence decays at a faster rate (i.e., n > 6/5). However, the exponent n = 6/5 is more sensitive to the domain size than to the Reynolds number. From the simulations, we find that n remains close to 6/5 as long as R-lambda greater than or similar to 10 and l less than or similar to 0.3L(box); here, R-lambda is the Reynolds number based on the Taylor microscale and L-box is the domain size. We also notice that n becomes slightly lower than 6/5 for a part of the decay period. Interestingly, this trend n < 6/5 is also observed earlier in freely decaying grid-generated turbulence.