Extended mild-slope equations for the propagation of small-amplitude water waves over variable bathymetry regions, recently proposed by Massel (1993) and Porter & Staziker (1995), are shown to exhibit an inconsistency concerning the sloping-bottom boundary condition, which renders them non-conservative with respect to wave energy. In the present work, a consistent coupled-mode theory is derived from a variational formulation of the complete linear problem, by representing the vertical distribution of the wave potential as a uniformly convergent series of local vertical modes at each horizontal position. This series consists of the vertical eigenfunctions associated with the propagating and all evanescent modes and, when the slope of the bottom is different from zero, an additional mode, carrying information about the bottom slope. The coupled-mode system obtained in this way contains an additional equation, as well as additional interaction terms in all other equations, and reduces to the previous extended mild-slope equations when the additional mode is neglected. Extensive numerical results demonstrate that the present model leads to the exact satisfaction of the bottom boundary condition and, thus, it is energy conservative. Moreover, it is numerically shown that the rate of decay of the modal-amplitude functions is improved from O(n⁻²), where n is the mode number, to O(n⁻⁴), when the additional sloping-bottom mode is included in the representation. This fact substantially accelerates the convergence of the modal series and ensures the uniform convergence of the velocity field up to and including the boundaries.

The consistent coupled-mode theory (Athanassoulis & Belibassakis, J. Fluid Mech. vol. 389, 1999, p. 275) is extended and applied to the hydroelastic analysis of large floating bodies of shallow draught or ice sheets of small and uniform thickness, lying over variable bathymetry regions. A parallel-contour bathymetry is assumed, characterized by a continuous depth function of the form $h( {x,y}) {=} h( x )$, attaining constant, but possibly different, values in the semi-infinite regions $x {<} a$ and $x {>} b$. We consider the scattering problem of harmonic, obliquely incident, surface waves, under the combined effects of variable bathymetry and a floating elastic plate, extending from $ x {=} a$ to $x {=} b$ and $ {-} \infty {<} y{<}\infty $. Under the assumption of small-amplitude incident waves and small plate deflections, the hydroelastic problem is formulated within the context of linearized water-wave and thin-elastic-plate theory. The problem is reformulated as a transition problem in a bounded domain, for which an equivalent, Luke-type (unconstrained), variational principle is given. In order to consistently treat the wave field beneath the elastic floating plate, down to the sloping bottom boundary, a complete, local, hydroelastic-mode series expansion of the wave field is used, enhanced by an appropriate sloping-bottom mode. The latter enables the consistent satisfaction of the Neumann bottom-boundary condition on a general topography. By introducing this expansion into the variational principle, an equivalent coupled-mode system of horizontal equations in the plate region ($a {\leq} x {\leq} b)$ is derived. Boundary conditions are also provided by the variational principle, ensuring the complete matching of the wave field at the vertical interfaces ($x{=}a$ and $x{=}b)$, and the requirements that the edges of the plate are free of moment and shear force. Numerical results concerning floating structures lying over flat, shoaling and corrugated seabeds are presented and compared, and the effects of wave direction, bottom slope and bottom corrugations on the hydroelastic response are presented and discussed. The present method can be easily extended to the fully three-dimensional hydroelastic problem, including bodies or structures characterized by variable thickness (draught), flexural rigidity and mass distributions.