** Interplanetary flight: **

*1) Introduction*

*2) **Heliocentric Transfer Orbit*

*4) Determining Orbital Elements*

*5) Hyperbolic Departure and Approach*

On New Year’s Day 1801, Giuseppe Piazzi of Italy observed for the first time an object that turned out to be the minor planet Ceres. He was only able to observe the asteroid for about one month before it was lost in the glare of the Sun. The challenge of rediscovering Ceres when it reappeared from behind the Sun seduced the intellect of young German mathematician Carl Friedrich Gauss. Exactly one year later, on New Year’s Day 1802, Ceres was rediscovered precisely where the ingenious and detailed calculations of Gauss had predicted she must be found.

The method that Gauss used is just as pertinent today as it was in 1802, but for a different reason. The data that Gauss used to determine the orbit of Ceres consisted of the right ascension and declination at three observation times. His method is much simplified if the original data consists to two position vectors and the time-of-flight between them. The technique of determining an orbit from two position vectors and time is of considerable interest to modern astrodynamics since it has direct application in the solution of intercept and rendezvous. We may define the Gauss problem as follows: Given position vectors **r _{1}** and

**r**, the time of flight from

_{2}**r**to

_{1}**r**, which we will call

_{2}*t*, and the direction of motion, find velocity vectors

**v**and

_{1}**v**. By “direction of motion” we mean whether the satellite is to go from

_{2}**r**to

_{1}**r**the “short way”, through an angular change () of less than radians, or the “long way”, through an angular change greater than radians. Obviously, there are an infinite number of orbits passing through

_{2}**r**and

_{1}**r**, but only two that have the specified time-of-flight – one for each possible direction of motion.

_{2}

One thing is immediately obvious from Figure 5.7; the two vectors **r _{1}** and

**r**uniquely define the plane of the transfer orbit. If the vectors

_{2}**r**and

_{1}**r**are collinear and in opposite directions (=), the plane of the transfer orbit is not determined and a unique solution for

_{2}**v**and

_{1}**v**is not possible. If the two position vectors are collinear and in the same direction (=0 or 2), the orbit is a degenerate conic, but a unique solution is possible of

_{2}**v**and

_{1}**v**. The relationship between the four vectors

_{2}**r**,

_{1}**r**,

_{2}**v**and

_{1}**v**is contained in the

_{2}*f*and

*g*expressions given below. It is not surprising, therefore, that nearly every known method for solving the Gauss problem may be derived from the

*f*and

*g*relations.

where,

As previously defined, is true anomaly, *E* is eccentric anomaly, *a* is semi-major axis, and *GM* is the standard gravitational parameter, while *p* is a geometrical constant of the conic section called the *parameter* or *semi-latus rectum*.

Consider equations (5.5), (5.6) and (5.7). There are seven variables — **r _{1}**,

**r**, , t, p, a and E — but the first four are known, so what we have is three equations in three unknowns. The only trouble is that the equations are transcendental in nature, so a trial-and-error solution is necessary. We may outline the general method of solution as follows:

_{2}- Guess a trial value for one of the three unknowns, p, a or E directly or indirectly by guessing some other parameter of the transfer orbit that in turn establishes p, a or E.
- Use equations (5.5) and (5.7) to compute the remaining two unknowns.
- Test the result by solving equation (5.6) for
*t*and check it against the given value of time-of-flight. - If the computed value of
*t*does not agree with the given value, adjust the trial value of the iteration variable and repeat the procedure until it does agree.

The last step is perhaps the most important of all, since the method used to adjust the trial value of the iteration variable is what determines how quickly the procedure converges to a solution. There are several methods for solving the Gauss problem, however the only one we will present here is the p-iteration technique. The method consists of guessing a trial value of *p*from which we can compute *a* and E. The first step in the solution is to find an expression for *a* as a function of *p* and the given information. We will find it convenient to define three constants that may be determined from the given information:

where r_{1} = |**r _{1}**| and r

_{2}= |

**r**|.

_{2}Once *p* is specified, a unique value of *a* is determined from

