Gravitation Maneuver Using the Families of Super-Unstable Orbits around the Libration Points

The families of periodic solutions to an autonomous Hamiltonian system in that part where the solutions are unstable have their specific field of influence. Under strong instability, the orbits that have fallen in such a field of influence are drawn into the family as in a whirlpool and then are thrown away from it. In the particular case of the restricted three-body problem, the orbits around the libration points L ₁ and L ₂ are super-unstable and the Keplerian elements in motion over these orbits change drastically. When the orbits fall into such a domain, for some time the motion is close to the motion around the libration point, and after being thrown out of this domain, the Keplerian orbital elements also change essentially.     

Collision Risk Mitigation in Geostationary Orbit

The short- and long-term effects of spacecraft explosions, as a function of the end-of-life re-orbit altitude above the geostationary orbit (GEO), were analyzed in terms of their additional contribution to the debris flux in the GEO ring. The simulated debris clouds were propagated for 72yrs, taking into account all the relevant orbital perturbations.The results obtained show that 6–7 additional explosions in GEO would be sufficient, in the long term, to double the current collision risk with sizable objects in GEO. Unfortunately, even if spacecraft were to re-orbit between 300 and 500km above GEO, this would not significantly improve the situation. In fact, an altitude increase of at least 2000km would have to be adopted to reduce by one order of magnitude the long-term risk of collision among geostationary satellites and explosion fragments. The optimal debris mitigation strategy should be a compromise between the reliability and effectiveness of spacecraft end-of-life passivation, the re-orbit altitude and the acceptable debris background in the GEO ring. However, for as long as the re-orbit altitudes currently used are less than 500km above GEO, new spacecraft explosions must be avoided in order to preserve the geostationary environment over the long term.     

Stationary Nose Structures of Protons in the Inner Magnetosphere: Observations by the ION Instrument onboard the Interball-2 Satellite and Modeling

Nose structures are objects formed by H⁺ particles penetrating into the inner magnetosphere [1, 2]. We present the results of experimental studies and numerical modeling of the nose structures. Statistical processing of the observations of nose structures in 1997 by the ION instrument onboard the Interball-2 satellite at heights of 10000–15000 km demonstrates that the probability of formation of the nose structures under quiet magnetic conditions (with current values K _p = 0–1) in the nighttime sector of the magnetosphere is 90%. The probability of observation of the nose structures in the daytime sector equals 50% at the current value K _p = 0–1, and the correlation between the observations of nose structures and K _p can be improved (up to 75%) if the K _p index is taken 6 h before the observed events. It is shown that nose structures are a characteristic feature not only of the substorm processes but also of quasi-stationary phenomena in the quiet magnetosphere. The nose structures observed in magnetically quiet periods are called stationary nose structures in this work. By modeling drift trajectories for protons, it is shown that the stationary nose structures are formed in all sectors of the MLT. The stationary nose structures observed by the ION instrument are modeled in the night, morning, and daytime sectors of the MLT. The relation between the stationary nose structures and ion spectral gaps is considered.     

Formation of a Polarization Jet during the Expansion Phase of a Substorm: Results of Ground-Based Measurements

The ground-based polarization jet measurements at the Yakutsk ionosphere station (L= 3.0) for the years 1989–1991 (110 events) are compared with variations of the AE-index and with parameters of the local magnetic activity. It is shown that polarization jet development in the near midnight sector can occur during a period of no longer than 10 min on the expansion phase of a substorm. The formation of the polarization jet is accompanied by a specific magnetic field variation corresponding in shape to a fast passage of the Harang Discontinuity above the station. Statistical data are given on ground level observations of the polarization jet, which are close to those measured from satellites. The mean delay (averaged over the full data bank) between the onset of a substorm with AE 500 nT and the moment of the polarization jet appearance at L= 3.0 is equal to 0.5 h near midnight and to 1.0 – 1.5 h in the evening sector. Estimations show that the duration of the polarization jet formation when energetic ions are injected into the Harang Discontinuity region above the ground station can last for about 10 min, and during this time the Harang Discontinuity can be shifted to the west. This is in qualitative agreement with the described observations.     

Assessing the Impact Risk of Orbital Debris on Space Tethers

