Бюджеты Delta-v и модифицированная гравитация

Я строю мягкий научно-фантастический мир, в котором гравитация работает немного по-другому: ускорение, вызванное гравитацией, падает со скоростью, намного превышающей наш обратный квадрат; и ускорение из-за гравитации «переворачивается» и колеблется через равные промежутки времени, становясь скорее отталкивающим, чем привлекательным, и наоборот (похоже, я думаю, на MOND , но в меньших планетарных масштабах).

Ускорение г опыт объекта незначительной массы можно смоделировать с помощью:

г С е λ д ты потому что ( д ты + ф )

где

  • С - ускорение свободного падения на поверхности планеты 1 м / с 2 ),
  • λ - безразмерная постоянная затухания (описывает постепенное падение силы),
  • д - расстояние от поверхности в метрах,
  • ты - скалярная константа 2 (преобразует метры в радианы угловой меры),
  • ф - фазовый угол, почти всегда нулевой и им можно пренебречь.

1 Если вы будете следовать математике, гравитационные силы под поверхностью колеблются с экспоненциально большей величиной, что, очевидно, разорвет планету на куски. Итак, формула хороша только для гравитации над поверхностью. Ниже применяется другой (и неважный) набор правил.

2 ты "=" 0,5 π / д 0 , где д 0 расстояние между двумя точками невесомости или 1 / 2 длина волны косинуса.

е λ д ты и потому что ( д ты ) безразмерны, оставляя результат в единицах м / с / с . Умножение результата на массу объекта дает силу, приложенную к объекту, которая может быть положительной или отрицательной.

Поскольку это построение мира, многие из этих параметров переключаются и не возникают (похоже) из фундаментальных взаимодействий. Это нормально, я думаю. Мне не нужно полностью конкретизировать и обосновывать, почему эти константы именно такие, какими я их выбрал (или даже почему гравитация отличается), я могу просто спрятать их в черный угол и работать с их выводом.


У этой формы гравитации есть некоторые странные свойства. Вот график:

гравитационный график2

λ "=" 0,5 , С "=" 1

Когда объект поднимается над поверхностью, ускорение свободного падения падает до тех пор, пока на некотором расстоянии π / ты , ускорение становится равным нулю и инвертируется, превращаясь в отталкивающее. Векторы силы направлены в сторону от поверхности планеты, и объект ускоряется, пока не пересечет повторяющуюся точку невесомости и не начнет замедляться.
Однако, поскольку расстояние, на котором происходит ускорение/замедление, постоянно, а величина ускорения/замедления всегда уменьшается, недвигательный объект всегда будет ускоряться от планеты — навсегда — после того, как он пересечет начальную точку невесомости. над поверхностью планеты. Чтобы даже приблизиться к планете, объект должен иметь некоторую минимальную «скорость входа».


Как выглядят орбиты?

гравитационный график 1

Гравитация колеблется между положительным и отрицательным, притягательным и отталкивающим. (На графике притягивающие области выделены сплошным цветом.) Кажется, что устойчивые орбиты бывают двух разновидностей: орбиты, которые остаются исключительно внутри притягивающих областей (в основном круговые), и орбиты, которые колеблются между притягивающими и отталкивающими областями (в основном синусоидальные). Высокоэллиптические орбиты невозможны; если у объекта достаточно орбитальной энергии, чтобы пересечь более одной точки невесомости, он будет выброшен из системы. Траектории, которые почти пересекают вторую точку невесомости, являются «критическими».

критическая траектория 1 критическая траектория 2

(Привлекательные области выделены черным цветом, отталкивающие — темно-серым цветом. В центре — траектории относительно невращающегося тела.)
Здесь мы видим объект (второе изображение) на критически колеблющейся траектории, падающий на планету, но не обладающий достаточной кинетической энергией для преодоления область отталкивания и ускоряется обратно, хотя и недостаточно сильно, чтобы войти во вторую область отталкивания, и так до бесконечности. Тангенциальная скорость кажется почти независимой от радиальной, а это означает, что два объекта могут иметь общую орбитальную область, но при этом иметь совершенно разные скорости по сравнению друг с другом. Отсюда более толстая полоса прослеживаемых орбит (представляющих более быстрое движение) в пределах большей, более медленной и разреженной траектории.

Стабильные орбиты существуют только внутри областей притяжения, где векторы силы указывают на планету. Круговые траектории в этих регионах могут иметь совершенно разные орбитальные скорости из-за большого гравитационного градиента. В пике области, где вектор силы самый сильный, орбитальная скорость самая высокая, тогда как по обе стороны от пика скорости уменьшаются и приближаются к нулю. Геосинхронные траектории не ограничиваются пространством областей невесомости, но также могут колебаться между областями притяжения и отталкивания. На поверхности планеты это выглядело бы как приближение и удаление спутника (в одном и том же положении или по траектории в виде восьмерки) в течение дня.

