The figures show velocity magnitude and pressure distributions of the driven cavity flow with Reynolds number Re=1000. The flow was assumed incompressible. The ElmerSolver uses a stabilized finite element formulation to solve the incompressible Navier-Stokes equations. In this simulation bilinear quadrilateral finite elements were used.

If you want to reproduce the results, download the file Cavity.tar.gz and follow the instructions in the file README. You should have the Elmer package installed in your computer, however.

The velocity results with four different Reynolds number (100, 400, 1000 and 3200) are compared below with the computed results by Ghia et al. in Journal of Computational Physics (48), pp 347-411, 1982. The mesh size in Elmer computations was 128x128 of 4-nodes elements. The upper figure below presents the values of u-velocity along the vertical line through the geometric center of the cavity. The following figure present v-velocity values on a corresponding horizontal line. The computed values (red lines) match very well with the reference values (the circles).