International Journal of Aerospace Engineering

Volume 2012 (2012), Article ID 480320, 9 pages

http://dx.doi.org/10.1155/2012/480320

## Solving Coplanar Power-Limited Orbit Transfer Problem by Primer Vector Approximation Method

Department of Mechanical and Aerospace Engineering, University of Missouri, Columbia, MO 65211, USA

Received 1 November 2011; Revised 12 January 2012; Accepted 16 January 2012

Academic Editor: Alessandro A. Quarta

Copyright © 2012 Weijun Huang. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

The coplanar orbit transfer problem has been an important topic in astrodynamics since the beginning of the space era. Though many approximate solutions for power-limited orbit transfer problem have been developed, most of them rely on simplifications of the dynamics of the problem. This paper proposes a new approximation method called primer vector approximation method to solve the classic power-limited orbit transfer problem. This method makes no simplification on the dynamics, but instead approximates the optimal primer-vector function. With this method, this paper derives two approximate solutions for the power-limited orbit transfer problem. Numerical tests show the robustness and accuracy of the approximations.

#### 1. Introduction

Most trajectory optimization problems are nonlinear problems with no analytic solutions. However, to the coplanar power-limited orbit transfer in the classical inverse-square gravity field, many researchers have proposed approximate solutions, for example, Edelbaum [1–4], Zee [5], Marinescu [6], Marec and Vinh [7], Haissig et al. [8], Kechichian’s [9], and Casalino and Colasurdo’s [10]. These proposed solutions are built on assumptions about the transfer scenarios. For example, the solution of [6] assumes the transfer is in a close range, the solutions of [2, 4–6] assume the transfer happens in a long duration, and the solutions of [1, 3, 7–11] assume the admissible control to be within a rather limited set.

An approximate solution to an optimal control problem implies both the approximation of the control policy and the approximation of the dynamics. However, we automatically approximate the optimal control policy as well, when we approximate the dynamics. The mentioned references more or less approximate the dynamics. One potential problem with approximating dynamics is that, once the assumptions are violated the obtained control will be infeasible. From a software point of view, an infeasible control produces unexpected results and might cause the software to crash. Therefore, a more robust approximation method for optimal control problem should be to purely make approximations on the optimal control, not the dynamics.

A popular method for generating feasible transfer trajectories is the shape-based method [12–15]. In its essence, this method generates feasible control without any compromise on the dynamics. Though there are different variates of the shape-based method, as far as the author knows, there is no theoretical research to address the connection between the real optimal solution and the trajectory generated by the shape-based method.

This paper proposes an innovative approximation method—the primer vector approximation method, which combines the advantages of the approximation method and the shape-based method and uses feasible control to approximate the optimal control. The method reformulates the classic transfer problem with a nonlinear transformation from Carter and Humi [16]. The purpose of this reformulation is to put all the nonlinear terms to the coefficients of the control variables. Thus the optimal control vector, called primer vector, of the new formulation can be analyzed and approximated without affecting the dynamics. To demonstrate this method, this paper derives two approximate solutions for both low-thrust close-range transfer and low-thrust long-duration transfer.

Four transfer scenarios are designed to test the two approximate solutions numerically. The tests show that the approximations are close to the optimal solutions when the scenarios are within the assumptions of the approximations. More importantly, the solutions are feasible even when the scenarios are far away from the assumptions of the approximations. Theoretically, both approximations can generate feasible transfer control between any two kinds of orbits, including hyperbolic and parabolic orbit.

The remaining part of this paper is organized in the following way. The second section is about the formulations of the coplanar power-limited orbit transfer problem. The third section introduces the primer vector approximation method and derives two approximate solutions. The fourth section numerically tests the two solutions. A conclusion is made in the fifth section.

#### 2. Formulations of Coplanar Power-Limited Orbit Transfer Problem

##### 2.1. Polar Coordinate Formulation

The polar coordinate system used in this paper is presented in Figure 1. The meanings of the symbols are given in Nomenclature.

Mathematically, the dynamics of a controlled satellite in an inverse-square gravity field can be described as

Set and . The two second-order ordinary differential equations (ODEs), (1), are equivalent to the following four first-order ODEs:

Set . If a transfer start from time to time , the objective function of the coplanar power-limited orbit transfer problem is
The initial orbit at is defined by a vector and the final orbit at is defined by a vector . A *free-time rendezvous problem* is to find the optimal control function , which minimizes the objective function, (3), and transfers a spacecraft from the original position at time to the final position at the optimal arriving time . This problem is the simplest problem to introduce the primer vector approximation method. A more common power-limited transfer problem unbounded thrust is the *fixed-time rendezvous problem*, the solution to which can be found by adding a time constraint on the approximate solution to the free-time rendezvous problem. We can adopt the method in Novak’s [15] paper to use a solution to the free-time rendezvous problem as a basis for the solution to the related fixed-time rendezvous problem. This paper, however, focuses only on the free-time rendezvous problem.

To numerically solve the free-time rendezvous problem, we can use the indirect shooting method, which first transforms the problem into a two-point boundary value problem (TPBVP) and then solves the TPBVP with the shooting technique. To setup the TPBVP, we write down the augmented Hamiltonian of the optimal problem as follows:

The adjoint ODEs of the problem are

The optimal control is , where and is called primer vector. A common shooting method will shoot the optimal arrival time and the values of the four variables , , , and until the shooting functions and transversality condition are satisfied. At each iteration of the shooting process, there are eight ODEs, (2) and (5), to solve. The shooting process will be time-consuming if the initial guess is poor or the number of the transfer revolutions is large.

In this paper, the indirect shooting method with the polar coordinate formulation will be used to generate optimal solutions for the four tests in the fourth session.

##### 2.2. A New Formulation

The polar formulation of the power-limited orbit transfer problem has nonlinear terms in the dynamic equations and adjoint equations. A new formulation of the power-limited orbit transfer problem will be derived here with a nonlinear transformation from Carter and Humi [16]. This transformation will transfer the state variables with respect to time to new state variables with respect to polar angle .

Use and to represent the first and second derivative operations with respective to the polar angle . Set . Thus and . The new state variables are defined as , , and . The mapping from the new variables back to the variables in the polar coordinate is

From the chain rule, we have

Therefore, with (7) and (1), the dynamic equations of the new variables are

For further simplification, write (8) in matrix form, set and define

Thus, state space representation (SSR) of the ODEs in (8) is

The new boundaries of a transfer in the new dynamic equations will be defined by the initial polar angle , the final polar angle and two constant vectors of the new states

Under the new state vector and new control vector , the objective function, (3), becomes

The objective function, (12), dynamic equations, (10) and the boundary conditions, (11), compose the new formulation of the coplanar power-limited orbit transfer problem.

Mathematically, a coplanar power-limited orbit transfer problem can be described as

To solve this problem, set to be the new augmented Hamiltonian and to be the new adjoint vector corresponding to the new state variable

Apply the theory of optimal control and the optimal control of the new formulation becomes , in which is the primer vector of the new formulation and

The mapping from the new primer vector to the original primer vector is

Then the adjoint equations are

In conclusion, the transfer problem given by (13) is an alternative formulation of the coplanar power-limited orbit transfer problem. The optimal solution for this new formulation needs to satisfy a TPBVP defined by the dynamic equation, (10), adjoint equation, (17), the boundary conditions, (11).

##### 2.3. Properties of the New Formulation

The new formulation has unique properties that facilitate analysis. Firstly, the system can be solved explicitly when there is no control () on the state equation, (10). The solutions are

Secondly, the nonlinear term in the adjoint vector equation (17) has the following interesting property.

Theorem 1. *For an optimal trajectory, the solution for the adjoint vector has the following form:
**
where is the initial polar angle of the trajectory, is a 3 by 3 unit matrix, and is a diagonal matrix function.*

*Proof. *Since (17) is a linear nonhomogeneous equation with constant coefficients, its solution is the sum of the general solution for the related homogeneous equation and the particular integral. The solution can be written down as
Set
and (use index notation for )
With (20), (21), and (22), we obtain (19) of Theorem 1.

Up to this point, no assumption has been made about the new formulation of the transfer problem. It is, however, easy to see that, if the thrust level is very small, will become very small. Therefore, for a low-thrust transfer, the adjoint vector can be approximated by setting .

#### 3. Primer Vector Approximation Method

##### 3.1. Equivalent Optimal Control Problem