Tethers are being proposed for a growing number of space applications. However, they may be particularly vulnerable to orbital debris and meteoroid impacts. In order to provide useful reference data for tether systems design, detailed analytical and numerical computations were carried out to assess the average impact rate of artificial debris and meteoroids. The specific geometric properties of tethers as debris targets, when compared to typical satellites, are discussed, and the results obtained are presented in tabular form, as a function of debris size and tether diameter.The computations were carried out for six circular orbits, spanning three altitudes (600, 800 and 1000km) and two inclinations (30° and 50°). Tether diameters in between 1mm and 2cm and debris larger than 0.1mm were considered in the analysis. The collision risk of tethers with spacecraft and upper stages in orbit was estimated as well.In the debris interval and orbital regimes considered, artificial debris represent the dominant contributor to the impact rate. At 600km and in the 0.1–10mm size range, the meteoroid and orbital debris impact rates are still comparable; however, at higher altitudes and in the 1–10cm size range, meteoroids contribute 20–30 times less to the collision probability.The results obtained confirm that for single-strand tethers in low Earth orbit the probability to be severed by orbital debris and meteoroid impacts is quite significant, making necessary the adoption of innovative designs for long duration missions.     

Tether Satellite System Collision Study

A study was performed to determine the probability of collision with resident space objects and untrackable debris for the tether component of the Tethered Satellite System (TSS) after it broke away from the Space Shuttle orbiter (mission STS-75) in February 1996. Both an analytical and a numerical approach were used in this study, and the results obtained with these two methods were found to be in good agreement. These results show that the deployed tether is expected to have been impacted by several particles 0.1mm or larger in size. The probability of collision with objects 10cm in size or larger was on the order of 10^–3 per month. Since the severed tether reentered within one month after deployment, the collision hazard to other objects while in orbit was small. The analytical methods used in this study are useful for tether collision evaluations in general.     

基于碰撞概率的交会对接近距离导引段的轨迹安全

提出了一种利用碰撞概率思想来描述交会对接近距离导引段飞行器轨迹安全的方法,这种方法采用最大瞬时碰撞概率来定量描述两个飞行器的碰撞危险程度。推导了最大瞬时碰撞概率的计算公式,其过程可分为三步:首先基于CW方程外推初始相对状态和误差协方差矩阵,然后利用外推后的相对状态和误差协方差矩阵计算瞬时碰撞概率,最后求出感兴趣区间内的最大瞬时碰撞概率。数值算例表明最大瞬时碰撞概率是描述交会对接近距离段两个飞行器碰撞危险程度的一个较好的指标。     

Dynamics of the 0.3–100 MeV Proton Energy Spectra in the Inner Heliosphere in Quiet Periods of 1974–1991

Variations of the proton spectra in the 0.3–100 MeV energy range based on the data of various instruments installed onboard the IMP-8 satellite are studied for very quiet, quiet, and quasi-stable solar activity periods during the years 1974–1991. As many as 118 spectra were approximated by two power laws: the left-hand and galactic branches of the spectrum were fitted by the AE ^– function and a dependence of the CE type, respectively, the sum J(E) = AE ^– + CE providing the total spectrum. It is shown that the spectra vary within a solar cycle with a shift of the minimum energy (E _min) to higher energies with increasing solar activity. It follows from the relations between the spectrum parameters thus obtained that, in particular periods of time, an increase (decrease) of the particle flux in the low-energy branch of the spectrum and an intensification (depression) of the GCR particle flux modulation take place simultaneously. This is manifested in a shift of the spectrum parallel to the energy axis. The study of the spectra in the most quiet time during three successive solar minima have shown that low-energy (0.3–10 MeV) protons, as well as GCR, are subject to the 22-year variation in the solar magnetic cycle.     

Empirical Modeling of the Statistical Structure of Radio Signals from Satellites Moving over Mid- and High-Latitude Trajectories in the Southern Hemisphere

The results of developing the empirical model of parameters of radio signals propagating in the inhomogeneous ionosphere at middle and high latitudes are presented. As the initial data we took the homogeneous data obtained as a result of observations carried out at the Antarctic Molodezhnaya station by the method of continuous transmission probing of the ionosphere by signals of the satellite radionavigation Transit system at coherent frequencies of 150 and 400 MHz. The data relate to the summer season period in the Southern hemisphere of the Earth in 1988–1989 during high (F > 160) activity of the Sun. The behavior of the following statistical characteristics of radio signal parameters was analyzed: (a) the interval of correlation of fluctuations of amplitudes at a frequency of 150 MHz ( _kA ); (b) the interval of correlation of fluctuations of the difference phase ( _k); and (c) the parameter characterizing frequency spectra of amplitude (P _A) and phase (P ) fluctuations. A third-degree polynomial was used for modeling of propagation parameters. For all above indicated propagation parameters, the coefficients of the third-degree polynomial were calculated as a function of local time and magnetic activity. The results of calculations are tabulated.     

