State and disturbance estimation for test masses of drag-free satellites based on self-recurrent wavelet neural network

In recent years, the drag-free satellites have been widely used for some fundamental physical experiments, such as checking short-range effects of general relativity, geopotential determination and the exploration of static ocean current. And the space-borne detector of gravitational waves is one of the important applications for drag-free satellites in the future. In this study, the estimation of relative motion state and disturbance for test masses of drag-free satellite after release are researched. Firstly, the relative motion model between the test masses is established based on corresponding reference frames. Secondly, a self-recurrent wavelet neural network estimator is designed to estimate the pure gravitational relative motion state and disturbance for the test masses, and the sliding mode controller is used to transfer the relative motion state to the ideal state. Finally, the proposed estimation methods are verified by the simulation results.     

Modeling and identification of an electro-mechanical system: The LISA grabbing positioning and release mechanism case

In the frame of space missions, mechanisms often constitute critical systems whose functionality and performance need to be tested on ground before the mission launch. The LISA scientific space mission will detect gravitational waves by measuring the relative displacement of pairs of free-floating test masses set into geodesic motion onboard of three spacecrafts. Inside each satellite, the injection of the test masses from the caged configuration into the geodesic trajectory will be performed by the grabbing positioning and release mechanism. To provide a successful injection, the test masses must be dynamically released with a minimal residual velocity against adhesion with the holding device. A parameter that determines the test mass residual velocity is the quickness of the retraction of the holding device. The need arises then to characterize the dynamic response of the release mechanism in order to predict its behaviour in the in-flight conditions. Once a validated model of the mechanism is available, the compliance of the system to the tight requirement on the maximum allowed residual velocity of the test mass may be verified. Starting from an electro-mechanical model of the mechanism dynamics, this paper presents the results of the experimental identification of its relevant parameters.     

Detecting gravitational waves from accreting neutron stars

The gravitational waves emitted by neutron stars carry unique information about their structure and composition. Direct detection of these gravitational waves, however, is a formidable technical challenge. In a recent study we quantified the hurdles facing searches for gravitational waves from the known accreting neutron stars, given the level of uncertainty that exists regarding spin and orbital parameters. In this paper we reflect on our conclusions, and issue an open challenge to the theoretical community to consider how searches should be designed to yield the most astrophysically interesting upper limits. With this in mind we examine some more optimistic emission scenarios involving spin-down, and show that there are technically feasible searches, particularly for the accreting millisecond pulsars, that might place meaningful constraints on torque mechanisms. We finish with a brief discussion of prospects for indirect detection.     

Gravitational wave detection on the Moon and the moons of Mars

The Moon and the moons of Mars should be extremely quiet seismically and could therefore become sensitive gravitational wave detectors, if instrumented properly. Highly sensitive displacement sensors could be deployed on these planetary bodies to monitor the motion induced by gravitational waves. A superconducting displacement sensor with a 10-kg test mass cooled to 2 K will have an intrinsic instrument noise of 10⁻¹⁶ m Hz^−1/2. These sensors could be tuned to the lowest two quadrupole modes of the body or operated as a wideband detector below its fundamental mode. An interesting frequency range is 0.1–1 Hz, which will be missed by both the ground detectors on the Earth and LISA and would be the best window for searching for stochastic background gravitational waves. Phobos and Deimos have their lowest quadrupole modes at 0.2–0.3 Hz and could offer a sensitivity h_min ? 10⁻²² Hz^−1/2 within their resonance peaks, which is within two orders of magnitude from the goal of the Big Bang Observer (BBO). The lunar and Martian moon detectors would detect many interesting foreground sources in a new frequency window and could serve as a valuable precursor for BBO.     

Free fall tests of the accelerometers of the MICROSCOPE mission

The MICROSCOPE mission is fully dedicated to the in-orbit test of the Universality of free fall, the so-called Weak Equivalence Principle (WEP), with an expected accuracy better than 10⁻¹⁵. The test principle consists in comparing the accelerations of two proof masses of different composition in the Earth gravitational field. The payload embarks two pairs of test-masses made of Platinum Rhodium and Titanium alloys at the core of two dedicated coaxial electrostatic accelerometers. These instruments are under qualification for a launch in 2016.     

Dark energy,co-evolution of massive black holes with galaxies,and ASTROD-GW

The detection of low frequency band (100 nHz–100 mHz) and very low frequency band (300 pHz–100 nHz) gravitational waves (GWs) is important for exploration of the equation of state of dark energy and the co-evolution of massive black holes (MBHs) with galaxies. Most galaxies are believed to have a massive black hole in the galactic core. In the formation of these black holes, merging and accretion are the two main processes. Merging of massive black holes generate GWs which could be detected by space GW detectors and Pulsar Timing Arrays (PTAs) to cosmological distances. LISA (Laser-Interferometric Space Antenna) is most sensitive to the frequency band 1 mHz–100 mHz, ASTROD-GW (ASTROD [Astrodynamical Space Test of Relativity using Optical Devices] optimized for Gravitational Wave detection) is most sensitive to the frequency band 100 nHz–1 mHz and PTAs are most sensitive to the frequency band 300 pHz–100 nHz. In this paper, we discuss the sensitivities and outlooks of detection of GWs from binary massive black holes in these frequency bands with an emphasis on ASTROD-GW. The GWs generated by the inspirals, merging and subsequent ringdowns of binary black holes are standard sirens to the cosmological distance. Using GW observations, we discuss the methods for determining the equation of state of dark energy and for testing the co-evolution models of massive black holes. ASTROD-GW is an optimization of ASTROD to focus on the goal of detection of GWs. The mission orbits of the 3 spacecraft forming a nearly equilateral triangular array are chosen to be near the Sun-Earth Lagrange points L3, L4 and L5. The 3 spacecraft range interferometrically with one another with arm length about 260 million kilometers. With 52 times longer in arm length compared to that of LISA, the strain detection sensitivity is 52 times better toward larger wavelength. The scientific aim is focused for gravitational wave detection at low frequency. The science goals include detection of GWs from MBHs, and Extreme-Mass-Ratio Black Hole Inspirals (EMRI), and using these observations to find the evolution of the equation of state of dark energy and to explore the co-evolution of massive black holes with galaxies.     

Performance assessment of RTPPP positioning with SSR corrections and PPP-AR positioning with FCB for multi-GNSS from MADOCA products

The state-space representation (SSR) product of satellite orbit and clock is one of the most essential corrections for real-time precise point positioning (RTPPP). When it comes to PPP ambiguity resolution (PPP-AR), the fractional cycle bias (FCB) matters. The Japan Aerospace Exploration Agency (JAXA) has developed a multi-GNSS (i.e., global navigation satellite system) advanced demonstration tool for orbit and clock analysis (MADOCA), providing free and precise orbit and clock products. Because of the shortage of relevant studies on performance evaluation, this paper focuses on the performance assessment of RTPPP and PPP-AR by real-time and offline MADOCA products. To begin with, the real-time MADOCA products are evaluated by comparing orbit and clock with JAXA final products, which gives an objective impression of the correction. Second, PPP tests in static and simulated kinematic mode are conducted to further verify the quality of real-time MADOCA products. Finally, the offline MADOCA products are assessed by PPP and PPP-AR comparisons. The results are as follows: (1) Orbit comparisons produced an average error of about 0.04–0.13 m for the global positioning system (GPS), 0.14–0.16 m for the global navigation satellite system (GLONASS), and 0.07–0.08 m for the quasi-zenith satellite system (QZSS). The G15 satellite had the most accurate orbit, with a difference of 0.04 m between the JAXA orbit products and MADOCA’s counterpart, while the R07 satellite had the least accurate orbit with a difference of 0.16 m. Clock products had an accuracy of 0.4–1.3 ns for GPS, 1.4–1.6 ns for GLONASS, and 0.7–0.8 ns for QZSS in general. The G15 satellite had the most accurate clock with a difference of only 0.40 ns between the JAXA clock products and MADOCA products, and the R07 satellite had the least accurate clock with a difference of 1.55 ns. The orbit and clock products for GLONASS performed worse than those of GPS and QZSS. (2) After convergence, the positioning accuracy was 3.0–8.1 cm for static PPP and 8.1–13.7 cm for kinematic PPP when using multi-GNSS observations and precise orbit and clock products. The PFRR station performed the good performance both in static and kinematic mode with an accuracy of 2.99 cm and 8.08 cm, respectively, whereas the CPNM station produced the worst static performance with an error of 8.09 cm, and the ANMG station produced the worst kinematic performance with a counterpart of 13.69 cm. (3) The PPP-AR solution was superior to the PPP solution, given that, with respect to PPP, post-processing PPP-AR improved the positioning accuracy and convergence time by 13–32 % (3–89 %) in GPS-only mode by 2–15 % (5–60 %) in GPS/QZSS mode. Thus, we conclude that the current MADOCA products can provide SSR corrections and FCB products with positioning accuracy at the decimeter or even centimeter level, which could meet the demands of the RTPPP and PPP-AR solutions.     

In-flight testing of the injection of the LISA Pathfinder test mass into a geodesic

LISA Pathfinder is a technology demonstrator space mission, aimed at testing key technologies for detecting gravitational waves in space. The mission is the precursor of LISA, the first space gravitational waves observatory, whose launch is scheduled for 2034. The LISA Pathfinder scientific payload includes two gravitational reference sensors (GRSs), each one containing a test mass (TM), which is the sensing body of the experiment. A mission critical task is to set each TM into a pure geodesic motion, i.e. guaranteeing an extremely low acceleration noise in the sub-Hertz frequency bandwidth. The grabbing positioning and release mechanism (GPRM), responsible for the injection of the TM into a geodesic trajectory, was widely tested on ground, with the limitations imposed by the 1-g environment. The experiments showed that the mechanism, working in its nominal conditions, is capable of releasing the TM into free-fall fulfilling the very strict constraint imposed on the TM residual velocity, in order to allow its capture on behalf of the electrostatic actuation.However, the first in-flight releases produced unexpected residual velocity components, for both the TMs. Moreover, all the residual velocity components were greater than maximum value set by the requirements. The main suspect is that unexpected contacts took place between the TM and the surroundings bodies. As a consequence, ad hoc manual release procedures had to be adopted for the few following injections performed during the nominal mission. These procedures still resulted in non compliant TM states which were captured only after impacts. However, such procedures seem not practicable for LISA, both for the limited repeatability of the system and for the unmanageable time lag of the telemetry/telecommand signals (about 4400 s). For this reason, at the end of the mission, the GPRM was deeply tested in-flight, performing a large number of releases, according to different strategies. The tests were carried out in order to understand the unexpected dynamics and limit its effects on the final injection. Some risk mitigation maneuvers have been tested aimed at minimizing the vibration of the system at the release and improving the alignment between the mechanism and the TM. However, no overall optimal release strategy to be implemented in LISA could be found, because the two GPRMs behaved differently.     

The 16-day waves in the mesosphere and lower thermosphere over Wuhan (30.6°N, 114.5°E) and Adelaide (35°S, 138°E)

Winds from a meteor radar at Wuhan (30.6°N, 114.5°E) and a MF radar at Adelaide (35°S, 138°E) are used to study the 16-day waves in the mesosphere and lower thermosphere (MLT). The height range is 78–98 km at Wuhan and 70–98 km at Adelaide. By comparison, it is found that the zonal components at both sites are generally larger than the meridional ones, and eastward motion of the zonal background winds is favorable for the 16-day waves penetration to the MLT region. The zonal maximum amplitude appears in the autumn (September–October) around 86–98 km at Wuhan and in the winter months and early spring (July–October) around 72–82 km at Adelaide. Differences are found in wave amplitudes and time of appearance between the two years of 2002 and 2003. In 2003, the intensity of the wave amplitudes is relatively smaller than that for 2002 at both sites. The summer 16-day waves are comparatively weaker at Adelaide in both years, but stronger in 2002 at Wuhan near the mesopause and the lower thermosphere (86–98 km). The strong summer waves at Wuhan may come from the winter southern hemisphere.     

Utilization of CAM,CAF, MAX,and FAX for space radiation analyses using HZETRN

To estimate astronaut health risk due to space radiation, one must have the ability to calculate various exposure-related quantities that are averaged over specific organs and tissue types. Such calculations require computational models of the ambient space radiation environment, particle transport, nuclear and atomic physics, and the human body. While significant efforts have been made to verify, validate, and quantify the uncertainties associated with many of these models and tools, relatively little work has focused on the uncertainties associated with the representation and utilization of the human phantoms. In this study, we first examine the anatomical properties of the Computerized Anatomical Man (CAM), Computerized Anatomical Female (CAF), Male Adult voXel (MAX), and Female Adult voXel (FAX) models by comparing the masses of various model tissues used to calculate effective dose to the reference values specified by the International Commission on Radiological Protection (ICRP). The MAX and FAX tissue masses are found to be in good agreement with the reference data, while major discrepancies are found between the CAM and CAF tissue masses and the reference data for almost all of the effective dose tissues. We next examine the distribution of target points used with the deterministic transport code HZETRN (High charge (Z) and Energy TRaNsport) to compute mass averaged exposure quantities. A numerical algorithm is presented and used to generate multiple point distributions of varying fidelity for many of the effective dose tissues identified in CAM, CAF, MAX, and FAX. The point distributions are used to compute mass averaged dose equivalent values under both a galactic cosmic ray (GCR) and solar particle event (SPE) environment impinging isotropically on three spherical aluminum shells with areal densities of 0.4 g/cm², 2.0 g/cm², and 10.0 g/cm². The dose equivalent values are examined to identify a recommended set of target points for each of the tissues and to further assess the differences between CAM, CAF, MAX, and FAX. It is concluded that the previously published CAM and CAF point distributions were significantly under-sampled and that the set of point distributions presented here should be adequate for future studies involving CAM, CAF, MAX, or FAX. It is also found that the errors associated with the mass and location of certain tissues in CAM and CAF have a significant impact on the mass averaged dose equivalent values, and it is concluded that MAX and FAX are more accurate than CAM and CAF for space radiation analyses.     

数学规划的物理模型及算法

利用力学原理，建立了用微分方程描述的带不等式约束优化问题的物理模型，用人工释能和摄动的思想，推广了现行的一些可行方向算法，并对算法的收敛性进行了讨论。     

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.     

OPTIS – A Satellite test of Special and General Relativity

OPTIS has been proposed as a small satellite platform in a high elliptical orbit (apogee 40,000 km, perigee 10,000 km) and is designed for high precision tests of foundations of Special and General Relativity. The experimental set-up consists of two ultrastable Nd:YAG lasers, three crossed optical resonators (monolithic cavities), an ensemble of atomic clocks, an optical comb generator, laser tracking devices and a drag-free control system. OPTIS enables improved tests of (1) the isotropy and (2) constancy of the speed of light, (3) special relativistic time dilation, (4) the universality of the gravitational redshift by comparison of various clocks, can measure (5) the absolute value of the gravitational redshift, (6) the Lense–Thirring effect and (7) the perigee advance and (8) can make a test of a hypothetical Yukawa part in the gravitational potential. To avoid any influence from atmospheric drag, solar radiation, or Earth albedo, the satellite needs drag-free control to depress the residual acceleration down to 10⁻¹⁴ m/s² in the frequency range between 10⁻² and 10⁻³ Hz. Precise thermal control must be used to stabilize the cavity temperature to within one part in 10⁷ at time scales of 100 s and to one part in 10⁵ on the orbit time scale.     

Low-mass bare strange stars

Strange stars with low masses are suggested to exist in reality, the origin of which could be via accretion-induced collapse of white dwarfs. Such a strange star is likely bare, and would thus spin very fast, even to a period of <0.1 ms. Strange stars with low masses may differ from those with solar masses in various astrophysical appearances. Observations to test this “low-mass” idea are proposed.     

The day-to-day variability in the mesosphere and lower thermosphere in low latitudes: A study using MF radar

This study presents the analysis of planetary waves (PWs) using daily mean wind velocities for four years (August 2013 to July 2017) of continuous measurements using MF radar over the low latitude Indian region Kolhapur (16.8° N; 74.2° E). The MF radar at Kolhapur was upgraded in 2013. These are the first results of PWs after the upgradation of MF radar. The seasonal and intra-seasonal variabilities of East-West (EW) traveling PWs in the MLT region have been studied. In the present work, the data was analyzed to study the waves with various periodicities (e.g. 3–4, 5–8, 15–17, and 30–60 days). The 3.5 day [Ultra-Fast Kelvin (UFK)] wave shows semiannual variability with burst like wave activity observed during the summer months and December solstice. In addition, it is observed to be stronger in the spring equinoctial period. A strong semiannual oscillation (SAO) has been observed in a 6.5-day wave with peaks near the equinoxes. Similar to SAO over the low latitude MLT region, the wave activity is stronger in April/May than in September/October. The 6.5-day waves are observed to be stronger when the background mean wind is westward. From the analysis, it has been seen that the period before and after the equinoctial period is favorable for the 6.5-day wave propagation. The 16-day wave has no significant seasonal dependence; instead, the waves spread to almost all seasons. The Madden-Julian Oscillations (MJOs) have been observed to be propagating with an average wind speed of ～ 5 m/s when the background mean wind is eastward. The occurrence of MJO is observed during the summer and winter months. These results are the first of their kind in two aspects: first, they show the PWs with enhanced altitude coverage covering up to 110 km, and second, they show the PWs not contaminated due to equatorial electro jet influence.     

Multi-wavelength study of a coronal mass ejection: a flare event from AR#9393

The active region, AR#9393, produced a number of intense flares during March–April 2001. In this paper, we report the analysis of an X1.1 flare event of April 2, 2001 and its associated coronal mass ejection. The timing and location of the H_α eruption, radio burst activities, and the onset of mass ejection suggest an energy release that occurred close to the surface of the sun. At this region, as shown by the magnetogram, X-ray and EUV images, the field configuration was complex and the 3-D extrapolation revealed the presence of a magnetic null point. Results also suggest that the energy release is followed by the magnetic reconnection between the low-lying loops near the separator point and outlying loops. This study provides the support for the magnetic break-out process to trigger the energy release in eruptive flare event.     

Linear stability of ring systems around oblate central masses

In this paper, we consider the stability of a ring of bodies of equal mass uniformly distributed around a large oblate central mass. The purpose of this and previous papers is to shed light on the stability of Saturn’s rings. Previous papers have been limited by the assumptions that (1) all ring bodies are at the same distance from the central body, (2) the central body acts like a point mass (i.e., is a perfect sphere), and (3) the ring bodies all have the same mass and are evenly spaced around the ring. The third limitation is probably the least important; as long as there are a large number of masses and the mass distribution is approximately uniform then the system should behave as a system of equally-spaced, equal-mass bodies. The first main purpose of this paper is to remove the second limitation. But, the paper also aims to address limitation (1). Recent computer simulations of single-ring systems have shown that the threshold for stability, as determined by a linear stability analysis, matches precisely the stability threshold predicted by simulation. In other words, a linear stability analysis while presumably just a mathematical abstraction actually tells us something quite real. Furthermore, simulations of multi-ring systems suggest that instability comes from azimuthal perturbations; small azimuthal changes are more destabilizing that small radial perturbations. Hence, in this paper, we also consider the situation where the central body consists not just of an oblate central mass but also incorporates a flat ring representing in aggregate all ring bodies at radii other than the one under consideration. The central oblate body together with a flat ring is modeled simply by introducing two oblateness terms to the gravitational potential associated with the central mass. The subsequent analysis is almost identical to the case of a single oblateness term. For Saturn, the oblateness of the central mass is six orders of magnitude more significant than the rings at other radii as a destabilizing influence.     