The primer vector approximation method aims to find an approximation for the optimal control (primer vector function), (15), without affecting the accuracy of the state dynamics. Thus, it is important to ensure that the approximate control will lead to a precise integration of the state equation, (10). This is the core of the primer vector approximation method.

To introduce the strategy for primer vector approximation, first define a family of simpler optimal control problems where is a set of matrix functions and all of its members are continuous while . The following theorem gives the relationship between the family and the optimal control problem, defined in (13).

Theorem 2. *Given a transfer problem , defined by (13), among the family of optimal control problem , defined by (23), there is at least one member problem with that has the same solution, that is, .*

*Proof. *Suppose that the optimal solution for the given problem , defined by (13), is known and set . With Theorem 1, the optimal primer vector function can be written as an explicit function of
Set and plug it into the state equation of (13). We obtain the following relationships for an optimal solution to problem :
Meanwhile, any member problem in can be solved by finding the solution for the following TPBVP:
This solution is
Select a problem with
Put (28) into (27). After simplification, we obtain
Therefore, the optimal control and optimal trajectory of problems and are the same, that is,

If and a transfer problem has the same solution, we call the *equivalent optimal control problem* of the original problem in (13).

Theorem 3. *If is an equivalent optimal control problem of a transfer problem , defined by (13), then the optimal control of another member problem is a feasible control for the transfer problem .*

*Proof. *From the general solution equation (27), we know that no matter what functional form of we pick, the obtained control will always generate a trajectory that satisfies the two end point boundaries. Thus, is a feasible control for the power-limited transfer problem.

Theorem 2 gives us a new way to solve any optimal power-limited transfer problem . That is, we first find its equivalent problem in , and then the optimal solution is given by , , and (27). In fact, from the proof of Theorem 2, we can construct such a problem with ,
Though we do not really know until we solve the original problem , we can always “shape” the matrix function ) based on our assumptions and knowledge of the transfer scenario. Suppose the “shaped” matrix function is and it corresponds to a problem . Theorem 3 indicates that, even though is not exactly , it gives a feasible solution for the original problem as long as . Thus, if we can find a proper satisfying , then the solution for , (27), is a feasible approximate solution for , essentially . In other words, once we find (an approximation of ), the approximate solution to is given by (27) with . Since any corresponds to a unique primer vector function in (27), this process of finding an approximation of is named as *primer vector approximation* (PVA) in this paper. The following section uses this PVA to derive two explicit approximate solutions to the coplanar power-limited orbit transfer problem.

##### 3.2. Primer Vector Approximation under the Low-Thrust Assumption

In this section, we assume the thrusting magnitude is very low (i.e., is very small), such that , becomes a symmetric matrix function and

The thrusting magnitude is usually very low in two kinds of transfers. One is the close range transfer and the other the long duration multirevolution transfer.

###### 3.2.1. Approximate Solution for Close Range Transfer

Because the transfer happens in the vicinity of the initial orbit, a reasonable approximation of is where represents the initial orbit. Put (33) into (27). We obtain the approximate optimal solution for a close range transfer.

Though (27) generally requires numerical integration to get , the approximate optimal trajectory and control are still explicit functions of . Moreover, is independent of the target location . Therefore, for all possible close-range transfers of the initial orbit , we only need to compute one time for one period of the initial orbit. This property is greatly useful when we schedule transfers for a satellite formation around an elliptic reference orbit specified by .

When the initial orbit is circular, the computation of can be done analytically. Using index notation and setting , and , we have With (34), the approximate optimal solution, (27), becomes analytic.

###### 3.2.2. Approximate Solution for Long Duration Multirevolution Transfer

A long duration multirevolution transfer tends to build up its orbital energy and angular momentum monotonically. Thus, it is reasonable to approximate the characteristic matrix function with a linear matrix function

With (35), the matrix function in (27) can be expanded analytically. With the analytic , the approximate optimal solution given by (27) is analytic. (The complete formula for is too lengthy to show here.)

#### 4. Numerical Test

##### 4.1. Cases of Free-Time Rendezvous Problem

Without loss of generality, we can set the gravitation constant and use normalized distance in the test cases.

