Calibration of the in-orbit center-of-mass of TaiJi-1
Abstract
Taiji program is a space mission aiming to detect gravitational waves in the low frequency band. Taiji-1 is the first technology demonstration satellite of the Taiji Program in Space, with the gravitational reference sensor (GRS) serving as one of its key scientific payloads. For accurate accelerometer measurements, the test-mass center of the GRS must be positioned precisely at the center of gravity of the satellite to avoid measurement disturbances caused by angular acceleration and gradient. Due to installation and measurement errors, fuel consumption during in-flight phase, and other factors, the offset between the test-mass center and the center-of-mass (COM) of the satellite can be significant, degrading the measurement accuracy of the accelerometer. Therefore, the offset needs to be estimated and controlled within the required range by the center-of-mass adjustment mechanism during the satellite’s lifetime. In this paper, we present a novel method, the Extended Kalman Filter combined with Rauch-Tung-Striebel Smoother, to estimate the offset, while utilizing the chi-square test to eliminate outliers. Additionally, the nonlinear Least Squares estimation algorithm is employed as a crosscheck to estimate the offset of COM. The two methods are shown to give consistent results, with the offset estimated to be mm, mm, and mm. The results indicate a significant improvement on the noise level of GRS after the COM calibration, which will be of great help for the future Taiji program.
I Introduction
In 2016, the LIGO collaboration made a groundbreaking discovery by detecting gravitational waves (GW) [1], thereby providing a direct verification of the predictions made by Albert Einstein in his general theory of relativity a century ago [2]. This discovery has had a profound impact on basic scientific research worldwide.
Space-based gravitational wave detection represents an intriguing frontier for future studies of the gravitational universe, as it can extend the reach of gravitational wave astronomy beyond that of the ground-based detectors, thereby allowing for a wider range of gravitational radiation sources to be observed [3], which will provide invaluable information to deepen our understanding of the evolution of early universe and the nature of gravity. Several space-borne gravitational-wave observatories have been proposed, such as LISA [4, 5, 6], DECIGO [7], ASTROD [8], Taiji [3, 9], and Tianqin [10].
The Taiji program, initiated by the Chinese Academy of Sciences (CAS), is a space mission that aims to detect gravitational waves (GW) in the frequency band between 0.1 mHz and 1.0 Hz, which is important in the fields of astronomy and cosmology [9, 11, 12]. The Taiji program proposes to detect GW signals using the Michelson laser interferometer principle, where each end of the interferometer contains a test mass (TM) serving as a reference body. This reference body is required to be free from spurious accelerations relative to its local inertial frame, and any spurious accelerations will affect the detection of tidal deformations caused by gravitational waves.
To facilitate the development of technology for the Taiji program, a three-step road map has been proposed [13, 14]. As the first step, a technology demonstrator satellite, Taiji-1 [14, 15], was launched on August 31, 2019. One of the key technologies validated by Taiji-1 is the Gravity Reference Sensor (GRS), which serves as an accelerometer and consists of sensors and electronic components [16]. The sensor comprises an electrode housing and a TM surrounded by the sensing electrode, as shown in Fig. 1.
GRS has three axes, including one non-sensitive axis and two sensitive axes, the directions of , , and in Fig. 1 correspond to the non-sensitive axis (radial direction), the first sensitive axis (flight direction), and the second sensitive axis, respectively. The sensor utilizes capacitive sensing technology to measure the disturbance acceleration of Taiji-1. The resulting data is sent to the drag-free controller, which instructs the thruster to apply force to compensate for the disturbing force. As the reference body of the future Taiji program interferometer, GRS needs to effectively mitigate the impact of all non-gravitational accelerations so that to reach the desired sensitivity of Taiji. Moreover, the accuracy of GRS is crucial for the drag-free control system, as the GRS readout results are used as inputs for issuing commands to control the spacecraft. GRS is susceptible to a range of noise sources, including Brownian noise arising from the surrounding air near TM, charge-induced noise due to TM charge accumulation, readout noise from voltage signals, temperature gradients, circuit noise, magnetic noise, self-gravity noise.
A precise positioning of the test-mass of the accelerometer at the center-of-mass (COM) of the Taiji-1 satellite is crucial to suppress non-gravitational and perturbed accelerations such as angular motion-related accelerations and accelerations due to gravity gradients [17, 18]. The COM of the accelerometer is adjusted to be at the COM of the satellite before launching. However, during the satellite’s operation, the consumption of propellant causes the COM of the satellite to change relative to the satellite frame, leading to a shift of the accelerometer COM relative to the satellite COM over time. Therefore, it is crucial to regularly measure the deviation of the COM position of the two during the entire life cycle of the satellite, and use the COM adjustment mechanism to perform in-orbit adjustments so that the deviation is within a certain range [19, 20, 21].
In Section II, the principle of the COM calibration is presented. The accelerometer measurement model is described in Section III. We discuss the use of the Extended Kalman filter model and Rauch-Tung-Striebel Smoother for COM calibration in Section IV, while the detection and removal of outliers principle is explained in Section V. The performance of the COM calibration is evaluated in Section VI, and the results and conclusions are summarized in Section VII.
II Principle of COM calibration during in-orbit
As shown in Fig. 1, the primary components of the high-precision electrostatic levitation accelerometer include a free fall TM, a set of capacitive electrode plate that surrounds the TM (together forming a sensitive probe), and a peripheral capacitive sensing and electrostatic feedback control circuit. This circuit enables the detection of position and attitude changes between the TM and the electrodes, as well as the measurement of acceleration through its feedback voltage. There are six sensing, control, and feedback circuits that use the same principle to measure three translational accelerations and three angular accelerations of the TM concurrently.
An offset of the COM of the electrostatic accelerometer from the COM of the satellite can cause measurement disturbances, primarily attributed to the angular motion of the satellite. Therefore, an estimate of the offset can be achieved using appropriate algorithms based on the relationship between the angular motion and the measurement disturbance.
The on-orbit center-of-mass calibration experiment is conducted using the Taiji-1 attitude control system, which consists of the attitude sensors and the actuators, as shown in Fig. 2. The additude sensors includes star sensors, three-axis magnetometers, and gyroscopes. The actuators include cold-gas micro-thrusters, Hall-effect micro-thrusters, magnetic torquers, and momentum wheel. The X-axis magnetic torquer is mounted on the side panel of the satellite, the Y-axis magnetic torquer is mounted on the bottom panel, and the Z-axis magnetic torquer is mounted on the top panel. The direction of the positive magnetic moment aligns with the satellite’s three axes, as shown in the left figure of Fig. 2.
In the calibration experiment, periodic torques, generated by the magnetic torquers, is applied to the satellite, inducing a cyclic oscillation by the Earth’s magnetic field, which will cause disturbance to the accelerometer measurements and ensure its magnitude is sufficiently large to disregard other disturbance effect, such as solar pressure torque and aerodynamic torque. Then, the offset can be modeled from readout of accelerometer and star tracker. Moreover, the accuracy of GRS is crucial for the drag-free control system, as the GRS readout results are used as inputs for issuing commands to control the spacecraft.
The TM is made of titanium alloy TC4, with a magnetic susceptibility of , therefore the effects of the coupling between the TM and the magnetic torquers is much smaller than the present GRS accuracy and can be neglected.
III Accelerometer measurement Model
The accelerometer measurement output model of the relative acceleration of TM and the electrode cage for Taiji-1 is presented as follows [15, 17, 18],
(1) |
where represents the theoretical measurements of the accelerometer, represent the COM offset between the accelerometer and the satellite, and denote the first and second partial derivation of relative to time. and are the angular velocity and the angular acceleration respectively. is the acceleration due to the gravitational gradient, the accelerometer measurement model ignores the perturbations due to the solid Earth tides, ocean tides, rotational deformations, the planets including the Sun and Moon, and general relativity which can be neglected as for the calibration span is short enough. Furthermore, represents the non-gravitational acceleration on the satellite, such as atmospheric drag, solar radiation pressure, and the Earth radiation pressure.
During the offset calibration period, the deviation of the COM of the TM from the COM of the satellite can be approximated as a constant offset, given the short measurement time. Therefore, can be expressed as,
(2) |
Considering the relatively smooth change of acceleration caused by non-conservative forces and gravity gradient, it can be approximated as a linear change within a limited time span. Therefore, the calibration interval is designed to be several minutes.
The attitude control system of the Taiji-1 satellite is equipped with a star tracker and a gyroscope to respectively measure the attitude angle of the satellite relative to the inertial coordinate system and the angular velocity . Additionally, the angular accelerations are obtained by a second-order polynomial fitting method.
Assuming negligible scale factors and misalignment errors, after substituting the and , the accelerometer measurement model can be expressed as follows,
(3) |
where can be expressed as
(4) |
The corresponding axis of the spacecraft can be represented by and () for angular acceleration and angular velocity, respectively. The linear slope is represented by and the constant bias by . In this study, we removed the linear effect by detrending . Therefore, the model of utilized in this study is,
(5) |
where is the measurement noise.
IV Extended Kalman filter model and Rauch-Tung-Striebel Smoother for COM calibration
Kalman filter is a high-efficiency recursive filter [23]. The filtering theory proposed by Kalman is only applicable to linear systems. An Extended Kalman Filter (EKF) was proposed in Ref. [24, 25], and can be applied to the nonlinear field.
The accelerometer measurement model used in this study is given by Eq. 5, and one can define the state variable as the COM offset and use the Kalman filter to estimate the offset. Here, the state vector is denoted as with for . The state equation, as derived in Ref. [20], is,
(6) |
Here, signify the step of filter, represents the state transition matrix from step to step . As no manipulation is performed on the state vector, the state transition matrix can be written as , where is the identity matrix.
The output of the accelerometer can be defined as the observation equation,
(7) |
where , and is the discrete measurement noise that satisfies,
(8) |
Here, denotes the variance matrix of the measurement noise, assuming that all are independent, unbiased, and possess finite variance which implies that is a diagonal matrix.
The predicted observation equation is given by,
(9) |
The estimated covariance matrix is given by,
(10) |
where denotes the system noise variance matrix, which is assumed to be zero in this study.
The Kalman gain is
(11) |
and the state update value after Kalman filter is,
(12) |
Here, represents the prior estimate value, and is the posterior estimate.
The update of the error covariance matrix is,
(13) |
After the Kalman filter is applied, the difference between the predicted and observed measurements is known as the filtered residual. The filtered residual can be used to evaluate the accuracy of the filter. The covariance matrix of the filtered residual provides information on the uncertainty of the estimation. Therefore, it is important to analyze both the filtered residuals and the covariance matrix of the filtered residuals , as follows,
(14) |
(15) |
The characteristic property of the Kalman filter is that the filtered residual vectors are uncorrelated, and even independent in the case of Gaussian distribution.
After applying the aforementioned Kalman Filter, the state estimator of the dynamic system can be further refined using the Rauch-Tung-Striebel (RTS) smoother, as described in Ref. [26]. The smoothing equations are given by the following formulas,
(16) |
(17) |
(18) |
(19) |
(20) |
(21) |
(22) |
The RTS smoother provides smoothed estimates of the state mean and state covariance at time step , denoted as and , respectively. The smoother gain on time step , denoted as , corrects the RTS smoother estimate. The recursion is initialized at the last time step T of Kalman filter with and . The smoothed residuals are denoted as , and the covariance matrix of smoothed residuals is denoted as .
To obtain a more accurate result, a combined method based on Kalman filter and the RTS Smoother (KF-RTS) and outliers removal is proposed. The data is filtered with the Kalman filter, and smoothed by the RTS smoother, which takes into account all the available valid data points. A chi-square confidence test is performed during the smoothing process for assessing the data quality and removing outliers. This KF-RTS process is iterated until there are no outliers.
A possible drawback of the filter algorithm is the fact that one needs an initial value of the state vector together with its covariance matrix. This can be obtained by fitting a small number of measurements at the start of the track by a conventional least-squares fit, but this is not an elegant solution. The other possibility is to start with an arbitrary state vector and an infinite covariance matrix, i.e. a large multiple of the identity matrix. This is completely in the spirit of the filtering approach, but may lead to numerical instabilities in the computation of the gain matrix, since the infinities have to cancel in order to give a finite gain matrix. This may be difficult on a computer with a short word length.
In this article the initial value of the state vector is set to zero and its covariance is chosen to be 0.001, which is large enough. The measurement noise is calculated by selecting a segment of stationary data from the corresponding data.
Here, we use the nonlinear least squares (NLLS) method [27] for parameter estimation, as a crosscheck of the KF-RTS method. As is well-known, the key challenge for NLLS is to find the value of that minimizes the function
(23) |
where . The Levenberg-Marquardt (LM) algorithm [28, 29], which is an efficient method to minimize , is used to find the optimal parameters in this article.
V Detection and removal of outliers
During the Kalman filtering and RTS smoothing, each data point was utilized to obtain an optimal estimate of the state value, and this process also allows for an assessment of the quality of each data point.
The chi-square value [30] is commonly used for assessing the quality of data and detecting outliers, which can be caused by spacecraft maneuvering, non-Gaussian noise, electronic noise, or other coupled noise, and deviate significantly from the normal sequence of measurements.
The residuals of the global fit can be utilized to identify measurements with large residuals as potential outliers. Kalman filters and smoother enable the exploitation of complete information locally, to determine the validity of a measurement with high probability. The smoothed residual chi-square value can serve as a useful decision criterion for the data quality of the measurement,
(24) |
It is demonstrated that the test on smoothed residual chi-square consistently outperforms the test on filtered residual chi-square [31]. Thus, searching for possible outliers should be carried out during smoothing, as it allows the utilization of complete parameter information.
During the smoothing process, the measurement point can be removed from the smoothed estimate to obtain , which represents the optimal estimate of the system state at step using all data information except . This optimal estimate can be utilized for the detection and removal of outliers. To remove from the estimate , an inverse Kalman filter can be applied with the covariance matrix of taken as negative. This step of the filter is described in Ref. [31], and the smoothed estimate of without using , , can be calculated as,
(25) |
in which
(26) |
and
(27) |
If is a valid measurement and the covariance matrix of its Gaussian readout error is known, the quantity follows a chi-square distribution with degrees of freedom, where is the dimension of .
The measurement can be identified as an outlier if the value of exceeds a certain threshold . This threshold is chosen as the quantile of the corresponding distribution, where represents the probability of rejecting a valid measurement and is chosen to be in this paper. Other choices of value, 0.005, 0.01, 0.05 and 0.1 have been tried, and the effects on the estimated COM offset is found to be negligible.
The measurement can be removed permanently from the list as an outlier, and the RTS smoother can continue with and instead of and for updating the estimates when . To remove all outliers, this Kalman filter and RTS smoother must be recomputed without the outliers and iterated until convergence is achieved.
VI Results of COM calibration
The COM bias estimation in this paper is performed using the two algorithms discussed earlier. Figure 3 presents a comparison between the linear acceleration results obtained from the NLLS estimation and the original data, while Fig. 4 illustrates the comparison between the linear acceleration results obtained from the extended Kalman filter algorithm and the original data. Both methods exhibit excellent agreement with the original data. The results of COM calibration using the KF-RTS algorithm before removing outliers are shown in Fig. 5, with the shaded area indicating the range of one standard deviation.
Figure 6 illustrates the comparison between the fit results obtained using the NLLS algorithm and the experimental data with outliers removed. The final round results of the offset obtained using the KF-RTS algorithm are presented in Fig. 7, with the shaded area indicating the range of one standard deviation. Figure 8 displays the comparison between the final round results obtained using the KF-RTS algorithm and the experimental data.
Table 1 presents the estimation results for the COM offset obtained using the two methods. Despite the use of different estimation algorithms, the results obtained from both methods are highly consistent, thus increasing the confidence in the accuracy of the results. It is worth mentioning that, the value of as the goodness of a fit is 2.33 for the first round and is 0.88 for the final round, where represents the number of degrees of freedom.
Axis | X (m) | Y (m) | Z (m) |
---|---|---|---|
NLLS with outliers | 175 41 | 610 101 | 777 31 |
First KF-RTS | 173 27 | 606 66 | 777 20 |
NLLS without outliers | 191 27 | 643 66 | 818 20 |
Final KF-RTS | 189 27 | 638 68 | 818 20 |
Figure 9 and Figure 10 depict the comparison of the amplitude spectral density (ASD) between the calibration experimental data and flight data before and after the COM calibration. Figure 9 shows that the peak observed in the experimental data is caused by a modulation signal, and its influence is highly suppressed after the COM calibration, which is consistent with our expectations. Figure 10 demonstrates a significant reduction in the acceleration noise level of GRS readings in the frequency range of 0.001–0.1 Hz after the COM calibration.
VII Discussion and Conclusions
Detecting and reducing the deviation between the COM of the inspection load and the COM of the satellite is crucial for high-precision accelerometers, as it can improve their accuracy. It is also a significant step in the development of space-based gravitational wave observatory to achieve their scientific objectives. Furthermore, the calibration can help obtain a high-precision gravity field, which is valuable for conducting geoscience research with greater accuracy.
In this study, the offset of the COM between the inspection load and the satellite is estimated using the KF-RTS smoother. Outliers are detected using the chi-square test, and the inverse Kalman filter is applied to remove them. The LM algorithm, as a cross check, is used to find the optimal offset parameters for the NLLS method. The results obtained with both methods are in very good agreement, and the offset of COM between the inspection load and the satellite is estimated with an accuracy of . After obtaining the COM offset, one can reduce it using the COM adjustment mechanism, and the effects of the COM offset can also be suppressed in the data-processing. The COM calibration is crucial for improving the accuracy of the accelerometer which will directly impact the detection sensitivity of the final space-based gravitational wave observatory.
For the space-borne gravitational-wave observatories, such as Taiji-3, there are three satellites and each satellite is equipped with two TMs, which means that the COM of TMs does not coincide with that of the satellite. However, as a reference body for the satellite, it would cause residual acceleration if the TMs are away from their nominal positions. Therefore, it is necessary to periodically estimate or monitor the deviation of the COM of TM from a fixed point and make adjustments when necessary. The same calibration principle and methods as discussed in this paper can be used. The Taiji-3 satellite will be equipped with higher-precision star-sensitive instruments, which will enable more accurate results.
Acknowledgements
This work is partially supported by the Strategic Priority Research Program of the Chinese Academy of Sciences Grant No. XDA15020700 and XDA15021100, and by the Fundamental Research Funds for the Central Universities. We acknowledge support from the National Space Science Data Center, National Science and Technology Infrastructure of China.
References
- Abbott et al. [2016] B. P. Abbott et al. (LIGO Scientific Collaboration and Virgo Collaboration), Phys. Rev. Lett. 116, 061102 (2016).
- Einstein [1916] A. Einstein, Sitzungsber. Preuss. Akad. Wiss. Berlin (Math. Phys.) 1916, 1 (1916).
- Hu and Wu [2017] W.-R. Hu and Y.-L. Wu, Natl. Sci. Rev. 4, 685 (2017).
- Bender et al. [1998] P. Bender et al., Max-Planck-Institut für Quantenoptik, Garching (1998).
- Jennrich et al. [2011] O. Jennrich et al., ESA/SRE (2011) 19 (2011).
- Amaro-Seoane et al. [2017] P. Amaro-Seoane et al., arXiv 1702.00786 (2017).
- Seto et al. [2001] N. Seto, S. Kawamura, and T. Nakamura, Phys. Rev. Lett. 87, 221103 (2001).
- Ni [2013] W.-T. Ni, Int. J. Mod. Phys. D 22, 1341004 (2013).
- Cyranoski et al. [2016] D. Cyranoski et al., Nature 531, 150 (2016).
- Luo et al. [2016] J. Luo et al., Classical and Quantum Gravity 33, 035010 (2016).
- Sathyaprakash and Schutz [2009] B. S. Sathyaprakash and B. F. Schutz, Living reviews in relativity 12, 1 (2009).
- Schutz [1999] B. F. Schutz, Classical and Quantum Gravity 16, A131 (1999).
- Luo et al. [2020a] Z. Luo, Z. Guo, G. Jin, Y. Wu, and W. Hu, Results in Physics 16, 102918 (2020a).
- Luo et al. [2020b] Z. Luo, Y. Wang, Y. Wu, W. Hu, and G. Jin, Progress of Theoretical and Experimental Physics 2021, 05A108 (2020b).
- Cai and Deng [2021] Z. Cai and J. Deng, Int. J. Mod. Phys. A 36, 2140020 (2021).
- Touboul et al. [1999] P. Touboul, E. Willemenot, B. Foulon, and V. Josselin, Boll. Geof. Teor. Appl 40, 321 (1999).
- Wang et al. [2010] F. Wang, S. Bettadpur, H. Save, and G. Kruizinga, Journal of Spacecraft and Rockets 47, 371 (2010).
- Huang et al. [2022] Z. Huang, S. Li, L. Cai, D. Fan, and L. Huang, Remote Sensing 14 (2022).
- Armano et al. [2018] M. Armano et al., Phys. Rev. D 97, 122002 (2018), arXiv:1806.08581 [astro-ph.IM] .
- Dong et al. [2009] F. Dong, H. Liao, C. Jia, and X. Xia, Science in China Series E: Technological Sciences 52, 1446 (2009).
- Guzman Cervantes et al. [2008] F. Guzman Cervantes, F. Steier, G. Wanner, G. Heinzel, and K. Danzmann, Applied Physics B 90, 395 (2008).
- Peng et al. [2021] X. Peng et al., Int. J. Mod. Phys. A 36, 2140026 (2021).
- Kalman [1960] R. E. Kalman, Journal of Basic Engineering 82, 35 (1960).
- Sunahara [1970] Y. Sunahara, Journal of Basic Engineering 92, 385 (1970).
- Bucy and Senne [1971] R. S. Bucy and K. D. Senne, Automatica 7, 287 (1971).
- Rauch et al. [1965] H. E. Rauch, F. Tung, and C. T. Striebel, AIAA journal 3, 1445 (1965).
- Al-Baali and Fletcher [1985] M. Al-Baali and R. Fletcher, Journal of the Operational Research Society 36, 405 (1985).
- Levenberg [1944] K. Levenberg, Quarterly of applied mathematics 2, 164 (1944).
- Marquardt [1963] D. W. Marquardt, Journal of the Society for Industrial and Applied Mathematics 11, 431 (1963).
- Cochran [1952] W. G. Cochran, The Annals of mathematical statistics , 315 (1952).
- Fruhwirth [1987] R. Fruhwirth, Nucl. Instrum. Meth. A 262, 444 (1987).