Predicting spacecraft surface degradation under atomic oxygen in Low Earth Orbit

Technische Samenvatting: Voorspellen van Oppervlakte-degradatie van Ruimtevaartuigen onder Invloed van Atoomzuurstof in een Lage Aardbaan

Probleemstelling
Nauwkeurige aerodynamische modellering voor thermosferisch onderzoek en satellietoperaties in een Lage Aardbaan (LEO) wordt gehinderd door gebrekkige kennis over de oppervlakte-degradatie veroorzaakt door atoomzuurstof (O1). Hoewel in-situ experimenten zoals de Materials International Space Station Experiments (MISSE) 2–8 de erosieopbrengsten en microscopische kegelvormige ruwheid van diverse polymeren hebben gedocumenteerd, ontbrak een theoretisch kader gebaseerd op eerste principes dat de chemische samenstelling koppelt aan deze macroscopische morfologische veranderingen. Bestaande empirische modellen passen in-situ gegevens aan, maar kunnen het gedrag van niet-geteste materialen niet voorspellen of de vorming van textuur verklaren. Omgekeerd schieten standaard reactieve Moleculaire Dynamica (MD) simulaties vaak tekort in het matchen met orbitale metingen vanwege beperkingen in de tijdschaal, oververhitting-artefacten en het negeren van de evoluerende topografie van het oppervlak, die de lokale invalshoeken verandert.

Methodologie
De auteurs stellen een multiscale framework gebaseerd op eerste principes voor die de brug slaat tussen chemie op atomaire schaal en macroscopische oppervlakte-evolutie. De aanpak bestaat uit twee gekoppelde componenten:

1. Reactieve Moleculaire Dynamica (MD) Simulaties:

   • Opzet: Uitgebreide simulaties werden uitgevoerd op vier polymeren (Kapton H, Teflon ETFE, Teflon FEP en Teflon PTFE) met behulp van het ReaxFF reactieve krachtveld. De simulaties bestreken lokale invalshoeken ($0^\circ$–$75^\circ$), impactsnelheden (3–10 km/s) en oppervlakte-temperaturen (200–400 K).
   • Thermostatering: Om de veelvoorkomende problemen met oververhitting bij korte MD-tijdschalen aan te pakken, werd een Nosé-Hoover chain thermostat toegepast tijdens de integratie van oppervlakte-atomen. Dit handhaafde een constante oppervlakte-temperatuur (het nabootsen van de trage thermische re-equilibratie van echte satellieten) zonder de hyperthermische gas-oppervlakte botsingsdynamica te verstoren.
   • Transiënte Filtering: Simulaties identificeerden een transiënt regime waarin losjes gebonden fragmenten worden uitgestoten. Gegevens uit deze initiële fase werden uitgesloten om het steady-state oxidatieve erosie-regime te isoleren.
   • Outputs: De MD-simulaties leverden drie lokale responsfuncties op die als input dienen voor het macro-model: lokale erosieopbrengst ($E_y$), reactiefractie ($\Psi$) en energie-accommodatiecoëfficiënt ($\alpha_E$).

2. Statistisch Macro-schaal Erosiemodel:

   • Framework: De lokale responsfuncties werden geïntegreerd in een statistisch erosiemodel gebaseerd op de poly-Gaussiaanse oppervlaktebeschrijving en de Level-Set Methode (LSM).
   • Mechanica: Het model behandelt oppervlakte-evolutie op twee schalen: de interactie op atomaire schaal en de macroscopische ruwheid gedreven door "as"-deeltjes (erosiebestendige residuen). Het houdt rekening met zowel eerste-reflectie als multi-reflectie deeltjesfluxen, waarbij de geometrie van het oppervlak in de loop van de tijd wordt bijgewerkt om langdurige blootstelling te simuleren (ongeveer 4 jaar onder MISSE 2-condities).
   • Validatie: Het model werd gevalideerd tegenover MISSE 2 vluchtgegevens, inclusclusief Field Emission Scanning Electron Microscopy (FESEM) beelden en gemeten erosieopbrengsten.

