Study of highly perturbed spacecraft formation dynamics via approximation

This paper explores methods for approximating and analyzing the dynamics of highly perturbed spacecraft formations with an emphasis on computationally efficient approaches. This facilitates on-board computation or rapid preliminary mission design analysis. Perturbed formation dynamics are often approximated as linear time-varying (LTV) systems, for which Floquet theory can be used to analyze the degree of system instability. Furthermore, the angular momentum of the relative orbital state can be computed with the approximate dynamics to provide additional insight. A general methodology is developed first and then applied to the problem of unstable formation dynamics in asteroid orbits. Here the dominant perturbative effects due to low-order gravitational harmonics and solar radiation pressure are modeled. Numerical simulations validate the approach and illustrate the approximation accuracy achieved.     

Estimating the spin axis orientation of the Echostar-2 box-wing geosynchronous satellite

For the first time, the spin axis orientation of an inactive box-wing geosynchronous satellite has been estimated from ground-based optical photometric observations of Echostar-2’s specular reflections. Recent photometric light curves obtained of Echostar-2 over four years suggest that unusually bright and brief specular reflections were occurring twice within an observed spin period. These bright and brief specular reflections suggested two satellite surfaces with surface normals separated by approximately 180°. The geometry between the satellite, the Sun, and the observing location at the time of each of the brightest observed reflections, was used to estimate Echostar-2’s equatorial spin axis orientation coordinates. When considering prograde and retrograde rotation, Echostar-2’s spin axis orientation was estimated to have been located within 30° of either equatorial coordinate pole. Echostar-2’s spin axis was observed to have moved approximately 180° in right ascension, within a time span of six months, suggesting a roughly one year spin axis precession period about the satellite’s angular momentum vector.     

Mapping mean lake surface from satellite altimetry and GPS kinematic surveys

Lake water height is a key variable in water cycle and climate change studies, which is achievable using satellite altimetry constellation. A method based on data processing of altimetry from several satellites has been developed to interpolate mean lake surface (MLS) over a set of 22 big lakes distributed on the Earth. It has been applied on nadir radar altimeters in Low Resolution Mode (LRM: Jason-3, Saral/AltiKa, CryoSat-2) in Synthetic Aperture Radar (SAR) mode (Sentinel-3A), and in SAR interferometric (SARin) mode (CryoSat-2), and on laser altimetry (ICESat). Validation of the method has been performed using a set of kinematic GPS height profiles from 18 field campaigns over the lake Issykkul, by comparison of altimetry’s height at crossover points for the other lakes and using the laser altimetry on ICESat-2 mission. The precision reached ranges from 3 to 7 cm RMS (Root Mean Square) depending on the lakes. Currently, lake water level inferred from satellite altimetry is provided with respect to an ellipsoid. Ellipsoidal heights are converted into orthométric heights using geoid models interpolated along the satellite tracks. These global geoid models were inferred from geodetic satellite missions coupled with absolute and regional anomaly gravity data sets spread over the Earth. However, the spatial resolution of the current geoid models does not allow capturing short wavelength undulations that may reach decimeters in mountaineering regions or for rift lakes (Baikal, Issykkul, Malawi, Tanganika). We interpolate in this work the geoid height anomalies with three recent geoid models, the EGM2008, XGM2016 and EIGEN-6C4d, and compare them with the Mean Surface of 22 lakes calculated using satellite altimetry. Assuming that MLS mimics the local undulations of the geoid, our study shows that over a large set of lakes (in East Africa, Andean mountain and Central Asia), short wavelength undulations of the geoid in poorly sampled areas can be derived using satellite altimetry. The models used in this study present very similar geographical patterns when compared to MLS. The precision of the models largely depends on the location of the lakes and is about 18 cm, in average over the Earth. MLS can serve as a validation dataset for any future geoid model. It will also be useful for validation of the future mission SWOT (Surface Water and Ocean Topography) which will measure and map water heights over the lakes with a high horizontal resolution of 250 by 250 m.     

Advanced Pulsed Plasma Thrusters and Their Application as a Part of Small Spacecraft Propulsion Systems

