Next Article in Journal
Factors Impacting Chinese and European Vertical Fight Efficiency
Previous Article in Journal
Assessment of the Potential of Electric Propulsion for General Aviation Using Model-Based System Engineering (MBSE) Methodology
Previous Article in Special Issue
Design of Orbit Controls for a Multiple CubeSat Mission Using Drift Rate Modulation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Variation of Osculating Orbit Elements Using Low-Thrust Photonic Laser Propulsion in the Two-Body Problem

Department of Aerospace Engineering, Tamkang University, Tamsui, New Taipei City 25137, Taiwan
Submission received: 1 January 2022 / Revised: 20 January 2022 / Accepted: 27 January 2022 / Published: 30 January 2022
(This article belongs to the Special Issue Spacecraft Trajectory Design and Optimization)

Abstract

:
This study investigated the variation of the osculating orbit elements of a spacecraft propelled by photonic laser propulsion (PLP) under the two-body problem assumption. The PLP thrusting system can produce continuous and constant thrust. This paper first reviewed its basics and then studied its influences on the variation of osculating orbit elements given a small PLP thrust. Gauss’s equations, perturbation theory, and normalization were introduced to investigate this problem. Our work approached the problem by studying the influences of small planar and out-of-plane PLP thrusts, respectively. Bounds on the variation of orbit elements were derived, and a sufficient condition that traps the mission spacecraft in the vicinity of the mother ship was also found. Numerical simulations are also presented to verify our results, including the bounds and the sufficient conditions. The results obtained in this paper are directly applicable to the usage of PLP thrust, a new type of thrusting system, in the future, and are potentially helpful to various space missions, especially interplanetary travel.

1. Introduction

This paper studied the trajectory variation of spacecraft propelled by a low-thrust photonic laser propulsion (PLP) system under the two-body problem environment. In the past few decades, many researchers have shown great interest in the photon thruster and its applications to interplanetary travel, because of its highly efficient and continuous thrusting properties [1,2,3,4,5]. In 2002, Thomas R. Meyer et al. [6] studied the idea of the laser elevator by momentum transfer using an optical resonator. Young K. Bae [7] proposed the concept of photonic laser propulsion in 2008. In that study, the PLP thruster was modeled as a continuous, low-thrust, and high-specific-impulse ( I s p ) engine. A potential application of the PLP-driven missions proposed by these studies is interplanetary travel. A journey to the Alpha Centauri star system was even proposed in 2016 by the Breakthrough Starshot project by Yuri Milner, Stephen Hawking, and Mark Zuckerberg [8].
However, these research works focused more on the PLP thruster itself, and less on the influence on the trajectory. A common assumption is seen in these research works: since the laser continuously pushes and accelerates the mission spacecraft, the mission spacecraft will eventually reach a very high speed, and interplanetary or interstellar travel can be accomplished in a very short traveling time, regardless of the practical problems of the aiming and the relativity effects. It was not until 2011 that Wang and Hsiao [9] investigated the trajectory of the PLP spacecraft under the two-body-problem assumption. Wang and Hsiao realized that some constraints still exist in the applications of PLP thrust [9]. Equilibrium points exist in the planet–mother ship–mission spacecraft system. The mission spacecraft will be trapped in the vicinity of the mother ship if the PLP force is smaller than a threshold, as proposed in [9].
On the other hand, much attention has been focused on a continuous, low-thrust engine, the efficiency of which was proven in the Deep Space 1 mission by NASA [10] and in a lunar mission by the European Space Agency (ESA) [11]. Several interplanetary missions demonstrated the feasibility of using a low-thrust engine, such as by electric propulsion, as the main propulsion system of the spacecraft. The idea of using the power of light has been studied for decades. Many space missions using light propulsion have been implemented using solar sails. The drawbacks of solar light are that it is passive and its magnitude is inversely proportional to the square of the distance. As a result, the propulsive force decays sharply as long as the spacecraft is far away from the Sun. A continuous and constant PLP force will suppress these drawbacks.
This article is a continuation of the work of [9]. As mentioned earlier, a criterion is imposed on the magnitude of the PLP thrust for practical applications. We were interested more in the influences of a small PLP thrust on the trajectory, because in the current technology, it is infeasible to generate a very large laser power from a space-borne base. By studying the spacecraft trajectory driven by a small PLP thrust, we are able to propose a reasonable laser power criterion for the application of PLP thrust for interplanetary or interstellar travel.
A preliminary work was presented in [12]. The results described here, as an extension of the preliminary work, derived more bounds on the variation of the planar orbit elements, adding a full section discussing the out-of-plane cases and providing many more numerical simulations for its verification. First, this paper briefly reviewed some facts and theories regarding the PLP system. The Gauss equations were employed to understand the variation of the osculating orbit elements. Normalization was also employed to generalize the results. Perturbation theory was applied to study the zeroth- and first-order approximations provided that the mission spacecraft is driven by a small PLP thrust. With certain assumptions, we were able to obtain bounds on the variation of the osculating orbit elements. Numerical simulations are provided to verify the proposed algorithms.
The paper is arranged as follows: Section 2 reviews some facts about PLP thrust; Section 3 reviews the assumptions of the investigated system and the normalization process; Section 4 briefly reviews the Gauss equations and perturbation theory; Section 5 and Section 6 study the variation of the osculating orbit elements by applying planar and out-of-plane thrusting forces, respectively; Section 7 provides some numerical simulations to verity the proposed algorithms; the last section of the paper concludes our work.

2. Photonic Laser Propulsion

Many scientists have discussed the concept of using a laser to provide thrust. However, lasers are very inefficient at generating thrust. Although photonic engines have a larger specific impulse compared to conventional ones, they have a smaller thrust-to-power ratio [7]. The specific impulse is approximately I s p = 3.06 × 10 7 s, whereas the thrust-to-power ratio is approximately T / P = 3.34 × 10 9 N/W.
Bae proposed an active resonant optical cavity between two space platforms. In this design, the photon thrust F produced on each mirror is given by [7]:
F = E c t ,
where E is the energy of each photon, c is the speed of light, and t is the interaction time. By some manipulations, the thrust can be expressed as [7]:
F = 2 P R m S c ,
where R m is the mirror reflectance and S is the apparent photon-thrust-amplification factor, defined as the ratio of the intracavity laser power to the extracavity laser power P. The term S is approximated by:
S = 1 1 R m .
The roof-of-concept demonstration can be found in [7]. Figure 1 presents some experimental data. According to Equations (2) and (3), the thrust is a continuous and stable force. Even though the force is very small, the continuous force keeps driving the spacecraft until it reaches the desired velocity. Figure 2 illustrates the application of PLP to a spacecraft. The launching process starts with a mother ship, which emits a laser beam at the mission ship to generate thrust. Because of the conservation of momentum, however, the mother ship moves in the opposite direction of the mission ship. Thus, a conventional thruster installed on the mission ship must act against the momentum caused by PLP to prevent the mother ship from falling out of orbit.

3. Equations of Motion

Figure 3 shows the relative locations of the central body, mother ship, and mission spacecraft. Let r be the position vector of the mission spacecraft, R the position vector of the mother ship, and r R be the relative position of the spacecraft with respect to the mother ship. Assume the masses of the mother ship and the spacecraft are negligibly small so that they do not produce any gravitational force. According to Newton’s gravitation law and Newton’s second law of motion, the equations of motion (EOMs) of the mother ship and spacecraft are given by:
R ¨ = μ R 3 R ,
r ¨ = μ r 3 r + F L ^ ,
where μ is the gravitational parameter of the central body, L = r R , L ^ = L / | L | , and F is the PLP force given by Equation (2). Note that Equations (4) and (5) are described in the inertial frame (IF).
Two coordinate systems were introduced in the inertial frame in the study. First of all, the J2000 Coordinate System, represented by ( I , J , K ) , was introduced to specify the orbit elements. In this system, I points toward the vernal equinox, K points toward the North Pole of the Earth, and J completes the triad by the right-hand rule. Although J2000 is defined by the Earth, it is also widely used in the exploration of other celestial bodies. The other coordinate system is represented by ( i , j , k ) . In this system, i points toward the initial position of the mother ship, k points toward the angular momentum of the mother ship, and J = k × i .
This study did not consider the reaction force by the PLP on the mother ship because in practical applications, we assumed that this reaction is counteracted by a traditional propulsion system. Moreover, this study also assumed that F is constant and acts along the relative position of the mission spacecraft and the mother ship.

