Mecánica orbital: ¿se estrellará un satélite?

Para un satélite con una velocidad y una posición iniciales sobre un planeta, me gustaría saber si experimentará una trayectoria elíptica, circular, parabólica o hiperbólica, o si chocará contra el planeta. La posición inicial ( X , y ) en relación con el centro de Marte, y la velocidad ( v X , v y ) son conocidos.

Siento que la órbita será circular si el potencial efectivo tu ( r ) = L 2 / 2 metro R 2 GRAMO METRO metro / R está en un mínimo (con respecto a r), donde L=momento angular. diferenciando, d tu d r = GRAMO METRO metro r 2 L 2 metro r 3 = 0 . Y por lo tanto, si se somete a una órbita circular, el radio es L 2 GRAMO METRO metro , ¿es esto correcto?

Si se ignora la resistencia del aire, ¿cuánta información se necesita para determinar si un satélite chocará contra un planeta o no? ¿Existe una energía mínima por debajo de la cual siempre entrará en espiral y chocará?

Nunca "entrará en espiral y chocará" mientras "se ignora la resistencia del aire" a menos que invoque la relatividad general o los efectos de múltiples cuerpos. En la mecánica newtoniana todas las "órbitas" de dos cuerpos en un r 1 potenciales---abiertas o cerradas---son secciones cónicas.
@dmckee: Bueno... todavía podría chocar contra el planeta si el pequeño eje de la órbita elíptica es más pequeño que el radio del cuerpo alrededor del cual está orbitando.
@ACuriousMind: sí, pero eso no es "en espiral".

Respuestas (2)

La ecuación de vis viva es

v 2 = GRAMO METRO ( 2 r 1 a )

Esto nos dice que si conocemos la velocidad y la distancia radial de un satélite en cualquier punto de la órbita, podemos conocerla en cualquier otro punto. Podemos usar esto para responder a su pregunta.

Primero: el satélite escapará si el parámetro a (el semieje mayor de la órbita) diverge (así que 1 a 0 ), que sucede cuando:

v > 2 GRAMO METRO r

La trayectoria será parabólica cuando la velocidad sea exactamente igual a la velocidad de escape (excentricidad = 1):

v = 2 GRAMO METRO r

Y para cualquier velocidad menor que esta, la trayectoria será elíptica. Que la órbita se sumerja por debajo de la superficie del planeta (es decir, que se estrelle) depende del valor de a y la excentricidad mi . Lo sabemos a está relacionado con la energía específica ϵ de la órbita por

ϵ = 2 GRAMO METRO a = 1 2 v 2 GRAMO METRO r

Así, basta conocer la velocidad y el radio instantáneo para conocer el parámetro a :

a = 2 GRAMO METRO GRAMO METRO r 1 2 v 2 = 2 r 1 v 2 r 2 GRAMO METRO

Esto solo será válido para v < v mi s C a pag mi .

Ahora necesitamos usar esto para encontrar la distancia de acercamiento más cercano. Aquí necesitamos usar el momento angular. Tomando prestada la inspiración de esta respuesta , escribimos

| METRO | = | r × v | = GRAMO METRO a ( 1 mi 2 )

Podemos resolver esto para mi . Entonces usamos

r pag r a = 1 mi 1 + mi

con el que te combinas r pag + r a = a para resolver r pag . Y si eso es menor que el radio del planeta, se estrellará.

Dejaré los detalles de ese cálculo como ejercicio. Probablemente hay un atajo que olvidé.

estas usando mi y ϵ arriba. Para evitar confusiones, debe usar uno u otro, pero no ambos.
@ ja72: no son lo mismo. mi es la excentricidad, ϵ es la energía característica. Estos son los símbolos comúnmente utilizados. ¿Tiene alguna sugerencia de símbolos alternativos?
@ ja72 Sí, en el párrafo justo después de la ecuación para la velocidad de escape (comenzando con "Y para cualquier velocidad menor que esta..."). A menos que la "excentricidad mi " no es una definición suficiente para ti (también doy la relación entre excentricidad y perigeo/apogeo).
Otra forma de encontrar la excentricidad es tomar la magnitud del vector de excentricidad .

Para un enfoque más conceptual: si la energía mecánica total ( k + tu ) es positivo, la órbita será hiperbólica. (Recuerda eso tu será negativo.)

Si k + tu = 0 , la órbita será parabólica.

Si k + tu < 0 , la órbita será algún tipo de elipse.

Si k + tu = tu 2 , la órbita elíptica será circular (excentricidad, mi = 0).

Si la órbita es elíptica y no circular, necesitas encontrar el periapsis:

r pag = a ( 1 mi ) .
Use el enfoque que da Floris para ese cálculo.