¿Cuándo (aproximadamente) volverá a entrar el satélite Kwangmyongsong-4 de Corea del Norte?

Nadie puede decirlo con certeza, pero se puede usar el estándar ISO apropiado para hacer una estimación:

Sistemas espaciales — Estimación de la vida útil en órbita (ISO 27852) (ver página 20)

La figura B5 dice que un objeto en órbita polar con un perigeo de unos 500 km (tomado de su imagen) volverá a entrar en unos 25 años. Dependiendo de la cantidad de arrastre que produzca, podría ser más largo.

Si desea esforzarse un poco para obtener una respuesta más precisa, pruebe UTIAS Review 43 . Te dice que a 500 km de altitud, la longitud del camino libre de una molécula es de 77 km, por lo que no es tanto un gas, sino unas pocas moléculas perdidas. Solo por esto, podemos ver que no volverá a entrar "pronto". Para hacer un cálculo, necesitaría estimar el área de las superficies expuestas a la dirección orbital y la masa. La revisión habla específicamente de satélites giratorios. Si no puede hacer una buena estimación, puede intentar obtener los parámetros orbitales con la mayor precisión posible ahora y luego nuevamente en un año y modificar el coeficiente de arrastre para que se ajuste al cambio. Con esto puedes hacer una mejor estimación de la vida útil orbital.

Si haces esto, ¡por favor dinos lo que encontraste! Si lo publica como un comentario, lo incluiré en la respuesta.

Aquí hay una simulación para predecir la fecha de reingreso.

La simulación incluye las aceleraciones newtoniana y relativista de todos los planetas, Sol y Luna.
El campo de gravedad de la Tierra se modela con el modelo de gravedad SGG-UGM-1 (calculado utilizando la anomalía de gravedad derivada de EGM2008 y los datos de observación de GOCE) truncado al grado y orden 15 (para ahorrar el tiempo de ejecución, manteniendo una buena precisión en comparación con el total). modelo).
Para el cálculo de la densidad del aire, utilizo el modelo NRLMSISE-00 junto con un archivo de datos actualizado para los índices solar y geomagnético. Los índices reales se pueden encontrar aquí: www.celestrak.com/spacedata/SW-All.txt, mientras que el modelo utilizado para la predicción a largo plazo de los índices se puede encontrar aquí: https://www.nasa.gov/msfcsolar .

El primer paso consiste en determinar el mejor coeficiente balístico para minimizar un parámetro de simulación en particular. Después de 46 minutos, el programa encuentra un coeficiente balístico de aproximadamente 118 kg/m^2 (no es fijo, porque el coeficiente de arrastre varía con la composición del aire).
El siguiente gráfico muestra el resultado de los dos últimos meses:

vemos que el vector de radio medio integrado se ajusta muy bien al vector de radio medio obtenido con la biblioteca SGP4 de CSpOC para la época TLE.

Ahora se puede iniciar la simulación:
1) con un TLE y el propagador SGP4 calcular el estado inicial (posición y velocidad) del satélite para la época TLE;
2) propagar ese estado inicial con un propagador especialmente diseñado (mi propagador se basa en el integrador 8(5,3) Dormand-Prince);
3) cuando la altitud del satélite desciende por debajo de 70 km, detener la simulación; esta es la fecha de reingreso.

Aquí está el resultado obtenido con 7 TLE desde 20003.47248329 hasta 20010.45318683:

El gráfico muestra:
Vector de radio medio de 1 órbita (gráfico azul): vector de radio promediado sobre la anomalía excéntrica (es el semieje mayor). No debe confundirse con el semieje mayor osculador.

Vector de radio mínimo de 1 órbita (gráfico rojo): el vector de radio más pequeño en una órbita. No debe confundirse con el perigeo osculador.

Vector de radio máximo de 1 órbita (gráfico verde): el vector de radio más grande en una órbita. No debe confundirse con el apogeo osculador.

excentricidad real = (Ra - Rp) / (Ra + Rp) , donde Ra es el vector de radio máximo de 1 órbita y Rp es el vector de radio mínimo de 1 órbita.

La fecha de reingreso promedio es 2025-05-20 y la diferencia entre la última y la anterior es de solo 3,5 días (lo que significa que los TLE son muy precisos para este satélite). La única incógnita es la densidad del aire.

La respuesta que dio Rikki-Tikki-Tavi es correcta, con una pequeña excepción. El periapsis, no el apoasis, es la clave para usar tales gráficos con fines de reingreso. Debe utilizarse el periápside de 466 km. Usando la misma fuente, estimaría que la vida útil no es más de 5 años, siendo 2-3 más probable. Sospecho que el satélite estaba un poco bajo, ciertamente se comparó con el satélite anterior lanzado por Corea del Norte , y también sería consistente con su caída.