Mercury Ion Analyzer (MIA) onboard Mercury Magnetospheric Orbiter: MMO

The Mercury Magnetopsheric Orbiter (MMO) is one of the spacecraft of the BepiColombo mission; the mission is scheduled for launch in 2014 and plans to revisit Mercury with modern instrumentation. MMO is to elucidate the detailed plasma structure and dynamics around Mercury, one of the least-explored planets in our solar system. The Mercury Plasma Particle Experiment (MPPE) on board MMO is a comprehensive instrument package for plasma, high-energy particle, and energetic neutral particle atom measurements. The Mercury Ion Analyzer (MIA) is one of the plasma instruments of MPPE, and measures the three dimensional velocity distribution of low-energy ions (from 5 eV to 30 keV) by using a top-hat electrostatic analyzer for half a spin period (2 s). By combining both the mechanical and electrical sensitivity controls, MIA has a wide dynamic range of count rates for the proton flux expected around Mercury, which ranges from 10⁶ to 10¹² cm⁻² s⁻¹ str⁻¹ keV⁻¹, in the solar wind between 0.3 and 0.47 AU from the sun, and in both the hot and cold plasma sheet of Mercury’s magnetosphere. The geometrical factor of MIA is variable, ranging from 1.0 × 10⁻⁷ cm² str keV/keV for large fluxes of solar wind ions to 4.7 × 10⁻⁴ cm² str keV/keV for small fluxes of magnetospheric ions. The entrance grid used for the mechanical sensitivity control of incident ions also work to significantly reduce the contamination of solar UV radiation, whose intensity is about 10 times larger than that around Earth’s orbit.     

Spectral characteristics of electron fluxes at L < 2 under the Radiation Belts

The paper presents the analysis of experimental data on electron fluxes with energies 10 keV–10 MeV. Data were obtained during 1978–2005 years in different space experiments (COSMOS-900, MIR Space Station, ACTIVE, SAMPEX, CORONAS-I, CORONAS-F, NOAA POES-17, TATYANA and others). Two areas of electron flux enhancements are studied in the paper: the near-equatorial (L < 1.2) zone and the middle-latitude (1.2 < L < 1.9) zone. It is shown that electron flux enhancements are regularly registered at L < 2 and the observed formations have some typical features. Electron peaks at L < 1.2 appear sporadically while peaks at 1.2 < L < 1.9 are observed regularly. The approximations of spectra by several functions including kappa-function are presented.     

Optimal guidance based on receding horizon control for low-thrust transfer to libration point orbits

This paper addresses the design and computation of a guidance law for a transfer mission from an orbit near the Earth to a halo orbit around the libration point L₂ in the Sun–Earth system. The guidance law, which is designed based on receding horizon control and compensates for launch velocity errors that are introduced by inaccuracies of the launch vehicle, is solved using the generating function method. During the design of the closed-loop guidance law, the entire transfer mission, which is considered a nonlinear optimal control problem, is evaluated to obtain a nominal reference trajectory. Using the launch velocity errors and the uncertainty of the model, a spacecraft controlled by the proposed guidance law tracks the reference trajectory. Furthermore, the original Riccati differential equation in the receding horizon control algorithm is replaced by an equivalent convenient form of the Riccati differential equation that is based on the generating function. The high-efficiency solution of the equivalent equation avoids the online direct integration of the original Riccati differential equation, which significantly increases the computational efficiency for the receding horizon control problem. Numerical simulations using a nonlinear bicircular four-body model demonstrate the capabilities of the proposed receding horizon guidance law for the transfer mission. In addition, the generating function method improves the computational efficiency by at least one order of magnitude over the backward sweep method in solving the receding horizon control problem.     

Design and analysis of a direct transfer trajectory from a near rectilinear halo orbit to a low lunar orbit

