Next Article in Journal
Defect Detection inside a Rail Head by Ultrasonic Guided Waves
Next Article in Special Issue
A Gain Scheduling Design Method of the Aero-Engine Fuel Servo Constant Pressure Valve with High Accuracy and Fast Response Ability
Previous Article in Journal
Convergence of Inexact Iterates of an Algorithm Based on Unions of Nonexpansive Mappings in Metric Spaces
Previous Article in Special Issue
Lifetimes of an Exomoon Orbiting a Jupiter-Like Planet in a Double Star System with the Mass of the Sun
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Assessing Parameters of the Coplanar Components of Perturbing Accelerations Using the Minimal Number of Optical Observations

1
Keldysh Institute of Applied Mathematics of Russian Academy of Science, 125047 Moscow, Russia
2
Astronomical Scientific Center JSC, 111123 Moscow, Russia
3
Peoples’ Friendship University of Russia (RUDN University), 117198 Moscow, Russia
*
Author to whom correspondence should be addressed.
Symmetry 2022, 14(12), 2564; https://doi.org/10.3390/sym14122564
Submission received: 7 October 2022 / Revised: 24 November 2022 / Accepted: 30 November 2022 / Published: 5 December 2022
(This article belongs to the Special Issue Advances in Mechanics and Control)

Abstract

:
The method presented in this paper is developed to assess the parameters (the application moment and the magnitude of a velocity impulse) of a maneuver-like perturbation of motion of the center of mass of a spacecraft in a near-circular orbit. The assessment is based on the information on the spacecraft’s trajectory before the maneuver and the optical observations of the spacecraft’s angular position (right ascension and declination angles) after the maneuver. This study considers the cases of solely transversal (in-track) or transversal and radial components of the velocity increment vector. A single pair of the values of angles is used for the assessment of the single transversal maneuver parameters and two pairs are used in the other cases. The method also makes it possible to estimate the parameters of a continuous maneuver performed with low-thrust engines. For this case the property of its symmetry is used. The approach described in this article makes it possible to determine the spacecraft’s orbit after the maneuver much faster and more accurately in comparison to traditional methods.

1. Introduction