4. Variation of Osculating Orbit Elements and Perturbation Theory

4.1. Osculating Orbit Elements and Normalized Gauss Equations

The motion described by Equation (5) is Keplerian if F 0 . A Keplerian trajectory can be described by six canonical orbit elements: a the semi-major axis, e the eccentricity, i the inclination, Ω the longitude of ascending node, ω the argument of periapsis, and M the mean anomaly. These elements are functions of the spacecraft position r ( t ) and velocity v ( t ) in the IF. Without external perturbation forces, the orbit elements remain constant over time.
When an external force is exerted on a spacecraft, the variation of its osculating orbit elements is described by the Gauss-type planetary equations [13]. Moreover, normalization is helpful to generalize applications. In this paper, the only perturbation results from the PLP thrust between the mother ship and the mission ship. Hence, it is natural to normalize this system with the mother ship’s parameters.
Assume a spacecraft is suffering from a continuous force F = F r e ^ r + F f e ^ f + F h e ^ h , where e ^ r and e ^ h are the unit vectors along the position and angular momentum, respectively, and e ^ f = e ^ h × e ^ r . Define n R = μ / R 3 to be the mean motion of the mother ship, and some fundamental parameters are normalized as follows:
r ˜ = r R ,
τ = n R t ,
F ˜ = F R 2 μ = F r R 2 μ e ^ r + F f R 2 μ e ^ f + F h R 2 μ e ^ h = F ˜ r e ^ r + F ˜ f e ^ f + F ˜ h e ^ h ,
a ˜ = a R ,
n ˜ = n n R = 1 a ˜ 3 .
Since n ˜ = μ ˜ / a ˜ 3 , we conclude that μ ˜ = 1 . The elements e, i, ω , Ω , and f are already dimensionless. Hence, e ˜ = e , i ˜ = i , and so on. Moreover, the derivative of a differential function X ( t ) with respect to dimensionless time τ can by found by the chain rule, given by,
d X d τ = d X d t d t d τ = 1 n R d X d t .
According to these definitions, the normalized parameters and elements of the mother ship are determined as follows:
R ˜ ( t ) = R ( t ) R = R ^ ( t ) ,
V ˜ ( t ) = d R ˜ d τ = 1 n R R d R d t = V ( t ) V .
Consequently, R ˜ ( t ) = e ^ r , V ˜ ( t ) = e ^ f , and μ ˜ = 1 , and the values of the following parameters can be found by:
E ˜ m = 1 2 V ˜ 2 μ ˜ R ˜ = 1 2 ,
a ˜ m = a m R = 1 ,
n ˜ m = n R n R = 1 ,
where the subscript m denotes the mother ship. Moreover, e m = 0 , i m = i 0 , Ω m = Ω 0 , and ω m = ω 0 . Without lost of generality, we set f m = 0 for the mother ship at t = 0 .
Accordingly, the normalized Gauss equations are given as follows:
d Ω d τ = n ˜ a ˜ r ˜ 1 e 2 F ˜ h sin u csc i ,
d i d τ = n ˜ a ˜ r ˜ 1 e 2 F ˜ h cos u ,
d e d τ = n ˜ a ˜ 2 1 e 2 F ˜ r sin f + F ˜ f cos f + cos E ,
d ω d τ = n ˜ a ˜ 2 e 1 e 2 F ˜ r cos f + F ˜ f 1 + r ˜ p ˜ sin f cos i d Ω d τ ,
d a ˜ d τ = 2 n ˜ a ˜ 2 F ˜ r a ˜ e 1 e 2 sin f + F ˜ f a ˜ 2 1 e 2 r ˜ ,
d n ˜ d τ = 3 2 n ˜ a ˜ d a ˜ d τ .
where f is the true anomaly and u = f + ω .
The remaining part of the paper focuses on the normalized Gauss equations. To simplify future derivations, the tilde signs in Equations (17) to (22) are dropped in the following part of the paper.
Notable, perturbations may cause osculating or non-osculating orbits, depending on the types of perturbations [14,15]. In the problem discussed in this article, orbits are osculating because the PLP thrust is not velocity dependent. Consequently, the term “orbit elements” refers to “osculating orbit elements” in the remainder of this paper.

4.2. Review of Perturbation Theory

This paper focused on the influence of small PLP thrust. Hence, perturbation theory can be employed to understand and approximate the variation of the orbit elements. The following is a brief review of perturbation theory.
Consider a differentiable function x ( t ) that satisfies the following differential equation:
d x d t = ϵ f ( x ) ,
where ϵ 1 . Suppose the solution to the differential equation is:
x ( t ) = x ( 0 ) ( t ) + ϵ x ( 1 ) ( t ) + ϵ 2 x ( 2 ) ( t ) + O ( ϵ 3 ) ,
where x ( 0 ) ( t ) , x ( 1 ) ( t ) , x ( 2 ) ( t ) , etc., are assumed differential.
Substituting Equation (24) into Equation (23) yields:
x ˙ ( 0 ) + ϵ x ˙ ( 1 ) + ϵ 2 x ˙ ( 2 ) + = ϵ f ( x ( 0 ) + ϵ x ( 1 ) + ϵ 2 x ( 2 ) + ) = ϵ f ( x ( 0 ) ) + f ϵ ϵ = 0 ϵ + .
Comparing the coefficients between two sides leads to:
d x ( 0 ) d t = 0 ,
d ϵ x ( 1 ) d t = ϵ f ( x ( 0 ) ) ,
Equation (26) determines a constant x ( 0 ) , and ϵ x ( 1 ) can be found by:
ϵ x ( 1 ) = ϵ f ( x ( 0 ) ) d t .
If ϵ is small enough, the solution x ( t ) can be approached by the first-order approximation as follows:
x ( t ) x ( 0 ) + ϵ x ( 1 ) ( t ) .
Note that the notations for the zeroth-order approximation are changed to x ¯ in the following derivations to make the equations neat, that is,
x ( t ) x ¯ ( t ) + ϵ x ( 1 ) ( t ) .

5. Influence of Planar PLP Thrust

If the PLP force is planar, i.e., F = F r e ^ r + F f e ^ f , it is apparent that d Ω / d τ = 0 and d i / d τ = 0 . The remaining equations to solve are Equations (19)–(22). Notably, d Ω / d τ in Equation (20) equals zero.

5.1. Zeroth-Order Approximation

The geometric parameters of the mother ship and mission spacecraft are depicted in Figure 4. As determined in Equation (26), the zeroth-order approximation is the unperturbed orbital elements of the mission spacecraft. Here, we assumed that the mission spacecraft deviates from the mother ship by a normalized offset | | L ( 0 ) | | = l 0 1 , but keeps a velocity identical to the mother ship, i.e., v = V = j at t = 0 . A potential scenario to implement this assumption is that the mother ship holds the mission spacecraft with robotic arms to place it at the launch site in the vicinity of the mother ship. When the laser beam is initiated, the robotic arms release and the spacecraft is driven away.

5.1.1. Semi-Major Axis

The zeroth-order approximation of semi-major axis a ¯ can be obtained from the initial energy. According to the assumption, the initial energy of the mission spacecraft is obtained by:
E = 1 2 × 1 2 1 r .
Since 1 l 0 r 1 + l 0 , it is conclusive that:
E min = 1 2 1 1 l 0 = 1 2 ( 1 l 0 ) 1 1 2 ( 1 + l 0 ) = 1 2 l 0 .
Similarly, E max 1 / 2 + l 0 . As a result,
a ¯ min = 1 2 E min = ( 1 2 l 0 ) 1 1 2 l 0 .
Similarly, a ¯ max = 1 + 2 l 0 . Consequently, we conclude that a ¯ 1 + 2 l 0 .

5.1.2. Eccentricity