несколько геосинхронных 1 несколько геосинхронных 2

Вблизи более высоких орбит проходят несколько синхронных траекторий, почти «подпрыгивающих на месте».


С чем у меня проблемы

Я думаю, что моя проблема в основном математическая (я знаю столько математики, сколько должен программист среднего уровня, но недостаточно, чтобы проанализировать и вывести уравнения, которые мне нужны, чтобы ответить на мой собственный вопрос, я думаю).
Я бы хотел, чтобы мои путешественники исследовали такие объекты, и хотя свойства этой гравитационной системы магические, сами путешественники скованы уравнением ракеты Циолковского .
Что мне нужно, так это метод вычисления дельта-v-бюджетов с учетом модифицированной гравитации.

Скажем, мои путешественники хотят взлететь с планеты на некоторую высоту над поверхностью (предположим, что там нет атмосферы). Им нужен метод, чтобы вычислить, сколько дельта-v требуется, чтобы добраться туда (что, вероятно, включает в себя замедление «там», поскольку отталкивающие области будут работать, чтобы ускорить их выход из системы), а также сколько дельты -v им нужно вернуться . Помните, что отталкивающие области помогут им выбраться с планеты (и, возможно, даже немного скруглить свою орбиту), но они будут сражаться с теми же самыми областями на обратном пути.
Возможно, проблема, выраженная в простейшем виде, будет состоять в том, чтобы выяснить, какая дельта-v требуется для радиального перемещения из одной области — отталкивающей или притягательной — в любую другую.
Ради простоты нас беспокоят только копланарные изменения траектории.

Я думаю, вам придется переформулировать уравнения vis-viva (которые представляют собой компромисс между потенциальной и кинетической энергией). В твоем мире все еще есть половина метра в квадрате?
@JamesMcLellan Да, я так думаю. Мое намерение состояло в том, чтобы изменить только то, как ведет себя гравитация. (На самом деле это неверно, но я не думаю, что это различие было слишком важным, чтобы включать его в вопрос. В моем вымышленном мире все еще существует обычная гравитация. Эта модифицированная гравитация либо является истинной эйнштейновской гравитацией, ведущей себя по-другому (каким-то образом ) , или дополнительное, накладывающееся поле, которое действует подобно гравитации, только в этом измененном виде.)
Вы предполагаете, что из-за нарушения закона сохранения механической энергии?
Просто хотел подтвердить
похоже на электронные волны ;D
Какое приложение вы используете для построения графика?
Просто чтобы убедиться, что мы имеем дело исключительно с проблемой 1 тела, верно? Нет центральной звезды или других планет, о которых нужно заботиться?
@Алиса, верно. Никаких других возмущающих систем.

Ответы (2)

Я думаю, вам нужно переформулировать уравнения vis-viva. Давайте попробуем это и посмотрим, что произойдет:

Круговая орбита:

  • г "=" г М р 2
  • в 2 "=" г М р в 2 "=" г р
  • г С е λ д ты потому что ( д ты + ф )

С д является р р , u — частотный член, С является г М р 2 и ф можно игнорировать, просто заменив:

  • г г М р 2 е λ ( р р ) ты потому что ( ( р р ) ты + 0 )
  • в 2 "=" [ г М р 2 е λ ( р р ) ты потому что ( ( р р ) ты + 0 ) ] р

Давайте попробуем! Используя Землю в качестве эталона (M = 5,9 × 10 24 , R=6378 км), попробуем с синусоидальными эффектами далеко ты "=" π р , и r на поверхности (6378 км ~ 6 378 000 м); Я получаю 9,67 м/с/с для g, это то же самое значение, которое я получаю, просто вычисляя г М р 2 . Я получаю орбитальную скорость 7855 м/с; что соответствует г М р

Если я перейду к синусоидальному эффекту при r = 1,5 R, g ~ 0 и v ~ 0, чего я и ожидал.

Итак, с большей уверенностью, давайте попробуем применить это к...

Эллиптическая орбита:

  • в 2 "=" г М ( 2 р 1 а ) в 2 "=" ( 2 г М р г М а ) в 2 "=" ( 2 г р г р 2 а )
  • г г М р 2 е λ ( р р ) ты потому что ( ( р р ) ты )
  • в 2 "=" ( 2 г М р 2 е λ ( р р ) ты потому что ( ( р р ) ты ) р г М р 2 е λ ( р р ) ты потому что ( ( р р ) ты ) р 2 а )

Тестирование снова, все еще с использованием Земли, и попытка эллиптической орбиты, где большая ось равна R (круговая орбита): 7855 м / с ... как раз то, что я ожидаю от круговой орбиты.

Как насчет синусоидального эффекта внутри при r = 1,5 R? v ~ 0 м/с; чего я ожидал.

Трансферы:

Разгон с более низкой орбиты р 1 к более высокому р 2

  • Δ в 1 "=" мю р 1 ( 2 р 2 р 1 + р 2 1 ) Δ в 1 "=" г р 1 ( 2 р 2 р 1 + р 2 1 )
  • Δ в 1 "=" ( г М р 2 е λ ( р р ) ты потому что ( ( р р ) ты ) ) р 1 ( 2 р 2 р 1 + р 2 1 )

И круговая новая орбита

  • Δ в 2 "=" мю р 2 ( 1 2 р 1 р 1 + р 2 ) Δ в 2 "=" г р 2 ( 1 2 р 1 р 1 + р 2 )
    Δ в 2 "=" ( г М р 2 е λ ( р р ) ты потому что ( ( р р ) ты ) ) р 2 ( 1 2 р 1 р 1 + р 2 )

Скорость убегания:

Не могу просто заменить здесь новый g.

  • 1 2 м в е 2 "=" г дельта р

Заимствование BMFs отличный ответ. г "=" С е λ д ты ( грех ( д ты + ф ) λ потому что ( д ты + ф ) ) ( λ 2 + 1 ) ты

Подставляя d и S и опуская фи

  • г дельта р "=" г М р 2 е λ ( р р ) ты ( грех ( ( р р ) ты ) λ потому что ( ( р р ) ты ) ) ( λ 2 + 1 ) ты

Затем, возвращаясь к уравнению скорости убегания

  • 1 2 м в е 2 "=" г М р 2 е λ ( р р ) ты ( грех ( ( р р ) ты ) λ потому что ( ( р р ) ты ) ) ( λ 2 + 1 ) ты

Отбрасывая m, и решая для в е

  • в е 2 "=" 2 [ г М р 2 е λ ( р р ) ты ( грех ( ( р р ) ты ) λ потому что ( ( р р ) ты ) ) ( λ 2 + 1 ) ты ]
Вау, это выглядит многообещающе. У меня будет время проанализировать это позже, но первые два графика г и в 2 выглядеть правильно
В уравнениях переноса Хомана, как вы пришли к замене мю р н с г р н ? (Также, с каким радиусом 1 или 2 г рассчитано?)
Я не расширил "g"= ... . Я полагал, что это было сделано дважды в круглом и эллиптическом примерах. Я чувствовал, что это было грязно. Я добавлю это.
Я думаю, что меня смущает то, что g принимает параметр радиуса немного r, а с уравнениями переноса у нас есть 2 маленьких r
@BMF, могу я одолжить твой интеграл?
Да, конечно, я написал это на случай, если это может помочь с ответом :)
В уравнениях переноса Хомана г можно сказать, что это функция р , и р 1 и р 2 и более низкие и более высокие орбиты, делая Δ в мгновенная дельта-v, необходимая для достижения р при изменении траектории от р 1 к р 2 ?

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

г С е λ д ты потому что ( д ты + ф )

по размеру д , то есть расстояние в метрах. Я не мог понять смысл этого результата (я подумал, что это может быть полезно для определения скорости), поэтому я задал этот вопрос, но теперь я думаю, что немного лучше понимаю, к чему мог привести этот результат. был.

г С е λ д ты потому что ( д ты + ф ) д д "=" С е λ д ты ( грех ( д ты + ф ) λ потому что ( д ты + ф ) ) ( λ 2 + 1 ) ты + в я

в я является «начальной скоростью». При установке минимума графика по положительной оси x весь график скорости лежит над осью x, что, как я полагаю, означает, что минимальное значение — это «скорость входа», необходимая для объекта, начиная с бесконечности, чтобы достичь д "=" 0 .

ускорение над скоростью

Уравнение гравитационного ускорения выделено красным (значения y увеличены в 20 000 раз для видимости), интегрирование этого уравнения выделено синим.

График интегральной формулы имеет минимальное значение вдоль положительной оси x, которое я интерпретировал как минимальную скорость, необходимую для приближения к планете (при условии, что интегральная формула представляет мгновенную скорость). Это минимальное значение существует на графике в точке д "=" д 0 , где д 0 это расстояние в метрах между двумя точками невесомости (я думаю, они называются точками перегиба).
Таким образом, оценивая интегральную формулу для д "=" д 0 и сдвинув график вверх по оси у на эту величину, мы должны найти график скорости падения объекта из бесконечности на поверхность планеты, упомянутую выше «скорость входа».

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

Какой программой вы там пользуетесь?
@ДжеймсМаклеллан Десмос . В частности, инструмент для построения графиков.
г можно найти из гравитационного потенциала В к г "=" В , так что вы нашли (с точностью до знака минус) гравитационный потенциал.
@ HDE226868 Да, я думаю, ты прав. п "=" г д , который очень похож на п Е г "=" м г час Вам нужно как-то экстраполировать уравнение для ускорения во времени (вместо высоты над землей) и интегрировать его. Кажется, это не простой способ сделать это.