- #1

- 267

- 2

**Flattening of a planet**

Hi all,

Does anyone know how to calculate the flattening of a planet, say a gaseous one for simplicity, simply from it's radius and rotational period?

Natski

Last edited:

You are using an out of date browser. It may not display this or other websites correctly.

You should upgrade or use an alternative browser.

You should upgrade or use an alternative browser.

- Thread starter natski
- Start date

- #1

- 267

- 2

Hi all,

Does anyone know how to calculate the flattening of a planet, say a gaseous one for simplicity, simply from it's radius and rotational period?

Natski

Last edited:

- #2

- 12,134

- 161

Grav. potential:

PE_grav = - G M1 M2 / r

where

M1 is the planet's mass and M2 is a "test mass"

r is the distance from the planet's center to the surface.

Centrifugal force:

M2 w^2 r cos(latitude)

where

M2 is the same test mass as before

w is the rotation rate in rad/sec

From the cent. force, we get the Centifugal potential is

PE_cent = -(1/2) M2 w^2 r^2 cos(latitude)^2

Set the total potential equal to a constant, and solve for r as a function of latitude:

PE_grav + PE_cent = constant.

Some hints:

1. To find the constant, at least approximately, neglect PE_cent and plug in the known mass and average r in the PE_grav expression.

2. Looks like you'll get a cubic equation to solve. Ouch! This approximation might be helpful:

|PE_cent| << |PE_grav|

Not that you'd neglect PE_cent. But you might be able to do some sort of series expansion, and drop higher-order terms.

- #3

stevebd1

Gold Member

- 749

- 41

[tex]\textit{f}=\frac{a-b}{a}=1-b:a \approx\frac{3\pi}{2GT^{2}\rho}[/tex]

where a is the long radius (equatorial), b is the short radius (polar) G is the gravitational constant, T is the rotation period (in seconds) and [tex]\rho[/tex] is density

http://en.wikipedia.org/wiki/Equatorial_bulge

regards

Steve

- #4

- 267

- 2

Thanks for that Steve, it is a shame Wikipedia has no reference for that equation...

- #5

- 12,134

- 161

I can post it if anybody is interested, but it'll be ascii math as I'm not yet proficient in Latex.

- #6

stevebd1

Gold Member

- 749

- 41

https://www.physicsforums.com/showthread.php?t=231579

I got it to a level where it was a relatively straight forward trial and error process, then it turned into a cubic equation.

Redbelly, feel free to post your method though I suggest you try latex, it's a lot better than using type set and quicker to understand.

Steve

- #7

- 12,134

- 161

I'm about halfway through a Latex version of my derivation, but have to log off for the night and will pick it up again tomorrow.

By the way, I inadvertently posted an incomplete version that appeared here briefly; I've deleted that post until the complete derivation is ready.

By the way, I inadvertently posted an incomplete version that appeared here briefly; I've deleted that post until the complete derivation is ready.

Last edited:

- #8

- 267

- 2

- #9

- 12,134

- 161

The density need not be uniform, it can vary with distance from the planet's center.

I use a a reference frame rotating with the planet, so that the forces on a test mass**m** are the gravitational force and a centrifugal force. I find the potential energy (of the test mass) due to these forces, which is a function of:

*r*, the distance from the planet's center

The latitude.

The planet's rotation rate and mass.

Setting the potential equal to a constant defines an equipotential surface, which must be the case for the planet's surface (if it's not, then material would flow or fall "downhill", in order to make it so).

By solving the potential=constant equation we get the center-to-surface distance as a function of latitude, and can use that to determine the ratio of polar/equatorial radii.

As I'm using Latex for the first time, for me it's slow going typing in the symbols.

edit added:

Time for me to head in to work soon, tonight I should have time to finish it up.

I use a a reference frame rotating with the planet, so that the forces on a test mass

The latitude.

The planet's rotation rate and mass.

Setting the potential equal to a constant defines an equipotential surface, which must be the case for the planet's surface (if it's not, then material would flow or fall "downhill", in order to make it so).

By solving the potential=constant equation we get the center-to-surface distance as a function of latitude, and can use that to determine the ratio of polar/equatorial radii.

As I'm using Latex for the first time, for me it's slow going typing in the symbols.

edit added:

Time for me to head in to work soon, tonight I should have time to finish it up.

Last edited:

- #10

- 12,134

- 161

This is Part 1 (out of 2). Here I'll derive an exact relation between radius and latitude. In part 2 I'll make approximations to get a solution for r.

In a reference frame rotating with the planet, the forces on a test mass

[tex]

F_{grav} = -\frac{GMm}{r^2}[/tex], directed inward toward the planet's

[tex]

F_{cent} = m\omega^2 r \cos(\theta)[/tex], directed outward from the planet's

where

[tex]M[/tex] is the planet's mass

[tex]r[/tex] is the distance from the planet's center to a surface point

[tex]\omega[/tex] is the planet's rotation rate (rad/sec)

[tex]\theta[/tex] is the latitude at a surface point (0 at the equator, 90 degrees at the poles)

