Estoy trabajando en una herramienta para simular trayectorias de cometas. Solo el método habitual en el que coloca el sol y los planetas en sus órbitas y luego resuelve las trayectorias de los cometas utilizando la ley de gravitación (quizás teniendo en cuenta la sublimación, etc. para una mejor precisión). Hay sitios web que te darán los elementos orbitales correspondientes a un cometa, con alguna época. Sin embargo, no se especifica dónde se encuentra el cometa en la época. Pregunta:
Me gustaría colocar un cometa conocido en sus efemérides observadas a una distancia considerable del sol en su entrada, para poder observarlo entrar, girar y salir. ¿Cómo puedo encontrar ese puesto usando datos en las bases de datos disponibles?
( Este es el algoritmo que uso para determinar la posición de los planetas, creo que también podría ser útil en este caso).
El grupo de dinámica del sistema solar del JPL tiene una herramienta en línea llamada HORIZONTES que se ha convertido en la herramienta de facto para buscar o calcular las efemérides de los planetas, los satélites naturales y varios miles de planetas menores. También proporcionan el kit de herramientas SPICE , que puede usar en su propio software. Esta herramienta necesita algunos archivos de datos generados previamente (como DE405/DE406 para planetas) y le permite generar efemérides precisas para todos y cada uno de los cuerpos en esos archivos de datos.
Luego está el Centro de Planetas Menores , el organismo de la IAU responsable de todos los planetas menores (clasificación, denominación, etc.). Ofrecen una herramienta llamada Servicio de Efemérides de Planetas Menores (MPES) , que le permite generar efemérides para una lista de hasta 100 planetas menores a la vez. Luego ofrecen la base de datos MPC Orbit , que proporciona archivos MPCORB.dat (o diarios), que contienen los elementos orbitales actuales de todos los planetas menores conocidos (y/o los recién descubiertos). Tenga en cuenta que, a diferencia del sitio al que se vinculó, el archivo MPCORB.dat incluye los 6 elementos orbitales y la época, lo que significa que puede determinar de forma única su posición y velocidad en el espacio.
Si mal no recuerdo, SPICE o MPC también proporcionan valores de deriva, es decir, las tasas de cambio de primer y segundo orden de los elementos orbitales en la época. Significado: más precisión proveniente de la herramienta de efemérides, se requieren menos cálculos de su parte. En general: use el trabajo que ya se ha hecho, en lugar de hacerlo usted mismo :)
Si realmente necesita o quiere propagar los asteroides: le sugiero que use los métodos de Encke en combinación con un integrador simpléctico o un integrador de Runge-Kutta-Nyström de alto orden .
Tenga en cuenta que estos integradores no se pueden usar si desea tener en cuenta cosas que dependen de la velocidad instantánea, como la presión de radiación. Pero bueno, si vas a entrar en ese tipo de detalles, bueno... tal vez sea mejor que hagas una nueva pregunta sobre qué libros leer primero :)
marca adler
TPel tiempo en el periapsis. Por ejemplo, para el cometa Halley, veráTP= 2446467.3953170511yTP= 1986-Feb-05.8953170511. Con el tiempo en el periapsis, puedes encontrar la posición del cometa en cualquier momento. También puede solicitar "vectores" en ese sistema y obtener la posición y la velocidad directamente en los tiempos solicitados. Hoy, el cometa Halley está en-2.044172910106513E+01 2.492823871168252E+01 -9.730954540464577E+00, en AU.