Comments:
The computational grid, named phriv for the Great Bay Estuary system is created using linear equilateral triangles fitted to the shoreline boundary. There are 12113 nodes and 21235 elements. The bathymetric depth contour reveals a center channel running south to north.

M2 tidal forcing is specified as a Dirichlet elevation boundary condition across the southern and northern inlet boundaries. The amplitude is 1.35m at the southern boundary. At the northern boundary the values obtained from gbriv mesh run are used as a boundary condition. The simulation was started with fluid at rest and was terminated after five M2 tidal period (62.10hrs). The simulation setup is given below. All the simulation parameters are summarized in Table below.

Residual transport:

Comments:

Comments:
The above transient solutions are numerically well-behaved and appear to describe the realistic dynamics of the tidal flooding and dewatering process in this complex domain.