The potential energy of the test mass is then

[tex]

U = -\int\vec{F}\cdot\vec{dr}

[/tex]

[tex]

= -\frac{GMm}{r} - \frac{1}{2} m \omega^2 r^2 \cos^2(\theta)

[/tex]

Note we picked up an extra cosine factor in the centrifugal term, because force and radial vector make an angle to one another:

Work = force x displacement x [tex]\cos(\theta)[/tex]

The planet's surface conforms to an equipotential surface. In particular, the potential energy at a pole equals the potential energy anywhere else on the surface:

[tex]

-\frac{GMm}{r_{pole}} = -\frac{GMm}{r} - \frac{1}{2} m \omega^2 r^2\cos^2(\theta)

[/tex]

Multiply both sides of this equation by the factor

[tex]-\: \frac{r}{GMm}[/tex] :

[tex]

\frac{r}{r_{pole}}= 1 + \frac{r^3 \omega^2}{GM}\cos^2(\theta)

[/tex]

or

[tex]

\frac{r}{r_{pole}}= 1 + (\frac{r}{r_{pole}})^3 \cdot\frac{r_{pole}^3\omega^2}{GM}\cos^2(\theta)

[/tex]

This is an exact equation, using dimensionless variables and parameters:

[tex]

\frac{r}{r_{pole}}[/tex] ,

[tex]

\frac{r_{pole}^3\omega^2}{GM}[/tex] ,

and

[tex]

\theta

[/tex]

While cubic equations are analytically solvable, the formula involved is quite cumbersome. In part 2 I'll use the approximation [tex]r \approx r_{pole}[/tex] to get a solution for r.

- #11

- 12,134

- 161

- #12

- 12,134

- 161

[tex]

\frac{r}{r_{pole}}= 1 + (\frac{r}{r_{pole}})^3 \cdot \frac{1}{2} \left[\frac{r_{pole}^3\omega^2}{GM}\right] \cdot \cos^2(\theta)

[/tex]

Part 2 of derivation

We're after solutions for r in the above equation. We assume that r does not vary too much over the planet's surface, in other words

[tex]

\frac{r}{r_{pole}} = 1 + \delta

[/tex]

where

[tex] \delta \ll 1 [/tex]

Our equation becomes, to 1st order in [tex] \delta [/tex] :

[tex]

(1+\delta) \approx 1 + (1+3\delta) \cdot \frac{1}{2} \left[\frac{r_{pole}^3\omega^2}{GM}\right] \cdot \cos^2(\theta)

[/tex]

or

[tex]

\delta \approx (1+3\delta) \cdot \frac{1}{2} \left[\frac{r_{pole}^3\omega^2}{GM}\right] \cdot \cos^2(\theta)

[/tex]

Evidently the term in square brackets is of order [tex] \delta [/tex] , which becomes of order [tex] \delta^2 [/tex] upon multiplication by [tex] 3\delta [/tex] . As we are neglecting terms of order [tex] \delta^2 [/tex] and higher, we may drop the [tex] 3\delta [/tex] from the equation so that

[tex]

\delta \approx \frac{1}{2} \left[\frac{r_{pole}^3\omega^2}{GM}\right] \cdot \cos^2(\theta)

[/tex]

and

[tex]

\frac{r}{r_{pole}} = 1 + \delta \approx 1 + \frac{1}{2} \left[\frac{r_{pole}^3\omega^2}{GM}\right] \cdot \cos^2(\theta)

[/tex]

.

Part 2a. Comparison with the Wikipedia result.

Wikipedia gives the ratio of polar to equatorial radii as

[tex]

1-\frac{r_{pole}}{r_{equator}} \approx \frac{3\pi}{2GT^{2}\rho}

[/tex]

or

[tex]

\frac{r_{pole}}{r_{equator}} \approx 1 - \frac{3\pi}{2GT^{2}\rho}

[/tex]

where [tex] \rho [/tex] is the planet's density and T is the rotation period.

Our equation gives, using [tex] \theta = 0 [/tex] degrees for the equatorial radius,

[tex]

\frac{r_{equator}}{r_{pole}} \approx 1 + \frac{1}{2} \left[\frac{r_{pole}^3\omega^2}{GM}\right]

[/tex]

We make substitutions

[tex]

\omega = 2 \pi / T \text{ and } \frac{M}{r_{pole}^3} \approx \frac{4\pi\rho}{3}

[/tex]

(The expression involving density is valid in the approximation [tex]r \approx r_{pole}[/tex] everywhere on the surface.)

and get

[tex]

\frac{r_{equator}}{r_{pole}} \approx 1 +\frac{3\pi}{2GT^{2}\rho}

[/tex]

So that

[tex]

\frac{r_{pole}}{r_{equator}} \approx \frac{1}{1 +\frac{3\pi}{2GT^{2}\rho}}

[/tex]

or

[tex]

\frac{r_{pole}}{r_{equator}} \approx 1 - \frac{3\pi}{2GT^{2}\rho}

[/tex] ,

in agreement with the result given in Wikipedia.

Share: