Estuarine Flooding and Dewatering Simulation



Simulation Results of Piscataqua River, New Hampshire





Computational grid and simulation setup

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

Simulation parameters for the Piscataqua River.


Mass balance

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.


Tidally averaged residual analysis

Residual transport:

pr6-resq.gif (10.9kb) Upper and Lower Piscataqua River

Residual velocity:
pr6-resv.gif (11.3kb) Upper and Lower Piscataqua Rive

Residual bottom stress:
pr6-ress.gif (10.7kb) Upper and Lower Piscataqua River

Residual sediment transport:
gbay6-rest.gif (9.7kb) Upper and Lower Piscataqua River

Comments:


Transient solution

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.

[ Introduction ] [ Great Bay ]
[ Bellamy River ] [ Oyster River ] [ Piscataqua River ] [ Portsmouth Harbour ]

[Oyster,Bellamy and Piscataqua Rivers] [Great Bay with Oyster and Bellamy Rivers]
[Portsmouth Harbour with Piscataqua River] [Whole Great Bay Estuary System]


[back]



[email protected]

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

1
Hosted by www.Geocities.ws