4.8. Results

The velocity field of the fluid flow through the midplane on the z-direction of the cavity at Re=100 is shown in Figure 4.6: Velocity field through the cavity for Re=100., while Figure 4.7: Streamlines through the cavity for Re = 100: Ghia at alshows the streamline contours for the cavity flow on the same midplane. Figure 4.7: Streamlines through the cavity for Re = 100: Ghia at al(a) shows the streamlines obtained by [7], while Figure 4.7: Streamlines through the cavity for Re = 100: Ghia at al(b) shows the streamlines from the Rocky simulation.

A qualitative comparison can be done and one can see that the solutions present similar behavior. It is possible to observe that a vortex is generated by the shear force of the moving wall (lid) and the core of the vortex places at the upper half of the cavity.

Figure 4.6: Velocity field through the cavity for Re=100.

Velocity field through the cavity for Re=100.

Figure 4.7: Streamlines through the cavity for Re = 100: Ghia at al

Streamlines through the cavity for Re = 100: Ghia at al

Figure 4.8: Streamlines through the cavity for Re = 100: Rocky simulation.

Streamlines through the cavity for Re = 100: Rocky simulation.

Relevant verification parameters used in the literature for the lid-driven cavity problem are the horizontal fluid velocity (u) along the cavity vertical centerline, as well as the vertical fluid velocity (v) along the horizontal cavity centerline. Comparisons of the u and v fluid velocities along the lines through the cavity center for Re=100 are shown in Figure 4.9: Comparison of u-velocity along the vertical line through the cavity center for Re = 100. and Figure 4.10: Comparison of v-velocity along the horizontal line through the cavity center for Re=100., respectively. The fluid velocities calculated by Rocky presents strongly correlated values to the reference ones.

Figure 4.9: Comparison of u-velocity along the vertical line through the cavity center for Re = 100.

Comparison of u-velocity along the vertical line through the cavity center for Re = 100.

Figure 4.10: Comparison of v-velocity along the horizontal line through the cavity center for Re=100.

Comparison of v-velocity along the horizontal line through the cavity center for Re=100.

The comparisons of the absolute and relative errors between the velocities results for Re=100 are shown in Figure 4.11: Absolute and relative errors for u-velocity along the vertical line through the cavity center for Re = 100. and Figure 4.12: Absolute and relative errors v-velocity along the horizontal line through the cavity center for Re = 100.. The relative errors are mostly below 30% and higher differences can be observed at the points close to the stationary walls and other points where the velocities values are close to zero, in which high relative differences are expected.

Figure 4.11: Absolute and relative errors for u-velocity along the vertical line through the cavity center for Re = 100.

Absolute and relative errors for u-velocity along the vertical line through the cavity center for Re = 100.

Figure 4.12: Absolute and relative errors v-velocity along the horizontal line through the cavity center for Re = 100.

Absolute and relative errors v-velocity along the horizontal line through the cavity center for Re = 100.