Predicting spacecraft surface degradation under atomic oxygen in Low Earth Orbit

Resumen Técnico: Predicción de la degradación de superficies de naves espaciales bajo oxígeno atómico en órbita terrestre baja

Planteamiento del problema
El modelado aerodinámico preciso para la investigación termosférica y las operaciones de satélites en órbita terrestre baja (LEO) se ve obstaculizado por el escaso conocimiento de la degradación de superficies causada por el oxígeno atómico (O1). Si bien los experimentos in situ, como los de los Experimentos de la Estación Espacial Internacional de Materiales (MISSE) 2–8, han documentado los rendimientos de erosión y el rugosidad cónica microscópica de varios polímeros, carecía de un marco teórico de primeros principios que vinculara la composición química con estos cambios morfológicos a macroescala. Los modelos empíricos existentes ajustan los datos in situ, pero no pueden predecir el comportamiento de materiales no probados ni explicar la formación de texturas. Por el contrario, las simulaciones estándar de Dinámica Molecular (MD) reactiva suelen fallar al intentar igualar las mediciones orbitales debido a las limitaciones de escala temporal, artefactos de sobrecalentamiento y la negligencia de la evolución de la topografía superficial, la cual altera los ángulos de incidencia locales.

Metodología
Los autores proponen un marco multiescala de primeros principios que une la química a escala atómica con la evolución de la superficie a macroescala. El enfoque consta de dos componentes acoplados:

1. Simulaciones de Dinámica Molecular (MD) reactiva:

   • Configuración: Se realizaron simulaciones extensas en cuatro polímeros (Kapton H, Teflon ETFE, Teflon FEP y Teflon PTFE) utilizando el campo de fuerza reactivo ReaxFF. Las simulaciones cubrieron ángulos de incidencia local ($0^\circ$–$75^\circ$), velocidades de impacto (3–10 km/s) y temperaturas superficiales (200–400 K).
   • Termostatización: Para abordar los problemas de sobrecalentamiento comunes en las MD de escalas temporales cortas, se aplicó un termostato de cadena de Nosé-Hoover durante la integración de los átomos superficiales. Esto mantuvo una temperatura superficial constante (imitando el lento reequilibrio térmico de los satélites reales) sin perturbar la dinámica de colisión gas-superficie hipertérmica.
   • Filtrado Transitorio: Las simulaciones identificaron un régimen transitorio donde se expulsan fragmentos débilmente ligados. Los datos de esta fase inicial se excluyeron para aislar el régimen de erosión oxidativa en estado estacionario.
   • Salidas: Las simulaciones de MD generaron tres funciones de respuesta local como entradas para el modelo macro: rendimiento de erosión local ($E_y$), fracción de reacción ($\Psi$) y coeficiente de acomodación de energía ($\alpha_E$).

2. Modelo de Erosión Estadístico a Macroescala:

   • Marco: Las funciones de respuesta local se integraron en un modelo de erosión estadística basado en la descripción de superficie poli-gaussiana y el Método de Conjunto de Niveles (LSM).
   • Mecánica: El modelo trata la evolución de la superficie en dos escalas: la interacción a escala atómica y la rugosidad macroscópica impulsada por partículas de "ceniza" (residuos resistentes a la erosión). Tiene en cuenta tanto el flujo de partículas de primera reflexión como el de reflexiones múltiples, actualizando la geometría de la superficie a lo largo del tiempo para simular la exposición prolongada (aprox. 4 años de condiciones de MISSE 2).
   • Validación: El modelo se validó con los datos de vuelo de MISSE 2, incluyendo imágenes de Microscopía Electrónica de Barrido de Emisión de Campo (FESEM) y rendimientos de erosión medidos.

Contribuciones Clave

• Predicción de Primeros Principios: El estudio establece un vínculo directo entre la composición química de un material y su rendimiento de erosión y morfología superficial predichos, sin depender del ajuste empírico a datos in situ.
• Protocolo de MD Termostatizado: La introducción de un protocolo de termostatización específico permite duraciones de bombardeo más largas en MD, evitando el sobrecalentamiento artificial y permitiendo estadísticas de erosión más realistas.
• Predicción de la Morfología: El marco predice con éxito el desarrollo de características superficiales cónicas y sus relaciones de aspecto, capturando la transición de superficies lisas a rugosas.
• Correlación Rugosidad-Erosión: Los autores identifican y cuantifican una fuerte correlación entre el rendimiento de erosión global y la rugosidad de estado estacionario (específicamente la desviación estándar de las pendientes de la superficie).

Resultados

• Precisión del Rendimiento de Erosión: El marco reprodujo los rendimientos de erosión global medidos para Kapton H, Teflon ETFE, FEP y PTFE con un error relativo máximo del 19% (errores promedio entre el 4% y el 19%). Este rendimiento supera significativamente a los modelos empíricos previos (que tenían errores del 29–36% para los mismos materiales).
• Acuerdo Morfológico: Las morfologías superficiales simuladas coincidieron cualitativamente con las imágenes FESEM de MISSE 2, reproduciendo correctamente las características cónicas, las relaciones de aspecto y las texturas específicas (por ejemplo, el aspecto "derretido" de ETFE y la textura intrínseca de PTFE).
• Regímenes Transitorios: El estudio identificó dos regímenes de erosión distintos:
  1. Materiales de Alta Cristalinidad (ej. PTFE, FEP): Exhiben un rendimiento de erosión inicial alto que disminuye a medida que la rugosidad se desarrolla y protege la superficie, alcanzando finalmente un estado estacionario.
  2. Materiales de Baja Cristalinidad/Reactivos (ej. Kapton H): Muestran dependencias diferentes respecto al ángulo de incidencia, lo que conduce a diferentes patrones de evolución de la rugosidad.
• Predicción de la Rugosidad: Utilizando la correlación derivada entre el rendimiento de erosión y la rugosidad, los autores estimaron los parámetros de rugosidad para los 41 materiales probados en MISSE 2. Predicen que incluso los materiales con menores rendimientos de erosión desarrollan una rugosidad significativa ($\sigma_s \geq 1.5$), lo que sugiere que la mayoría de las superficies en LEO promoverán predominantemente la reflexión difusa del gas.

Significancia
El artículo afirma que este marco proporciona una herramienta robusta y basada en la física para predecir la degradación a largo plazo de las superficies de las naves espaciales. Al capturar con precisión tanto el rendimiento de la erosión como la rugosidad geométrica evolutiva, el modelo aborda una fuente crítica de incertidumbre en la investigación termosférica: el modelado de la interacción gas-superficie. Los autores sostienen que la rugosidad predicha implica que los coeficientes de arrastre aerodinámico para superficies erosionadas deben modelarse utilizando supuestos de reflexión difusa (modelo DRIA) con variaciones limitadas entre diferentes recubrimientos poliméricos. Este trabajo ofrece una vía para mejorar la predicción de la órbita de los satélites y la selección de materiales para misiones de Órbita Terrestre Muy Baja (VLEO) sin depender únicamente de experimentos de exposición in situ costosos y prolongados.