Keplerian Orbits: Does a=v^2/r Still Apply?

[DiamondGeezer]

In a Keplerian orbit (an ellipse), the acceleration vector always points to the focus containing the central mass. Since that means the acceleration vector is no longer orthogonal to the orbit (except at two points), does the relation a = \frac{v^2}{r} still apply if we consider v to be the velocity orthogonal to the radius vector?

 

[Ben Niehoff]

There is no a priori reason that it should, although it might by coincidence.

Keplerian orbits follow the areal law: the radius vector sweeps out equal areas in equal times.

It should take only a bit of algebra to show how this relates to the velocity orthogonal to the radius.

 

[D H]

does the relation a = \frac{v^2}{r} still apply if we consider v to be the velocity orthogonal to the radius vector?

In general, no. At two very specific locations on the orbit, yes.

 

[DiamondGeezer]

What I really require is an idea of how tangential velocity and radial acceleration are related in Keplerian orbits. Any clues?

 

[D H]

Is this a homework problem, or you you simply curious?

 

[DiamondGeezer]

I am curious, and no its not homework. I'm working on a problem of my own using data from JPL Horizons and I'm a little stuck on this issue.

 

[D H]

OK.

First, some nomenclature. You mentioned tangential velocity in post #4. The velocity vector is tautologically tangent to the orbit, so the tangential velocity is just the velocity. You meant the transverse (not tangential) component of the velocity vector.

(I'm having troubles with LaTeX, so I'll have to make do with (yuck) HTML math. Bear with me.)

In polar coordinates, the radial vector from the Sun to a planet at a radial distance r from the Sun and a true anomaly θ is

r = r rhat

where rhat is the unit vector

rhat = cosθ xhat + sinθ yhat

Here xhat points in the direction of perihelion and yhat is normal to xhat such that the velocity vector is directed along the yhat axis when the planet is at perihelion.

The velocity vector in polar coordinates is

v = r θdot θhat + rdot rhat

You are looking for the transverse component of velocity, v_t = r θdot. This is just h/r, where h=r² θdot. h is the specific orbital angular momentum (angular momentum divided by the reduced mass), which is a constant of the orbit:

h² = α G (M_s+M_p)

where α is the semi-latus rectum of the orbit. The semi-latus rectum α, the semi-major axis a, and eccentricity e are related via α=a(1-e²).

Dividing by r²,

h²/r² = α G(M_s+M_p)/r²

The left hand side is the square of the transverse component of velocity while the right hand side is the product of the semi-latus rectum and the relative acceleration. In short,

v_t² = α acc

 

[DiamondGeezer]

I need to clarify three things. In Keplerian motion:

1. The acceleration a is the radial acceleration.

2. The transverse velocity v_{trans} is normal to the orbit at any point in the orbit

3. At any point in the orbit, the fraction \frac{v_{trans}^2}{a} is always equal to the semi-latus rectum.

Correct?

 

[D H]

The transverse velocity is the component of the velocity vector normal to the radius vector, not to the orbit. Other than that, yes, you have it correct.

 

[DiamondGeezer]

I'm just recopying your reply into tex so I can see the math in its proper form

In polar coordinates, the radial vector from the Sun to a planet at a radial distance r from the Sun and a true anomaly θ is

\hat r = \frac{\vec{r}}{|r|}

where \hat r is the unit vector

\hat r = \cos(\theta ) \hat x + \sin(\theta ) \hat y

Here \hat x points in the direction of perihelion and \hat y is normal to \hat x such that the velocity vector is directed along the \hat y axis when the planet is at perihelion.

The velocity vector in polar coordinates is

\vec{v} = \vec{r}\dot{\theta} \cdot \hat{\theta} + \vec{r} \cdot \hat r

You are looking for the transverse component of velocity, v_t = \vec{r} \dot{\theta}. This is just h/r, where h=r^2 \dot{\theta}. h is the specific orbital angular momentum (angular momentum divided by the reduced mass), which is a constant of the orbit:

h^2 = \alpha G (M_s+M_p)