非线性相对运动的飞行器碰撞概率研究

发展了一种非线性相对运动情况下的飞行器碰撞概率描述和计算方法。采用拟最大瞬时碰撞概率和总碰撞概率两个主要指标来定量描述飞行器碰撞的危险性，并针对这两个指标发展了相应的计算方法。对于拟最大瞬时碰撞概率的计算，首先需要在相对轨迹上寻找最小相对距离或最大概率密度的点，然后在这个点上计算拟最大瞬时碰撞概率。对于总碰撞概率的计算，由于总碰撞概率的增量可以表示成一个时间的函数，因此总碰撞概率的计算就是一个在时间域上的积分。所建立的这种描述和计算方法适用于相对速度和误差协方差矩阵随时间任意变化的情况。数值算例表明所提出的指标比单纯采用总碰撞概率能够更好的描述飞行器的碰撞危险性，并且碰撞概率的计算方法和线性相对运动情况下的计算方法是兼容的。     

Measurement of Vibration Microaccelerations onboard the Mir Orbital Station Using the VM-09 System of Accelerometers

The results of processing the data of measurements of microaccelerations, carried out onboard the Mir orbital station using the Russian VM-09 system of accelerometers, are described. The system was developed by the Composite Research-Production Association. The sensitivity of this system was 10^–4 m/s²; its frequency band had limits from a few tenths of a hertz up to 100 Hz. The measurements were carried out in the real-time mode of data transmission to the Earth, when the orbital station flew over the telemetry data receiving point. The instrument's sampling rate was 200 measurements per second, and the length of a continuous run of measurements did not exceed 10 min. The following problems are considered in the paper: (1) isolation of cyclic trends from the measurement data; (2) estimation of spectral density of the data component with a continuous spectrum; and (3) low-frequency filtration of the measurement data     

Developing a Low-Velocity Collision Model Based on the Nasa Standard Breakup Model

We have conducted a series of low-velocity impact experiments to understand the dispersion properties of fragments newly created by low-velocity impacts possible in space, especially in geostationary Earth orbit. The test results are utilized to establish a mathematical prediction model to be used in debris generation and propagation codes. Since the expected collision velocity between catalogued objects in geostationary Earth orbit shows a peak at a few hundreds meters per second, these impact experiments were conducted at a velocity range lower than 300m/s. As a typical structure of satellites in geostationary Earth orbit, thin aluminum honeycomb sandwich panels with carbon fiber reinforced plastics face sheets were prepared, while the projectile was a stainless steel ball of 9mm diameter. The data collected through these impact experiments have been re-analyzed based on the method used in the National Aeronautics and Space Administration (NASA) standard breakup model 1998 revision. The results indicate that the NASA standard breakup model derived from hypervelocity impacts could be applied to low-velocity collision possible in geostationary Earth orbit with some modifications.     

Application of LB-Planning for Designing the Programs of Trajectory Measurements

The LB-planning is a construction of the dependence between three variables k, m, and , where k is the number of parameters of the linear dynamic system (LDS) under study, m is the specified number of measurements (communication sessions), and (defect) is the maximum value of failure (set of all criterion functions) in a chosen m-point measurement program relative to the optimal program. The problem is reduced to the solution of an extremal problem, which, at m = k + 1, is found rather simply by maximization of the Tchebyshev determinant (TD-planning).     

Ray Based Diffraction Tomography of the Ionosphere and Laboratory Inhomogeneous Plasma

A new procedure for restoration of the plasma inhomogeneities with improved resolution is suggested. The procedure deals with the double weighted Fourier transform (DWFT) of the observed wavefield in coordinates of both receivers = (x, y) and sources ₀ = (x ₀, y ₀) [1]. Phase increments between the sources and receivers, being found from DWFT representation, can be used for extracting information on small perturbations of the dielectric constant ~(, z) in a way similar to traditional radio tomography. The resulting resolution of the method is close to the diffraction limit = h/D in the horizontal direction and z = (h/D)² in the vertical direction, where h is the height of inhomogeneities and D is the length of the ground-based receiving system.     

New insights on the orbital debris collision hazard at GEO

