Next Article in Journal
Decoupling Analysis of Water Consumption and Economic Growth in Tourism in Arid Areas: Case of Xinjiang, China
Previous Article in Journal
Sustainable Higher Education: From E-Learning to Smart Education
Previous Article in Special Issue
Exergy Transfer Analysis of Biomass and Microwave Based on Experimental Heating Process
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Target Selection for a Space-Energy Driven Laser-Ablation Debris Removal System Based on Ant Colony Optimization

1
School of Mechatronics Engineering, Harbin Institute of Technology, Harbin 150001, China
2
Beijing Institute of Spacecraft Environment Engineering, Beijing 100091, China
*
Author to whom correspondence should be addressed.
Sustainability 2023, 15(13), 10380; https://doi.org/10.3390/su151310380
Submission received: 10 May 2023 / Revised: 27 June 2023 / Accepted: 27 June 2023 / Published: 30 June 2023
(This article belongs to the Special Issue Renewable Energy: Sources, Conversion and Utilization)

Abstract

:
The space-energy driven laser-ablation debris removal technology can remove or detach multiple centimeter-level space debris in a single mission. However, the space-energy driven platform can only rely on its own equipment capabilities to detect and identify space debris. It is necessary to select multiple potentially removable debris targets to improve the removal efficiency. In this paper, target selection for a space-energy driven laser-ablation debris removal system is analyzed based on ant colony optimization. The intersection and interaction periods were given by the optimal driving sequence calculation for multiple debris. Parameters such as the detection range, pulsed energy, repetition frequency of the laser and trajectory of debris have been considered as inputs of the simulation. Target selection and optimal action time have been calculated when a single debris entered the detection range of the laser system. This optimization can significantly improve the overall efficiency and laser energy utilization of the space-based laser platform for the same randomly generated debris group, compared to the mode driven sequentially according to the order of entering the laser action range. The results showed that after being filtered by the ant colony algorithm, the number of removable debris doubled, and the de-orbit altitude increased by 15.9%. The energy utilization rate of the laser removal system has been improved by 74.6%. This optimization algorithm can significantly improve the overall work efficiency and laser energy utilization rate of the space-energy driven system. It can remove more debris or have a larger effective orbit reduction distance value for all debris.

1. Introduction

Space debris is a non-functional artificial object and its components operate in orbit or re-enter the atmosphere on Earth. With the increasing frequency of space activities, the amount of space debris has rapidly increased, posing a great threat to the safe operation of in-orbit assets and the conduct of space activities [1]. The international community has continued to pay attention to this issue, and countries around the world have reached a broad consensus on jointly carrying out space debris environmental governance [2]. The measures taken for post-mission disposal and collision avoidance of space debris are insufficient to suppress the continuous deterioration of the space debris environment. The current international consensus is that the active removal of space debris is an inevitable choice and the only way to ensure the long-term sustainability of space [3,4,5].
Currently, most active space debris removal technologies, including robotic arm capture, fly net capture, harpoons, etc., are mainly aimed at large-sized debris above the meter level and are difficult to apply to the removal of space debris on a scale of 1–10 cm [6,7,8]. However, there is an urgent need to remove debris on a scale of 1–10 cm. On the one hand, 1–10 cm debris cannot be actively avoided or passively protected at the current technological level, which is the main factor causing fatal damage to spacecrafts [9,10,11]. Targeted active removal is needed to ensure the safety of spacecraft in orbit. On the other hand, due to the large amount and the frequent collisions with each other, debris can generate secondary debris when they collide [12,13,14]. It is necessary to actively remove them to suppress the continuous deterioration of the space debris environment.
Space-energy driven laser-ablation debris removal involves using a strong laser beam to irradiate the surface of the debris, causing the melting, vaporization, and ionization of the material in the irradiation zone and forming of a plasma plume [15,16]. The impulse coupling of the plasma plume causes the debris to obtain a reverse velocity increment, thereby driving the movement of the debris, changing its orbit, and lowering it into the atmosphere to achieve removal. This technology has the characteristics of simple operation, long action distance, rapid response, low cost, and high reliability [17]. It is especially suitable for the removal of 1–10 cm scale debris.
Due to the difficulty in tracking and cataloging 1–10 cm debris, the potential target cannot be predicted and must be detected independently by the laser platform [18]. It is hard to develop removal strategies in advance like the mechanical capture debris removal method. During the removal period, the laser platform may discover multiple space debris at the same time. Each of the detected debris may enter the ablation range, with a limited time window and potential conflicts [19]. It is necessary to analyze, select and plan tasks in real-time based on the detected debris.
The optimization of the selection and removal sequence of multiple debris during laser removal is essentially the path optimization of the removal system. There are many mature algorithms for path optimization. Ant colony optimization algorithms have good robustness and excellent global search ability, which can quickly achieve the global optimization [20,21,22].The core idea of ant colony optimization is the self-organizing collective behavior of ants during foraging. It means that the shortest distance between the food source and the nest can be found through collective search behavior. When ants move between the food source and the nest, they will release a volatile pheromone along the way, and other ants will follow the path of the released pheromone, and eventually the ants can converge to the same path. At present, the max/min ant colony optimization can avoid premature convergence and achieve better performance, making it the best ant colony optimization algorithm to solve problems such as the traveling salesman and quadratic allocation [23].
In this paper, mission planning based on the typical task flow of space-based laser driven space debris removal is introduced. Single and multi debris de-orbit optimization control algorithms were established based on the ant colony algorithm. The computational effectiveness was verified using the debris generated near the platform orbit based on the platform orbit information as an input.

2. Target Selection Model

Due to the requirement of tracking accuracy, space-based laser platforms cannot perform laser irradiation under the guidance of other platforms. Even with the support of additional detection data, the laser ablation process still needs to rely on the platform’s own tracking and aiming equipment to detect and capture the target. When debris is detected, it is necessary to determine whether it will enter the platform’s range of action based on its trajectory, in order for the debris to be ablated. During the ablation process, the impulse determines whether the interaction can effectively lower the orbit or propel the debris away. The impulse direction generated by laser ablation is determined by the relative positions of the debris and laser platform. Debris that enters the scope of ablation has a limited time window, during which effective trajectory changes can be achieved through propulsion. The process is illustrated in Figure 1. When multiple debris enter the scope of action simultaneously or sequentially, challenges arise regarding the order of ablation and selection in the presence of time conflicts. This requires further optimization of the optimal action sequence based on their respective time windows and expected driving effects.
Based on the introduction above, the process of a space-energy driven laser-ablation debris removal system is as follows:
  • n targets T1−Tn were found at time t0−t1.
  • Calculate and filter the discovered debris to eliminate those that do not fall within the scope of the platform.
  • Single target selection optimization. Calculate the relative position change of all targets that will enter the scope of action one by one according to their estimated tracks. Analyze their optimal action time when they can be effectively cleared (perigee is less than 200 km) to obtain the optimal driving time.
  • Multi objective sequence optimization. With the goal of achieving the best overall driving effect throughout the entire time period, analyze the time conflict situation of the action window for all targets entering the range of action, resolve conflicts, and form a driving time sequence.
  • After the target enters the range, it is sequentially driven according to the sequence.

2.1. Centimeter Space Debris Removal Scene

Due to the lack of a cataloging database for centimeter-level space debris, it is necessary to construct a simulation calculation scenario that is as close to the actual situation as possible based on the centimeter-level space debris data near the laser platform orbit. The generated debris is the calculation input for sequence optimization. In this article, based on the orbital parameters of the debris removal platform, the Monte Carlo random sampling method is used to construct the orbital elements of the initial space debris near the production platform orbit in Equations (1) and (2) [24].
R m = R s c + d R m = R s c + d R m [ cos E l e cos A z cos E l e sin A z sin E l e ]
V m = V m [ cos α sin β cos α cos β sin α ]
Az is the azimuth angle of the debris relative to the removal system, Ele is the elevation angle of the debris relative to the removal system, dRm is the relative position vector of debris m relative to the laser platform in the inertial frame, where dRm meets the detection distance requirement [25]. Vm is the velocity vector of debris m in the inertial frame, where
V m [ μ R max , μ R min ]
Based on Equation (3), the initial space debris dataset can be obtained by using the Monte Carlo random sampling method. It can be sorted according to the distribution characteristics of space debris; for example, the orbital inclination of debris is less than 100° and the eccentricity is less than 0.05. If the perigee of debris is greater than 200 km, the space debris database can be obtained.
The coordinate transformation matrix from the inertial coordinate system to the orbital coordinate system (RSW) is in Equation (4) [26]:
R o i = R z ( u ) R x ( i ) R z ( Ω ) = [ cos u cos Ω sin u cos i sin Ω cos u sin Ω + sin u cos i cos Ω sin u sin i sin u cos Ω cos u cos i sin Ω sin u sin Ω + cos u cos i cos Ω cos u sin i sin i sin Ω sin i cos Ω cos i ] [ x i y i z i ] = R i o [ r 0 0 ] = a ( 1 e 2 ) 1 + e cos f [ cos u cos Ω sin u cos i sin Ω cos u sin Ω + sin u cos i cos Ω sin u sin i ] = r [ cos δ cos α cos δ sin α sin δ ]
δ and α are the declination and right ascension of the debris, respectively. If the positions of two space debris are equivalent, it is obvious in Equation (5) that:
a s ( 1 e s 2 ) 1 + e s cos f s = a ( 1 e 2 ) 1 + e cos f [ cos u s cos Ω s sin u s cos i s sin Ω s cos u s sin Ω s + sin u s cos i s cos Ω s sin u s sin i s ] = [ cos u cos Ω sin u cos i sin Ω cos u sin Ω + sin u cos i cos Ω sin u sin i ]
If the initial position of the laser platform is near the ascending/descending intersection point, Monte Carlo random sampling results in a probability of generating debris that meet the requirements of around 4%. If there is a chance of a close intersection for two space debris with different orbital inclinations, the time window of the leap is near to the ascending or descending intersection point. According to Formula (3), it is necessary to meet the following conditions in Equation (6) [27]:
a s ( 1 e s 2 ) 1 + e s cos f s = a ( 1 e 2 ) 1 + e cos f [ cos u s cos Ω s sin u s cos i s sin Ω s cos u s sin Ω s + sin u s cos i s cos Ω s sin u s sin i s ] = [ cos Ω sin Ω 0 ]

2.2. Single Debris Selection

The single debris optimization control algorithm needs to target single space debris. Based on the input of space debris, laser system and optimal de-orbit effect, the laser ablation sequence can be determined. The laser ablation sequence includes start point, end point, total period of the ablation and total pulse number of the laser system. The optimal sequence is the shortest working time of the laser when two sequences can achieve the same de-orbit reduction effect.

2.2.1. Laser Impulse Calculation

According to the information of the space debris and laser system, the maximum de-orbit effect of debris can be calculated. If the perigee of space debris is lower than 200 km, it will be seen as successful removal.
Simplify space debris as a planar target for laser impulse calculation, with the following formula in Equation (7) [28]:
f L τ L = I 0 m = C m I L τ L A L m cos ( n ^ , L ^ R ) n ^
where Cm is the impulse coupling coefficient, IL is the laser power density, τ L is the laser pulse width, AL is the laser irradiation area, n is the normal vector irradiated by a laser, L R is the direction of laser irradiation. The conditions for the target material to be irradiated by the laser are in Equation (8):
cos ( n ^ , L ^ R ) < 0
In this study, it is assumed that the space debris is a circular aluminum alloy thin plate. The plane normal vector irradiated by the laser is in the opposite direction to the laser irradiation direction. We make this assumption because in current detection methods it is difficult to obtain centimeter-level debris dimensions and material information at a long distance.

