Trayectorias hiperbólicas

Supongamos un asteroide de masa m que se aproxima a un planeta de masa M con velocidad v0 en el infinito siguiendo una dirección (recta de color rojo) que dista b del centro del planeta tal como se muestra en la figura. A medida que el asteroide se acerca al planeta, va cambiando su velocidad tanto en módulo como en dirección (en color azul), hasta que se aleja hacia el infinito siguiendo una dirección simétrica respecto del eje X, que es la asíntota de la otra rama de la hipérbola.

Para analizar el problema supondremos que la masa del asteroide es muy pequeña comparada con la del planeta, m<<M y no tendremos en cuenta la influencia del Sol ni la del resto del los planetas.

La ecuación de la trayectoria es

r= d 1+εcosθ ε= 1+ 2 L 2 E m 3 G 2 M 2 d= L 2 GM m 2

r es la distancia del centro del planeta al asteroide y θ es el ángulo que hace el radio vector que une el asteroide con el centro del planeta con el eje X, tal como se muestra en la figura.

El parámetro ε denominado excentricidad define el tipo de trayectoria. La trayectoria es una hipérbola si ε>1 es decir, si la energía del asteroide E>0.

La energía del asteroide en el infinito es solamente cinética

E= 1 2 m v 0 2

El momento angular es

L=mv0b

Los parámetros ε y d valen

 ε= 1+ b 2 v 0 4 G 2 M 2 d= b 2 v 0 2 GM

La distancia más cercana del asteroide al centro del planeta ocurre cuando θ=0

r p = d 1+ε = b 2 v 0 2 GM 1+ 1+ b 2 v 0 4 G 2 M 2

Calculamos esta distancia, aplicando la conservación de la energía y del momento angular en el punto P más cercano al centro del planeta.

L=m v 0 b=m v p r p E= 1 2 m v 0 2 = 1 2 m v p 2 G Mm r p

En este sistema de dos ecuaciones con dos incógnitas despejamos vp y rp

v p = GM b v 0 ( 1+ 1+ b 2 v 0 4 G 2 M 2 ) r p = b 2 v 0 2 GM 1+ 1+ b 2 v 0 4 G 2 M 2

Dependiendo de las magnitudes M, v0, b y el radio R del planeta el asteroide choca rp<R o pasa el planeta. Cuando rp=R el asteroide choca tangencialmente con el planeta.

Asíntotas

El ángulo para el cual r→∞, es

α=arccos( 1 ε )

La dirección de la velocidad del asteroide cambia a medida que se acerca al planeta y luego, cuando se aleja. El cambio en la dirección de la velocidad es φ=2α-180 tal como se aprecia en la figura.

tan(φ/2)=-cotα. Utilizando la relación trigonométrica

1+ tan 2 α= 1 cos 2 α tan φ 2 = GM b v 0 2

La asíntota de una de las ramas de la hipérbola hace un ángulo α con el eje X y dista b del origen.

Donde x es la intersección de la asíntota con el eje X

x= b sin(180α) = bε ε 2 1

Ejemplo

Para el planeta Júpiter cuya masa es M=1.90·1027 kg y radio es R=6.98·106 m.

G=6.67e-11; %constante G
M=1.90e27; %masa de Júpiter
R=6.98e7; %radio de Júpiter

b=2*R; %parámetro de impacto
v0=14.6*1000; %velocidad
ex=sqrt(1+b^2*v0^4/(G*M)^2); %excentricidad
d=(b*v0)^2/(G*M);

rp=d/(1+ex); %máximo acercamiento
vp=v0*b/rp; %velocidad
fprintf('Máximo acercamiento %1.3f, velocidad (km/s) %3.2f\n'
,rp/R, vp/1000)
Máximo acercamiento 0.232, velocidad (km/s) 126.048

Hipérbola girada

Vamos a determinar la ecuación de la trayectoria de una partícula que se mueve a una distancia muy grande del centro de fuerzas paralelamente al eje X, con velocidad v0, tal como se muestra en la figura.

El momento angular L y la energía E valen

L=mb v 0 E= 1 2 m v 0 2

La ecuación de la trayectoria es

r= d 1+εcos( θ θ 0 )

Calculamos el ángulo θ0 sabiendo que cuando la partícula está muy alejada del origen θ

cos( θ θ 0 )= 1 ε ( d r 1 ) rθπ cos( π θ 0 )= 1 ε cos θ 0 = 1 ε

La ecuación de la trayectoria en términos del parámetro adimensional p es

r= d 1+cosθ+ ε 2 1 sinθ r= bp 1+cosθ+psinθ p= b v 0 2 GM

Donde π<θ>π-2α, siendo 2α-π el cambio en la dirección de la velocidad

α=arccos( 1 ε )=arccos( 1 p 2 +1 )

p=1.5; %parámetro adimensional
b=3; %parámetro de impacto
r=@(x) b*p./(1+cos(x)+p*sin(x));
angulo=2*acos(-1/sqrt(1+p^2));
theta=pi-20*pi/180:-pi/180:pi-angulo+20*pi/180;
%trayectoria
hold on
plot(r(theta).*cos(theta),r(theta).*sin(theta))
%origen, centro de fuerzas
plot(0,0,'ro', 'markersize',8,'markerfacecolor','y')
hold off
axis equal
grid on
xlabel('x')
ylabel('y')
title('Hipérbola girada')

Actividades

Se elige un planeta en el control titulado Planeta

Se introduce

Se pulsa el botón titulado Nuevo

La trayectoria del asteroide empieza a dibujarse a una distancia r=10·R, diez veces el radio del planeta. En la parte superior derecha, se nos proporciona los datos del tiempo en horas y la distancia en radios del planeta elegido.

Se sugiere al lector:

Datos

Planeta GM Radio (km)
Tierra 3.986·1014 6378
Marte 4.283·1013 3397
Júpiter 1.267·1017 71 492
Saturno 3.793·1016 60 268

Referencias

Van Allen J. A., Encounter of an asteriod with a planet. Am. J. Phys. 74 (8) August 2006, pp. 717-719