Deriving an equation for orbital period

In summary: I was hoping you could derive from your equation Kepler's.In summary, the conversation discusses an equation for determining orbital period given altitude, speed, and mass of the orbiting bodies. The equation is derived using the specific orbital energy and properties of a circular orbit. It is then compared to Kepler's 3rd law and found to be valid for circular orbits, but not for highly eccentric orbits. The conversation also highlights the importance of Kepler's genius in creating an equation that is applicable to all planets in the solar system, regardless of their eccentricities.
  • #1
cephron
124
0
Hi guys, I derived an equation for determining orbital period, given an altitude, speed, and mass of the primary and the object orbiting it. I think it makes sense, but I'd welcome anyone who is willing to check it for conceptual error or nonsensical math.

Here is the equation:

[tex]P = \frac{2 \pi ~ \mu }{ \Big( \frac{2 \mu}{r} - v^2 \Big) ^{ \frac{3}{2}} }[/tex]
Where
[itex]P[/itex] is the orbital period in seconds,
[itex]\mu[/itex] is the standard gravitational parameter [itex]G(m_1 + m_2)[/itex]
([itex]G[/itex] is the gravitational constant in [itex]m^3~kg^{-1}~s^{-2}[/itex],
[itex]m_1[/itex] and [itex]m_2[/itex] are the masses of the orbiting bodies),
[itex]r[/itex] is the distance between them in meters, and
[itex]v[/itex] is their relative speed in meters per second.

Concept
The idea is that, if an orbiting body has a certain speed at a certain instant, its orbital period is not affected by the direction in which it is headed at that moment in time. It could be taking a circular orbit or a heavily elliptical one, and its period would be the same because the orbit's specific energy would be the same.
So, what this equation does is:
1) Calculate the specific orbital energy using the givens (altitude, speed, m1, m2);
2) From this specific energy, use the properties of the circular orbit to find [itex]v_c[/itex] and [itex]r_c[/itex], the speed and altitude/distance values for the equivalent circular orbit;
3) Calculate the orbital period of this equivalent circular orbit.

Derivation - the pieces
So, we start with the equation to calculate specific orbital energy:
[tex]\epsilon = \frac{1}{2}v^2 - \frac{\mu}{r}[/tex]
And then we will use the properties of a circular orbit to find [itex]v_c[/itex] and [itex]r_c[/itex]:
[tex]\frac{1}{2}v_c^2 = -\epsilon[/tex][tex]\frac{-\mu}{r_c} = 2\epsilon[/tex]which simplify to:
[tex]v_c = \sqrt{-2\epsilon}[/tex][tex]r_c = \frac{-\mu}{2\epsilon}[/tex]
And finally, we'll use [itex]v_c[/itex] and [itex]r_c[/itex] to get the orbital period:
[tex]P = \frac{2\pi r_c}{v_c}[/tex]

Derivation - putting them together
With the circular orbit equations all plugged together, it looks like this:
[tex]P = \frac{-2\pi \frac{\mu}{2\epsilon}}{\sqrt{-2\epsilon}}[/tex]which simplifies to:
[tex]P = \frac{2\pi \mu}{(-2\epsilon)^{\frac{3}{2}}}[/tex]
And then, plugging in the specific orbital energy equation, we get:
[tex]P = \frac{2\pi \mu}{\Big( -2(\frac{1}{2}v^2 - \frac{\mu}{r}) \Big)^{\frac{3}{2}}}[/tex]which simplifies to the final equation at the head of this post.