Stacked recurrent neural network based high precision pointing coupled control of the spacecraft and telescopes

In some space missions especially in the field of space gravitational wave detection, the telescope needs to point to a certain target through attitude movement and pointing control. In several mainstream gravitational wave detection missions, the detector usually consists of a cluster of three identical satellites, flying in a quasi-equilateral triangular formation with a big edge length, so every satellite needs two telescopes to point each other and constitute three giant Michelson-Type interferometers. Therefore, a satellite platform system with two telescopes is researched in this paper. This research helps to characterize the attitude motion of a telescope for space astronomical observation or space gravitational wave detection, provides new method on the telescope’s high-precision pointing control. For this purpose, we derive a satellite-telescope coupling attitude model, design the sliding mode controller for satellite and the stacked recurrent neural network adaptive controller for telescope. In the stacked recurrent neural network adaptive controller design, a sliding mode control technology is adopted. In addition, we propose a combinatorial optimization method for network weights in the stacked recurrent neural network training process, that is, the output layer is corrected by the adaptive law, and the correction of other layers adopt the error backpropagation method. Finally, a numerical simulation method verifies the effectiveness of the controller design.     

Analytic and numerical treatment of motion of dust grain particle around triangular equilibrium points with post-AGB binary star and disc

This paper investigates the motion around the triangular equilibrium points, of a passively gravitating dust particle in the gravitational field of a low-mass post-AGB binary system, surrounded by circumbinary disc. The two bodies of the binary are modeled as a triaxial star and a radiating-oblate star. Due to small deviation of disc stars on circular orbits, we have assumed that the Coriolis and centrifugal forces of the stars are slightly perturbed. The triangular equilibrium points of the particle are found. These points are defined by, triaxiality of the primary star, oblateness and radiation of the secondary one and the gravitational potential from the disc mass. Further, when the disc mass increases, the particle moves nearer to the stars and farther away from the disc. In general, these equilibrium points are linearly stable when μ < μ_C; where μ is the mass ratio and μ_C is the critical mass function, defined by the parameters of the system. The effects of each of these parameters on the size of the stability region are stated, and the periodic motion around the stable points is examined. It is seen that the orbits are ellipses, and the orientation, eccentricities, lengths of the semi-major and semi-minor axes are influenced by the parameters of the problem. In particular, for our numerical linear stability analysis, we have taken an extremely depleted pulsating star, IRAS 11472-0800 as the post-AGB triaxial star, with a weakly-radiating young white dwarf star; G29-38 as the secondary. For this system, the stability result of the triangular points comes out different. Here, μ_C < μ throughout the entire range of the mass ratio and the critical mass function. Hence, the triangular equilibrium points are unstable. The stability of the orbits is tested using the Poincaré surfaces of section (Pss). The region of stability is controlled by the introduced parameters and the Jacobi constant.     

Detection of thermal X-ray emission in the halo of the plerionic supernova remnant G21.5-0.9

The detection of a soft thermal X-ray component in the spectrum of a bright knot in the halo of the plerion G21.5-0.9 is reported. Using a collisional ionization equilibrium model for an hot optically thin plasma, a temperature of kT 0.12–0.24 keV, a mass of 0.3–1.0 M_⊙ and a density of 1.6–6 cm⁻³ is derived. The spectral analysis suggests a possible overabundance of Silicon with respect to the solar value in the knot; if this will be confirmed this object may be a clump of shocked ejecta.