Attitude equilibria of dual spin satellites subject to gravitational torques ofn bodies

We consider the motion of a dual spin satellite placed in the gravitational field of n material points, assuming that the satellite has no influence on the motion of these points. The main bodies are located at the libration points of the classical n bodies problem. We investigate the set of relative equilibria of the satellite. As in the elementary case of a gyrostat attracted by a single point, we show that this problem is equivalent to find the extremum of a quadratic function. We obtain all possible equilibria of the satellite by solving two algebraic equations. Sufficient conditions of stability of these relative equilibria are given.     

Stochastic wobble of accretion discs and jets from turbulent rocket torques

Models of accretion discs and their associated outflows often incorporate assumptions of axisymmetry and symmetry across the disc plane. However, for turbulent discs these symmetries only apply to averaged quantities and do not apply locally. The local asymmetries can induce local imbalances in outflow power across the disc mid-plane, which can in turn induce local tilting torques. Here we calculate the effect of the resulting stochastic torques on disc annuli that are a consequence of standard mean field accretion disc models. The torques induce a random walk of the vector perpendicular to the plane of each averaged annulus. This random walk is characterized by a radially dependent diffusion coefficient which we calculate for small angle tilt. We use the coefficient to calculate a radially dependent time-scale for annular tilt and associated jet wobble. The wobble time depends on the square of the wander angle so the age of a given system determines the maximum wobble angle. We apply this to examples of blazars, young stellar objects and binary engines of pre-planetary nebulae and microquasars. It is noteworthy that for an averaging time   t _w∼ 3 d  , we estimate a wobble angle for jets in SS 433 of  θ∼ 0.8°  , not inconsistent with observational data. In general the non-periodic nature of the stochastic wobble could distinguish it from faster periodic jet precession.     

Secular resonance,solar spin down,and the orbit of Mercury

A mechanism capable of accounting for the large mean eccentricity (0.175) and inclination (7°.2) of Mercury is discussed. Provided the gravitational field of the rapidly rotating primordial Sun had a sufficiently large second degree harmonic (i.e., J₂ ? order 10^?3), subsequent solar spin down would drive the orbit of Mercury through two secular resonances with Venus, one involving the precession of the line of apsides, the other one involving the regression of the nodal line. Resonance passage generates contributions to the eccentricity and inclination that are proportional to the square root of the characteristic solar spin down time. We find that an initial solar rotation l period of $\text{P ? 5}\text{1}\text{2}\text{hr}$ guarantees passage through resonance and that a spin down time of $\text{τ = Ω|dΩ/dt|}{}^{\mathrm{?1}}$ of order 10⁶ years could have produced the observed eccentricity and inclination. Such a primordial rotation rate is comparable to the measured rotations of very young stars and the spin down time appears consistent with the time scale derived for magnetic braking of the Sun's rotation by an intense solar wind during a T-Tauri stage of solar evolution.     

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.     

Stochastic Hill's equations for the study of errant rocket burns in orbit

The problem of errant rocket burns in low Earth orbit is of growing interest, especially in the area of safety analysis of nuclear powered spacecraft. The development of stochastic Hill's equations provides a rigorous mathematical tool for the study of such errant rocket maneuvers. These equations are analyzed within the context of a theory of linear dynamical systems driven by a random white noise. It is established that the trajectories of an errant rocket are realizations of a Gauss-Markov process, whose mean vector is given by the solution of a deterministic rocket problem. The time-dependent covariance matrix of the process is derived in an explicit form.     

The orbit of Cosmos 307 rocket and its use in atmospheric research

About 1500 observations from 46 observing stations were used to determine the orbit of Cosmos 307 rocket (1969-94B) at 25 epochs spread throughout its nine-month orbital life. The determination was made using the RAE computer program for the refinement of orbital parameters, PROP.     

The mutual potential and gravitational torques of two bodies to fourth order

The mutual gravitational potential of two bodies of arbitrary shape is expressed to fourth order in an extension of MacCullagh's Formula for a single body. The expressions for the gravitational torques acting on each body are derived in a form convenient for use in the differential equations describing the rotational dynamics.     

Lumped geopotential harmonics of order 29, from analysis of the orbit of Cosmos 837 rocket

