A02. Waves propagating towards a beach

The bathymetry in this set of tests is identical to the previous one (A01); this time source terms are activated. The purpose of this test is to verify the correct implementation of the complete model for various grids (1-d, regular, curvilinear, repeated). Plot output includes output quantity ′Force′.

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 an incident wave field is assumed which propagates in a direction 30° from S. The peak period is 10 s.

In this series of tests there are no analytical tests to compare with, so the 1-d results are used as a basis for comparison in the graphs (see plotfile A02.ps).

11 cases are considered:

1d computation
regular 2d grid in standard orientation
regular grid rotated over 10 degr.
curvilinear 'wavy' grid, the grid can be seen in plot A02d.hpg
curvilinear 'jagged' grid, the grid can be seen in plot A02e.hpg
curvilinear 'jagged' grid, the grid can be seen in plot A02f.hpg
spherical coordinates
repeating grid
unstructured grid
unstructured grid, uses hotfile generated by A02j as Initial state
nonstationary 1d computation
nonstationary computation on a regular 2d grid in standard orientation
The results are plotted in file A02.ps; values of Hs along the transect normal to the coast are plotted on the first page; values of the average direction on the second. The first panel on each page compares the analytic solution (with zero source term) with the 1D computation. The other panels Compare the 1D solution with various 2D grid solutions, regular as well as curvilinear grids.

In case A02h a table A02h.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.
In the nonstationary computations a time-dependent table is made for the point located at 3000 m from the incident wave boundary, in order to show how the nonstationary computation converges to the stationary solution.

The central part of the curvilinear grid in case A02d
The central part of the jagged grid in case A02e
The central part of the jagged grid in case A02f

© 2012: Nico Booij