Dispersion Relation for a Floating Elastic Plate

Separation of Variables

The dispersion equation arises when separating variables subject to the boundary conditions for a Floating Elastic Plate of infinite extent. The same equation arises when separating variables in two or three dimensions and we present here the two-dimensional version. The equations are described in detail in the Floating Elastic Plate page and we begin with the equations The equations of motion for the Frequency Domain Problem with radial frequency [math]\displaystyle{ \,\omega }[/math] in terms of the potential alone is

[math]\displaystyle{ D\frac{\partial^4}{\partial x^4} \frac{\partial \phi}{\partial z} + \left(\rho g - \omega^2 \rho_i h \right) \frac{\partial \phi}{\partial z} = - \rho \omega^2 \phi, \, z=0 }[/math]

plus the equations within the fluid

[math]\displaystyle{ \nabla^2\phi =0 }[/math]

[math]\displaystyle{ \frac{\partial \phi}{\partial z} = 0, \, z=-H }[/math]

where [math]\displaystyle{ \,g }[/math] is the acceleration due to gravity, [math]\displaystyle{ \,\rho_i }[/math] and [math]\displaystyle{ \,\rho }[/math] are the densities of the plate and the water respectively, [math]\displaystyle{ \,h }[/math] and [math]\displaystyle{ \,D }[/math] are the thickness and flexural rigidity of the plate.

We then look for a separation of variables solution to Laplace's Equation and obtain the following expression for the velocity potential

[math]\displaystyle{ \phi(x,z) = e^{ikx} \cosh k(z+H) \, }[/math]

If we then apply the condition at [math]\displaystyle{ z=0 }[/math] we see that the constant [math]\displaystyle{ \,k }[/math] (which corresponds to the wavenumber) is given by

[math]\displaystyle{ -D k^5 \sinh(kH) - k \left(\rho g - \omega^2 \rho_i h \right) \sinh(kH) = -\rho \omega^2 \cosh(kH) \,\,\,(1) }[/math]

This is the dispersion equation for a floating elastic plate.

Solution of the dispersion equation

The dispersion equation was first solved by Fox and Squire 1994 and they determined that the solution consists of one real, two complex, and infinite number of imaginary roots plus their negatives. Interestingly the eigenfunctions form an over complete set for [math]\displaystyle{ L_2[-H,0]\, }[/math]. Also, there are some circumstances (non-physical) in which the complex roots become purely imaginary. The solution of this dispersion equation is far from trivial and the optimal solution method has been developed by Tim Williams and is described below.

Non-dimensional form

The dispersion equation (1) is often given in non-dimensional form. The form used by Tayler 1986 is to scale length with respect to [math]\displaystyle{ L\, }[/math] and time with respect to [math]\displaystyle{ \sqrt{L/g}\, }[/math]. The non-dimensional equations then become

[math]\displaystyle{ -\beta k^5 \sinh(kH) - k \left(1 - \alpha \gamma \right) \sinh(kH) = -\alpha \cosh(kH) \, }[/math]

where [math]\displaystyle{ H }[/math] is the nondimensional water depth, [math]\displaystyle{ \alpha = L\omega^2/g,\, }[/math] [math]\displaystyle{ \beta = D/(\rho g L^4)\, }[/math] and [math]\displaystyle{ \gamma = \rho_i h/(\rho L)\, }[/math]. This equation has four parameters [math]\displaystyle{ \alpha,\beta,\gamma,H\, }[/math]. We can trivially remove one of these by setting the [math]\displaystyle{ L=H\, }[/math]. Alternatively we can set the length parameter to the Characteristic Length [math]\displaystyle{ L = (D/(\rho g))^{1/4}\, }[/math] (so that in the notation above [math]\displaystyle{ \beta = 1\, }[/math] and we also have three parameters).

Tim Williams developed the following non-dimensionalisation

[math]\displaystyle{ \cosh(k H)-(k^4+\varpi)k\sinh(k H)=0, }[/math]

where [math]\displaystyle{ \varpi=(1-k_\infty\sigma)/(k_\infty L),\quad k_\infty=\omega^2/g,\quad\sigma=\rho_ih/\rho,\quad L^5=D/(\rho\omega^2). }[/math]

[math]\displaystyle{ k_\infty }[/math] is the wave number for a wave traveling in open water of infinite depth, and [math]\displaystyle{ \sigma\, }[/math] is the amount that the plate would be submerged relative to a region of open water. Note that if [math]\displaystyle{ \sigma\, }[/math] is greater than [math]\displaystyle{ 1/k_\infty, }[/math] or [math]\displaystyle{ 2\pi\, }[/math] times the infinite depth open water wavelength, then the parameter [math]\displaystyle{ \varpi }[/math] becomes negative.

The dispersion relation relates the wavenumber [math]\displaystyle{ k/L\, }[/math] and thus wave speed [math]\displaystyle{ L\omega/k\, }[/math] to the above parameters.