Нова система координат підвищує надійність відстеження апаратів біля Місяця

У міру того як космічні агентства та приватні компанії готуються до тривалої присутності людини біля Місяця, на перший план виходить малопомітна, але критична задача: як надійно відстежувати космічні апарати у величезному та хаотичному просторі між Землею й Місяцем. Нова робота Мааніні Ґупти (Maaninee Gupta) та Кайла Дж. ДеМарса (Kyle J. DeMars) з Університету Пердью, опублікована в журналі Celestial Mechanics and Dynamical Astronomy і описана виданням Scienmag, пропонує перспективне рішення.

Дослідники показали, що спеціалізований набір орбітальних координат — Modified Generalized Equinoctial Orbital Elements (M-GEqOE), тобто модифіковані узагальнені еквідиктичні орбітальні елементи, може суттєво поліпшити прогнозування невизначеності положення та швидкості апаратів у цислунарному просторі — об’ємі, обмеженому орбітою Місяця навколо Землі.

Чому цислунарний простір такий складний для навігації

Традиційно система космічної обізнаності (Space Domain Awareness) зосереджувалася на навколоземному просторі, де супутники добре спостерігаються радарами й телескопами. У цислунарній зоні ситуація зовсім інша. Цей регіон у багато разів більший за околиці Землі, а рух об’єктів визначається конкуруючим тяжінням Землі та Місяця, до якого додається вплив Сонця.

Через таку складну динаміку невеликі помилки в оцінці стану апарата можуть перетворюватися на кардинально різні прогнозовані траєкторії. Це створює труднощі для точного знання того, де апарат був, де він є зараз і де перебуватиме в майбутньому.

До того ж більшість сучасних систем навігації спирається на припущення, що невизначеність стану апарата описується гаусовим (дзвоноподібним) розподілом. Поблизу Землі це працює задовільно. У цислунарному просторі через сильну нелінійність динаміки така «дзвоноподібна» хмара швидко перетворюється на викривлені, витягнуті й навіть багатогорбкові структури, які однією гаусівською кривою вже не опишеш.

Коли системи навігації та оцінки ризиків продовжують покладатися на гаусові припущення там, де вони вже не виконуються, вносяться помилки в каталоги орбіт, в оцінку зіткнень і виявлення маневрів. Це може погіршити безпеку місій, зокрема й таких, як майбутня місячна станція Gateway.

Перехід від декартових координат до орбітальних елементів

Ґупта та ДеМарс підійшли до проблеми, змінивши «мову», якою описується рух апарата. Замість звичних декартових координат (положення та швидкість у тривимірному просторі) вони застосували орбітальні елементи — величини, що задають розмір, форму й орієнтацію орбіти.

У своїй роботі вони розвинули ідею узагальнених еквідиктичних орбітальних елементів, запропонованих раніше Бауу (Baù) та колегами. Ці елементи формулюються через повну енергію орбітальної системи. Такий енергетичний підхід дає змогу одразу «вшити» збурювальні сили — як консервативні (наприклад, гравітація Місяця й Сонця), так і неконсервативні — безпосередньо в рівняння еволюції орбітальних елементів, а не додавати їх окремо через складні варіаційні рівняння.

Оригінальна схема мала суттєве обмеження: вона працювала лише для траєкторій із негативною повною енергією (тобто зв’язаних орбіт). Модифікація, запропонована авторами, долає цю проблему.

Замість великої півосі вони вводять узагальнений параметр півшироти (semi-latus rectum), а замість середньої довготи — істинну довготу. У такій формі набір M-GEqOE залишається коректним для значно ширшого спектра траєкторій, типових для цислунарного простору.

Шість елементів M-GEqOE включають: узагальнений параметр півшироти, два параметри узагальненого вектора ексцентриситету, два параметри, що задають орієнтацію орбітальної площини, та істинну довготу, яка швидко змінюється вздовж орбіти. Разом вони описують неосцилюючу еліптичну траєкторію, в самій геометрії якої вже враховано вплив збурень.

Перевірка в моделях із реалістичною динамікою

Щоб перевірити підхід, дослідники моделювали рух апаратів у високоточній гравітаційній моделі, що включала вплив Землі, Місяця й Сонця, використовуючи планетні ефемериди Лабораторії реактивного руху NASA.