The orbit of Cosmos 837 rocket (1976-62E) has been determined at 36 epochs between January and September 1978, using the RAE orbit refinement program PROP 6 with about 3000 observations. The inclination was 62.7° and the eccentricity 0.039. The orbital accuracy achieved was between 30m and 150m, both radial and crosstrack. The orbit was near 29:2 resonance in 1978 (exact resonance occurred on 14 May) and the values of orbital inclination obtained have been analysed to derive lumped 29th-order geopotential harmonic coefficients, namely:

$\text{10}{}^9\text{C}{}^{\mathrm{0,2}}{}_{29}\text{= ? 10 ± 15}$

and

$\text{10}{}^9\text{S}{}^{\mathrm{0,2}}{}_{29}\text{= ?76 ± 12}$

. These will be used in future, when enough results at different inclinations have accumulated, to determine individual coefficients of order 29. The values of lumped harmonics obtained from analysis of the values of eccentricity were not well defined, because of the high correlations between them and the errors in removing the very large perturbation (31 km) due to odd zonal harmonics.     

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 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.     

Analysis of the orbit of cosmos 268 rocket (1969-20B)

The orbit of Cosmos 268 rocket (1969-20B) has been determined at 28 epochs during its 342-day life, with the aid of the PROP5 orbit refinement program. All available observations were used, including 16 from the Hewitt camera at Malvern, 28 from the 200-mm camera at Meudon, 56 from the kinetheodolite at the Cape Observatory, 700 visual observations from volunteer observers, 500 US Navy observations and 200 British radar observations. The orbits are of very good accuracy for such a high-drag satellite, most of the values of inclination having standard deviations less than 0.002°. The most accurate orbits are those utilizing photographic observations, and the best of these has standard deviations of 0.00001 in eccentricity and 0.0001° in inclination.

The values of inclination obtained, after correction to allow for the effects of other perturbing forces, have been analysed to determine zonal wind speeds in the upper atmosphere at heights a little above perigee (230–250 km) averaged over latitudes up to about 25°. The results show a clear distinction between the wind at night (21 to 03 hr local time), which is west-to-east with an average speed of 140 ± 50 m/sec, and the wind by day (08 to 17 hr), which is east-to-west with an average speed of 110 ± 50 m/sec.     

Parametric study of near-wake structure of spherical and cylindrical bodies in the laboratory

Some aspects of the interaction between metal bodies and streaming rarefied plasmas were studied in a newly constructed Plasma Wind Tunnel as part of an attempt to investigate (via simulation) phenomena relevant to the spacecraft/space plasma interaction. Detailed near-wake ion current profiles for both spherical and cylindrical bodies at different body potentials (φ_S) and at different plasma flow parameters are presented. Various features of the profiles can be correlated, at least qualitatively, with both plasma and body characteristics. For example, the width of the wake zone appears proportional to the Debye length (λ_D) and depends on the potential of the target body although it appears to be relatively insensitive to the ratio $\text{S = V}{}_{\mathrm{<mtext>flow</mtext>}}\text{/(}\text{2kT}{}_{\mathrm{e}}\text{M}{}_{\mathrm{+}}\text{)}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}$. The amplitude of the ion current peak(s) also appears proportional to λ_D while, for fixed φ_S, the location of the peak is directly related to S and possibly dependent upon body geometry. The general importance of body geometry is qualitatively demonstrated. In addition, a discussion of the relevance of the above studies to previous in situ data obtained from the Ariel I and Gemini/Agena missions is given.     

On the spin motion of the rocket 1970–97 B

The spin period of COSMOS 378 rocket was studied during all its lifetime. An exponential increase of the spin period values is pointed out.     

Features of the upper atmosphere revealed by analysis of the orbit of cosmos 1009 rocket, 1978-50B

COSMOS 1009 rocket was launched on 19 May 1978 into an orbit with initial perigee height 150 km and apogee 1100 km: its lifetime was only 17 days. The orbit has been determined daily during the final 14 days of its life, using the RAE orbit refinement program PROP6,with about 1100 observations supplied by NORAD. An average accuracy of about 60 m, radial and cross-track, was achieved.The orbits were analysed to reveal three features of the upper atmosphere at heights between 125 and 175 km. From the decrease in perigee height, five values of density scale height, accurate to ±4%, were obtained. The first three were within 10% of those from CIRA 1972; the fourth, after a magnetic storm, was higher than expected; the fifth gave evidence of the decrease in drag coefficient at heights below 130 km.Atmospheric oblateness produced a change of 4° in perigee position during the last four days of the life. Analysis showed that the ellipticity of the upper atmosphere was approximately equal to that of the Earth, f, for the first two of the four days, and about $\text{1}\text{2}\text{f}$ in the last two.The orbital inclination decreased during the 14 days by about 50 times its standard deviation, and the observed variation was analysed to determine zonal winds at heights of 150–160 km at latitudes near 47° north. The zonal wind was very weak (0±30 m/s) for 23–28 May at local times near 03h; and 90±30 m/s east-to-west for 29 May to 4 June at local times near 01 h.     