2.2.2. Orbital Motion Model

Based on the relative motion between the space debris and laser platform, the intersection process between debris and the platform can be obtained. For space debris (perigee below 400 km), the long-term impact of atmospheric drag on orbit must be considered, and the formula for calculating atmospheric drag is in Equation (9) [29]:
f D = C D A 2 m ρ V V V = v ω e × r
where CD ≈ 2.2 is the resistance coefficient, A/m is the surface to mass ratio of space debris, ωe is the rotational velocity of the Earth. According to the fitting formula of American standard atmosphere model, the atmospheric density of 200–600 km can be calculated by the formula below in Equation (10):
ln ( ρ ρ m ) = 9.595 9.7875 × 10 3 ( h h u 200 ) + 7.0725 × 10 6 ( h h u 200 ) ( h h u 400 )
Under the combined action of laser ablation force and atmospheric resistance, the movement of space debris can be described by position velocity in Equation (11) [30]:
r ¨ s = g ( r s ) + f D , s r ¨ d = g ( r d ) + f D , d + f L ( r d , r s ) = g ( r d ) + f D , d + F max u n
Compared with the average laser ablation force, the impact of atmospheric drag can be ignored in sequence planning when the altitude of the space debris is higher than 400 km. At this point, the equation can be expressed as in Equations (12)–(14) [31]:
{ r ˙ s = v s v ˙ s = μ r s r s 3 r ˙ d = v d v ˙ d = μ r d r d 3 + F max u n
F max = { 0 i f r d r s > r min C m E f L / M e l s e
n = r d r s r d r s
where u [ 0 , 1 ] is the on/off condition of the laser system.
Under the σ ( t ) = h p , min ϕ ( t ) = h p , min h p ( r d , v d ) 0 condition, an appropriate optimal control rate to minimize the debris perigee, or minimize the startup time when the perigee height is less than 200 km can be given. The minimum performance index function can be described as in Equations (15) and (16) [32]:
J = h p ( t f ) + λ 0 t 0 t f u F max d t = [ a ( 1 e ) R E ] + t 0 t f λ 0 u F max d t = ϕ ( r d ( t f ) , v d ( t f ) ) + t 0 t f L ( r d , v d , u F max ) d t
where
a = ( 2 r d ( t f ) v d ( t f ) 2 μ ) 1 , e = [ 1 r d ( t f ) a ] 2 + r d ( t f ) 2 v d ( t f ) 2 μ a
This is a typical optimal control problem with state space constraints. When solving based on the indirect method in optimal control, the Lagrangian function is in Equations (17) and (18) [33]:
J α = J + t 0 t f κ σ d t = ϕ + t 0 t f ( L + κ σ ) d t
ψ = H + κ σ = L + λ x ˙ + κ σ
The Hamiltonian function is in Equation (19):
H = λ x ˙ + L ( x , u , t ) = λ r v d + λ v ( μ r d r d 3 + F max u n ) + λ 0 u F max = λ r v d λ v ( μ r d r d 3 ) + ( λ 0 + λ v n ) F max u λ ˙ = H x κ σ x λ ˙ r = μ r d 3 ( λ v 3 ( λ v r d ) r d 2 r d ) κ σ r d = μ r d 3 ( λ v 3 ( λ v r d ) r d 2 r d ) + κ ϕ r d λ ˙ v = λ r κ σ v d = λ r + κ ϕ v d
λ0 is the weight coefficient of perigee height and laser power on time.
According to the optimal control rate, the Hamilton function can be minimized as follows. When u is a continuous function in Equation (20),
H u = 0
when u is discontinuous in Equation (21),
u * = arg u [ 0 , 1 ] min H
That means u = 1 , i f ( ρ < 0 ) , the switch function is ρ = ( λ 0 + λ v n ) F max , n = r d r s r d r s . At the same time, the relaxation conditions κ σ = 0 , κ 0 , σ 0 are met. This indicates that the clearance height constraint σ = 0 is reached, the perigee will not change and no force u = 0 will be applied. Otherwise, κ = 0 . At the same time, the terminal covariate variables need to meet Equation (22):
{ λ r ( t f ) = ϕ r d ( t f ) λ v ( t f ) = ϕ v d ( t f )
The solution is divided into two steps. First, solve the optimal control problem of the lowest perigee altitude κ = 0 to obtain the time window for orbit descent. For space debris that can be cleared, the laser on/off time is optimized to minimize the laser working period.

2.3. Multiple Debris Selection

To solve the combinatorial optimization problem, expressing it in a standard format is the first step. Then, use the ant colony algorithm to determine the decision point according to the pheromone as the feedback carrier between exploration and utilization. Construct the pheromone of each ant individual incrementally according to the corresponding pheromone-updated rules. Plan the behavior direction of the ant colony activity from the overall perspective. Work out the optimal solution of the combinatorial optimization problem from the cycle. The process is shown in Figure 2.
The continuous ant colony algorithm (CACA) is a method that extends traditional ant colony algorithms to continuous domains. Theoretically, the global optimal obtained through this algorithm. Due to its proximity to the original formula of ant colony algorithm, ACOR has the potential to solve mixed discrete–continuous optimization problems [34]. But, this article mainly studies its use to solve pure continuous optimization problems.
The core idea of ant colony algorithms is to select solutions based on a biased probability and gradually build a complete solution set. In each step, ants choose the composition of the solution from the feasible solution by the probability property according to the formula in Equation (23) [35]:
p ( c i j | s p ) = τ i j α η ( c i j ) β c i l N ( s p ) τ i l α η ( c i l ) β , c i j N ( s p ) P i ( j ) = c i l N ( s p ) τ i l α η ( c i l ) β
Each element of the feasible solution and its probability constitute the discrete probability distribution of the ant sampling, so as to select the elements to be added to the current solution space Cij is the j-th element of the solution constructed by ant i, which is heuristic information and is generally chosen as the reciprocal of the objective function. The pheromone τ i j is updated as follows in Equation (24) [36]:
τ i j = { ( 1 ρ ) τ i j + ρ Δ τ if τ i j s c h ( 1 ρ ) τ i j otherwise
Pheromone Q / L k is usually updated with global information Δ τ , which has a good performance in solving TSP problems.
Referring to the idea of discrete probability distribution in the ant colony algorithm, ACOR proposes a continuous probability distribution for variable optimization in a continuous domain. It is further transformed into sampling in the Gaussian kernel probability density function GI in Equation (25) [37]:
G i ( x ) = l = 1 k ω l g l i ( x ) = l = 1 k ω l 1 σ l i e ( x μ l i ) 2 2 σ l i 2 μ l i = s l i ω l = 1 q k 2 π e ( l 1 ) 2 2 q 2 k 2 σ l i = ξ e = 1 k | s e i s l i | k 1
where k is the number of stored solutions (which cannot be less than the number of variables), q and ξ are the setting parameters of the algorithm.
The multi debris optimization algorithm needs to calculate the intersection process between debris and the laser system. Time to enter the ablation range, time to leave, time of the intersection window and global de-orbit effect all have to be considered as the optimization goal. The core issue is to resolve conflicts between space debris with conflicting time windows based on single space debris optimization. During the process, it is necessary to consider the rotation time of the laser system attitude. There should be a sufficient time interval between adjacent space debris based on their different positions.
Tabu lists provide important information when making choices in ant colony algorithms [38]. They can be used to solve conflict resolution. In order to achieve decoupling of sequence optimization and on-off time optimization for better optimization efficiency, it is necessary to determine whether the time series between adjacent space debris conflict. The resolution strategy needs to be deterministic to ensure that the sequence and the on/off time can be determined.
The time series relationships between adjacent debris can be divided into four types.
  • s t o p i 1 + T s e i b s e i , if there is no conflict, debris i should be added to the feasible solution.
  • b s e i 1 < b s e i T s e i < s t o p i 1 , while s t o p i T s e i b s e i T s e i s t o p i 1 , a conflict resolution algorithm needs to be used.
  • b s e i T s e i < b s e i 1 ( < s t o p i 1 ) a n d s t o p i 1 < s t o p i T s e i , some conflicts require the use of conflict resolution algorithms.
  • b s e i T s e i < b s e i 1 a n d s t o p i 1 s t o p i T s e i , if there is a complete conflict and no solution, then debris i should be added to the tabu table.
Among them, numerical integration is required to calculate the benefits g a i n i = f ( b s e i , s t o p i ) . Steering time T s e i = g ( s t o p i 1 , a r r i ) is required to be calculated according to the inertial pointing of the two debris and the attitude response speed of the laser system, which needs to meet the following requirements in Equation (26):
w a i t i = bse i a r r i = bse i ( s t o p i 1 + T s e i ) = bse i s t o p i 1 g ( s t o p i 1 , bse i ) 0
At the initial stage w a i t i < 0 , it is necessary to adjust b s e i (Strategy 1) or s t o p i 1 (Strategy 2) to meet T s e i = g ( s t o p i 1 , b s e i ) = b s e i s t o p i 1 . An iterative algorithm is used to solve the problem. Among them, the strategy selection adopts a greedy strategy. When the profit of the current debris is more than that of the adjustment, the constraint is satisfied, so select strategy 1.Otherwise, choose strategy 2.
When selecting the next debris according to the state transition probability, this paper comprehensively considers decision factors such as pheromone concentration, debris removal income, removal window length, and waiting time. The calculation formula is as follows in Equation (27) [3]:
p ( c i j | s p ) = τ i j α g a i n ( c i j ) β [ 1 / w i d t h ( c i j ) ] γ [ 1 / w a i t ( c i j ) ] δ c i l N ( s p ) τ i l α g a i n ( c i l ) β [ 1 / w i d t h ( c i l ) ] γ [ 1 / w a i t ( c i l ) ] δ , c i j N ( s p ) P i ( j ) = c i l N ( s p ) τ i l α g a i n ( c i l ) β [ 1 / w i d t h ( c i l ) ] γ [ 1 / w a i t ( c i l ) ] δ
where α , β , γ , δ are the weights of various decision-making factors.

3. Results and Discussion

3.1. Debris Removal Scenarios

The inputs of the simulation are in Table 1. The debris are simplified into circular aluminum alloy sheets with a diameter of 10 cm and a thickness of 1 mm.
There are 20,000 debris generated that can enter the detection range of the platform using the laser system orbit as the input. The distribution of these debris is shown in Figure 3 and Figure 4. After target selection, 118 debris can enter the range of the laser system and have a time window for de-orbit.
The distribution of the relative distance and velocity between the space debris and the laser system are shown in Figure 5 and Figure 6. The results showed that the relative velocity between most of the space debris and the laser system is large, making it difficult to remove.
The results indicated that the intersection of space debris is approximately in the same orbit as the laser system. The orbital inclination and the right ascension of the ascending node are approximately the same. The relative speed is low and the window of time to meet the operating distance is long. Debris can be successfully cleared after a period of laser ablation, but the number of such debris is relatively small. Most of the debris and the laser system have different orbits, with a difference of about 180° in right ascension at the ascending node as shown in Figure 7 and Figure 8. This results in a relatively high relative velocity and a shorter window of time to meet the operating distance, making it difficult to clear.
If the initial position of the platform is not near to the ascending/descending intersection point, Monte Carlo random sampling results in the probability of generating debris that meet the requirements of about 1%.
Taking the initial latitude angle of the laser system u = 120°, the distribution of debris data is shown in Figure 9 and Figure 10. It is evident that when space debris and the laser system are coplanar, the relative velocity is low. In most cases, space debris and the laser system are not coplanar and have relatively high velocities. At this point, sin u s sin i s = sin u sin i and arcsin ( sin u sin i ) i s , u s 180 arcsin ( sin u sin i ) .
Taking the initial latitude angle of the laser platform u = 90°, the distribution of space debris is shown in Figure 11 and Figure 12. Similar to the previous example, in most cases, the relative speed is relatively high.

3.2. Single Debris selection

For debris with a descent window that cannot be cleared, the descent window is the on-off time corresponding to the maximum change in the height of the debris orbit. For debris that can be cleared, ACOR optimization is used to obtain the clearing window corresponding to the shortest ablation time. The output of some of the debris is shown in Table 2. Debris no. 64, 234, and 370 can be effectively removed when the perigee is less than 200 km. Most of the debris-driven windows have conflicts, as shown in Figure 13.

3.3. Multi Debris Selection

For 20 space debris with orbit reduction windows, the ant colony algorithm shown in Table 3 was used to obtain the optimal orbit reduction sequence as shown in Figure 14, and the optimization results are shown in Table 4. The driving sequence is 7-16-18-19. It can be seen that four debris were effectively lowered and two of them were successfully removed.
If optimization is not entered, the regular driving sequence of the laser system is to sequentially drive the debris according to the entry time. It means continuously driving the first debris that enters the scope of the laser, and then moving on to the next debris when it leaves the scope o fthe laser or is effectively removed. The de-orbit effect obtained is shown in Table 5, with a driving sequence of 1-8-15-18. As shown in Figure 14, it can be seen that this scheme only effectively removes one debris from the orbit. From this, it is shown that after target selection and sequence optimization, the number of removable debris has doubled, and the de-orbit effect has increased by 15.9%. The energy utilization rate of the laser system has increased by 74.6%.

4. Conclusions

This study conducts research on the problem of target selection for a space-energy driven laser-ablation debris removal system based on ant colony optimization. Optimization algorithms for single target selection, optimal driving time windows, and multi-objective optimal driving sequences have been established. Calculation and verification have been conducted using the platform’s peripheral debris information, produced by the Monte Carlo method, as samples. The results show that this algorithm can effectively meet the requirements of target selection and sequence optimization for the space-energy driven laser-ablation debris removal system. Compared with the situation without optimization, the number of removable debris is doubled, and the de-orbit effect is increased by 15.9%. The energy utilization rate of the space-energy driven laser-ablation debris removal system has increased by 74.6%. This algorithm is of great significance for the in-orbit engineering application of the space-energy driven laser-ablation debris removal system, especially for the process control of debris removal tasks of large-scale, high-energy, and high-frequency space-based pulse laser platforms.
In this paper, only ant colony optimization algorithms are used to study the screening and path planning of laser removal of space debris, but there are slow rates of convergence and local optimization problems. In the process of laser removal of space debris, there are various requirements such as maximizing the energy utilization rate of the removal system, optimizing the de-orbit effect, fast response ability, and orbital transfer, as well as removal scenarios under various distributions of space debris. In subsequent research, considering the initial state of debris and removal systems, hybrid and meta-acoustics algorithms could be used to study scenarios where a single removal system responds to multiple debris [39,40]. Considering the mutual cooperation and optimization among multiple systems, adaptive algorithms and self-adaptive algorithms could be used to study the scenario of removing space debris from a constellation composed of multiple removal systems [41,42]. Considering the characteristics of debris distribution intervals, island algorithms could be used to study the removal of multiple debris within different orbital ranges by relay [43,44]. In future research, the proposed approaches above could be effectively compared to each other to improve the efficiency of the last removal system.

Author Contributions

W.Y., H.F. and Z.S. designed the study and simulation; W.Y. and C.C. wrote the code; W.Y., Q.W. and C.C. conducted the data analysis; W.Y. and C.C. provided the mathematical methods; W.Y., Z.S. and C.C. drafted the paper; W.Y., H.F., Z.S., Q.W. and C.C. edited the paper. All authors have read and agreed to the published version of the manuscript.

Funding

This work was financially supported by the Special Fund for Space Debris Project (KJSP2020010304).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

A special thanks is expressed to Zizheng Gong for his support with the simulation algorithm.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Alexander, L.; Vladimir, A. Review of contact and contactless active space debris removal approaches. Prog. Aerosp. Sci. 2022, 134, 100858. [Google Scholar]
  2. Shan, M.; Shi, L. Comparison of Tethered Post-Capture System Models for Space Debris Removal. Aerospace 2022, 9, 33. [Google Scholar] [CrossRef]
  3. Barea, A.; Urrutxua, H.; Cadarso, L. Large-scale object selection and trajectory planning for multi-target space debris removal missions. Acta Astronaut. 2020, 170, 289–301. [Google Scholar] [CrossRef]
  4. Leomanni, M.; Bianchini, G.; Garulli, A.; Giannitrapani, A.; Quartullo, R. Orbit control techniques for space debris removal missions using electric propulsion. J. Guid. Control. Dyn. 2020, 43, 1259–1268. [Google Scholar] [CrossRef]
  5. Phipps, C.R.; Boustie, M.; Chevalier, J.M.; Baton, S.; Brambrink, E.; Berthe, L.; Schneider, M.; Videau, L.; Boyer, S.A.R.; Scharring, S. Laser impulse coupling measurements at 400 fs and 80 ps using the LULI facility at 1057 nm wavelength. J. Appl. Phys. 2017, 122, 193103. [Google Scholar] [CrossRef]
  6. Hakima, H.; Bazzocchi, M.C.F. Low-Thrust Trajectory Design for Controlled Deorbiting and Reentry of Space Debris. In Proceedings of the 2021 IEEE Aerospace Conference, Big Sky, MT, USA, 6–13 March 2021; Volume 50100, pp. 1–10. [Google Scholar]
  7. Jankovic, M.; Yüksel, M.; Babr, M.M.; Letizia, F.; Braun, V. Space debris ontology for ADR capture methods selection. Acta Astronaut. 2020, 173, 56–68. [Google Scholar] [CrossRef]
  8. Phipps, C.R.; Loktionov, E.Y.; Bonnal, C.; Boyer, S.A.E.; Sharaborova, E.; Tahan, G. Is laser space propulsion practical?: Review. Appl. Opt. 2021, 60, 1–11. [Google Scholar] [CrossRef]
  9. Chang, H.; Ye, J.; Zhou, W. Effects of Typical Laser Wavelength on Impulse Coupling Characteristics Ablated by Nanosecond Pulsed Laser with Al Target. J. Propuls. Technol. 2015, 36, 1754–1760. [Google Scholar]
  10. Stefan, S.; Jascha, W.; Eckel, H.A. Laser-based removal of irregularly shaped space debris. Opt. Eng. 2016, 56, 011007. [Google Scholar]
  11. Phipps, C.R.; Luke, J.; Funk, D.; Moore, D.; Glownia, J.; Lippert, T. Laser impulse coupling at 130 fs. Appl. Surf. Sci. 2006, 252, 4838–4844. [Google Scholar] [CrossRef]
  12. Mark, C.P.; Kamath, S. Review of Active Space Debris Removal Methods. Space Policy 2019, 47, 194–206. [Google Scholar] [CrossRef]
  13. Jiang, T.; Cao, Y.F.; Ding, M. Progress of Space Debris Detection Technology. Laser Optoelectron. Prog. 2022, 59, 1415010. [Google Scholar]
  14. Hakima, H.; Emami, M.R. Assessment of active methods for removal of LEO debris. Acta Astronaut. 2018, 144, 225–243. [Google Scholar] [CrossRef]
  15. Scharring, S.; Eisert, L.; Lorbeer, R.A. Momentum predictability and heat accumulation in laser-based space debris removal. Opt. Eng. 2018, 58, 011004. [Google Scholar] [CrossRef] [Green Version]
  16. Tsuno, K.; Wada, S.; Ogawa, T. Laser ablation induced impulse study for removal of space debris mission using small satellite. Appl. Phys. 2022, 128, 121–130. [Google Scholar] [CrossRef]
  17. Zhao, S.; Steindorfer, M.; Kirchner, G. Attitude analysis of space debris using SLR and light curve data measured with single-photon detector. Adv. Space Res. 2020, 65, 1518–1527. [Google Scholar] [CrossRef]
  18. Daisuke, S.; Yasuhiro, Y.; Toshiya, H.; Yuki, I.; Tadanori, F. Contactless attitude control of an uncooperative satellite by laser ablation. Acta Astronaut. 2022, 196, 275–281. [Google Scholar]
  19. Sagnières, L.B.M.; Sharf, I.; Deleflie, F. Simulation of long-term rotational dynamics of large space debris: A TOPEX/Poseidon case study. Adv. Space Res. 2020, 65, 1182–1195. [Google Scholar] [CrossRef]
  20. Ali, M. Interactive fuzzy Bayesian search algorithm: A new reinforced swarm intelligence tested on engineering and mathematical optimization problems. Expert Syst. Appl. 2022, 187, 115954. [Google Scholar]
  21. Ali, M.; Mahsa, M. Differential evolution method integrated with a fuzzy decision-making mechanism and Virtual Mutant agent: Theory and application. Appl. Soft Comput. 2021, 112, 107808. [Google Scholar]
  22. Ngoc, T.H.; Tien, V.T.N.; Nguyen, T.T.; Quoc-Manh, N. Optimizing Magnification Ratio for the Flexible Hinge Displacement Amplifier Mechanism Design. In Proceedings of the 2nd Annual International Conference on Material, Machines and Methods for Sustainable Development, Nha Trang, Vietnam, 12–15 November 2020; pp. 769–778. [Google Scholar]
  23. Al-Nimr, M.; Milhem, A.; Al-Bishawi, B.; Al Khasawneh, K. Integrating Transparent and Conventional Solar Cells TSC/SC. Sustainability 2020, 12, 7483. [Google Scholar] [CrossRef]
  24. Nazarenko, A.I.; Usovik, I.V. Space debris in low earth orbits region: Formation and reduction process analysis in past decade. Acta Astronaut. 2022, 194, 210231. [Google Scholar] [CrossRef]
  25. Sarkar, M.; Barad, S.; Mohit, A. Development of a novel autonomous space debris collision avoidance system for uncrewed spacecraft. Proc. Inst. Mech. Eng. Part G J. Aerosp. Eng. 2022, 236, 2940–2952. [Google Scholar] [CrossRef]
  26. Driedger, M.; Asgari, A.; Ferguson, P. Feasibility of Gathering Resident Space Object Range Measurements Using In-Orbit Observers. IEEE J. Radio Freq. Identif. 2022, 6, 250–257. [Google Scholar] [CrossRef]
  27. Obukhov, V.A.; Kirillov, V.A.; Petukhov, V.G. Control of a service satellite during its mission on space debris removal from orbits with high inclination by implementation of an ion beam method. Acta Astronaut. 2022, 194, 390–400. [Google Scholar] [CrossRef]
  28. Carvalho, J.; Moraes, R.D.; Prado, A. Analysis of the orbital evolution of space debris using a solar sail and natural forces. Adv. Space Res. 2022, 70, 125–143. [Google Scholar] [CrossRef]
  29. Herrera-Jaramillo, D.A.; Henao-Bravo, E.E.; González Montoya, D.; Ramos-Paja, C.A.; Saavedra-Montes, A.J. Control-Oriented Model of Photovoltaic Systems Based on a Dual Active Bridge Converter. Sustainability 2021, 13, 7689. [Google Scholar] [CrossRef]
  30. Miao, H.; Zhang, L.; Liu, S.; Zhang, S.; Memon, S.; Zhu, B. Laser Sealing for Vacuum Plate Glass with PbO-TiO2-SiO2-RxOy Solder. Sustainability 2020, 12, 3118. [Google Scholar] [CrossRef] [Green Version]
  31. Tang, J.; Liu, G.; Pan, Q.T. A Review on Representative Swarm Intelligence Algorithms for Solving Optimization Problems: Applications and Trends. IEEE/CAAJ 2021, 8, 1627–1643. [Google Scholar] [CrossRef]
  32. Jorgensen, M.K.; Sharf, I. Optimal planning for a multiple space debris removal mission using high-accuracy low-thrust transfers. Acta Astronaut. 2020, 172, 56–69. [Google Scholar] [CrossRef]
  33. Viavattene, G.; Devereux, E.; Snelling, D.; Payne, N.; Wokes, S.; Ceriotti, M. Design of multiple space debris removal missions using machine learning. Acta Astronaut. 2022, 193, 277–286. [Google Scholar] [CrossRef]
  34. Ajay, K.; Seema, B. Generalized Ant Colony Optimizer: Swarm-based meta-heuristic algorithm for cloud services execution. Computing 2019, 101, 1609–1632. [Google Scholar]
  35. Liu, S.; Zhang, Q.; Zhou, D. Obstacle Avoidance Path Planning of Space Manipulator Based on Improved Artificial Potential Field Method. J. Inst. Eng. 2014, 95, 31–39. [Google Scholar] [CrossRef]
  36. Zhang, T.; Shen, H.; Yang, Y. Ant Colony Optimization-based Design of Multiple-target Active Debris Removal Mission. Trans. Jpn. Soc. Aeronaut. Space Sci. 2018, 61, 201–210. [Google Scholar] [CrossRef] [Green Version]
  37. Stuart, J.; Howell, K.; Wilson, R. Application of multi-agent coordination methods to the design of space debris mitigation tours. Adv. Space Res. 2016, 57, 1680–1697. [Google Scholar] [CrossRef]
  38. Ntagiou, E.V.; Armellin, R.; Iacopino, C.; Policella, N.; Donati, A. Ant-based automated mission planning for data relay space missions. J. Aerosp. Inf. Syst. 2019, 16, 249–262. [Google Scholar] [CrossRef]
  39. Zhao, H.T.; Zhang, C.S. An online-learning-based evolutionary many-objective algorithm. Inf. Sci. 2020, 509, 1–21. [Google Scholar] [CrossRef]
  40. Maxim, A.D. An Adaptive Polyploid Memetic Algorithm for scheduling trucks at a cross-docking terminal. Inf. Sci. 2021, 565, 390–421. [Google Scholar]
  41. Junayed, P.; Arriana, L.N.; Amir, M.; Tian, G.D.; Li, Z.W.; Wang, H.; Maxim, A.D. Exact and metaheuristic algorithms for the vehicle routing problem with a factory-in-a-box in multi-objective settings. Adv. Eng. Inform. 2022, 52, 101623. [Google Scholar]
  42. Gholizadeh, H.; Fazlollahtabar, H.; Amir, M.F.F.; Maxim, A.D. Preventive maintenance for the flexible flowshop scheduling under uncertainty: A waste-to-energy system. Environ. Sci. Pollut. Res. 2021, 1–20. [Google Scholar] [CrossRef]
  43. Dulebenets, M.A.; Kavoosi, M.; Abioye, O.; Pasha, J. A Self-Adaptive Evolutionary Algorithm for the Berth Scheduling Problem: Towards Efficient Parameter Control. Algorithms 2018, 11, 100. [Google Scholar] [CrossRef] [Green Version]
  44. Masoud, R.; Nastaran, O.A.; Niloofar, A.S. Ambulance routing in disaster response considering variable patient condition: NSGA-II and MOPSO algorithms. J. Ind. Manag. Optim. 2022, 18, 1035–1062. [Google Scholar]
Figure 1. Typical task flow for removing debris on a space-based laser platform.
Figure 1. Typical task flow for removing debris on a space-based laser platform.
Sustainability 15 10380 g001
Figure 2. Flowchart of basic ant colony algorithm.
Figure 2. Flowchart of basic ant colony algorithm.
Sustainability 15 10380 g002
Figure 3. Perigee distribution of space debris.
Figure 3. Perigee distribution of space debris.
Sustainability 15 10380 g003
Figure 4. Inclination distribution of space debris.
Figure 4. Inclination distribution of space debris.
Sustainability 15 10380 g004
Figure 5. Relative distance distribution between space debris and laser platform.
Figure 5. Relative distance distribution between space debris and laser platform.
Sustainability 15 10380 g005
Figure 6. Relative velocity distribution between space debris and laser platform.
Figure 6. Relative velocity distribution between space debris and laser platform.
Sustainability 15 10380 g006
Figure 7. Distribution of declination/latitude argument and orbital inclination of right ascension of the ascending node.
Figure 7. Distribution of declination/latitude argument and orbital inclination of right ascension of the ascending node.
Sustainability 15 10380 g007
Figure 8. Distribution of relative speed and declination/latitude argument of right ascension of the ascending node.
Figure 8. Distribution of relative speed and declination/latitude argument of right ascension of the ascending node.
Sustainability 15 10380 g008
Figure 9. Distribution of inclination and argument of latitude/right ascension of the ascending node.
Figure 9. Distribution of inclination and argument of latitude/right ascension of the ascending node.
Sustainability 15 10380 g009
Figure 10. Distribution of relative velocity and orbital inclination/declination of ascending node.
Figure 10. Distribution of relative velocity and orbital inclination/declination of ascending node.
Sustainability 15 10380 g010
Figure 11. Distribution of declination/latitude argument and orbital inclination of ascending node.
Figure 11. Distribution of declination/latitude argument and orbital inclination of ascending node.
Sustainability 15 10380 g011
Figure 12. Distribution of relative velocity and orbital inclination/declination of ascending node.
Figure 12. Distribution of relative velocity and orbital inclination/declination of ascending node.
Sustainability 15 10380 g012
Figure 13. Ablation time window for some debris.
Figure 13. Ablation time window for some debris.
Sustainability 15 10380 g013
Figure 14. Ablation duration before and after optimization.
Figure 14. Ablation duration before and after optimization.
Sustainability 15 10380 g014
Table 1. Laser system configurations.
Table 1. Laser system configurations.
Parameters NameValue
Altitude650 km
Inclination97.7°
Right ascension120°
Detection range200 km
Laser driving range100 km
Attitude pointing turning speed15°·s−1
Laser single pulse energy300 J
Spot diameter15 cm
Pulse width50 nm
Repetition frequency100 Hz
Table 2. Results of single debris selection.
Table 2. Results of single debris selection.
Debris No.Start Point (s)End Point (s) Perigee before Ablation (km)Perigee in the End (km)Total Ablation Time (s)Total Impulses (Counts)
53.3669.1304547.13538.435.7644577
90.40064.0593646.80644.983.6587366
120.00251.0425560.64560.401.0400104
130.11100.2911551.62540.510.180118
160.05591.0536515.45509.860.9977100
200.03203.8829485.59482.393.8509385
220.01982.6900732.91730.932.6702267
232.44856.8648578.77333.934.4163442
243.721864.6199597.27225.0960.89816090
371.116124.2846685.82574.2123.16852317
6411.8880104.4701509.83199.8592.58219258
23410.0590167.2413474.85199.74157.182315,718
37024.1154140.9113520.92199.78116.795911,680
Table 3. Parameter settings for ant colony algorithm.
Table 3. Parameter settings for ant colony algorithm.
Parameters NameVariableValue
Number of antsmAnt20
Transfer Rule Parametersstr0.5
Importance of pheromonealpha1
Importance of heuristic functionbeta3
Importance of waiting timegama1
Importance of window width for orbit reductiondelta1
Pheromone volatilization factorrho0.85
Maximum number of iterationsiter_max20
Stop condition for deviation of objective functionMaxStallGen5
Table 4. Results of multiple debris selection with optimization.
Table 4. Results of multiple debris selection with optimization.
Sequence No.Debris No.Start Point (s)End Point (s)Perigee before Ablation (km)Perigee after Ablation (km)Total Ablation Time (s)Total Impulses (Counts)
76411.8880104.4701509.8295199.885092.58199259
16269108.8915162.5621595.7943389.822653.67025367
18234167.2413281.6858474.8490199.7381114.444511,445
1999286.8095468.8938524.2445450.2049183.084018,309
Table 5. Results of multiple debris selection without optimization.
Table 5. Results of multiple debris selection without optimization.
Sequence No.Debris No.Start Point (s)End Point (s)Perigee before Ablation (km)Perigee after Ablation (km)Total Ablation Time (s)Total Impulses (Counts)
13060.82381.1715554.2266553.25880.347634
83518.693850.4794556.0903455.080641.78604179
152464.6202162.1677597.2702227.399197.54789755
18234167.2413281.6858474.8490199.7381114.444511,445
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yang, W.; Fu, H.; Shao, Z.; Wu, Q.; Chen, C. Target Selection for a Space-Energy Driven Laser-Ablation Debris Removal System Based on Ant Colony Optimization. Sustainability 2023, 15, 10380. https://doi.org/10.3390/su151310380

AMA Style

Yang W, Fu H, Shao Z, Wu Q, Chen C. Target Selection for a Space-Energy Driven Laser-Ablation Debris Removal System Based on Ant Colony Optimization. Sustainability. 2023; 15(13):10380. https://doi.org/10.3390/su151310380

Chicago/Turabian Style

Yang, Wulin, Hongya Fu, Zhongxi Shao, Qiang Wu, and Chuan Chen. 2023. "Target Selection for a Space-Energy Driven Laser-Ablation Debris Removal System Based on Ant Colony Optimization" Sustainability 15, no. 13: 10380. https://doi.org/10.3390/su151310380

APA Style

Yang, W., Fu, H., Shao, Z., Wu, Q., & Chen, C. (2023). Target Selection for a Space-Energy Driven Laser-Ablation Debris Removal System Based on Ant Colony Optimization. Sustainability, 15(13), 10380. https://doi.org/10.3390/su151310380

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