An analysis is performed of the orbital debris collision hazard to operational spacecraft at geosynchronous orbit (GEO). As part of the examination, the contribution of individual components of the population are considered and presented to provide a clearer linkage between object characteristic and resulting risk. Our examination of GEO collision risk reveals several critical new insights: (1) the current probability of collision in GEO is relatively low, yet the future is difficult to predict due to our limited ability to observe objects in GEO and the uncertainty in past and future debris-generating events in GEO; (2) the probability of collision in GEO is not uniform by longitude — it is seven times greater in regions centered about the geopotential wells; (3) the probability of a mission-terminating collision is greatly dependent upon the approximately 2200 objects in the 10 cm–1 m range observed in GEO but not yet cataloged; (4) hardware relocated to GEO “graveyard” disposal orbits pose a potential additional, but not fully understood, collision hazard to operational GEO satellites; and (5) the collision hazard throughout the course of a day or year is highly episodic (i.e. non-uniform).     

Advanced technologies available for future solid propellant grains

Significant advances have been made during the last decade in several fields of solid propulsion: the advances have enabled new savings in the motor development phase and recurring costs, because they help limit the number of prototypes and tests.The purpose of the paper is to describe the improvements achieved by SNPE in solid grain technologies, making these technologies available for new developments in more efficient and reliable future SRMs: new energetic molecules, new solid propellants, new processes for grain manufacturing, quick response grain design tools associated with advanced models for grain performance predictions.Using its expertise in chemical synthesis, SNPE develops new molecules to fit new energetic material requirements.Tests based on new propellant formulations have produced good results in the propellant performance/safety behavior ratio. New processes have been developed simultaneously to reduce the manufacturing costs of the new propellants.In addition, the grain design has been optimized by using the latest generation of predictive theoretical tools supported by a large data bank of experimental parameters resulting from over 30 years' experience in solid propulsion:

• Computer-aided method for the preliminary grain design

• Advanced models for SRM operating and performance predictions

Analysis and implementation of in-plane stationkeeping of continuously perturbed Walker constellations

The stationkeeping of symmetric Walker constellations is analyzed by considering the perturbations arising from a high order and degree Earth gravity field and the solar radiation pressure. These perturbations act differently on each group of spacecraft flying in a given orbital plane, causing a differential drift effect that would disrupt the initial symmetry of the constellation. The analysis is based on the consideration of a fictitious set of rotating reference frames that move with the spacecraft in the mean sense, but drift at a rate equal to the average drift rate experienced by all the vehicles over an extended period. The frames are also allowed to experience the $J_2$-precession such that each vehicle is allowed to drift in 3D relative to its frame. A two-impulse rendezvous maneuver is then constructed to bring each vehicle to the center of its frame as soon as a given tolerance deadband is about to be violated. This paper illustrates the computations associated with the stationkeeping of a generic Walker constellation by maneuvering each leading spacecraft within an orbit plane and calculating the associated velocity changes required for controlling the in-plane motions in an exacting sense, at least for the first series of maneuvers. The analysis can be easily extended to lower flying constellations, which experience additional perturbations due to drag.     

Reconstruction of spacecraft rotational motion using a Kalman filter

Quasi-static microaccelerations of four satellites of the Foton series (nos. 11, 12, M-2, M-3) were monitored as follows. First, according to measurements of onboard sensors obtained in a certain time interval, spacecraft rotational motion was reconstructed in this interval. Then, along the found motion, microacceleration at a given onboard point was calculated according to the known formula as a function of time. The motion was reconstructed by the least squares method using the solutions to the equations of satellite rotational motion. The time intervals in which these equations make reconstruction possible were from one to five orbital revolutions. This length is increased with the modulus of the satellite angular velocity. To get an idea on microaccelerations and satellite motion during an entire flight, the motion was reconstructed in several tens of such intervals. This paper proposes a method for motion reconstruction suitable for an interval of arbitrary length. The method is based on the Kalman filter. We preliminary describe a new version of the method for reconstructing uncontrolled satellite rotational motion from magnetic measurements using the least squares method, which is essentially used to construct the Kalman filter. The results of comparison of both methods are presented using the data obtained on a flight of the Foton M-3.     

碰撞概率对轨道误差和接近几何条件的灵敏度分析

基于显式表达式对碰撞概率的灵敏度进行解析分析,给出其对接近距离的RSW分量、轨道误差、接近角度和目标大小的两类灵敏度（绝对灵敏度和相对灵敏度）。以美俄卫星碰撞事件为例进行计算和分析,得到灵敏度的变化曲线和数值。该解析方法可以克服数值分析方法中普遍存在的使用不便、计算量大、精度不高等问题。