A01 shoaling and refraction


The purpose of this test is to verify the correct implementation of shoaling and refraction.

A domain with flat Inclined bottom is considered. The depth decreases from 10 m at y=0 to depth -1 m at y=4000 m.

At the Southern boundary (y=0) an incident wave field is assumed which propagates in a direction 30° from S. In order to come close to the analytical solution this incident wave field has a narrow frequency spectrum with a narrow directional distribution. The peak period is 10 s.

The analytical results are computed as follows; the wave direction θ (measured from the normal to the coast) results from Snel's law:

sin θ / cphase = sin θ0 / cphase,0

and the wave height from: cg cosθ H2 = cg,0 cosθ0 H02

In plotfile A01.plt computed values for Hs and average direction are compared with the analytical results.

9 cases are considered:

A01a
1d computation
A01b
regular 2d grid in standard orientation
A01c
regular grid rotated over 10 degr.
A01d
curvilinear 'wavy' grid, the grid can be seen in plot A01d.hpg
A01e
curvilinear 'jagged' grid, the grid can be seen in plot A01e.hpg
A01f
curvilinear 'jagged' grid, the grid can be seen in plot A01f.hpg
A01g
spherical coordinates
A01h
repeating grid
A01j
unstructured grid

In case A01h a table A01h.tbl is written which shows values of a number of variables along a straight line parallel to the coast; all variables should be constant along this line.

The central part of the curvilinear grid in case A01d
The central part of the jagged grid in case A01e
The central part of the jagged grid in case A01f

© 2012: Nico Booij