Presently, the orbital information on approximately 26,000 space objects in near-Earth orbits is available from public sources accessible via the Internet. It is relatively easy to propagate motion for most of them since all considerable natural perturbing factors which define this motion are described with acceptable accuracy by the known motion models. However, there are space objects whose motion is heavily influenced by additional perturbations. These perturbations are hard to account for directly in the process of the determination and propagation of the parameters of motion of the space object’s center of mass. First of all, it concerns the maneuvering spacecraft. At present, there are more than five thousand such objects in orbit and their number is increasing fast due to very large spacecraft constellation deployments. These objects perform maneuvers to maintain their orbits or change them for the fulfillment of new tasks. As a rule, only the operators of these maneuvering spacecraft have a priori and actual data about the orbit parameter alterations. Furthermore, maneuvering spacecraft of some design have the ability to calculate their maneuvers on board. In the course of the space object catalogue maintenance process, one has to wait for the accumulation of the sufficient amount of observation data for the accurate determination of new orbits of these objects after the maneuver is performed.
The use of the observation information for the assessment of parameters of motion of the center of mass of maneuvering spacecraft, simultaneously with the assessment of parameters of the performed maneuver, allows speeding up the process of orbits updating in the space object catalogue. This catalogue is necessary for solving different tasks. For example, it allows to predict and assess the potentially hazardous close approaches between space objects with the maneuvering spacecraft being one of them.
Numerous modern approaches to the determination and assessment of the maneuver parameters of spacecraft are based on the improvement and modification of Kalman type filters, which were actively developed at the end of the last century for atmospheric aircraft [1,2]. Some of the first works to take into account the peculiarities of near-Earth space flight were [3,4]. The work [3] shows the feasibility of the application of a sequential filter for the procession of radiolocation measurements for the determination of long-duration and impulsive maneuver parameters. Woodburn et al. in [4] considers the maneuvers of small duration and uses the radiolocation measurements as preliminary ones. The solution is based on the application of the sequential filter for the assessment of the parameters of motion from the initial point after the moment of the impulse application to the moment of the last measurement. Then, the state vectors before and after the maneuver allow assessing the parameters of the impulsive maneuvers. As the approaches to the solution of the cataloging problem of the maneuvering objects developed, the methods allowing distinguishing the fact of maneuvering were divided into separate groups [5,6,7,8,9]. The methods of the maneuver parameter assessment can be divided into the following criteria: the processing of data on state vectors and covariance matrices exclusively [10,11] or the inclusion of trajectory measurements [3,4,12,13,14,15,16,17]—the usage of radiolocation measurements [12,13,14] or optical measurements [15,16,17]. The methods based on radar measurements are aimed to raise the resistance of filters to the occurrence of unaccounted accelerations. This can be achieved through the combination and modification of Kalman filters: extended Kalman filter (EKF), state smoothing, limited batch least-squares data reduction (BLSQ) in [13] and predictive corrective iterations for the EKF and unscented Kalman filter (UKF) in [14]. The work [14] describes the concept of the variable structure estimator (VSE), which also uses the EKF and belongs to the multiple model adaptive estimation (MMAE) family. The works [15,16,17] should be noted as they present the problem statement most close to the statement of this research. G. Escribano et al. [15] suggests a new method based on the presentation of the space of states as a stochastic hybrid system, which allows implementing the methods of successive Monte Carlo filtration, with which the maneuver parameters are determined during the process of Bayesian inference. The lack of dynamical constrictions and effectiveness of the problem solution from the fuel expenditures point of view can be treated as a disadvantage of this method. In the meantime, the work [16] contains an approach for narrowing the space of permissible solutions on the basis of the limitation of the finite energy expenditures. This approach is complemented with the algorithm of analyses of the previous spacecraft maneuvers which allows propagating the probability density of the alteration of the orbital parameters for the given moments of time after the last realized maneuver. K. Hill [17] suggests the two angle pairs initial orbit with conjunction analysis method (TAPIOCA), which consists of the preliminary assessment of the orbit on the basis of two angular measurements after the maneuver realization and successive search for the point of the closest approach between the maneuvering spacecraft before and after the maneuver. The obtained maneuver parameters are used in the least squares method for the refinement of the orbit after the maneuver. The inability of the assessment of the long-duration maneuver parameters during measurement gathering in the process of the dynamic operation fulfillment can be treated as a disadvantage of this method. The simultaneous maneuver assessment along with the orbit determination after the maneuver are also used in [18]. Although many methods have been proposed, they all come down to filtering in one form or another, which has significant computational costs.
The maneuvers of space objects in near-circular orbits are assessed in this paper and the preceding works of the authors. This simplifies the problem and allows suggesting semianalytical methods of its solution, which decreases the problem solution time substantially. The crucial importance of the performance arises from the necessity of regular operative assessments of maneuvers of a large number of real maneuvering spacecraft.
Especially important is the problem of the maneuver parameter assessment [19] when the orbit shaped by the maneuver is determined with the use of observation data with considerable uncertainties. These uncertainties arise due to the small number of used observations with low accuracy. An example of solving a similar problem using the method developed in [19] as applied to a real maneuvering spacecraft is provided in Table 1.
Table 1 contains the results of the assessment of a single-impulse coplanar maneuver of a spacecraft in a geostationary orbit with errors in the orbit determination [19]. The second column indicates the result of the traditional method of maneuver assessment (the velocity difference in the point of maximum proximity) obtained without the accountancy of orbit determination errors. The third one shows the assessment result obtained with the accountancy of the determination error. The fourth one provides the magnitude of the real implemented velocity impulse. In the case when the orbit after the maneuver was determined with the use of the observation information of the short observation interval with substantial uncertainties, the traditional method provides an assessment of the maneuver parameters which is far from the actual values. The method described in [19] allows increasing considerably the accuracy of the maneuver parameters assessment (the magnitudes and the application moments). Furthermore, one can obtain the terminal orbit considerably closer to the real orbit than the orbit obtained with the help of observations on the short observation interval if the estimated velocity impulse was applied to the initial orbit.
The use of the refined orbital parameters after performing the maneuver makes it possible to significantly increase the speed of the reliable recognition and evaluation of potentially dangerous encounters of space objects (one of them being a maneuvering spacecraft).
The next step is taken in this work. The maneuver parameters are assessed directly from the minimal number of observations without waiting for the sufficient number of observations needed for the initial approximate orbit determination to accumulate. The developed method of the maneuver assessment allows substantially speeding up the process of the acquisition of reliable and accurate enough assessments of the motion parameters of these objects. This result is possible thanks to the use of the minimal number of pairs α and δ (the angles of the right ascension and declination), which define the direction from the observer on the Earth’s surface to the points on the celestial sphere. At these points, the maneuvering spacecraft is detected at the corresponding moments of time during the optical monitoring session.
The suggested method of the maneuver assessment allows determining the parameters of the single-impulse maneuvers with different attitudes of the thrust vector in the case of the use of the high- and low-thrust engines.

2. Materials and Methods

2.1. Variants of the Problems Being Solved

2.1.1. Basic Problem Statement

Let the state vector of the maneuvering space object X = X(t0) be accurately known at the moment t0. The observers’ positions on the Earth’s surface Y = Y(ti) at the moment ti (ti > t0) and the angles αi and δi, which set the direction from the observer to the maneuvering space object at this moment, are known. It is necessary to determine the moment of application, the magnitude and the attitude of the velocity impulse for maneuvers close to impulses, and for maneuvers performed with a low-thrust engine, to evaluate the moments of switching the propulsion system on and off and the attitude of the thrust vector, as well as the acceleration caused by the operation of the propulsion system.

2.1.2. Equations for the Velocity Impulse Influence