Behaviour
If the orbiting object is traveling at a speed greater than its escape velocity, the kinetic energy of the orbit is greater than the potential energy. The term inside the (3/2) exponent works out to be negative, and the result is undefined--not a closed orbit, so no period.
If the orbiting object is traveling at exactly escape velocity, the kinetic and potential energy cancel out, leaving a denominator of zero--undefined again.
If the potential energy just barely outweighs the kinetic energy, we have a tiny positive denominator, and so we get a huge orbital period. This seems sane--if we consider an object orbiting at almost its escape velocity, it would probably have a large orbital period.
So yeah...this seems to be a pretty simple derivation, but I'm very noob at this stuff, so please point out if I botched something. And if you look through it and it checks out, please indicate your approval, because I'll probably be wanting to use this in further questions of mine soon. Many thanks!
 
Last edited:
Physics news on Phys.org
  • #2
One simple way to check out your equation, Cephron;

It must replicate Kepler's 3rd law.
IOW's you must be able to transform your equation back to the Kepler eqn. namely; the period must transform back to ...

P = sq.rt.[(4(pi)^2) (R^3)/ GM]

Give it a try ,
... for simplicity assume small orbital mass so you can use u = GM instead of G(M+m), and also substitute v = sq.rt.(GM/R), and then solve for P.

It appears that you may recover Kepler...although I haven't gone through the algebra ...if not then your eqn. is not acceptable classically.

Creator
 
Last edited:
  • #3
Given the substitutions you offered, it appears to check out! Thanks!

I'm not sure how you knew to substitute v for sqrt(GM/r), though...what assumptions are needed for that substitution to work? Would the equation still hold for elliptical orbits if it uses that substitution to get back to Kepler?
 
  • #4
cephron said:
Given the substitutions you offered, it appears to check out! Thanks!

I'm not sure how you knew to substitute v for sqrt(GM/r), though...what assumptions are needed for that substitution to work?

The same assumption applies as we used for u = GM, namely, that the orbiting mass is almost negligible compared to the central mass, but in addition, there is the more important assumption of very little eccentricity (which was the initial constraint with your derivation using circular orbit.)

Why could I use that substitution? Because for circular orbits the kinetic orbital energy always equals -1/2 the potential energy.
IOW;
K.E. =-(1/2)U
1/2mv^2 = - (1/2)GMm/r

and solving for orbital vel.(ignoring the sign convention)...
v = sqrt(GM/r) ...for CIRCULAR orbits. (v is actually the mean orbital speed).

Would the equation still hold for elliptical orbits if it uses that substitution to get back to Kepler?

Not sure, but your equation would probably not be applicable because mean orbital speed actually decreases with increased eccentricity and that would alter the substitution.

Kepler's equation allows you to substitute the semi-major axis of elliptical orbits for the radius (in the eqn.) and you can still use his same equation.
But in your case you derived with the initial assumption of circular orbit, and unfortunately the mean orbital speed v changes with eccentricity.

Fundamentally, the problem arises since you are using a velocity dependent term whereas Kepler uses only period and radius which sweeps out equal areas in equal time periods in spite of a changing radius and thus circumvents the eccentricity problem. Now you know the genius of Kepler.

Furthermore, the ratio of R^3/P^2 is the same value for ALL planets in the solar system, namely GM/(2pi)^2...(where M is solar mass)...in spite of their different orbital eccentricities as long as R is used as the semi-major axis.

Good questions, BTW.

Creator
 
Last edited:
  • #5
http://burtleburtle.net/bob/physics/orbit101.html

I'm not certain how reliable this website is, but I'm guessing from his apparently successful gravity simulations that he knows something about all this.

Anyways, the key thing I noticed on this webpage was where he says:
The position and absolute velocity determine the period of an orbit. Direction does not matter.
and had the little java applet showing the simulation of a bunch of objects departing a point in different directions, but arriving back there at the same time. This is essentially what prompted me to build the above equation: if period is a function of position and speed (and of course mass), then there should be an equation which gives us that.

