next up previous contents
Next: Indicadores de caos Up: Tratamiento teórico Previous: Reducción del sistema   Contents

Formalismo del problema de Sitnikov

Como se había mencionado en la introducción (sec. 2.1), el problema de Sitnikov es una reducción del caso general. Para hacer más sencillo el cálculo numérico se considerarán unidades tales que \bgroup\color{Black}$\mathrm{G}=1$\egroup, y \bgroup\color{Black}$m=1$\egroup. Considerando que las masas de los cuerpos dominantes sean \bgroup\color{Black}$m_1=m_2=m/2$\egroup, tienen un movimiento elíptico en torno a su centro de masa, con excentricidad \bgroup\color{Black}$e$\egroup, y período \bgroup\color{Black}$2\pi$\egroup. Entonces el movimiento de un tercer cuerpo situado en el eje perpendicular al plano de movimiento de las dos masas (plano \bgroup\color{Black}$xy$\egroup) y que pase por el centro masa (por el origen), estará restringido sólo a ese eje (el eje \bgroup\color{Black}$z$\egroup).

Si ponemos \bgroup\color{Black}$r(t)$\egroup como la distancia de las dominantes al centro de masa, y con las consideraciones anteriores, podemos definir el hamiltoniano que describe a la masa en el eje \bgroup\color{Black}$z$\egroup en términos de unidades de masa de la siguiente forma:

\begin{displaymath}
H\left(z,v,t\right)=\frac{v^2}{2}-\frac{1}{\sqrt{z^2+r^2(t)}}
\end{displaymath} (9)

Del resultado de los dos cuerpos, se tiene que (considerando las unidades escogidas):
\begin{displaymath}
2r(\varphi)=\frac{1-e^2}{1+e\cos \varphi}=1-e\cos \varphi+\mathrm{O}\left(e^2\right)
\end{displaymath} (10)

Y el ángulo polar \bgroup\color{Black}$\varphi$\egroup se puede poner en función del tiempo \bgroup\color{Black}$t$\egroup, va la integral
\begin{displaymath}
t=\int_0^t\frac{4r^2(\varphi')}{\sqrt{1-e^2}}\mathrm{d}\varphi'
\end{displaymath} (11)

Entonces con las ecuaciones (10) y (11) y descartando términos de orden \bgroup\color{Black}$e^2$\egroup y mayores, obtenemos:

\begin{displaymath}
r(t)=\frac{1}{2}\left(1-e\cos \varphi\right)+\mathrm{O}\left(e^2\right)
\end{displaymath} (12)

as, el hamiltoniano (9) podrá escribirse a primer orden en \bgroup\color{Black}$e$\egroup, como:
\begin{displaymath}
H\left(z,v,t\right)=\frac{v^2}{2}-\frac{1}{\sqrt{z^2+1/4}}-e\frac{\cos t}{4\left(z^2+1/4\right)^{3/2}}
\end{displaymath} (13)

Como se puede ver, este hamiltoniano depende del tiempo y tiene perodo \bgroup\color{Black}$2\pi$\egroup.

Las ecuaciones de movimiento para este hamiltoniano son:

$\displaystyle \dot{z}$ $\textstyle =$ $\displaystyle v$ (14)
$\displaystyle \dot{v}$ $\textstyle =$ $\displaystyle -\frac{8z}{\left(4z^2+1\right)^{3/2}}-e\frac{24z}{\left(4z^2+1\right)^{5/2}}\cos t$ (15)

Para \bgroup\color{Black}$e=0$\egroup el sistema se vuelve independiente del tiempo y por lo tanto es integrable, pero para \bgroup\color{Black}$e>0$\egroup la dependencia explcita del tiempo hace que no rompa la simetra temporal y no habrán cantidades integrables convirtiendo al problema en no integrable.


next up previous contents
Next: Indicadores de caos Up: Tratamiento teórico Previous: Reducción del sistema   Contents
Francisco Javier Rodríguez Arias 2004-12-16