Astronomy & Celestial Mechanics
the sky, computed
Generator

The tisserand-map generator

The contours a planetary encounter cannot cross
The contours a planetary encounter cannot cross. Contours of the Tisserand parameter with respect to Jupiter in the plane of semi-major axis and eccentricity, drawn for a coplanar orbit. An encounter with Jupiter moves a comet along one of these curves and never across one, because T is what the encounter conserves. The heavy contours are at T = 3 and T = 2, and they are the boundaries the comet families are defined by: T > 3 means no encounter is possible at all, since v∞²/v_J² = 3 − T and a negative squared speed is not a trajectory; 2 < T < 3 is the Jupiter family; below 2 the approach speed exceeds Jupiter's own orbital speed and the orbits are the nearly isotropic ones. The shaded boundary is the crossing condition — an orbit whose pericentre is outside Jupiter's, or whose apocentre is inside it, never meets the planet whatever its T. The five comets are placed at their JPL elements and labelled with the T the literature quotes; all five with inclination sit off the coplanar contours by exactly the cos i in the definition, which is why 1P/Halley's is negative — a retrograde orbit meets Jupiter at nearly twice Jupiter's speed.

Contours of the Tisserand parameter with respect to Jupiter in the plane of semi-major axis and eccentricity, drawn for a coplanar orbit. An encounter with Jupiter moves a comet *along* one of these curves and never across one, because T is what the encounter conserves. The heavy contours are at T = 3 and T = 2, and they are the boundaries the comet families are defined by: *T > 3 means no encounter is possible at all*, since v∞²/v_J² = 3 − T and a negative squared speed is not a trajectory; 2 < T < 3 is the Jupiter family; below 2 the approach speed exceeds Jupiter's own orbital speed and the orbits are the nearly isotropic ones. The shaded boundary is the crossing condition — an orbit whose pericentre is outside Jupiter's, or whose apocentre is inside it, never meets the planet whatever its T. The five comets are placed at their JPL elements and labelled with the T the literature quotes; all five with inclination sit off the coplanar contours by exactly the cos i in the definition, which is why 1P/Halley's is negative — a retrograde orbit meets Jupiter at nearly twice Jupiter's speed.

5 essays call tisserand-map. The drawing above is what it returns with no arguments at all; every call below passes it something, because a placement that passes nothing draws whichever member of the family the generator happens to default to rather than the one its essay argues about.

Where it is called

Every figure listed here is the same construction drawn at different numbers, so a correction to one is a correction to all of them.

After one revolution the error is 3π times longer than it is wide, and after 300 it is 2827. The two semi-axes of a fitted orbit's position uncertainty, against elapsed revolutions, for a solution whose semi-major axis is uncertain by 12 kilometres. The radial extent does not grow at all: a body on a slightly larger orbit is slightly further out and stays so. The along-track extent grows linearly, because δn/n = −(3/2)δa/a makes a semi-major-axis error into a mean-motion error and a mean-motion error into a phase that runs away — a·δM = 3πN·δa after N revolutions. The ratio is 3π ≈ 9.42 after a single revolution and 2827 after 300, which is why an asteroid recovered after one apparition is found within a few arcseconds of its predicted place along its own track and could be a long way from it in time. Every consequence of this in practice — that an impact probability is a one-dimensional integral rather than a volume, that a keyhole is an interval, that the next observation worth taking is the one across the track rather than the one that fits best — is a restatement of these two lines diverging. Orbits

An error that is nearly all in one direction

A fitted orbit's uncertainty is not a ball. Within a few revolutions it has collapsed onto a line along the track, because an error in the size of an orbit is an error in its period and an error in period is a phase that runs away — which is why an impact probability is an integral along a curve rather than over a volume.

 The contours a planetary encounter cannot cross. Contours of the Tisserand parameter with respect to Jupiter in the plane of semi-major axis and eccentricity, drawn for a coplanar orbit. An encounter with Jupiter moves a comet along one of these curves and never across one, because T is what the encounter conserves. The heavy contours are at T = 3 and T = 2, and they are the boundaries the comet families are defined by: T > 3 means no encounter is possible at all, since v∞²/v_J² = 3 − T and a negative squared speed is not a trajectory; 2 < T < 3 is the Jupiter family; below 2 the approach speed exceeds Jupiter's own orbital speed and the orbits are the nearly isotropic ones. The shaded boundary is the crossing condition — an orbit whose pericentre is outside Jupiter's, or whose apocentre is inside it, never meets the planet whatever its T. The five comets are placed at their JPL elements and labelled with the T the literature quotes; all five with inclination sit off the coplanar contours by exactly the cos i in the definition, which is why 1P/Halley's is negative — a retrograde orbit meets Jupiter at nearly twice Jupiter's speed. Gravitation

The number that survives the encounter

