Dispersion relation for water waves (Scripts) Publisher's description
from Frederic Moisy
This set of functions simply provides an easy way to work with the dispersion relation of surface waves
This set of functions simply provides an easy way to work with the dispersion relation of surface waves, given by
omega(k) = sqrt ( tanh(k*h0) * (g*k + gamma*k^3/rho))
where omega is the pulsation (in rad/s), k the wavenumber (in 1/m), h0 the depth, g the gravity, gamma the surface tension and rho the density.
The function kfromw allows one to invert the dispersion relation, i.e. to give the value of omega for a given value of k.
(For infinite depth, kfromw simply inverts the cubic polynomial. For finite depth, a zero-finding method is used, starting from the infinite depth solution).
By default, the physical parameters (liquid densities, surface tension, etc.) are set for an air-water interface under usual conditions, with a water layer of infinite depth ("deep water waves"). Use the function wave_parameter to change those properties.
See the published file "demo" to learn more about this package.
System Requirements:MATLAB 7 (R14)
Program Release Status: New Release
Program Install Support: Install and Uninstall