Keplera Орбіти: Три Закони, Одна Інверсія-Квадратна Сила

Від спостережень до законів, потім від законів до сили

Йоганн Кеплер десятиліттями працював над тим, щоб узгодити спостереження Тихо Брахе за Марсом та вивести три емпіричні правила руху планет на початку 1600-х років, за десятиліття до того, як хтось знав, чому вони правдиві. Це було потрібно Ісааку Ньютону, щоб показати, що всі три закони Кеплера не є незалежними фактами про сонячну систему, а є єдиним математичним наслідком однієї закону сили, застосованого до проблеми двох тіл.

Закон всесвітнього тяжіння Ньютона, F = -GMm/r², що спрямований по лінії з’єднання двох тіл, пояснив, чому всі три закони Кеплера є нерозривно пов'язаними.

Перший закон: еліпти, а не кола

Розв’язуючи рівняння Ньютона для об’єкта, що рухається по замкнутій орбіті під дією сили тяжіння зворотно-пропорційної квадрату відстані, отримуємо конічну перерізну гіперболу, у якій центральна маса знаходиться в одному фокусі, а не в центрі. Для закритого об’єкта (негативна загальна енергія) ця конічна гіпербола завжди є еліптичною; коло є лише особливим точним випадком з ексцентриситетом e = 0, який вимагає певного співвідношення між швидкістю запуску та радіусом. Зміна швидкості запуску навіть незначно призведе до розтягування орбіти в еліпс, де планета найближча до зірки знаходиться на перигелії, а найдальша – на апогеї, точно відповідно до форми, яку Кеплер вивів з непередбачувано нерухомого еліптичного шляху Марса.

Другий закон: рівні площі, наслідок збереження омосектного моменту

Закон Кеплера другого стверджує, що лінія від зірки до планети витікає рівними площами за рівні проміжки часу – планета рухається найшвидше біля перицентру та найповільніше біля афелу. Оскільки гравітаційна сила завжди напрямлена вздовж цієї ж лінії («центральна сила»), вона здійснює нульовий момент сил навколо зірки, тому збереження омосектного моменту L = m·r×v протягом орбіти. Швидкість, з якою витікає площа, дорівнює точно dA/dt = L/(2m), постійній – отже, рівні площі за рівні проміжки часу є нічим іншим, як збереженням омосектного моменту, вираженим у геометричній формі.

Третій закон: період в квадраті, напіввелике положення кубічне

T² = (4π² / G·M) · a³ T = орбітальний період a = напіввелике положення (середнє між найближчим і найвіддаленішим відстанню) M = маса центрального тіла (власна маса планети випадає, коли вона значно менша за M) Третій закон Кеплера випливає з балансування гравітаційного впливу проти вимог щодо центripetal сили для замкненої орбіти та інтегрування протягом одного періоду. Це закон, який найчастіше використовують астрономи: виміряйте період і напіввелике положення планети (або супутника, або зірки в системі подвійних зірок), і центральна маса M випливає безпосередньо - це буквально те, як визначали маси Сонця, планет з місяцями та незліченної кількості зір-господарів екзопланет, ніколи не ставлячи їх на ваги.

T² = (4π² / G·M) · a³

T  = orbital period
a  = semi-major axis (average of closest and farthest distance)
M  = mass of the central body (the planet's own mass drops out
     when it is much smaller than M)

Моделювання: чому важливий інтегратор

Рівняння руху самі по собі прості – прискорення дорівнює -GM r̂/r² – але стабільно інтегрувати їх чисельно протягом багатьох орбіт не автоматично. Наївний явний розклад Ейлера (просунути позицію за допомогою старої швидкості, а потім оновити швидкість) систематично вводить енергію в систему на кожному кроці; кругова орбіта, змодельована таким чином, помітно згортається зовні, і цей ефект погіршується, а не покращується, чим довше триває симуляція. Четвертий порядок Runge-Kutta (RK4) вирішує це для масштабів, які тут мають значення, відбираючи прискорення в чотири точки протягом кожного кроку – на початку, двічі в середині з пробними півкроками та в прогнозованій кінцевій точці – і об'єднуючи їх в оновлення, локальна помилка якого масштабується до п’ятого степеня часу кроку, що значно краще поводиться протягом десятків орбіт, які потрібно відобразити переконливо.

Енергія, що визначає форму орбіти

Один останній елемент об’єднує всю картину: загальна механічна енергія орбіти, E = ½mv² - GMm/r, визначає, яка конїчне секторіа визначиться. E < 0 дає обмежене еліпс (кеплюхове середовище); E = 0 дає параболічну траєкторію, яка ледь-ледь виривається в нескінченність з нульовою кінцевою швидкістю; E > 0 дає гіперболічний прохід, який виривається з надлишком швидкості. Коригування кутового моменту запуску у симуляторі, такому як цей, насправді полягає в підвищенні або зниженні E, і спостерігати, як орбіта змінюється від тісної еліптичної траєкторії через майже-круглу, широку еліптичну та врешті-решт розривається на параболічну та гіперболічну, – це найшвидший спосіб побачити всі три закони Кеплера – і їхні межі – в одному безперервному перемикачі.

Часті запитання

Спробуйте наживо

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