View source for Astronomy:Radial trajectory

In astrodynamics and [[Physics:Celestial mechanics|celestial mechanics]] a '''radial trajectory ''' is a [[Astronomy:Kepler orbit|Kepler orbit]] with zero [[Physics:Angular momentum|angular momentum]]. Two objects in a radial trajectory move directly towards or away from each other in a straight line.
{{Astrodynamics}}
== Classification ==

There are three types of radial trajectories (orbits).<ref>Thomson, William Tyrrell; ''Introduction to Space Dynamics'', Dover, 1986</ref>

* Radial elliptic trajectory: an orbit corresponding to the part of a degenerate ellipse from the moment the bodies touch each other and move away from each other until they touch each other again. The relative speed of the two objects is less than the [[Astronomy:Escape velocity|escape velocity]]. This is an elliptic orbit with semi-minor axis = 0 and eccentricity = 1. Although the eccentricity is 1, this is not a parabolic orbit. If the [[Physics:Coefficient of restitution|coefficient of restitution]] of the two bodies is 1 (perfectly elastic) this orbit is periodic. If the coefficient of restitution is less than 1 (inelastic) this orbit is non-periodic.
* Radial parabolic trajectory, a non-periodic orbit where the relative speed of the two objects is always equal to the escape velocity. There are two cases: the bodies move away from each other or towards each other.
* Radial hyperbolic trajectory: a non-periodic orbit where the relative speed of the two objects always exceeds the escape velocity. There are two cases: the bodies move away from each other or towards each other. This is a hyperbolic orbit with semi-minor axis = 0 and eccentricity = 1. Although the eccentricity is 1 this is not a parabolic orbit.

Unlike standard orbits which are classified by their [[Astronomy:Orbital eccentricity|orbital eccentricity]], radial orbits are classified by their [[Astronomy:Specific orbital energy|specific orbital energy]], the constant sum of the total kinetic and potential energy, divided by the [[Physics:Reduced mass|reduced mass]]:
: <math> \epsilon = \frac{v^2}{2} - \frac{\mu}{x} </math>

where ''x'' is the distance between the centers of the masses, ''v'' is the relative velocity, and <math>\mu = {G}\left(m_1 + m_2\right)</math> is the [[Astronomy:Standard gravitational parameter|standard gravitational parameter]].

Another constant is given by:

: <math> w = \frac{1}{x} - \frac{v^2}{2\mu} = \frac{-\epsilon}{\mu} </math>

*For elliptic trajectories, w is positive. It is the inverse of the [[Astronomy:Apsis|apoapsis distance]] (maximum distance).
*For parabolic trajectories, w is zero.
*For hyperbolic trajectories, w is negative, It is <math> \textstyle \frac{-v_\infty^2}{2\mu} </math> where <math> \textstyle v_\infty </math> is the velocity at infinite distance.

== Time as a function of distance ==

Given the separation and velocity at any time, and the total mass, it is possible to determine the position at any other time.

The first step is to determine the constant w. Use the sign of w to determine the orbit type.
: <math> w = \frac{1}{x_0} - \frac{v_0^2}{2\mu} </math>
where <math> \textstyle x_0 </math> and <math> \textstyle v_0 </math> are the separation and relative velocity at any time.

=== Parabolic trajectory  ===
:<math> t(x) = \sqrt{ \frac{2x^3}{9\mu} }</math>

where ''t'' is the time from or until the time at which the two masses, if they were point masses, would coincide, and ''x'' is the separation.

