Eksperymenty przeprowadzono zarówno w konfiguracji pionowej, jak i w unikalnym poziomym tunelu aerodynamicznym na University of Maryland, przedstawionym na Rysunku 1. Zamiast tradycyjnego tunelu o obiegu zamkniętym lub z systemem zasysania, instalacja na University of Maryland wykorzystuje dmuchawę o zmiennej prędkości do wytworzenia nadciśnienia w komorze plenumerowej o wymiarach 100 x 75 x 100 cm, która wymusza przepływ powietrza przez kanał na przeciwległym końcu. Taka konfiguracja umożliwia prowadzenie ciągłych eksperymentów spalania, ponieważ dym nie jest recyrkulowany, tunel aerodynamiczny nie ulega uszkodzeniu i nie jest oddziaływany przez ogień, a termopary mogą swobodnie przemieszczać się w sekcji próbkowania. Kanał wylotowy składa się z zbieżnej sekcji o długości 122 cm i szerokości 30,5 cm, połączonej z komorą plenumerową. Aby wyprostować przepływ i zmniejszyć intensywność turbulencji wlotowych, przy wejściu i wyjściu sekcji zbieżnej umieszczono gęste ekrany siatkowe, a 110 cm przed wylotem tunelu zainstalowano strukturę plastra miodu o grubości 5 cm z otworami o średnicy 0,3 cm. Prędkość przepływu opuszczającego tunel aerodynamiczny jest regulowana poprzez zmianę prędkości wentylatora za pomocą sterownika szerokości impulsu (PWM), a próbki paliwa umieszcza się przy wylocie tunelu, gdzie prędkości przepływu zostały zweryfikowane przy użyciu anemometru drutowego.
Próbki paliwa przy wylocie tunelu aerodynamicznego umieszczono na ogniwie obciążenia, które w sposób ciągły mierzy ubytek masy próbki w czasie. Aby uniknąć zakłóceń wywołanych przez wiatr w obrębie ogniwa obciążenia, próbka została uniesiona na arkuszu aluminium (30,5 x 61,0 cm x 1,5 mm grubości) za pomocą dwóch uchwytów w kształcie litery U i otoczona izolacyjną płytą z włókna ceramicznego o grubości 1,27 cm, aby zapewnić gładką powierzchnię wokół palącej się próbki. Górna powierzchnia płyty została pokryta czarną matową farbą wysokotemperaturową o emisyjności około 98%, aby zapewnić odpowiednie tło do wizualnej obserwacji płomienia oraz uszczelnić izolację, która zawiera również spoiwa organiczne. Ponieważ płyta izolacyjna stanowi względnie tępy korpus dla nadchodzącego przepływu, umieszczenie układu z próbką bezpośrednio w wylocie tunelu aerodynamicznego prowadziło do separacji przepływu i znacznych turbulencji obserwowanych w płomieniach. Wcześniejsze prace Ha et al. wykazały, że przymocowanie płyty przedłużającej do przedniej części próbki paliwa zapobiega tej separacji przepływu i zapewnia laminarny profil przepływu dopływającego do próbki. W związku z tym od przedniej krawędzi próbki do wylotu tunelu aerodynamicznego zamontowano cienką, metalową wargę o szerokości 10 cm i długości 40,6 cm, co pozwoliło uzyskać laminarny płomień dyfuzyjny, który ostatecznie okazał się zgodny z istniejącą teorią7.
Do testowania paliw ciekłych potrzebny był porowaty, niepalny knot. Do eksperymentów z przepływem wymuszonym wybrano arkusz wełny krzemianowej metali alkalicznych o wymiarach 10 cm x 10 cm x 1,27 cm ze względu na wysoką porowatość i niską przewodność cieplną. Aby zapobiec wyciekowi paliwa z próbki, zastosowano klej krzemian sodu do przymocowania folii aluminiowej do wszystkich ścian poza powierzchnią przednią. Próbkę poddano również „wypalaniu” w celu usunięcia spoiw organicznych, przesuwając palnikiem nad próbką przez około 20 minut, aż do momentu zmiany barwy płomienia z żółtej na niebieską (co wskazywało na usunięcie spoiw z próbki). Podczas testów knoty nasączano około 120 ml paliwa ciekłego (etanolu lub metanolu), co okazało się punktem nasycenia dla knotów o szerokości 10 cm.
Masową szybkość spalania paliwa określono poprzez pomiar utraty masy próbki w czasie podczas spalania przy częstotliwości 1 Hz. Układ próbki był umieszczony na precyzyjnej wadze o maksymalnym zakresie 32,2 kg i rozdzielczości 0,1 g, co było wystarczające do pomiaru tej szybkości ubytku masy z wysoką precyzją. Po zapaleniu próbki palnikiem gazowym szybkość ubytku masy skondensowanego paliwa wzrasta w funkcji czasu, aż do osiągnięcia stałej wartości, która następnie stopniowo spada pod koniec badania w miarę wypalania się paliwa. Ten obszar „stacjonarny”, w którym procesem spalania dominuje parowanie paliwa, a nie dyfuzja przez knot, jest obszarem zainteresowania, z którego pobierane są dane. W przypadku knota cieczowego stwierdzono, że próbki spalają się ze stałą szybkością ubytku masy przez około 400 s, co stanowi w przybliżeniu środkowe 80% czasu trwania badania. Wszystkie przedstawione szybkości spalania są średnimi z co najmniej sześciu powtórzonych testów w określonych warunkach, przy czym powtarzalność pomiarów znajdowała się w granicach 1,2% wartości średniej.
Do badań próbki paliwa stałego wybrano polimetakrylan metylu (PMMA), ponieważ spala się on stosunkowo stabilnie i nie ulega zwęglaniu. Aby zapalić próbkę, nad jej powierzchnią przesunięto palnik acetylenowy przez 50-60 sec, aż cała powierzchnia została równomiernie zapalona. Ze względu na niewielki rozmiar próbki paliwa oraz stwierdzoną wysoką powtarzalność wyników eksperymentalnych, metodę tę uznano za wystarczającą do inicjacji zapłonu. W przeciwieństwie do paliw ciekłych nasączonych w niepalnym knocie, paliwa stałe ulegają regresji w funkcji czasu, przez co nigdy nie osiągają w pełni ustalonego reżimu spalania. Zamiast tego wybrano pobieranie próbek z wczesnych etapów spalania, gdy paliwo pozostawało stosunkowo płaskie, co na podstawie pomiarów eksperymentalnych określono jako pierwsze 150 sec po zapłonie.
W przypadku paliw w stanie ciekłym i stałym rozkład temperatur nad powierzchnią paliwa mapowano w fazie gazowej przy użyciu termopar cienkodrutowych. Dla PMMA pomiary temperatury przeprowadzono w 6 punktach nad powierzchnią, zaczynając od warstwy stopionej w głąb fazy gazowej w odstępach 0,25 mm (dla testów konwekcji wymuszonej). W przypadku paliw ciekłych pomiary te wykonywano od cienkiej warstwy paliwa na powierzchni aż do 6 punktów z tą samą rozdzielczością. Profile te pobierano w 12 lokalizacjach wzdłuż długości powierzchni paliwa, w ciągu 400 s od zapłonu dla próbek ciekłych i w ciągu 150 s dla PMMA.
Wspomniane pomiary temperatury przeprowadzono przy użyciu mikrotermopar typu R Pt/Pt-13% Rh (zgrzewanych punktowo) o dwóch średnicach drutu, 50 μm (0.002 in) i 75 μm (0.003 in), z odpowiednio przybliżonymi średnicami złącz pomiarowych wynoszącymi 100 μm i 150 μm. Rozmiar termopar dobrano tak, aby były one jak najmniejsze, unikając jednocześnie częstego pękania (w celu zminimalizowania konieczności stosowania korekt radiacyjnych), niemniej niektóre korekty radiacyjne były nadal niezbędne. Zastosowano dwie termopary o różnych średnicach, aby lepiej określić odpowiednią korektę radiacyjną (opisaną poniżej). Mikrotermopary przesuwano za pomocą zestawu sterowanych komputerowo stolików X-Y unislides o maksymalnej rozdzielczości przestrzennej 1.5 μm. Sygnały napięciowe były następnie pozyskiwane, warunkowane i digitalizowane za pomocą modułu akwizycji danych o czułości pomiarowej do 0.02 °C. Oprogramowanie LabVIEW wykorzystano do synchronizacji ruchu obu termopar o średnicach drutu 50 μm i 75 μm z pomiarem temperatury w obrębie próbki.
Aby określić stosunkowo dokładną korektę promieniowania, dwa opisane rozmiary termopar przesuwano nad tym samym obszarem podczas powtórnych testów. W celu obliczenia strat ciepła z próbki zastosowano korelację Collisa i Williamsa5-6,8,
(1)
gdzie Nu to liczba Nusselta, a Re = Udw / v to liczba Reynoldsa, wyznaczona dla 0,02 < Re < 44, przy właściwościach ocenianych w temperaturze filmu, Τm, będącej średnią z temperatury gazu, Τg, oraz temperatury termopary, Τtc. Tutaj liczba Reynoldsa Re jest zdefiniowana zgodnie z powyższym dla lokalnej prędkości przepływu gazu U i lepkości kinematycznej v. dw w równaniu (1) oznacza średnicę drutu termopary.
W przypadku pomiarów w stanie ustalonym, takich jak opisany tutaj, bilans energii na złączu termopary sprowadza się do konwekcyjno-radiacyjnego bilansu ciepła (pomijając błędy wynikające z przewodnictwa i efektów katalitycznych), wyrażonego wzorem
(2)
(3)
gdzie Τg jest rzeczywistą temperaturą gazu, Τtc jest temperaturą złącza (lub kuleczki) termopary, Τsurr jest temperaturą otoczenia, εtc jest emisyjnością złącza termopary, σ jest stałą Stefana-Boltzmanna, a h jest współczynnikiem konwekcyjnego przekazywania ciepła przepływu nad złączem termopary, zdefiniowanym jako h = k Nu/d. k to przewodność cieplna gazu, Nu to liczba Nusselta, a d to średnica drutu termopary. Wybór korelacji liczby Nusselta ma kluczowe znaczenie przy obliczaniu korekty radiacyjnej do zmierzonej temperatury termopary, ponieważ, jak pokazano w równaniu (3), korekta radiacyjna jest odwrotnie proporcjonalna do liczby Nusselta. Wybór ten jest jednak utrudniony ze względu na istnienie wielu „odpowiednich” korelacji liczby Nusselta oraz trudność w oszacowaniu właściwości mieszaniny gazów otaczającej termoparę, w szczególności jej przewodności cieplnej. Większość dowodów w literaturze wyraźnie wskazuje jednak, że korelacja liczby Nusselta dla cylindra jest najbardziej odpowiednia do opisu konwekcyjnego przekazywania ciepła dla niemal wszystkich praktycznych termopar5-6, najlepiej korelacja Collisa i Williamsa8.
Korelację liczby Nusselta należy podstawić do stanu ustalonego bilansu konwekcyjno-promienistego (równanie 3), a pomijając niewielką zależność od temperatury, powstaje układ dwóch równań z dwiema niewiadomymi (a mianowicie Τg oraz U),
(4)
i
(5)
Równania (4) i (5) muszą być rozwiązywane iteracyjnie wspólnie w każdym punkcie, ponieważ przewodność elektryczna fazy gazowej oraz lepkość kinematyczna są funkcjami temperatury. Temperaturę kuleczki należy przyjąć jako pierwszą iterację temperatury gazu do wyznaczenia przewodności cieplnej i lepkości kinematycznej, a wartość iteracyjną należy aktualizować aż do osiągnięcia niskich błędów. Podczas rozwiązywania równań wynika, że poprawka na promieniowanie (tj. różnica między odczytem termopary a temperaturą rzeczywistą) zwiększa się dla termopar o większej średnicy i maleje wraz ze wzrostem prędkości przepływu nad kuleczką. dw1 oraz dw2 w równaniach (4) i (5) reprezentują średnice drutów termopary użytych w naszej badaniu.
Emisyjność kuleczki (εtc) może być również wyznaczona jako funkcja temperatury przy użyciu metody opisanej przez Jakoba9. W swojej analizie Jakob rozwiązuje równania falowe Maxwella dla zespolonych współczynników załamania na powierzchni metalicznej jako funkcję jej oporności elektrycznej. Przyjęto założenie dla granicy niskiej oporności i dużych współczynników załamania, co jest prawdziwe dla metali, co pozwala uzyskać prostą korelację dla całkowitej emisyjności półsferycznej platyny (Pt) w postaci:
(6)
gdzie dla platyny re ≈ re,273T / 273, przy czym T jest podane w K, a re,273 = 11x10-6 Ω–cm.
Zatem emisyjność platyny wynosi5-6
(7)
dla 0 < T < 2,330 K. Emisyjność kuleczki termopary lub złącza, występująca w równaniach (4) i (5), może zatem zostać obliczona przy użyciu powyższego wyrażenia. W przypadku równań (6) i (7) iteracja nie jest konieczna, ponieważ rzeczywista wartość temperatury kuleczki jest znana; jedynie temperaturę i prędkość gazu w równaniach (4) i (5) należy wyznaczyć iteracyjnie.
Podczas eksperymentów dwie termopary wprowadzono dokładnie do tych samych punktów pomiarowych, a dane pobrano w celu uwzględnienia korekty radiacyjnej w pomiarach temperatury. Korekty zastosowane w wyniku iteracji równań (4) i (5) były niewielkie, na przykład wyniosły jedynie +79 K dla termopary o średnicy drutu 50 μm przy 1 700 K i mniej niż 5 K w pobliżu powierzchni paliwa6. Ponieważ termopary przechodzą również przez obszary wysokich gradientów temperatury, należy rozważyć straty wynikające z przewodzenia ciepła przez drut; jednak ze względu na małe powierzchnie przekroju drutów termopar, błędy te obliczono na poziomie < 1%, w związku z czym żadne korekty nie były konieczne5-6.
Dzięki umieszczeniu powierzchni paliwa w centrum strumienia powietrza przy wylocie tunelu aerodynamicznego zapewniono łatwy dostęp do powierzchni paliwa w celu wykonania pomiarów za pomocą mikrotermopary i anemometru drutowego. Podczas prób zimnego przepływu w tunelu aerodynamicznym (bez spalania) prędkość strumienia swobodnego, U∞ tunelu aerodynamicznego była kalibrowana przy użyciu anemometru drutowego z częstotliwością próbkowania 50 000 próbek/s przez łączny czas 10 s na punkt. Wyznaczono profil prędkości wzdłuż wylotu całego tunelu, co wykazało występowanie spójnego przepływu tłokowego emanującego z centrum wylotu tunelu. Jest to spodziewane w przypadku kanału kwadratowego, takiego jak wylot naszego tunelu aerodynamicznego. Wcześniejsze pomiary przeprowadzone przez Sforza et al.10 wykazały, że długość rdzenia potencjalnego kwadratowego strumienia przy liczbie Reynoldsa Red w zakresie od 2,6 do 8,8 x 104 powinna wynosić około 5 d w dół strumienia od wylotu, gdzie d jest wysokością kanału. Dla d = 30,48 cm, czyli szerokości wylotu tunelu aerodynamicznego, Red mieści się w zakresie od 1,5 x 104 do 3,9 x 104, co oznacza, że próbka pozostaje w odległości 1d (20 cm) od wylotu tunelu. Powtarzalność tych pomiarów mieściła się w granicach 3% wartości średniej.
Temperatury mierzono na powierzchni zapalonej płyty z PMMA o wymiarach 10 cm x 10 cm x 1,27 cm, umieszczonej przy wylocie z tunelu aerodynamicznego pracującego przy U∞ = 0,79 m/sec i 2,06 m/sec. Opisane powyżej procedury wykorzystano do przeprowadzenia pomiarów temperatury, które poddano bezwymiarowaniu w odniesieniu do długości normalnej y*=y/L oraz temperatury, T* = (T - Tw,p / Tfl,ad - Tw,p), gdzie Τw,p oraz Τfl,ad reprezentują odpowiednio temperaturę ścianki i adiabatyczną temperaturę płomienia dla danego paliwa, y to pozycja normalna do powierzchni paliwa, w której mierzona jest temperatura, a L to długość powierzchni paliwa. Następnie obliczono bezwymiarowe gradienty temperatury normalne do powierzchni, (∂T* / ∂y*)y*=0, dopasowując wielomian piątego stopnia do bezwymiarowych temperatur i wyznaczając nachylenie na powierzchni paliwa, y*=0.
Rycina 2 (a) przedstawia te bezwymiarowe gradienty temperatury wzdłuż powierzchni paliwa. Są one wyraźnie najwyższe przy krawędzi natarcia powierzchni paliwa, gdzie płomień znajduje się najbliżej powierzchni paliwa, i maleją w kierunku krawędzi spływu (x = 100 mm), gdzie płomień jest najdalej od powierzchni paliwa. Bezwymiarowe gradienty temperatury można wykorzystać do wyznaczenia lokalnej szybkości spalania masy poprzez zastosowanie korelacji4,6,
(8)
gdzie B to liczba przenikania masy dla danego paliwa, kw przewodność cieplna powietrza wyznaczona w temperaturze ścianki, cp ciepło właściwe powietrza wyznaczone w adiabatycznej temperaturze płomienia paliwa, a L długość powierzchni pirolizującego paliwa. Stwierdzono wówczas, że lokalna szybkość spalania masy zmienia się w sposób podobny do bezwymiarowych gradientów temperatury, przedstawionych na Rysunku 2 (b).
W przeciwieństwie do paliw ciekłych, w przypadku PMMA lokalną szybkość spalania masy można również przybliżyć a posteriori poprzez pomiar lokalnego ubytku powierzchni w stałych odstępach czasu2,11. Próbki PMMA spalano w reprezentatywnych warunkach przez okresy czasu zaczynając od 50 sec i zwiększając je w interwałach 50 sec, po czym następowało wygaszenie próbki. Strumień masy pirolizy dla PMMA oblicza się dla każdego położenia x wzdłuż osi symetrii centralnej, korzystając z przybliżenia pierwszego rzędu podanego przez Pizzo et al.11, omówionego w innych publikacjach4-6. Średnią gęstość PMMA, ρs = 1,190 kg/m3, wykorzystano wraz z pomiarami ubytku powierzchni wzdłuż powierzchni paliwa do wyznaczenia szybkości utraty masy podczas każdego 50 sec interwału na całej długości próbki paliwa. Chociaż pożądany byłby krótszy krok czasowy, błędy pomiarowe sprawiają, że staje się on niepraktyczny, gdy kroki czasowe są krótsze niż 50 sec5.
Aby porównać lokalne szybkości ubytku masy mierzone za pomocą termopar z tymi uzyskanymi z profili regresji, wykorzystano dane z czasów wypalania paliwa wynoszących 100 i 150 s w celu porównania lokalnych szybkości spalania masy przedstawionych na Ryc. 2 (b). Czasy te odpowiadają w przybliżeniu momentom, w których dokonano pomiarów. Jak widać na rycinie, obie metody pomiaru lokalnej szybkości spalania masy dają bardzo zbliżone wyniki, co sugeruje, że zastosowana metodologia sprawdza się w przypadku tego typu płomieni.
W przypadku płomieni zdominowanych przez konwekcję, takich jak te małe płomienie laminarne, gradienty temperatury na powierzchni paliwa mogą być również wykorzystane do wyznaczenia konwekcyjnych strumieni ciepła, ponieważ są one w istocie bezpośrednio powiązane z gradientem temperatury na powierzchni. Wykorzystując pomiary szybkości ubytku masy, można również wyznaczyć składowe strumienia ciepła płomienia wzdłuż strefy pirolizy. Stosując kilka przybliżeń bilansu cieplnego na powierzchni paliwa, opisanych w literaturze2-3, składowe te można określić dla powierzchni płonącej płyty z PMMA. Rysunek 3 przedstawia ten wynik dla płomienia PMMA stabilizowanego przy prędkości swobodnego strumienia otoczenia U∞= 2,06 m/s. Technika ta może być zatem niezwykle użyteczna w ocenie wielu parametrów opisujących spalanie małych próbek paliw, co prowadzi do lepszego zrozumienia procesu spalania, a w szczególności zależności między fazą stałą a gazową.

