Two-dimensional `turbulent' ยฏow around a rectangular cylinder has been simulated at Re 10,000 using a sixthorder-accurate ยฎnite volume method for the discretization of convection and diffusion. The spatial discretization consists of a combination of a seventh-order upwind-biased method for the convective terms and an eighth-order central method for the diffusive terms, discretized on a stretched and staggered grid. To cope with the stretching of the grid, Lagrange interpolations are used.
The method applied to obtain a boundary condition for the velocity in the x-direction at the outยฏow boundary is shown not to affect the ยฏow in the interior of the computational domain in a way that is visible in various snapshots of the vorticity ยฎeld. The variation in the velocity in the x-direction with time is itself found to be relatively small near the outยฏow boundary.
Several turbulence statistics have been gathered from a simulation of the ยฏow developed during 77 dimensionless time units. Snapshots of the vorticity ยฎeld of the developed ยฏow show the presence of a vortexstreet-like structure. Typical 2D turbulent behaviour, such as the appearance of monopolar, dipolar and tripolar vortices due to the amalgamation of vorticity in the wake and the x ร1a2 scaling of the velocity defect in the wake, has been obtained.