Orbit parameters change instantaneously if the interval of the engines’ work is small with respect to the maneuvering space object’s orbit period.
The spacecraft moves along an orbit that is close to a circular orbit (geostationary orbit), which we call the reference orbit. This allows one to write down the equations of motion in deviations of the real orbit from a circular one in a cylindrical coordinate frame [20]. The solution of the linearized equations allows one to write down the influence of velocity impulses.
Each of these velocity impulses, applied at the point with the angles φi defined as an argument of latitude on the reference orbit at time ti (i = 1, …, N), causes the alteration of orbit elements at the point defined by the angle φf. The sum of these deviations caused by N velocity impulses can be written as [21]:
i = 1 N r 0 ( Δ V r i V 0 sin ( φ f φ i ) + 2 Δ V t i V 0 ( 1 cos ( φ f φ i ) ) ) = Δ r ,
i = 1 N ( Δ V r i cos ( φ f φ i ) + 2 Δ V t i sin ( φ f φ i ) ) = Δ V r ,
i = 1 N ( Δ V r i sin ( φ f φ i ) Δ V t i ( 1 2 cos ( φ f φ i ) ) ) = Δ V t ,
i = 1 N r 0 ( 2 Δ V r i V 0 ( 1 cos ( φ f φ i ) ) Δ V t i V 0 ( 3 ( φ f φ i ) 4 sin ( φ f φ i ) ) ) = Δ n ,
i = 1 N r 0 Δ V z i V 0 sin ( φ f φ i ) = z ,
i = 1 N Δ V z i cos ( φ f φ i ) = V z .
Here, r0 and V0 are the radius and velocity of the reference circular orbit in the vicinity of which the motion occurs; ΔVri, ΔVti and ΔVzi are the radial, transversal and lateral components of the i-th velocity impulse, correspondingly; Δr, ΔVr, ΔVt and Δn are the deviation along the radius vector, deviations of the radial and transversal components, and deviation along the orbit caused by the velocity impulses; and the angles φi and φf are measured from the maneuvering space object’s position at the moment t0 to the direction of the maneuvering space object’s motion.

2.1.3. Determination of Deviations Along the Radius and Along the Orbit Using Observations

In order to solve the stated problem, one should first calculate the deviations caused by the maneuver at the moment ti along the radius Δri and along the orbit Δni using the information about the known angles αi and δi at the moment ti.
For this purpose, at first, it is necessary to find the intersection point of the beam originating at the observer’s position on the Earth’s surface with the right ascension and declination angles αi and δi and the plane of the initial maneuvering space object’s orbit. This plane’s orientation can be determined with the help of numerical integration with the use of the initial conditions vector X(t0) at the moment ti. Note that the numerical propagation of the orbit from time t0 to ti should be performed with a fairly accurate model that incorporates major disturbances. Then, Δ r i = r p ( t i ) r o r b i   is calculated, where rp(ti) is the magnitude of the vector rp(ti), which is pointed out of the Earth’s center to the point of intersection (this vector lies in the plane of the initial maneuvering space object’s orbit) and rorbi is the magnitude of the radius vector of the initial orbit aimed along the vector rp(ti). Therefore, Δ r i indicates the difference in the radius vector magnitude of the real observed orbit, which exhibits the yet undefined impulse and propagates the initial orbit without the impulse. In addition, the distance along the orbit Δni between the position of the maneuvering space object in the initial orbit propagated at the moment ti (it is determined by the vector r(ti)) and its position in the orbit shaped by the maneuver at this moment (it is set by the vector rp(ti)) is calculated. To calculate Δni, the angle between the vectors r(ti) and rp(ti) is multiplied by r0 in a framework of linearized equations of relative motion.
Further on, we work with Δr and Δn directly. The technical part of the work described in this section requires the main efforts in writing programs.

2.1.4. Possible Variants of the Problems Being Solved

As the different variants of maneuver realizations are possible, there are several problems that need to be solved. The methods of solving these problems and the number of Δr and Δn (α and δ) pairs used also differ. Only the problem statements relating to coplanar maneuvers are provided.
These problems are:
  • A transversal velocity impulse is performed. Its application angle φ and the magnitude of the transversal component ΔVt should be determined. There are two unknown variables in the problem, hence, one pair of Δr and Δn (α and δ) is enough for its solution. It is the most frequently met problem;
  • The same as above but the impulse is radial. The problem in this statement is solved in some rare special cases;
  • A velocity impulse with the transversal and radial components is performed. Its application angle φ and the magnitudes of the transversal and radial components ΔVt and ΔVr should be determined. There are three unknown variables in the problem, hence, two pairs of deviations Δr and Δn are needed for its solution;
  • A long-duration transversal maneuver is performed. The angular duration of the maneuver Δφ, the angle which corresponds to its medium point φm and the constant acceleration of the maneuvering space object during the maneuver wt should be found. There are three unknown variables in the problem, hence, two pairs of deviations Δr and Δn are needed for its solution;
  • A long-duration maneuver with the transversal and radial components is performed. The angular duration of the maneuver Δφ, the angle which corresponds to its medium point φm and the transversal wt and radial wr constant accelerations are to be derived. There are four unknown variables in the problem, hence, two pairs of deviations Δr and Δn (α and δ) are needed for its solution.
There are only three unknown variables in the third and the fourth problems. Two pairs of Δr and Δn are used which provide four equations using Equations (1) and (4). Three equations out of four are directly used to find problem parameters, while the fourth equation can be used for a selection of a more accurate solution. This is discussed in detail when solving specific problems.

2.2. Small-Duration Maneuver Parameters Assessment

2.2.1. Transversal and Radial Velocity Impulse Assessment

