next up previous contents
Next: Internal waves Up: Cavity flow Previous: Boundary conditions   Contents

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\includegraphics{figures/cavity}


next up previous contents
Next: Internal waves Up: Cavity flow Previous: Boundary conditions   Contents
2014-08-06