6 - Asteroid Orbits and Encounters

In order to estimate the likely nature of an encounter between the probe and a NEO a number of close encounters with known NEOs were modelled to provide a set of sample missions. A list of all currently-predicted NEO close encounters for the next several decades was obtained from the IAU Minor Planet Center web site; from this every encounter within 0.07 AU (slightly larger than the planned maximum encounter range) for the next 10 years was extracted (see Table 5.1). Each encounter was then modelled in detail in order to determine its geometry.

Encounter Modelling. For each encounter it was assumed that the probe would intercept the target at the point of closest approach whilst travelling radially outwards from Earth at a notional velocity of 5 km/s (see Section 7). To locate the intercept point and determine the encounter relative velocity and angle the asteroid's position and velocity vectors at the time of closest approach were then needed. This was found by making use of the Jet Propulsion Laboratory's 'Horizons' on-line ephemeris generation system. 'Horizons' allows users to generate an ephemeris covering a large number of possible parameters for any body or pair of bodies in the solar system. For example, it can be used to generate a list of position and velocity vectors of an asteroid relative to Earth on an hourly basis around the time of predicted closest encounter. The output can be referred to one of several co-ordinate systems (ecliptic plane is most useful for this calculation) and presented in various sets of units. An example output run from 'Horizons' is shown at Appendix 1 to this section.

Using the range output the time of closest approach was confirmed and the position and velocity vector for that time noted. The position vector relative to Earth then gave the direction (and hence velocity vector, normalized to 5 km/s) of the probe. With position and velocity vectors for both asteroid and probe the circumstances of the encounter could be modelled.

Fig 6.1 shows the encounter geometry, projected onto the plane of the ecliptic. The vector difference between the asteroid and probe velocity vectors gives the relative encounter speed. Also of interest though are the two angles illustrated in more detail in Fig 6.2.

Figure 6.1. Asteroid Encounter Geometry

Figure 6.2. Asteroid Encounter Geometry - Detail

The Probe/Sun Angle is the angle between the probe flight path and the Sun. This is important for probe configurations where the high-gain antenna is fixed, as it determines the Sun's aspect relative to the probe body assuming that the probe cruises with its antenna in Earth-pointing mode. The Asteroid/Sun Angle (also known as the Phase Angle indicates the illumination of the asteroid as the probe approaches it. A low value indicates that the asteroid approaches almost fully illuminated, 90° that it is half-full during approach and values approaching 180° that the asteroid is coming out of the Sun relative to the probe. Although the probe flightpath can usually be adjusted to ensure that the encounter is over the Sunward side of the asteroid a high Asteroid/Sun angle should be avoided as it makes acquisition of the target for tracking purposes difficult. Once the asteroid's position and velocity vectors are known then basic vector algebra can be used to find these angles via the probe's velocity vector and the Sun's direction vector from the Earth (calculated as a function of date).

The flyby parameters for the encounters listed at Table 5.1 are given at Table 6.1.

Asteroid
JD
Date
Rel V (km/s)
Sun/Flt Angle (°)
Rel V /Sun Angle (°)
1988 EG2450873.417
28-Feb-98
12.036
99.276
145.735
1996 FG32451143.25
25-Nov-98
11.437
96.039
147.577
1992 SK2451263.75
26-Mar-99
12.606
44.947
73.155
Golevka2451332.292
02-Jun-99
10.923
127.937
22.698
Mithra2451770.833
14-Aug-00
18.279
85.547
165.271
Nereus2452297.042
22-Jan-02
7.365
18.101
144.744
1994 PM2452868.083
16-Aug-03
26.891
98.303
18.946
1996 GT2452957.667
14-Nov-03
30.61
80.19
19.833
Minos2453038.375
02-Feb-04
12.7
83.118
31.2
Toutatis2453278.083
29-Sep-04
12.036
59.737
62.172
1992 BF2453433.208
03-Mar-05
10.071
106.131
26.048
1992 UY42453590.875
08-Aug-05
11.632
132.832
17.303
1991 VK2454122
21-Jan-07
10.009
15.576
135.319
Hathor2454395.75
22-Oct-07
15.869
66.992
166.103

Table 6.1. Flyby parameters for forthcoming asteroid encounters

From Table 6.1 it is apparent that not all encounters are suitable for flyby missions. The encounters with 1994 PM and 1996 GT both feature very high relative velocities, whilst six encounters (1996 EG, 1996 FG3, Mithra, Nereus, 1991VK and Hathor) all involve poor asteroid/Sun angles. Note that the encounter with Nereus, whilst very favourable in relative velocity terms, involves the target being poorly illuminated. The encounters with 1992 SK, Golevka, Minos, Toutatis, 1992 BF and 1992 UY4 are, however, all favourable. Thus 6 out of 14 encounters over the next 10 years are potentially suitable for probe flyby missions. Assuming that this is a representative sample, some 40% of all close NEO encounters should be suitable candidates for such flyby missions.

