Dostępny powszechnie materiał masowy stanowił wysoce porowaty preformowany węgiel związany wiązkami węglowymi (CBCF), składający się z krótkich włókien izolacyjnych pochodzących z bawełny (włókno celulozowe wykonane z oczyszczonej celulozy). Pokrojone, nieciągłe włókna węglowe są ze sobą połączone w matrycy powstałej w wyniku karbonizacji żywicy fenolowej. W trakcie tego procesu włókna ulegają uporządkowaniu, a mikrostruktura i właściwości materiału stają się anizotropowe. Następnie materiał poddawany jest obróbce w próżni w temperaturach powyżej 2300 K, aby zagwarantować jego stabilność temperaturową oraz brak wydzielania gazów. Materiał został obrabiany wewnętrznie na półkuliste (HS) próbki testowe o promieniu 25 mm i długości 50 mm. Próbki mają początkową gęstość około 180 kg/m3 oraz początkową porowatość wynoszącą 90%.
Do wszystkich eksperymentów mających na celu odtworzenie aero-termodynamicznego środowiska strumieni plazmy podczas ponownego wejścia w atmosferę wykorzystano stanowisko VKI Plasmatron, generujące strumień gazowy o wysokiej entalpii i silnym stopniu dysocjacji, przepływający z prędkością subsoniczną. Gaz jest ogrzewany indukcyjnie za pomocą cewki, co pozwala uzyskać strumień plazmy o wysokiej czystości. Przegląd komory testowej oraz schemat aparatury do pomiarów ablacji in-situ przedstawiono na rys. 1(a) i 1(b). Warunki testowe oraz ogólne wyniki eksperymentów, takie jak uśredniona szybkość cofania się powierzchni próbki uzyskana na podstawie obrazowania HSC i utrata masy, zamieszczono w tabeli 1. Zastosowano pirometr dwukolorowy, wykorzystujący szerokie (0,75–1,1 µm) oraz wąskie (0,95–1,1 µm) pasmo spektralne do wyznaczania temperatury z częstotliwością próbkowania 1 Hz (1300–3300 K). Użycie dwóch wąskich pasm falowych oraz założenie, że emisyjność nie zależy od długości fali, pozwala oszacować temperaturę powierzchni bez konieczności znajomości jej emisyjności. Pirometr skierowano i sfokusowano na strefę stagnacji próbki przez 1 cm grubą kwarcową szybkę, pod kątem 35° względem linii stagnacji. Kalibrację urządzenia przeprowadzono do temperatury 3300 K za pomocą źródła ciała doskonale czarnego.
Ubytek powierzchniowy mierzono za pomocą HSC z rozdzielczością 0,2 mm. Wyraźnie widać, że pomiary ubytku linijką suwmiarkową dawały zazwyczaj większe wartości niż te przeprowadzone za pomocą obrazowania HSC, przy różnicy całkowitego ubytku między obiema metodami w zakresie od 0,45 do 0,9 mm. Największą niepewność pomiaru wprowadzało uciskanie kruchej warstwy węglowej linijką suwmiarkową. Szybkości ubytku w powietrzu wahają się od 44,6 do 58,4 µm/s. Ponadto widoczne jest, że szybkości ubytku określone metodą HSC w plazmie powietrznej nie różnią się znacząco, co najprawdopodobniej wynika z reżimu ablacji kontrolowanej dyfuzją. W tym reżimie temperatura powierzchni jest wystarczająco wysoka, aby spowodować całkowite zużycie dostępnego tlenu na powierzchni, wskutek czego ablacja jest ograniczona przez dyfuzję tlenu przez warstwę graniczną32,33. Z kolei w środowisku, w którym reakcja utleniania jest kontrolowana kinetycznie, tlen dyfunduje przez warstwę graniczną szybciej niż jest zużywany na powierzchni, a ablacja wzrasta wraz ze wzrostem temperatury powierzchni. Szybkości ubytku materiału CBCF w warunkach wysokiej entalpii zostały również zgłoszone przez MacDonalda et al. (56 µm/s)22 oraz Löhle’go et al. (50 µm/s)34. Te wartości mieszczą się pomiędzy naszymi pomiarami, choć MacDonald et al. użyli cylindrycznej próbki testowej, a Löhle et al. próbki testowej osadzonej w chłodzonej wodą sondzie.
Do obserwacji fazy gazowej wykorzystano trzy spektrometry o niskiej rozdzielczości. Zaletą tego instrumentu jest szybkie skanowanie szerokiego zakresu spektralnego (200–1 000 nm), które pozwala na wykrywanie wielu cząsteczek i atomów obecnych w analizie ablacji.
Zintegrowane natężenia emisji CN, przedstawione w funkcji odległości od powierzchni ablacyjnej, wykazują bardzo dobrą zgodność między sobą (Rys. 2). Dane są oznaczone zgodnie z ich odpowiednimi położeniami względem powierzchni próbki jako „bliskie”, „średnie” i „dalekie”. Trzy spektrometry zmierzyły takie samo natężenie emisji fioletowej linii CN, gdy stała ścieżka optyczna zbierała światło z tej samej odległości przed powierzchnią. Zintegrowane natężenia wszystkich trzech spektrometrów niemal się pokrywają w odległości 3,4 mm przed powierzchnią ablacyjną. Oba przypadki pokazują, że zarejestrowana emisja fioletowa CN osiągała maksimum tuż przed próbką testową, a następnie malała w warstwie granicznej. Na podstawie tych wyników wnioskujemy, że spalanie materiału w powietrzu przez cały czas trwania testu było bardzo stabilne, a sygnał emisji zmniejszył się o około 90% w odległości 5 mm przed powierzchnią. Następnie eksperymentalne widma fioletowej emisji CN zostały wykorzystane do porównania z widmami symulowanymi w celu wyznaczenia temperatur gazów. Widma syntetyczne uzyskano za pomocą programu SPECAIR 2.2, zakładając rozkład Boltzmanna poziomów wzbudzonych, a do oszacowania temperatur rotacyjno-postępowych Trot oraz temperatur wibracyjno-elektronowych Tvib zastosowano metodę najmniejszych kwadratów (Rys. 3). Przedstawiono dwa przypadki: przy niskim (a) i wysokim (b) ciśnieniu, przy czym widma zostały zarejestrowane blisko ściany w warstwie granicznej. Oszacowane temperatury wskazują na dużą rozbieżność z równowagą termiczną przy niskim ciśnieniu (Rys. 3(a)). Tę samą analizę przeprowadzono dla kilku odległości od powierzchni, co lepiej ilustruje odchylenie od równowagi termicznej w pobliżu ściany przy niskim ciśnieniu (Rys. 4(a), 15 hPa), która ulega wyrównaniu w miarę przechodzenia przez warstwę graniczną. Uzyskane temperatury miały wartości rzędu 8200 K dla Trot i 21 000 K dla Tvib w pobliżu ściany, przy czym Tvib zmniejszała się do wartości 8200 K w miarę przechodzenia przez warstwę graniczną. Jest to w przeciwieństwie do warunku równowagi panującego w całej warstwie granicznej przy wyższym ciśnieniu (Rys. 4(b), 200 hPa). Granice temperatur ustalono na podstawie niepewności pomiaru natężenia emisji spektrometru na poziomie 10%, co pozwala na teoretyczną zmienność widma w ramach tych granic podczas procedury dopasowania.
Przy niskim ciśnieniu przeniesienie wzbudzenia między cząsteczkami jest ograniczone ze względu na mniejszą liczbę zderzeń, co może tłumaczyć efekt wyrównywania w kierunku krawędzi warstwy granicznej. Zakładamy silny wpływ azotu cząsteczkowego przy niskich entalpiach plazmy na produkcję CN, po którym następuje wzbudzenie wibracyjne CN. Przyjmuje się, że dysocjacyjna adsorpcja silnie wibrującego azotu prowadzi do powstawania reaktywnych centrów na powierzchni, co z kolei skutkuje produkcją CN. Boubert i Vervisch opisują ten proces w plazmie azotu/dwutlenku węgla przy niskim ciśnieniu35. Proces ten może prowadzić do gromadzenia się atomów azotu na powierzchni, przy czym reakcje egzotermiczne powodują przekształcenie nadmiaru energii w wzbudzenie rotacyjne i wibracyjne CN.
Mikrografy wykazały, że utlenianie węgla w plazmie powietrznej prowadzi do powstawania kształtu soplek u ablacyjnych włókien z głębokością utlenienia około 0,2 mm (ryc. 5(a)). Taki rodzaj formowania się soplek w wyniku ablacji jest szeroko opisywany w literaturze dotyczącej materiałów kompozytowych węgiel-węgiel36-38. Kształt sopelki (kąt rozwarcia) zależy od konkurencji między reakcją a dyfuzją na powierzchni materiału porowatego i w związku z tym zmienia się wraz z dyfuzją tlenu. Przyjmuje się, że ta długość odpowiada średniej głębokości dyfuzji tlenu. Kształt sopelki dodatkowo potwierdza ablację kontrolowaną przez dyfuzję. W przeciwieństwie do tego, ablację ograniczoną przez reakcję charakteryzuje możliwość przepływu tlenu w głąb struktury włókien, co prowadzi do lokalnego powstawania ubytków w węglowych włóknach.
Podczas niektórych testów ablacji obserwowano intensywne iskrzenie (ryc. 5(b)), które mogło być spowodowane odrywaniem się gorących skupisk włókien od powierzchni. Ablacja w plazmie azotowej prowadziła do silnego degradowania włókien na całej ich powierzchni, co skutkowało powolnym cofaniem się materiału wskutek azotowania (ryc. 5(c)). Ponieważ reaktywność węgla wobec azotu jest znacznie niższa niż wobec tlenu, azot może przenikać głębiej w materiał, powodując degradację wzdłuż całego włókna.

