6. Internal waves
Consider a continuously stratified fluid with ρo(z) the vertical density profile.
z ρ p' ξ p
ρο(z)
Figure by MIT OpenCourseWare. Figure 1.
At a point P raise a parcel of water by small amount ξ from its equilibrium position P
adiabatically. The change in pressure experienced by the parcel is
dp = -ρo gξ
dp if ρ (z) is the medium density, while the change in density is dρ= [Remember the o 2 cs
⎛ ⎞ ∂p 2 definition of adiabatic compressibility ⎜ ⎟ = cs ]. The buoyancy force acting on the ⎝ ∂ρ ⎠
particle induces an acceleration because the density of the particle at the displaced
position is different from the background density of the medium
dρo ρo |P'=ρo |P + ξ and dz
−ρogξ ρ | = ρ |particle at + ρ ρo particle at P' P d = 2 cs
1 d2ξ ρ gξ ρ = g[(ρ +ρ ξ) − (ρ − o )] o 2 o oz o 2 dt cs
d2ξ g ρ g = [(ρ + o )ξ] 2 ρ oz 2 dt o cs
d2ξ gρ g2 + (− oz − )ξ=0 2 ρ 2 dt o cs
−g dρ g2 Define N2(z) = o − the Brünt-Vaisala frequency or buoyancy frequency ρ 2 o dz cs
2 d2ξ + N2ξ=0 dt2 is the equation of the harmonic oscillator with solution e±iNt.
Thus the parcel oscillates about its equilibrium position and the Brünt-Vaisala frequency is the natural frequency of oscillation determined by the local density stratification and the fluid’s compressibility. As we saw, in the ocean compressibility effects are negligible and we can assume
−g ∂ρ N2 = o ρo ∂z
The motion is the oscillation of the particle around the equilibrium position with frequency N. Let us compare the frequency of surface gravity waves with the buoyancy frequency and use deep water gravity waves
ω N 1 ∂ρ 1 ∂ρ Δρ int . = = (− o )1/2~(− o λ)1/2 ≈ o( o )1/2 ωsurf. gk kρo ∂z ρo ∂z ρo
1 ∂ρ k~ ; λ o =Δρ density change over λ λ ∂z o
2 −3 Δρo ~10 over total depth:
ωint << ωsurface .
The frequency of internal waves is much less than the frequency of surface gravity waves. The restoring force for surface gravity waves is g; the restoring force for
Δρ internal waves is o g, reduced gravity. ρo
Internal waves in an unbounded fluid
Consider an incompressible, stratified non-rotating fluid characterized by a basic state of rest and hydrostatic equilibrium:
dp u ≡ o, ρ=ρ (z), p = p (z) and o =−gρ o o dz o
The fluid experiences small amplitude perturbations around the basic state, so we can use the linearized equations of motion, assuming also that the motion is frictionless and adiabatic:
∂u ∂p (1) ρ =− ρtotal = ρo+ρ; ρ<<ρo o ∂t ∂x
∂v ∂p (2) ρ =− ptotal = po+ρ; p< ∂w ∂p (3) ρ =− −ρg ρ = perturbation density o ∂t ∂z ∂u ∂v ∂w (4) + + = 0 p = perturbation pressure ∂x ∂y ∂z ∂p ∂ρ (5) + w o = 0 The basic reference state cancels out in 3 ∂t ∂z If u = v = o and p = o, we have the special solution 3 ∂w ∂ρ ∂ρ ρ =−ρg; =−w o o ∂t ∂t ∂z ∂2w ∂ρ g∂ρ ∂2w −g ∂ρ ρ =−g =+ o w → + ( o )w = o o 2 2 ∂t ∂t ∂z ∂t ρo ∂z General equation ∂ ∂u ∂ ∂v ∂ ∂w Differentiate (4): + + = 0 ∂t ∂x ∂t ∂y ∂t ∂z ∂u ∂v Substitute into it ( , ) given by (1) and (2) ∂t ∂t ∂ 1 ∂p ∂ 1 ∂p ∂ ∂w (− ) + (− ) + ( ) = 0 ∂x ρo ∂x ∂y ρo ∂y ∂z ∂t ∂2w 1 ∂2 ∂2 (I) = ∇2 p where ∇2 = + H H 2 2 ∂z∂t ρo ∂x ∂y Eliminate ρ between (3) and (5) 2 2 2 ∂ w ∂ p ∂ρ ∂ p dρo ρo =− − g =− + gw ∂t2 ∂z∂t ∂t ∂z∂t ∂z ∂2w g ∂ρ 1 ∂2p or + (− o )w =− 2 ∂t ρo ∂z ρo ∂z∂t ∂2w 1 ∂2p (II) + N2w =− 2 ∂t ρo ∂z∂t 2 Eliminate p between (I) and (II) taking ∇H of (II) ∂2 1 ∂2 1 ∂2 ∂ ∂w ∇2 ∇2 w + N2∇2 w =− ∇2 p =− (ρ ) H 2 H H 2 H 2 o ∂t ρo ∂z∂t ρo ∂t ∂z ∂z ∂2 Take to the LHS and factorize ∂t2 4 ∂2 ⎡ 1 ∂ ∂w ⎤ Internal wave equation ∇2 w + (ρ ) + N2∇2 w = 0 2 ⎢ H o ⎥ H ∂t ⎣ ρo ∂z ∂z ⎦ Consider the last term in square bracket: 1 ∂ ∂w 1 ∂ρ ∂w ∂2w (ρ ) = o + o 2 ρo ∂z ∂z ρo ∂z ∂z ∂z Ratio ∂ρo ∂w 1 ∂ρ 1 d ∂ρ d ∂z ∂z = o • 2 = o << = 2 d 1 ∂ w ρo ∂z d ρo ∂z D ρo ∂z2 d = vertical scale of w ~ thickness of thermocline and D is the ocean depth or atmospheric height -> valid also in the atmosphere. So we can ignore the first term and we simplify to: 2 ⎡ 2 2 2 ⎤ ⎛ 2 2 ⎞ ∂ ∂ ∂ ∂ 2 ∂ ∂ ⎢ + + ⎥w + N ⎜ + ⎟w = 0 ∂t2 ⎣⎢ ∂x2 ∂y2 ∂z2 ⎦⎥ ⎝ ∂x2 ∂y2 ⎠ total Laplacian horizontal Laplacian 2 ∂ 2 2 2 ∇ w + N ∇Hw = 0 internal wave equation. ∂t2 Notice that if N = 0, no stratification, ∇2 ω = 0 the motion is irrotational -> surface gravity waves. Consider now a plane wave solution in three-dimensions w = wocos(kx+ly+mz-ωt) (k,l,m) wave number Substituting into the internal wave equation, we obtain the dispersion relationship N2(k2 + l2) ω2 = (k2 + l2 + m2) 5 N2k2 or ω2 = H K2 k ω=±N H K (m,z) K m θ (l,y) ϕ (k,x) KH Figure by MIT OpenCourseWare. Figure 2. ω = ±N cos θ And the wave numbers are k = Kcosθcoφ,l= Kcosθsinφ; m = Ksinθ The dispersion relation for internal waves is of a quite different character from the dispersion relation for surface waves: ω of surface waves depends only on the magnitude of K and not on its direction ω of internal waves is independent of the magnitude of K and depends only on the orientation of the wave vector, i.e. on the angle θ of the wavenumber K with the horizontal. 6 Consider the displacement ζ along a line of constant phase by definition perpendicular to the phase line z K θ x θ Fz θ dz θ Figure 3. Figure by MIT OpenCourseWare. Thus particle motion is along the wave crests i.e., along the lines of constant phase d2ζ + N2ζ=0 dt2 We can see this also considering i(kx+ly+mz−ωt) (u,v,w) = (uo,v o,w o)e Continuity equation gives kuo+lvo+mwo = 0 => K •u = 0 Fluid velocity is perpendicular to the wave vector, i.e. is along the crests of the waves => wave motion is transverse 7 K u θ u 3 high u θ2 low θ 1 high phase propagation Figure 4. Figure by MIT OpenCourseWare. When the wave vector is horizontal, m = o, the particle motion is purely vertical and ωmax = N A displacement ζ along the phase line gives a vertical displacement dz = ζcosθ The buoyancy force in the vertical corresponding to the displacement dz is 2 2 Fz = -N dz = -N ζcosθ The component of this force along the phase line is 2 2 Fζ = -N ζcos θ and the motion of the particle along the phase line is governed by: 2 2 d ζ d ζ 2 2 = Fζ → + (N cos θ)ζ=0 dt2 dt2 Thus ±Ncosθ is the frequency of oscillation of a particle along the phase line. If θ = 0, we have vertical oscillations with frequency N. 8 Dispersion effects For internal waves, the surfaces of constant frequency in wavenumber space are the cones θ = constant, i.e. ω = constant ω N The phase velocity |c|= = cosθ is directed along the wave vector and therefore lies K K ∂ω ∂ω ∂ω on the cone. The group velocity c = i + ˆj + kˆ is the gradient of ω in g ∂k ∂l ∂m wavenumber space and therefore by definition is perpendicular to the surfaces ω = constant. It follows that c g is at right angles with ( K,c ). When the group velocity has an upward component the phase speed has a downward component. Let us explicitly evaluate the group velocity k2 + l2 N2k2 Write ω2 = N2 = H k2 + l2 + m2 K2 ∂ω N m2 k N = = sinθ(sinθcosφ) 2 ∂k K K kH K ∂ω N m2 l N = = sinθ(sinθsinφ) 2 ∂k K K kH K ∂ω k m N =−N H =− sinθ(cosθ) ∂m K3 K Note ω ∂ω N2 =− cos2 θ<0 m ∂m K2 9 The vertical phase velocity is always opposite to the vertical group velocity. Waves propagating their phase upwards will be propagating their energy downwards and viceversa, N or c = sinθ(sinθcosφ,sinθsinφ,−cosθ). g K Explanation of the Figure For internal waves, wave crest AA’ has moved to 4 wavelengths downward to the left, the wave group has moved upward parallel to the crests, at right angle to the phase propagation.