Por qué importa
En órbita baja, la vida de un satélite pequeño la fija la atmósfera más que el diseño. Sin propulsión, la vida es a la vez una restricción, porque la misión debe terminar antes de que la nave caiga, y una obligación, porque la nave debe abandonar la órbita en un plazo razonable tras el fin de la misión. Esta nota deriva el decaimiento por arrastre de las órbitas circulares con el modelo coherente más sencillo, el que usa la revisión sobre PocketQubes[1], compara factores de forma y contrasta el resultado con reentradas observadas.
Arrastre y parámetro balístico
La deceleración por arrastre crece con la densidad y con el cuadrado de la velocidad respecto al aire:
con la densidad , el coeficiente de arrastre , la sección transversal y la masa . La combinación es el parámetro balístico; su inverso, , el . En el flujo molecular libre de se suele suponer [2]. Una nave sin control de actitud presenta una sección cambiante; para un cuerpo convexo en rotación aleatoria, el área proyectada media es la cuarta parte de su superficie total (teorema de Cauchy).1Un cubo de 50 mm tiene 150 cm² de superficie y un área proyectada media de 37,5 cm². De frente al flujo presentaría 25 cm²; a lo largo de una diagonal del cuerpo, 43,3 cm².
Con estas hipótesis, el parámetro balístico de un 1P a plena carga es 0,033 m²/kg y el de un 1U a plena carga, 0,018 m²/kg. Los PocketQubes tienen poca masa por unidad de área[3], y esa es toda la razón de la diferencia de vida.
La ecuación del decaimiento
La energía orbital específica de una órbita circular es . El arrastre extrae energía al ritmo , de modo que
y con el semieje mayor decae como
El ritmo es lineal en la densidad y en el parámetro balístico: reducir a la mitad duplica aproximadamente la vida en la misma atmósfera.
Una atmósfera gobernada por el Sol
Las tablas exponenciales estáticas, como las de Vallado[2], ignoran el principal motor de la densidad termosférica: el Sol. El modelo de Kennewell[4] conserva una sola , pero la liga a la actividad solar y geomagnética:
| Altitud (km) | Baja | Moderada | Alta | Altura de escala, moderada (km) | Modelo exponencial estático | Alta ÷ baja |
|---|---|---|---|---|---|---|
| 300 | 1,7 × 10⁻¹¹ | 3,1 × 10⁻¹¹ | 4,8 × 10⁻¹¹ | 42 | 2,4 × 10⁻¹¹ | 3× |
| 400 | 1,4 × 10⁻¹² | 3,8 × 10⁻¹² | 7,9 × 10⁻¹² | 44 | 3,7 × 10⁻¹² | 6× |
| 500 | 1,4 × 10⁻¹³ | 5,7 × 10⁻¹³ | 1,6 × 10⁻¹² | 47 | 7,0 × 10⁻¹³ | 11× |
| 600 | 1,9 × 10⁻¹⁴ | 1,1 × 10⁻¹³ | 3,7 × 10⁻¹³ | 49 | 1,5 × 10⁻¹³ | 20× |
| 700 | 3,2 × 10⁻¹⁵ | 2,5 × 10⁻¹⁴ | 1,1 × 10⁻¹³ | 52 | 3,6 × 10⁻¹⁴ | 33× |
El modelo estático (Vallado, ajuste de CIRA-72) ignora la actividad solar. La diferencia entre actividad baja y alta crece con la altitud y supera un orden de magnitud por encima de unos 450 km. Actividad solar: baja F10.7 = 70 sfu, Ap = 6; moderada 140 sfu, Ap = 12; alta 210 sfu, Ap = 18.
Con actividad moderada, la altura de escala es de unos 42 km a 300 km y de 47 km a 500 km; entre actividad baja ( sfu, ) y alta (210 sfu, 18) varía entre unos 35 y 55 km, y la densidad a 500 km cambia en un factor diez.
Si la altura de escala y se mantienen constantes, la ecuación se integra analíticamente: el tiempo para descender una altura de escala desde es, en la forma más simple del análisis clásico de King-Hele[5],
La mayor parte de la vida transcurre en las primeras alturas de escala por debajo de la altitud inicial; por eso esta estimación tosca queda a unas decenas por ciento del resultado numérico.
Integración
Integramos la ecuación del decaimiento con un esquema de punto medio (Runge–Kutta de segundo orden), limitando cada paso a un día y a 100 m de pérdida de altitud, y suponemos la reentrada a 180 km, por debajo de la cual la vida restante se mide en días. Con estos ajustes se reproducen las vidas de la revisión con la precisión que allí se indica[1].
- 1P: 20,4 meses
- 3P: 22,6 meses
- 1U: 3,1 años
- 3U: 3,9 años
Mostrar los datos como tabla
| Tiempo (años) | 1P: 20,4 meses | 3P: 22,6 meses | 1U: 3,1 años | 3U: 3,9 años |
|---|---|---|---|---|
| 0,0 | 500 | 500 | 500 | 500 |
| 0,5 | 482 | 484 | 491 | 493 |
| 1,0 | 454 | 461 | 480 | 485 |
| 1,5 | 394 | 420 | 466 | 475 |
| 2,0 | — | — | 447 | 463 |
| 2,5 | — | — | 418 | 448 |
| 3,0 | — | — | 343 | 426 |
| 3,5 | — | — | — | 389 |
| 4,0 | — | — | — | — |
| Factor de forma | Masa (kg) | Parámetro balístico (m²/kg) | Desde 350 km | Desde 400 km | Desde 450 km | Desde 500 km | Desde 550 km | Desde 600 km |
|---|---|---|---|---|---|---|---|---|
| 1P | 0,25 | 0,0330 | 28 días | 3 meses | 8 meses | 20,4 meses | 4,3 años | 10,1 años |
| 2P | 0,50 | 0,0306 | 31 días | 3 meses | 8 meses | 22,0 meses | 4,6 años | 10,9 años |
| 3P | 0,75 | 0,0298 | 32 días | 3 meses | 9 meses | 22,6 meses | 4,7 años | 11,2 años |
| 1U | 2,00 | 0,0180 | 52 días | 5 meses | 14,2 meses | 3,1 años | 7,8 años | 18,6 años |
| 3U | 6,00 | 0,0143 | 2 meses | 6 meses | 17,8 meses | 3,9 años | 9,8 años | 23 años |
Modelo de densidad de Kennewell, coeficiente de arrastre 2,2, área media en rotación aleatoria A_tot/4, reentrada a 180 km, nave en el límite de masa de su estándar. Actividad solar: baja F10.7 = 70 sfu, Ap = 6; moderada 140 sfu, Ap = 12; alta 210 sfu, Ap = 18.
Destacan dos resultados. Los parámetros balísticos de los PocketQubes 1P, 2P y 3P son casi idénticos (0,030–0,033 m²/kg), porque masa y área crecen juntas a lo largo del eje largo. Y los CubeSats en su límite de masa viven mucho más que los PocketQubes. Un 1U más ligero, como un diseño típico de 1,33 kg, queda más cerca de los PocketQubes; las filas de CubeSat son cotas superiores.
Vida frente a altitud
- 1P, 0,033 m²/kg
- 3P, 0,030 m²/kg
- 1U, 0,018 m²/kg
- 3U, 0,014 m²/kg
- 1P, de actividad baja a alta
- 1U, de actividad baja a alta
Mostrar los datos como tabla
| Altitud inicial (km) | 1P, 0,033 m²/kg | 3P, 0,030 m²/kg | 1U, 0,018 m²/kg | 3U, 0,014 m²/kg |
|---|---|---|---|---|
| 350 | 0,08 | 0,09 | 0,14 | 0,18 |
| 400 | 0,23 | 0,26 | 0,42 | 0,53 |
| 450 | 0,64 | 0,71 | 1,18 | 1,49 |
| 500 | 1,70 | 1,89 | 3,12 | 3,92 |
| 550 | 4,26 | 4,72 | 7,82 | 9,82 |
| 600 | 10,13 | 11,23 | 18,59 | 23,35 |
| 650 | 22,87 | 25,35 | 41,97 | 52,72 |
Como la vida es aproximadamente inversa al parámetro balístico, desde la misma altitud un PocketQube vive más o menos la mitad que un 1U. Ambas curvas suben con fuerza con la altitud: cada 50 km multiplican la vida por dos o tres.
El Sol domina la incertidumbre
| Altitud (km) | Baja | Moderada | Alta | Baja ÷ alta |
|---|---|---|---|---|
| 400 | 6 meses | 3 meses | 47 días | 4,2× |
| 450 | 22,0 meses | 8 meses | 4 meses | 5,9× |
| 500 | 5,8 años | 20,4 meses | 9 meses | 8,1× |
| 550 | 17,4 años | 4,3 años | 19,2 meses | 10,9× |
| 600 | 49 años | 10,1 años | 3,4 años | 14,5× |
Modelo de densidad de Kennewell, coeficiente de arrastre 2,2, área media en rotación aleatoria A_tot/4, reentrada a 180 km, nave en el límite de masa de su estándar. Actividad solar: baja F10.7 = 70 sfu, Ap = 6; moderada 140 sfu, Ap = 12; alta 210 sfu, Ap = 18.
El mismo 1P liberado a 500 km vive casi seis años con actividad baja y nueve meses con actividad alta. Esta dispersión es mucho mayor que la incertidumbre de la masa o del coeficiente de arrastre. Una previsión de vida debe seguir la fase del ciclo solar en el lanzamiento y durante la misión, y las propias previsiones del ciclo cambian: el panel NOAA/NASA de 2019 esperaba que el ciclo 25 alcanzara su máximo con un número de manchas de 115 en julio de 2025[6], mientras que la actualización de 2023 adelantó el máximo a 2024 con 137–173[7]. Los modelos empíricos como NRLMSISE-00[8], alimentados con esas previsiones, son las herramientas de los análisis de cumplimiento; los estudios exhaustivos del decaimiento de CubeSats llegan a la misma conclusión sobre la actividad solar[9].
Contraste con el vuelo
Transporter-3 liberó dos PocketQubes bien documentados cerca de 525 km el 13 de enero de 2022. Grizu-263A, un 1P de 237 g, reentró el 3 de octubre de 2023 tras 628 días; Delfi-PQ, un 3P de 545 g, el 6 de enero de 2024 tras 723 días[10].
| Escenario de parámetro balístico | Parámetro balístico (m²/kg) | 100 sfu constantes | 141 sfu constantes | Flujo observado (revisión) |
|---|---|---|---|---|
| Solo el cubo en rotación libre (C_D 2,2) | 0,035 | 1795 d | 920 d | 845 d |
| Cubo con cinta de antena de 5 mm (nominal) | 0,039 | 1611 d | 825 d | 774 d |
| Cota superior: C_D 2,4, cinta de 10 mm | 0,047 | 1337 d | 685 d | 666 d |
| Valor efectivo que reproduce la reentrada (B_eff ≈ 0,050) | 0,050 | 1252 d | 641 d | 628 d |
Las columnas de flujo constante se recalculan aquí con Ap = 12; 141 sfu es la media de la misión. La última columna es la reconstrucción de la revisión con el F10.7 mensual observado (sección 5.8); el valor efectivo reproduce por construcción la reentrada del 3 de octubre de 2023.
Para Grizu-263A, un flujo medio constante de 141 sfu con el parámetro balístico nominal da unos 825 días, un 31 % más de lo observado; la reconstrucción de la revisión con el flujo mensual observado lo reduce a 774 días, y un parámetro balístico efectivo de ≈0,050 m²/kg (0,0502 m²/kg) reproduce la reentrada. Para Delfi-PQ, el valor nominal en rotación aleatoria de 0,041 m²/kg a 141 sfu da unos 790 días frente a los 723 observados. Ambos residuos tienen el mismo signo y causas plausibles: una sola altura de escala, la órbita exacta de liberación, apéndices como las cintas de antena que aumentan el área de arrastre y una actitud que no es una rotación aleatoria uniforme. Los historiales de que alimentan estas comparaciones también se degradan en los objetos ligeros[11].
Flujo de cálculo
flowchart LR A["Factor de forma: dimensiones, masa"] --> B["Área media A = A_tot / 4"] B --> C["Parámetro balístico C_D A / m"] D["Densidad de Kennewell: F10.7, Ap"] --> F["da/dt = −√(μa) ρ C_D A / m"] C --> F F --> G["RK2, paso ≤ 1 día y ≤ 100 m"] G --> H["Reentrada a 180 km: vida"]
Conclusión
Un modelo de primer orden con una densidad gobernada por el Sol recoge lo que importa para el diseño: desde la misma altitud, los PocketQubes viven aproximadamente la mitad que los CubeSats, y gran parte de la incertidumbre la fija el Sol, no la nave. Con flujo constante, el modelo sobrestima las vidas observadas de Transporter-3 en el sentido que cabe esperar de sus simplificaciones. La calculadora de vida orbital ejecuta el mismo modelo para cualquier nave, y el artículo sobre la regla de cinco años lo convierte en límites de altitud.