Rysunek 1. Układ eksperymentalny. (a) Schemat układu eksperymentalnego wykorzystanego do pomiaru szybkości ubytku masy i profilu temperatury w płomieniu dyfuzyjnym warstwy przyściennej z wymuszonym konwekcją. (b) Układ eksperymentalny do badania płomieni dyfuzyjnych warstwy przyściennej w warunkach wymuszonego przepływu. Kliknij tutaj, aby zobaczyć powiększoną wersję tego rysunku.

Rysunek 2. Wyniki gradientu temperatury i lokalnej szybkości spalania. (a) Zmienność normalnych bezwymiarowych gradientów temperatury wzdłuż powierzchni paliwa dla płomienia dyfuzyjnego warstwy przyściennej PMMA przy U∞= 0.79 m/sec i 2.06 m/sec, odpowiednio. (b) Zmienność lokalnych szybkości spalania masy dla płomieni dyfuzyjnych warstwy przyściennej PMMA przy różnych warunkach strumienia swobodnego. Lokalne szybkości spalania masy uzyskane poprzez bezwymiarowe gradienty temperatury zestawiono z danymi eksperymentalnymi uzyskanymi poprzez regresję powierzchni PMMA. Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.

Rysunek 3. Wyniki strumienia ciepła w warunkach przepływu wymuszonego. Rozkład poszczególnych składowych strumienia ciepła płomienia w strefie pirolizy dla dyfuzyjnego płomienia warstwy przyściennej PMMA przy U∞= 2.06 m/sec. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.