Общие рекомендации по моделированию ионной спирали с малой тягой?

Мне сказали, что для ионных спиралей с малой тягой дельта V будет разницей между скоростями стартовой и целевой орбит. Например, дельта V между LEO со скоростью 7,7 км/с и GEO со скоростью 3,1 км/с будет равна 4,6 км/с. Это правильно? Как это получается?

Для ионного двигателя, разгоняющего космический корабль до 1 миллиметра в секунду ^ 2, я получаю 11,6 дня, чтобы разогнаться до 1 км/с.

Чтобы рассчитать время, необходимое для перехода от LEO к GEO, я взял 11,6 * 4,6, чтобы получить 54 дня. Опять же, я не знаю, правильно ли это.

Вот попытка нарисовать спираль LEO-GEO с ускорением 1 мм/с^2:

введите описание изображения здесь

Это был мой ход мыслей: конечная орбита составляет около 24 часов, начальная орбита составляет 1,5 часа, поэтому средняя орбита составляет около 12,75 часов. !2,75 часа уходят в 54 дня примерно в 100 раз, поэтому я сделал логарифмическую спираль, которая 100 раз поворачивается от LEO к GEO. (ярко окрашенные области — пояса Ван Аллена)

Сильно подозреваю, что это неправильно. Мне кажется, что на НОО ионная спираль будет более туго накручена и постепенно расслабится по мере подъема космического корабля. Но на данный момент я понятия не имею, как смоделировать ионную спираль в электронной таблице Excel.

В этих вопросах я использую конкретный пример перехода от LEO к GEO, но я надеюсь на рекомендации для более общих сценариев.

Если вы хотите поиграть с подобными вещами, я настоятельно рекомендую вам закончить работу с Excel. Мне нравится Mathematica, у которой в настоящее время есть распродажа Pi Day для домашней и студенческой версий. (Я не имею никакого отношения к Wolfram и не имею к нему никакого отношения.)
@MarkAdler Крысы! Сейчас июнь (и 2016 тоже). Слишком поздно. Мне интересно, у них тоже есть дневная распродажа Тау? cf Tau объяснил Ви Харт , официальный сайт Tau Day .

Ответы (2)

Правило, которое у вас есть для общего Δ В спирали с малой тягой — это верхний предел, достигаемый, когда тяга становится равной нулю. Однако на это уходит бесконечное количество времени. Общая Δ В спирали с ненулевой тягой меньше, а время конечно. Но это хорошее эмпирическое правило для быстрых расчетов при попытке установить осуществимость.

Вывод эмпирического правила довольно прост. Посмотрите на бесконечно малый перенос Хомана. Вы обнаружите, что Δ В сумма двух бесконечно малых вспышек на начальной орбите и в апоапсисе переходной орбиты равна разнице орбитальных скоростей. Затем, если вы суммируете их для конечного повышения орбиты, вы получите разницу в Δ В начальной и конечной орбиты.

Чтобы узнать реальную сумму Δ В и для построения фактической траектории, которая не совершает бесконечное количество оборотов, прежде чем куда-нибудь попасть, лучше всего использовать численное интегрирование.

Вот пример спирали с круговой орбиты на побег ( С 3 знак равно 0 ):

спираль

Это нормализовано к начальной круговой орбите, где расстояния указаны в единицах начального радиуса орбиты, а ускорение постоянно при 10 3 гравитационного ускорения тела на начальном радиусе орбиты. Общая Δ В для побега составляет 0,856 начальной орбитальной скорости по сравнению с 1,0 для эмпирического правила. Общее время до побега составляет 136 начальных периодов обращения. Он обходит тело около 40 раз, прежде чем сбежать.

Первые несколько орбит достаточно близки, чтобы вы не могли их разобрать при показанном разрешении. Это становится еще хуже для меньших ускорений. 10 3 на самом деле довольно высок. Я выбрал его, чтобы вы могли лучше видеть спираль. Это время с низкой околоземной орбиты составляет около 8,5 суток. Типичный спиральный выход может быть больше похож на месяцы с ускорением в 10 4 начального ускорения свободного падения или меньше. Попытки нарисовать сплошной диск почти до конца, где вы видите спиральный выход.

Вот пример спирали с НОО (400 км) на ГСО с той же нормировкой и нормированным постоянным ускорением 10 4 . Это занимает около двух месяцев над 945 витками. В этом случае общая Δ В очень близко к эмпирическому правилу. Это упрощенно, так как конечный угол траектории полета здесь составляет около половины градуса. Так что есть время и Δ В осталось сделать круговую орбиту.

спираль к гео

Вы можете аппроксимировать этот график, продвигаясь по одной орбите за раз, используя период обращения, умноженный на ускорение, как Δ В и поднимая орбиту на соответствующую величину, соединяя каждую с линейно возрастающей спиралью.

Благодарю вас! Я вижу нормальный угол траектории полета к концу спирали от окружности до C3=0. Тяга горизонтальная на всем протяжении? Или выровнен с вектором скорости?
Всегда выравнивается с вектором скорости. Это наиболее эффективное использование тяги для увеличения удельной энергии. Финал γ составляет 31°.
@uhoh Извините, должно быть, в комментарии была опечатка. Я вернулся к исходному блокноту, в котором делал графики, да и финальные γ для случая ускорения 0,001 было 39,2°.
@MarkAdler, могу я использовать ваши иллюстрации траектории ионов в книжке-раскраске, которую я делаю? Я бы отдал тебе должное, конечно.
Конечно...........

Просто чтобы добавить к вашему первоначальному вопросу "это правильно?" - да, но только для спиралей между двумя круговыми орбитами.

Отвечая на ваш запрос о «рекомендациях для более общих сценариев»: если вы хотите найти ∆V для перехода с эллиптической орбиты, такой как GTO, и спиральной на GEO, вы можете следовать методу, такому как во втором абзаце Марка Адлера выше - с существенная корректировка. Для такого переноса лучше всего придерживаться короткой дуги вокруг апогея и принять, что часть дуги, удаленная от точного апогея, понесет потери. Оценка этой потери требует численного интегрирования, хотя, если вам действительно нужна приблизительная цифра, вы можете принять косинусную зависимость, то есть эффективность 1 в апогее, 0 на концах малой полуоси (1/4 орбиты).

Я думал о круговых орбитах. Но теперь, когда вы упомянули об этом, я хотел бы смоделировать орбиты от эллиптических до эллиптических. Итак, я проголосовал за ваш ответ. хотя я должен признаться, что я не уверен, что вы подразумеваете под косинусной зависимостью.
Со своей стороны, я признаюсь, что о «косинусной зависимости» я подумал двадцать лет назад, когда у меня не было времени изучать ее глубже. Я хотел использовать косинусоидальную форму, чтобы понять, какую длину дуги использовать для полной тяги. Будет методом проб и ошибок выбирать компромисс между эффективностью топлива и временем переброски. Вы также можете решить модулировать уровень или направление тяги (чтобы оставаться ортогональным радиус-вектору) вдали от апогея. Я не мог сразу сказать, что было бы лучше.
Вы могли бы получить приблизительное представление об эффектах добавления небольшого ∆V через уравнение vis-viva, т.е. проверить воздействие на большую полуось в местах, постепенно удаляющихся от апогея. v ^ 2 = ГМ (2 / г - 1 / а)
Одна последняя мысль, которая только что пришла ко мне: пример, который я описал, не даст вам точную круговую орбиту, вам нужно будет привести ее в порядок, как только вы пройдете большую часть пути.