Fine-Tuning Encounters. The above encounters were all modelled by assuming that the probe was sent to intercept the asteroid at its point of closest approach to Earth. If that approach is significantly closer than the 0.05 AU baseline range of the probe then there is scope for bringing forward or delaying the flyby to make its parameters more favourable. Delaying the flyby in particular ensures that the probe is 'chasing' the asteroid and so reduces the encounter velocity, although the need to reach a more distant intercept point may mean that the probe actually has to be launched earlier than for a minimum-distance encounter, despite the later flyby date. Nonetheless such a move can have beneficial effects; for example, calculations show that intercepting Toutatis at a range of 0.05 AU on 7 Oct 2004, rather than at 0.01 AU on 29 Sep 2004, has the following effects:

In the event of a newly-discovered NEO being a suitable mission target, trajectory analysis will thus be required to identify the optimum launch and encounter dates, subject of course to constraints on the former.


Previous Contents Next

MSc Report Home Page

Appendix 1

Example output from JPL 'Horizons' Ephemeris Generator

(data truncated to fit page)

Date: Sat, 2 Aug 1997 04:06:49 -0700 (PDT)

From: Horizons On-Line Ephemeris System <horizons@ssd.jpl.nasa.gov>

Subject: SMALL BODY #C399_T(6037 (1988 EG)) (1/1)

Apparently-To: sjbradshaw@cix.compulink.co.uk

Automated mail xmit by PORT_LOGIN, PID = 28792 Sat Aug 2 04:06:48 1997

++++++++++++++++++++++++++++++++ (part 1 of 1) +++++++++++++++++++++++++++++++

$$SOH

*******************************************************************************

JPL On-Line Ephemerides Horizons

*******************************************************************************

TARGET BODY : 6037 (1988 EG) {Source : MPC23772}

COORD. CENTER : Earth (399) {Source : DE-0403LE-0403}

*******************************************************************************

EXECUTION DATE : Sat Aug 2 04:06:13 1997 (Pasadena time)

REQUESTED BY : PORT_LOGIN

REFERENCE FRAME: ICRF/J2000.0 rotated to ecliptic of Epoch J2000.0

START TIME : 1998 FEB 28 18:00 TDB

STOP TIME : 1998 MAR 01 06:00 TDB

STEP (MINUTES) : 60

OUTPUT UNITS : KM-S

OUTPUT FORMAT : 03

OUTPUT TYPE : GEOMETRIC cartesian states

*******************************************************************************

Initial heliocentric osc. elements wrt ecliptic and mean equinox of J2000.0:

EPOCH= 2450600.5 != 1997-Jun-01.0000000 (TDB)

EC= .499176323 QR= .635620962 TP= 2450817.8183546

OM= 182.947887 W= 241.5103969 IN= 3.4793031

Asteroid physical parameters:

GM= n.a. RAD= n.a. ROTPER= n.a.

H= 18.7 G= .150 B-V= n.a.

*******************************************************************************

JDTDB, , X, Y, Z, Vx, Vy, Vz, OWLT, R, dR/dt,

*******************************************************************************

$$EOH

$$SOE

2450873.250000000, 1998 FEB 28 18:00:00.000, 0.1153899407528818E+07…

2450873.291666667, 1998 FEB 28 19:00:00.000, 0.1103102400232166E+07…

2450873.333333333, 1998 FEB 28 20:00:00.000, 0.1052304271454364E+07…

2450873.375000000, 1998 FEB 28 21:00:00.000, 0.1001504974935174E+07…

2450873.416666667, 1998 FEB 28 22:00:00.000, 0.9507044615356922E+06…

2450873.458333333, 1998 FEB 28 23:00:00.000, 0.8999026822274327E+06…

2450873.500000000, 1998 MAR 01 00:00:00.000, 0.8490995910672247E+06…

2450873.541666667, 1998 MAR 01 01:00:00.000, 0.7982951392307281E+06…

2450873.583333333, 1998 MAR 01 02:00:00.000, 0.7474892779943049E+06…

2450873.625000000, 1998 MAR 01 03:00:00.000, 0.6966819617020786E+06…

2450873.666666667, 1998 MAR 01 04:00:00.000, 0.6458731418143213E+06…

2450873.708333333, 1998 MAR 01 05:00:00.000, 0.5950627698817551E+06…

2450873.750000000, 1998 MAR 01 06:00:00.000, 0.5442508005053997E+06…

$$EOE

$$SOD

*******************************************************************************

Symbol meaning

JDTDB Epoch Julian Date, Barycentric Dynamical Time

X x-component of position vector (km)

Y y-component of position vector (km)

Z z-component of position vector (km)

Vx x-component of velocity vector (km/sec)

Vy y-component of velocity vector (km/sec)

Vz z-component of velocity vector (km/sec)

OWLT One-way Newtonian light-time (sec)

R Range; distance from coordinate center (km)

dR/dt Range-rate; radial velocity wrt coord. center (km/sec)

Geometric states/elements have no aberration corrections applied.

Computations by ...

Solar System Dynamics Group, Horizons On-Line Ephemeris System

4800 Oak Grove Drive, Jet Propulsion Laboratory

Pasadena, CA 91109 USA

information: http://ssd.jpl.nasa.gov/

connect : telnet ssd.jpl.nasa.gov 6775

*******************************************************************************

$$EOD

$$EOF