A comet that passes Jupiter comes away with every orbital element changed. One combination of them is not changed, and it is enough to recognise the comet afterwards, to sort the comet families, and to say where a spacecraft can and cannot go.

 Along a contour is free; across one costs, and 6 resonances sit on this one. Perihelion against aphelion, in units of Jupiter's orbit, with contours of constant v∞ — which is the Tisserand parameter through v∞² = (3 − T)v_p². Every contour crosses the line r_p = r_a = 1, because a spacecraft has to be at the planet's orbit to have an encounter there at all, and every point on one contour is reachable from every other point on it with no propellant: a flyby rotates the v∞ vector without changing its length, which moves the pump angle and slides the spacecraft along the curve. A burn is the only thing that moves it between curves, and that is the whole economy of a gravity-assist tour. The diagonals are resonant orbits — 1:1, 3:2, 2:1, 5:2, 3:1, 4:1 with the planet — and each is a straight line of slope −1 because a resonance fixes the semi-major axis and r_p + r_a = 2a. They matter because a spacecraft on one comes back to the same place at the same time as the planet, which is what makes a second encounter possible without waiting for a chance alignment. The marked points are where the v∞ = 0.3 contour meets each: a tour is a walk along the highlighted curve from one to the next, and the arithmetic that decides whether it can be walked is the turn one flyby delivers. At 1.35 Jupiter radii and 3.9 km s⁻¹ that turn is 163°, so the longest step drawn here needs 0.1 encounters — which is why a real tour has dozens of flybys and why the ones with the largest steps are the ones that need a deep-space manoeuvre in between. Spaceflight

The same planet, three times

A flyby cannot change the encounter speed, only its direction, so a tour has to be designed in the space of what is conserved. One pass moves a spacecraft along a single curve and no further than the planet can bend it — and reaching a distant target means walking that curve, returning to the same planet again and again.

 Along a contour is free; across one costs, and 6 resonances sit on this one. Perihelion against aphelion, in units of Jupiter's orbit, with contours of constant v∞ — which is the Tisserand parameter through v∞² = (3 − T)v_p². Every contour crosses the line r_p = r_a = 1, because a spacecraft has to be at the planet's orbit to have an encounter there at all, and every point on one contour is reachable from every other point on it with no propellant: a flyby rotates the v∞ vector without changing its length, which moves the pump angle and slides the spacecraft along the curve. A burn is the only thing that moves it between curves, and that is the whole economy of a gravity-assist tour. The diagonals are resonant orbits — 1:1, 3:2, 2:1, 5:2, 3:1, 4:1 with the planet — and each is a straight line of slope −1 because a resonance fixes the semi-major axis and r_p + r_a = 2a. They matter because a spacecraft on one comes back to the same place at the same time as the planet, which is what makes a second encounter possible without waiting for a chance alignment. The marked points are where the v∞ = 0.3 contour meets each: a tour is a walk along the highlighted curve from one to the next, and the arithmetic that decides whether it can be walked is the turn one flyby delivers. At 1.35 Jupiter radii and 3.9 km s⁻¹ that turn is 163°, so the longest step drawn here needs 0.1 encounters — which is why a real tour has dozens of flybys and why the ones with the largest steps are the ones that need a deep-space manoeuvre in between. Spaceflight

A map of the transfers that are free

Drawn as contours of perihelion against aphelion, the invariant that survives an encounter becomes a map. A flyby slides a spacecraft along its own contour and costs nothing; a burn is the only thing that moves it between contours — so tour design is reading a graph.

 An invariant that moves by 8.0e-3 once the planet's orbit is real. The Tisserand parameter of a comet on an orbit of semi-major axis 5 and eccentricity 0.8, followed through a close passage of Jupiter, integrated twice. The lower trace has Jupiter on a perfect circle, which is the problem the parameter is an exact constant of: it survives the encounter having moved by 1.0e-4, which is the integrator's own error and not a physical change, and the spike at the moment of closest approach is the osculating elements being briefly meaningless while the comet is inside Jupiter's sphere of influence rather than the constant failing. The upper trace is the same encounter with Jupiter on its real orbit, eccentricity 0.0489. The parameter comes out changed by 8.0e-3, 79 times as much, because the Jacobi constant exists only when the rotating frame is uniformly rotating and a planet on an ellipse does not provide one. That number is small and it is not negligible: comet families are separated by boundaries in this parameter placed to two decimal places, and a comet that drifts across one over several encounters has changed class without anything having happened to it that a single encounter could account for. Orbits

An invariant that is only almost one

The Tisserand parameter survives a close encounter with Jupiter exactly, and comet families are separated by boundaries in it drawn to two decimal places. The exactness holds for a Jupiter on a circle. Jupiter's eccentricity is 0.0489, and integrating the same encounter twice shows what that costs.

The whole library · All essays