To find the initial eccentricity of the mission spacecraft, we had to compute the initial angular momentum. Let ψ be the angle between L and R . Figure 4 indicates that r = ( 1 + l 0 cos ψ ) i + l 0 sin ψ j and v = j and:
h = | | r × v | | = 1 + l 0 cos ψ ,
and r = 1 + 2 l 0 cos ψ + l 0 2 . The eccentricity is computed by:
e = 1 + 2 E h 2 μ 2 = 1 + 2 1 2 1 r ( 1 + l 0 cos ψ ) 2 = 1 + 1 2 ( 1 + 2 l 0 cos ψ + l 0 2 ) ( 1 / 2 ) ( 1 + l 0 cos ψ ) 2 l 0 2 + k 1 l 0 3 + k 2 l 0 4 l 0 ,
where k 1 = 2 cos ψ 4 cos 3 ψ , k 2 = cos 2 ψ 3 cos 4 ψ , and the approximation of 1 / r is detailed in Appendix A. As a result, e ¯ = l 0 .

5.1.3. Longitude of the Ascending Node and Argument of Periapsis

Subject to an initial planar offset, only the magnitude of the angular momentum of the mission spacecraft may change, but the axis of the angular momentum will not. This indicates that the longitude of ascension remains identical to the mother ship, i.e., Ω ¯ = Ω m .
However, the initial offset of the mission spacecraft may cause the shift of the eccentricity axis. Since the argument of periapsis is defined as the angle between the eccentricity axis and the axis of the ascending node, the zeroth-order of the argument of periapsis may change. An easier approach is to find the angle difference from the mother ship, i.e., ω ¯ = ω m + Δ ω , where Δ ω satisfies:
e · R = e cos Δ ω .
Recall that r = ( 1 + l 0 cos ψ ) ) i + l 0 sin ψ j and h = r × v = ( 1 + l 0 cos ψ ) ) k . Then,
e = ( v × h ) / μ r ^ = ( 1 + l 0 cos ψ ) i ( 1 + l 0 cos ψ ) i + l 0 sin ψ j 1 + l 0 2 + 2 l 0 cos ψ .
Consequently,
e · R = ( 1 + l 0 cos ψ ) 1 + l 0 cos ψ 1 + l 0 2 + 2 l 0 cos ψ = ( 1 + l 0 cos ψ ) 1 1 1 + l 0 2 + 2 l 0 cos ψ 2 l 0 cos ψ ( 1 + l 0 cos ψ ) 2 l 0 cos ψ .
Since e ¯ l 0 , we conclude that:
Δ ω = cos 1 2 cos ψ .

5.1.4. Mean Motion

The mean motion is a function of the semi-major axis, given by n = μ / a 3 . In the normalized system,
n = 1 a 3 = a 3 2
Then, the zeroth-order approximation can be found by:
n ¯ = a ¯ 3 2 = ( 1 + Δ a ) 3 2 1 3 2 Δ a
Since a ¯ [ 1 2 l 0 , 1 + 2 l 0 ] , we conclude that:
1 3 l 0 n ¯ 1 + 3 l 0

5.2. First-Order Approximation and Bounds of Elements under Small Thrust

Suppose the PLP force is small enough so that the right-hand-side functions in Equations (17) to (22) are very small. The perturbation theory can be applied to approximate the solutions. We only consider the influence on planar motion in this section.

5.2.1. Semi-Major Axis

According to the conclusion in Section 4.2, the evolution of the semi-major axis can be approximated as:
a a ¯ + ϵ a ( 1 ) + ,
where a ¯ is the solution of the unperturbed system. Section 5.1.1 concluded that a ¯ [ 1 2 l 0 , 1 + 2 l 0 ] , depending on the initial offset of the mission spacecraft.
The first-order approximation ϵ a ( 1 ) can be found by:
ϵ a ( 1 ) = ϵ f ( a ¯ , e ¯ , n ¯ ) d τ = 2 n ¯ a ¯ 2 F r a ¯ e ¯ 1 e ¯ 2 sin f + F f a ¯ 2 1 e ¯ 2 r d τ
The integration in the time domain is very complicated. Instead, we integrated this system in the eccentric-anomaly domain. According to Kepler’s law, n τ = E e sin E . We can change the variables by:
n d τ = ( 1 e cos E ) d E
Moreover, the following identities were employed:
r = a ( 1 e 2 ) 1 + e cos f = a ( 1 e cos E )
cos f = cos E e 1 e cos E
sin f = 1 e 2 sin E 1 e cos E
Then, the integration can be rewritten as:
ϵ a ( 1 ) = 0 E 2 a ¯ 3 F r e ¯ sin E + F f 1 e ¯ 2 d E ,
where E 0 . The assumption of E 0 holds for the whole article.
As shown in Figure 4, F r = F cos ϕ and F f = F sin ϕ , where ϕ satisfies:
l 2 + r 2 2 l r cos ϕ = 1
Moreover,
1 sin ϕ = l sin ( f τ )
r = a ( 1 e 2 ) 1 + e cos f
Solving Equations (50) to (52) yields ϕ as a function of f, i.e., ϕ = ϕ ( f ) . Since the true anomaly f is a function of E, described in Equations (47) and (48), by the change of variables, we eventually obtained ϕ = ϕ ( E ) .
The integration can be written as:
ϵ a ( 1 ) = 0 E 2 a ¯ 3 F e ¯ cos ϕ sin E + sin ϕ 1 e ¯ 2 d E .
Define:
cos α a = e ¯ sin E e ¯ 2 sin 2 E + 1 e ¯ 2 = e ¯ sin E 1 e ¯ 2 cos 2 E
sin α a = 1 e ¯ 2 1 e ¯ 2 cos 2 E
Then, Equation (53) can be rewritten as:
ϵ a ( 1 ) = 0 E 2 a ¯ 3 F 1 e ¯ 2 cos 2 E cos ( ϕ α a ) d E
Since 1 e ¯ 2 cos 2 E 1 and 1 cos ( ϕ α a ) 1 , the integration satisfies,
ϵ a ( 1 ) = 0 E 2 a ¯ 3 F 1 e ¯ 2 cos 2 E cos ( ϕ α a ) d E 0 E 2 a ¯ 3 F d E = 2 E F a ¯ 3
If we integrate over one period, i.e., E = 2 π , then ϵ a ( 1 ) 4 π F a ¯ 3 . Accordingly,
a 1 ± 2 l 0 + 4 π F ( 1 ± 2 l 0 ) 3 1 ± 2 l 0 + 4 π F ( 1 ± 6 l 0 ) = ( 1 + 4 π F ) ± ( 2 ± 24 π F ) l 0
However, Equation (58) holds only if the PLP force is small enough. In other words, ϵ a ( 1 ) a ¯ . A reasonable assumption is that:
ϵ a ( 1 ) 4 π F a ¯ 3 0.1 a ¯
Then, we conclude:
F 1 40 π a ¯ 2 1 40 π ( 1 4 l 0 ) 0.008
Equation (60) provides a sufficient condition. If F is smaller than 0.008, the trajectory will not leave the vicinity of the mother ship for a while. Actually, approximations and assumptions were made during the derivation. The upper bound of F for a small perturbation can be larger. The result also agrees with the upper bound found in Hsiao’s paper [9], where a more rigorous bound was applied by F l 0 / 16 .

5.2.2. Eccentricity