The paper presents the research results of the effect of a capacitor energy storage device configuration on the specific characteristics of advanced modern propulsion systems based on the ablative pulsed plasma thrusters (APPT). These thrusters are designed to perform specific tasks within the small spacecrafts with the onboard power capacity up to 200 W.     

Optical Design for a Cubesat: Unobscured Telescope,Using Freeform Mirrors and a Curved Detector

The possibility of constructing an optical scheme of a telescope for a small satellite based on mirrors with freeform surfaces and a curved photodetector is considered. It is shown that the use of a new element base allows the effective area of the main mirror to be maximized and the field of view up to 9.4 times to be increased in comparison with the classical scheme of the Ritchey–Chretien scheme in ensuring the high image quality.     

Features of whistling electromagnetic wave propagation descending from above upon the night ionosphere

The results of numerical solution of the wave equations for the oblique incidence of whistling electromagnetic waves upon the night ionosphere from above have been obtained and analyzed. In the studied region of altitudes, within the wavelength scale, charged particle concentration varies drastically, and damping caused by collisions between the charged and neutral particles decreases considerably. Below, the sharp lower boundary of the ionosphere, the refractive index of the whistler wave approaches unity, and plasma turbulence transform into atmospheric electromagnetic waves. The dependences of the whistler reflection factor are found in terms of energy and horizontal magnetic component of the electromagnetic wave near the Earth’s surface on the frequency and the wave vector transverse component for the plain-layered medium model at two values of latitude. Strong dependences have been found on the wave angle of incidence and frequency. At rather small angles of incidence, the wave disturbances reach the Earth’s surface, and the module of reflection coefficient logarithm is in the range of 0.4–1. At large angles of incidence, the reflection coefficient module varies over a wide range depending on specific conditions. The obtained results explain the absence of oscillation modes of plasma magnetosphere maser in the night magnetosphere.     

The Role and Contributions of Energetic Neutral Atom (ENA) Imaging in Magnetospheric Substorm Research

Energetic Neutral Atom (ENA) imaging has contributed substantially to substorm research. This technique has allowed significant advances in areas such as observation and quantification of injected particle drift as a function of energy, observation of dynamics in the tail that are directly related to the effects of imposed (growth phase) and induced (expansion phase) electric fields on the plasma, the prompt extraction of oxygen from the ionosphere during substorms, the relationship between storms and substorms, and the timing of substorm ENA signatures. We present discussion of the advantages and shortcomings of the ENA technique for studying space plasmas. Although the technique is in its infancy, it is yielding results that enrich our understanding of the substorm process and its effects.     

Microthrusters in silicon for aerospace application

This paper reviews the results of the thermal and static analysis of small motor aerospace technology (SMART) propulsion system, constituted of a microthrusters array realised by MEMS technology on silicon wafers. This system has been studied using FEM (NASTRAN) and the results have been verified by the electro-thermic analogy and the FDM method, using, respectively, SPICE and MATLAB codes. The simulation results demonstrated the feasibility of SMART systems for aerospace applications such as attitude control and deorbiting missions for small satellite station-keeping. A theoretical impulse of 20 mNs has been calculated for the SMART system.     

Elemental and Isotopic Abundances of Carbon and Nitrogen in Meteorites

We have taken an inventory of the elemental and isotopic abundances of major carbon- and nitrogen-bearing components in different groups of meteorites. Primary phases, inherited from the solar nebula, are frequently isotopically heterogeneous, and surprisingly resistant to modification through parent body processing. Even melted and recrystallised meteorites retain primordial carbon and nitrogen isotopic signatures. This revised version was published online in August 2006 with corrections to the Cover Date.     

Influence of the 2:1 resonance in the orbits of GPS satellites

Due to the characteristics of their orbits the GPS satellites are submitted to the following main perturbations: terrestrial gravitational field, luni-solar gravitational attraction and solar radiation pressure (including the effects of the Earth's shadow). An additional perturbation arises due to the 2:1 commensurability of the orbital period of the satellite with the period of the Earth's rotation. An analytical theory is briefly presented to solve the equations of motion including the previously mentioned effects. The analytical solution, based on the Lie-Hori method, is compared with a numerical integration of the equations.