Adaptive Node Distribution via Probability Density Functions for Birkhoff Pseudospectral Trajectory Optimisation

Pseudospectral methods are a powerful framework for spacecraft trajectory optimisation that rely on fixed, problem-independent collocation nodes such as Legendre-Gauss-Lobatto (LGL) nodes. These nodes are located according to polynomial root distributions, which can result in insufficient resolution in regions of high dynamical activity, such as those encountered during collision avoidance manoeuvres. This paper presents an adaptive node distribution strategy for Birkhoff pseudospectral trajectory optimisation, where the collocation nodes are placed according to a probability density function derived from the thrust magnitude profile of the spacecraft. By concentrating nodes in regions of high dynamical variation, the proposed method simultaneously achieves feasibility preservation and improved convergence with fewer nodes in a collision avoidance case study when compared to the universal Birkhoff method, with Model Predictive Control serving as an independent validation comparator. The case study highlights the practical consequences of inter-node constraint violations and validates the adaptive node scheme as a robust and efficient alternative for fuel-optimal spacecraft manoeuvre design.