Whitham equation

From HandWiki
Revision as of 16:27, 6 February 2024 by NBrush (talk | contribs) (change)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Short description: Non-local model for non-linear dispersive waves

In mathematical physics, the Whitham equation is a non-local model for non-linear dispersive waves. [1][2][3]

The equation is notated as follows:

[math]\displaystyle{ \frac{\partial \eta}{\partial t} + \alpha \eta \frac{\partial \eta}{\partial x} + \int_{-\infty}^{+\infty} K(x-\xi)\, \frac{\partial \eta(\xi,t)}{\partial \xi}\, \text{d}\xi = 0. }[/math]

This integro-differential equation for the oscillatory variable η(x,t) is named after Gerald Whitham who introduced it as a model to study breaking of non-linear dispersive water waves in 1967.[4] Wave breaking – bounded solutions with unbounded derivatives – for the Whitham equation has recently been proven.[5]

For a certain choice of the kernel K(x − ξ) it becomes the Fornberg–Whitham equation.

Water waves

Using the Fourier transform (and its inverse), with respect to the space coordinate x and in terms of the wavenumber k:

  • For surface gravity waves, the phase speed c(k) as a function of wavenumber k is taken as:[4]
[math]\displaystyle{ c_\text{ww}(k) = \sqrt{ \frac{g}{k}\, \tanh(kh)}, }[/math]   while   [math]\displaystyle{ \alpha_\text{ww} = \frac{3}{2} \sqrt{\frac{g}{h}}, }[/math]
with g the gravitational acceleration and h the mean water depth. The associated kernel Kww(s) is, using the inverse Fourier transform:[4]
[math]\displaystyle{ K_\text{ww}(s) = \frac{1}{2\pi} \int_{-\infty}^{+\infty} c_\text{ww}(k)\, \text{e}^{iks}\, \text{d}k = \frac{1}{2\pi} \int_{-\infty}^{+\infty} c_\text{ww}(k)\, \cos(ks)\, \text{d}k, }[/math]
since cww is an even function of the wavenumber k.
[math]\displaystyle{ c_\text{kdv}(k) = \sqrt{gh} \left( 1 - \frac{1}{6} k^2 h^2 \right), }[/math]   [math]\displaystyle{ K_\text{kdv}(s) = \sqrt{gh} \left( \delta(s) + \frac{1}{6} h^2\, \delta^{\prime\prime}(s) \right), }[/math]   [math]\displaystyle{ \alpha_\text{kdv} = \frac{3}{2} \sqrt{\frac{g}{h}}, }[/math]
with δ(s) the Dirac delta function.
[math]\displaystyle{ K_\text{fw}(s) = \frac12 \nu \text{e}^{-\nu |s|} }[/math]   and   [math]\displaystyle{ c_\text{fw} = \frac{\nu^2}{\nu^2+k^2}, }[/math]   with   [math]\displaystyle{ \alpha_\text{fw}=\frac32. }[/math]
The resulting integro-differential equation can be reduced to the partial differential equation known as the Fornberg–Whitham equation:[6]
[math]\displaystyle{ \left( \frac{\partial^2}{\partial x^2} - \nu^2 \right) \left( \frac{\partial \eta}{\partial t} + \frac32\, \eta\, \frac{\partial \eta}{\partial x} \right) + \frac{\partial \eta}{\partial x} = 0. }[/math]
This equation is shown to allow for peakon solutions – as a model for waves of limiting height – as well as the occurrence of wave breaking (shock waves, absent in e.g. solutions of the Korteweg–de Vries equation).[6][3]

Notes and references

Notes

  1. (Debnath 2005)
  2. (Naumkin Shishmarev)
  3. 3.0 3.1 (Whitham 1974)
  4. 4.0 4.1 4.2 4.3 (Whitham 1967)
  5. (Hur 2017)
  6. 6.0 6.1 6.2 (Fornberg Whitham)

References