Once we have selected a trial value of *p* and computed *a* from equation (5.12), we are ready to solve for *t* and check it against the given time-of-flight. First, however, we need to determine E, or F in case *a* is negative (i.e. the orbit is hyperbolic).

From the trail value of *p* and the known information, we can compute f, g and from equations (5.5), (5.6) and (5.7). If *a* is positive, we can determine E from equations (5.5) and (5.7), rearranged as follows:

If *a* is negative, the corresponding f and g expressions involving F yield the following. Since we always assume F is positive, there is no ambiguity in determining F from this one equation.

The time-of-flight may now be determined from equation (5.16) or the corresponding equation involving F:

The limiting values of *p* correspond to two parabolic orbits passing through **r _{1}** and

**r**. The values of

_{2}*p*that specify the parabolic orbits we’ll call p

_{i}and p

_{ii}. For less than radians,

*p*must lie between p

_{i}and infinity. For greater than radians,

*p*must lie between 0 and p

_{ii}. Since it is important that the first trial value, as well as all subsequent guesses for

*p*, lie within the prescribed limits, we should first compute p

_{i}or p

_{ii}.

The method used to adjust the trial value of *p* to give the desired time-of-flight is crucial in determining how rapidly *p* converges to a solution. Several simple methods may be used successfully, such as the bisection or linear interpolation techniques.

In the *bisection* method we must find two trial values of *p*, one that gives too small a value for *t* and one that gives too large a value. The solution is then bracketed and, by choosing our next trial value half way between the first two, we can keep it bracketed while reducing the interval of uncertainty to some arbitrarily small value.

In the *linear interpolation* method we choose two trial values of *p*, which we will call p_{n-1} and p_{n}. If t_{n-1} and t_{n} are the times-of-flight corresponding to these trial values of *p*, then we select a new value from,

This scheme can be repeated, always retaining the latest two trial values of *p* and their corresponding times-of-flight for use in computing a still better trial value from equation (5.20). It is not necessary that the initial two trial values bracket the answer.

We can summarize the steps involved in solving the Gauss problem via the p-iteration technique as follows:

- Evaluate the constants k, and m from r
_{1}, r_{2}and using equations (5.9) through (5.11). - Determine the limits on the possible values of
*p*by evaluating p_{i}and p_{ii}from equations (5.18) and (5.19). - Pick a trial value of
*p*within the appropriate limits. - Using the trial value of
*p*, solve for*a*from equation (5.12). The type conic orbit will be known from the value of*a*. - Solve for f, g and from equations (5.5), (5.6) and (5.7).
- Solve for E or F, as appropriate, using equations (5.13) and (5.14) or equation (5.15).
- Solve for
*t*from equation (5.16) or (5.17) and compare it with the desired time-of-flight. - Adjust the trial value of
*p*using one of the iteration methods discussed above until the desired time-of-flight is obtained. - Evaluate from equation (5.8) and then solve for
**v**and_{1}**v**using equations (5.3) and (5.4)._{2}

The p-iteration method converges in all cases except when **r _{1}** and

**r**are collinear. Its main disadvantage is that separate equations are used for the ellipse and hyperbola.

_{2}**Selecting a Transfer Orbit**

Each time the Gauss problem is solved, the result gives just one of an infinite number of possible transfer orbits. It was previously stated that is generally desirable that the transfer orbit be tangential to Earth’s orbit at departure. This is true only in that it minimizes the V required to inject the spacecraft into its transfer orbit; however, it likely results in a less than optimum condition at target intercept. A one-tangent burn produces a trajectory that crosses the orbit of the target planet with a relatively large flight path angle, resulting in a large relative velocity between the spacecraft and planet. This relative velocity can be significantly reduced by selecting a transfer orbit that reduces the angle between the velocity vectors of the spacecraft and target at the moment of intercept. Improving the intercept condition (1) increases the duration of a close flyby encounter, (2) reduces the V required for orbit insertion, or (3) lowers the spacecraft’s velocity at atmospheric entry.

