The mild-slope equation models the propagation and transformation of water waves, as they travel through waters of varying depth and interact with lateral boundaries such as cliffs, beaches, seawalls and breakwaters. As a result, it describes the variations in wave amplitude, or equivalently wave height. From the wave amplitude, the amplitude of the flow velocity oscillations underneath the water surface can also be computed. These quantities—wave amplitude and flow-velocity amplitude—may subsequently be used to determine the wave effects on coastal and offshore structures, ships and other floating objects, sediment transport and resulting geomorphology changes of the sea bed and coastline, mean flow fields and mass transfer of dissolved and floating materials. Most often, the mild-slope equation is solved by computer using methods from numerical analysis.
A first form of the mild-slope equation was developed by Eckart in 1952, and an improved version—the mild-slope equation in its classical formulation—has been derived independently by Juri Berkhoff in 1972. Thereafter, many modified and extended forms have been proposed, to include the effects of, for instance: wave–current interaction, wave nonlinearity, steeper sea-bed slopes, bed friction and wave breaking. Also parabolic approximations to the mild-slope equation are often used, in order to reduce the computational cost.
In case of a constant depth, the mild-slope equation reduces to the Helmholtz equation for wave diffraction.
The last equation shows that wave energy is conserved in the mild-slope equation, and that the wave energy is transported in the -direction normal to the wave crests (in this case of pure wave motion without mean currents). The effective group speed is different from the group speed
The first equation states that the effective wavenumber is irrotational, a direct consequence of the fact it is the derivative of the wave phase , a scalar field. The second equation is the eikonal equation. It shows the effects of diffraction on the effective wavenumber: only for more-or-less progressive waves, with the splitting into amplitude and phase leads to consistent-varying and meaningful fields of and . Otherwise, κ2 can even become negative. When the diffraction effects are totally neglected, the effective wavenumber κ is equal to , and the geometric optics approximation for wave refraction can be used.
Details of the derivation of the above equations
When is used in the mild-slope equation, the result is, apart from a factor :
Now both the real part and the imaginary part of this equation have to be equal to zero:
The effective wavenumber vector is defined as the gradient of the wave phase:
Luke's Lagrangian formulation gives a variational formulation for non-linear surface gravity waves. For the case of a horizontally unbounded domain with a constant density, a free fluid surface at and a fixed sea bed at Luke's variational principle uses the Lagrangian
where is the velocity potential, with the flow velocity components being and in the , and directions, respectively. Luke's Lagrangian formulation can also be recast into a Hamiltonian formulation in terms of the surface elevation and velocity potential at the free surface. Taking the variations of with respect to the potential and surface elevation leads to the Laplace equation for in the fluid interior, as well as all the boundary conditions both on the free surface as at the bed at
In case of linear wave theory, the vertical integral in the Lagrangian density is split into a part from the bed to the mean surface at and a second part from to the free surface . Using a Taylor series expansion for the second integral around the mean free-surface elevation and only retaining quadratic terms in and the Lagrangian density for linear wave motion becomes
The term in the vertical integral is dropped since it has become dynamically uninteresting: it gives a zero contribution to the Euler–Lagrange equations, with the upper integration limit now fixed. The same is true for the neglected bottom term proportional to in the potential energy.
The waves propagate in the horizontal plane, while the structure of the potential is not wave-like in the vertical -direction. This suggests the use of the following assumption on the form of the potential
with normalisation at the mean free-surface elevation
Here is the velocity potential at the mean free-surface level Next, the mild-slope assumption is made, in that the vertical shape function changes slowly in the -plane, and horizontal derivatives of can be neglected in the flow velocity. So:
Now is first taken equal to and then to As a result, the evolution equations for the wave motion become:
with ∇ the horizontal gradient operator: ∇ ≡(∂/∂x ∂/∂y)T where T denotes the transpose.
The next step is to choose the shape function and to determine and
Vertical shape function from Airy wave theory
Since the objective is the description of waves over mildly sloping beds, the shape function is chosen according to Airy wave theory. This is the linear theory of waves propagating in constant depth The form of the shape function is:
with now in general not a constant, but chosen to vary with and according to the local depth and the linear dispersion relation:
Here a constant angular frequency, chosen in accordance with the characteristics of the wave field under study. Consequently, the integrals and become:
The following time-dependent equations give the evolution of the free-surface elevation and free-surface potential 
From the two evolution equations, one of the variables or can be eliminated, to obtain the time-dependent form of the mild-slope equation:
and the corresponding equation for the free-surface potential is identical, with replaced by The time-dependent mild-slope equation can be used to model waves in a narrow band of frequencies around
Consider monochromatic waves with complex amplitude and angular frequency
with and chosen equal to each other, Using this in the time-dependent form of the mild-slope equation, recovers the classical mild-slope equation for time-harmonic wave motion:
Applicability and validity of the mild-slope equation
The standard mild slope equation, without extra terms for bed slope and bed curvature, provides accurate results for the wave field over bed slopes ranging from 0 to about 1/3. However, some subtle aspects, like the amplitude of reflected waves, can be completely wrong, even for slopes going to zero. This mathematical curiosity has little practical importance in general since this reflection becomes vanishingly small for small bottom slopes.
Dingemans, M. W. (1997), Water wave propagation over uneven bottoms, Advanced Series on Ocean Engineering 13, World Scientific, Singapore, ISBN981-02-0427-2, OCLC36126836, 2 Parts, 967 pages.
Liu, P. L.-F. (1990), "Wave transformation", in B. Le Méhauté and D. M. Hanes, Ocean Engineering Science, The Sea 9A, Wiley Interscience, pp. 27–63, ISBN0-471-52856-0
Mei, Chiang C. (1994), The applied dynamics of ocean surface waves, Advanced Series on Ocean Engineering 1, World Scientific, ISBN9971-5-0789-7, 740 pages.
Porter, D.; Chamberlain, P. G. (1997), "Linear wave scattering by two-dimensional topography", in J. N. Hunt, Gravity waves in water of finite depth, Advances in Fluid Mechanics 10, Computational Mechanics Publications, pp. 13–53, ISBN1-85312-351-X