A near rectilinear halo orbit (NRHO) is regarded as a potential orbit for a future deep space station that can effectively support sustainable crewed lunar exploration missions. In this paper, the L₂ southern NRHOs are selected as the research object, and a direct transfer trajectory from an NRHO to a low lunar orbit (LLO) is designed and analyzed. First, based on the circular restricted three-body problem, the characteristics of NRHOs are discussed including geometric behaviour, stability and synodic resonance. Second, optimal direct transfer trajectories are obtained by combining a local gradient optimization with a numerical continuation strategy. A reachable domain calculation model is established. Finally, simulations are carried out to verify the feasibility of the trajectory design method. The relationship between the velocity change and the reachable domain is further analysed through simulation calculations.     

A neutral gas mass spectrometer to measure the chemical composition of the stratosphere

The Polar Balloon Atmospheric Composition Experiment (P-BACE) is a new generation of neutral gas mass spectrometer based on the time-of-flight principle. P-BACE is the scientific experiment on the Mars Environment Analog Platform (MEAP) flown successfully on a balloon mission in summer 2008. The MEAP mission was flown with a 334,000 m³ helium balloon in the stratosphere on a semicircular trajectory from northern Sweden around the North Pole to Canada using the summer northern hemispheric wind current. The atmospheric conditions at an atmospheric altitude of 35–40 km are remarkably similar to those on the surface of Mars and thus the balloon mission was an ideal testbed for our mass spectrometer P-BACE. Originally this instrument was designed for in situ measurements of the chemical composition of the Martian atmosphere.P-BACE has a unique mass range from 0 to 1000 amu/q with a mass resolution m/Δm (FWHM) > 1000, and the dynamic range is at least six orders of magnitude. During this experiment, the acquisition of one mass spectrum is a sum of 65,535 single spectra, recorded in a time frame of 66 s.The balloon mission lasted 5 days and had successfully demonstrated the functionality of the P-BACE instrument during flight conditions. We had recorded more than 4500 mass spectra. With little modifications, P-BACE can be used on a planetary mission for Mars, but for example also for Venus or Mercury, if placed on a satellite.     

Propagation and origin of ultra high-energy cosmic rays

Propagation of UHE protons through CMB radiation leaves the imprint on energy spectrum in the form of Greisen–Zatsepin–Kuzmin (GZK) cutoff, bump (pile-up protons) and dip. The dip is a feature in energy range 1 × 10¹⁸–4 × 10¹⁹ eV, caused by electron–positron pair production on CMB photons. Calculated for power-law generation spectrum with index γ_g = 2.7, the shape of the dip is confirmed with high accuracy by data of Akeno-AGASA, HiRes, Yakutsk and Fly’s Eye detectors. The predicted shape of the dip is robust: it is valid for the rectilinear and diffusive propagation, for different discreteness in the source distribution, for local source overdensity, deficit, etc. This property of the dip allows us to use it for energy calibration of the detectors. The energy shift λ for each detector is determined by minimum χ² in comparison of observed and calculated dip. After this energy calibration the absolute fluxes, measured by AGASA, HiRes and Yakutsk detectors remarkably coincide in energy region 1 × 10¹⁸–1 × 10²⁰ eV. Below the characteristic energy E_c ≈ 1 × 10¹⁸ eV the spectrum of the dip flattens for both diffusive and rectilinear propagation, and more steep galactic spectrum becomes dominant at E < E_c. The energy of transition E_tr < E_c approximately coincides with the position of the second knee E_2kn, observed in the cosmic ray spectrum. The dip-induced transition from galactic to extragalactic cosmic rays at the second knee is compared with traditional model of transition at ankle, the feature observed at energy 1 × 10¹⁹ eV.     

ETALON spin period determination from kHz SLR data

Using kHz Satellite LASER Ranging (SLR) data of the SLR station Graz only, we determined the spin periods of the two ETALON satellites – launched into high orbits of about 20,000 km – and their spin period increase during 3 years. The determined spin period values and spin period increase rates at 2004-01-01 are: T_ET1 = 63 s + 0.484 s/year, and T_ET2 = 65.5 s + 0.401 s/year.     

