Расчет бюджета дельты V от Земли до Меркурия

Контекст

При проверке базовых вычислений для вычисления необходимого бюджета дельты V, чтобы добраться до Меркурия, я испытываю некоторые трудности с интерпретацией числа 8650 м / с от пересечения Земли с Меркурием, как указано на карте метро Солнечной системы, как представлено в этом вопросе .Карта метро до планет с бюджетом дельта V.

Кроме того, я заметил, что у меня нет четкого понимания того, что подразумевается под «перехватом Земли». Следовательно, я хотел бы представить свои расчеты дельты V вместе с интерпретацией «(Земли) Пересечения» и спросить, не допустил ли я каких-либо ошибок в своих расчетах и/или предположениях.

Предположения

  • «Перехват Земли» означает, что космический корабль достиг космической скорости относительно Земли.
  • Сферический радиус Земли составляет примерно 6371000 м.
  • Земля считается однородным шаром.
  • Высота НОО составляет 250 км над поверхностью Земли.
  • мю Е а р т час "=" г М Е а р т час "=" 3,98 * 10 1 4
  • Космический аппарат может достигать НОО с дельтой V 9400 м/с, из которых 1600 м/с приходится на трение, что дает круговую скорость 7800 м/с вокруг НОО . Проверяя это с помощью уравнения Vis-Viva для р "=" а (круговая орбита).
    В С я р с л Е О "=" мю Е а р т час р л Е О "=" 3,98 * 10 1 4 6371000 + 250000 "=" 7753.18.. м с
    Это считается подтвержденным, поскольку на самом деле Лео может находиться на высоте чуть ниже или выше 250 км.
  • Затем вычисляется дельта V, необходимая для достижения скорости убегания Земли с НОО, сначала вычисляя скорость убегания на высоте НОО, используя :
    В л Е О е с с "=" 2 г М р л Е О "=" 2 мю Е а р т час р л Е О "=" 2 * 3,98 * 10 1 4 6371000 + 250000 "=" 10964,6 м с
    И вычитая из этого ранее вычисленную круговую скорость НОО:
    Δ В л Е О Е а р т час е с с а п е "=" В л Е О е с с В С я р с л Е О "=" 10964,6 7753,18 "=" 3211,46 м с
    Что, кажется, достаточно близко к бюджету, чтобы достичь точки перехвата Земли с НОО, как показано на карте метро Солнечной системы. Следовательно, это считается проверенным.
  • Этот V-побег интерпретируется как космический корабль, способный «свободно» перемещаться в любое положение на земной орбите вокруг Солнца. По сути, у него недостаточно энергии, чтобы двигаться выше или ниже по своей орбите вокруг Солнца, но он «свободен» от хватки Земли/сферы Хилла. (Обратите внимание, что я думаю, что это не совсем точно, потому что для фактического движения по орбите Земли вокруг Солнца потребовалась бы незначительная дельта V).
  • Поскольку предполагается, что «Перехват Земли» означает, что космический корабль просто находится где-то на той же высоте орбиты вокруг Солнца, что и Земля, предполагается, что для достижения Меркурия достаточно выйти на круговую орбиту, равную орбите Земли. большая полуось (0,387 а.е.) эллиптической орбиты Меркурия. Можно ударить Меркурия при ударе ради расчета.

Расчеты

На основе этих предположений выполняется фактический расчет от точки пересечения Земли до точки пересечения Меркурия. Сначала уравнения Vis-Viva переписываются для вычисления круговой скорости Земли и Меркурия:

В Е а р т час "=" мю ( 2 р Е а р т час С ты н 1 а Е а р т час С ты н ) В Е а р т час "=" мю ( 2 р Е а р т час С ты н 1 р Е а р т час С ты н ) В Е а р т час "=" мю ( 1 р Е а р т час С ты н ) В Е а р т час "=" мю С ты н р Е а р т час С ты н
Заполнение цифр:

  • мю С ты н "=" 1,33 10 20 м 3 с 2
  • р Е а р т час С ты н "=" 1,496 10 11 м Урожайность:
    В Е а р т час "=" мю С ты н р Е а р т час С ты н В Е а р т час "=" 1,33 10 20 1,496 10 11 "=" 29816.73075900643 м с
    Затем можно вычислить орбитальную скорость Меркурия вокруг Солнца, используя р М е р с ты р у С ты н "=" 0,387 1,496 10 11 :