My idea for how to calculate the period of any elliptical orbit was to calculate the period of a circular orbit with the same specific orbital energy. It looks to me like specific orbital energy is dependant on the exact same things, no more, no less, as orbital period: distance, absolute velocity, and mass of the bodies. This makes me believe--please correct me if I'm wrong!--that there is a one-to-one mapping between orbital period and specific orbital energy. Thus, if you've calculated the orbital period for one orbit of energy [itex]\epsilon[/itex], you've calculated them all, right? This is the assumption upon which my equation is based...if the one-to-one mapping between specific orbital energy and orbital period is true, then I think my equation would work for calculating the period of any elliptical orbit.
 
  • #6
The orbital period of an object can be expressed as [tex]T=2\pi\sqrt{\frac{a^3}{\mu}}[/tex]
where a is the semi-major axis (see Wikipedia for example). This can be related to the orbital energy via [itex]a=\frac{-\mu}{2\varepsilon}[/itex] with the energy [itex]\varepsilon=\frac{v^2}{2}-\frac{\mu}{r}[/itex]. Putting everything together, I think you get your formula.
 
  • #7
Yes, your formula works and it does eliminate the need to calculate the semi-major axis.

I would note one thing. The geocentric gravitational constant (mu) is just G * M, with M being the mass of the Earth and G being the universal gravitational constant.

You're interested in the motion of orbiting objects, and you only need the acceleration due to gravity to figure that out (and since F=ma, all objects accelerate at the same rate regardless of their mass).

That's also why specific energy (actually, specific energy per unit of mass) is used instead of your actual energy level (which does require mass).

Your mu works for a satellite, since the mass of the satellite is miniscule compared to the mass of the Earth, but the mass of the satellite really shouldn't be there.
 
  • #8
cephron said:
[ This makes me believe--please correct me if I'm wrong!--that there is a one-to-one mapping between orbital period and specific orbital energy.

I see what you are trying to do Cephron, but there is no need to derive a new equation. The correspondance between orbital period of an elliptical orbit and specific orbital energy is already established through this other eqn. for total specific orbital energy, e ...

e = -u / 2a ...where a is the semi-major axis.

And since by Kepler's law the Period squared must be proportional to semi-major axis cubed, then the period must also be proprotional to 'e'.

if the one-to-one mapping between specific orbital energy and orbital period is true, then I think my equation would work for calculating the period of any elliptical orbit.

Again , I don't think so. The reason is because in your derivation you constrained the velocity to a circular orbit, using v = sqrt(-2Gm/r), . However, for elliptical orbits the velocity goes as the more general vis-viva equation which is different and must include the semi-major axis...see here...(first equation at top of page)..
http://en.wikipedia.org/wiki/Vis-viva_equation

However, I could be wrong; I haven't done the substitution to see how it comes out.
Sorry it took so long to respond..

...
 
Last edited:

1. What is an orbital period?

An orbital period is the time it takes for an object to complete one full orbit around another object, such as a planet around a star.

2. How do you derive an equation for orbital period?

The equation for orbital period can be derived using Kepler's Third Law, which states that the square of a planet's orbital period is directly proportional to the cube of its average distance from the Sun.

3. What is Kepler's Third Law?

Kepler's Third Law, also known as the "harmonic law", is a mathematical relationship that describes the orbital period of a planet in relation to its distance from the Sun.

4. What factors affect the orbital period of a planet?

The orbital period of a planet is primarily affected by its distance from the central object it is orbiting and the mass of that object. Other factors such as the shape and size of the orbit can also have an impact.

5. Can the equation for orbital period be applied to all objects in orbit?

Yes, the equation for orbital period can be applied to any object in orbit around another object, as long as the mass and distance of the two objects are known.

Similar threads

Replies
3
Views
986
  • Special and General Relativity
Replies
5
Views
1K
  • Advanced Physics Homework Help
Replies
5
Views
1K
Replies
86
Views
4K
Replies
62
Views
5K
Replies
1
Views
2K
Replies
3
Views
951
  • Special and General Relativity
2
Replies
42
Views
3K
  • Advanced Physics Homework Help
Replies
26
Views
3K
Replies
8
Views
1K
Back
Top