Neutralization of negative electric charge of a satellite by ionospheric plasma ions

The dynamics of the process of neutralization of a satellites negative electric charge by the charge of a layer of positive ions of ionospheric plasma is investigated. The changes in the electric field strength, as well as in the velocity and density of plasma ions, are obtained depending on the initial position of a particle and the distance to the satellites surface. It is found that under the conditions of accepted approximations, as the time grows unlimitedly, the satellites charge tends to zero but does not reach this value. The time during which the satellites electric charge decreases in magnitude as many times as specified is determined. It is shown that the negative charge of a satellite and the charge of a layer of plasma ions can form a structure of a double charged layer type whose lifetime is much longer than the time of its formation.Translated from Kosmicheskie Issledovaniya, Vol. 43, No. 1, 2005, pp. 9–18.Original Russian Text Copyright © 2005 by Fedorov.     

Stochastic models of the earth’s rotation taking the effects of the moon and planets into account

A generalization of kinetic and spectral and correlation equations of fluctuations in the Earths rotation is made taking into account the effect of the Moon and planets (these equations form the basis for corresponding stochastic models). Special attention is paid to a parameterized version of kinetic equations based on the method of normal approximation. The application of linear stochastic models for obtaining the spectral and correlation characteristics of the Earths rotation is considered, including the case of random parameters of the tensor of inertia of a deformable Earth. The results can be used for constructing stochastic models of the Earths rotation over time intervals of 6–10 years.Translated from Kosmicheskie Issledovaniya, Vol. 43, No. 1, 2005, pp. 54–66.Original Russian Text Copyright © 2005 by Markov, Dasaev, Perepelkin, I. N. Sinitsyn, V. I. Sinitsyn.     

Evolution of Rotational Motion of a Symmetrical Satellite with Viscoelastic Rods

We study the motion of a symmetrical satellite with a pair of flexible viscoelastic rods in a central Newtonian gravitational field. A restricted problem formulation is considered, when the satellite's center of mass moves along a fixed circular orbit. A small parameter is introduced which is inversely proportional to the stiffness of flexible elements. Another small parameter is equal to the ratio of the squared orbital angular velocity and the squared magnitude of the initial angular velocity of the satellite. In order to describe the satellite rotational motion relative to the center of mass, we use the canonical Andoyer variables. In the undisturbed formulation of the problem, i.e., at = 0 and = 0, these variables are the action–angle variables. Equations describing the evolution of motion are derived by an asymptotic method which combines the method of separating motions for systems with an infinite number of degrees of freedom and the Krylov–Bogolyubov method for systems with fast and slow variables. The manifolds of stationary motions are found, and their stability is investigated on the basis of equations in variations. Phase portraits are constructed which describe the rotational motion of a satellite at the stage of slow dissipative evolution.     

Some Technical Issues of an Optically Focused Small Space Debris Tracking and Cataloging System

Congressional language in the 1998 US Senate Armed Services Committee authorization bill directed ... the Secretary of the (United States) Air Force to undertake a design study of a system that could catalog and track debris down to one centimeter in size out to 1000kilometer in altitude. The US Air Force Research Laboratory, in conjunction with other US National Laboratories and the National Aeronautics and Space Administration (NASA) conducted a study that examined what technical systems and operations would be required to perform such a mission. This paper outlines the study process, details the findings, draws conclusions, and makes recommendations as to what would be needed to develop an optically based system capable of cataloging and tracking small debris in low Earth orbit.     

Quasistatic Microaccelerations onboard the Foton Satellite during Its Flight at High-Altitude Orbits

The level of quasistatic microaccelerations onboard the Foton-M satellite is predicted for its flights in two orbits: the planned orbit with the altitudes in perigee h = 262 km and in apogee h = 304 km and the orbit with h = 262 km and h = 350 km. The prediction is based on mathematical simulation of the satellite motion with respect to its center of mass under the action of gravitational and aerodynamic moments. The model is represented by the system of equations of the satellite rotational motion. Parameters of this system are chosen from the condition of coincidence of the motion of preceding Foton satellites (h 220 km and h 400 km) calculated using this model with the results of determination of actual rotational motion of the Foton-11 and Foton-12 satellites. With the help of the model thus calibrated, a calculation is made of the rotational motion of the Foton-M satellite and of the quasistatic microaccelerations onboard it. As is shown by the results of simulation, the use of the first and the second orbits will result in reductions of microaccelerations by 30% and 60%, respectively.     

The role of horizontal components of the neutral wind velocity in anthropogenic impacts on the earth’s ionosphere

The influence of zonal and meridional components of the neutral wind under anthropogenic impacts of water molecules on the ionospheres electron concentration distributed along a magnetic field tube is considered. Computer simulation method based on a mathematical model of the F2 region along a field tube is used.Translated from Kosmicheskie Issledovaniya, Vol. 43, No. 1, 2005, pp. 76–80.Original Russian Text Copyright © 2005 by Medvedev, Ishanov, Zenkin.     

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.     

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.     

Energy Spectrum of Galactic 10–100 MeV Protons in Quiet Sun Periods

The dynamics of the proton energy spectrum during the solar cycle is studied. The spectra were determined by 1–100 MeV particle fluxes measured by different instruments mounted aboard the Earth's IMP-8 satellite for more than one hundred quiet-time intervals in the period between 1974 and 1991. The galactic branch of the spectra (E _p > 10 MeV) constructed for every quiet interval was fitted by a power law function, J =CE . The theory predicts that in the 1–100 MeV energy range, where the adiabatic cooling of particles is dominant, = 1, while we have derived a double-peak distribution. The main maximum has the mean value = 1.35. The mean value of the second, much weaker maximum, is = 0.95. Within the main maximum, values are distributed in accordance with the Gaussian law with a standard deviation D/ = 0.12. The substantial difference of from unity requires the elaboration of a new model of modulation processes in the inner heliosphere. The values corresponding to the second maximum show that modulation processes correspond sometimes to theoretical conceptions. It is shown that correlates weakly with parameters A and describing the solar branch of the spectrum (J(E) = AE ^–). At the same time, a more significant correlation is observed between and the solar activity index, R _z, the counting rate of the Deep River neutron monitor, and the energy value in the minimum of the energy spectrum flux, E _min.     

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.     

Multisatellite systems with linear structure and their application for continuous coverage of the earth

Multisatellite systems with linear structure (SLS) are defined, and their application for a continuous global or zonal coverage of the Earths surface is justified. It is demonstrated that in some cases these systems turned out to be better than usually recommended kinematically regular systems by G.V. Mozhaev, delta systems of J.G. Walker, and polar systems suggested by F.W. Gobets, L. Rider, and W.S. Adams. When a comparison is made using the criterion of a minimum radius of one-satellite coverage circle, the SLS beat the other systems for the majority of satellite numbers from the range 20–63, if the global continuous single coverage of the Earth is required. In the case of a zonal continuous single coverage of the latitude belt ±65°, the SLS are preferable at almost all numbers of satellites from 38 to 100, and further at any values up to 200 excluding 144.Translated from Kosmicheskie Issledovaniya, Vol. 43, No. 1, 2005, pp. 36–53.Original Russian Text Copyright © 2005 by Saulskiy.     

Global geomagnetic response to a sharp compression of the magnetosphere and IMF variations on October 29, 2003

We study the response of the ionosphere and magnetosphere to a sudden commencement (SC) on October 29, 2003, at 06:11 UT. It is shown that the geomagnetic response had the form of two successive stages. In the first 5 min after the SC, a strong intensification of a two-vortex current system of the DP2 type was observed in latitudes 67°-65°, with variations of H -4000 nT (+700 nT). At the same time, energetic electrons were injected without dispersion to geosynchronous orbits simultaneously in the sectors 16, 04, and 07 MLT. In the subsequent 5–15 min, a new intensification of the western electrojet took place in all time sectors at latitudes 70°. Around midnight, this electrojet was extended in the poleward direction up to the polar cap latitudes ( 75°-83°). It had an unusually high velocity of extension (up to 5.0 km/s) and was accompanied by typical dispersionless substorm injections, but only at meridians 04 and 07 MLT. From comparing the development of electrojets with the data of satellite observations in the solar wind and magnetosphere, we suggest that 3–5 min after the SC onset a dipolization of the magnetic field at the geosynchronous orbit occurred. It was connected with the decay of the current flowing across the magnetotail. The subsequent extension of the region of current decay into the tail up to 150 R_E proceeded with a velocity of 1000 km/s, which exceeds the known velocities of such an extension by a factor of 5.Translated from Kosmicheskie Issledovaniya, Vol. 42, No. 6, 2004, pp. 622–631.Original Russian Text Copyright © 2004 by Solovyev, Moiseyev, Mullayarov, Du, Engebretson, Newitt.     

Radial gradient of low-energy (∼0.5–2 MeV) protons at 20–80 AU in the solar activity minima of the 22nd and 23rd cycles

The value of the radial gradient of low-energy (0.5–2 MeV) protons in the heliosphere at distances of 20–80 AU in the periods of solar activity minima in 1985–1987 and 1994–1997 was estimated using the data of the Voyager-1 and Voyager-2 spacecraft (s/c). Preliminary results on the dependence of the radial gradient on the distance were obtained for protons of these energies. The value of the radial gradient varies from –3% (AU)^–1 to –1% (AU)^–1 at distances from the Sun of 20–60 AU, reaching +0.7% (AU)^–1 at maximum considered distances (80 AU). The sign reversal of the proton radial gradient at a distance of 60–70 AU is interpreted as the appearance of a new component: up to the point of inversion there are mainly particles of the solar origin and/or accelerated in the inner heliosphere, while after the reversal of the gradients sign the fluxes of particles prevail whose source is located far from the Sun (maybe in the vicinity of the heliosphere boundary in the region of existence of the termination shock).Translated from Kosmicheskie Issledovaniya, Vol. 43, No. 1, 2005, pp. 3–8.Original Russian Text Copyright © 2005 by Logachev, Zeldovich, Surova.     

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.     

Theory of Singular Optimal Controls with Applications to Space Flight Mechanics

The problem of existence of solutions to the problems of optimum control with a singularity condition H(x, , u) const, u U is studied. The necessary conditions of appearance of the singular control, as well as time location and the numbers of singular sections of the trajectories of a controlled object under consideration, are obtained. The numerical algorithms of solving the problems with the condition of singularity of control are constructed. The widely known problem of optimum control of the motion of a high-altitude rocket probe is also considered, and it is proved that in this problem no intermediate section can exist in the law of variation of the thrust. The results of numerical experiments are presented together with their relative evaluation.     

A new method of determining the spectral index for protons in the energy range about 10<Superscript>13</Superscript> eV

We consider a relationship between the difference in spectral indices of the spectra of single hadrons and all hadrons (_sngl – _h) and the difference in the indices of the spectra of galactic cosmic ray (GCR) protons and nuclei. It is demonstrated that at the mountain level the ratio (_p – _Z)/(_sngl – _h) is always larger than unity, if (_sngl – _h) > 0.1. From the experimental value _sngl – _h = 0.4 ± 0.05 we derive that, in the vicinity of E = 10 TeV, _p – _Z 0.49 ± 0.06 , i.e., _p 3.09 ± 0.06.__________Translated from Kosmicheskie Issledovaniya, Vol. 43, No. 2, 2005, pp. 83–87.Original Russian Text Copyright © 2005 by Grigorov, Tolstaya.     

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.     

Families of Periodic Solutions to the Beletsky Equation

A review of achievements in the investigation of the planar periodic oscillations and rotations of a satellite around its center of mass, which moves along an elliptic orbit in the central gravitational field, is presented. These oscillations and rotations of the satellite are described by an ordinary differential equation (Beletsky equation) of the second order with periodic coefficients and two parameters. The equation is equivalent to a periodic Hamiltonian system with a single degree of freedom and has a singularity. It turned out that two-parameter families of the generalized periodic solutions to this simple equation form complicated structures of a new type. A comparison of numerous separate results made it possible to outline a sufficiently unified and complete picture of the location and structure of the families of generalized 2-periodic solutions with an integer number of rotation. These families are compared with known data related to the resonance rotations of the natural satellites in the solar system.     

Families of periodic solutions to Hamiltonian systems: Nonsymmetrical periodic solutions for a planar restricted three-body problem

A monodromy matrix calculated at a single arbitrary point of the periodic solution to a Hamiltonian system allows one to obtain both the direction of continuation for the family of solutions of the first (in Poincarés sense) kind and the multiplicity and direction of branching for periodic solutions of the second kind. In case of resonances 1 : 1 and 1 : 2 one needs to take into account the structure of elementary divisors of the monodromy matrix. Using the planar circular restricted three-body problem as an example, the infiniteness of the process of branching for a nonintegrable system and its finiteness for an integrable system are demonstrated. It is proved that periodic solutions of both first and second kinds which are obtained by continuation of symmetric periodic solutions of a restricted problem are also symmetric. The only exception is the case of resonance 1 : 1 and two second-order cells of the monodromy matrix in the Jordanian form. In this case, all periodic solutions of the second kind turned out to be nonsymmetrical. Examples of the families of nonsymmetrical periodic solutions are given.__________Translated from Kosmicheskie Issledovaniya, Vol. 43, No. 2, 2005, pp. 88–110.Original Russian Text Copyright © 2005 by Kreisman.     

Regularization of the Limit Problem of Satellite Librations in a Keplerian Orbit Taking Light Pressure into Account

The planar librations of a satellite whose center of mass moves along an elliptic orbit are considered. It is assumed that not only the gravitational moment but also the forces of light pressure act upon the satellite. Account is taken of the fact that the right-hand sides of the differential equations are nonanalytic functions of the phase variables. When e 1, e being the orbit's eccentricity, the deformations of solutions are considered for the case of a satellite moving along a highly elongated orbit. Such transformation of the initial system of differential equations is carried out so that the new system becomes regular up to the value e = 1. A limit problem corresponding to the case e = 1 is considered. When the azimuth angle of the light source coincides with the direction to the pericenter, the dynamical system is reversible. In this case, the known families of the periodic solutions to the problem can be continued up to the limit case.