The equations for the determination of the velocity transversal impulse application angle φ and the magnitude ΔVt can be described as:
2 r 0 Δ V t V 0 ( 1 cos ( φ f φ ) ) = Δ r , r 0 Δ V t V 0 ( 3 ( φ f φ ) 4 sin ( φ f φ ) ) = Δ n .
Here, the angle φf corresponds to the point in which the deviations Δr and Δn are calculated.
Figure 1 shows a change in the shape of the orbit, accompanied by a deviation along the radius caused by the transversal component of the velocity impulse. Figure 2 shows the deviation along the orbit depending on the angular distance from the moment of application of the velocity impulse.
The deviations along the radius Δr and along the orbit Δn = Δu · r0 shown in Figure 1 and Figure 2 can be calculated in the program using the angles αi and δi. Then, they themselves are used to determine the maneuver parameters.
The parameter ΔVt in Equation (7) is determined from the first equation and is used in the second one. This leads to Equation (8) for the determination of φ, which is in turn used in Equation (7) so ΔVt can be determined.
Δ V t = V 0 Δ r 2 r 0 ( 1 cos ( φ f φ ) )   ,
Δ r 3 ( φ f φ ) 4 sin ( φ f φ ) 2 2 cos ( φ f φ ) = Δ n .
This is the simplest problem in which the unique solution can be found. The accuracy of the deviation calculations Δr and Δn effect the solution quality.
Similarly, one can obtain formulas for determining the parameters of the radial velocity impulse.
Δ V r = V 0 Δ r r 0 sin ( φ f φ )   ,
Δ r ( 1 cos ( φ f φ ) ) sin ( φ f φ ) = Δ n .
This is also a relatively easy problem, in which the unique solution can also be found fast. The radial velocity impulses can rarely be met, but the possibility of obtaining the parameters of such a maneuver should be provided.

2.2.2. Assessment of the Velocity Impulse with Transversal and Radial Components

