Estuarine Flooding and Dewatering Simulation |
| pr6-mesh.gif (8.3kb) | Finite element grid (3417 nodes and 6043 elements) |
| pr6-bat2d.gif (18.3kb) | Color map of bathymetric depth |
Comments:
The computational grid, named
pr6 for the
Piscataqua River is created using linear equilateral triangles
fitted to the shoreline boundary. There are 3417 nodes and 6043 elements.
The bathymetric depth contour
reveals a center channel running north to
south.
M2 tidal forcing is specified as a Dirichlet elevation boundary conditions
across the southern inlet boundary at the bottom and across the northern inlet boundary at the top . The
amplitude is 1.00m at the bottom boundary and 0.93m at the top boundary. The
simulation was started with fluid at rest and was terminated after
six M2 tidal period (74.52hrs). The simulation setup is
given below. All the simulation parameters are summarized in Table below.
| Description | Parameters |
| Bathymetry range | 0 - 18.93m |
| Porous layer thickness | 0.25m |
| Hydraulic conductivity | 0.0003162 |
| Drag coefficient | 0.0025 |
| Time increment | 111.78sec |
| Time steps per tidal period | 400 |
| Tidal periodicity | 12.42hrs |
| Duration of simulation | 74.52hrs (6 cycles) |
| Numerical implicity | 1 (fully implicit) |
| Number of nonlinear iterations | 10 |
| pr6-mb2.gif (11.9kb) | Time series of the total fluid volume and the cumulative transport across the four transects |
| pr6-ts4.gif (13.9kb) | Time series of normal velocity, normal transport, normal bottom stress, and normal sediment transport across the four transects |
Comments:
Time history of the total fluid volume and the cumulative transport across
four selected transects as shown in
pr6
show the system establishes a dynamic equilibrium rapidly, and mass
conservation is thereafter maintained throughout the simulation after
the initial tidal period. The maximum normal velocity across the transects
increases then decreases as one moves northward with maximum level across
the narrowest transect, while the normal transport across the same transects
decreases to maintain mass conservation. It is noticed that asymmetry exists
between flood and ebb.
Residual transport:
| pr6-resq.gif (10.9kb) | Upper and Lower Piscataqua River |
| pr6-resv.gif (11.3kb) | Upper and Lower Piscataqua Rive |
| pr6-ress.gif (10.7kb) | Upper and Lower Piscataqua River |
| gbay6-rest.gif (9.7kb) | Upper and Lower Piscataqua River |
Comments:
| pr6-ebb.gif (17.9kb) | Maximum ebb at dynamic equilibrium. |
| pr6-fld.gif (19.4kb) | Maximum flood at dynamic equilibrium. |
| pr6_hv.fli (0.42Mb) | FLC animation of tidal time transport (H*v) |
| pr6-v.fli (0.34Mb) | FLC animation of tidal time velocity (v) |
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.


Last modified: September 21, 1998 (Safak Nur ERTURK)