Tuesday, June 26, 2012

Derivation of the day: scale height

Background
The scale height, H, refers to the height above the plane of a disk of particles surrounding a larger body.  Due to eccentricities or inclinations of orbits from pressure within the system, these disks are not completely flat and each particle has an associated "random velocity" defined as u = ivorb or u = evorb.  We find this height by balancing forces within the disk.

Units
u = random velocity
Ω = angular frequency
G = Newtonian gravitation constant
a = radius of orbit

Derivation
Let us consider a particle with an inclination given by i = arcsin(H/a).  Starting with the ideal gas law,

P = nkT ~ nmv2 = mv2/H(Area)

where n is a volume density. Force = (Pressure)(Area) so we can write that the force pushing the particle up is approximately:

mv2/H

We also consider the downward force due to the gravitation of the star:

Fg = (GMm/a2)sin(i) = (GMm/a2)(H/a)

From Kepler's Law, we can substitute GM = Ω2a3 so

mv2/H = Ω2mH

Rearranging, we find an expression for H:

H = u/Ω

No comments:

Post a Comment