In Table 1, four coplanar rendezvous cases are chosen to test the accuracy of the two analytic approximations. Cases A and B test the approximate solutions for close range transfers, with case A featuring a circular reference while case B featuring a high elliptic reference. Cases C and D test the approximate solutions for long duration transfers. Case C is a circular-to-elliptic transfer, while case D is a circular-to-circular transfer. We use “Approx. C.R.” to identify the approximate solution for close range transfer, (27) and (23), and use “Approx. L.D.” to identify the approximate solution for long range transfer, (27) and (35).

##### 4.2. Results and Discussions

The results of cases A and B show that “Approx. C.R” captures the close-range-transfer primer vector dynamics. Where the circular reference is used (case A), the percentage error of the cost and the time of flight (TOF) are about 0.5% and 0.1%, respectively. As the eccentricity of the reference orbit increases, the primer vector dynamics becomes more complicated. However, even when the eccentricity of the reference orbit is as high as 0.8 (case B), “Approx. C.R” still captures the primer vector dynamics well. The percentage error of the cost and the TOF are about 1.5% and 0.16%, respectively.

The results of cases C and D show that “Approx. L.D” captures the long-duration-transfer primer vector dynamics. The percentage error of the cost and the time of flight (TOF) are around 1% and 2.3%, respectively, for case C, and around 0.26% and 1.4%, respectively, for case D.

It is interesting to see that “Approx. C.R”, intended for close-range transfers, works very well for long-duration transfers too. It is, however, worth pointing out that, since “Approx. C.R” requires a numerical integration, when the initial orbit is elliptic, it is computationally slower than “Approx. L.D.” But when the initial orbit is circular, the analytic "Approx. C.R" actually becomes computationally faster than “Approx. L.D.”

A significant advantage of the primer vector approximation method is that it precisely follows the dynamic equations of motion and gives a feasible solution. Table 2 shows the missed target errors of case D, which are obtained by simulating Equations (2) with the generated control profiles of the two approximations. The results verify the advantage of the proposed method.

Figures 2 and 3 show the thrust angle and thrust magnitude histories of case B, while Figures 5 and 6 show those of case D. Figures 4 and 7 present the trajectories for case B and case D, respectively. In the figures, the black solid line represents the optimal solution, while the black dashed line is from “Approx. LD” and the gray solid line from “Approx. CR”. Only the “Approx. CR” solution (gray line) and the true optimal solution are shown in Figure 3, because the “Approx. LD” solution is too far away from the optimum to effectively display.

#### 5. Conclusions

This paper uses a new method to approximate solutions for a nonlinear optimal control problem. This method begins with a transformation to push all the nonlinearity to the coefficients of the control terms. Then it analyzes of the adjoint equations and embeds the process of finding solutions into the process of approximating the primer vector curve of the equivalent linear quadratic optimal control problem. This method is powerful and leads to an extremely simple and accurate explicit approximate solution. At the cases of circular close-range transfers and long-duration transfers, analytic approximate solutions exist.

This paper focuses on the display of the primer vector approximation method itself and includes no numerical comparison of this new method with other approximation methods. However, as far as the author knows, the proposed method is quite unique and radically different from other approximation methods. This method has a precise integration of the state dynamics, while other approximation methods more or less approximate the state dynamics. As a result, even though the given transfer scenario largely violates the assumptions of the approximation, the obtained solutions are still feasible. This is a major advantage of using the primer vector approximation method.

#### Nomenclature

: | Standard gravitational parameter of the Earth |

: | Distance from the attraction center to the spacecraft |

: | Polar angle of the spacecraft in the polar coordinate |

: | Radial direction component of the velocity of the spacecraft |

: | Normal direction component of the velocity of the spacecraft |

: | Thrust angle of the spacecraft |

: | Adjoint variable corresponding to |

: | Adjoint variable corresponding to |

: | Adjoint variable corresponding to |

: | Adjoint variable corresponding to |

: | Control acceleration in the radial direction of the spacecraft |

: | Control acceleration in the normal direction of |

: | Augmented Hamiltonian of the polar coordinate formulation of trajectory optimization |

: | Augmented Hamiltonian of the new formulation of trajectory optimization |

: | Objective function of the power-limited transfer problem |

: | Objective function of the equivalent problem |

: | State vector of a space vehicle in the new formulation |

: | Adjoint vector of the the new formulation |

: | Control variable in the new formulation |

: | Primer vector in the new formulation |

: | State vector of the equivalent problem |

: | Adjoint vector of the equivalent problem |

: | Control variable of the equivalent problem |

