Hamilton--Jacobi Reachability for Spacecraft Collision Avoidance

Résumé technique : Atteignabilité Hamilton–Jacobi pour l'évitement de collisions spatiales

Énoncé du problème
L'expansion rapide des constellations en orbite basse (LEO), telles que Starlink de SpaceX, remet en cause les architectures existantes de gestion du trafic spatial (STM). Les méthodes traditionnelles de désencombrement au sol et de planification stratégique des manœuvres sont insuffisantes pour l'échelle des futures constellations (potentiellement 42 000 satellites). L'article aborde le besoin de couches de sécurité automatisées et fondées mathématiquement pour gérer l'évitement de collisions entre un satellite contrôlé et un satellite non coopératif ou secondaire au sein d'une même orbite circulaire. Le défi central consiste à déterminer quand initier des manœuvres d'évitement pour garantir la sécurité malgré l'intention inconnue et les perturbations bornées du satellite non coopératif.

Méthodologie
Les auteurs proposent un cadre formel basé sur l'analyse d'atteignabilité Hamilton–Jacobi (HJ), formulée comme un jeu différentiel à somme nulle.

1. Dynamique et modélisation :

   • Le mouvement relatif est modélisé dans le repère Radial–Tangente–Normale (RTN) en utilisant la dynamique planaire de Hill–Clohessy–Wiltshire (HCW), fournissant une approximation localement linéaire du mouvement orbital.
   • L'état du système $X = [x, y, \dot{x}, \dot{y}]^\top$ représente les positions et vitesses tangentes et radiales.
   • Joueur 1 (Satellite contrôlé) : Applique des commandes bornées (accélérations $u$) pour éviter les collisions.
   • Joueur 2 (Satellite non coopératif) : Modélisé comme une perturbation adverse bornée $d$ à intention inconnue, représentant l'accélération relative au pire cas.

2. Architecture de système hybride :
   Le processus de prise de décision est modélisé comme un automate hybride comportant trois modes discrets :

   • $q_{nom}$ : Maintien de position nominal ou orbite.
   • $q_{eva}$ : Exécution de manœuvres actives d'évitement de collisions.
   • $q_{ret}$ : Phase de récupération de trajectoire.
     Les transitions entre ces modes sont régies par des « conditions de garde » basées sur l'état du système par rapport à l'Ensemble Atteignable en Retour (BRS) et aux régions de récupération sûre.

3. Formulation de l'atteignabilité :

   • Ensemble cible ($T$) : Défini comme l'ensemble des configurations relatives dangereuses violant les normes de séparation minimales de la Commission fédérale des communications (FCC) (une zone de collision physique).
   • Ensemble Atteignable en Retour (BRS) : Calculé via l'équation aux dérivées partielles de Hamilton–Jacobi–Isaacs (HJI). Le BRS caractérise tous les états initiaux à partir desquels une collision est inévitable sous les perturbations au pire cas, indépendamment de la commande optimale du Joueur 1.
   • Garantie de sécurité : Les états situés en dehors du BRS admettent des trajectoires prouvées sans collision. La frontière du BRS sert de surface de commutation critique ; si l'état du système s'approche de cette frontière, l'automate hybride déclenche une transition vers le mode d'évitement ($q_{eva}$).

4. Stratégies de contrôle :

   • Nominal/Récupération : Utilise des lois de commande proportionnelle-dérivée (PD) pour amener le satellite aux positions et vitesses relatives souhaitées.
   • Manœuvres d'évitement : Quatre modes d'évasion spécifiques sont définis basés sur la géométrie RTN :
     1. Poussée radiale positive (dérive vers l'arrière).
     2. Poussée radiale négative (accélère vers l'avant).
     3. Freinage purement tangentiel.
     4. Accélération purement tangentielle.

Contributions clés

• Cadre formel de sécurité : L'article intègre l'atteignabilité HJ avec une logique de contrôle de supervision hybride pour fournir des garanties de sécurité rigoureuses et au pire cas pour l'évitement de collisions spatiales.
• Modélisation adverse : Contrairement aux méthodes probabilistes de filtrage, cette approche traite le satellite non coopératif comme une perturbation bornée, assurant la sécurité même sous les manœuvres les plus agressives (dans les limites) de l'autre partie.
• Conception d'automate hybride : Les auteurs définissent une logique de commutation spécifique déterminant exactement quand abandonner les opérations nominales pour des manœuvres d'évitement et quand revenir au maintien de position, basée sur les ensembles atteignables calculés.
• Considérations d'évolutivité : Le cadre est conçu avec des architectures décentralisées à l'esprit, adapté au désencombrement au niveau individuel des véhicules dans de grandes constellations.

Résultats
Les auteurs ont validé le cadre en utilisant des simulations numériques avec une boîte à outils modifiée des méthodes de niveaux (Level Set Methods Toolbox).

• Paramètres de simulation : Les simulations supposaient une orbite LEO à 500 km, des limites de poussée symétriques de $\pm 0,1 \, \text{m/s}^2$ pour le satellite contrôlé, et des perturbations bornées de $\pm 0,05 \, \text{m/s}^2$ pour le satellite non coopératif.
• Vérification : L'étude a calculé avec succès des ensembles de capture 2D et 3D et des tubes atteignables en retour. Les résultats démontrent que l'automate hybride identifie correctement les frontières de commutation entre les modes nominal, d'évitement et de récupération.
• Validation de la sécurité : Les ensembles calculés confirment que les états en dehors du BRS permettent des trajectoires sans collision sous les lois de commande définies et les bornes de perturbation.

Signification et revendications
L'article revendique offrir une alternative mathématiquement rigoureuse aux méthodes d'évitement de collisions purement probabilistes. En utilisant l'atteignabilité HJ, le cadre fournit des marges de sécurité explicites au pire cas, ce qui est crucial pour la durabilité des futures constellations de satellites à grande échelle. Les auteurs soulignent que leur approche permet des « comportements prouvés sûrs » et des « décisions de supervision fondées mathématiquement ».

L'article reste modeste concernant sa portée actuelle, notant que :

• L'analyse est restreinte à la dynamique HCW planaire (2D), omettant le mouvement hors plan, les perturbations séculaires (comme $J_2$) et les effets orbitaux non linéaires.
• Les ensembles atteignables calculés sont des garanties formelles spécifiques au modèle linéarisé adopté.
• La formulation actuelle est conservative ; elle suppose que toutes les réalisations de perturbations admissibles sont possibles, ce qui peut surestimer les états dangereux.

Les travaux futurs identifiés par les auteurs incluent l'extension du cadre au mouvement relatif complet en 3D, l'intégration de perturbations orbitales plus réalistes, et le développement d'approximations évolutives (telles que la décomposition ou des fonctions de valeur apprises) pour répondre à la complexité computationnelle des calculs d'ensembles atteignables de haute dimension.