Study of correlations between waves and particle fluxes measured on board the DEMETER satellite

The topic of relativistic electron dynamics in the outer radiation belt has received considerable attention for many years. Nevertheless, the problem of understanding the physical phenomenon involved is far from being resolved. In this paper, we use DEMETER observations to examine the variations of the energetic electron fluxes and ELF/VLF wave intensities in the inner magnetosphere during the intense 8 November 2004 magnetic storm. Electron flux spectra and associated wave intensity spectra are analysed throughout the magnetic storm and common characteristics or differences to other storm events are retained. The overall objective of this study is to identify and derive parameters that are relevant for particle flux modelling; the time constant characterizing the persistent decay after particle enhancement was found to be one of these important model parameters.The analysis of the 8 November 2004 event reveals that for L-shell parameter higher than 4, an electron flux dropout is observed during the storm’s main phase for electrons in the energy range 0.1–1 MeV, as has been reported from other measurements. Characteristic wave spectra accompanying this phase are analysed. They show a typical enhancement in the frequency range 0.3–10 kHz at onset for all L-shell values under consideration (2 < L < 5). During the first stage of the recovery phase, the electron fluxes are increased to a level higher than the pre-storm level, whereas the level of wave intensity in the frequency range observed below 300 Hz is at its highest. In the second stage, the particle flux decrease goes hand in hand with a global wave activity decline, the relaxation time of the latter being smaller than the former’s one. In some other cases, long-lasting electron enhancement associated with constant wave activity has been observed during this latter stage. For the above mentioned storm, while at low L values the decay time constants are higher for low energy electrons than for high energy electrons, this order is reversed at high L values. At about L = 3.6 the time constant is independent of electron energy.     

Low-energy near Earth asteroid capture using Earth flybys and aerobraking

Since the Sun-Earth libration points L₁ and L₂ are regarded as ideal locations for space science missions and candidate gateways for future crewed interplanetary missions, capturing near-Earth asteroids (NEAs) around the Sun-Earth L₁/L₂ points has generated significant interest. Therefore, this paper proposes the concept of coupling together a flyby of the Earth and then capturing small NEAs onto Sun–Earth L₁/L₂ periodic orbits. In this capture strategy, the Sun-Earth circular restricted three-body problem (CRTBP) is used to calculate target Lypaunov orbits and their invariant manifolds. A periapsis map is then employed to determine the required perigee of the Earth flyby. Moreover, depending on the perigee distance of the flyby, Earth flybys with and without aerobraking are investigated to design a transfer trajectory capturing a small NEA from its initial orbit to the stable manifolds associated with Sun-Earth L₁/L₂ periodic orbits. Finally, a global optimization is carried out, based on a detailed design procedure for NEA capture using an Earth flyby. Results show that the NEA capture strategies using an Earth flyby with and without aerobraking both have the potential to be of lower cost in terms of energy requirements than a direct NEA capture strategy without the Earth flyby. Moreover, NEA capture with an Earth flyby also has the potential for a shorter flight time compared to the NEA capture strategy without the Earth flyby.     

Asymmetric precipitation in a coronal loop as explanation of a singular observed spectrum

Almost 10 years of solar submillimeter observations have shown new aspects of solar activity, such as the presence of rapid solar spikes associated with the launch of coronal mass ejections and an increasing submillimeter spectral component in flares. We analyse the singular microwave–submillimeter spectrum of an M class solar flare on 20 December, 2002. Flux density observations measured by Sun patrol telescopes and the Solar Submillimeter Telescope are used to build the radio spectrum, which is fitted using Ramaty’s code. At submillimeter frequencies the spectrum shows a component different from the microwave classical burst. The fitting is achieved proposing two homogeneous sources of emission. This theoretical fitting is in agreement with differential precipitation through a magnetically asymmetric loop or set of loops. From a coronal magnetic field model we infer an asymmetric magnetic structure at the flare location. The model proposed to quantify the differential precipitation rates due to the asymmetry results in a total precipitation ratio Q₂/Q₁≈10⁴–10⁵, where Q₁(Q₂) represents the total precipitation in the loop foot with the high (low) magnetic field intensity. This ratio agrees with the electron total number ratio of the two sources proposed to fit the radio spectrum.     