Tables 2 and 3 below provide sample data for a hypothetical mission to Mars in the year 2020. Table 2 gives the V required for Trans-Mars Injection (TMI) for a variety of different departure dates and times of flight. TMI is the maneuver that places the spacecraft into a trajectory that will intercept Mars at the desired place and time. In this sample, it is assumed that TMI is performed from an Earth parking orbit with an altitude of 200 km. Table 3 gives the V required for Mars-Orbit Insertion (MOI) for the same departure dates and times of flight found in Table 2. MOI, as it names implies, is the maneuver that slows the spacecraft to a velocity that places it into the desired orbit around Mars. In this sample, it is assumed that MOI is performed at periapsis of a insertion orbit with a periapsis altitude of 1,000 km and an apoapsis altitude of 33,000 km. Placing a spacecraft into a high eccentricity orbit such as this is common, as it provides for a MOI burn with a relatively low V.

Table 2 | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|

Trans-Mars Injection DV (m/s), launch altitude = 200 km | |||||||||||

Departure Date, 2020 |
Time of Flight (days) | ||||||||||

180 | 185 | 190 | 195 | 200 | 205 | 210 | 215 | 220 | 225 | 230 | |

7/7 | 3876 | 3862 | 3854 | 3851 | 3853 | 3863 | 3881 | 3912 | 3962 | 4043 | 4180 |

7/12 | 3841 | 3830 | 3824 | 3823 | 3826 | 3835 | 3851 | 3877 | 3917 | 3978 | 4074 |

7/19 | 3819 | 3812 | 3808 | 3808 | 3811 | 3819 | 3833 | 3853 | 3882 | 3925 | 3988 |

7/26 | 3834 | 3829 | 3826 | 3826 | 3829 | 3836 | 3846 | 3862 | 3883 | 3913 | 3956 |

8/2 | 3892 | 3887 | 3885 | 3884 | 3886 | 3890 | 3897 | 3908 | 3923 | 3943 | 3972 |

8/9 | 3999 | 3994 | 3990 | 3987 | 3987 | 3987 | 3991 | 3996 | 4005 | 4017 | 4034 |

8/16 | 4162 | 4154 | 4147 | 4141 | 4137 | 4133 | 4131 | 4131 | 4133 | 4138 | 4146 |

8/23 | 4386 | 4373 | 4362 | 4351 | 4341 | 4332 | 4325 | 4318 | 4313 | 4310 | 4309 |

Table 3 | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|

Mars Orbit Insertion DV (m/s), insertion orbit = 1000 × 33000 km | |||||||||||

Departure Date, 2020 |
Time of Flight (days) | ||||||||||

180 | 185 | 190 | 195 | 200 | 205 | 210 | 215 | 220 | 225 | 230 | |

7/7 | 1371 | 1258 | 1163 | 1086 | 1025 | 982 | 959 | 957 | 984 | 1052 | 1187 |

7/12 | 1290 | 1188 | 1102 | 1033 | 979 | 940 | 918 | 915 | 933 | 982 | 1074 |

7/19 | 1186 | 1097 | 1024 | 965 | 920 | 888 | 870 | 866 | 879 | 911 | 970 |

7/26 | 1093 | 1019 | 957 | 909 | 872 | 847 | 833 | 830 | 840 | 864 | 905 |

8/2 | 1016 | 954 | 904 | 865 | 837 | 818 | 808 | 808 | 817 | 836 | 867 |

8/9 | 957 | 907 | 868 | 838 | 817 | 804 | 799 | 801 | 811 | 828 | 853 |

8/16 | 920 | 881 | 852 | 830 | 816 | 809 | 808 | 813 | 823 | 839 | 862 |

8/23 | 910 | 881 | 860 | 846 | 838 | 836 | 839 | 846 | 857 | 873 | 893 |

- Calculate and
_{o},

( and _{o} are always less than 360^{o})

