Running the test

This test case is run with the command

make test

which runs the testXY and testXZ cases in turn. These cases can be individually run with the commands

make testXY

and

make testXZ

Central differencing for advection of momentum (suntans.dat: nonlinear = 2) is employed and nonhydrostatic pressure used for the testXZ case (suntans.dat: nonhydrostatic = 1). Small time steps of 0.001 and 0.005 are employed to resolve eddies over a total time of 40, corresponding to a rough approximation for steady state. Viscosity controls the Reynolds number and is set to yield $Re=3200$ (suntans.dat: nu = 0.0003125 and suntans.dat: nu_H = 0.0003125) for comparisons with Zang et al.  [6]. Results from the simulation can be compared with these literature data values via the Matlab scripts CompareResultXY.m and CompareResultXZ.m. A more resolved case than the examples (13,500 cells) is shown in Figure 13, where $L=1$ is the square domain's length, $U=1.0$ is lid velocity, and $u$ and $v$ are the zonal and meridional velocities, respectively.

Figure 13: Velocity profiles for refined XY cavity flow (lines) compared with literature data from Zang et al.  [6](squares)

0.8