ERDC/CHL CHETN-IV-32
June 2001
Example Applied to Grays Harbor, Washington: The Grays Harbor sample problem
represents an intermediate-sized domain situated in the Pacific Ocean off of the State of
Washington and the province of British Columbia. The domain and finite element grid are
shown in Figure 13 and bathymetry/topography are shown in Figure 15. Grays Harbor is highly
resolved in the domain and Willapa Bay is represented with medium resolution. Lower-
resolution representations of the Columbia Estuary, the Strait of Juan de Fuca, the Strait of
Georgia and Queen Charlotte Sound are included.
The input files are grays_s03_g04_r11.grd (unit 14) and grays_s03_g04_r11.inp (unit 15) which
submerged jetty to the south of the inlet at Grays Harbor is represented as an internal-barrier
boundary (Figures 14 and 16). This barrier has been set up to allow for both overtopping and
leakage through the barrier itself. The crown of this internal barrier ranges between 1.5 m below
the geoid along the western portion of the jetty to 3.0 m above the geoid along its eastern portion.
In addition to barrier overtopping, this barrier incorporates the leaky internal-barrier boundary to
account for the possibility of holes in barriers that can lead to cross-barrier through flow. Two
nodes have been set up to allow leaky through flow in the region where the jetty rises above the
geoid. This boundary has been designated as IBTYPE = 25 so that the computed cross-barrier
through and over flow will be implemented as a natural boundary condition. This flow boundary
has been implemented as a leaky boundary by specifying two node pairs as having a pipe crown
height located at PIPEHT = 1.0 m, a pipe friction factor PIPECOEF = 0.1, and a pipe diameter
PIPEDIAM = 0.25 m. At all boundary node pairs where no cross-barrier pipe flow is desired, a
pipe crown height equal to 100 m has been specified.
Figure 13. Domain definition and finite element grid
for the Grays Harbor case
16