: | Primer vector of the equivalent problem. |

*Subscripts*

0: | Initial value |

: | Final value. |

*Superscripts*

: | Transpose of a matrix |

*: | Optimal value |

: | First derivative with respect to the polar angle |

: | Second derivative with respect to the polar angle. |

#### References

- T. N. Edelbaum, “Theory of maxima and minima,” in
*Optimization techniques: with applications to aerospace systems*, G. Leitmann, Ed., vol. 5, chapter 1, Academic Press, 1962. View at Google Scholar - T. N. Edelbaum, “Optimum low-thrust rendezvous and station keeping,”
*Journal of Spacecraft and Rockets*, vol. 40, no. 6, pp. 960–965, 2003. View at Publisher · View at Google Scholar - T. N. Edelbaum, “Optimum power-limited orbit transfer in strong gravity fields,”
*AIAA Journal*, vol. 3, pp. 921–925, 1965. View at Publisher · View at Google Scholar - T. N. Edelbaum, “An asymptotic solution for optimum power limited orbit transfer,”
*AIAA Journal*, vol. 4, pp. 1491–1494, 1966. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - C. Zee, “Low constant tangential thrust spiral trajectories,”
*AIAA Journal*, vol. 1, no. 7, pp. 1581–1583, 1963. View at Google Scholar - A. Marinescu, “Optimal low-thrust orbital rendezvous,”
*Journal of Spacecraft and Rockets*, vol. 13, no. 7, pp. 385–392, 1976. View at Publisher · View at Google Scholar · View at Scopus - J. P. Marec and N. X. Vinh, “Optimal low-thrust, limited power transfers between arbitrary elliptical orbits,”
*Acta Astronautica*, vol. 4, no. 5-6, pp. 511–540, 1977. View at Google Scholar · View at Scopus - C. M. Haissig, K. D. Mease, and N. X. Vinh, “Minimum-fuel, power-limited transfers between coplanar elliptical orbits,”
*Acta Astronautica*, vol. 29, no. 1, pp. 1–15, 1993. View at Google Scholar · View at Scopus - J. A. Kechichian, “Reformulation of Edelbaum's low-thrust transfer problem using optimal control theory,”
*Journal of Guidance, Control, and Dynamics*, vol. 20, no. 5, pp. 988–994, 1997. View at Google Scholar · View at Scopus - L. Casalino and G. Colasurdo, “Improved Edelbaum's approach to optimize low earth/geostationary orbits low-thrust transfers,”
*Journal of Guidance, Control, and Dynamics*, vol. 30, no. 5, pp. 1504–1510, 2007. View at Publisher · View at Google Scholar · View at Scopus - Y. Gao and C. A. Kluever, “Analytic orbital averaging technique for computing tangential-thrust trajectories,”
*Journal of Guidance, Control, and Dynamics*, vol. 28, no. 6, pp. 1320–1323, 2005. View at Google Scholar · View at Scopus - N. Markopoulos, “Explicit, near-optimal guidance for power-limited transfers between coplanar circular orbits,”
*Journal of Guidance, Control, and Dynamics*, vol. 19, no. 6, pp. 1317–1325, 1996. View at Google Scholar · View at Scopus - A. E. Petropoulos and J. M. Longuski, “Shape-based algorithm for automated design of low-thrust, gravity-assist trajectories,”
*Journal of Spacecraft and Rockets*, vol. 41, no. 5, pp. 787–796, 2004. View at Google Scholar · View at Scopus - B. J. Wall and B. A. Conway, “Shape-based approach to low-thrust rendezvous trajectory design,”
*Journal of Guidance, Control, and Dynamics*, vol. 32, no. 1, pp. 95–102, 2009. View at Publisher · View at Google Scholar · View at Scopus - D. M. Novak and M. Vasile, “Improved shaping approach to the preliminary design of low-thrust trajectories,”
*Journal of Guidance, Control, and Dynamics*, vol. 34, no. 1, pp. 128–147, 2011. View at Publisher · View at Google Scholar · View at Scopus - T. Carter and M. Humi, “Three-impulse rendezvous near circular orbit determined through a new primer vector analysis,” in
*Proceedings of the AAS/AIAA Spaceflight Mechanics Meeting*, vol. 140 of*Advances in the Astronautical Sciences Series*, 2011.