J.L. Hench et al. / Continental Shelf Research 22 (2002) 26152631
2628
for inlets based on a dynamical analysis. The four
qV
qV
qV
V
U
fU
idealized inlets analyzed here provide baseline
qy |ffl{zffl}
qt
qx
|{z} |fflfflfflfflfflfflfflfflfflfflfflfflffl{zfflfflfflfflfflfflfflfflfflfflfflfflffl} M
cases for the dynamics one might expect at natural
y3
My1
My2
inlets.
pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi!
Cf U 2 V 2
qZ
g
V 0:
A:2
H
qy
|fflffl{zfflffl} |fflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflffl{zfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflffl}
Acknowledgements
My4
My5
We thank two anonymous reviewers for con-
For simplicity, horizontal diffusion terms are
structive comments that significantly improved the
omitted from (A.1) and (A.2) as they did not
manuscript. Funding for this study was provided
contribute significantly to the momentum balances
by the Coastal Inlets Research Program (CIRP) of
of our inlet simulations. At each time-step, we
define a two-dimensional, orthogonal, curvilinear
Corps of Engineers, and the Naval Research
coordinate system such that at each grid point, one
Laboratory, Office of Naval Research (award
coordinate direction points in the along-stream
N00014-97-C-6010).
direction s and the other points in the across-
stream direction n in the right hand sense (see
Fig. 8). The orientation of the s2n coordinate
Appendix A. Rotation of momentum equations into
system relative to the original fixed x2y coordi-
a streamwise-normal coordinate system
nate system, as expressed by the streamline angle
ax; y; t; varies such that at all points the along-
The use of a streamwise-normal coordinate
stream velocity Us is equivalent to the speed, and
system allows a more intuitive physical interpreta-
the across-stream velocity Un 0: From Fig. 8, it
tion for strongly curving flow fields. Gill (1982)
should be apparent that momentum terms in s2n
elegantly derives the steady frictionless form of the
two-dimensional momentum equations directly in
the s2n coordinate system. Here we derive
corresponding transient frictional equations by
transforming the familiar x2y equations. This
transformation is useful in that standard x2y
momentum terms can be used to directly compute
term values in the s2n coordinate system. The
procedure is as follows. U x; y; t; V x; y; t and
Zx; y; t are computed on a fixed x2y grid for all
time steps. The elevation and velocity fields are
used to reconstruct each term in the x2y
momentum equations at each node in the x2y
grid. We designate the individual x2y momentum
terms Mxi and Myi as
qU
qU
qU
V
U
fV
qy |ffl{zffl}
t
qx
qz}
Fig. 8. Definition sketch for the streamline coordinate system.
|{ |fflfflfflfflfflfflfflfflfflfflfflfflffl{zfflfflfflfflfflfflfflfflfflfflfflfflffl} M
At each computational node a s2n coordinate system is
x3
Mx1
Mx2
established and aligned with the local velocity vector U : The
pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi!
local axes rotation angle is ax; y; t: Mi represents the ith force
Cf U 2 V 2
qZ
g
U 0;
A:1
H
qfflx
components Mxi and Myi (e.g. g qZ=qx and g qZ=qy). Mi is
|fflffl{zffl }
|fflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflffl{zfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflffl}
rotated on to the s2n axes to determine local streamwise and
Mx4
Mx5
normal components Msi and Mni (e.g. g qZ=qs and g qZ=qn).