Movimiento orbital · FísicaTercera ley de Kepler y rapidez orbital
Introducción
El movimiento orbital es la trayectoria cerrada y repetitiva que un cuerpo pequeño traza alrededor de uno mucho más grande cuando la gravedad aporta todo el giro lateral que el cuerpo pequeño necesita. La simulación fija una estrella en el centro del lienzo y lanza un planeta desde un radio elegido con la rapidez tangencial exacta que requiere una órbita circular, así que la trayectoria que se dibuja sola es el caso más simple que predice la ley de gravitación de Newton.
El tema sostiene cada satélite, cada sonda espacial y cada predicción sobre dónde estará Júpiter o Marte en el cielo el próximo mes. Los planificadores de misiones que eligen una ventana de lanzamiento para un rover de Marte usan la misma fórmula del periodo de Kepler que la simulación muestra en sus indicadores, y un satélite de comunicaciones geoestacionario se ubica exactamente en el radio donde el periodo orbital coincide con un día terrestre.
Cuesta no asumir que un planeta más pesado debe orbitar más rápido porque la gravedad lo atrae con más fuerza. La simulación revela que la masa propia del planeta no aparece en absoluto en la rapidez de órbita circular: con Radio de la Órbita = 100 m y Masa de la Estrella = 50, la lectura de Rapidez se asienta en 22,36 m/s sin importar cuánto pese el planeta, porque la atracción gravitatoria y la fuerza centrípeta requerida escalan juntas con esa masa.
La física explicada
Para que un planeta mantenga una trayectoria circular alrededor de la estrella central, la atracción gravitatoria hacia adentro G·M·m/r² debe igualar exactamente la demanda centrípeta m·v²/r. La masa orbitante m se cancela en ambos lados, dejando v² = G·M/r y por lo tanto v = sqrt(G·M/r). La simulación trabaja en unidades escaladas donde el valor s del deslizador Masa de la Estrella representa G·M = s · 10³ m³/s². En los valores por defecto de Radio de la Órbita = 100 m y Masa de la Estrella = 50, eso da G·M = 50000 m³/s² y una rapidez circular de sqrt(50000/100) = sqrt(500) ≈ 22,36 m/s, que es exactamente lo que la lectura de Rapidez reporta a lo largo de la corrida.
Cerrar la órbita lleva tiempo. Una revolución cubre una circunferencia de 2π·r, así que el periodo es T = 2π·r/v = 2π·sqrt(r³/(G·M)), la tercera ley de Kepler en su forma newtoniana. Para los valores por defecto, T = 2π·sqrt(100³/50000) = 2π·sqrt(20) ≈ 28,10 s. La simulación corre hasta 85 s, así que la lectura de Ángulo recorre desde 0° hasta 360° unas tres veces antes de que el límite de tiempo detenga el bucle, y la lectura de Tiempo en cada regreso a 0° coincide con el periodo predicho dentro del ruido de integración.
El radio y la masa central tiran del periodo en direcciones opuestas. Duplicar el radio de 100 m a 200 m con Masa de la Estrella = 50 multiplica T por 2^1,5 ≈ 2,83, así que el periodo crece de 28,10 s a unos 79,48 s y solo cabe una órbita completa dentro del límite de 85 s. Duplicar la Masa de la Estrella de 50 a 100 con Radio de la Órbita = 100 m eleva G·M a 100000 m³/s², encogiendo T a 2π·sqrt(100³/100000) ≈ 19,87 s, el mismo círculo, recorrido más rápido porque la atracción central es más fuerte.
La gravedad también fija un techo al movimiento ligado. La rapidez de escape al radio r es vesc = sqrt(2·G·M/r), que es exactamente sqrt(2) ≈ 1,414 veces la rapidez circular al mismo radio. Con los valores por defecto eso es sqrt(1000) ≈ 31,62 m/s, así que los 22,36 m/s que muestra la lectura de Rapidez dejan al planeta cómodamente ligado, con cerca del 71 % de la rapidez que necesitaría para liberarse.
Ecuaciones clave
Esta es la fuerza hacia adentro que la estrella central ejerce sobre el planeta a separación r. Con G·M = 50000 m³/s² y r = 100 m, la fuerza por unidad de masa del planeta es 50000 / 100² = 5 m/s², la aceleración centrípeta que la simulación sostiene durante la corrida completa de 85 s, y el origen de la flecha roja hacia adentro dibujada en el planeta en cada cuadro.
Igualando la atracción gravitatoria a la demanda centrípeta m·v²/r y cancelando la masa del planeta resulta esta expresión. En el Radio de la Órbita = 100 m y Masa de la Estrella = 50 por defecto (es decir, G·M = 50000 m³/s²), v = sqrt(50000 / 100) = sqrt(500) ≈ 22,36 m/s. La lectura de Rapidez se mantiene dentro de unas centésimas de ese valor en cada revolución, ya que el integrador simpléctico de Euler conserva la energía en órbitas keplerianas cerradas.
Una revolución completa traza una circunferencia 2π·r a rapidez constante v, así que el periodo es T = 2π·r / v, que se simplifica a esta forma con el radio al cubo. Para los valores por defecto, T = 2π · sqrt(100³ / 50000) = 2π · sqrt(20) ≈ 28,10 s. La lectura de Ángulo regresa a su 0° inicial con esa cadencia aproximada, y la lectura de Tiempo en ese momento confirma la predicción dentro del presupuesto de deriva del integrador.
Esta es la rapidez de lanzamiento al radio r por encima de la cual la energía mecánica total se vuelve positiva y la órbita ya no se cierra. En los valores por defecto, vesc = sqrt(2 · 50000 / 100) = sqrt(1000) ≈ 31,62 m/s, exactamente sqrt(2) ≈ 1,414 veces la rapidez circular de 22,36 m/s mostrada por la lectura de Rapidez. Cualquier lanzamiento por encima de vesc escaparía en una trayectoria hiperbólica en lugar de cerrarse en una órbita.
Variables clave
| Símbolo | Nombre | Unidad | Significado |
|---|---|---|---|
| G | Constante gravitatoria | m³/(kg·s²) | Constante universal 6,674·10⁻¹¹ en SI |
| M | Masa central | kg | Masa de la estrella en el centro del lienzo |
| G·M | Parámetro gravitatorio | m³/s² | Valor del deslizador Masa de la Estrella × 10³ |
| r | Radio orbital | m | Distancia centro a centro entre planeta y estrella |
| v | Rapidez orbital | m/s | Rapidez tangencial mostrada por la lectura de Rapidez |
| T | Periodo orbital | s | Tiempo de una revolución de regreso a Ángulo = 0° |
| vesc | Rapidez de escape | m/s | Umbral de trayectoria no ligada al radio r |
Ejemplos del mundo real
¿Cómo se elige el radio de un satélite geoestacionario?
Un satélite geoestacionario debe completar exactamente una órbita por día sidéreo, 86164 s, para que su posición angular sobre el ecuador nunca cambie respecto del suelo. Despejando r en la tercera ley de Kepler con ese periodo y el G·M de la Tierra ≈ 3,986·10¹⁴ m³/s² se obtiene r ≈ 42164 km desde el centro de la Tierra, o cerca de 35.786 km sobre el ecuador. Elige una altitud distinta y el satélite se desplaza hacia el este o el oeste por el cielo cada día.
La simulación captura el mismo acoplamiento de un solo sentido entre radio y periodo en unidades escaladas. Manteniendo Masa de la Estrella = 50 y subiendo Radio de la Órbita de 100 m a 159 m, el periodo predicho pasa de 28,10 s a 2π·sqrt(159³/50000) ≈ 56,20 s, exactamente el doble. El radio no es un parámetro libre una vez que el operador se compromete a una huella estacionaria, y la misma ley v = sqrt(G·M/r) que la lectura de Rapidez obedece a 22,36 m/s para la órbita por defecto de 100 m fuerza una rapidez única cercana a 3,07 km/s en cada satélite geoestacionario.
¿Por qué la Estación Espacial Internacional orbita tan rápido?
La ISS se sitúa aproximadamente a 400 km sobre la superficie de la Tierra, en r ≈ 6771 km desde el centro terrestre. La fórmula de órbita circular da v = sqrt(G·M/r) ≈ 7,66 km/s y un periodo cercano a 92 min, razón por la cual las tripulaciones ven cerca de dieciséis amaneceres y atardeceres al día. Los astronautas no están sin peso porque la gravedad se haya apagado (a esa altitud, la gravedad superficial sigue siendo aproximadamente el 89 % de su valor al nivel del mar) sino porque la estación y todo lo que hay dentro están en caída libre continua alrededor del planeta.
La simulación ilustra el mismo acoplamiento radio-rapidez a una escala más cómoda. Reducir el Radio de la Órbita de 200 m a 50 m con Masa de la Estrella = 50 eleva la Rapidez predicha de 15,81 m/s a sqrt(50000/50) ≈ 31,62 m/s (exactamente el doble de rápido a un cuarto del radio) y encoge el periodo de 79,48 s a unos 9,93 s. Cuanto más cercana es la órbita, más curvatura debe imponer la gravedad por segundo, y más rápido debe viajar el cuerpo en órbita para mantener el ritmo. Las órbitas bajas reales también se desgastan: la atmósfera tenue a 400 km ejerce suficiente arrastre para que la ISS necesite reimpulsos periódicos, mientras que la lectura de Distancia sin arrastre de la simulación se mantiene cerca de 100,00 m durante toda la corrida de 85 s.
¿Cómo envían los planificadores una sonda a Marte?
La transferencia más barata de la Tierra a Marte es una elipse de Hohmann cuyo perihelio se sitúa en la órbita de la Tierra y cuyo afelio apenas toca la órbita de Marte. El semieje mayor es el promedio de los dos radios (unas 1,262 UA) y la tercera ley de Kepler fija el tiempo de viaje de media órbita en aproximadamente 8,5 meses. Las ventanas de lanzamiento se abren cada 26 meses, cuando la Tierra y Marte se alinean para que la sonda llegue al mismo lugar donde Marte estará.
El mismo escalamiento periodo-radio aparece en la simulación. Con Masa de la Estrella = 50 fija, subir el Radio de la Órbita de 100 m a 200 m hace crecer T de 28,10 s a unos 79,48 s (un factor de 2^1,5) exactamente la pendiente de 1,5 que los astrónomos planetarios ven cuando grafican log T contra log r para los planetas del Sistema Solar. La lectura de Rapidez baja correspondientemente de 22,36 m/s a 15,81 m/s, la misma ley v ∝ 1/sqrt(r) que rige cada cálculo de órbita de transferencia, y la predicción de rapidez de escape de la simulación de vesc ≈ 31,62 m/s en los valores por defecto es el mismo umbral, escalado, que una sonda interplanetaria debe cruzar al dejar el pozo gravitatorio de la Tierra rumbo a Marte.
Lecturas adicionales
- Velocidad de escape: el umbral de rapidez sqrt(2)·vcirc en el cual una órbita ligada se convierte en una trayectoria hiperbólica y el cuerpo orbitante se aleja al infinito.
- Movimiento circular: el marco general de fuerza centrípeta que la gravedad satisface en el caso orbital, con ejemplos resueltos para cualquier ley de fuerza hacia adentro.
- Péndulo simple: otro sistema periódico de pequeña amplitud cuyo periodo depende de la geometría y no de la masa del cuerpo que oscila, en paralelo con la cancelación de m en la deducción de la rapidez orbital.