Latitudinal distribution of solar coronal active regions

We studied the cyclic evolution of the latitudinal distribution of solar coronal active regions based on daily images from SOHO EIT for the period 1995–2017. Fully automated software was used, which included the following steps: initial preparation of images in the data series, normalization of histograms and correction of limb brightening, segmentation of images using threshold intensity values obtained from their histograms, scanning of segmented images in heliographic coordinates and obtaining profiles of latitudinal distribution of coronal active regions for each image of the data series. From the output data, we obtained a temporary change in the latitudinal distribution profiles and the migration of activity centers on the solar disk. From the period of minimum activity to the next minimum in both hemispheres, activity centers begin to migrate from high latitudes towards the equator. At the same time, the general center of activity repeatedly changes the direction of migration. The latitudinal distribution of the so-called presence factor of coronal active regions closely resembles the magnetic butterfly diagram, which proves their direct causal relationships. Variations in the presence factor of coronal active regions are correlated with cyclic variations in the sunspot daily numbers.     

Plasma polarization electric field derived from radio sounding of solar wind acceleration region with spacecraft signals

Presented is the analytical approximation of averaged solar wind velocity radial dependence in the solar wind acceleration region at heliolatitudes below 60° under low and moderate solar activity. This empirical approximation is based on the data of radio sounding of the solar corona with radio signals from various spacecraft. Deduced is an equation connecting the solar wind velocity radial dependence and the radial dependence of solar wind plasma polarization electric field intensity. This allows constructing a semi-empirical radial dependence of plasma polarization electric field corresponding to the empirical radial dependence of solar wind velocity. Main properties of the semi-empirical dependence, which is based on radio sounding data, are described.     

Equatorial anomaly according to the Interkosmos-19 data and IRI model: A comparison

The comparison of the IRI model with the foF2 distribution in the equatorial anomaly region obtained by topside sounding onboard the Interkosmos-19 satellite has been carried out. The global distribution of foF2 in terms of LT-maps was constructed by averaging Intercosmos-19 data for summer, winter, and equinox. These maps, in fact, represent an empirical model of the equatorial anomaly for high solar activity F10.7 ~ 200. The comparison is carried out for the latitudinal foF2 profiles in the characteristic longitudinal sectors of 30, 90, 210, 270, and 330°, as well as for the longitudinal variations in foF2 over the equator. The largest difference between the models (up to 60%) for any season was found in the Pacific longitudinal sector of 210°, where there are a few ground-based sounding stations. Considerable discrepancies, however, are sometimes observed in the longitudinal sectors, where there are many ground-based stations, for example, in the European or Indian sector. The discrepancies reach their maximum at 00 LT, since a decay of the equatorial anomaly begins before midnight in the IRI model and after midnight according to the Interkosmos-19 data. The discrepancies are also large in the morning at 06 LT, since in the IRI model, the foF2 growth begins long before sunrise. In the longitudinal variations in foF2 over the equator at noon, according to the satellite data, four harmonics are distinguished in the June solstice and at the equinox, and three harmonics in the December solstice, while in the IRI model only two and one harmonics respectively are revealed. In diurnal variations in foF2 and, accordingly, in the equatorial anomaly intensity, the IRI model does not adequately reproduce even the main, evening extremum.     

Solar sailing in realistic mode,based on a locally-optimal control laws

Ballistic design of solar sailing missions in the solar system is composed of defining the design parameters, the control programs, and the trajectories that provide performance goals of a flight. The use of a solar sail spacecraft imposes specific restrictions on mission parameters that include the degradation limit on the flight duration, the maximum temperature of solar sail's surface, the minimum distance from the Sun, the maximum angular velocity of the spacecraft's rotation and others.Many authors considered the impact of these restrictions on the design of the mission separately, but they used a sophisticated method of finding the exact optimal motion control or applied the most straightforward laws of motion control. This paper uses local-optimal control laws at the complete mathematical models of motion and functioning of solar sail spacecraft to describe a technique of designing interplanetary missions. The described method avoids the need to obtain an accurate optimal solution to the control problem and does not cause significant computational difficulties.     

Solar sails for perturbation relief: Application to asteroids

A major cause of spacecraft orbital variation comes from natural perturbations, which, in close proximity of a body, are dominated by its non-spherical nature. For small bodies, such as asteroids, these effects can be considerable, given their uneven (and uncertain) mass distribution. Solar sail technology is proposed to reduce or eliminate the net secular effects of the irregular gravity field on the orbit. Initially, a sensitivity analysis will be carried out on the system which will show high sensitivity to changes in initial conditions. This presents a challenge for optimisation methods which require an initial guess of the solution. As such, the Genetic Algorithm (GA) is proposed as the preferred optimisation method as this requires no initial guess from the user. A multi-objective optimisation is performed which aims to achieve a periodic orbit whilst also minimising the effort required by the sail to do so. Given the system sensitivity, the control law for one orbit is not necessarily applicable for any subsequent orbit. Therefore, a new method of updating the control law for subsequent orbits is presented, based on linearisation and use of a Control Transition Matrix (CTM). The techniques will later find application in a multiple asteroid rendezvous mission with a solar sail as the primary propulsion system.     

Orbit determination and prediction of GEO satellite of BeiDou during repositioning maneuver

