Satellite fault diagnosis using a bank of interacting Kalman filters

The main objective of this work is development and testing of a detection, isolation, and diagnosis algorithm based on interacting multiple model (IMM) filters for both partial (soft) and total (hard) reaction wheels faults in a spacecraft. This is shown to be accomplished under a number of different faulty mode scenarios for these actuators associated with the attitude control system (ACS) of a satellite. Various operating and faulty conditions due to changes and anomalies in the temperature, the power supply line voltage, and the loss of effectiveness of the torque and the current are considered in each reaction wheel associated with the three axes of the satellite. Once a fault mode is detected and isolated the recovery procedure can subsequently be engaged by invoking appropriate switching control strategies for the ACS. The application of a bank of interacting multiple Kalman filters for detection and diagnosis of anticipated reaction wheel failures in the ACS is described and developed. Compared with other model-based fault detection, diagnosis and isolation(FDDI) strategies developed in the control systems literature, our FDDI strategy is shown, through extensive numerical simulations, to be more accurate and robust with potential for extension to a number of other application areas.     

Magnetosphere Imaging Instrument (MIMI) on the Cassini Mission to Saturn/Titan

The magnetospheric imaging instrument (MIMI) is a neutral and charged particle detection system on the Cassini orbiter spacecraft designed to perform both global imaging and in-situ measurements to study the overall configuration and dynamics of Saturn’s magnetosphere and its interactions with the solar wind, Saturn’s atmosphere, Titan, and the icy satellites. The processes responsible for Saturn’s aurora will be investigated; a search will be performed for substorms at Saturn; and the origins of magnetospheric hot plasmas will be determined. Further, the Jovian magnetosphere and Io torus will be imaged during Jupiter flyby. The investigative approach is twofold. (1) Perform remote sensing of the magnetospheric energetic (E > 7 keV) ion plasmas by detecting and imaging charge-exchange neutrals, created when magnetospheric ions capture electrons from ambient neutral gas. Such escaping neutrals were detected by the Voyager l spacecraft outside Saturn’s magnetosphere and can be used like photons to form images of the emitting regions, as has been demonstrated at Earth. (2) Determine through in-situ measurements the 3-D particle distribution functions including ion composition and charge states (E > 3 keV/e). The combination of in-situ measurements with global images, together with analysis and interpretation techniques that include direct “forward modeling’’ and deconvolution by tomography, is expected to yield a global assessment of magnetospheric structure and dynamics, including (a) magnetospheric ring currents and hot plasma populations, (b) magnetic field distortions, (c) electric field configuration, (d) particle injection boundaries associated with magnetic storms and substorms, and (e) the connection of the magnetosphere to ionospheric altitudes. Titan and its torus will stand out in energetic neutral images throughout the Cassini orbit, and thus serve as a continuous remote probe of ion flux variations near 20R _S (e.g., magnetopause crossings and substorm plasma injections). The Titan exosphere and its cometary interaction with magnetospheric plasmas will be imaged in detail on each flyby. The three principal sensors of MIMI consists of an ion and neutral camera (INCA), a charge–energy–mass-spectrometer (CHEMS) essentially identical to our instrument flown on the ISTP/Geotail spacecraft, and the low energy magnetospheric measurements system (LEMMS), an advanced design of one of our sensors flown on the Galileo spacecraft. The INCA head is a large geometry factor (G ∼ 2.4 cm² sr) foil time-of-flight (TOF) camera that separately registers the incident direction of either energetic neutral atoms (ENA) or ion species (≥5^∘ full width half maximum) over the range 7 keV/nuc < E < 3 MeV/nuc. CHEMS uses electrostatic deflection, TOF, and energy measurement to determine ion energy, charge state, mass, and 3-D anisotropy in the range 3 ≤ E ≤ 220 keV/e with good (∼0.05 cm² sr) sensitivity. LEMMS is a two-ended telescope that measures ions in the range 0.03 ≤ E ≤ 18 MeV and electrons 0.015 ≤ E≤ 0.884 MeV in the forward direction (G ∼ 0.02 cm² sr), while high energy electrons (0.1–5 MeV) and ions (1.6–160 MeV) are measured from the back direction (G ∼ 0.4 cm² sr). The latter are relevant to inner magnetosphere studies of diffusion processes and satellite microsignatures as well as cosmic ray albedo neutron decay (CRAND). Our analyses of Voyager energetic neutral particle and Lyman-α measurements show that INCA will provide statistically significant global magnetospheric images from a distance of ∼60 R _S every 2–3 h (every ∼10 min from ∼20 R _S). Moreover, during Titan flybys, INCA will provide images of the interaction of the Titan exosphere with the Saturn magnetosphere every 1.5 min. Time resolution for charged particle measurements can be < 0.1 s, which is more than adequate for microsignature studies. Data obtained during Venus-2 flyby and Earth swingby in June and August 1999, respectively, and Jupiter flyby in December 2000 to January 2001 show that the instrument is performing well, has made important and heretofore unobtainable measurements in interplanetary space at Jupiter, and will likely obtain high-quality data throughout each orbit of the Cassini mission at Saturn. Sample data from each of the three sensors during the August 18 Earth swingby are shown, including the first ENA image of part of the ring current obtained by an instrument specifically designed for this purpose. Similarily, measurements in cis-Jovian space include the first detailed charge state determination of Iogenic ions and several ENA images of that planet’s magnetosphere.This revised version was published online in July 2005 with a corrected cover date.     

Orbits and Integrals in Self-Consistent Systems

We find the forms of the orbits in a self-consistent galactic model generated by a N-body simulation of the collapse of a protogalaxy. The model represents a stationary elliptical galaxy of type E5, which is approximately axisymmetric around its longest axis. The orbits are of three main types, box orbits (including box-like orbits), tube orbits and chaotic orbits. The box or box-like and tube orbits are represented by closed invariant curves on a Poincaré surface of section. The forms of the orbits and of the invariant curves can be explained by a third integral of motion I, that is given by the Giorgilli (1979) computer program. The nonresonant form of the third integral explains the box orbits, while a resonant form of this integral explains both the box orbits and the 1:1 tube orbits. The N-body model gives the distribution of velocities F, which is an exponential of the third integral.     

Spectral-domain covariance estimation with a priori knowledge

A knowledge-aided spectral-domain approach to estimating the interference covariance matrix used in space-time adaptive processing (STAP) is proposed. Prior knowledge of the range-Doppler clutter scene is used to identify geographic regions with homogeneous scattering statistics. Then, minimum-variance spectral estimation is used to arrive at a spectral-domain clutter estimate. Finally, space-time steering vectors are used to transform the spectral-domain estimate into a data-domain estimate of the clutter covariance matrix. The proposed technique is compared with ideal performance and to the fast maximum likelihood technique using simulated results. An investigation of the performance degradation that can occur due to various inaccurate knowledge assumptions is also presented     

Photoelectron Counting Statistics and Error Performance of m-Distributed Optical Signal

By the term "m-distributed optical signal" we mean a noise-like optical signal whose envelope (or intensity) fluctuation probability is modeled by Nakagami's "m-distribution." Using the m-distribution which has been widely used as an analytical model of the fading envelope in radio communications, it is shown that one can generally analyze the statistical properties such as the photoelectron count probabilities and error probabilities for the wider class of noise-like optical signals; some numerical results are also given.     

A relative navigation system for formation flight

A relative navigation system based on both the Inertial Navigation System (INS) and the Global Positioning System (GPS) is developed to support situational awareness during formation flight. The architecture of the system requires an INS/GPS integration across two aircraft via a data link. A fault-tolerant federated filter is used to estimate the relative INS errors based on relative GPS measurements and a range measurement obtained from the data link. The filter is constructed based on a reduced-order model of the relative INS error process. A method for analyzing the filter performance is presented. A case involving two helicopters in formation flight is studied under three different night trajectories to account for the effect of vehicle motion on the INS state transition matrix. The results of the covariance analysis are compared with actual night results over an instrumented test range.     

Instrument Failure Detection in Partially Observable Systems

Instrument failure detection using the dedicated observer scheme (DOS) depends on partial state observability through each instrument which is monitored. For instrument fault detection by the DOS technique, a quantitative measure of partial state observability is established for each instrument and used to determine a necessary condition on the output structure of the system. This measure, called the internal redundancy of the instrument, indicates the complexity of the logic required for failure detection, and it also indicates where some hardware redundancy can be introduced into the system to improve the fault detection capability of the DOS. The principles developed are applied to a simulation of the pitch axis autopilot of the A7 jet aircraft.     

Pure proportional navigation against time-varying target manoeuvres

Capturability of the pure proportional navigation (PPN) guidance law against a target executing bounded piecewise continuous time-varying manoeuvres is investigated. A qualitative analysis is carried out to obtain a set of sufficient conditions for capture defined on the engagement parameters and initial conditions. These conditions are significantly less restrictive than the ones obtained previously by others using the Lyapunov method. It is shown that the actual capture region for time-varying target manoeuvres, obtained by using the conditions derived, is much larger than that obtained from the Lyapunov technique. We also show that though a bounded tine varying target manoeuvre does change the constant target manoeuvre capture region to some extent, it does not reduce it drastically. Further, we show that the worst case capture region is obtained when the target executes a constant manoeuvre equal to the bound on the manoeuvre level. Some bounds on the missile lateral acceleration are also obtained for certain regions in the engagement plane. These results are generalizations and extensions of existing results on the capturability of the PPN guidance law against targets executing constant or time-varying target manoeuvres     

The plasma mantle: Composition and other characteristics observed by means of the PROGNOZ-7 satellite

A brief summary of the main results of magnetospheric ion composition measurements in general is first presented. PROGNOZ-7 measurements in the nightside plasma mantle are then described and analyzed. Some of the results are the following: In the nightside mantle not too far from midnight the properties of the mantle are sometimes consistent with the open magnetosphere model. However during most magnetic storm situations O⁺ ions appear in the mantle in large proportions and with high energies. The acceleration process affecting the ions has been found in several cases to give equal amounts of energy to all ions independent of mass. Along the flanks of the magnetosphere the flow of the plasma is often low or absent. The O⁺ content is high (up to 20%) and the energy spectrum of both ions and electrons may be very hot, even up to the level of the ring current plasma in the keV range.The O⁺ content in the plasma mantle is positively correlated with the magnetospheric activity level. The mantle, however, does not appear to be the dominating source for the storm time ring current. Direct acceleration of ionospheric ions onto the closed field lines of the plasma sheet and ring current is most likely the main source. The magnetopause on the nightside and along the flanks of the magnetosphere appears to be a fairly solid boundary for mantle ions of ionospheric origin. This is especially evident during periods with high geomagnetic activity, when the mantle is associated with fairly strong fluxes of O⁺ ions.An interesting observation in most of the mantle passages during geomagnetically disturbed periods is the occurrence of intense, magnetosheath like, regions deep inside the mantle. In some cases these regions with strong antisunward flow and with predominant magnetosheath ion composition was observed in the innermost part of the mantle, i.e. marking a boundary region between the lobe and the mantle. These magnetosheath penetration events are usually associated with strong fluxes of accelerated ionospheric ions in nearby parts of the mantle. Evanescent penetration regions with much reduced flow properties are frequently observed in the flank mantle.     

Asymptotically Stable Motions of a Gyro in a Spinning Vehicle

An analysis of the motion of a single-axis rate gyroscope mounted in a space vehicle which is spinning with uncertain constant angular velocity xz about the spin axis of the gyro is presented. The nonlinearity in the equation of motion of the gimbal is retained. Using the Lyapunov approach, regions are obtained in the k-xz plane where k is the torsional spring constant, for asymptotic stability and for global asymptotic stability, and an estimate of the region of stability in state space is presented. Analytical relations for the selection of the gyro parameters are derived.