Geomagnetic effects on mid-latitude railways: A statistical study of anomalies in the operation of signaling and train control equipment on the East-Siberian Railway

Intense geomagnetically induced currents (GIC) can hamper rail traffic by disturbing signaling and train control systems. GIC threats have been a concern for technological systems at high-latitude locations due to geomagnetic disturbances driven by substorm expansion electrojet or convection electrojet intensifications. However, other geomagnetic storm processes such as storm sudden commencement (SSC) and geomagnetic pulsations can also cause GIC concerns for technological systems. We present in this paper the first evidence based on statistical data for links between geomagnetic disturbances and faulty operations (anomalies) in the functioning of railway automatics and telemetry. We analyze anomalies of automatic signaling and train control equipment which occurred in 2004 on the East-Siberian Railway (corrected geomagnetic latitude _m = 46–51°N and longitude λ_m = 168–187°E). Our results reveal a seasonal effect in the number of anomalies per train similar to the one observed in geomagnetic activity (Kp, Ap, Dst indices). We also found an increase by a factor of 3 in the total duration of daily anomalies during intense geomagnetic storms (local geomagnetic index specific to Siberian Observatory A_max > 30), with a significant correlation between the daily sum of durations of anomalies with geomagnetic activity. Special attention was paid to failures not related to recognized technical malfunctions. We found that the probability of these failures occurring in geomagnetically disturbed periods was 5–7 times higher than the average anomaly occurrence.     

Single-event upset in geostationary transfer orbit during solar-activity maximum period measured by the Tsubasa satellite

This paper reports single-event upset (SEU) occurrence related to the space radiation environment in geostationary transfer orbit during solar-activity maximum period measured by the Tsubasa satellite. Most SEUs are measured in the inner radiation belt, indicating that they are mainly caused by trapped protons. Thus, the spatial distribution and the temporal variation of the SEU count correlate well with those of trapped protons. The peak SEU rate appears around L = 1.4. The transition point from SEUs caused by trapped protons to those caused by galactic cosmic rays is around L = 2.6. During the experiment period, increased SEU count was sometimes detected due to solar and geomagnetic events outside the inner radiation belt.     

Optical performance of PHEBUS/EUV detector onboard BepiColombo

BepiColombo, a mission of ESA (European Space Agency) in cooperation with JAXA (Japan Aerospace Exploration Agency), will explore Mercury, the planet closest to the Sun. BepiColombo will launch in 2014 on a journey lasting up to six and a half years; the data gathering phase should occupy a one year nominal mission, with a possible extension of another year. The data which will be brought back from the orbiters will tell us about the Hermean surface, atmospheric composition, and magnetospheric dynamics; it will also contribute to understanding the history and formation of terrestrial planets. The PHEBUS (Probing of Hermean Exosphere by Ultraviolet Spectroscopy) instrument will be flown on MPO: Mercury Planetary Orbiter, one of the two BepiColombo orbiters. The main purpose of the instrument is to reveal the composition and the distribution of the exosphere of Mercury through EUV (Extreme Ultraviolet: 55–155 nm) and FUV (Far Ultraviolet: 145–315 nm) measurements. A consortium composed of four main countries has been formed to build it. Japan provides the two detectors (EUV and FUV), Russia implements the scanning system, and France and Italy take charge of the overall design, assembly, test, integration, and also provide two small NUV (Near Ultraviolet) detectors (for the light from calcium and potassium molecules). An optical prototype of the EUV detector which is identical to the flight configuration has been manufactured and evaluated. In this paper, we show the first spectra results observed by the EUV channel optical prototype. We also describe the design of PHEBUS and discuss the possibility of detecting noble gases in Mercury’s exosphere taking the experimental results so far into account.     

Launch and commissioning of the PAMELA experiment on board the Resurs-DK1 satellite

PAMELA is a satellite borne experiment designed to study with great accuracy cosmic rays of galactic, solar, and trapped nature in a wide energy range (protons: 80 MeV–700 GeV, electrons 50 MeV–400 GeV). Main objective is the study of the antimatter component: antiprotons (80 MeV–190 GeV), positrons (50 MeV–270 GeV) and search for antimatter (with a precision of the order of 10⁻⁸). The experiment, housed on board the Russian Resurs-DK1 satellite, was launched on June, 15th 2006 in a 350 × 600 km orbit with an inclination of 70°. The detector consists of a permanent magnet spectrometer core to provide rigidity and charge sign information, a Time-of-Flight system for velocity and charge information, a silicon–tungsten calorimeter and a neutron detector for lepton/hadron identification. An anticounter system is used off-line to reject false triggers coming from the satellite. In self-trigger mode the calorimeter, the neutron detector and a shower tail catcher are capable of an independent measure of the lepton (e⁺ + e⁻) component up to 2 TeV. In this work we focus on the first months of operations of the experiment during the commissioning phase.     

The PAMELA experiment: A space-borne observatory for heliospheric phenomena

PAMELA is a multi-purpose apparatus composed of a series of scintillator counters arranged at the extremities of a permanent magnet spectrometer to provide charge, time-of-flight and rigidity information. Lepton/hadron identification is performed by a silicon–tungsten calorimeter and a Neutron detector placed at the bottom of the device. An Anticounter system is used offline to reject false triggers coming from the satellite. The device was put into orbit on June 15th 2006 in a pressurized container on board the Russian Resurs-DK1 satellite. The satellite is flying along a high inclination (70°), low Earth orbit (350–600 km), allowing to perform measurements in different points and conditions of the geomagnetosphere. PAMELA main goal is a precise measurement of the antimatter ( 80 MeV–190 GeV, e⁺ 50 MeV–270 GeV) and matter (p 80–700 GeV, e⁻ 50 MeV–400 GeV) component of the galactic cosmic rays. In this paper we focus on the capabilites of observations of heliospheric cosmic rays: trapped and semi-trapped particles in the proton and electron belts, solar particle events, Jovian electrons will be studied in the three years of expected mission.     

A linear stability analysis of large-Prandtl-number thermocapillary liquid bridges

A linear stability analysis is applied to determine the onset of oscillatory thermocapillary convection in cylindrical liquid bridges of large Prandtl numbers (4  Pr  50). We focus on the relationships between the critical Reynolds number Re_c, the azimuthal wave number m, the aspect ratio Γ and the Prandtl number Pr. A detailed Re_c–Pr stability diagram is given for liquid bridges with various Γ. In the region of Pr > 1, which has been less studied previously and where Re_c has been usually believed to decrease with the increase of Pr, we found Re_c exhibits an early increase for liquid bridges with Γ around one. From the computed surface temperature gradient, it is concluded that the boundary layers developed at both solid ends of liquid bridges strengthen the stability of basic axisymmetric thermocapillary convection at large Prandtl number, and that the stability property of the basic flow is determined by the “effective” part of liquid bridge.     

有限推力变轨的月球探测器发射轨道设计

月球探测器的发射是一个多体问题,也是严格的轨道交会问题。发射轨道方案从理论上有很多种形式,但在工程应用中,受到火箭能力、发射场和测控条件等多种因素制约,方案选择限制非常严格。利用现有火箭,采用有限推力从停泊轨道转移发射是最常用的一种发射方式。本文以某典型火箭为例,对采用有限推力变轨方案的奔月轨道设计流程、发射方案、轨道设计方法、精度分析、发射窗口设计进行了详细分析和研究。     

