Code covered by the BSD License  

Highlights from
Dispersion relation for water waves

Be the first to rate this file! 31 Downloads (last 30 days) File Size: 40.5 KB File ID: #27391

Dispersion relation for water waves

by

 

28 Apr 2010 (Updated )

Dispersion relation, and its inverse, for surface waves (eg, finding wavenumber from frequency).

| Watch this File

File Information
Description

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.

MATLAB release MATLAB 7 (R14)
Tags for This File   Please login to tag files.
Please login to add a comment or rating.
Updates
30 Apr 2010

typing errors

Contact us