Під час розрахунків вони зіткнулися з тонкою числовою проблемою: узагальнений кутовий момент, який залежить від ефективної потенціальної енергії, за фізично коректних умов має бути невід’ємним. Однак третє тіло (наприклад, Місяць) може зробити цей вираз від’ємним, навіть коли реальна фізика лишається коректною.

Вирішенням стало додавання постійного зсуву до потенціалу, що обчислюється один раз перед початком моделювання. Це змінює лише нульовий рівень потенціальної енергії, але не впливає на рівняння руху, оскільки вони залежать тільки від градієнта потенціалу.

Валідація охопила чотири характерні орбіти:

• Орбіта Near Rectilinear Halo 9:2 — базова орбіта місячної станції Gateway. Для неї M-GEqOE протягом приблизно шести з половиною днів добре узгоджувалися з класичною декартовою N-body моделлю, а помилки залишалися обмеженими, за винятком помірного зростання біля зближення з Місяцем (перилун).
• Еліптична місячна «заморожена» орбіта — сильно нахилена ексцентрична орбіта, придатна для тривалого спостереження південного полюса Місяця та майбутніх ретрансляційних супутників; її відстежували протягом 30 днів.
• Резонансні орбіти 4:1 та 3:1 із центром на Землі — з періодами близько 27 днів на оберт, які враховують комбінований вплив Землі, Місяця та Сонця. Одна з них подібна до орбіти, яку використовував апарат NASA Interstellar Boundary Explorer під час продовженої місії.

Як нові координати впливають на невизначеність

Найяскравіші результати стосувалися саме невизначеності. Для кожного сценарію вчені виконали Монте-Карло моделювання з 10 тисячами «частинок», які еволюціонували в однаковій динаміці, але описувалися або через M-GEqOE, або через декартові координати. Так можна було оцінити, як сам вибір координат впливає на форму розподілу похибок.

Гаусівськість розподілів перевіряли за допомогою критерію Генце–Цірклера (Henze–Zirkler) — статистичного тесту багатовимірної нормальності, що оцінює відповідність гаусовому розподілу прямо з вибіркових статистик. Це важливо, бо альтернативні показники, як-от дивергенція Кульбака–Лейблера, вимагають додаткових припущень про апроксимувальний розподіл, які можуть приховати справжні відхилення.

У всіх чотирьох випадках картина була схожою: невизначеність у координатах M-GEqOE довше залишалася близькою до гаусової, ніж у декартовому описі.

На орбіті Near Rectilinear Halo декартове представлення демонструвало різкий відхід від гаусівської форми саме під час проходження перилуну, де нелінійності найбільші, тоді як у M-GEqOE розподіл залишався гаусівським упродовж усього оберту.

Для резонансної орбіти 4:1 гаусівськість зберігалася до третього проходу перигею (приблизно 14-й день) у M-GEqOE, тоді як у декартових координатах вона «ламалася» вже на наступному проходженні перигею.

Навіть там, де обидва представлення зрештою втрачали гаусівські властивості, як на орбіті 3:1, відхилення в узагальнених елементах були помітно меншими й зникали швидше.

Візуалізація «хмар» невизначеності в просторі власних векторів коваріаційної матриці додатково підкреслила різницю: у декартових координатах з’являлися вигнуті, «хвостаті» структури, що явно виходили за рамки простої дзвоноподібної форми, тоді як ансамблі в M-GEqOE зберігали більш гладкий, наближено гаусівський вигляд.

Значення для майбутнього трафіку біля Місяця

Точне прогнозування невизначеності — фундамент для всіх задач оцінювання, відстеження й зіставлення даних про об’єкти на орбіті. Автори наголошують, що зменшення відхилень від гаусівської форми в M-GEqOE може спростити й негаусівські методи опису невизначеності, наприклад зменшити кількість компонент у фільтрах на основі гаусівських сумішей.

У міру зростання руху в цислунарному просторі — завдяки запланованим місячним станціям, комерційним місіям і науковим обсерваторіям — здатність точно передбачати, де перебувають апарати та наскільки впевнені ці прогнози, стане ключовою умовою безпечних та узгоджених операцій.

Робота Ґупти й ДеМарса демонструє, що для приборкання хаотичного середовища не завжди потрібні лише потужніші комп’ютери чи складніші моделі сил. Іноді вирішальним кроком стає розумніший вибір координат, у яких фізика задачі трьох тіл вбудована в математику з самого початку.