On the dispersion relation for trapped internal waves

Abstract

An analysis is constructed in order to estimate the dispersion relation for internal waves trapped in a layer and propagating linearly in a fluid of infinite depth with a rigid surface. The main interest is in predicting the structure of internal wave wakes, but the results are applicable to any internal waves. It is demonstrated that, in general 1/cp = 1/CpO + k/ωmax + ∈(k) where cp is the wave phase speed for a particular mode, CpO is the phase speed at k = 0, ωmax is the maximum possible wave angular frequency and ωmax ≤ Nmax where Nmax is the maximum buoyancy frequency. Also, ∈(0) = 0, ∈(k) = o(k) for k large, and is bounded for finite k. In particular, when ∈(k) can be neglected, the dispersion relation for a lowest mode wave is approximately 1/cp ≈ (∫∞0N2(y)ydy)-½ + k/ωmax. The eigenvalue problem is analysed for a class of buoyancy frequency squared functions N2(x) which is taken to be a class of realvalued functions of a real variable x where O ≤ x ∞ such that N2(x) = O(e-βx) as x → ∞ and 1/β is an arbitrary length scale. It is demonstrated that N2(x) can be represented by a power series in e-βx. The eigenfunction equation is constructed for such a function and it is shown that there are two cases of the equation which have solutions in terms of known functions (Bessel functions and confluent hypergeometric functions). For these two cases it is shown that ∈(k) can be neglected and that, in addition, ωmax = Nmax. More generally, it is demonstrated that when k → ∞ it is possible to approximate the equation uniformly in such a way that it can be compared with the confluent hypergeometric equation. The eigenvalues are then, approximately, zeros of the Whittaker functions. The main result which follows from this approach is that if N2(x) is O(e-βx) as x → ∞ and has a maximum value N2max then a sufficient condition for 1/cp ∼ k/Nmax to hold for large k for the lowest mode is that N2(t)/t is convex for O ≤ t ≤ 1 where t = e-βx.