Hardware-in-the-loop simulations of GPS-based navigation and control for satellite formation flying

A relative navigation and formation control algorithm for satellite formation flying was developed, and a hardware-in-the-loop (HIL) simulation testbed was established and configured to evaluate this algorithm. The algorithm presented is a relative navigation estimation algorithm using double-difference carrier-phase and single-difference code measurements based on the extended Kalman filter (EKF). In addition, a state-dependent Riccati equation (SDRE) technique is utilized as a nonlinear controller for the formation control problem. The state-dependent coefficient (SDC) form is formulated to include nonlinearities in the relative dynamics. To evaluate the relative navigation and control algorithms developed, a closed-loop HIL testbed is configured. To demonstrate the performance of the testbed, a test formation flying scenario comprising formation acquisition and keeping in a low earth orbit (LEO) has been established. The relative navigation results from the closed-loop simulations show that a 3D RMS of 0.07 m can be achieved for position accuracy. The targeted leader–follower formation flying in the along-track separation of 100 m was maintained with a mean position error of approximately 0.2 m and a standard deviation of 0.9 m. The simulation results show that the HIL testbed is capable of successful demonstration of the GPS-based satellite autonomous formation flying mission.     

Integrated orbit and attitude hardware-in-the-loop simulations for autonomous satellite formation flying

Development and experiment of an integrated orbit and attitude hardware-in-the-loop (HIL) simulator for autonomous satellite formation flying are presented. The integrated simulator system consists of an orbit HIL simulator for orbit determination and control, and an attitude HIL simulator for attitude determination and control. The integrated simulator involves four processes (orbit determination, orbit control, attitude determination, and attitude control), which interact with each other in the same way as actual flight processes do. Orbit determination is conducted by a relative navigation algorithm using double-difference GPS measurements based on the extended Kalman filter (EKF). Orbit control is performed by a state-dependent Riccati equation (SDRE) technique that is utilized as a nonlinear controller for the formation control problem. Attitude is determined from an attitude heading reference system (AHRS) sensor, and a proportional-derivative (PD) feedback controller is used to control the attitude HIL simulator using three momentum wheel assemblies. Integrated orbit and attitude simulations are performed for a formation reconfiguration scenario. By performing the four processes adequately, the desired formation reconfiguration from a baseline of 500–1000 m was achieved with meter-level position error and millimeter-level relative position navigation. This HIL simulation demonstrates the performance of the integrated HIL simulator and the feasibility of the applied algorithms in a real-time environment. Furthermore, the integrated HIL simulator system developed in the current study can be used as a ground-based testing environment to reproduce possible actual satellite formation operations.     

Response of the mid-latitude D-region ionosphere to the total solar eclipse of 22 July 2009 studied using VLF signals in South Korean peninsula

In this paper, we analyze VLF signals received at Busan to study the the D-region changes linked with the solar eclipse event of 22 July 2009 for very short (∼390 km) transmitter–receiver great circle path (TRGCP) during local noon time 00:36–03:13 UT (09:36–12:13 KST). The eclipse crossed south of Busan with a maximum obscuration of ∼84%. Observations clearly show a reduction of ∼6.2 dB in the VLF signal strength at the time of maximum solar obscuration (84% at 01:53 UT) as compared to those observed on the control days. Estimated values of change in Wait ionospheric parameters: reflection height (h′) in km and inverse scale height parameter (β) in km⁻¹ from Long Wave Propagation Capability (LWPC) model during the maximum eclipse phase as compared to unperturbed ionosphere are 7 km and 0.055 km⁻¹, respectively. Moreover, the D-region electron density estimated from model computation shows 95% depletion in electron density at the height of ∼71 km. The reflection height is found to increase by ∼7 km in the D-region during the eclipse as compared to those on the control days, implying a depletion in the Lyman-α flux by a factor of ∼7. The present observations are discussed in the light of current understanding on the solar eclipse induced D-region dynamics.     

Investigating suitable orbits for the Swarm constellation mission – The frozen orbit

This paper proposes a suitable orbit design for the lower pair of ESA's Swarm constellation mission, flying side-by-side in near-polar and circular orbits with a separation of only 1.4° at ascending node. Both orbits are suggested to be frozen orbits to minimize the evolution, and an along-track separation strategy is applied to avoid collision risk. The characteristics of the proposed orbit type are examined through numerical techniques including high-fidelity perturbation models. The prime change from the initial configuration is an along-track separation. The perturbations causing the along-track drift are analyzed by switching on/off certain perturbations. The results indicate that the tesseral harmonics and the atmospheric drag yield dominant effects. The atmospheric drag effect shows a dependence on the local time of the ascending node. From two months of orbit propagation for the altitude 300 km the maximum along-track drift we obtain is about 80 km, which is still within the measurement requirement range. Several maneuver strategies for maintaining the proposed orbit design are suggested. The results analyzed for the proposed orbit design show that collision risk can be avoided by along-track separation within the frozen orbit design. Consequently, this combination is considered as a suitable approach for Swarm's lower pair.