Numerical simulation of flows past periodic arrays of cylinders

Abstract

We present a detailed numerical investigation of three unsteady incompressible flow problems involving periodic arrays of staggered cylinders. The first problem is a uniperiodic flow with two cylinders in each cell of periodicity. The second problem is a biperiodic flow with two cylinders in each cell, and the last problem is a uniperiodic flow with ten cylinders. Both uniperiodic flows are periodic in the direction perpendicular to the main flow direction. In all three cases, the Reynolds number based on the cylinder diameter is 100, and initially the flow field has local symmetries with respect to the axes of the cylinders parallel to the main flow direction. Later on, these symmetries break, vortex shedding is initiated, and gradually the scale of the shedding increases until a temporally periodic flow field is reached.

We furnish extensive flow data, including the vorticity and stream function fields at various instants during the temporal evolution of the flow field, time histories of the drag and lift coefficients, Strouhal number, initial and mean drag coefficients, amplitude of the drag and lift coefficient oscillations, and the phase relationships between the drag and lift oscillations associated with each cylinder. Our data confirms that, at this Reynolds number, there are no stable steady-state solutions with local symmetries. Of course, one can obtain such unphysical solutions by assuming symmetry conditions along the axes of the cylinders parallel to the main flow direction and taking half of the computational domain needed normally. In such cases, the “steady-state” flow fields obtained would be identical to the flow fields observed at the initial stages of our computations. However, we show that such flow fields do not represent the temporally periodic flow fields even in a time-averaged sense, because, in all three cases, the initial drag coefficients are different from the mean drag coefficients. Therefore, we conclude that stability studies involving periodic arrays of cylinders should be carried out, as it is done in this work, with the true implementation of the spatial periodicity.