Solution for Nonlinear Three-Dimensional Intercept Problem with Minimum Energy

Classical orbit intercept applications are commonly formulated and solved as Lambert-type problems, where the time-of-flight (TOF) is prescribed. For general three-dimensional intercept problems, selecting a meaningful TOF is often a difficult and an iterative process. This work overcomes this limitation of classical Lambert’s problem by reformulating the intercept problem in terms of a minimum-energy application, which then generates both the desired initial interceptor velocity and the TOF for the minimum-energy transfer. The optimization problem is formulated by using the classical Lagrangian f and g coefficients, which map initial position and velocity vectors to future times, and a universal time variable x. A Newton-Raphson iteration algorithm is introduced for iteratively solving the problem. A generalized problem formulation is introduced for minimizing the TOF as part of the optimization problem. Several examples are presented, and the results are compared with the Hohmann transfer solution approaches. The resulting minimum-energy intercept solution algorithm is expected to be broadly useful as a starting iterative for applications spanning: targeting, rendezvous, interplanetary trajectory design, and so on.


Introduction
A fundamental problem of astrodynamics is concerned with computing intercept trajectories or interplanetary mission orbit for objects in space [1,2].These calculations are often performed assuming a predetermined time-of-flight (TOF).This is the well-known Lambert's problem [3][4][5][6].Applications of Lambert's problem are common in interplanetary trajectory design, spacecraft intercept, rendezvous, ballistic missile targeting, and so on.These problems are formulated by specifying the initial position vectors of an interceptor and target satellite.When the TOF is specified, the initial velocity vector for the interceptor is an unknown implicit function of the local gravity field and can only be recovered by a successive approximation strategy.Other authors have considered alternative problem formulations for a specified TOF that have included minimum-fuel and multiple-impulse strategies [7,8].This work overcomes the limitations of these approaches by introducing a minimum-energy approach that simultaneously generates both the TOF and the initial velocity vector for the interceptor.
For the special case of coplanar orbits, the Hohmann transfer algorithm generates a two-impulse minimumenergy orbit transfer by using tangential burns [3][4][5].This technique provides a reference orbit transfer for various space applications.For direct applications of the Hohmann transfer to interplanetary orbit transfer, the position vectors of the target planet and initial departure planet are specified assuming a prescribed TOF.When the spacecraft is far from the initial position, one must be alert to the possibility that a multiorbit maneuver may be required.
Clearly, the TOF is a critical parameter for various applications.Once a TOF is determined, the rest of the procedure is solved readily by the solution of Lambert's problem.This work addresses the problem that there are no adequate methods available for determining a TOF, especially, in general three-dimensional (3D) cases.The problem of finding an optimal TOF only becomes well defined when one specifies a minimization criteria.To this end, a minimum-energy version of classical interceptor problem is formulated for recovering the TOF for a 3D orbit transfer.The results of this calculation are useful as a reference value for interplanetary trajectory design, spacecraft intercept, rendezvous, ballistic missile targeting, and so on.Of course, one can also bound the range of achievable transfer trajectory times by solving for the minimum TOF consistent with the maximum energy that can be generated.Yielding a mission design space that spans the range of TOF consists in the range [TOF min , TOF max ].
The design goal for the optimization problem is to simultaneously recover the required initial interceptor velocity and the TOF for the intercept.The mathematical advantage of this approach is that the problem has a unique optimal solution, rather than the family of solutions that characterize the classical Lambert's problem.Mathematically, the problem is defined by a constrained optimization algorithm.Particular care is exercised in formulating the problem for handling the near-parabolic orbits that arise in intercept applications.Analytically, this is handled in a comprehensive way by introducing a universal variable that permits a single TOF equation to be developed that is valid for all conic orbits.This work is organized in three sections.First, Kepler's equation is used to define the TOF equation.This is followed by a description of the universal variable used for the problem formulation.For completeness, Lambert's problem is briefly described.Second, the minimum-energy problem for the intercept problem is introduced and solved.Third, simulation results are presented which compare the TOF obtained for an interplanetary trajectory design with a trajectory developed using the Hohmann transfer methodology and interceptor design solution approaches.

