36 Integrals with Coalescing SaddlesApplications

§36.13 Kelvin’s Ship-Wave Pattern

A ship moving with constant speed V on deep water generates a surface gravity wave. In a reference frame where the ship is at rest we use polar coordinates r and ϕ with ϕ=0 in the direction of the velocity of the water relative to the ship. Then with g denoting the acceleration due to gravity, the wave height is approximately given by

36.13.1 z⁡(ϕ,ρ)=∫−π/2π/2cos⁡(ρ⁢cos⁡(θ+ϕ)cos2⁡θ)⁢dθ,

where

36.13.2 ρ=g⁢r/V2.

The integral is of the form of the real part of (36.12.1) with y=ϕ, u=θ, g=1, k=ρ, and

36.13.3 f⁡(θ,ϕ)=−cos⁡(θ+ϕ)cos2⁡θ.

When ρ>1, that is, everywhere except close to the ship, the integrand oscillates rapidly. There are two stationary points, given by

36.13.4 θ+⁢(ϕ) =12⁢(arcsin⁡(3⁢sin⁡ϕ)−ϕ),
θ−⁢(ϕ) =12⁢(π−ϕ−arcsin⁡(3⁢sin⁡ϕ)).

These coalesce when

36.13.5 |ϕ|=ϕc=arcsin⁡(13)=19∘⁢.47122.

This is the angle of the familiar V-shaped wake. The wake is a caustic of the “rays” defined by the dispersion relation (“Hamiltonian”) giving the frequency ω as a function of wavevector 𝐤:

36.13.6 ω⁡(𝐤)=g⁢k+𝐕⋅𝐤.

Here k=|𝐤|, and 𝐕 is the ship velocity (so that V=|𝐕|).

The disturbance z⁡(ρ,ϕ) can be approximated by the method of uniform asymptotic approximation for the case of two coalescing stationary points (36.12.11), using the fact that θ±⁢(ϕ) are real for |ϕ|<ϕc and complex for |ϕ|>ϕc. (See also §2.4(v).) Then with the definitions (36.12.12), and the real functions

36.13.7 u⁡(ϕ) =Δ1/2⁢(ϕ)2⁢(1f+′′⁢(ϕ)+1−f−′′⁢(ϕ)),
v⁡(ϕ) =12⁢Δ1/2⁢(ϕ)⁢(1f+′′⁢(ϕ)−1−f−′′⁢(ϕ)),

the disturbance is

36.13.8 z⁡(ρ,ϕ)=2⁢π⁢(ρ−1/3⁢u⁡(ϕ)⁢cos⁡(ρ⁢f~⁢(ϕ))⁢Ai⁡(−ρ2/3⁢Δ⁢(ϕ))⁢(1+O⁡(1/ρ))+ρ−2/3⁢v⁡(ϕ)⁢sin⁡(ρ⁢f~⁢(ϕ))⁢Ai′⁡(−ρ2/3⁢Δ⁢(ϕ))⁢(1+O⁡(1/ρ))),
ρ→∞.

See Figure 36.13.1.

See accompanying text
Figure 36.13.1: Kelvin’s ship wave pattern, computed from the uniform asymptotic approximation (36.13.8), as a function of x=ρ⁢cos⁡ϕ, y=ρ⁢sin⁡ϕ. Magnify

For further information see Lord Kelvin (1891, 1905) and Ursell (1960, 1994).