Predicting spacecraft surface degradation under atomic oxygen in Low Earth Orbit

Resumo Técnico: Predição da Degradação de Superfícies de Naves Espaciais sob Oxigênio Atômico em Órbita Terrestre Baixa

Definição do Problema
A modelagem aerodinâmica precisa para pesquisa termosférica e operações de satélites em Órbita Terrestre Baixa (LEO) é dificultada pelo conhecimento precário da degradação de superfícies causada pelo oxigênio atômico (O1). Embora experimentos in-situ, como os experimentos Materials International Space Station Experiments (MISSE) 2–8, tenham documentado os rendimentos de erosão e o rugosidade cônica microscópica de vários polímeros, faltava um arcabouço teórico de primeiros princípios que vinculasse a composição química a essas mudanças morfológicas de macroescala. Os modelos empíricos existentes ajustam dados in-situ, mas não conseguem prever o comportamento de materiais não testados ou explicar a formação de textura. Por outro lado, simulações padrão de Dinâmica Molecular (MD) reativa frequentemente falham em corresponder às medições orbitais devido a limitações de escala de tempo, artefatos de superaquecimento e a negligência da evolução da topografia superficial, que altera os ângulos de incidência locais.

Metodologia
Os autores propõem um arcabouço multiescala de primeiros princípios que faz a ponte entre a química em escala atômica e a evolução da superfície em macroescala. A abordagem consiste em dois componentes acoplados:

1. Simulações de Dinâmica Molecular (MD) Reativa:

   • Configuração: Simulações extensivas foram realizadas em quatro polímeros (Kapton H, Teflon ETFE, Teflon FEP e Teflon PTFE) utilizando o campo de força reativo ReaxFF. As simulações cobriram ângulos de incidência local ($0^\circ$–$75^\circ$), velocidades de impacto (3–10 km/s) e temperaturas de superfície (200–400 K).
   • Termostatação: Para abordar os problemas de superaquecimento comuns em MD de curto prazo, um termostato de cadeia de Nosé-Hoover foi aplicado durante a integração dos átomos da superfície. Isso manteve uma temperatura de superfície constante (imitando o lento reequilíbrio térmico de satélites reais) sem perturbar a dinâmica de colisão gás-superfície hipertérmica.
   • Filtragem Transiente: As simulações identificaram um regime transiente onde fragmentos fracamente ligados são ejetados. Dados desta fase inicial foram excluídos para isolar o regime de erosão oxidativa em estado estacionário.
   • Saídas: As simulações de MD geraram três funções de resposta local que servem como entradas para o modelo macro: rendimento de erosão local ($E_y$), fração de reação ($\Psi$) e coeficiente de acomodação de energia ($\alpha_E$).

2. Modelo de Erosão Estatística em Macroescala:

   • Arcabouço: As funções de resposta local foram integradas em um modelo de erosão estatística baseado na descrição de superfície poli-Gaussiana e no Método de Nível de Conjunto (Level-Set Method - LSM).
   • Mecânica: O modelo trata a evolução da superfície em duas escalas: interações em escala atômica e rugosidade em macroescala impulsionada por partículas de "cinzas" (resíduos resistentes à erosão). Ele considera tanto o fluxo de partículas de primeira reflexão quanto de múltiplas reflexões, atualizando a geometria da superfície ao longo do tempo para simular a exposição prolongada (aprox. 4 anos de condições MISSE 2).
   • Validação: O modelo foi validado contra dados de voo do MISSE 2, incluindo imagens de Microscopia Eletrônica de Varredura de Emissão de Campo (FESEM) e rendimentos de erosão medidos.

Principais Contribuições

• Predição de Primeiros Princípios: O estudo estabelece um vínculo direto entre a composição química de um material e seu rendimento de erosão e morfologia de superfície previstos, sem depender de ajuste empírico a dados in-situ.
• Protocolo de MD Termostatizado: A introdução de um protocolo de termostatação específico permite durações de bombardeio mais longas em MD, prevenindo o superaquecimento artificial e permitindo estatísticas de erosão mais realistas.
• Predição de Morfologia: O arcabouço prevê com sucesso o desenvolvimento de características superficiais cônicas e suas razões de aspecto, capturando a transição de superfícies lisas para rugosas.
• Correlação Rugosidade-Erosão: Os autores identificam e quantificam uma forte correlação entre o rendimento de erosão global e a rugosidade da superfície em estado estacionário (especificamente o desvio padrão das inclinações da superfície).

Resultados

• Precisão do Rendimento de Erosão: O arcabouço reproduziu os rendimentos de erosão globais medidos para Kapton H, Teflon ETFE, FEP e PTFE com um erro relativo máximo de 19% (erros médios variando de 4% a 19%). Este desempenho excede significativamente os modelos empíricos anteriores (que apresentavam erros de 29–36% para os mesmos materiais).
• Concordância Morfológica: As morfologias de superfície simuladas corresponderam qualitativamente às imagens FESEM do MISSE 2, reproduzindo corretamente as características cônicas, as razões de aspecto e texturas específicas (ex: o aspecto "derretido" do ETFE e a textura intrínseca do PTFE).
• Regimes Transientes: O estudo identificou dois regimes de erosão distintos:
  1. Materiais de Alta Cristalinidade (ex: PTFE, FEP): Exibem um rendimento de erosão inicial alto que diminui conforme a rugosidade se desenvolve e protege a superfície, atingindo eventualmente um estado estacionário.
  2. Materiais de Baixa Cristalinidade/Reativos (ex: Kapton H): Mostram dependências diferentes do ângulo de incidência, levando a diferentes padrões de evolução de rugosidade.
• Predição de Rugosidade: Usando a correlação derivada entre rendimento de erosão e rugosidade, os autores estimaram os parâmetros de rugosidade para todos os 41 materiais testados no MISSE 2. Eles preveem que mesmo materiais com menores rendimentos de erosão desenvolvem rugosidade significativa ($\sigma_s \geq 1,5$), sugerindo que a maioria das superfícies em LEO promoverá predominantemente a reflexão difusa de gases.

Significância
O artigo afirma que este arcabouço fornece uma ferramenta robusta e baseada em física para prever a degradação de longo prazo de superfícies de espaçonaves. Ao capturar com precisão tanto o rendimento de erosão quanto a rugosidade geométrica evolutiva, o modelo aborda uma fonte crítica de incerteza na pesquisa termosférica: a modelagem da interação gás-superfície. Os autores afirmam que a rugosidade prevista implica que os coeficientes de arrasto aerodinâmico para superfícies erodidas devem ser modelados usando suposições de reflexão difusa (modelo DRIA) com variação limitada entre diferentes revestimentos poliméricos. Este trabalho oferece um caminho para melhorar a previsão de órbita de satélites e a seleção de materiais para missões de Órbita Terrestre Muito Baixa (VLEO) sem depender exclusivamente de experimentos de exposição in-situ caros e demorados.