Mathematical Review
A fundamental approach for determining the TOF for spacecraft starts with Kepler's equation that is given by where  is the mean anomaly,  and  denote the eccentric anomaly and the eccentricity, respectively,  is the time of periapsis passage,  is the TOF,  is the mean motion defined as √/ 3 ,  denotes the gravitational constant, and  is the semimajor axis of orbit.
As  ∼ 1, the solution for Kepler's equation becomes more difficult to obtain.This problem is overcome by introducing the universal variable given by [3] where  is the position of spacecraft.As shown in [3], by introducing the universal variable defined by (2), one can express Kepler's equation and the radial spacecraft coordinate in the following form: where  is assumed to be zero without loss of generality and r 0 and v 0 are the initial position and velocity vectors of spacecraft, respectively.These necessary conditions describe the position and velocity of an orbiting object as a function of time.If the value of the universal variable from ( 3) is known, the position of the spacecraft at that time is evaluated.Even though ( 3) is transcendental in , a Newton's iteration technique is used to successfully solve for  when the TOF, , is given.
Assuming that there are no external forces, then the four vectors r 0 , v 0 , r, and v are assumed to be governed by Keplerian motion.To compute v and r in terms of v 0 , r 0 , and , the position and velocity vectors of spacecraft at time  are described as [9] where , , ḟ , and ġ are scalar time-dependent constants, which are subject to the following constraint: where The energy minimum form of Lambert's problem is solved by introducing the classical Lagrangian coefficients and universal variable in the problem necessary conditions.

Time-of-Flight for Minimum-Energy Orbit Transfer
The major objective in this paper is to compute (i) the TOF and (ii) the initial velocity for an interceptor object for two arbitrary given position vectors so that the transfer orbit energy is a minimum.This problem differs from the classical Lambert's problem, which fixes a TOF and only recovers a solution for the initial velocity for the interceptor object, given initial and final position vectors of spacecraft.
The optimization problem is formulated by assuming that a target and an interceptor exist in arbitrary orbits, respectively.The problem geometry is illustrated in Figure 1, where r 0 and v 0 denote the initial position and velocity vectors of the target, respectively, and r0 and v0 represent the initial position and velocity vectors of interceptor, respectively.
The unknowns for the problem are the TOF and initial velocity correction for the interceptor.The goal of the trajectory optimization is to reduce the displacement position vector locating the interceptor relative to the target to zero values at the TOF, while minimizing the orbit energy of the interceptor.The problem is formulated as a nonlinear optimization problem.

Constrained Optimization Problem
For given r 0 , v 0 , and r0 , find  and v0 by minimizing the performance index defined as the interceptor's orbit energy, J, defined as subject to where  and x (9) denote the universal variables for the target and the interceptor, respectively, and  is the TOF to be determined.The displacement vectors for the target and interceptor are expressed using  and  as follows: As a constraint vector,  ∈ R 2 , (3) for  and x is rewritten as Note that the augmented variables to be obtained are v0 , , x, and .

Optimal Necessary Conditions
Since the second term of the energy is constant, it does not affect the performance index so that the index is redefined, without loss of generality, as [10] J The Hamiltonian is formed by appending the constraints of ( 9) with Lagrange multipliers as follows: where To minimize the performance index with respect to the augmented variables, the necessary conditions provide the following [11]: where η/v 0 ∈ R 1×3 represents a row vector, r/v 0 ∈ R 3×3 is a matrix (refer to the Appendix for detail derivation), and The necessary conditions of ( 14)-( 17) are simplified by the following manipulations.First, from (15), the Lagrange multiplier  2 is obtained as Second, substituting (19) into (17) yields which can be solved for , leading to Third, by using ( 14), the Lagrange multiplier  1 follows as Collecting the Lagrange multiplier solutions from ( 19) and ( 22), introducing the results into (16), one obtains This equation is further simplified by recalling the terminal constraint  = r, leading to Substituting (21) into (3) yields the final necessary condition required for finding the TOF for the intercept problem: where the new matrix is defined for simplicity as Satisfaction of the new equation implies that the interceptor can meet the target with minimum energy in a time provided by the computed TOF, not by a predetermined time.

Summary
The approach for obtaining the nonlinear 3D intercept problem by using the classical Lagrangian  and  coefficients is summarized.Computing the TOF and the interceptor's initial velocity is the goal to meet the requirement that the final distance between the two spacecrafts becomes zero.Firstly, one can simply obtain  and  0 using the following: where the orbital energy is given by Then, one can find the universal variables, the initial velocity for the interceptor, and the TOF using the following equations: where the semimajor axis â of the interceptor can be iteratively computed with estimated v0 .The Newton-Raphson iteration algorithm is applied to solve the previous equations.
Next, one can compute all of the  and  expressions using (7).Then, ( 5) is applied to obtain the final position and velocity vectors.
There are many feasible performance indices to specify a TOF.For example, consider the candidate performance index where  is nonnegative weight.By adding the time as one part of the performance index, the TOF is expected to be shortened with respect to the variation of .In a similar manner with the minimum-energy procedure in the previous section, the optimization solution to this problem is readily determined.The partial derivative of the Hamiltonian  with respect to  is given by Finally, a cost-effective equation weighted to the time is obtained as Numerical convergences based on different methods and their overall computational cost depend on the chosen parameterization, the initial guess, and the numerical technique used for solving the resulting equation.It is known that singularities exist when solving Lambert's problem that prevent some algorithms from converging for particular cases or make convergence extremely slow.For example, Lambert's method fails when the transfer angle is 180 degrees [9].Therefore, the features of the suggested method must be analyzed.However, this is out of the scope of this paper, which is focused on approaches to determine the TOF and initial velocity of the interceptor.

Application Examples
The specification of a TOF for an intercept problem is generally not unique, and a family of solutions are possible when the initial trust level is variable.As a result, the process of determining a useful TOF requires experimentation and iteration.The minimum-energy optimization approach of this work finds a unique value for the TOF.The solution for the intercept problem simultaneously determines the initial interceptor velocity vector and TOF.Numerical examples are presented that compare and contrast the classical Hohmann transfer with the proposed method.
Let us briefly review the Hohmann transfer and compare the minimum-energy problem with it.The geometry of the Hohmann transfer is illustrated in Figure 2. The distances of the departure and arrival orbits are denoted as  1 and  2 , respectively.
The semimajor axis for the elliptic orbit and the energy are given by [3] and, then, the departure velocity of the transfer orbit is readily obtained as Since the velocity of the departure orbit is given by the velocity change for the Hohmann transfer is calculated as and the TOF of the Hohmann transfer is written as where  ℎ is obtained from (33).Two circular orbits are assumed with the radii of  1 = 4000 km and  2 = 6000 km, respectively.Then, the velocity of the departure circular orbit is V 1 = 9.9825 km/s, and the remaining parameters for the Hohmann transfer are obtained as  = 1759.3sec, V ℎ0 = 10.9353km/s, and ΔV = 0.9528 km/s.To navigate to the final position of the arrival orbit by the Hohmann transfer, the initial position and velocity vectors are assumed to be given by  Four cases including the Hohmann case are analyzed.Initial velocities, velocity changes, and TOF obtained by the solution of the proposed minimum-energy problem are arranged in Table 1.Also, the initial positions of target spacecraft and their resultant transfer trajectories are displayed in Figure 3.
As shown in Table 1, the result for the Hohmann case is nearly identical to the output from the classical approach in (35)-(38) with a small numerical error.It proves that the proposed approach provides optimal solutions we are looking for.Moreover, it is obvious that it gives the minimum velocity change, which is tangential with the trajectory, compared with the other cases.If the target is positioned at case 1, 2, or 3, relative to the interceptor's initial position and it is required to start the orbit transfer mission at this time, it would be a great advantage to have a reference minimum-energy trajectory to accomplish the mission.Fortunately, the results in Table 1 can be utilized, since they represent the minimum velocity in each case.This means that there are no more efficient trajectories in these cases than the transfer orbit listed in Table 1.When the target is positioned forward compared to the Hohmann transfer, the phase angle, sometimes called flight-path-angle, at departure should be negative to meet the optimal trajectory requirement.When the target, on the contrary, is positioned backward, the flight-path-angle should be positive.
Even if a circular orbit is selected for the comprehensive analysis by comparing with the Hohmann transfer, the application of the proposed approach is not limited.Therefore, an illustrative example in Figure 4 is conducted to demonstrate the performance of space maneuver of the interceptor.There are two arbitrary elliptic orbits, which are not coplanar.The initial positions of the target and interceptor orbit are depicted in Figure 4.By solving the nonlinear 3D  intercept problem through minimizing energy, the TOF and initial velocity is computed.Then, applying minimum-energy velocity obtained for the interceptor confirms that the final distance between the two orbits is zero at the computed TOF.
Note that the problems formulated by the universal variable and  and  expressions in this paper are solved in 3D space for supporting the design of arbitrary intercept problems with minimum energy.The intercept problem is easily generalized by introducing a time weighting factor in the definition of the optimal control problem, which allows a systematic exploration of the intercept design space as a function of the allowed transfer energy.In this example, the outer orbit is selected as a target orbit to be captured as illustrated in Figure 5.The initial position of the interceptor is at the inner orbit.The TOF obtained from the optimization problem is shortened when  increases, and the results are illustrated in Figure 5.As expected, the longest TOF is obtained when  is zero.If  approaches one in this simulation case, the interceptor can hit the target in a very short time.It means that shortening the TOF is the optimal solution to minimize the chosen cost function.

Figure 1 :
Figure 1: Geometry of the minimum-energy problem for a TOF.

Figure 3 :Figure 4 :
Figure 3: Trajectory generation by initial positions of target.

Figure 5 :
Figure 5: TOF for intercept due to variation of .