Radial trajectory explained

In astrodynamics and celestial mechanics a radial trajectory is a Kepler orbit with zero angular momentum. Two objects in a radial trajectory move directly towards or away from each other in a straight line.

Classification

There are three types of radial trajectories (orbits).^[1]

• 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 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 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 orbital eccentricity, radial orbits are classified by their specific orbital energy, the constant sum of the total kinetic and potential energy, divided by the reduced mass:$$\backslash varepsilon\; =\; \backslash frac\; -\; \backslash frac$$where is the distance between the centers of the masses, is the relative velocity, and

\mu = G \left(m _{1} + m _{2\right)} is the standard gravitational parameter.

Another constant is given by:$$w\; =\; \backslash frac\; -\; \backslash frac\; =\; \backslash frac$$

• For elliptic trajectories, w is positive. It is the inverse of the apoapsis distance (maximum distance).
• For parabolic trajectories, w is zero.
• For hyperbolic trajectories, w is negative, It is

style (- v ^{2} _{inf ty})/(2 \mu) where style v _{infty} 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 . Use the sign of to determine the orbit type.$$w\; =\; \backslash frac\; -\; \backslash frac$$where $x\_0$ and $v\_0$ are the separation and relative velocity at any time.

Parabolic trajectory

$$t(x)\; =\; \backslash sqrt$$where is the time from or until the time at which the two masses, if they were point masses, would coincide, and is the separation.

This equation applies only to radial parabolic trajectories, for general parabolic trajectories see Barker's equation.

Elliptic trajectory

$$t(x,\; w)\; =\; \backslash frac$$where is the time from or until the time at which the two masses, if they were point masses, would coincide, and is the separation.

This is the radial Kepler equation.^[2]

Hyperbolic trajectory

$$t(x,\; w)\; =\; \backslash frac$$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):$$t(x,\; w)\; =\; \backslash lim\_\; \backslash frac$$or by expanding in a power series:$$t(x,\; w)\; =\; \backslash frac\; \backslash left.\; \backslash left(\backslash frac\; x^\backslash frac\; +\; \backslash frac\; w\; x^\backslash frac\; +\; \backslash frac\; w^2\; x^\backslash frac\; +\; \backslash frac\; w^3\; x^\backslash frac\; +\; \backslash frac\; w^4\; x^\backslash frac\; \backslash cdots\; \backslash right)\; \backslash right|\_$$

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 Kepler problem. This section solves the Kepler problem for radial orbits.

The first step is to determine the constant $w$. Use the sign of $w$ to determine the orbit type.$$w\; =\; \backslash frac\; -\; \backslash frac$$Where $x\_0$ and $v\_0$ are the separation and velocity at any time.

Parabolic trajectory

$$x(t)\; =\; \backslash left(\backslash frac\; \backslash mu\; t^2\; \backslash right)^\backslash frac$$

Universal form (any trajectory)

Two intermediate quantities are used:, and the separation at time the bodies would have if they were on a parabolic trajectory, .$$w\; =\; \backslash frac\; -\; \backslash frac\; \backslash quad\; \backslash text\; \backslash quad\; p\; =\; \backslash left(\backslash frac\; \backslash mu\; t^2\; \backslash right)^$$

Where is the time,

x _{0} is the initial position, v _{0} is the initial velocity, and \mu = G \left(m _{1} + m _{2\right)} .

The inverse radial Kepler equation is the solution to the radial Kepler problem:$$x(t)\; =\; \backslash sum\_^\backslash infty\; \backslash left(\backslash lim\_\; \backslash left[\backslash frac\{w^\{n-1\}\; p^n\}\{n!\}\; \backslash frac\{\backslash mathrm\{d\}^\{n-1\}\}\{\backslash mathrm\{d\}r^\{n-1\}\}\; \backslash left(r^n\; \backslash left[\; \backslash frac\{3\}\{2\}\; \backslash left(\backslash arcsin\backslash left[\backslash sqrt\{r\}\backslash right]\; -\; \backslash sqrt\; \backslash right)\; \backslash right]^\; \backslash right)\; \backslash right]\; \backslash right)$$

Evaluating this yields:$$x(t)\; =\; p\; -\; \backslash frac\; w\; p^2\; -\; \backslash frac\; w^2\; p^3\; -\; \backslash frac\; w^3\; p^4\; -\; \backslash frac\; w^4\; p^5\; -\; \backslash frac\; w^5\; p^6\; -\; \backslash fracw^6\; p^7\; \backslash cdots$$

Power series can be easily differentiated term by term. Repeated differentiation gives the formulas for the velocity, acceleration, jerk, snap, etc.

Orbit inside a radial shaft

The orbit inside a radial shaft in a uniform spherical body^[3] would be a 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

• Kepler's equation
• Kepler problem
• List of orbits

References

• Cowell, Peter (1993), Solving Kepler's Equation Over Three Centuries, William Bell.

External links

• Kepler's Equation at Mathworld http://mathworld.wolfram.com/KeplersEquation.html

Notes and References

1. Thomson, William Tyrrell; Introduction to Space Dynamics, Dover, 1986
2. Brown, Kevin; MathPages
3. Strictly this is a contradiction. However, it is assumed that the shaft has a negligible influence on the gravity.