This paper describes the mathematical formulation, the numerical solution,
and the validation of a linear refraction-diffraction model for steep bathy
metry. The model involves two coupled governing equations derived from, res
pectively, the exact seabed boundary condition and the Laplace equation. It
reduces to the extended and the original mild-slope model, when the seabed
slope is small. Although the present approach is based on depth-integratio
n of flow characteristics, it correctly accounts for the vertical component
of the seabed fluid velocity. The formulation is based on the weighted-res
idual method, and the hybrid element solution is derived from a Galerkin ap
proach. The capability of the present model to simulate flow velocity and w
ave amplitude over three-dimensional bedforms is examined in a parametric s
tudy. The computed results are compared with the original and extended mild
-slope solutions and verified with those of a three-dimensional wave model.
The present depth-integrated model has the same data requirements as other
two-dimensional models, but provides accurate three-dimensional results wi
th only a fraction of the CPU time that would be required by a three-dimens
ional model.