¿Métodos o técnicas algorítmicas para encontrar conjunciones en conjuntos de elementos Keplerianos de alto N?

Supongamos que quisiera responder a la pregunta ¿Pasará Starman/Roadster particularmente cerca de algún asteroide en los próximos años? o tratar de predecir conjunciones de satélites alrededor de la Tierra (por ejemplo, SÓCRATES de Celestrak ), y tenía elementos orbitales keplerianos simples, no necesariamente tablas de elementos osculantes ni derivados por primera vez de los elementos (¡aunque una respuesta sobre eso sería genial!)

Podría propagar esas órbitas keplerianas simples de un cuerpo en pequeños pasos de tiempo, calcular todo norte puestos y todo norte ( norte 1 ) / 2 distancias y buscar cualquiera por debajo de una distancia d C o norte j , pero esa podría no ser la forma más eficiente de hacerlo.

Pregunta: ¿Cuáles son los métodos o técnicas algorítmicas para hacer este tipo de búsqueda de manera más eficiente? Necesito una explicación o una referencia autorizada, no solo un nombre.

¿Esta pregunta es distinta de los métodos o técnicas algorítmicas para encontrar conjunciones en grandes conjuntos de vectores de estado? porque pregunta específicamente sobre métodos que operan sobre elementos keplerianos. Por ejemplo, en el caso de elementos keplerianos simples, puede descartar inmediatamente algunas combinaciones al principio si el periápside de uno es más de d C o norte j más grande que la apoapsis del otro.

Estoy llamando a una efemérides más que simplemente una "tabla interpolable de vectores de estado". Por lo general, una efemérides se optimizará para un algoritmo de interpolación específico o incluso solo será útil dentro de un entorno de software específico, mientras que una tabla de vectores de estado es solo eso.

Respuestas (1)

Propagador norte objetos y control norte ( norte 1 ) / 2 distancias en cada pequeño paso de tiempo es computacionalmente ineficiente hasta el punto de ser poco práctico. Si no recuerdo mal, actualmente hay ~17000 objetos catalogados en SPACE-TRAK . Suponiendo que verifiquemos las colisiones cada segundo,* tendríamos que calcular 10 13 Normas euclidianas en un solo día. Algún tipo de filtrado es absolutamente necesario.

Creo que el trabajo de Hoots, Crawford y Roerich (1984) es uno de los primeros y más citados enfoques del problema. No estoy seguro de los detalles, pero creo que consta de 3 pasos, con un grado creciente de refinamiento:

  1. Un filtro de apogeo-perigeo, en el que se excluyen las órbitas concéntricas,
  2. Un filtro geométrico, en el que se comprueba que la distancia mínima entre las órbitas está por debajo de un determinado umbral,
  3. Un filtro de tiempo, en el que se verifica que los objetos en las órbitas que pasan los filtros 1 y 2 estén cerca del mismo punto al mismo tiempo.

Es probable que un filtro de este tipo siga dando muchos falsos positivos, y habrá que comprobarlos calculando la distancia euclidiana entre ellos. Otros factores a considerar son A) que las órbitas keplerianas son una aproximación válida solo por un período de tiempo limitado desde la última determinación de la órbita, y generalmente se necesita una trayectoria precisa obtenida a través de la propagación numérica; y B) que se debe tener en cuenta la incertidumbre en la posición/velocidad de los objetos, determinando en última instancia la probabilidad de colisión.

De todos modos, creo que esta es una línea de investigación activa todavía hoy. Puede encontrar más enfoques en las citas de Hoots et al. paper , pero creo que la mayoría de ellos todavía se reducen a enfoques de filtrado secuencial mejorados (que se pueden realizar en paralelo para una gran cantidad de objetos).


* Además, 1 segundo es demasiado grande para comprobar si hay colisiones. Las naves espaciales LEO tienen velocidades relativas del orden de 10 km/s, por lo que usar un intervalo de 1 segundo perdería la mayoría de las colisiones.

¡Gran! Voy a darles un buen vistazo hoy, gracias!
@uhoh, ¿puede dejarnos saber sus comentarios sobre esto? Estoy muy interesado en la respuesta, pero actualmente no tengo tiempo para investigar. Gracias
@ChrisR Estuve ocupado los últimos días y fui negligente en mis "tareas de intercambio de pilas". Con suerte mañana...
No he podido ir a la biblioteca para obtener una copia de Hoots (1984), así que no lo aceptaré todavía, pero esta es una gran respuesta, ¡gracias de nuevo!
Tengo Hoots (1984) pronto lo comentaré pero esto tiene buena pinta. ¡Gracias!