Growth rate and decay of magnetospheric ring current

The rate of energy input to the ring current is studied as a function of solar wind parameters. The ring current dissipation rate is also examined. The decay constant τ in the main phase of a storm has been shown to be independent of its intensity and to equal (4 ± 2) h. In the recovery phase τ rises with increasing storm intensity.     

Determination and geophysical interpretation of the orbit of China 2 rocket (1971-18B)

The orbit of China 2 rocket, 1971-18B, has been determined at 114 epochs throughout its 5-yr life, using the RAE orbit refinement program PROP 6, with more than 7000 radar and optical observations from 83 stations.The rocket passed slowly enough through the resonances 14:1, 29:2, 15:1 and 31:2 to allow lumped geopotential harmonic coefficients to be calculated for each resonance, by least-squares fittings of theoretical curves to the perturbation-free values of inclination and eccentricity. These lumped coefficients can be combined with values from satellites at other inclinations, to obtain individual harmonic coefficients.The rotation rate of the upper atmosphere, at heights near 300 km, was estimated from the decrease in orbital inclination, and values of 1.15, 1.05, 1.10 and 1.05 rev/day were obtained between April 1971 and January 1976. From the variation in perigee height, 25 values of density scale height were calculated, from April 1971 to decay. Comparison with values from the COSPAR International Reference Atmosphere 1972 shows good agreement between April 1971 and October 1975, but the observational values are 10% lower, on average, than CIRA thereafter.A further 1400 observations, made during the final 15 days before decay, were used to determine 15 daily orbits. Analysis of these orbits reveals a very strong West-to-East wind, of 240 ± 40 ms^?1, at a mean height of 195 km under winter evening conditions, and gives daily values of density scale height in the last 7 days before decay.     

The determination and analysis of the orbit of NIMBUS 1 rocket, 1964-52B: Variations in orbital inclination

The influence of aerodynamic drag and the geopotential on the motion of the satellite 1964-52B is considered. A model of the atmosphere is adopted that allows for oblateness, and in which the density behaviour approximates to the observed diurnal variation. A differential equation governing the variation of the orbital inclination combining the effects of air drag with those of the Earth's gravitational field is given.The 310 observed values of inclination are modified by the removal of perturbations due to luni-solar attraction, solid Earth and ocean tides, solar radiation pressure, low-order long-periodic tesseral harmonic perturbations and changes due to precession. The method of removal of these effects is given in some detail.The variations in inclination due to drag are analysed to give four values of the average atmospheric rotation rate at heights of 296–476 km at latitude 0–54°. These values are as expected from previous analyses.The analysis of the change in inclination due to solar radiation pressure shows that this rapidly tumbling cylindrical satellite may be considered as equivalent to a spherical satellite of a given area-to-mass ratio.Analysis of the inclination near 15:1 resonance with the geopotential yields values of lumped geopotential harmonics of order 15 and 30, namely, $\text{10}{}^9\text{C}\text{?}{}^{0.1}{}_{15}\text{= ?31.2 ± 2.3 10}{}^9\text{S}\text{?}{}^{0.1}{}_{15}\text{= ?4.4 ± 3.2 10}{}^9\text{C}\text{?}{}^{0.2}{}_{30}\text{= 39.0 ± 10.7 10}{}^9\text{S}\text{?}{}^{0.2}{}_{30}\text{= 51.8 ± 10.0}$     

The effect of solar radiation pressure on the orbit of a cylindrical satellite

The effect of solar radiation pressure on the orbits of cylindrical satellites is considered. The cylinder is assumed to reflect radiation both specularly and diffusely. The resultant forces on a stationary cylindrical satellite are given. The evolution of the satellite's orbit is described for two particular modes of rotation. In both cases the satellites are assumed to be in circular Sun-synchronous orbits.