applet-magic.com
Thayer Watkins
Silicon Valley
& Tornado Alley
USA

The Angular Momentum of
a Rankine Vortex on a Sphere

Rankine Vortices

One approach to modeling tropical cyclones is as rankine vortices. A Rankine vortex is one in which the inner part turns as a disk and in the outer part the velocity is inversely proportional to the distance from the center; i.e.,


v(r) = va(r/a) for r≤a
v(r) = va/(r/a) for r≥a
 

where va is the tangential velocity at r=a. In a cyclone a would be the radius of the eyewall.

In computing the angular momentum in a Rankine vortex each band of width Δr in the outer portion has the same angular momentum so a Rankine vortex that extends to infinity would have an infinite angular momentum. Or, stated another way, the angular momentum of a Rankine vortex is directly proportional to the outer limit of the vortex even though the velocity at the outer limit may be negligibly small. But tropical cyclones are on the sphere of the Earth and therefore necessarily of finite extent.

A Rankine Vortex on a Sphere

Let θ be the angle from the center of the vortex to a point a distance r from the center and let R be the radius of the sphere. The distance r is then Rθ and hence the formula for the velocities in a Rankine vortex on a sphere is:


v(θ) = va(Rθ/a) for Rθ≤a
v(θ) = va/(Rθ/a) for Rθ≥a
 

The angular momentum of an infinitesimal element with respect to the spin axis of the vortex is the areal density ρ times the velocity times the distance s of the element from the spin axis. The distance s of a point at angle θ from the spin axis is Rsin(θ). The area of a band of with dθ centered at an angle θ is (2πs)Rdθ=(2πRsin(θ)dθ.

The angular momentum L is then the integral of the angular momenta in bands of width Rdθ; i.e.,


L = ∫0π(ρv(θ)Rsin(θ))(2πRsin(θ))Rdθ
= 2πρR30πv(θ)sin2(θ)dθ
 

In meteorology quantities are usually stated per unit mass. The mass of the vortex is ρ(4πR2) so the angular momentem of the vortex per unit mass is


L/M = (R/2)∫0πv(θ)sin2(θ)dθ
 

The integral must be evaluated for the inner and outer parts of the vortex separately. For the inner part of the vortex the angle θ ranges from 0 to a/R so the integral is


0a/Rva(Rθ/a)sin2(θ)dθ
 

The indefinite integral of


θsin2(θ)
is
θ2/4 - (θ/4)sin(2θ)-(1/8)cos(2θ)
 

so the angular momentum of the inner portion of the vortex evaluates to


(R/2)va(R/a)[θa2/4 - (θa/4)sin(2θa)-(1/8)cos(2θa) + (1/8)]
where θa=a/R.
 

The evaluation of the angular momentum for the outer part of the vortex involves the integral of sin2(θ)/θ. The indefinite integral of this quantity is


ln(θ)/2 - [cos(2θ) - 1]
 

Thus the angular momentum per unit mass for the outer portion of the vortex evaluates to


(R/2)va(a/R)[ln(π)/2 - ln(θa)/2 + (cos(2θa)-1)]
= va(a/2)[ln(π/θa)/2 + (cos(2θa)-1)]
 

where again θa is the angular distance of the eyewall of the vortex from the center of the sphere; i.e., a/R.

Representative Values
of the Parameters and
Values for Cyclones

In SI units the typical values for a hurricane are:


HOME PAGE OF applet-magic
HOME PAGE OF Thayer Watkins