A similar procedure can be applied to analyze the evolution of eccentricity. With the change of variables, Equation (19) can be rewritten as:
ϵ e ( 1 ) = 0 E a ¯ 2 1 e ¯ 2 F r 1 e ¯ 2 sin E + F f ( 2 cos E e ¯ e ¯ cos 2 E ) d E
Defining K e = 1 e ¯ 2 and applying the identity cos 2 E = ( 1 + cos 2 E ) / 2 yield:
ϵ e ( 1 ) = 0 E a ¯ 2 K e F r K e sin E + F f 2 cos E 3 e ¯ 2 e ¯ cos 2 E 2 d E = 0 E a ¯ 2 F K e cos ϕ K e sin E + sin ϕ 2 cos E 3 e ¯ 2 e ¯ cos 2 E 2 d E
Define
K = ( K e sin E ) 2 + 2 cos E 3 e ¯ 2 e ¯ cos 2 E 2 2 = e ¯ 2 cos 4 E 4 e ¯ cos 3 E + ( e ¯ 4 + 4 e ¯ 2 + 3 ) cos 2 E 4 e ¯ cos E + ( e ¯ 2 + 1 + e ¯ 4 ) 0
cos α e = K e sin E K
sin α e = 2 cos E 3 e ¯ / 2 e ¯ cos 2 E / 2 K
Equation (62) becomes:
ϵ e ( 1 ) = 0 E a ¯ 2 F K e K cos ( ϕ α e ) d E
The parameter K in Equation (63) is a function of e ¯ and cos E . The mission spacecraft is assumed to launch very close to the mother ship. It is reasonable to assume 0 e ¯ 1 . Moreover, 1 cos E 1 .
Figure 5 gives the simulation of the maximum value versus different e ¯ ’s. We can see that, over all possible e ¯ ’s, K 2 . As a result, we conclude that Equation (62) leads to:
ϵ e ( 1 ) = 0 E a ¯ 2 F K e K cos ( ϕ α e ) d E 0 E 2 a ¯ 2 F K e cos ( ϕ α e ) d E 2 E a ¯ 2 F K e
If the system is integrated over one period, i.e., E = 2 π , we obtain:
ϵ e ( 1 ) 4 π a ¯ 2 F 1 e ¯ 2
To satisfy the assumption of small perturbation, we apply F ( 1 / 40 π ) and obtain:
ϵ e ( 1 ) a ¯ 2 1 e ¯ 2 10
Thus,
e e ¯ + a ¯ 2 1 e ¯ 2 10 l 0 + a ¯ 2 10
Since 1 2 l 0 a ¯ 1 + 2 l 0 , the upper bound of e varies between:
l 0 + 1 10 ( 1 2 l 0 ) 2 l 0 + 1 10 ( 1 4 l 0 ) = 0.1 + 0.6 l 0   or   0.1 + 1.4 l 0
As a result, we conclude that e 0.1 + 1.4 l 0 .

5.2.3. Argument of Periapsis

The first-order approximation of the argument of periapsis is governed by Equation (20). In this case, only a planar PLP force is applied. It is conclusive that d Ω / d τ = 0 , implying that F r and F f are the only two terms that influence the variation of periapsis. The first-order approximation is attainable by integrating Equation (20), given by:
ϵ ω ( 1 ) = 0 2 π n ¯ a ¯ 2 e ¯ 1 e ¯ 2 F ¯ r cos f + F ¯ f 1 + r p sin f d τ = 0 2 π a ¯ 3 K e e ¯ F ¯ r ( cos E e ) + F ¯ f 1 + 1 e cos E K e sin E d E
Recall that F r = F cos ϕ and F f = F sin ϕ . The integral can be further written as:
ϵ ω ( 1 ) = a ¯ 3 F e ¯ 0 2 π K 1 ( e ¯ , E ) cos ϕ + K 2 ( e ¯ , E ) sin ϕ d E
where K 1 ( e ¯ , E ) = K e ( cos E e ¯ ) and K 2 ( e ¯ , E ) = K e + sin E e ¯ sin E cos E . Define K ω = K 1 2 + K 2 2 , and:
sin α ω = K 1 K ω
cos α ω = K 2 K ω
The integration can be shortened as:
ϵ ω ( 1 ) = a ¯ 3 F e ¯ 0 2 π K ω sin ( ϕ + α ω ) d E a ¯ 3 F e ¯ 0 2 π K ω d E ,
where K ω is a function of e ¯ and E. Because K ω is highly nonlinear, it is difficult to integrate Equation (76) analytically. Let K ω max ( e ¯ ) denote the maximum value of K ω ( e ¯ , E ) over E = [ 0 , 2 π ] given e ¯ . Then, the upper bound of Equation (76) can be found by:
ϵ ω ( 1 ) 2 π K ω max ( e ¯ ) a ¯ 3 F e ¯ .
The value of K ω max ( e ¯ ) for a given e ¯ is presented in Figure 6 numerically. One can see that K ω ( e ¯ , E ) 2.1343 for all potential eccentricities. Since e ¯ = l 0 , it is reasonable to assume that e ¯ 0.1 and K ω 2.02 according to Figure 6. Moreover, to satisfy the assumption of small perturbation, we let F 0.008 , or equivalently, 4 π a ¯ 3 F 0.1 a ¯ , as in our previous discussion. This assumption leads to 2 π a ¯ 3 F 0.05 a ¯ , and the first-order approximation can be concluded as:
ϵ ω ( 1 ) 2.02 · 2 π a ¯ 3 F e ¯ = 0.1 a ¯ l 0 = 0.1 1 ± 2 l 0 l 0 0.1 l 0 ± 0.2
with l 0 0.1 . Note that this result may have a large error for l 0 0 .

5.2.4. Mean Motion

The evolution of mean motion is described in Equation (22). Similarly, under small perturbations, the mean motion can be approximated by n n ¯ + ϵ n ( 1 ) with ( a , e ) substituted by their zeroth-order approximations.
By handingthe integration variable to E, Equation (22) is rewritten as:
d ϵ n ( 1 ) d E = 3 2 n ¯ a ¯ d a d E = 3 2 n ¯ a ¯ d ( a ¯ + ϵ a ( 1 ) ) d E = 3 2 n ¯ a ¯ d ϵ a ( 1 ) d E
Integrating the two sides over 2 π , with the result of Equation (57), yields:
ϵ n ( 1 ) = 3 2 n ¯ a ¯ d ϵ a ( 1 ) d E d E 3 2 n ¯ a ¯ ( 4 π F a ¯ 3 )
The lower bound of ϵ n ( 1 ) depends on the value of 4 π F a ¯ 3 . According to Equation (59), we may let 4 π F a ¯ 3 = K n a ¯ with 0 K n 0.1 . Plugging the result into Equation (80) yields:
ϵ n ( 1 ) 3 K n 2 n ¯
Accordingly, the approximation of n is obtained by:
n = n ¯ + ϵ n ( 1 ) 1 3 K n 2 n ¯
Moreover, 1 3 l 0 n ¯ 1 + 3 l 0 , and the lower bound of the mean motion varies between:
1 3 K n 2 ( 1 3 l 0 ) 1 3 K n 2 3 l 0
As a result, we conclude that n 1 3 K n 2 3 l 0 with 0 K n 0.1 .

6. Influence of Out-of-Plane PLP Thrust

The out-of-plane PLP thrust applies along the direction perpendicular to the orbital plane of the the mother ship initially. All parameters are shown in Figure 7. The initial position offset is assumed L ( 0 ) = l 0 k . Hence, the position and velocity of the mission spacecraft are, respectively, given by r = i + l 0 k and v = j .
Although the out-of-plane cases were studied in this paper, in practice, a pure out-of-plane PLP thrust for the whole mission is infeasible. When the mission spacecraft is initially driven by an out-of-plane PLP thrust, in order to aim at the mission spacecraft, the motion of the spacecraft will naturally cause planar components of the PLP force. As a result, a practical motion subject to initial out-of-plane PLP thrust will eventually be described by the combination of planar and out-of-plane effects.

6.1. Zeroth-Order Approximation

6.1.1. Semi-Major Axis, Eccentricity, and Mean Motion

According to Figure 7, it can be found that:
r = 1 + l 0 2 1
h = r × v = l 0 i + k .
Consequently, h = | | h | | = 1 + l 0 2 1 . Moreover, the value of the mission ship energy, with μ = 1 , is computed as:
E = 1 2 × v 2 μ r = 1 2
The approximated semi-major axis is:
a ¯ = μ 2 E = 1
The approximated eccentricity is:
e ¯ = 1 + 2 E h 2 μ 2 = 0
The mean motion n ¯ is a function of semi-major axis a ¯ . Since the semi-major axis does not change in the zeroth-order approximation, the mean motion remains unchanged as well, i.e., n ¯ = 1 .

6.1.2. Inclination