In order to establish a continuous GEO satellite orbit during repositioning maneuvers, a suitable maneuver force model has been established associated with an optimal orbit determination method and strategy. A continuous increasing acceleration is established by constructing a constant force that is equivalent to the pulse force, with the mass of the satellite decreasing throughout maneuver. This acceleration can be added to other accelerations, such as solar radiation, to obtain the continuous acceleration of the satellite. The orbit determination method and strategy are illuminated, with subsequent assessment of the orbit being determined and predicted accordingly. The orbit of the GEO satellite during repositioning maneuver can be determined and predicted by using C-Band pseudo-range observations of the BeiDou GEO satellite with COSPAR ID 2010-001A in 2011 and 2012. The results indicate that observations before maneuver do affect orbit determination and prediction, and should therefore be selected appropriately. A more precise orbit and prediction can be obtained compared to common short arc methods when observations starting 1 day prior the maneuver and 2 h after the maneuver are adopted in POD (Precise Orbit Determination). The achieved URE (User Range Error) under non-consideration of satellite clock errors is better than 2 m within the first 2 h after maneuver, and less than 3 m for further 2 h of orbit prediction.     

Direct position determination using single moving rotating linear array: Noncoherent and coherent processing

To improve the resolution and accuracy of Direct Position Determination (DPD), this paper investigates the problem of positioning multiple emitters directly with a single moving Rotating Linear Array (RLA). Firstly, the geometry of the RLA is formulated and analysed. According to its geometry, the intercepted noncoherent signals in multiple interception intervals are modeled. Correspondingly, the MUltiple SIgnal Classification (MUSIC) based noncoherent DPD approach is proposed. Secondly, the synchronous coherent pulse signals are individually considered and formulated. And the coherent DPD approach which aims for localizing this special type of signal is presented by stacking all array responses at different interception intervals. Besides, we also derive the constrained Cramér-Rao Lower Bound (CRLB) expression for both noncoherent and coherent DPD with RLA under the constraint that the altitudes of the emitters are known. At last, computer simulations are included to examine the performance of the proposed approach. The results demonstrate that the localization accuracy and resolution of DPD with single moving linear array can be significantly improved by the array rotation. In addition, coherent DPD with RLA further improves the resolution and increases the maximum emitter number that can be localized compared with the noncoherent DPD with RLA.     

加筋壁板大开口结构稳定性分析及试验验证

为了研究加筋壁板大开口结构的压缩承载极限问题,规划了两组具有不同设计参数的试验件进行了静力试验,并同时采用了线性特征值法以及非线性的弧长法对其压缩极限进行了对比分析。在采用弧长法分析时,引入了模态扰动技术,而且考虑了材料的弹塑性性质。计算结果与试验的对比表明:针对复杂结构稳定性问题,屈曲特征值法具有较大局限性,而结合模态扰动的弧长法对这种复杂的加筋壁板类结构压缩极限的分析精度较好。通过分析给出了加筋壁板大开口结构设计的一些建议。     

Predictive fault-tolerant control of an all-thruster satellite in 6-DOF motion via neural network model updating

The problem of controlling an all-thruster spacecraft in the coupled translational-rotational motion in presence of actuators fault and/or failure is investigated in this paper. The nonlinear model predictive control approach is used because of its ability to predict the future behavior of the system. The fault/failure of the thrusters changes the mapping between the commanded forces to the thrusters and actual force/torque generated by the thruster system. Thus, the basic six degree-of-freedom kinetic equations are separated from this mapping and a set of neural networks are trained off-line to learn the kinetic equations. Then, two neural networks are attached to these trained networks in order to learn the thruster commands to force/torque mappings on-line. Different off-nominal conditions are modeled so that neural networks can detect any failure and fault, including scale factor and misalignment of thrusters. A simple model of the spacecraft relative motion is used in MPC to decrease the computational burden. However, a precise model by the means of orbit propagation including different types of perturbation is utilized to evaluate the usefulness of the proposed approach in actual conditions. The numerical simulation shows that this method can successfully control the all-thruster spacecraft with ON-OFF thrusters in different combinations of thruster fault and/or failure.     

Direction finding for two-dimensional incoherently distributed sources with Hadamard shift invariance in non-uniform orthogonal arrays

This paper proposes a novel algorithm for Two-Dimensional (2D) central Direction-of-Arrival (DOA) estimation of incoherently distributed sources. In particular, an orthogonal array structure consisting of two Non-uniform Linear Arrays (NLAs) is considered. Based on first-order Taylor series approximation, the Generalized Array Manifold (GAM) model can first be established to separate the central DOAs from the original array manifold. Then, the Hadamard rotational invariance relationships inside the GAMs of two NLAs are identified. With the aid of such relationships, the central elevation and azimuth DOAs can be estimated through a search-free polynomial rooting method. Additionally, a simple parameter pairing of the estimated 2D angular parameters is also accomplished via the Hadamard rotational invariance relationship inside the GAM of the whole array. A secondary but important result is a derivation of closed-form expressions of the Cramer-Rao lower bound. The simulation results show that the proposed algorithm can achieve a remarkably higher precision at less complexity increment compared with the existing low-complexity methods, which benefits from the larger array aperture of the NLAs. Moreover, it requires no priori information about the angular distributed function.