This equation applies only to radial parabolic trajectories, for general parabolic trajectories see [[Astronomy:Parabolic trajectory#Barker's equation|Barker's equation]].

=== Elliptic trajectory ===

: <math> t(x, w) = \frac{ \arcsin\left( \sqrt{w\,x} \right) - \sqrt{w\,x\ (1 - w\,x)} }{ \sqrt{2\mu\, w^3} } </math>
where ''t'' is the time from or until the time at which the two masses, if they were point masses, would coincide, and ''x'' is the separation.

This is the radial Kepler equation.<ref>Brown, Kevin; [http://www.mathpages.com/rr/s4-03/4-03.htm ''MathPages'']</ref>

See also [[Physics:Equations for a falling body#Overview|equations for a falling body]].

=== Hyperbolic trajectory ===
:<math> t( x, w ) =  \frac{ \sqrt{ (|w|x)^2 + |w|x} - \ln\left(\sqrt{|w|x} + \sqrt{1 + |w|x}\right) }{ \sqrt{2\mu \, |w|^3} }</math>
where ''t'' is the time from or until the time at which the two masses, if they were point masses, would coincide, and ''x'' is the separation.

=== Universal form (any trajectory) ===

The radial Kepler equation can be made "universal" (applicable to all trajectories):

: <math> t(x, w) = \lim_{u \to w} \frac{ \arcsin\left(\sqrt{u\, x}\right) - \sqrt{u\,x\ (1 - u\,x)} }{ \sqrt{2\mu\, u^3} } </math>
or by expanding in a power series:

: <math> t( x, w ) = \frac{1}{\sqrt{2\mu}} \left. \left(\frac{2}{3} x^\frac{3}{2} + \frac{1}{5} w x^\frac{5}{2} + \frac{3}{28} w^2 x^\frac{7}{2} + \frac{5}{72} w^3 x^\frac{9}{2} + \frac{35}{704} w^4 x^\frac{11}{2} \cdots \right) \right|_{-1 < w \cdot x < 1}</math>

== The radial Kepler problem (distance as function of time) ==

The problem of finding the separation of two bodies at a given time, given their separation and velocity at another time, is known as the [[Physics:Kepler problem|Kepler problem]]. This section solves the Kepler problem for radial orbits.

The first step is to determine the constant <math> \textstyle w </math>. Use the sign of <math> \textstyle w </math> to determine the orbit type.
: <math> w = \frac{1}{x_0} - \frac{v_0^2}{2\mu} </math>
Where <math> \textstyle x_0 </math> and <math> \textstyle v_0 </math> are the separation and velocity at any time.

=== Parabolic trajectory ===

::<math> x(t) = \left( \frac{9}{2} \mu t^2 \right)^\frac{1}{3} </math>

See also position as function of time in a straight escape orbit.

=== Universal form (any trajectory) ===

Two intermediate quantities are used: w, and the separation at time t the bodies would have if they were on a parabolic trajectory, p.

: <math> w = \frac{1}{x_0} - \frac{v_0^2}{2\mu} \quad \text{and} \quad p = \left( \frac{9}{2} \mu t^2 \right)^\frac{1}{3} </math>

Where t is the time, <math> x_0 </math> is the initial position, <math> v_0 </math> is the initial velocity, and <math>\mu ={G}\left(m_1 + m_2\right)</math>.

The [[Astronomy:Kepler's equation#Inverse radial Kepler equation|inverse radial Kepler equation]] is the solution to the radial Kepler problem:
: <math>x(t) =
  \sum_{n=1}^\infty \left(
    \lim_{r \to 0} \left[
      \frac{w^{n-1} p^n}{n!}
      \frac{\mathrm{d}^{n-1}}{\mathrm{d}r^{n-1}} \left(
        r^n \left[
          \frac{3}{2} \left( \arcsin\left[\sqrt{r}\right] - \sqrt{r - r^2} \right)
        \right]^{-\frac{2}{3}n}
      \right)
    \right]
  \right)
</math>

Evaluating this yields:
:<math>x(t) =  p - \frac{1}{5} w p^2 - \frac{3}{175} w^2 p^3 - \frac{23}{7875} w^3 p^4 - \frac{1894}{3031875} w^4 p^5 - \frac{3293}{21896875} w^5 p^6 - \frac{2418092}{62077640625}w^6 p^7 \cdots</math>

Power series can be easily differentiated term by term. Repeated differentiation gives the formulas for the velocity, acceleration, jerk, snap, etc.
<!--
t(x)= \frac{ E - \sin(E) }{ \sqrt{ 8 \mu } \cdot w^\frac{3}{2}} | { E = 2 \arcsin \left( \sqrt{ \frac{x}{d} } \right) }
t(x)= \frac{\sinh(E) - E }{ \sqrt{ 8 \mu } \cdot w^\frac{3}{2}} | { E = 2 \arcsin( \sqrt{ x \cdot v_\inf  } )
t(x)= \frac{ E - \sin(E) }{ \sqrt{ 8 \mu } \cdot w^\frac{3}{2}} | { E = 2 \arcsin( \sqrt{ x \cdot w } ) }
-->

==Orbit inside a radial shaft==
The orbit inside a radial shaft in a uniform spherical body<ref>Strictly this is a contradiction. However, it is assumed that the shaft has a negligible influence on the gravity.</ref> would be a [[Physics:Simple harmonic motion|simple harmonic motion]], because gravity inside such a body is proportional to the distance to the center. If the small body enters and/or exits the large body at its surface the orbit changes from or to one of those discussed above. For example, if the shaft extends from surface to surface a closed orbit is possible consisting of parts of two cycles of simple harmonic motion and parts of two different (but symmetric) radial elliptic orbits.

==See also==
*[[Astronomy:Kepler's equation|Kepler's equation]]
*[[Physics:Kepler problem|Kepler problem]]
*[[Astronomy:List of orbits|List of orbits]]

== References ==
* Cowell, Peter (1993), Solving Kepler's Equation Over Three Centuries, William Bell.
{{Reflist}}

== External links ==
* Kepler's Equation at Mathworld [http://mathworld.wolfram.com/KeplersEquation.html]

{{DEFAULTSORT:Radial Trajectory}}
[[Category:Orbits]]
[[Category:Astrodynamics]]

{{Sourceattribution|Radial trajectory}}