Circular Basin

UNDER CONSTRUCTION

Analytical Solution

Dupont (2001) presented an analytical solution for a closed circular domain on an f-plane, with radius $R$, a linear bottom friction, and a spatially variable wind stress equal to $\tau _{Wx}=Wy/R$, $\tau _{Wy}=0$ where $W$ is the gradient of the wind forcing and $y$ is the vertical coordinate. The water surface elevation solution is given by

$\eta ={\begin{cases}{\frac {W{{r}^{2}}\sin 2\theta }{4gHR}},&{\mbox{if }}f_{c}=0\\{\frac {W{{f}_{c}}}{RgH\kappa }}\left[{\frac {{R}^{2}}{8}}+{\frac {{r}^{2}}{4}}\left({\frac {\kappa }{{f}_{c}}}\sin 2\theta -1\right)\right],&{\mbox{if }}f_{c}\neq 0\\\end{cases}}$

(1)

The current velocities are independent of the Coriolis parameter and are given by

$u={\frac {Wy}{2R\kappa }}$

(2)

$v=-{\frac {Wx}{2R\kappa }}$

(3)

Setup

Figure 1. Computational grid.

The model is run to steady state from zero current and water level initial conditions with $W=10^{-}4{\text{m}}^{2}{\text{s}}^{-2}$, $\kappa =10^{-3}{\text{s}}^{-1}$ , and $f_{c}=0$. Table 1 shows the general settings used for CMS-Flow. Figure 1 shows the computational grid with 5 levels of refinement from 2000 m to 125 m.

Table 1. General Settings for Wind-driven flow in a circular basin

Parameter	Value
Time step	3600 s
Simulation Duration	72 hrs
Ramp Period	24 hrs
Initial Water Depth	10 m
Mixing Terms	Off
Wall Friction	Off
Linear Bottom Friction Coefficient	0.001

Results

Table 2. Goodness of fit statistics for the current velocity and water level

Variable	RRMSE, %	RMAE, %	R^2	Bias
U-Velocity	3.88	0.64	0.997	-4.06e-5
V-Velocity	3.87	0.64	0.997	4.06e-5
Water Level	0.16	0.13	1.000	-3.56e-6

References

• Dupont, F., 2001. Comparison of numerical methods for modelling ocean circulation in basins with irregular coasts. Ph.D. thesis, McGill University, Montreal.

---

Test Cases

Documentation Portal