The inclination of the mission spacecraft is influenced by the orbital elements of the mother ship. Let i m , ω m , and Ω m denote the inclination, argument of periapsis, and longitude of the ascending node of the mother ship, respectively. The transformation of the orbit reference to the J2000 coordinate system is given by [13]:
p J = R x ( f m ) R z ( Ω m ) R x ( i m ) R z ( ω m ) p O R
where p O R = x i + y j + z k is a vector with respect to the orbital reference and p J = X I + Y J + Z K is with respect to the J2000 system. Moreover, R x ( θ ) and R z ( θ ) are the rotation matrices with respect to the i- and k-axis, respectively.
Recall that h = l 0 i + k . Therefore,
cos i ¯ = h · K | | h | | = 1 1 + l 0 2 ( cos i m l 0 sin i m sin ω m )
with f m assumed to be zero. The inclination of mission spacecraft can be found by:
i ¯ = cos 1 cos i m l 0 sin i m sin ω m 1 + l 0 2 cos 1 cos i m + d cos 1 ( χ ) d χ d χ d l 0 l 0 = 0 l 0 = i m 1 1 χ 2 ( l 0 sin i m sin ω m ) i m + sin i m sin ω m | sin i m | l 0
where:
χ = cos i m l 0 sin i m sin ω m 1 + l 0 2
Note that Equation (91) is applicable only if i m 0 . Otherwise, this function is singular. In the case that i m = 0 , we have to approach the problem in a different way. When i m = 0 , Equation (90) degeneratesto:
cos i ¯ = 1 1 + l 0 2
The Taylor expansion of cos i ¯ is:
cos i ¯ = 1 i ¯ 2 2 ! +
The second-order approximation of the right-hand side gives:
1 1 + l 0 2 = ( 1 + l 0 2 ) 1 / 2 = 1 l 0 2 2 ! +
Equating the two sides yields i ¯ l 0 .

6.1.3. Longitude of Ascending Node and Argument of Periapsis

Recall that the angular momentum of the mission spacecraft is h = l 0 i + k . As a result, the ascending node n can be found by:
n = K × h | | K × h | | = n x I + n y J n x 2 + n y 2 ,
where:
n x = l 0 ( sin Ω m cos ω m + cos Ω m cos i m sin ω m ) + cos Ω m sin i m
n y = l 0 ( cos Ω m cos ω m sin Ω m cos i m sin ω m ) + sin Ω m sin i m
As a result, the zeroth-order approximation of the longitude of ascending node is obtained by:
Ω ¯ = tan 1 n y n x
Given l 0 = 0 , Ω = Ω m .
The transverse of the node, denoted as n , can be found by:
n = h × n | | h × n | | = n x I + n y J + n z K n x 2 + n y 2 + n z 2
where:
n x = 1 + l 0 2 cos i ¯ sin Ω m sin i m l 0 cos Ω m cos ω m + l 0 sin Ω m cos i m sin ω m
n y = 1 + l 0 2 cos i ¯ cos Ω m sin i m + l 0 sin Ω m cos ω m + l 0 cos Ω m cos i m sin ω m
n z = l 0 2 ( cos 2 i m cos 2 i m cos 2 ω m + cos 2 ω m ) + 2 l 0 cos i m sin i m sin ω m + sin 2 i m
and the quantity 1 + l 0 2 cos i ¯ is presented in Equation (90). Consequently,
tan ω ¯ = r · n r · n = sin i m sin ω m + l 0 cos i m cos ω m sin i m
Given l 0 = 0 , ω = ω m .

6.2. First-Order Approximation

6.2.1. Inclination

The variation of inclination subject to a small continuous thrust is described in Equation (18). Consider an out-of-plane, small PLP thrust F acting on the mission ship, i.e., F h = F . Integrating Equation (18) over 2 π :
ϵ i ( 1 ) = 0 2 π n ¯ a ¯ r ¯ 1 e ¯ 2 F cos ( f + ω ¯ ) d τ
The trigonometric function in the integrant can be further expanded using the identity cos ( f + ω ¯ ) = cos f cos ω ¯ sin f sin ω ¯ . Along with the change of variables in Equations (45) to (48), the integral can be written as:
ϵ i ( 1 ) = 0 2 π a ¯ r ¯ 1 e ¯ 2 F ( cos E e ¯ ) cos ω ¯ 1 e ¯ 2 sin E sin ω ¯ d E = a ¯ 2 e ¯ F cos ω ¯ 1 e ¯ 2 0 2 π cos 2 E d E = π a ¯ 2 e ¯ F cos ω ¯ 1 e ¯ 2 = 0
since e ¯ = 0 .

6.2.2. Longitude of Ascending Node

The variation of inclination subject to a small continuous thrust is described in Equation (18). Consider an out-of-plane, small PLP thrust F acting on the mission ship, i.e., F h = F . Integrating Equation (18) over 2 π :
ϵ Ω ( 1 ) = 0 2 π n ¯ a ¯ r ¯ 1 e ¯ 2 F sin ( f + ω ¯ ) csc i ¯ d τ
Similarly, the trigonometric identity and the change of variable bring:
ϵ Ω ( 1 ) = 0 2 π a ¯ r ¯ 1 e ¯ 2 F 1 e ¯ 2 sin E cos ω ¯ + ( cos E e ¯ ) sin ω ¯ csc i ¯ d E = a ¯ 2 e ¯ F sin ω ¯ csc i ¯ 1 e ¯ 2 0 2 π cos 2 E d E = π a ¯ 2 e ¯ F sin ω ¯ csc i ¯ 1 e ¯ 2 = 0

6.2.3. Argument of Periapsis

The variation of inclination subject to a small continuous thrust is described in Equation (20). With an out-of-plane, small PLP force F acting on the mission ship, i.e., F h = F , as the only thrust, the differential equation of the first-order term degenerates to:
d ϵ ω ( 1 ) d τ = cos i ¯ d ϵ Ω ( 1 ) d τ = n ¯ a ¯ r ¯ 1 e ¯ 2 F sin ( f + ω ¯ ) cot i ¯
It is obvious that the integral of Equation (109) over 2 π is:
ϵ ω ( 1 ) = 0 2 π n ¯ a ¯ r ¯ 1 e ¯ 2 F sin ( f + ω ¯ ) cot i ¯ d τ = π a ¯ 2 e ¯ F sin ω ¯ cot i ¯ 1 e ¯ 2 = 0

7. Numerical Simulations

Figure 8, Figure 9, Figure 10, Figure 11, Figure 12, Figure 13, Figure 14, Figure 15, Figure 16, Figure 17, Figure 18, Figure 19, Figure 20, Figure 21, Figure 22, Figure 23 and Figure 24 present numerical simulations of several selected cases. The numerical simulations were performed to verify the algorithms and the results proposed in the proceeding sections. Five types of simulations are presented here.
First, a planar PLP thrust was applied to the mission spacecraft for an orbit period. Since the original launching direction may be arbitrary in real missions, we present four examples in this part, corresponding to four initial directions. Secondly, one of the aforementioned four cases was selected to integrate over a longer time period. This simulation was performed to qualitatively demonstrate that a PLP thrust smaller than the proposed criteria is not powerful enough to expel a mission spacecraft for interplanetary travel, as concluded in our algorithm. Thirdly, an out-of-plane example is presented to qualitatively show the variation of the orbit elements. However, the results derived in the preceding section cannot be applied to determine the bound of variation, because the motions in the out-of-plane cases are always coupled with planar components. Fourthly, a dimensional case is presented to show how our results can be applied. Finally, an interesting example shows that a PLP thrust can control a spacecraft formation by the Clohessy–Wiltshire (C-W) equations.
Throughout all the numerical simulations, the magnitude of the normalized PLP thrust was set as F = 0.008 . The normalized orbit elements of the mother ship were randomly given as ( a m , e m , i m , Ω m , ω m ) = ( 1 , 0 , 30 , 40 , 50 ) . Since the mother ship is in a circular orbit, without lost of generality, the initial true anomaly can be set as zero, i.e., f m ( 0 ) = 0 . All trajectories that examine the variations of the orbit elements were obtained by integrating Equations (4) and (5) directly and then converted to the osculating elements.

