The theory of asynchronous relative motion I: time transformations and nonlinear corrections

Using alternative independent variables in lieu of time has important advantages when propagating the partial derivatives of the trajectory. This paper focuses on spacecraft relative motion, but the concepts presented here can be extended to any problem involving the variational equations of orbital motion. A usual approach for modeling the relative dynamics is to evaluate how the reference orbit changes when modifying the initial conditions slightly. But when the time is a mere dependent variable, changes in the initial conditions will result in changes in time as well: a time delay between the reference and the neighbor solution will appear. The theory of asynchronous relative motion shows how the time delay can be corrected to recover the physical sense of the solution and, more importantly, how this correction can be used to improve significantly the accuracy of the linear solutions to relative motion found in the literature. As an example, an improved version of the Clohessy-Wiltshire (CW) solution is presented explicitly. The correcting terms are extremely compact, and the solution proves more accurate than the second and even third order CW equations for long propagations. The application to the elliptic case is also discussed. The theory is not restricted to Keplerian orbits, as it holds under any perturbation. To prove this statement, two examples of realistic trajectories are presented: a pair of spacecraft orbiting the Earth and perturbed by a realistic force model; and two probes describing a quasi-periodic orbit in the Jupiter-Europa system subject to third-body perturbations. The numerical examples show that the new theory yields reductions in the propagation error of several orders of magnitude, both in position and velocity, when compared to the linear approach.     

Efficient design of direct low-energy transfers in multi-moon systems

In this contribution, an efficient technique to design direct (i.e., without intermediate flybys) low-energy trajectories in multi-moon systems is presented. The method relies on analytical two-body approximations of trajectories originating from the stable and unstable invariant manifolds of two coupled circular restricted three-body problems. We provide a means to perform very fast and accurate computations of the minimum-cost trajectories between two moons. Eventually, we validate the methodology by comparison with numerical integrations in the three-body problem. Motivated by the growing interest in the robotic exploration of the Jovian system, which has given rise to numerous studies and mission proposals, we apply the method to the design of minimum-cost low-energy direct trajectories between Galilean moons, and the case study is that of Ganymede and Europa.     

Study and application of the resonant secular dynamics beyond Neptune

We use a secular representation to describe the long-term dynamics of transneptunian objects in mean-motion resonance with Neptune. The model applied is thoroughly described in Saillenfest et al. (Celest Mech Dyn Astron, doi: 10.1007/s10569-016-9700-5, 2016). The parameter space is systematically explored, showing that the secular trajectories depend little on the resonance order. High-amplitude oscillations of the perihelion distance are reported and localised in the space of the orbital parameters. In particular, we show that a large perihelion distance is not a sufficient criterion to declare that an object is detached from the planets. Such a mechanism, though, is found unable to explain the orbits of Sedna or \(2012\text {VP}_{113}\), which are insufficiently inclined (considering their high perihelion distance) to be possibly driven by such a resonant dynamics. The secular representation highlights the existence of a high-perihelion accumulation zone due to resonances of type 1:k with Neptune. That region is found to be located roughly at \(a\in [100;300]\) AU, \(q\in [50;70]\) AU and \(I\in [30;50]^{\circ }\). In addition to the flux of objects directly coming from the Scattered Disc, numerical simulations show that the Oort Cloud is also a substantial source for such objects. Naturally, as that mechanism relies on fragile captures in high-order resonances, our conclusions break down in the case of a significant external perturber. The detection of such a reservoir could thus be an observational constraint to probe the external Solar System.     

Time distributions in satellite constellation design

The aim of the time distribution methodology presented in this paper is to generate constellations whose satellites share a set of relative trajectories in a given time, and maintain that property over time without orbit corrections. The model takes into account a series of orbital perturbations such as the gravitational potential of the Earth, the atmospheric drag, the Sun and the Moon as disturbing third bodies and the solar radiation pressure. These perturbations are included in the design process of the constellation. Moreover, the whole methodology allows to design constellations with multiple relative trajectories that can be distributed in a minimum number of inertial orbits.     

On the co-orbital motion in the planar restricted three-body problem: the quasi-satellite motion revisited

In the framework of the planar and circular restricted three-body problem, we consider an asteroid that orbits the Sun in quasi-satellite motion with a planet. A quasi-satellite trajectory is a heliocentric orbit in co-orbital resonance with the planet, characterized by a nonzero eccentricity and a resonant angle that librates around zero. Likewise, in the rotating frame with the planet, it describes the same trajectory as the one of a retrograde satellite even though the planet acts as a perturbator. In the last few years, the discoveries of asteroids in this type of motion made the term “quasi-satellite” more and more present in the literature. However, some authors rather use the term “retrograde satellite” when referring to this kind of motion in the studies of the restricted problem in the rotating frame. In this paper, we intend to clarify the terminology to use, in order to bridge the gap between the perturbative co-orbital point of view and the more general approach in the rotating frame. Through a numerical exploration of the co-orbital phase space, we describe the quasi-satellite domain and highlight that it is not reachable by low eccentricities by averaging process. We will show that the quasi-satellite domain is effectively included in the domain of the retrograde satellites and neatly defined in terms of frequencies. Eventually, we highlight a remarkable high eccentric quasi-satellite orbit corresponding to a frozen ellipse in the heliocentric frame. We extend this result to the eccentric case (planet on an eccentric motion) and show that two families of frozen ellipses originate from this remarkable orbit.     

Tides in a body librating about a spin–orbit resonance: generalisation of the Darwin–Kaula theory

The Darwin-Kaula theory of bodily tides is intended for celestial bodies rotating without libration. We demonstrate that this theory, in its customary form, is inapplicable to a librating body. Specifically, in the presence of libration in longitude, the actual spectrum of Fourier tidal modes differs from the conventional spectrum rendered by the Darwin–Kaula theory for a nonlibrating celestial object. This necessitates derivation of formulae for the tidal torque and the tidal heating rate, that are applicable under libration. We derive the tidal spectrum for longitudinal forced libration with one and two main frequencies, generalisation to more main frequencies being straightforward. (By main frequencies we understand those emerging due to the triaxiality of the librating body.) Separately, we consider a case of free libration at one frequency (once again, generalisation to more frequencies being straightforward). We also calculate the tidal torque. This torque provides correction to the triaxiality-caused physical libration. Our theory is not self-consistent: we assume that the tidal torque is much smaller than the permanent-triaxiality-caused torque, so the additional libration due to tides is much weaker than the main libration due to the permanent triaxiality. Finally, we calculate the tidal dissipation rate in a body experiencing forced libration at the main mode, or free libration at one frequency, or superimposed forced and free librations.     

Non-resonant secular dynamics of trans-Neptunian objects perturbed by a distant super-Earth

We use a secular model to describe the non-resonant dynamics of trans-Neptunian objects in the presence of an external ten-Earth-mass perturber. The secular dynamics is analogous to an “eccentric Kozai mechanism” but with both an inner component (the four giant planets) and an outer one (the eccentric distant perturber). By the means of Poincaré sections, the cases of a non-inclined or inclined outer planet are successively studied, making the connection with previous works. In the inclined case, the problem is reduced to two degrees of freedom by assuming a non-precessing argument of perihelion for the perturbing body. The size of the perturbation is typically ruled by the semi-major axis of the small body: we show that the classic integrable picture is still valid below about 70 AU, but it is progressively destroyed when we get closer to the external perturber. In particular, for \(a>150\) AU, large-amplitude orbital flips become possible, and for \(a>200\) AU, the Kozai libration islands at \(\omega =\pi /2\) and \(3\pi /2\) are totally submerged by the chaotic sea. Numerous resonance relations are highlighted. The most large and persistent ones are associated with apsidal alignments or anti-alignments with the orbit of the distant perturber.     

Performance results of HESP physical model

As a continuation to the published work on model based calibration technique with HESP(Hanle Echelle Spectrograph) as a case study, in this paper we present the performance results of the technique. We also describe how the open parameters were chosen in the model for optimization, the glass data accuracy and handling the discrepancies. It is observed through simulations that the discrepancies in glass data can be identified but not quantifiable. So having an accurate glass data is important which is possible to obtain from the glass manufacturers. The model’s performance in various aspects is presented using the ThAr calibration frames from HESP during its pre-shipment tests. Accuracy of model predictions and its wave length calibration comparison with conventional empirical fitting, the behaviour of open parameters in optimization, model’s ability to track instrumental drifts in the spectrum and the double fibres performance were discussed. It is observed that the optimized model is able to predict to a high accuracy the drifts in the spectrum from environmental fluctuations. It is also observed that the pattern in the spectral drifts across the 2D spectrum which vary from image to image is predictable with the optimized model. We will also discuss the possible science cases where the model can contribute.     

Bandpass calibration of a wideband spectrometer using coherent pulse injection

We present a relatively simple time domain method for determining the bandpass response of a system by injecting a nanosecond pulse and capturing the system voltage output. A pulse of sub-nanosecond duration contains all frequency components with nearly constant amplitude up to 1 GHz. Hence, this method can accurately determine the system bandpass response to a broadband signal. In a novel variation on this impulse response method, a train of pulses is coherently accumulated providing precision calibration with a simple system. The basic concept is demonstrated using a pulse generator-accumulator setup realised in a Bedlam board which is a high speed digital signal processing unit. The same system was used at the Parkes radio telescope between 2–13 October 2013 and we demonstrate its powerful diagnostic capability. We also present some initial test data from this experiment.