The angle , *longitude of periapsis*, is sometimes used in place of argument of periapsis. As a substitute for the time of periapsis passage, any of the following may be used to locate the spacecraft at a particular time, t_{o}, known as the “epoch”: _{o}, *true anomaly at epoch*, u_{o}, *argument of latitude at epoch*, or _{o}, *true longitude at epoch*.

If there is no periapsis (circular orbit), then is undefined, and _{o} = + u_{o}. If there is no ascending node (equatorial orbit), then both and u_{o} are undefined, and _{o} = + _{o}. If the orbit is both circular and equatorial, _{o} is simply the true angle from **x** to **r _{o}**, both of which are always defined.

The procedure outlined above describes a spacecraft in a solar orbit, but the method works equally well for satellites in Earth orbit, or around another planet or moon, where the position and velocity vectors are known in the geocentric-equatorial reference plane. Note, however, that it is customary for the geocentric-equatorial coordinate system to use unit vectors **i**, **j** and **k**instead of **x**, **y** and **z** as used in the heliocentric-ecliptic system.

PROBLEM 5.3A flight to Mars is launched on 2020-7-20, 0:00 UT. The planned time of flight is 207 days. Earth's postion vector at departure is 0.473265X- 0.899215YAU. Mars' postion vector at intercept is 0.066842X+ 1.561256Y+ 0.030948ZAU. Calculate the parameter and semi-major axis of the transfer orbit. SOLUTION, Given: t = 207 daysr= 0.473265_{1}X- 0.899215YAUr= 0.066842_{2}X+ 1.561256Y+ 0.030948ZAU GM = 1.327124×10^{20}m^{3}/s^{2}= 1.327124×10^{20}/ (149.597870×10^{9})^{3}= 3.964016×10^{-14}AU^{3}/s^{2}From vector magnitude,r= SQRT[ 0.473265_{1}^{2}+ (-0.899215)^{2}]r= 1.016153 AU_{1}r= SQRT[ 0.066842_{2}^{2}+ 1.561256^{2}+ 0.030948^{2}]r= 1.562993 AU From vector dot product, = arccos[ (0.473265 × 0.066842 - 0.899215 × 1.561256) / (1.016153 × 1.562993) ] = 149.770967_{2}^{o}Equations (5.9), (5.10) and (5.11), k = r_{1}× r_{2}× (1 - cos ) k = 1.016153 × 1.562993 × (1 - cos(149.770967)) k = 2.960511 AU = r_{1}+ r_{2}= 1.016153 + 1.562993 = 2.579146 AU m = r_{1}× r_{2}× (1 + cos ) m = 1.016153 × 1.562993 × (1 + cos(149.770967)) m = 0.215969 AU Equations (5.18) and (5.19), p_{i}= k / ( + SQRT(2 × m)) p_{i}= 2.960511 / (2.579146 + SQRT(2 × 0.215969)) p_{i}= 0.914764 AU p_{ii}= k / ( - SQRT(2 × m)) p_{ii}= 2.960511 / (2.579146 - SQRT(2 × 0.215969)) p_{ii}= 1.540388 AU Since < , 0.914764 < p < Equation (5.12), Select trial value, p = 1.2 AU a = m × k × p / [(2 × m -^{2}) × p^{2}+ 2 × k × × p - k^{2}] a = 0.215969 × 2.960511 × 1.2 / [(2 × 0.215969 - 2.579146^{2}) × 1.2^{2}+ 2 × 2.960511 × 2.579146 × 1.2 - 2.960511^{2}] a = 1.270478 AU Equations (5.5), (5.6) and (5.7), f = 1 - r_{2}/ p × (1 - cos ) f = 1 - 1.562993 / 1.2 × (1 - cos(149.770967)) f = -1.427875 g = r_{1}× r_{2}× sin / SQRT[ GM × p ] g = 1.016153 × 1.562993 × sin(149.770967) / SQRT[ 3.964016×10^{-14}× 1.2 ] g = 3,666,240 = SQRT[ GM / p ] × tan(/2) × [(1 - cos ) / p - 1/r_{1}- 1/r_{2}] = SQRT[ 3.964016×10^{-14}/ 1.2 ] × tan(149.770967/2) × [(1 - cos(149.770967)) / 1.2 - 1/1.016153 - 1/1.562993 ] = -4.747601×10^{-8}Equation (5.13), E = arccos[ 1 - r_{1}/ a × (1 - f) ] E = arccos[ 1 - 1.016153 / 1.270478 × (1 + 1.427875) ] E = 2.798925 radians Equation (5.16), t = g + SQRT[ a^{3}/ GM ] × (E - sin E) t = 3,666,240 + SQRT[ 1.270478^{3}/ 3.964016×10^{-14}] × (2.798925 - sin(2.798925)) t = 21,380,951 s = 247.4647 days Select new trial value of p and repeat above steps, p = 1.300000 AU, a = 1.443005 AU, t = 178.9588 days Equation (5.20), p_{n+1}= p_{n}+ (t - t_{n}) × (p_{n}- p_{n-1}) / (t_{n}- t_{n-1}) p_{n+1}= 1.3 + (207 - 178.9588) × (1.3 - 1.2) / (178.9588 - 247.4647) p_{n+1}= 1.259067 AU Recalculate using new value of p, p = 1.259067 AU, a = 1.336197 AU, t = 201.5624 days Perform additional iterations, p = 1.249221 AU, a = 1.318624 AU, t = 207.9408 days p = 1.250673 AU, a = 1.321039 AU, t = 206.9733 days p = 1.250633 AU, a = 1.320971 AU, t = 206.9999 days <-- close enough

