DIRECT NUMERICAL-SIMULATION OF TRANSITION TO TURBULENCE FROM A HIGH-SYMMETRY INITIAL CONDITION

The three-dimensional (3-D) time evolution of a high-symmetry initial condition [J. Phys. Soc. Jpn. 54, 2132 (1985)] is simulated using a Fo urier pseudospectral method for Re = 1/nu = 500, 1000, 2000, and 5000 with an effective resolution of 1024(3) collocation points (171(3) ind ependent modes, maximum wave number k(max) = 340). It is found that mu ch before the peak enstrophy is reached, there is a short interval whe n the local quantities increase sharply. It is also found that during this interval, six vortex dipoles (at the origin) and three dipoles (a t the pi/2 comer) collapse toward two separate vorticity null points a t the opposite corners of the domain in a nearly self-similar fashion. The coherent vortices break up afterward, followed by a sharp decreas e in local quantities. The singularity analysis shows that, within the limits of the resolution, the maximum vorticity scales approximately as (T-T(c))-1, shortly before the breakup. However, the increase in pe ak vorticity stops at a certain time, possibly due to viscous dissipat ion effects. The temporal evolution of the width of the analyticity st rip shows that delta approaches zero at a rate faster than exponential , but reaches a minimum value and starts to increase. This suggests th at the solution remains uniformly analytic, as is the case in the visc ous Burgers equation.