В М е р с ты р у "=" мю С ты н р М е р с ты р у С ты н В М е р с ты р у "=" 1,33 10 20 0,387 1,496 10 11 "=" 47929.68129706198 м с

Отсюда требуемый Δ В можно вычислить как:

Δ В М е р с ты р у "=" В М е р с ты р у В Е а р т час "=" 47929.68129706198 29816.73075900643 "=" 18112.95053805555 м с
Однако это
18112.95053805555 м с
находится не по соседству
8650 м с
показан для перехвата Земли с перехватом Меркурия.

Код расчетов

Для полноты картины вот код Python, который выполнял фактические вычисления от перехвата Земли до перехвата Меркурия:

import math

# Initialize parameters:
mu_sun=1.33*10**20
r_earth_sun=1.496*10**11
r_mercury_sun_au=0.387
r_mercury_sun=r_mercury_sun_au*r_earth_sun

# Compute orbital velocities
v_earth=(mu_sun/r_earth_sun)**0.5
print(f'v_earth=\n{v_earth}')

v_mercury=(mu_sun/r_mercury_sun)**0.5
print(f'v_mercury=\n{v_mercury}')

dv_mercury=v_mercury-v_earth
print(f'dv_mercury=\n{dv_mercury}')

Вопрос

Что я сделал не так, рассчитывая расстояние от Земли до Меркурия?

Ответы (1)

Я вижу две основные вещи, которые сбивают с толку ваши расчеты:

  1. Вы не можете просто вычесть скорость одной планеты из скорости другой, чтобы получить затраты на межпланетный переход.

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

Таким образом, числа, которые вы должны попытаться получить,

  • разница в скорости между Землей и афелиевой скоростью такого эллипса.
  • разница в скорости между Меркурием и перигелийной скоростью такой орбиты.

Поскольку вы уже владеете уравнением Vis-Viva, вы должны быть в состоянии вычислить их оба.

  1. Ожидание побега с Земли (и Меркурия) для выполнения межпланетного перехода неэффективно.

То есть вы не просто складываете два числа из предыдущего шага. Вместо этого и импульс спасения, и переходы выполняются вместе на высоте LEO (и LMO).

Ключевое уравнение:

в 2 "=" в е 2 + в 2

Где в - скорость на НОО после выполнения импульса, в е скорость убегания и в - скорость относительно Земли после побега.


Рабочий пример земной стороны вещей. Вы должны были бы сделать то же самое на Меркурии, чтобы сократить пролёт до перехвата.

в а п час е л я о н "=" мю С ты н ( 2 а Е а р т час 2 а Е а р т час а М е р с ты р у ) 22250 м / с

Это относительная скорость к Земле 7533 м / с

Затем мы можем использовать:

в 2 "=" в е 2 + в 2

Введите скорость убегания LEO:

в 2 "=" ( 10964,6 м / с ) 2 + в 2

Введите желаемую относительную скорость после побега:

в 2 "=" ( 10964,6 м / с ) 2 + ( 7533 м / с ) 2

А затем решить:

в "=" 13303 м / с

Вычитая круговую скорость, мы имеем Δ в 1 "=" 5550 м / с

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

Является ли следующим точным пониманием общего значения вашего отзыва: я должен вычислить дельту V, используя перенос Хомана на Меркурий, где скорость в Афелии равна скорости вращения Меркурия вокруг Солнца, вместо того, чтобы вычислять разницу в круговых орбитах ?
@atYes, вам придется использовать правильный перевод Хохмана, хотя я старался быть осторожным, говоря, какие части перевода должны быть включены. (и «скорость в Афелии равна скорости вращения Меркурия вокруг Солнца» не имеет смысла).