Rycina 1. Palnik plazmowy i przegląd układu eksperymentalnego. (a) Ogólny widok komory testowej palnika plazmowego VKI, pokazujący próbkę testową poza systemem zatrzymywania, sondy strumienia ciepła i ciśnienia oraz dostęp optyczny dla radiometrów, HSC i optyki spektrometru. (b) Schemat układu eksperymentalnego. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Tabela 1. Warunki testów Plasmatronu i wyniki eksperymentalne próbek wstępnie uformowanego węgla. Odniesienie przypadku testowego, gaz testowy, ciśnienie statyczne ps, ciśnienie dynamiczne pd, moc generatora P, średnie natężenie strumienia cieplnego na zimnej ścianie qcw, czas ekspozycji próbki τ, średnia temperatura powierzchni Ts, szybkość cofania się powierzchni r/τ oraz szybkość ubytku masy m/τ.

Rysunek 2. Przestrzenny emisja fioletowa CN w warstwie granicznej. Profile emisji zarejestrowane przez trzy sąsiednie spektrometry podczas ablacji preforma w powietrzu dobrze się pokrywają, gdy stałe ścieżki optyczne zbierały światło z tej samej odległości przed powierzchnią ablacyjną: stabilne zużycie materiału oraz rozmiar reaktywnej warstwy granicznej wynoszący około 5 mm przed powierzchnią (warunek A1a). Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rycina 3. Temperatury pasma fioletowego CN oszacowane metodą dopasowania widmowego. Metoda najmniejszych kwadratów dla najlepszego dopasowania widm fioletowego CN, obliczonych za pomocą programu SPECAIR 2.2, dostarczyła temperatur translacyjno-obrotowych i wibracyjno-elektronowych Trot oraz Tvib: (a) Warunek A1a: Trot = 8 240 K ± 400 K, Tvib = 21 600 K ± 1 700 K, TLTE = 12 600 K ± 500 K (symulacja równowagowa TLTE podana dla porównania); (b) Warunek A1a: Trot = 6 880 K ± 200 K, Tvib = 7 120 K ± 180 K. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rysunek 4. Profile temperatury pasma fioletowego CN w warstwie przyściennej. Temperatury translacyjno-obrotowe oraz wibracyjno-elektronowe Trot i Tvib wyznaczone z przeskalowanych i dopasowanych widm pasma fioletowego CN, obliczonych za pomocą narzędzia do symulacji promieniowania na czterech odległościach od powierzchni ablacyjnej, wskazują na odchylenie od stanu równowagi termicznej w pobliżu ścianki przy niskim ciśnieniu 15 hPa (a), natomiast przy ciśnieniu 200 hPa równowaga ta występuje w całej warstwie przyściennej (b). Zakresy temperatur [K] oparto na niepewności pomiaru natężenia emisji spektrometru na poziomie 10%, co pozwala na teoretyczną zmienność widma w ramach tych granic podczas procedury dopasowania. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rysunek 5. Mikrografy elektronowe z mikroskopu skaningowego po ablacji powietrznej (a), w tym zdjęcie in-situ (b) oraz mikrografy po ablacji azotem (c). (a) Mikrografy po ablacji powietrznej wykonane na powierzchni czołowej w pobliżu punktu stagnacji wykazują cieniowanie włókien węglowych spowodowane utlenianiem od końca włókna, prowadzące do kształtu nacieku, głębokość dyfuzji tlenu wynosi około 200 µm (ablacja ograniczona dyfuzją); (b) Zdjęcie wykonane podczas testu ablacji cylindrycznej próbki (czas naświetlania: 1/200 s) ilustruje intensywne iskrzenie; (c) Silna korozja została zaobserwowana w azocie na całej długości włókna. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.