Formulation and Analysis for Integrated Spacecraft Routing and Trajectory Design Problem
Pith reviewed 2026-05-10 03:35 UTC · model grok-4.3
The pith
Two formulations for integrated spacecraft routing and trajectory optimization reveal distinct trade-offs in optimality, speed, and reliability.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The paper introduces an arc-based formulation that employs an iterative decoupling scheme alternating between mixed-integer linear programming and sequential nonlinear trajectory optimization, and a path-based formulation that defines variables on paths and uses column generation together with a labeling algorithm. A geosynchronous satellite servicing case study together with numerical experiments quantifies the trade-offs between the two in terms of solution optimality, computational time, and robustness against non-converging or trivial solutions.
What carries the argument
Arc-based formulation with iterative decoupling between mixed-integer linear programming and nonlinear trajectory optimization; path-based formulation with column generation and labeling algorithm for route identification.
If this is right
- Mission planners can select the arc-based method when solution quality is paramount and the path-based method when faster computation is needed.
- Partial en-route propellant replenishment can be incorporated into routing decisions without separating routing from trajectory design.
- Numerical benchmarks provide guidance on expected run times and solution gaps for similar servicing scenarios.
- Robustness metrics help identify problem sizes where one formulation is more likely to produce usable results.
Where Pith is reading between the lines
- Hybrid solvers that switch between the two formulations based on instance size could combine their strengths.
- The observed robustness issues point to a need for improved warm-starting techniques for the nonlinear subproblems.
- Similar integrated routing-trajectory models could be tested on low-Earth orbit constellations or asteroid rendezvous missions.
Load-bearing premise
Nonlinear trajectory optimization subproblems can be solved reliably and the iterative decoupling scheme converges to useful solutions without excessive computational cost or trivial outcomes.
What would settle it
A set of test instances where the nonlinear solver fails to converge for a large fraction of subproblems or the overall procedure repeatedly returns trivial routes would show the formulations are not yet practical for real missions.
read the original abstract
This paper studies the integrated spacecraft routing and trajectory optimization problem for satellite servicing missions involving partial en-route propellant replenishment. Unlike terrestrial routing problems, spacecraft operate in a dynamic environment, and we need to optimize the spacecraft routing over a network with nonlinear and time-dependent trajectory costs. In this paper, we tackle this problem using two different formulations. The first formulation, referred to as the arc-based formulation, defines variables based on arcs and employs an iterative decoupling scheme that alternates between mixed-integer linear programming and sequential nonlinear trajectory optimization. The second formulation, referred to as the path-based formulation, defines variables based on paths/routes and leverages column generation and a labeling algorithm to accelerate the identification of promising routes. Through a geosynchronous satellite servicing case study and numerical experiments, we quantify the computational trade-offs between these two formulations in terms of the solution optimality, computational time, and robustness against non-converging or trivial solutions.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies the integrated spacecraft routing and trajectory optimization problem for satellite servicing missions with partial en-route propellant replenishment. It presents an arc-based formulation that uses an iterative decoupling scheme alternating between mixed-integer linear programming and sequential nonlinear trajectory optimization, and a path-based formulation that employs column generation together with a labeling algorithm. The two approaches are compared on a geosynchronous satellite servicing case study, with numerical experiments quantifying trade-offs in solution optimality, computational time, and robustness against non-converging or trivial solutions.
Significance. If the reported numerical results hold, the work supplies concrete guidance on selecting between decomposition strategies for hybrid discrete-continuous optimization problems that arise in space logistics. The explicit quantification of robustness metrics directly addresses a practical concern for mission planners. The manuscript is strengthened by the provision of the two complete mathematical programs, the labeling algorithm, and the focused case-study evaluation.
minor comments (2)
- The abstract outlines the methods and evaluation approach but does not report any quantitative results (e.g., relative optimality gaps or run-time ratios) from the geosynchronous case study; adding one or two such figures would allow readers to gauge the claimed trade-offs immediately.
- Section describing the labeling algorithm: a compact pseudocode listing would improve clarity and reproducibility of the column-generation procedure.
Simulated Author's Rebuttal
We thank the referee for the positive assessment of our work on arc-based and path-based formulations for the integrated spacecraft routing and trajectory optimization problem, including the recognition of its practical guidance on decomposition strategies and robustness metrics. The recommendation for minor revision is noted.
Circularity Check
No significant circularity detected in formulations or claims
full rationale
The paper presents two distinct optimization formulations (arc-based with iterative MILP-NLP decoupling, path-based with column generation and labeling) for the integrated routing-trajectory problem. These are derived from standard mixed-integer and nonlinear programming techniques without any reduction of a claimed result to its own fitted inputs or self-referential definitions. No load-bearing self-citations, uniqueness theorems from prior author work, or ansatzes smuggled via citation appear in the text. The geosynchronous case study supplies direct numerical measurements of optimality, runtime, and robustness, keeping the central quantification externally grounded rather than tautological.
Axiom & Free-Parameter Ledger
Reference graph
Works this paper leans on
-
[1]
Space Logistics Modeling and Optimization: Review of the State of the Art,
Ho, K., “Space Logistics Modeling and Optimization: Review of the State of the Art,”Journal of Spacecraft and Rockets, Vol. 61, No. 5, 2024, pp. 1417–1427. https://doi.org/10.2514/1.A35982
-
[2]
Orbital Depot Location Optimization for Satellite Constellation Servicing with Low-Thrust Transfers,
Choi, E., and Ho, K., “Orbital Depot Location Optimization for Satellite Constellation Servicing with Low-Thrust Transfers,” Journal of Spacecraft and Rockets, 2025, pp. 1–10. https://doi.org/10.2514/1.A36386
-
[3]
In-space Servicing, Assembly, and Manufacturing (ISAM) State of Play - 2024 Edition,
Arney, D., Mulvaney, J., Williams, C., Hernandez, W. A. R., Friz, J., Stockdale, C., Nelson, J., and Vargas, R. R., “In-space Servicing, Assembly, and Manufacturing (ISAM) State of Play - 2024 Edition,” Tech. rep., NASA Langley Research Center, 2024
work page 2024
-
[4]
Framework for Modeling and Optimization of On-Orbit Servicing Operations Under Demand Uncertainties,
Sarton du Jonchay, T., Chen, H., Gunasekara, O., and Ho, K., “Framework for Modeling and Optimization of On-Orbit Servicing Operations Under Demand Uncertainties,”Journal of Spacecraft and Rockets, Vol. 58, No. 4, 2021, pp. 1157–1173. https://doi.org/10.2514/1.A34978
-
[5]
On-Orbit Servicing Optimization Framework with High- and Low-Thrust Propulsion Tradeoff,
Sarton du Jonchay, T., Chen, H., Isaji, M., Shimane, Y., and Ho, K., “On-Orbit Servicing Optimization Framework with High- and Low-Thrust Propulsion Tradeoff,”Journal of Spacecraft and Rockets, Vol. 59, No. 1, 2022, pp. 33–48. https://doi.org/10.2514/1.A35094. 28
-
[6]
Choi, E., and Ho, K., “Optimal Routing and Trajectory Planning for Asteroid Mining With Partial In-Situ Resource Utilization,” AIAA SCITECH 2026 Forum, 2026, pp. 1–14. https://doi.org/10.2514/6.2026-2787
-
[7]
Ahn, J., de Weck, O., Geng, Y., and Klabjan, D., “Column generation based heuristics for a generalized location routing problem with profits arising in space exploration,”European Journal of Operational Research, Vol. 223, No. 1, 2012, pp. 47–59. https://doi.org/https://doi.org/10.1016/j.ejor.2012.06.018
-
[8]
Optimal Planetary Surface Routing Considering Exploration Values and Synergies,
Choi, E., and Ahn, J., “Optimal Planetary Surface Routing Considering Exploration Values and Synergies,”Journal of Aerospace Information Systems, Vol. 19, No. 6, 2022, pp. 406–420. https://doi.org/10.2514/1.I011052
-
[9]
Optimal Planetary Surface Exploration with Heterogeneous Agents,
Choi, E., and Ahn, J., “Optimal Planetary Surface Exploration with Heterogeneous Agents,”Journal of Spacecraft and Rockets, Vol. 60, No. 4, 2023, pp. 1343–1354. https://doi.org/10.2514/1.A35632
-
[10]
Thevehicleroutingproblem: Stateoftheartclassificationandreview,
Braekers,K.,Ramaekers,K.,andVanNieuwenhuyse,I.,“Thevehicleroutingproblem: Stateoftheartclassificationandreview,” Computers & Industrial Engineering, Vol. 99, 2016, pp. 300–313. https://doi.org/https://doi.org/10.1016/j.cie.2015.12.007
-
[11]
A concise guide to existing and emerging vehicle routing problem variants,
Vidal, T., Laporte, G., and Matl, P., “A concise guide to existing and emerging vehicle routing problem variants,”European Journal of Operational Research, Vol. 286, No. 2, 2020, pp. 401–416. https://doi.org/https://doi.org/10.1016/j.ejor.2019.10.010
-
[12]
Archetti, C., Speranza, M. G., and Vigo, D.,Vehicle Routing Problems with Profits, MOS-SIAM Series on Optimization, Society for Industrial and Applied Mathematics, 2014, Chap. 10, pp. 273–297. https://doi.org/10.1137/1.9781611973594.ch10, 0
-
[13]
A Green Vehicle Routing Problem,
Erdoğan, S., and Miller-Hooks, E., “A Green Vehicle Routing Problem,”Transportation Research Part E: Logistics and Transportation Review, Vol. 48, No. 1, 2012, pp. 100–114. https://doi.org/https://doi.org/10.1016/j.tre.2011.08.001, select Papers from the 19th International Symposium on Transportation and Traffic Theory
-
[14]
Green vehicle routing problem: A state-of-the-art review,
Asghari, M., and Mirzapour Al-e-hashem, S. M. J., “Green vehicle routing problem: A state-of-the-art review,”International Journal of Production Economics, Vol. 231, 2021, p. 107899. https://doi.org/https://doi.org/10.1016/j.ijpe.2020.107899
-
[15]
Thegreenvehicleroutingproblem: Asystematicliteraturereview,
Moghdani,R.,Salimifard,K.,Demir,E.,andBenyettou,A.,“Thegreenvehicleroutingproblem: Asystematicliteraturereview,” Journal of Cleaner Production, Vol. 279, 2021, p. 123691. https://doi.org/https://doi.org/10.1016/j.jclepro.2020.123691
-
[16]
Optimization for Multitarget, Multispacecraft Impulsive Rendezvous Considering J2 Perturbation,
Chen, S., Jiang, F., Li, H., and Baoyin, H., “Optimization for Multitarget, Multispacecraft Impulsive Rendezvous Considering J2 Perturbation,”Journal of Guidance, Control, and Dynamics, Vol. 44, No. 10, 2021, pp. 1811–1822. https://doi.org/10.2514/ 1.G005602
work page 2021
-
[17]
Advances in Space Research61(5), 1273–1285 (2018) https://doi.org/ 10.1016/j.asr.2017.12.025
Bang, J., and Ahn, J., “Two-phase framework for near-optimal multi-target Lambert rendezvous,”Advances in Space Research, Vol. 61, No. 5, 2018, pp. 1273–1285. https://doi.org/10.1016/j.asr.2017.12.025
-
[18]
Multitarget Rendezvous for Active Debris Removal Using Multiple Spacecraft,
Bang, J., and Ahn, J., “Multitarget Rendezvous for Active Debris Removal Using Multiple Spacecraft,”Journal of Spacecraft and Rockets, Vol. 56, No. 4, 2019, pp. 1237–1247. https://doi.org/10.2514/1.A34344. 29
-
[19]
Optimal Multitarget Rendezvous Using Hybrid Propulsion System,
Lee, D., and Ahn, J., “Optimal Multitarget Rendezvous Using Hybrid Propulsion System,”Journal of Spacecraft and Rockets, Vol. 60, No. 2, 2023, pp. 689–698. https://doi.org/10.2514/1.A35540
-
[20]
Large-scale object selection and trajectory planning for multi-target space debris removal missions,
Barea, A., Urrutxua, H., and Cadarso, L., “Large-scale object selection and trajectory planning for multi-target space debris removal missions,”Acta Astronautica, Vol. 170, 2020, pp. 289–301. https://doi.org/https://doi.org/10.1016/j.actaastro.2020.01. 032
-
[21]
OptimizationofRefuelingStrategiesforElectricPropulsionSpaceMissions,
Bendimerad,R.,Savransky,D.,andPetro,E.M.,“OptimizationofRefuelingStrategiesforElectricPropulsionSpaceMissions,” Journal of Spacecraft and Rockets, 2025, pp. 1–13. https://doi.org/10.2514/1.A36441
-
[22]
https: //doi.org/10.1007/978-1-4614-9149-1
Simchi-Levi, D., Chen, X., and Bramel, J.,The Logic of Logistics: Theory, Algorithms, and Applications for Logistics Management, Springer Series in Operations Research and Financial Engineering, Springer, New York, NY, 2014. https: //doi.org/10.1007/978-1-4614-9149-1
-
[23]
A branch-and-cut-and-price approach for the pickup and delivery problem with shuttle routes,
Masson, R., Ropke, S., Lehuédé, F., and Péton, O., “A branch-and-cut-and-price approach for the pickup and delivery problem with shuttle routes,”European Journal of Operational Research, Vol. 236, No. 3, 2014, pp. 849–862. https://doi.org/https: //doi.org/10.1016/j.ejor.2013.08.042, URL https://www.sciencedirect.com/science/article/pii/S0377221713007261, ...
-
[24]
A column generation algorithm for vehicle scheduling and routing problems,
Faiz, T. I., Vogiatzis, C., and Noor-E-Alam, M., “A column generation algorithm for vehicle scheduling and routing problems,” Computers & Industrial Engineering, Vol. 130, 2019, pp. 222–236. https://doi.org/https://doi.org/10.1016/j.cie.2019.02.032
-
[25]
Izzo, D., Binns, W., dariomm098, Mereta, A., Sprague, C. I., dhennes, den Abbeele, B. V., Andre, C., Nowak, K., Guy, N., Murguía, A. I. B., Pablo, Chapoton, F., GiacomoAcciarini, v. Looz, M., dietmarwo, Heddes, M., Babenia, A., Fournier, B., Simon, J., Willitts, J., Polnik, M., Narayanaswamy, S., Badger, T. G., and Yarndley, J., “esa/pykep: Optimize,” , O...
-
[26]
Bertsimas, D., and Tsitsiklis, J.,Introduction to linear optimization, Athena Scientific, 1997
work page 1997
-
[27]
Orbital Facility Location Problem for Satellite Constellation Servicing Depots,
Shimane, Y., Gollins, N., and Ho, K., “Orbital Facility Location Problem for Satellite Constellation Servicing Depots,”Journal of Spacecraft and Rockets, Vol. 61, No. 3, 2024, pp. 808–825. https://doi.org/10.2514/1.A35691. [28]Gurobi Optimizer Reference Manual, Gurobi Optimization, LLC, 2025. 30
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.