PROBLEM 5.4For the Mars transfer orbit in Problem 5.3, calculate the departure and intecept velocity vectors. SOLUTION, Given:r= 0.473265_{1}X- 0.899215YAUr= 0.066842_{2}X+ 1.561256Y+ 0.030948ZAUr= 1.016153 AU_{1}r= 1.562993 AU p = 1.250633 AU a = 1.320971 AU = 149.770967_{2}^{o}Equations (5.5), (5.6) and (5.7), f = 1 - r_{2}/ p × (1 - cos ) f = 1 - 1.562993 / 1.250633 × (1 - cos(149.770967)) f = -1.329580 g = r_{1}× r_{2}× sin / SQRT[ GM × p ] g = 1.016153 × 1.562993 × sin(149.770967) / SQRT[ 3.964016×10^{-14}× 1.250633 ] g = 3,591,258 = SQRT[ GM / p ] × tan(/2) × [(1 - cos ) / p - 1/r_{1}- 1/r_{2}] = SQRT[ 3.964016×10^{-14}/ 1.250633 ] × tan(149.770967/2) × [(1 - cos(149.770967)) / 1.250633 - 1/1.016153 - 1/1.562993 ] = -8.795872×10^{-8}= 1 - r_{1}/ p × (1 - cos ) = 1 - 1.016153 / 1.250633 × (1 - cos(149.770967)) = -0.514536 Equation (5.3),v= (_{1}r- f ×_{2}r) / g_{1}v= [(0.066842 + 1.329580 × 0.473265) / 3,591,258]_{1}X+ [(1.561256 + 1.329580 × -0.899215) / 3,591,258]Y+ [(0.030948 + 1.329580 × 0) / 3,591,258]Zv= 0.000000193828_{1}X+ 0.000000101824Y+ 0.00000000861759ZAU/s × 149.597870×10^{9}v= 28996.2_{1}X+ 15232.7Y+ 1289.2Zm/s Equation (5.4),v= ×_{2}r+ ×_{1}v_{1}v= [-8.795872×10_{2}^{-8}× 0.473265 - 0.514536 × 0.000000193828]X+ [-8.795872×10^{-8}× -0.899215 - 0.514536 × 0.000000101824]Y+ [-8.795872×10^{-8}× 0 - 0.514536 × 0.00000000861759]Zv= -0.000000141359_{2}X+ 0.0000000267017Y- 0.00000000443406ZAU/s × 149.597870×10^{9}v= -21147.0_{2}X+ 3994.5Y- 663.3Zm/s