where \alpha is the semi-latus rectum of the orbit. The semi-latus rectum \alpha, the semi-major axis a, and eccentricity e are related via \alpha=a(1-e^2).

Dividing by r^2,

\frac{h^2}{r^2} = \alpha \frac{G(M_s+M_p)}{r^2}

The left hand side is the square of the transverse component of velocity while the right hand side is the product of the semi-latus rectum and the relative acceleration. In short,

v_t^2 = \alpha a

[Hopefully this has been transcribed correctly! Please check!]

 

[DiamondGeezer]

The transverse velocity is the component of the velocity vector normal to the radius vector, not to the orbit. Other than that, yes, you have it correct.

Doh! Because the velocity normal to the velocity vector would be zero

Let me capture your reply in vector notation

v_t = \vec{v} \times \hat r

 

[DiamondGeezer]

The left hand side is the square of the transverse component of velocity while the right hand side is the product of the semi-latus rectum and the relative acceleration. In short,

v_t^2 = \alpha a

[Hopefully this has been transcribed correctly! Please check!]

Sorry this last part should be

v_t^2 = \alpha a_{cc}

 

[D H]

In polar coordinates, the radial vector from the Sun to a planet at a radial distance r from the Sun and a true anomaly θ is

\hat r = \frac{\vec{r}}{|r|}

As written, this defines the unit vector rhat. You are given r and θ. You don't know the vector \vec r. This equation is much better expressed as

\vec r = r \hat r

\vec{v} = \vec{r}\dot{\theta} \hat \theta + \vec{r} \cdot \hat r

No. \hat{\theta} is a vector. That should be a scalar r coupled with \hat{\theta}, not the vector \vec r. This got you in trouble in several places. Note that it is customary to denote the scalar magnitude of some vector \vec q with the corresponding symbol sans the vector glyph: q\equiv||\vec q||. That is what is being done here with r.

The second term is even worse. It has units of length, not velocity.

Apply the product rule to \vec r = r \hat r. Doing so yields

\vec v = \dot r \hat r + r\frac{d}{dt}(\hat r)

Now differentiate the unit vector rhat:

\frac{d}{dt}\hat r =<br /> -\sin\theta\,\dot{\theta}\hat x + \cos\theta\,\dot{\theta}\hat y<br /> = \dot{\theta}\hat\theta

Collecting terms yields the expression for the velocity vector in polar coordinates:

\vec{v} = r\dot{\theta} \hat \theta + \dot r \hat rFinally, you used the symbol a to represent both the semi-major axis and the acceleration. Using one symbol to represent two different things is never a good idea.

Let me capture that in vector notation

v_t = \vec{v} \times \hat r

No!

That vector is normal to the orbital plane. Think of it this way: The specific angular momentum vector is

\mathbf h = \mathbf r \times \mathbf vTo get the component of some vector q normal to some unit vector u you need to subtract the projection of q onto u:

<br /> \mathbf q_{\perp} =<br /> \mathbf q - (\mathbf q\cdot \hat{\mathbf u}) \hat{\mathbf u}<br />

Alternatively, you can use the vector triple product:

<br /> \mathbf q_{\perp} =<br /> \hat{\mathbf u} \times (\mathbf q \times \hat{\mathbf u}) =<br /> (\hat{\mathbf u} \times \mathbf q) \times \hat{\mathbf u}<br />

Note that in general the vector triple product is not associative. It is in this case because the unit vector uhat appears twice.

Using the above, the transverse component of the velocity vector is

<br /> \aligned<br /> \mathbf v_{\text{trans}} &amp;=<br /> \mathbf v - (\mathbf v\cdot \hat{\mathbf r}) \hat{\mathbf r} \\<br /> &amp;= \hat{\mathbf r} \times (\mathbf v \times \hat{\mathbf r}) \\<br /> &amp;= (\hat{\mathbf r} \times \mathbf v) \times \hat{\mathbf r}<br /> \endaligned<br />

 

[DiamondGeezer]