Belangrijkste Bijdragen

• Voorspelling op basis van Eerste Principes: De studie legt een direct verband tussen de chemische samenstelling van een materiaal en de voorspelde erosieopbrengst en oppervlaktemorfologie, zonder afhankelijk te zijn van empirische aanpassing aan in-situ gegevens.
• Thermostaterend MD-protocol: De introductie van een specifiek thermostaterend protocol maakt langere bombardementsduur in MD mogelijk, wat kunstmatige oververhitting voorkomt en meer realistische erosiestatistieken mogelijk maakt.
• Morfologievoorspelling: Het framework voorspelt succesvol de ontwikkeling van kegelvormige oppervlaktekenmerken en hun aspectratio's, waarbij de overgang van gladde naar ruwe oppervlakken wordt gevangen.
• Correlatie tussen Ruwheid en Erosie: De auteurs identificeren en kwantificeren een sterke correlatie tussen de globale erosieopbrengst en de steady-state oppervlakte-ruwheid (specifiek de standaarddeviatie van de oppervlaktehellingen).

Resultaten

• Nauwkeurigheid van de Erosieopbrengst: Het framework reproduceerde de gemeten globale erosieopbrengsten voor Kapton H, Teflon ETFE, FEP en PTFE met een maximale relatieve fout van 19% (gemiddelde fouten variërend van 4% tot 19%). Deze prestatie overtreft aanzienlijk eerdere empirische modellen (die fouten van 29–36% hadden voor dezelfde materialen).
• Morfologische Overeenkomst: De gesimuleerde oppervlaktemorfologieën kwamen kwalitatief overeen met de FESEM-beelden van MISSE 2, waarbij de kegelvormige kenmerken, aspectratio's en specifieke texturen (bijv. de "gesmolten" look van ETFE en de intrinsieke textuur van PTFE) correct werden gereproduceerd.
• Transiënte Regimes: De studie identificeerde twee verschillende erosieregimes:
  1. Materialen met hoge kristalliniteit (bijv. PTFE, FEP): Vertonen een initiële hoge erosieopbrengst die afneemt naarmate de ruwheid zich ontwikkelt en het oppervlak afschermt, om uiteindelijk een steady-state te bereiken.
  2. Materialen met lage kristalliniteit/reactieve materialen (bijv. Kapton H): Vertonen andere afhankelijkheden van de invalshoek, wat leidt tot verschillende ruwheidsontwikkelingspatronen.
• Ruwheidsvoorspelling: Met behulp van de afgeleide correlatie tussen erosieopbrengst en ruwheid, schatten de auteurs de ruwheidsparameters voor alle 41 in de MISSE 2 geteste materialen. Ze voorspellen dat zelfs materialen met lagere erosieopbrengsten een significante ruwheid ontwikkelen ($\sigma_s \geq 1.5$), wat suggereert dat de meeste LEO-oppervlakken voornamelijk diffuse gasreflectie zullen bevorderen.

Betekenis
Het artikel beweert dat dit framework een robuuste, op fysica gebaseerde tool biedt voor het voorspellen van de langetermijndegradatie van ruimtevaartuigenoppervlakken. Door zowel de erosieopbrengst als de evoluerende geometrische ruwheid nauwkeurig te vangen, adresseert het model een cruciale bron van onzekerheid in thermosferisch onderzoek: de modellering van gas-oppervlakte-interacties. De auteurs stellen dat de voorspelde ruwheid impliceert dat aerodynamische weerstandscoëfficiënten voor geërodeerde oppervlakken gemodelleerd moeten worden met diffusie-reflectie aannames (DRIA-model) met beperkte variatie tussen verschillende polymeercoatings. Dit werk biedt een pad naar verbeterde voorspelling van satellietbanen en materiaalselectie voor Very Low Earth Orbit (VLEO) missies zonder uitsluitend afhankelijk te zijn van kostbare en tijdrovende in-situ blootstellingsexperimenten.