PROBLEM 5.5 For the Mars transfer orbit in Problems 5.3 and 5.4, calculate the orbital elements. SOLUTION, Problem can be solved using eitherr&_{1}vor_{1}r&_{2}v– we will use_{2}r&_{1}v. Given:_{1}r= (0.473265_{1}X- 0.899215YAU) × 149.597870×10^{9}m/AU = 7.079944×10^{10}X- 1.345206×10^{11}Ymr= 1.016153 × 149.597870×10_{1}^{9}= 1.520144×10^{11}m GM = 1.327124×10^{20}m^{3}/s^{2}From problem 5.4,v= 28996.2_{1}X+ 15232.7Y+ 1289.2Zm/s Also,v= SQRT[ 28996.2^{2}+ 15232.7^{2}+ 1289.2^{2}] = 32,779.2 m/s Equations (5.21) and (5.22),h= (r_{Y}v_{Z}- r_{Z}v_{Y})X+ (r_{Z}v_{X}- r_{X}v_{Z})Y+ (r_{X}v_{Y}- r_{Y}v_{X})Zh= (-1.345206×10^{11}× 1289.2 - 0 × 15232.7)X+ (0 × 28996.2 - 7.079944×10^{10}× 1289.2)Y+ (7.079944×10^{10}× 15232.7 + 1.345206×10^{11}× 28996.2)Zh= -1.73424×10^{14}X- 9.12746×10^{13}Y+ 4.97905×10^{15}Zn= -h_{Y}X+ h_{X}Yn= 9.12746×10^{13}X- 1.73424×10^{14}YAlso,h= SQRT[ (-1.73424×10^{14})^{2}+ (-9.12746×10^{13})^{2}+ (4.97905×10^{15})^{2}] = 4.98291×10^{15}n= SQRT[ (9.12746×10^{13})^{2}+ (1.73424×10^{14})^{2}] = 1.95977×10^{14}Equation (5.23),e= [(v^{2}- GM /r) ×r- (r • v) ×v] / GMv^{2}- GM /r= 32779.2^{2}- 1.327124×10^{20}/ 1.520144×10^{11}= 2.01451×10^{8}r • v= 7.079944×10^{10}× 28996.2 - 1.345206×10^{11}× 15232.7 + 0 x 1289.2 = 3.80278×10^{12}e= [2.01451×10^{8}× (7.079944×10^{10}X- 1.345206×10^{11}Y) - 3.80278×10^{12}× (28996.2X+ 15232.7Y+ 1289.2Z) ] / 1.327124×10^{20}e= 0.106639X- 0.204632Y- 0.000037ZEquations (5.24) and (5.25), a = 1 / ( 2 /r-v^{2}/ GM ) a = 1 / ( 2 / 1.520144×10^{11}- 32779.2^{2}/ 1.327124×10^{20}) a = 1.97614×10^{11}me= SQRT[ 0.106639^{2}+ (-0.204632)^{2}+ (-0.000037)^{2}]e= 0.230751 Equations (5.26) though (5.30), cos i =h/_{Z}hcos i = 4.97905×10^{15}/ 4.98291×10^{15}i = 2.255^{o}cos =n/_{X}ncos = 9.12746×10^{13}/ 1.95977×10^{14}= 297.76^{o}cos =n • e/ (n×e) cos = (9.12746×10^{13}× 0.106639 - 1.73424×10^{14}× (-0.204632) + 0 × (-0.000037)) / (1.95977×10^{14}× 0.230751) = 359.77^{o}cos_{o}=e • r/ (e×r) cos_{o}= (0.106639 × 7.079944×10^{10}- 0.204632 × (-1.345206×10^{11}) - 0.000037 × 0) / (0.230751 × 1.520144×10^{11})_{o}= 0.226^{o}cos u_{o}=n • r/ (n×r) u_{o}= 0 (launch point = ascending node) Equations (5.31) and (5.32), = + = 297.76 + 359.77 = 297.53^{o}_{o}= + +_{o}_{o}= 297.76 + 359.77 + 0.23_{o}= 297.76^{o}