Cylindrical isomorphic mapping applied to invariant manifold dynamics for Earth–Moon Missions

Several families of periodic orbits exist in the context of the circular restricted three-body problem. This work studies orbital motion of a spacecraft among these periodic orbits in the Earth–Moon system, using the planar circular restricted three-body problem model. A new cylindrical representation of the spacecraft phase space (i.e., position and velocity) is described, and allows representing periodic orbits and the related invariant manifolds. In the proximity of the libration points, the manifolds form a four-fold surface, if the cylindrical coordinates are employed. Orbits departing from the Earth and transiting toward the Moon correspond to the trajectories located inside this four-fold surface. The isomorphic mapping under consideration is also useful for describing the topology of the invariant manifolds, which exhibit a complex geometrical stretch-and-folding behavior as the associated trajectories reach increasing distances from the libration orbit. Moreover, the cylindrical representation reveals extremely useful for detecting periodic orbits around the primaries and the libration points, as well as the possible existence of heteroclinic connections. These are asymptotic trajectories that are ideally traveled at zero-propellant cost. This circumstance implies the possibility of performing concretely a variety of complex Earth–Moon missions, by combining different types of trajectory arcs belonging to the manifolds. This work studies also the possible application of manifold dynamics to defining a suitable, convenient end-of-life strategy for spacecraft placed in any of the unstable orbits. The final disposal orbit is an externally confined trajectory, never approaching the Earth or the Moon, and can be entered by means of a single velocity impulse (of modest magnitude) along the right unstable manifold that emanates from the Lyapunov orbit at \(L_2\) .     

Satellite resonance with longitude-dependent gravity—II Effects involving the eccentricity

For a satellite in a nominally circular orbit at arbitrary inclination whose mean motion is commensurable with the Earth's rotation, the dependence of gravity on longitude leads to a resonant variation in eccentricity as well as the long-period oscillation in longitude. Provided forces capable of processing perigee are present, it is shown that the change in eccentricity for a satellite captured in librational resonance is not secular but periodic.

There are corresponding resonance effects for a satellite in a nominally equatorial but eccentric orbit. Here the commensurability condition is that the longitudes of the apses shall be nearly repetitive relative to the rotating Earth. There will be a long-period oscillation in longitude which can take the form of either a libration (trapped) or a circulation (free), and there will also be an oscillation of the orbital plane having the same period as the precession of perigee relative to inertial space.     

Analysis of the orbit of 1970-97b (cosmos 378 rocket)

Cosmos 378 rocket, 1970-97B, entered orbit on 17 November 1970, with orbital inclination 74.0°, period 105 min and perigee height 230 km, and decayed on 30 September 1972 after 683 days in orbit. The RAE computer program PROP was used, with more than 1900 observations from 64 stations, to determine the orbit at 39 epochs between February 1971 and September 1972.The main aim of the analysis was to determine the atmospheric rotation rate from the decrease in orbital inclination, which was determined with a mean standard deviation of 0.0010° and a best standard deviation of 0.0003°. After removal of relevant perturbations, analysis of the variation in inclination between July 1971 and April 1972 yields the surprisingly low average atmospheric rotation rate of 0.75 ± 0.05 rev/day, at a mean height of 250 km. The local time at perigee is however strongly biassed towards daytime values (07–16 hr), so the results lend support to the picture of east-to-west winds by day and west-to-east winds by night.Values of scale height are obtained by analysis of the change in perigee height.     

Sun-synchronous orbits near critical inclination

The long period dynamics of Sun-synchronous orbits near the critical inclination 116.6° are investigated. It is known that, at the critical inclination, the average perigee location is unchanged by Earth oblateness. For certain values of semimajor axis and eccentricity, orbit plane precession caused by Earth oblateness is synchronous with the mean orbital motion of the apparent Sun (a Sun-synchronism). Sun-synchronous orbits have been used extensively in meteorological and remote sensing satellite missions. Gravitational perturbations arising from an aspherical Earth, the Moon, and the Sun cause long period fluctuations in the mean argument of perigee, eccentricity, inclination, and ascending node. Double resonance occurs because slow oscillations in the perigee and Sun-referenced ascending node are coupled through the solar gravity gradient. It is shown that the total number and infinitesimal stability of equilibrium solutions can change abruptly over the Sun-synchronous range of semimajor axis values (1.54 to 1.70 Earth radii). The effect of direct solar radiation pressure upon certain stable equilibria is investigated.     

Analysis of the orbit of 1965-11D (COSMOS 54 rocket)

The satellite 1965-11D was the final-stage rocket used to launch Cosmos 54, 55 and 56 into orbit on 21 February 1965. The orbit of 1965-11D was inclined at 56° to the Equator, with an initial perigee height of 280 km; the lifetime was nearly 5 yr, with decay on 23 December 1969. The orbit has been determined at 75 epochs during the life, using the RAE orbit determination program PROP with over 4000 observations, photographic, visual and radar. Observations from the Hewitt camera at Malvern were available for 34 of the 75 orbits and typical accuracies for these orbits are 0.0005° in inclination and 100 m in perigee height.The variations in perigee height have been analyzed to determine reliable values of density scale height, at heights between 240 and 360 km. The analysis also revealed a rapid decrease of 5 km in perigee distance early in 1966, attributed to the escape of residual propellants.The variations in orbital inclination have been analyzed to determine upper-atmosphere zonal winds and 15th-order harmonics in the geopotential. The region of the upper atmosphere traversed by 1965-11D near its perigee is found to have had an average rotation rate of 1.10 ± 0.05 rev/day in 1966–1967, and 1.00 ± 0.03 rev/day between March 1968 and May 1969. In late 1969 there were probably wide variations in zonal winds, with east-to-west winds of order 100 m/s followed by west-to-east winds of order 200 m/s. The changes in inclination at the 15th-order resonance in July 1969 have been analyzed to give the first accurate values of lumped 15th-order harmonics obtained from a high-drag satellite. This success points the way towards similar analyses of the many other high-drag satellites that pass through 15th-order resonance, to evaluate individual geopotential coefficients of order 15 and even degree.     

Analysis of the orbit of 1968-90A (Cosmos 248)

The satellite 1968-90A (Cosmos 248), was launched in October 1968 into an orbit inclined at 62.25° to the equator, with an initial perigee height of 475 km, apogee height 543 km, and orbital period 94.8 min. The orbit has been determined at 57 epochs over nearly one and a quarter cycles of the argument of perigee from January 1972 until December 1975 with the aid of the RAE orbit refinement program PROP, using nearly 3000 observations. For most of these orbits the standard deviations in inclination are less than 0.0009° (corresponding to about 100m in cross-track distance). The values of eccentricity give perigee heights accurate to between 30 and 120m.The main purpose of the orbit determination was to provide accurate values of the eccentricity for use in determining the odd zonal harmonics in the Earth's gravitational potential. These values have been analysed to determine the amplitude of the oscillation in eccentricity, which is found to be 0.00433 ± 0.00001.     

Solar sail time-optimal interplanetary transfer trajectory design

The fuel consumption associated with some interplanetary transfer trajectories using chemical propulsion is not affordable. A solar sail is a method of propulsion that does not consume fuel. Transfer time is one of the most pressing problems of solar sail transfer trajectory design. This paper investigates the time-optimal interplanetary transfer trajectories to a circular orbit of given inclination and radius. The optimal control law is derived from the principle of maximization. An indirect method is used...     

Search for and study of nearly periodic orbits in the plane problem of three equal-mass bodies

We analyze nearly periodic solutions in the plane problem of three equal-mass bodies by numerically simulating the dynamics of triple systems. We identify families of orbits in which all three points are on one straight line (syzygy) at the initial time. In this case, at fixed total energy of a triple system, the set of initial conditions is a bounded region in four-dimensional parameter space. We scan this region and identify sets of trajectories in which the coordinates and velocities of all bodies are close to their initial values at certain times (which are approximately multiples of the period). We classify the nearly periodic orbits by the structure of trajectory loops over one period. We have found the families of orbits generated by von Schubart’s stable periodic orbit revealed in the rectilinear three-body problem. We have also found families of hierarchical, nearly periodic trajectories with prograde and retrograde motions. In the orbits with prograde motions, the trajectory loops of two close bodies form looplike structures. The trajectories with retrograde motions are characterized by leafed structures. Orbits with central and axial symmetries are identified among the families found.     

Restriction of motions in wide pairs in the Galactic field

The motions of the components of wide binary stars in the solar neighborhood in the regular Galactic gravitational field on time scales ～10¹⁰ yr have been studied numerically. The regions of restricted motions of the components in wide pairs have been found depending on the initial conditions: the magnitude of the relative velocity of the components, their mutual distance, and the inclination of the relative velocity vector to the Galactic plane. The size of the main part of the region of restricted motions is approximately equal to the tidal radius. Profound changes in the eccentricity of the binary orbit occur at inclinations close to 90°, which can lead to close approaches of the stars with a pericenter distance less than 1 AU. In the case of retrograde motions (the binary rotates in a direction opposite to the Galactic rotation), there is a region of restricted motions extending at least to 10 pc. Examples of the trajectories of relative motion of the stars and the change in osculating orbital elements are given for systems with restricted motions.     

High-precision repeat-groundtrack orbit design and maintenance for Earth observation missions

The focus of this paper is the design and station keeping of repeat-groundtrack orbits for Sun-synchronous satellites. A method to compute the semimajor axis of the orbit is presented together with a station-keeping strategy to compensate for the perturbation due to the atmospheric drag. The results show that the nodal period converges gradually with the increase of the order used in the zonal perturbations up to \(J_{15}\). A differential correction algorithm is performed to obtain the nominal semimajor axis of the reference orbit from the inputs of the desired nodal period, eccentricity, inclination and argument of perigee. To keep the satellite in the proximity of the repeat-groundtrack condition, a practical orbit maintenance strategy is proposed in the presence of errors in the orbital measurements and control, as well as in the estimation of the semimajor axis decay rate. The performance of the maintenance strategy is assessed via the Monte Carlo simulation and the validation in a high fidelity model. Numerical simulations substantiate the validity of proposed mean-elements-based orbit maintenance strategy for repeat-groundtrack orbits.     

Analysis of the orbit of 1970-114F in its last 20 days

The satellite 1970-114F, the final-stage rocket of the Molniya 1S communications satellite, decayed in the atmosphere on 3 March 1973. During the last 20 days of its life the orbit suffered exceptionally rapid decay, with the apogee height decreasing from 7000 to 1000 km while the perigee height remained near 110 km. About 650 observations, made by visual observers in Britain and by U.S. Navy sensors, have been used with the PROP6 orbit refinement program to determine orbits at 14 epochs. Although the decay rate was more than ten times greater than in any previous orbit determination with PROP, good orbits were obtained, the standard deviation in inclination being less than 0.002° on eight orbits.The combination of high drag and good accuracy allows three techniques in orbital analysis to be successfully applied for the first time. Since zonal winds have little effect on the orbit, the changes in inclination are analysed to determine meridional winds near perigee, at heights of 110–120 km, latitudes of 63–65°S, and 6–12 hr LT. The changes in right ascension of the node are also successfully analysed for the same purpose. The two methods agree in indicating a south-to-north wind of 40 ± 30 m/sec from 11 to 21 February, a geomagnetically quiet period, and a south-to-north wind averaging 150 ± 30 m/sec from 22 February to 3 March, a geomagnetically disturbed period. Thirdly, the changes in the argument of perigee are analysed to determine atmospheric oblateness, which is found to be equal to the Earth's oblateness, to within ±20%. Lastly, the drag coefficient in transition flow is evaluated and found to be 0.85 ± 0.20.     

Navigating to the Moon along low-energy transfers

This paper presents a navigation strategy to fly to the Moon along a Weak Stability Boundary transfer trajectory. A particular strategy is devised to ensure capture into an uncontrolled relatively stable orbit at the Moon. Both uncertainty in the orbit determination process and in the control of the thrust vector are included in the navigation analysis. The orbit determination process is based on the definition of an optimal filtering technique that is able to meet accuracy requirements at an acceptable computational cost. Three sequential filtering techniques are analysed: an extended Kalman filter, an unscented Kalman filter and a Kalman filter based on high order expansions. The analysis shows that only the unscented Kalman filter meets the accuracy requirements at an acceptable computational cost. This paper demonstrates lunar weak capture for all trajectories within a capture corridor defined by all the trajectories in the neighbourhood of the nominal one, in state space. A minimum Δv strategy is presented to extend the lifetime of the spacecraft around the Moon. The orbit determination and navigation strategies are applied to the case of the European Student Moon Orbiter.     

Analysis of the orbit of cosmos 387 (1970-111A) near 15th-order resonance

Cosmos 387 (1970-111A) was launched on 16 December 1970 into a near-circular orbit with an average height of 540 km and an inclination of 74.0°. On 5 November 1971 the orbit, in its slow contraction under the influence of air drag, passed through 15th-order resonance, when the ground track repeats after 15 revolutions. The orbit has been determined with the aid of the RAE orbit refinement program PROP at 19 epochs between May 1971 and June 1972, using 1500 optical and radar observations. The average accuracy is about 70 m in perigee height and 0.001° in inclination.The variation of orbital inclination while the satellite was experiencing 15th-order resonance, as given by these 19 orbits and 55 U.S. Navy orbits, has been analysed to obtain equations accurate to 4 per cent for the geopotential coefficients of order 15 and odd degree (15, 17, 19 …). These equations have subsequently been used (with others) in determining individual coefficients of order 15 and odd degree.The variation of eccentricity with argument of perigee showed unexpected complexity, including a tight loop near resonance (Fig. 4). Analysis of the variation in eccentricity has yielded, for the first time, accurate equations for the geopotential coefficients of order 15 and even degree (16, 18 …), thus opening the way to the evaluation of individual coefficients of this type. The variations in the argument of perigee and right ascension of the node have also been analysed.     

Strategies for plane change of Earth orbits using lunar gravity and derived trajectories of family G

The dynamics of the circular restricted three-body Earth-Moon-particle problem predicts the existence of the retrograde periodic orbits around the Lagrangian equilibrium point L1. Such orbits belong to the so-called family G (Broucke, Periodic orbits in the restricted three-body problem with Earth-Moon masses, JPL Technical Report 32–1168, 1968) and starting from them it is possible to define a set of trajectories that form round trip links between the Earth and the Moon. These links occur even with more complex dynamical systems as the complete Sun-Earth-Moon-particle problem. One of the most remarkable properties of these trajectories, observed for the four-body problem, is a meaningful inclination gain when they penetrate into the lunar sphere of influence and accomplish a swing-by with the Moon. This way, when one of these trajectories returns to the proximities of the Earth, it will be in a different orbital plane from its initial Earth orbit. In this work, we present studies that show the possibility of using this property mainly to accomplish transfer maneuvers between two Earth orbits with different altitudes and inclinations, with low cost, taking into account the dynamics of the four-body problem and of the swing-by as well. The results show that it is possible to design a set of nominal transfer trajectories that require ΔV _Total less than conventional methods like Hohmann, bi-elliptic and bi-parabolic transfer with plane change.     

Analysis of the orbit of 1970-65D,cosmos 359 rocket

Cosmos 359 rocket 1970-65D, was launched on 22 August 1970 into an orbit inclined at 51·2° to the Equator, with an initial perigee height of 209 km: it decayed on 6 October 1971 after a lifetime of 410 days. The orbit has been determined at 42 epochs during the lifetime, using the RAE orbit refinement program, PROP, with over 2600 observations. Observations from the Hewitt cameras at Malvern and Edinburgh were available for 10 of the 42 orbits.Ten values of density scale height, at heights between 185 and 261 km, have been determined from analysis of the variations in perigee height.Upper-atmosphere zonal winds and 15th-order harmonics in the geopotential have been evaluated from the changes in orbital inclination. The average atmospheric rotation rate, for heights near 220 km, is found to be 1·04 rev/day; but there are striking departures from the average, with well-established values of 1·30, 0·75, 1·35 and 0·95 over four successive 75-day intervals. The changes in inclination at the 15th-order resonance in November 1970 give values of lumped 15th-order harmonics, which will provide equations for evaluating coefficients of order 15 and even degree (16,18,…) and also show that useful results on the geopotential can be obtained from satellites with perigee as low as 200 km.     

Stability of a planet in the HD 41004 binary system

The Hill stability criterion is applied to analyse the stability of a planet in the binary star system of HD 41004 AB, with the primary and secondary separated by 22 AU, and masses of 0.7 M_⊙ and 0.4 M_⊙, respectively. The primary hosts one planet in an S‐type orbit, and the secondary hosts a brown dwarf (18.64 M_J) on a relatively close orbit, 0.0177 AU, thereby forming another binary pair within this binary system. This star‐brown dwarf pair (HD 41004 B+Bb) is considered a single body during our numerical calculations, while the dynamics of the planet around the primary, HD 41004 Ab, is studied in different phase‐spaces. HD 41004 Ab is a 2.6 M_J planet orbiting at the distance of 1.7 AU with orbital eccentricity 0.39. For the purpose of this study, the system is reduced to a three‐body problem and is solved numerically as the elliptic restricted three‐body problem (ERTBP). The Hill stability function is used as a chaos indicator to configure and analyse the orbital stability of the planet, HD 41004 Ab. The indicator has been effective in measuring the planet's orbital perturbation due to the secondary star during its periastron passage. The calculated Hill stability time series of the planet for the coplanar case shows the stable and quasi‐periodic orbits for at least ten million years. For the reduced ERTBP the stability of the system is also studied for different values of planet's orbital inclination with the binary plane. Also, by recording the planet's ejection time from the system or collision time with a star during the integration period, stability of the system is analysed in a bigger phase‐space of the planet's orbital inclination, ≤ 90°, and its semimajor axis, 1.65–1.75 AU. Based on our analysis it is found that the system can maintain a stable configuration for the planet's orbital inclination as high as 65° relative to the binary plane. The results from the Hill stability criterion and the planet's dynamical lifetime map are found to be consistent with each other. (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Upper-atmosphere features revealed by the orbit of 1980-43A

The satellite NOAA-B (1980-43A) was launched in May 1980 into an orbit with perigee height near 260 km and apogee height 1440 km, at an inclination of 92.2°.The lifetime was 11 months. The orbit has been determined at 40 epochs between October 1980 and May 1981 from about 3000 radar and optical observations. The average orbital accuracy, radial and cross-track, was about 100 m, with rather better accuracy in the final 14 days.The variation of orbital inclination has been analysed to determine two good values of atmospheric rotation rate, namely 1.10 ± 0.10 rev day^?1 at 300 km (average local time) and 1.15 ± 0.06 rev day^?1 at 225 km (evening).The effect of atmospheric rotation on the precession of the orbital plane of an actual satellite has never previously been detected; it is clearly apparent for 1980-43A in its last days and conforms to the expected theoretical change.The variation of perigee height has been analysed to determine ten values of atmospheric density scale height, for heights of 280–370 km. These values, accurate to about 3%, exceed by 15% the values indicated by the COSPAR International Reference Atmosphere. Solar activity was higher in the years 1980–1981 than at any time since early 1958 and it appears that the CIRA model underestimates the density and density scale height at high levels of solar activity.     

Comparing analytical and numerical approaches to meteoroid orbit determination using Hayabusa telemetry

Fireball networks establish the trajectories of meteoritic material passing through Earth's atmosphere, from which they can derive pre‐entry orbits. Triangulated atmospheric trajectory data require different orbit determination methods to those applied to observational data beyond the Earth's sphere of influence, such as telescopic observations of asteroids. Currently, the vast majority of fireball networks determine and publish orbital data using an analytical approach, with little flexibility to include orbital perturbations. Here, we present a novel numerical technique for determining meteoroid orbits from fireball network data and compare it to previously established methods. The re‐entry of the Hayabusa spacecraft, with its known pre‐Earth orbit, provides a unique opportunity to perform this comparison as it was observed by fireball network cameras. As initial sightings of the Hayabusa spacecraft and capsule were made at different altitudes, we are able to quantify the atmosphere's influence on the determined pre‐Earth orbit. Considering these trajectories independently, we found the orbits determined by the novel numerical approach to align closer to JAXA's telemetry in both cases. Using simulations, we determine the atmospheric perturbation to become significant at ~90 km—higher than the first observations of typical meteorite dropping events. Using further simulations, we find the most substantial differences between techniques to occur at both low entry velocities and Moon passing trajectories. These regions of comparative divergence demonstrate the need for perturbation inclusion within the chosen orbit determination algorithm.     

Long-Period Variations of the Eccentricity Vector Valid also for Near Circular Orbits around a Non-Spherical Body

This paper studies the long period variations of the eccentricity vector of the orbit of an artificial satellite, under the influence of the gravity field of a central body. We use modified orbital elements which are non-singular at zero eccentricity. We expand the long periodic part of the corresponding Lagrange equations as power series of the eccentricity. The coefficients characterizing the differential system depend on the zonal coefficients of the geopotential, and on initial semi-major axis, inclination, and eccentricity. The differential equations for the components of the eccentricity vector are then integrated analytically, with a definition of the period of the perigee based on the notion of “free eccentricity”, and which is also valid for circular orbits. The analytical solution is compared to a numerical integration. This study is a generalization of (Cook, Planet. Space Sci., 14, 1966): first, the coefficients involved in the differential equations depend on all zonal coefficients (and not only on the very first ones); second, our method applies to nearly circular orbits as well as to not too eccentric orbits. Except for the critical inclination, our solution is valid for all kinds of long period motions of the perigee, i.e., circulations or librations around an equilibrium point.