There are three unknown variables in the problem: the application angle φ and the magnitudes of the transversal and radial components ΔVt and ΔVr, hence, two pairs of deviations Δr and Δn calculated at the given moments φ1 and φ2 should be used for its solution.
The equations for the determination of the velocity impulse application angle φ and the magnitudes of the transversal and radial components of the velocity impulse ΔVt and ΔVr are:
r 0 ( Δ V r V 0 sin ( φ 1 φ ) + 2 Δ V t V 0 ( 1 cos ( φ 1 φ ) ) ) = Δ r 1 ,
r 0 ( Δ V r V 0 sin ( φ 2 φ ) + 2 Δ V t V 0 ( 1 cos ( φ 2 φ ) ) ) = Δ r 2 ,
r 0 ( 2 Δ V r V 0 ( 1 cos ( φ 1 φ ) ) Δ V t V 0 ( 3 ( φ 1 φ ) 4 sin ( φ 1 φ ) ) ) = Δ n 1 ,
r 0 ( 2 Δ V r V 0 ( 1 cos ( φ 2 φ ) ) Δ V t V 0 ( 3 ( φ 2 φ ) 4 sin ( φ 2 φ ) ) ) = Δ n 2 .
We determine ΔVt in Equation (15) and ΔVr in Equation (16) from Equations (11) and (12) and use them in Equation (14). Equation (14) is preferred over Equation (13) because it has a greater deviation of Δn than in Equation (13), which reduces the effect of the error in the determination of the value of Δn. We obtain Equation (17) for determining φ. Having determined φ, we use it in Equations (15) and (16) and determine ΔVt and ΔVr.
Δ V t = V 0 ( Δ r 1 sin ( φ 2 φ ) Δ r 2 sin ( φ 1 φ ) ) 2 r 0 ( ( 1 cos ( φ 1 φ ) ) sin ( φ 2 φ ) ( 1 cos ( φ 2 φ ) ) sin ( φ 1 φ ) )   ,
Δ V r = V 0 ( Δ r 1 ( 1 cos ( φ 2 φ ) ) Δ r 2 ( 1 cos ( φ 1 φ ) ) ) r 0 ( sin ( φ 1 φ ) ( 1 cos ( φ 2 φ ) ) sin ( φ 2 φ ) ( 1 cos ( φ 1 φ ) ) )   ,
Δ n 2 = 2 ( Δ r 1 ( 1 cos ( φ 2 φ ) ) Δ r 2 ( 1 cos ( φ 1 φ ) ) ) r 0 ( sin ( φ 1 φ ) ( 1 cos ( φ 2 φ ) sin ( φ 2 φ ) ( 1 cos ( φ 1 φ ) ) ( 1 cos ( φ 2 φ ) ) ( Δ r 1 sin ( φ 2 φ ) Δ r 2 sin ( φ 1 φ ) ) ( 3 ( φ 2 φ ) 4 sin ( φ 2 φ ) ) 2 ( ( 1 cos ( φ 1 φ ) ) sin ( φ 2 φ ) ( 1 cos ( φ 2 φ ) ) sin ( φ 1 φ ) ) ) .
Equation (17) sometimes has several solutions for this problem. In this case, the velocity impulse parameters for each solution are used in Equation (13). The solution for which the right and the left parts of Equation (13) match with higher accuracy is taken as the final solution of the problem.
By applying the calculated velocity impulse to the initial orbit, we can obtain the orbit after the maneuver. This orbit can already be used to calculate the rendezvous with this object and for other purposes.

2.3. Long-Duration Maneuver Parameters Assessment

For determination of the influence of transversal and radial thrust vector components on deviations along the radius vector and along the orbit, it is assumed that the thrust and the acceleration created by it are constant over the entire interval of the engine operation time.
First, let us find the influence of the maneuver performed with the low-thrust engine, directed along the transversal component in the orbital coordinate frame, on the deviation along the radius vector at the moment set by the angle φf.
It is assumed that the maneuver has the angular duration Δφ, the middle of the maneuver is set by the angle φm and the constant acceleration of the maneuvering space object within the maneuver is wt. The key is the assumption that the maneuver is symmetric about its center. By using Equation (1), along with the assumption of maneuver symmetry, one can find the deviation along the radius produced after the ΔVt application (ΔVt = wtΔt and Δt is the operating time of the propulsion system), which is equally distributed on the interval of the latitude argument Δφ:
Δ r t = 2 r 0 V 0 Δ V t Δ ϕ Δ φ / 2 Δ φ / 2 ( 1 cos ( φ f φ m φ ) ) d φ = 2 r 0 V 0 Δ V t Δ φ ( Δ φ 2 sin Δ φ 2 cos ( φ f φ m ) ) .
By using the equation Δ φ = λ 0 Δ t = k Δ V V 0 = w c w Δ V V 0 , where λ 0 = V 0 r 0 ,   k = m V 0 2 P   r 0 = w c w , wc is the centripetal acceleration of the reference circular orbit ( w c = V 0 2 r 0 ) and w t is the acceleration produced by the spacecraft’s engine ( w = P m ), we obtain the final formula for the deviation:
Δ r t = 2 r 0 w t w c ( Δ φ 2 sin Δ φ 2 cos ( φ f φ m ) ) .
Similarly, the influence of the radial component distributed on the interval of the latitude argument Δφ on the deviation along the radius can be found:
Δ r r = r 0 V 0 Δ V r Δ ϕ Δ φ / 2 Δ φ / 2 sin ( φ f φ m φ ) d φ = 2 r 0 V 0 Δ V r Δ φ sin Δ φ 2 sin ( φ f φ m ) ,
Δ r r = 2 r 0 w r w c sin Δ φ 2 sin ( φ f φ m ) .
By using Equation (4), we can find the influence of the maneuver performed with the low-thrust engine, the orientation of which is fixed by the transversal component in the orbital coordinate frame, on the deviation along the orbit. Just as before it is assumed that the maneuver has the angular duration Δ φ , the middle of the maneuver is determined by the angle φm:
Δ n t = r 0 V 0 Δ V t Δ ϕ Δ φ / 2 Δ φ / 2 ( 3 ( φ f φ m φ ) 4 sin ( φ f φ m φ ) ) d φ = r 0 V 0 Δ V t Δ ϕ ( 3 ( φ f φ m ) Δ φ 8 sin Δ φ 2 sin ( φ f φ m ) ) ,
Δ n t = r 0 w t w c ( 3 ( φ f φ m ) Δ φ 8 sin Δ φ 2 sin ( φ f φ m ) ) .
Similarly, the influence of the radial component distributed on the interval of the latitude argument Δφ on the deviation along the orbit can be found:
Δ n r = 2 r 0 V 0 Δ V r Δ ϕ Δ φ 2 Δ φ 2 ( 1 cos ( φ f φ m φ ) ) d φ = 2 r 0 V 0 Δ V r Δ φ ( Δ φ + 2 sin Δ φ 2 cos ( φ f φ m ) ) ,
Δ n r = 2 r 0 w r w c ( Δ φ + 2 sin Δ φ 2 cos ( φ f φ m ) ) ,
where φf is the point in which the deviations are calculated.
In the problem of determining the parameters of the maneuver performed with a low-thrust engine, there are four unknown variables φm, Δφ, wt and wr, hence, it is necessary to use two pairs of deviations Δr and Δn (α and δ) for its solution. By solving the equation system of four equations, in which the influence of the engines’ work is determined by Equations (18)–(21), we obtain the values of all unknowns. If it is known that the object uses only transversal orientation while performing the maneuver, there are only three unknown variables φm, Δφ and wt. It is enough to use three deviations: Δn for the first point and Δr and Δn for the second one. The deviations along the orbit are always more preferable to use, as they can be determined in a more accurate fashion, also it is more preferable to use the deviations which correspond to the more distant point from the maneuver.
By changing the elements of the orbit at the second point by the magnitude of the calculated deviations [21], we obtain an orbit taking into account the influence of the long-term operation of the propulsion system.

3. Results

The examples with the known realized velocity impulses from Table 2 were taken in order to obtain the results which allowed us to assess the accuracy of the found solution.
Optical measurements were performed with errors, which, for a considered mission, did not exceed five arcseconds. It could be seen that due to these errors, the deviations measured at different time points differed, although the measurement times were very close to each other.
All observations were conducted within an interval of less than two minutes, therefore, ten observations were replaced by one average measurement.
We obtained the mean values Δr = −5927 (km) and Δn = 13,761 (km).
For comparison purpose, we calculated using Equations (1)–(4) the deviations caused by the real velocity impulse ΔVt = −0.112 m/s applied at Δφ = 178.381° to be Δr = −6144 (km) and Δn = 14,176 (km).
It can be seen that the averaged deviations calculated from Table 3 prove to be quite close to the theoretical ones.
The solution using the averages of Δr and Δn was Δφ = −178.858 and ΔVt = −0.108 m/s. The derived maneuver magnitude and application angle are very close to the real maneuver parameters in Table 2.
Table 4 provides parameters of the second example. In this case the maneuver was not exactly an impulsive one. Instead, it took approximately 8 min. However, the maneuver could be treated as an impulsive one due to the fact that the maneuver performance duration was much shorter in comparison with the orbit period.
The run contained ten observations, the duration of the run was 40 s, the average time of the run was 23 June 2021 20:57:52.89, the average angle was 496° (136°) and the mean values of the deviations were Δr = 0.521 km and Δn = −65.623 km.
The theoretical values of the deviations for the angle between the moment of application of the real velocity impulse and the moment of the optic observations φ = 325.031° and ΔVt = 0.246 m/s were Δr = 1218 km and Δn = −65,130 km.
The found solution was the velocity impulse application angle of 160.49° (φ = −345.84° from the observed point) and the velocity impulse magnitude ΔVt = 0.249 m/s.
The maneuver magnitude was determined with high accuracy in both examples, even for the second one with the continuous, although relatively short, maneuver. This meant that the semimajor axis of the orbit after the maneuver was now known practically with the same accuracy as the accuracy with which the semimajor axis of the initial orbit was determined. The error in the measurement of the semimajor axis was approximately Δa = 100 m and the error in the measurement of the eccentricity was Δe = 0.000002.
It can be noted that the deviations determined by the observations and the deviations determined analytically differ from each other especially in the second example. It is mainly related to the deviations along the radius. This is the main problem in the use of this method connected with the present errors of the observations. The given approach cannot be used when the orbit observations are situated closely to the moment of the maneuver performance. In this case the deviations caused by the velocity impulse were less than the observation errors themselves.
In the second example, the angle between the moment of the velocity impulse application and the average moment of the observations was substantial enough and comprised 325°. The considerable deviation along the orbit (the theoretical deviation was Δn = −65,130 km and the value measured with the observations was Δn = −65,623 km) was accumulated for this period of time. The 0.5 km error was not crucial compared to the deviation itself of 66.5 km. However, the error along the radius vector of magnitude 0.7 km was substantial (the theoretical deviation was Δr = 1218 km and the deviation calculated with the measurements was Δr = 0.521 km), which was commensurable with the determined value. This error caused the error in the determination of the velocity impulse application moment. Furthermore, if we just took the average value of all 10 observations from the orbit observation run, we obtained the value Δr = 0.214 km. However, the considerable negative deviation along the orbit indicated that the positive velocity impulse was performed. Hence, all deviations along the radius were positive. This allowed the removal of four observations with negative deviations. Figure 1 illustrates this discussion, as a positive transversal impulse led to the radius vector being larger than the initial radius vector at any point on the orbit. As a result, the mean deviation increased from Δr = 0.214 km to Δr = 0.521 km. While using the deviation along the orbit for the correction of the deviation along the radius vector, one should bear in mind that the arousing (after the positive transversal velocity impulse) shift along the orbit was positive at the start and then reached its maximum at Δ n 0.4776 r 0 V 0 Δ V t with φ = 41°24′35″ (Figure 2) [22], then the shift started to decrease and after φ = 73°05′32″ the retardation from the unperturbed motion started. The retardation per one revolution was Δ n = 6 π r 0 V 0 Δ V t . This possibility of the alteration of the sign of the deviation along the orbit needed to be accounted for. The number of deviations along the radius with different signs, as well as the magnitude and the sign of the deviation along the orbit itself were considered for the removal of the erroneous deviations along the radius.
Avoiding the use of the radial deviations was the simplest way to remove errors in the radial direction. In this case, two deviations along the orbit Δn were used. The solution for this option for Example 1 is provided below.
The real velocity impulse parameters derived from Table 2 were ΔVt = −0.112 m/s and φ = 235.92°.
The found solution was ΔVt = −0.112 m/s and φ = 235.69°.
The magnitudes of the velocity impulse coincided and the angles of application of the velocity impulse almost coincided as well. The angles φ of the velocity application differ from those provided in Table 1 (φ = 178.38°) because they are measured from the second compact group of observations additionally used in this example.
The new solution confirms the one observation series solution. The accuracy of the solution increased due to the fact that the used deviation along the orbit increased and the error in determining this value remained within the same limits. However, the solution was obtained approximately six hours later. However, even in this case, we obtained the new orbit significantly earlier in comparison to the use of the usual technique of determining the orbit with observations after the maneuver and with a higher accuracy.

4. Conclusions

In the geostationary orbit, a quite frequent situation occurs when several spacecraft which maintain their orbits via maneuvers operate at almost the same point. It is essential to quickly and accurately predict the motion of each spacecraft and its neighbors in order to firmly protect each spacecraft from collisions. A similar problem arises when the spacecraft is transferred to the operating point.
This paper describes the method that allows the assessment of the maneuvers performed in the plane of a geostationary orbit using no more than two compact measurements of the right ascension and declination angles. The provided examples demonstrate high accuracy of this method. The traditional technique usually uses four spaced revolutions in order to determine the orbit after the maneuver with tolerable accuracy. Thus, the method from this work helps to reduce the time of the spacecraft’s orbit determination after the maneuver by at least two times, while the new orbit is determined with almost the same high accuracy with which it is determined before the maneuver.
Unlike traditional methods in which the orbit after the maneuver is first determined and then the maneuver parameters are estimated, in this work, the maneuver parameters are estimated directly from the measurements, and then using this information, the parameters of the formed orbit are determined. This increases the speed of determining the parameters of the formed orbit and the accuracy of determining its parameters.
Maneuver evaluation provides a double effect. One can identify the possibility of a collision with this object along with the determination of the purpose of the maneuver (maintenance of the orbit, transfer to a new operating point or transfer to a disposal orbit). The evaluation of the maneuvers of your own spacecraft makes it possible to determine the health of the propulsion system and whether the implemented maneuver matches to the planned one.
It is also possible to determine the velocity impulse caused by an imbalance in the operation of the orientation motors. As a result of this imbalance, along with the given rotation of the spacecraft, there is an undesirable movement of its center of mass. The evaluation of the emerging velocity impulse allows predicting the undesirable movement of the center of mass.

Author Contributions

Conceptualization, A.B. and V.A.; methodology, A.B.; software, N.G. and M.K.; validation, A.B. and V.A.; investigation, N.G. and M.K.; writing—original draft preparation, A.B.; writing—review and editing, V.A. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Moscow Center of Fundamental and Applied Mathematics, in agreement with the Ministry of Science and Higher Education of the Russian Federation, No. 075-15-2022-283.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mcintyre, G.; Hintz, K. Comparison of several maneuvering target tracking models. In Proceedings of the SPIE—The International Society for Optical Engineering 3374, Orlando, FL, USA, 13–15 April 1998. [Google Scholar] [CrossRef]
  2. Li, X.R.; Jilkov, V.P. Survey of Maneuvering Target Tracking. Part I. Dynamic Models. IEEE Trans. Aerosp. Electron. Syst. 2003, 39, 1333–1364. [Google Scholar] [CrossRef]
  3. Hujsak, R.S. Orbit Determination During High Thrust and Low Thrust Maneuvers; AAS 05-1 36; Analytical Graphics Inc.: Exton, PA, USA, 2002. [Google Scholar]
  4. Woodburn, J.; Carrico, J.; Wright, J.R. Estimation of Instantaneous Maneuvers Using a Fixed Interval Smoother. Adv. Astronaut. Sci. 2003, 116, 243–260. [Google Scholar]
  5. Bai, X.; Liao, C.; Pan, X.; Xu, M. Mining Two-Line Element Data to Detect Orbital Maneuver for Satellite. IEEE Access 2019, 7, 129537–129550. [Google Scholar] [CrossRef]
  6. Mukundan, A.; Wang, H.-C. Simplified Approach to Detect Satellite Maneuvers Using TLE Data and Simplified Perturbation Model Utilizing Orbital Element Variation. Appl. Sci. 2021, 11, 10181. [Google Scholar] [CrossRef]
  7. Clark, R.; Lee, R. Parallel processing for orbital maneuver detection. Adv. Space Res. 2020, 66, 444–449. [Google Scholar] [CrossRef]
  8. Huang, J.; Hu, W.; Zhang, L. Maneuver detection of space object for space surveillance. In Proceedings of the 6th European Conference on Space Debris, Darmstadt, Germany, 22–25 April 2013. [Google Scholar]
  9. Folcik, Z.; Cefola, P.; Abbot, R. GEO maneuver detection for space situational awareness (AAS 07-285). Adv. Astronaut. Sci. 2008, 129, 523. [Google Scholar]
  10. Serra, R.; Yanez, C.; Frueh, C. Tracklet-to-orbit association for maneuvering space objects using optimal control theory. Acta Astronaut. 2021, 181, 271–281. [Google Scholar] [CrossRef]
  11. Jaunzemis, A.D.; Mathew, M.V.; Holzinger, M.J. Control cost and Mahalanobis distance binary hypothesis testing for spacecraft maneuver detection. J. Guid. Control Dyn. 2016, 39, 2058–2072. [Google Scholar] [CrossRef] [Green Version]
  12. Guang, Z.; Xingzi, B.; Hanyu, Z.; Bin, L. Non-cooperative maneuvering spacecraft tracking via a variable structure estimator. Aerosp. Sci. Technol. 2018, 79, 352–363. [Google Scholar] [CrossRef]
  13. Kelecy, T.; Jah, M. Detection and orbit determination of a satellite executing low thrust maneuvers. Acta Astronaut. 2010, 66, 79. [Google Scholar] [CrossRef]
  14. Bergmann, C.M.; Zollo, A.; Herzog, J.; Schildknecht, T. Integrated Manoeuvre Detection and Estimation Using Nonlinear Kalman Filters During Orbit Determination of Satellites. In Proceedings of the 8th European Conference on Space Debris, Darmstadt, Germany, 20–23 April 2021. [Google Scholar]
  15. Escribano, G.; Sanjurjo-Rivo, M.; Siminski, J.A.; Pastor, A.; Escobar, D. Automatic maneuver detection and tracking of space objects in optical survey scenarios based on stochastic hybrid systems formulation. arXiv 2021, arXiv:2109.07801. [Google Scholar] [CrossRef]
  16. Siminski, J.; Fiedler, H.; Flohrer, T. Correlation of Observations and Orbit Recovery Considering Maneuvers. In Proceedings of the AAS/AIAA Space Flight Mechanics, San Antonio, TX, USA, 5–9 February 2017. [Google Scholar]
  17. Hill, K. Maneuver detection and estimation with optical tracklets. In Proceedings of the Advanced Maui Optical and Space Surveillance Technologies Conference, Maui, HI, USA, 9–12 September 2014; p. E26. [Google Scholar]
  18. Kamensky, S.; Tuchin, A.; Stepanyants, V.; Alfriend, K.T. Algorithm of Automatic Detection and Analysis of non-Evolutionary Changes in Orbital Motion of Geocentric Objects. In Proceedings of the AAS/AIAA Astrodynamics Specialist Conference, Pittsburgh, PA, USA, 9–13 August 2009; Paper AAS 09-103. pp. 1–34. [Google Scholar]
  19. Baranov, A.A.; Vikhrachev, V.O.; Karatunov, M.O.; Razumnyi, Y.N. Estimates of Near-Circular Orbits after a Single Correction: A Geometrical Method. J. Comput. Syst. Sci. Int. 2017, 56, 137–145. [Google Scholar] [CrossRef]
  20. Baranov, A. Spacecraft Manoeuvring in the Vicinity of a Near-Circular Orbit; Cambridge Scholars Publishing: Cambridge, UK, 2022; 520p. [Google Scholar]
  21. Baranov, A.A.; De Prado, A.F.; Razumny, V.Y.; Baranov, A.A. Optimal Low-thrust Transfers Between Close Near-Circular Coplanar Orbits. Cosm. Res. 2011, 3, 269–279. [Google Scholar] [CrossRef]
  22. El’yasberg, P.E. Introduction to the Theory of Flight of Artificial Earth Satellites; Gurevich, U.G., Ed.; Israel Program of Scientific Translation: Jerusalem, Israel, 1967; 357p. [Google Scholar]
Figure 1. Orbit shape change after the application of the transversal component of the velocity impulse.
Figure 1. Orbit shape change after the application of the transversal component of the velocity impulse.
Symmetry 14 02564 g001
Figure 2. Along-the-orbit shift after the application of the transversal component of the velocity impulse.
Figure 2. Along-the-orbit shift after the application of the transversal component of the velocity impulse.
Symmetry 14 02564 g002
Table 1. Assessment of one single-impulse coplanar maneuver with the presence of errors in orbit determination.
Table 1. Assessment of one single-impulse coplanar maneuver with the presence of errors in orbit determination.
Maneuver
Parameters
Traditional Method With the Uncertainty AccountancyReference
(Actual Values)
Vt (m/s) −0.3920.4190.419
Maneuver application time18 h 7 min 33 s18 h 8 min 17 s18 h 9 min 54 s
Table 2. Example 1: Real maneuver parameters.
Table 2. Example 1: Real maneuver parameters.
ParameterValue
Time of the initial conditions settingt0 = 2022.04.11 16:37:08.950 (GMT+3)
Time of the maneuver2022.04.12 08:18:25.000
Maneuver magnitudeΔVt = −0.112 m/s
The angle between the maneuver
and the observation
Δφ = 178.381°
Table 3. Example 1: Series of very close measurements of the position of an object.
Table 3. Example 1: Series of very close measurements of the position of an object.
DateTimeΔrΔn
12 April 202220:08:18.543−6.0597513.73975
12 April 202220:08:29.543−5.5984713.83167
12 April 202220:08:40.543−5.6719713.78066
12 April 202220:08:51.543−5.8114213.64288
12 April 202220:09:02.543−5.9399713.86445
12 April 202220:09:13.543−6.3060513.82325
12 April 202220:09:24.543−5.6185713.73580
12 April 202220:09:35.543−6.1495413.73975
12 April 202220:09:46.543−6.1088213.66814
12 April 202220:09:57.543−6.0146213.79292
Table 4. Example 2: Real maneuver parameters.
Table 4. Example 2: Real maneuver parameters.
ParameterValue
Performance date 22 June 2021
Time of the maneuver start23:17:30
Duration of the burn465.5 s
Velocity impulse magnitude0.246 m/s
Time of the impulsive maneuver (the middle of the burn duration)23:21:22.75
Maneuver application angle (from the moment of the initial conditions setting)170.75°
Acceleration0.000529 m/s2
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Baranov, A.; Agapov, V.; Golosova, N.; Karatunov, M. Assessing Parameters of the Coplanar Components of Perturbing Accelerations Using the Minimal Number of Optical Observations. Symmetry 2022, 14, 2564. https://doi.org/10.3390/sym14122564

AMA Style

Baranov A, Agapov V, Golosova N, Karatunov M. Assessing Parameters of the Coplanar Components of Perturbing Accelerations Using the Minimal Number of Optical Observations. Symmetry. 2022; 14(12):2564. https://doi.org/10.3390/sym14122564

Chicago/Turabian Style

Baranov, Andrey, Vladimir Agapov, Natalya Golosova, and Maksim Karatunov. 2022. "Assessing Parameters of the Coplanar Components of Perturbing Accelerations Using the Minimal Number of Optical Observations" Symmetry 14, no. 12: 2564. https://doi.org/10.3390/sym14122564

APA Style

Baranov, A., Agapov, V., Golosova, N., & Karatunov, M. (2022). Assessing Parameters of the Coplanar Components of Perturbing Accelerations Using the Minimal Number of Optical Observations. Symmetry, 14(12), 2564. https://doi.org/10.3390/sym14122564

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop