MIDACO on MINLP space applications

A numerical study on two challenging mixed-integer non-linear programming (MINLP) space applications and their optimization with MIDACO, a recently developed general purpose optimization software, is presented. These applications are the optimal control of the ascent of a multiple-stage space launch vehicle and the space mission trajectory design from Earth to Jupiter using multiple gravity assists. Additionally, an NLP aerospace application, the optimal control of an F8 aircraft manoeuvre, is discussed and solved. In order to enhance the optimization performance of MIDACO a hybridization technique, coupling MIDACO with an SQP algorithm, is presented for two of these three applications. The numerical results show, that the applications can be solved to their best known solution (or even new best solution) in a reasonable time by the considered approach. Since using the concept of MINLP is still a novelty in the field of (aero)space engineering, the demonstrated capabilities are seen as very promising.     

Estimates of vertical land motion along the southwestern coasts of Turkey from coastal altimetry and tide gauge data

The differences between coastal altimetry and sea level time series of tide gauges in between March 1993 and December 2009 are used to estimate the rates of vertical land motion at three tide gauge locations along the southwestern coasts of Turkey. The CTOH/LEGOS along-track coastal altimetry retrieves altimetric sea level anomalies closer to the coast than the standard along-track altimetry products. However, the use of altimetry very close to the coast is not found to improve the results. On the contrary, the gridded and interpolated AVISO merged product exhibits the best agreement with tide gauge data as it provides the smoothest variability both in space and time compared with along track altimetry data. The Antalya gauge to the south (in the Mediterranean Sea) and the Mentes/Izmir gauge to the west (in the Aegean Sea) both show subsidence while the Bodrum tide gauge to the south (in the Aegean Sea) shows no significant vertical land motion. The results are compared and assessed with three independent geophysical vertical land motion estimates like from GPS. The GIA effect in the region is negligible. The VLM estimates from altimetry and tide gauge data are in good agreement both with GPS derived vertical velocity estimates and those inferred from geological and archaeological investigations.     

Issues of geologically-focused situational awareness in robotic planetary missions: Lessons from an analogue mission at Mistastin Lake impact structure,Labrador, Canada

Remote robotic data provides different information than that obtained from immersion in the field. This significantly affects the geological situational awareness experienced by members of a mission control science team. In order to optimize science return from planetary robotic missions, these limitations must be understood and their effects mitigated to fully leverage the field experience of scientists at mission control.     

Testing the three-dimensional IRI-SIRMUP-P mapping of the ionosphere for disturbed periods

This paper describes the three-dimensional (3-D) electron density mapping of the ionosphere given as output by the assimilative IRI-SIRMUP-P (ISP) model for three different geomagnetic storms. Results of the 3-D model are shown by comparing the electron density profiles given by the model with the ones measured at two testing ionospheric stations: Roquetes (40.8°N, 0.5°E), Spain, and San Vito (40.6°N, 17.8°E), Italy. The reference ionospheric stations from which the autoscaled foF2 and M(3000)F2 data as well as the real-time vertical electron density profiles are assimilated by the ISP model are those of El Arenosillo (37.1°N, 353.3°E), Spain, Rome (41.8°N, 12.5°E), and Gibilmanna (37.9°N, 14.0°E), Italy. Overall, the representation of the ionosphere made by the ISP model is better than the climatological representation made by only the IRI-URSI and the IRI-CCIR models. However, there are few cases for which the assimilation of the autoscaled data from the reference stations causes either a strong underestimation or a strong overestimation of the real conditions of the ionosphere, which is in these cases better represented by only the IRI-URSI model. This ISP misrepresentation is mainly due to the fact that the reference ionospheric stations covering the region mapped by the model turn out to be few, especially for disturbed periods when the ionosphere is very variable both in time and in space and hence a larger number of stations would be required. The inclusion of new additional reference ionospheric stations could surely smooth out this concern.     

Performance evaluation of IRI-2007 at equatorial latitudes and its Matlab version for GNSS applications

International Reference Ionosphere (IRI) model is the widely used empirical model for ionospheric predictions, especially TEC which is an important parameter for radio navigation and communication. The Fortran based IRI-2007 does not support real-time interactive visualization and debugging. Therefore, the source code is converted into Matlab and is validated for the purposes of this study. This facilitates easy representation of results and for near real-time implementation of IRI in the applications including spacecraft launching, now casting, pseudolite based navigation systems etc. In addition, the vertical delay results over the equatorial region derived from IRI and GPS data of three IGS stations namely Libreville (Garbon, Africa), Brasilia (Brazil, South America) and Hyderabad (India, Asia) are compared. As the IRI model does not account for plasmasphere TEC, the vertical delays are underestimated compared to vertical delays of GPS signals. Therefore, the model should be modified accordingly for precise TEC estimation.     

Active debris removal of multiple priority targets

Today’s space debris environment shows major concentrations of objects within distinct orbital regions for nearly all size regimes. The most critical region is found at orbital altitudes near 800 km with high declinations. Within this region many satellites are operated in so called sun-synchronous orbits (SSO). Among those, there are Earth observation, communication and weather satellites. Due to the orbital geometry in SSO, head-on encounters with relative velocities of about 15 km/s are most probable and would thus result in highly energetic collisions, which are often referred to as catastrophic collisions, leading to the complete fragmentation of the participating objects. So called feedback collisions can then be triggered by the newly generated fragments, thus leading to a further population increase in the affected orbital region. This effect is known as the Kessler syndrome.     

Evaluation of a combined electrostatic and magnetostatic configuration for active space-radiation shielding

Developing successful and optimal solutions to mitigating the hazards of severe space radiation in deep space long duration missions is critical for the success of deep-space explorations. A recent report (Tripathi et al., 2008) had explored the feasibility of using electrostatic shielding. Here, we continue to extend the electrostatic shielding strategy and examine a hybrid configuration that utilizes both electrostatic and magnetostatic fields. The main advantages of this system are shown to be: (i) a much better shielding and repulsion of incident ions from both solar particle events (SPE) and galactic cosmic rays (GCR), (ii) reductions in the power requirement for re-charging the electrostatic sub-system, and (iii) low requirements of the magnetic fields that are well below the thresholds set for health and safety for long-term exposures. Furthermore, our results show transmission levels reduced to levels as low as 30% for energies around 1000 MeV, and near total elimination of SPE radiation by these hybrid configurations. It is also shown that the power needed to replenish the electrostatic charges due to particle hits from the GCR and SPE radiation is minimal.     

Height of shock formation in the solar corona inferred from observations of type II radio bursts and coronal mass ejections

Employing coronagraphic and EUV observations close to the solar surface made by the Solar Terrestrial Relations Observatory (STEREO) mission, we determined the heliocentric distance of coronal mass ejections (CMEs) at the starting time of associated metric type II bursts. We used the wave diameter and leading edge methods and measured the CME heights for a set of 32 metric type II bursts from solar cycle 24. We minimized the projection effects by making the measurements from a view that is roughly orthogonal to the direction of the ejection. We also chose image frames close to the onset times of the type II bursts, so no extrapolation was necessary. We found that the CMEs were located in the heliocentric distance range from 1.20 to 1.93 solar radii (Rs), with mean and median values of 1.43 and 1.38 Rs, respectively. We conclusively find that the shock formation can occur at heights substantially below 1.5 Rs. In a few cases, the CME height at type II onset was close to 2 Rs. In these cases, the starting frequency of the type II bursts was very low, in the range 25–40 MHz, which confirms that the shock can also form at larger heights. The starting frequencies of metric type II bursts have a weak correlation with the measured CME/shock heights and are consistent with the rapid decline of density with height in the inner corona.     

Variations in statistical parameters of the NmF2 equinoctial asymmetry with latitude and solar activity near noon

Hourly values of the F2-layer peak density, NmF2, measured by 95 ionosondes near noon from 1957 to 2011 at low and middle geomagnetic latitudes of the northern and southern geographic hemispheres are used in a statistical study of the NmF2 equinoctial asymmetry. The ratios, R, of NmF2 measured during 61 days around the March equinox to NmF2 measured during 61 days around the September equinox at the same UT near noon during geomagnetically quiet daytime conditions for approximately the same solar activity conditions over the same ionosonde are analyzed. The conditional probability of the occurrence of R in an interval of R, the most probable value of R, and the mean expected value of R are calculated for the first time for the low, moderate, and high solar activity levels to study variations in these statistical parameters with latitude and solar activity. These statistical parameters are averaged over 5° geomagnetic latitude interval in the northern and southern geographic hemispheres to calculate and to study for the first time trends in latitude and solar activity of these averaged NmF2 equinoctial asymmetry statistical characteristics.     

Numerical modeling of inhomogeneous orographic wave influence on planetary waves in the middle atmosphere

Parameterization of dynamical and thermal effects of stationary orographic gravity waves (OGWs) generated by the Earth’s surface topography is incorporated into a numerical model of general circulation of the middle and upper atmosphere. Responses of atmospheric general circulation and characteristics of planetary waves at altitudes from the troposphere up to the thermosphere to the effects of OGWs propagating from the earth surface are studied. Changes in atmospheric circulation and amplitudes of planetary waves due to variations of OGW generation and propagation in different seasons are considered. It is shown that during solstices the main OGW dynamical and heat effects occur in the middle atmosphere of winter hemispheres, where changes in planetary wave amplitudes due to OGWs may reach up to 50%. During equinoxes OGW effects are distributed more homogeneously between northern and southern hemispheres.