I'm glad I copied and asked, otherwise I'd have been going in the wrong direction (or lots of them).

The vector triple product (perversely) makes more sense to me than "subtracting the projection of...". I must make more use of the right hand rule to prevent mistakes!

I know about the confusing of the symbols for semi-major axis and acceleration - hence my later post.

I must admit that my knowledge of Keplerian orbits is pretty rudimentary up to this point, so there's plenty more to learn.

DH - this is plenty for me to digest and work on. In case this thread goes a little quiet from here on in, can I thank you for your attention to my problem and for your prompt responses?

 

[D H]

The vector triple product (perversely) makes more sense to me than "subtracting the projection of...".

That is kinda perverse. Look at it from the perspective of computational cost. "Subtracting the projection of ..." involves an inner product (3 multiplies, 2 adds), a scalar times a vector (3 multiplies), and a vector subtraction (3 adds), for a total of 6 multiplies and 5 adds. One cross product requires 6 multiplies and 3 adds. The vector triple product involves two cross products or 12 multiplies and 6 adds. The projection method, with half the number of multiplication operations as the vector triple product method, is nearly twice as fast as using the vector triple product -- and is also a lot less prone to numerical error.

DH - this is plenty for me to digest and work on. In case this thread goes a little quiet from here on in, can I thank you for your attention to my problem and for your prompt responses?

You're welcome! Glad to help.

 

[DiamondGeezer]

DH

Sorry to bother but I need some clarification. I'm using MathCAD (hence I can crunch vector data quickly) but still need brains to do the calculations correctly.

I'm using data from JPL Horizons for the planet Mercury (as an example)

My calculations run like this:

At a point xi,yi,zi from the center of the Sun the velocity vectors for Mercury are vxi, vyi, vzi

Let \vec{r_i}=\vec{x_i},\vec{y_i},\vec{z_i} and \vec{v_i}=\vec{vx_i},\vec{vy_i},\vec{vz_i}

Therefore the unit radius vector is \hat r_i = \frac{r_i}{|r_i|}

This means that the radial velocity is vrad_i= \hat r_i \cdot \vec{v_i}

To get the vector normal to the plane of the orbit:

\vec{vn_i}= \hat r_i \times \vec{v_i} (This should be zero)

Then the transverse vector should be \vec{vt_i}= ( \hat r_i \times \vec{v_i}) \times \hat r_i ?

In which case I need the magnitude of \vec{vt_i}, square it, divide it by the radial acceleration and that should be the length of the semi-latus rectum (ie a constant)?

 

[D H]

This means that the radial velocity is vrad_i= \hat r_i \cdot \vec{v_i}

Correct.

To get the vector normal to the plane of the orbit:

\vec{vn_i}= \hat r_i \times \vec{v_i} (This should be zero)

That this cross product is normal to the instantaneous orbital plane is correct. That this should be zero is absolutely incorrect. This is the specific angular momentum vector. If the planets truly did obey Kepler's laws (and they don't), this vector would be a constant.

Then the transverse vector should be \vec{vt_i}= ( \hat r_i \times \vec{v_i}) \times \hat r_i ?

Correct.

In which case I need the magnitude of \vec{vt_i}, square it, divide it by the radial acceleration and that should be the length of the semi-latus rectum (ie a constant)?

That would be a constant if the planets truly did obey Kepler's laws. Simple Keplerian orbits can be characterized by six orbital elements at some epoch time. Only of of these six elements, the true anomaly, changes with respect to time. The other orbital elements are all constants of the orbit.

The planets do not exactly obey Kepler's laws. Kepler's laws are only approximately correct. The planets perturb each others' orbits. In addition, even Newton's law of gravitation is not quite correct; general relativity yields a better (but much, much hairier) model. The deviation between general relativity and Newtonian gravity is the greatest for Mercury because of Mercury's relatively high orbital velocity.

The argument of perhielion is one of those constant orbital elements in a simple Keplerian orbit. Mercury's perihelion advances by 531 arcseconds per century because of the influence of the other planets and by 43 arcseconds per century because of general relativistic effects.