Practical stabilization in attitude thruster control

The attitude stabilization of a spacecraft using thrusters is considered from a practical point of view, i.e., when actuator constraints, uncertainties and failures, measurement noise, fuel consumption, and inertia matrix uncertainty are considered. A variable structure controller based on sliding modes is obtained first, which guarantees global exponential stability under actuator constraints, uncertainties, and failures. The design is based on a continuous average-proportional torque selection function and a set of feasible sliding surfaces, which may be adapted to improve performance. Exponential convergence towards ultimate bounds is guaranteed when all other practical issues arise.     

Possible shock wave in the local interstellar plasma, very close to the heliosphere

We show, using the HST — GHRS data on velocity and temperature in the nearby interstellar medium, that the observed 3 – 4 km s^–1 relative velocity between the Local Interstellar Cloud (LIC) and the so-called G-cloud located in the Galactic Center hemisphere can be quite naturally explained assuming that the two clouds do interact with each other. In the proposed interpretation the two media are separated by a (quasiperpendicular) MHD shock front propagating from the LIC into the G-cloud. The LIC plasma is then nothing else but the shocked (compression 1.3 – 1.4) gas of the G-cloud. A 1-D single-fluid solution of the Rankine — Hugoniot equations can fit the most probable observed values of the relative velocity (3.75 km/s), LIC (6700 K) and G-cloud (5400 K) kinetic temperatures, if the plasma-beta of the LIC plasma is in the range 1.3 – 1.5 (Table 1). This corresponds to a super — fast magnetosonic motion of the heliosphere through the LIC, independently of LIC density. The LIC magnetic field strength is 1.9 (3.1) G for the LIC electron density n_e = 0.04 (0.10) cm^–3. In this case the shock is less than 30 000 AU away and moves at about 10 km s^–1 relative to the LIC plasma. The Sun is chasing the shock and should catch up with it in about 10⁴ years. If the heliospheric VLP emissions cutoff at 1.8 kHz is indicative of n_e (LIC) = 0.04 cm^–3 (Gurnett et al., 1993), the (pure plasma) bowshock ahead of the heliopause could be the source of quasi-continuous heliospheric 2-kHz emission band. We believe that with the expected increase in the performance of modern spectroscopic instrumentation the proposed method of magnetic field evaluation may in the future find wider application in the studies of the interstellar medium.     

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     

ExoMars 2016 Schiaparelli Module Trajectory and Atmospheric Profiles Reconstruction

On 19^th October 2016 Schiaparelli module of the ExoMars 2016 mission flew through the Mars atmosphere. After successful entry and descent under parachute, the module failed the last part of the descent and crashed on the Mars surface. Nevertheless the data transmitted in real-time by Schiaparelli during the entry and descent, together with the entry state vector as initial condition, have been used to reconstruct both the trajectory and the profiles of atmospheric density, pressure and temperature along the traversed path.The available data-set is only a small sub-set of the whole data acquired by Schiaparelli, with a limited data rate (8 kbps) and a large gap during the entry because of the plasma blackout on the communications.This paper presents the work done by the AMELIA (Atmospheric Mars Entry and Landing Investigations and Analysis) team in the exploitation of the available inertial and radar data. First a reference trajectory is derived by direct integration of the inertial measurements and a strategy to overcome the entry data gap is proposed. First-order covariance analysis is used to estimate the uncertainties on all the derived parameters. Then a refined trajectory is computed incorporating the measurements provided by the on-board radar altimeter.The derived trajectory is consistent with the events reported in the telemetry and also with the impact point identified on the high-resolution images of the landing site.Finally, atmospheric profiles are computed tacking into account the aerodynamic properties of the module. Derived profiles result in good agreement with both atmospheric models and available remote sensing observations.     

Definition of Wingtip Devices Shape and Dimensions for a Subsonic Airliner

The results of flow field numerical simulation on the typical wing-body prototype of the modern DLR-F4 airliner under sub- and transonic compressible air flow are presented. Using the DLR-F4 CAD model, the effect of the wingtip end plate area and of the cant angle of a typical Whitcomb winglet is studied. The dependencies of the model lift-to-drag ratio increment on the flat wingtip end plate relative area and on the cant angle of an airfoil Whitcomb winglet are obtained. The concept of an elliptic winglet with a variable cant angle that similar to the winglet used on Airbus A350 is studied. A technique is developed for solving the multi-parameter design optimization task for the Whitcomb winglet, taking the maximum lift-to-drag ratio of the wing as a criterion for optimization.     

Efficiency Estimation of Electrothermal Thrusters Use in the Control System of the Technological Spacecraft Motion

The influence of various control systems of the orbital motion of a technological spacecraft on the level of microacceleration of its internal environment is simulated. Conclusions are drawn about the effectiveness of control systems with different actuators for realization of certain gravitationally sensitive processes onboard a spacecraft.     

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.     

Isotopic Signatures of Presolar Materials in Interplanetary Dust

Interplanetary dust particles collected in the stratosphere frequently exhibit enrichments in deuterium (D) and ¹⁵N relative to terrestrial materials. These effects are most likely due to the preservation of presolar interstellar materials. While the elevated D/H ratios probably resulted from mass fractionation during chemical reactions at very low < 100 K temperatures, the origin of the N isotopic anomalies remains unresolved. The bulk of the N-bearing material may have obtained its isotopic signatures from low temperature chemistry, but a nucleosynthetic origin is also possible. This revised version was published online in August 2006 with corrections to the Cover Date.     

A Class of TVD Type Combined Numerical Scheme for MHD Equations With a Survey About Numerical Methods in Solar Wind Simulations

It has been believed that three-dimensional, numerical, magnetohydrodynamic (MHD) modelling must play a crucial role in a seamless forecasting system. This system refers to space weather originating on the sun; propagation of disturbances through the solar wind and interplanetary magnetic field (IMF), and thence, transmission into the magnetosphere, ionosphere, and thermosphere. This role comes as no surprise to numerical modelers that participate in the numerical modelling of atmospheric environments as well as the meteorological conditions at Earth. Space scientists have paid great attention to operational numerical space weather prediction models. To this purpose practical progress has been made in the past years. Here first is reviewed the progress of the numerical methods in solar wind modelling. Then, based on our discussion, a new numerical scheme of total variation diminishing (TVD) type for magnetohydrodynamic equations in spherical coordinates is proposed by taking into account convergence, stability and resolution. This new MHD model is established by solving the fluid equations of MHD system with a modified Lax-Friedrichs scheme and the magnetic induction equations with MacCormack II scheme for the purpose of developing a combined scheme of quick convergence as well as of TVD property. To verify the validation of the scheme, the propagation of one-dimensional MHD fast and slow shock problem is discussed with the numerical results conforming to the existing results obtained by the piece-wise parabolic method (PPM). Finally, some conclusions are made. This revised version was published online in August 2006 with corrections to the Cover Date.