7.1. Verification of Bounds with a Planar PLP Thrust

In this section, the PLP force was assumed to be planar. The initial normalized offsets of the mission spacecraft from the mother ship were given as L 0 = 0.1   i , 0.1   j , 0.1   i , and 0.1   j , respectively.
Figure 8, Figure 9, Figure 10, Figure 11, Figure 12, Figure 13, Figure 14 and Figure 15 present the simulation results. All the bounds for orbit elements are discussed in Section 5. In Figure 8 and Figure 9, the initial excursion was set as L 0 = 0.1   i . The PLP force always aims along the line of sight of the mother ship and mission spacecraft. In Figure 8, the trajectories of mother ship and mission spacecraft are presented with respect to the inertial frame, respectively. In the right-most plot, the relative trajectory of the mission spacecraft to the mother ship is provided, and it is also presented in the inertial frame. Figure 9 shows the variation of orbit elements of the mission spacecraft. Since the true anomaly has nothing to do with the variation of the orbit, it is not presented in the figure. Every bound over one normalized orbit period is provided in the title of the corresponded subplot and also presented as a dashed line. It was verified that the bounds were well predicted by our algorithm. Similar results can also be observed in Figure 10, Figure 11, Figure 12, Figure 13, Figure 14 and Figure 15, corresponding to different initial conditions.

7.2. Influence of a Small Planar PLP Thrust

This section verifies the sufficient condition of the lowest thrust required for interplanetary travel. According to the previous results, the mission spacecraft might be trapped in the vicinity of the mother ship if the normalized thrust force is less than 0.008.
Figure 16 and Figure 17 verify the result. In these two simulations, the initial excursion L 0 = 0.1   j was selected. The simulation was run for 200 normalized periods. It is obvious that the mission spacecraft stayed in the vicinity of the mother ship, either from the relative trajectory or the variation of the orbit elements. Notably, the bounds in the titles of Figure 17 are only valid for one normalized period.
Although the fourth case in the previous simulations was selected to demonstrate in this section, similar results can be obtained using different initial offsets.
A counter example is presented in Figure 18 and Figure 19. In this simulation, all parameters of the mission spacecraft were set identical to the previous simulation. There were only two differences. First, the thrust level was set as twice the previous, i.e., F = 0.016 . Secondly, the simulation only ran for 25 normalized period because the trajectory diverged very fast after that.

7.3. Influence of a Small Out-of-Plane PLP Thrust

Figure 20 and Figure 21 present the simulation results of an out-of-plane case. Most parameters were selected identical to the planar cases. The only difference was imposed on the initial offset, setting it as L 0 = 0.1   k . Notably, the results derived in Section 6 cannot be employed to predict the simulation directly. Once the mission spacecraft is expelled by an out-of-plane thrust, the motion will generate a planar component naturally. This enforces the PLP force to have planar components, and the variation of orbit elements should consider both planar and out-of-plane perturbations. As a result, the simulation in this section only qualitatively demonstrates how a small out-of-plane PLP thrust influences the orbit elements of the mission spacecraft. In the simulation, all parameters were identical to the previous simulations.

7.4. Practical Applications

This simulation, presented in Figure 22 and Figure 23, was to verify the proposed algorithms with full-dimensional parameters and study potential practical applications. In this simulation, a mother ship was assumed to be placed in a circular orbit of altitude 400 km above the Earth’s surface. Other orbit elements were identical to the previous simulations. As for the PLP force, we had a laser of power 10 W, and the reflectance of the mirror installed on the mission spacecraft was assumed to be R m = 0.999998 . The mass of the mission spacecraft was assumed to be 1 kg. These three parameters were cited from our previous investigation in [9]. Suppose the mission spacecraft is lifted by a robotic arm for launch. The robotic arm was set as 17 m long, which was inspired by the Canadarm2 robotic arm installed on the International Space Station.
According to these settings, the PLP thrust can be computed as:
F P L P = 2 P R m S c = 0.0334   ( N )
The corresponded normalized force is given by:
F = F P L P R m 2 μ E = 0.0334 × 6778137 2 3.986 × 10 14 = 0.0038
Moreover, the initial offset was assumed to be 17 m, corresponding to a normalized offset of:
l 0 = 17 6778137 = 2.5081 × 10 9
The simulation was run for 200 orbit periods, where one period was:
T = 2 π 6778137 3 3.986 × 10 14 = 5553.6   ( s )
It is obvious, from the simulation results, that the proposed algorithms provided a good prediction on the variation of the orbit elements of the mission spacecraft.
One thing to note is that the bound of ω was unreasonably high. If we look back at Equation (20), a singularity will exist in the variation equation of ω given e = 0 , and large error might exist for e 0 . In our simulation, e ¯ = 2.5081 × 10 9 0 , and the singularity dominated the result. Hence, the bound of ω was actually meaningless in this case.

7.5. Formation Control by C-W Equations

An interesting example presented in this section is to control a spacecraft formation by the C-W equations using PLP thrust. The C-W equations are well known for the linearly approximated relative trajectory to a nominal circular orbit in the rotational frame. Let L = ζ e ^ r + η e ^ θ + ξ e ^ h and L = ζ ˙ e ^ r + η ˙ e ^ θ + ξ ˙ e ^ h be the relative velocity in the rotational frame. The linearized equations of motion with small PLP thrust by the C-W equations are:
ζ ¨ = 3 n R 2 ζ + 2 n R η ˙ + F ζ l ,
η ¨ = 2 n R ζ ˙ + F η l ,
ξ ¨ = n R 2 ξ + F ξ l .
In the normalized system, the equations look similar by setting n R = 1 . The unperturbed C-W equations have closed-form solutions. One of the famous special cases is linear drift, which requires the initial conditions to be L 0 = ( ζ 0 , η 0 , 0 ) and L 0 = ( 0 , 1.5 ζ 0 , 0 ) .
In order to test the robustness of our result, the linear drift case was selected. In the simulation presented in Figure 24, the initial conditions were set as L 0 = ( 0.1 , 0 , 0 ) and L 0 = ( 0 , 0.15 , 0 ) . The simulation was run for five orbit periods. It is clear in the figure that an uncontrolled trajectory, presented by a dashed line, linearly drifts away. The controlled trajectory presented by a solid curve, however, is confined within a small region. One should realize that the controlled trajectory is not a closed curve, though it looks as an ellipse. This example verifies the proposed PLP trust threshold and implies that this type of trust is potentially applicable to the formation flight of spacecraft.

8. Conclusions

This study investigated the variation of the orbit elements of a spacecraft propelled by photonic laser propulsion (PLP) under the two-body problem assumption. This study first reviewed some facts of PLP thrust and background knowledge of celestial mechanics. Perturbation theory was also reviewed and introduced to investigate this problem. Since PLP thrust is continuous and constant, the Gauss equations were introduced to study the variation of osculating orbit elements. In order to generalize the study, all equations were normalized based on the dynamics of the mother ship. The mission spacecraft subject to planar and out-of-plane PLP thrusts was studied, respectively. One should note that in practice, a pure out-of-plane PLP thrust for the whole mission is infeasible. The motion of the spacecraft will naturally cause planar components of the PLP force. As a result, a practical motion subject to initial out-of-plane PLP thrust will eventually be described by the combination of planar and out-of-plane effects. In this study, bounds on the variation of the orbit elements of the mission spacecraft were derived for both small planar and out-of-plane thrusts by zeroth- and first-order approximations. A sufficient condition that traps the mission spacecraft in the vicinity of the mother ship was also found. The study suggested that it may take a very long time period for the mission spacecraft to leave the mother ship if the normalized PLP thrust is less than 0.008. A similar, but more rigorous criterion was proposed in [9]. Our new discovery pushes the bound higher. This threshold was also verified by the linearized dynamics described by the Clohessy–Wiltshire equations. All proposed results, including the bounds and the sufficient conditions, were verified by numerical simulations. Our work presented in this paper is directly applicable to the usage of PLP thrust in future space missions, such as the formation flight about a circular orbit or interplanetary travel.

Funding

This research was partly funded by the Department of Education of Taiwan and partly by the National Science Council of Taiwan through Project NSC-100-2221-E-032-032.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of the data; in the writing of the manuscript; nor in the decision to publish the results.

Appendix A. Approximation of 1/r

As shown in Figure 4, r = 1 + 2 l 0 cos θ + l 0 2 . Given θ = θ 0 , the approximation of a function f ( r ) = f ( l 0 ) can be obtained through the Taylor expansion about the operation point l 0 = 0 . Let f ( r ) = 1 / r = 1 / 1 + 2 l 0 cos θ 0 + l 0 2 = f ( l 0 ) , and the function can be approximated by:
f ( l 0 ) f ( l 0 ) l 0 = 0 + f l 0 l 0 = 0 l 0 + 1 2 ! 2 f l 0 2 l 0 = 0 l 0 2 1 l 0 cos θ 0 + ( 3 cos 2 θ 0 1 ) l 0 2 2
Hence, the first-order approximation is found by 1 / r 1 l 0 cos θ 0 .

References

  1. Tsander, F.K. From a Scientific Heritage, NASA Technical Translation TTF-541; NASA: Washington, DC, USA, 1967. [Google Scholar]
  2. Forward, R.L. Solar photon thrustor. J. Spacecr. Rocket. 1990, 27, 411–416. [Google Scholar] [CrossRef]
  3. Gulevich, A.V.; Ivanov, E.A.; Kukharchuk, O.F.; Poupko, V.Y.; Zrodnikov, A.V. Application of nuclear photon engines for deep-space exploration. AIP Conf. Proc. 2001, 552, 957–962. [Google Scholar] [CrossRef]
  4. Wickman, J.H. Technology Assessment of Photon Propulsion: How Close Are We? In Proceedings of the AIAA/SAE/ASME Conference, Boston, MA, USA, 11–13 July 1981. AIAA Paper 81-1532. [Google Scholar]
  5. Zhang, Q.; Lubin, P.M.; Hughes, G.B.; Melis, C.; Walsh, K.J. Orbital simulations of laser-propelled spacecraft. In Nanophotonics and Macrophotonics for Space Environments IX; Taylor, E.W., Cardimona, D.A., Eds.; International Society for Optics and Photonics, SPIE: Bellingham, WA, USA, 2015; Volume 9616, pp. 61–72. [Google Scholar] [CrossRef] [Green Version]
  6. Meyer, T.R.; Pryor, W.R.; McKay, C.P.; McKenna, P.M. Laser Elevator: Momentum Transfer Using an Optical Resonator. J. Spacecr. Rocket. 2002, 39, 258–266. [Google Scholar] [CrossRef]
  7. Bae, Y. Photonic Laser Propulsion: Proof-of-Concept Demonstration. J. Spacecr. Rocket. 2008, 45, 153–155. [Google Scholar] [CrossRef]
  8. Jessica, F. Stephen Hawking, Mark Zuckerberg, Yuri Milner Launch $ 100M Space Project Called Breakthrough Starshot. 2016. Available online: https://www.natureworldnews.com/articles/20799/20160414/stephen-hawking-mark-zuckerberg-and-russian-millionaire-yuri-milner-launch-100m-space-project-called-breakthrough-starshot.htm (accessed on 14 April 2016).
  9. Hsiao, F.Y.; Wang, T.H.; Lee, C.C.; Yang, P.C.; Wei, H.C.; Tsai, P.Y. Trajectory of spacecraft with photonic laser propulsion in the two-body problem. Acta Astronaut. 2013, 84, 215–226. [Google Scholar] [CrossRef]
  10. NASA. Deep Space 1 (1998-061ADS1). 2008. Available online: https://nssdc.gsfc.nasa.gov/nmc/spacecraft/display.action?id=1998-061A (accessed on 14 April 2016).
  11. ESA. ESA and ANU Make Space Propulsion Breakthrough. 2006. Available online: http://www.spaceref.com/news/viewpr.html?pid=18706 (accessed on 11 January 2006).
  12. Hsiao, F.Y. Trajectory Evolution Under Laser Photonic Propulsion in The Two-Body Problem. In Advances in the Astronautical Sciences; Univelt, Inc.: Escondido, CA, USA, 2013; Volume 148, pp. 13–377. [Google Scholar]
  13. Danby, J.M.A. Fundamentals of Celestial Mechanics, 2nd ed.; Willmann-Bell: New York, NY, USA, 1992. [Google Scholar]
  14. Batrakov, Y.; Efroimsky, M.; Shefer, V. Nonosculating, Osculating, and Superosculating Orbits in Celestial Mechanics. 2007. Available online: https://www.researchgate.net/publication/273966361_Nonosculating_osculating_and_superosculating_orbits_in_celestial_mechanics (accessed on 11 July 2007).
  15. Efroimsky, M.; Goldreich, P. Gauge Symmetry of the N-body Problem in the Hamilton-Jacobi Approach. J. Math. Phys. 2003, 44, 5958. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Photon thrust data obtained with an output coupler mirror with a reflectance of 0.99967 [7].
Figure 1. Photon thrust data obtained with an output coupler mirror with a reflectance of 0.99967 [7].
Aerospace 09 00075 g001
Figure 2. Diagram of the photonic laser propulsion system on a spacecraft [6].
Figure 2. Diagram of the photonic laser propulsion system on a spacecraft [6].
Aerospace 09 00075 g002
Figure 3. Relative positions between the central body, mother ship, and mission spacecraft.
Figure 3. Relative positions between the central body, mother ship, and mission spacecraft.
Aerospace 09 00075 g003
Figure 4. Parameters of the mother mother ship and mission spacecraft. All parameters are normalized in this figure.
Figure 4. Parameters of the mother mother ship and mission spacecraft. All parameters are normalized in this figure.
Aerospace 09 00075 g004
Figure 5. The maximum value of K over E in Equation (63) as a function of eccentricity e ¯ .
Figure 5. The maximum value of K over E in Equation (63) as a function of eccentricity e ¯ .
Aerospace 09 00075 g005
Figure 6. The maximum value of K ω over E as a function of eccentricity e ¯ .
Figure 6. The maximum value of K ω over E as a function of eccentricity e ¯ .
Aerospace 09 00075 g006
Figure 7. Parameters of the mother ship and mission spacecraft subject to out-of-plane PLP thrust. All parameters are normalized in this figure.
Figure 7. Parameters of the mother ship and mission spacecraft subject to out-of-plane PLP thrust. All parameters are normalized in this figure.
Aerospace 09 00075 g007
Figure 8. (Left) Trajectory of the mother ship; (middle) trajectory of the mission spacecraft; (right) relative trajectory of the mission spacecraft to the mother ship. The initial excursion is L 0 = 0.1   i .
Figure 8. (Left) Trajectory of the mother ship; (middle) trajectory of the mission spacecraft; (right) relative trajectory of the mission spacecraft to the mother ship. The initial excursion is L 0 = 0.1   i .
Aerospace 09 00075 g008
Figure 9. Variation of the orbit elements of the mission spacecraft subject to PLP thrust F = 0.008 with the initial excursion being L 0 = 0.1   i . The units for ( i , Ω , ω ) are degrees, whereas a is dimensionless.
Figure 9. Variation of the orbit elements of the mission spacecraft subject to PLP thrust F = 0.008 with the initial excursion being L 0 = 0.1   i . The units for ( i , Ω , ω ) are degrees, whereas a is dimensionless.
Aerospace 09 00075 g009
Figure 10. (Left) Trajectory of the mother ship; (middle) trajectory of the mission spacecraft; (right) relative trajectory of the mission spacecraft to the mother ship. The initial excursion is L 0 = 0.1   j .
Figure 10. (Left) Trajectory of the mother ship; (middle) trajectory of the mission spacecraft; (right) relative trajectory of the mission spacecraft to the mother ship. The initial excursion is L 0 = 0.1   j .
Aerospace 09 00075 g010
Figure 11. Variation of the orbit elements of the mission spacecraft subject to PLP thrust F = 0.008 with the initial excursion being L 0 = 0.1   j . The units for ( i , Ω , ω ) are degrees, whereas a is dimensionless.
Figure 11. Variation of the orbit elements of the mission spacecraft subject to PLP thrust F = 0.008 with the initial excursion being L 0 = 0.1   j . The units for ( i , Ω , ω ) are degrees, whereas a is dimensionless.
Aerospace 09 00075 g011
Figure 12. (Left) Trajectory of the mother ship; (middle) trajectory of the mission spacecraft; (right) relative trajectory of the mission spacecraft to the mother ship. The initial excursion is L 0 = 0.1   i .
Figure 12. (Left) Trajectory of the mother ship; (middle) trajectory of the mission spacecraft; (right) relative trajectory of the mission spacecraft to the mother ship. The initial excursion is L 0 = 0.1   i .
Aerospace 09 00075 g012
Figure 13. Variation of the orbit elements of the mission spacecraft subject to PLP thrust F = 0.008 with the initial excursion being L 0 = 0.1   i . The units for ( i , Ω , ω ) are degrees, whereas a is dimensionless.
Figure 13. Variation of the orbit elements of the mission spacecraft subject to PLP thrust F = 0.008 with the initial excursion being L 0 = 0.1   i . The units for ( i , Ω , ω ) are degrees, whereas a is dimensionless.
Aerospace 09 00075 g013
Figure 14. (Left) Trajectory of the mother ship; (middle) trajectory of the mission spacecraft; (right) relative trajectory of mission spacecraft to the mother ship. The initial excursion is L 0 = 0.1   j .
Figure 14. (Left) Trajectory of the mother ship; (middle) trajectory of the mission spacecraft; (right) relative trajectory of mission spacecraft to the mother ship. The initial excursion is L 0 = 0.1   j .
Aerospace 09 00075 g014
Figure 15. Variation of the orbit elements of the mission spacecraft subject to PLP thrust F = 0.008 with the initial excursion being L 0 = 0.1   j . The units for ( i , Ω , ω ) are degrees, whereas a is dimensionless.
Figure 15. Variation of the orbit elements of the mission spacecraft subject to PLP thrust F = 0.008 with the initial excursion being L 0 = 0.1   j . The units for ( i , Ω , ω ) are degrees, whereas a is dimensionless.
Aerospace 09 00075 g015
Figure 16. (Left) Trajectory of the mother ship; (middle) trajectory of the mission spacecraft; (right) relative trajectory of the mission spacecraft to the mother ship. The initial excursion is L 0 = 0.1   j . This simulation was run for 200 normalized periods.
Figure 16. (Left) Trajectory of the mother ship; (middle) trajectory of the mission spacecraft; (right) relative trajectory of the mission spacecraft to the mother ship. The initial excursion is L 0 = 0.1   j . This simulation was run for 200 normalized periods.
Aerospace 09 00075 g016
Figure 17. Variation of the orbit elements of the mission spacecraft subject to PLP thrust F = 0.008 with the initial excursion being L 0 = 0.1   j . This simulation was run for 200 normalized period. The units for ( i , Ω , ω ) are degrees, whereas a is dimensionless. Notably, the bounds in the titles are only valid for 1 normalized period.
Figure 17. Variation of the orbit elements of the mission spacecraft subject to PLP thrust F = 0.008 with the initial excursion being L 0 = 0.1   j . This simulation was run for 200 normalized period. The units for ( i , Ω , ω ) are degrees, whereas a is dimensionless. Notably, the bounds in the titles are only valid for 1 normalized period.
Aerospace 09 00075 g017
Figure 18. (Left) Trajectory of the mother ship; (middle) trajectory of the mission spacecraft; (right) relative trajectory of the mission spacecraft to the mother ship. The initial excursion is L 0 = 0.1   j . This simulation was run for 25 normalized periods. Notably, the bounds in the titles are only valid for 1 normalized period.
Figure 18. (Left) Trajectory of the mother ship; (middle) trajectory of the mission spacecraft; (right) relative trajectory of the mission spacecraft to the mother ship. The initial excursion is L 0 = 0.1   j . This simulation was run for 25 normalized periods. Notably, the bounds in the titles are only valid for 1 normalized period.
Aerospace 09 00075 g018
Figure 19. Variation of the orbit elements of the mission spacecraft subject to PLP thrust F = 0.016 with the initial excursion being L 0 = 0.1   j . This simulation was run for 25 normalized period. The units for ( i , Ω , ω ) are degrees, whereas a is dimensionless.
Figure 19. Variation of the orbit elements of the mission spacecraft subject to PLP thrust F = 0.016 with the initial excursion being L 0 = 0.1   j . This simulation was run for 25 normalized period. The units for ( i , Ω , ω ) are degrees, whereas a is dimensionless.
Aerospace 09 00075 g019
Figure 20. (Left) Trajectory of the mother ship; (middle) trajectory of the mission spacecraft; (right) relative trajectory of the mission spacecraft to the mother ship. The initial excursion is L 0 = 0.1   k .
Figure 20. (Left) Trajectory of the mother ship; (middle) trajectory of the mission spacecraft; (right) relative trajectory of the mission spacecraft to the mother ship. The initial excursion is L 0 = 0.1   k .
Aerospace 09 00075 g020
Figure 21. Variation of the orbit elements of the mission spacecraft subject to PLP thrust F = 0.008 with the initial excursion being L 0 = 0.1   k .
Figure 21. Variation of the orbit elements of the mission spacecraft subject to PLP thrust F = 0.008 with the initial excursion being L 0 = 0.1   k .
Aerospace 09 00075 g021
Figure 22. (Left) Trajectory of the mother ship; (middle) trajectory of the mission spacecraft; (right) relative trajectory of the mission spacecraft to the mother ship. The initial excursion is L 0 = 2.5081 × 10 9   i . This simulation was run for 200 orbit period.
Figure 22. (Left) Trajectory of the mother ship; (middle) trajectory of the mission spacecraft; (right) relative trajectory of the mission spacecraft to the mother ship. The initial excursion is L 0 = 2.5081 × 10 9   i . This simulation was run for 200 orbit period.
Aerospace 09 00075 g022
Figure 23. Variation of the orbit elements of the mission spacecraft subject to PLP thrust F = 0.0038 with the initial excursion being L 0 = 2.5081 × 10 9   i . This simulation was run for 200 orbit period. The units for ( i , Ω , ω ) are degrees, whereas a is in kilometer. Notably, the bounds in the titles are only valid for 1 orbit period.
Figure 23. Variation of the orbit elements of the mission spacecraft subject to PLP thrust F = 0.0038 with the initial excursion being L 0 = 2.5081 × 10 9   i . This simulation was run for 200 orbit period. The units for ( i , Ω , ω ) are degrees, whereas a is in kilometer. Notably, the bounds in the titles are only valid for 1 orbit period.
Aerospace 09 00075 g023
Figure 24. Variation of the orbit elements of the mission spacecraft subject to PLP thrust F = 0.008 with the initial excursion being L 0 = 0.1   i . This simulation was run for 5 orbit period. ζ and η are in normalized length.
Figure 24. Variation of the orbit elements of the mission spacecraft subject to PLP thrust F = 0.008 with the initial excursion being L 0 = 0.1   i . This simulation was run for 5 orbit period. ζ and η are in normalized length.
Aerospace 09 00075 g024
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hsiao, F.-Y. Variation of Osculating Orbit Elements Using Low-Thrust Photonic Laser Propulsion in the Two-Body Problem. Aerospace 2022, 9, 75. https://0-doi-org.brum.beds.ac.uk/10.3390/aerospace9020075

AMA Style

Hsiao F-Y. Variation of Osculating Orbit Elements Using Low-Thrust Photonic Laser Propulsion in the Two-Body Problem. Aerospace. 2022; 9(2):75. https://0-doi-org.brum.beds.ac.uk/10.3390/aerospace9020075

Chicago/Turabian Style

Hsiao, Fu-Yuen. 2022. "Variation of Osculating Orbit Elements Using Low-Thrust Photonic Laser Propulsion in the Two-Body Problem" Aerospace 9, no. 2: 75. https://0-doi-org.brum.beds.ac.uk/10.3390/aerospace9020075

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