Previous Up Next

9.2  PASS: Coastally-trapped waves

Author
Stéphane Popinet
Command
sh waves.sh waves.gfs
Version
1.0.0
Required files
waves.gfs (view) (download)
waves.sh solution.gfv correlation.ref
Running time
7 minutes 14 seconds

A simulation is initialised with an exact solution found by Lamb [11] for inertia-gravity waves in a circular, flat-bottomed basin on an f-plane.

The waves should rotate unchanged around the basin. This is a good test of the shallow-water solver with Coriolis terms and non-trivial coastlines.

Figure 76 illustrates the solution obtained after three wave periods.

Following [5], Table 13 presents the phase error and the maximum correlation between the computed solution and the exact solution. Note that these results are not directly comparable to [5], because the timestep used is much larger. The time-explicit codes tested in [5] have stronger stability requirements than Gerris.


Figure 76: Surface-height after a rotation of 3 wave periods for a resolution of 9.375 km.


Table 13: Maximum correlation and phase error as functions of resolution.
LevelMaximum CAngle of max C
50.86452.7
60.97170.6
70.99230.2

9.2.1  PASS: Coastally-trapped waves with adaptive refinement

Author
Stéphane Popinet
Command
sh ../waves.sh adaptive.gfs
Version
1.0.0
Required files
adaptive.gfs (view) (download)
solution.gfv correlation.ref
Running time
4 minutes 25 seconds

Same test case as before but with adaptive refinement (using the vorticity criterion).

Figure 77 illustrates the solution obtained after three wave periods.


Figure 77: Surface-height after a rotation of 3 wave periods. The maximum resolution is 9.375 km (dark red), the minimum resolution is 150 km (dark blue).


Table 14: Maximum correlation and phase error as functions of resolution.
LevelMaximum CAngle of max C
50.86482.6
60.97070.6
70.9910.4


Previous Up Next