火星探测转移轨道初始设计与分析

根据推进方式和是否采用金星借力,火星转移轨道分为大推力直接转移轨道、大推力金星借力转移轨道、小推力直接转移轨道和小推力金星借力转移轨道4类。传统的轨道设计方法只是针对某一类特定的转移方案进行轨道优化,而并未针对不同的转移方案进行详细对比分析。文章以2020/2022年发射窗口为例,针对4类基本火星转移轨道进行研究。首先,基于不同轨道初始设计方法,对4类轨道进行了初始设计,得到了每类转移方案的能量最优转移轨道。然后,基于设计结果和能耗对4类转移方案进行了横向对比分析,得到了不同策略下的转移轨道的特性。基于小推力的火星探测任务轨道对发射能量要求低;大推力直接转移和借力金星的发射窗口交替分布,可以互为备份;基于小推力推进的探测器采用金星借力转移策略相比直接转移能够减少10%的能耗,优势十分明显。     

USV test flight by stratospheric balloon: Preliminary mission analysis

The Unmanned Space Vehicle test flights will use a 7 m 1300 kg aircraft. The first three launches will take place at the Italian Space Agency ASI base in Trapani–Milo, Sicily, through a stratospheric balloon that will drop the aircraft at a predefined height. After free fall acceleration to transonic velocities, the parachute deployment will allow a safe splash down in the central Mediterranean Sea. The goal of this article is to show the preliminary analysis results for the first USV flight.We carried out a statistical study for the year 2000–2003, evaluating the typical summer and winter launch windows of the Trapani–Milo base.First, in the center Mediterranean, we define safe recovery areas. They cannot be reached during the balloon ascending phase so, after a sufficiently long floating part able to catch the open sea, the balloon will go down to the release height (24 km). The simulation foresees a 400,000 m³ balloon and 3 valves for the altitude transfer.A safe splash down must occur far enough from the nearest coast: the minimum distance is considered around 25 km. The vehicle should be released at a distance, from the nearest coast, greater than this minimum amount plus the USV model maximum horizontal translation, during its own trajectory from balloon separation to splash down. In this way we define safe release areas for some possible translations.Winter stratospheric winds are less stable. The winter average flight duration is 7 h and it is probably too long for the diurnal recovery requirement and its scheduled procedures.Comparing past stratospheric balloons flights and trajectories computed using measured meteorological data (analysis), with their predictions made using forecast models and soundings, we obtain the standard deviation of the trajectory forecast uncertainty at the balloon–aircraft separation. Two cases are taken into account: predictions made 24 and 6 h before the launch.Assuming a Gaussian latitudinal uncertainty distribution for the prediction 6 h before the launch, we are able to identify the forecast trajectories that have a probability greater than 97% to reach the safe release areas.Simulating the summer windows trajectories for the years from 2000 to 2003 and for the favorable ground wind days, we obtain the number of trajectories with the desired forecast probabilities.     

长征运载火箭发射地火直接转移轨道研究

针对发射场射向范围和低温入轨级在停泊轨道上最长允许滑行时间严格受限下，工程设计火星探测任务发射轨道时，优化保障2~3周发射日期窗口的困难，以及精确设计拼接探测器分离点时刻轨道6根数耗时周期长等问题，提出了加入滑行时间限制再精确双向微分修正的设计算法（有别于传统B平面矢量法）。由地球影响球边界速度直接解出长征运载火箭入轨点轨道根数，将火箭飞行星下点弧段表示成滑行时间的分析式，快速找到满足诸多设计约束的初始发射轨道；并解决了考虑探测器小幅深空机动后火箭发射轨道的深入优化问题。精确力学模型下数值微分改进后，获得高精度计算结果，同等可比条件下运算精度不低于STK软件。实现了从地面起飞到抵达近火点的飞行轨道整体优化，确保了首次探火工程的顺利实施。