A Novel Approach for Direct Measurement of the Stretch Factor in Laminar Premixed Hydrogen-Air Flames Affected by Thermodiffusive Instabilities

Deze studie introduceert een nieuwe experimentele methode met OH-PLIF-imaging om de rekfactor ($I_0$) in laminaire waterstof-luchtvlammen direct te meten, waarbij een overgang van een quasi-stabiele naar een thermodiffusief instabiele toestand wordt benut om de toename in verbrandingssnelheid en oppervlakte te kwantificeren.

De kern

De Brandende V-vorm: Een Reis door de Onrust van Waterstofvlammen

Stel je voor dat je een waterstofvlam laat branden in een stroom lucht. Normaal gesproken zou je denken aan een rustig, glad vlammend oppervlak, zoals een kalme rivier. Maar bij waterstof, vooral als er weinig brandstof in zit (een "mager" mengsel), gebeurt er iets magisch en chaotisch: de vlam begint te rillen, te krullen en cellen te vormen, alsof het water kookt of een onrustige menigte.

Dit fenomeen noemen wetenschappers thermodiffusieve instabiliteit. Het klinkt ingewikkeld, maar het is eigenlijk heel simpel: waterstofdeeltjes bewegen sneller dan warmte. Hierdoor ontstaan er kleine "hotspots" die de vlam sneller laten branden, waardoor het oppervlak van de vlam uit elkaar trekt en rimpels krijgt.

Het Probleem: Hoe meet je de "kracht" van die onrust?
De onderzoekers van deze studie wilden weten: Hoeveel sneller brandt deze vlam precies door die rimpels? Ze noemen dit de rekfactor (stretch factor). Het is alsof je wilt weten hoeveel extra snelheid een auto krijgt door een extra turbo, maar dan voor een vlam.

Vroeger was dit lastig te meten. Je moest vaak complexe 3D-simulaties draaien op supercomputers, wat duur en tijdrovend is.

De Oplossing: Een V-vormige Vlam als Meetinstrument
De onderzoekers hebben een slimme, nieuwe manier bedacht. Ze gebruiken een heel dun keramisch staafje om een V-vormige vlam te stabiliseren.

Stel je die vlam voor als een V-vormige weg:

1. De bovenkant (Stabiel): Dicht bij het staafje is de vlam glad en rustig. Dit is de "stille zone".
2. De onderkant (Onstabiel): Na een stukje begint de vlam plotseling te rimpelen en te krullen. Dit is de "onrustige zone".

Het geniale aan hun methode is dat ze deze twee zones direct naast elkaar kunnen vergelijken.

Hoe werkt de meting? (De Creatieve Analogie)
Stel je voor dat je twee fietsers hebt die dezelfde weg afleggen:

• Fietsster A rijdt op een gladde, rechte weg (de stabiele vlam).
• Fietsster B rijdt op een weg vol kuilen en kronkels (de onstabiele vlam).

Als je kijkt naar de hoek die hun weg maakt met de horizon, zie je iets interessants: Fietsster B moet een steilere weg nemen om vooruit te komen. Dat betekent dat ze harder moet trappen (sneller brandt) om dezelfde afstand in dezelfde tijd te overbruggen.

De onderzoekers gebruiken een speciale camera (OH-PLIF) die de vlam als een foto vastlegt. Ze meten twee dingen:

1. De Hoek: Hoe steiler de vlam staat, hoe sneller hij brandt.
2. De Oppervlakte: Hoeveel "rimpels" er in de vlam zitten. Een gerimpelde vlam heeft meer oppervlak dan een gladde vlam, net zoals een gekreukeld stuk papier meer oppervlak heeft dan een glad vel.

Door de hoek en de oppervlakte te combineren, kunnen ze precies berekenen hoeveel extra "kracht" (snelheid) de vlam krijgt door die onrust.

Wat hebben ze ontdekt?

• Hoe minder brandstof, hoe meer chaos: Bij een heel mager mengsel (weinig waterstof) is de onrust het grootst. De vlam wordt extreem snel en rimpelig.
• Hoe meer brandstof, hoe rustiger: Als je iets meer waterstof toevoegt, kalmeert de vlam iets. De extra snelheid die de rimpels geven, wordt kleiner.
• De theorie klopt: Hun nieuwe meetmethode gaf resultaten die perfect overeenkwamen met wat de theorie voorspelde en wat de computersimulaties lieten zien.

Waarom is dit belangrijk?
Waterstof wordt steeds belangrijker als schone brandstof voor de toekomst. Maar omdat waterstofvlammen zo snel en onvoorspelbaar kunnen worden (ze kunnen bijvoorbeeld terugslaan in de brander), is het cruciaal om ze goed te begrijpen.

Deze nieuwe methode is als een nieuwe bril voor ingenieurs. Het stelt hen in staat om zonder dure supercomputers direct te zien hoe waterstofvlammen zich gedragen. Dit helpt bij het bouwen van veiligere en schonere motoren en turbines in de toekomst.

Kortom:
Ze hebben een slimme manier bedacht om de "onrust" van een waterstofvlam te meten door te kijken naar de hoek en de rimpels in een V-vormige vlam. Het bewijst dat hoe minder brandstof je hebt, hoe wilder de vlam wordt, en dat je dit nu heel nauwkeurig kunt meten.

Technische samenvatting

Probleemstelling

Waterstof (H2) is een cruciale energiedrager voor de decarbonisatie van de industrie. Echter, bij mager voorverwarmde verbranding (lean premixed) zijn vlammen met een Lewis-getal kleiner dan 1 (zoals waterstof-lucht vlammen) zeer gevoelig voor thermodiffusieve instabiliteiten (TDI). Deze instabiliteiten worden veroorzaakt door preferentiële en differentiële diffusie, wat leidt tot de vorming van cellulaire structuren en "vingerachtige" patronen in de vlamfront.

Deze TDI's hebben een significant effect op:

• De lokale vlamfrontsnelheid en het oppervlak van de vlam.
• Het ontstaan van superadiabatische temperaturen.
• Het risico op vlamterugslag (flashback) en vlam-wandinteracties.

Hoewel de fysica achter TDI's bekend is, is het experimenteel bepalen van de rekfactor (stretch factor, $I_0$) – een maat voor de versterking van de lokale vlamfrontsnelheid door diffusieve effecten – zeer uitdagend. Bestaande methoden (zoals straalbranders of bolvormig uitdijende vlammen) worden vaak beïnvloed door door de configuratie opgelegde kromming en interacties tussen vlamtakken, wat de meting van de intrinsieke instabiliteiten verstoort. Er is een gebrek aan betrouwbare experimentele data voor $I_0$ in laminair voorverwarmde waterstofvlammen.

Methodologie

De auteurs introduceren een nieuwe experimentele configuratie en analysemethode om $I_0$ direct te meten zonder directe meting van de brandstofverbruiksnelheid.

1. Experimentele Opstelling:

• Configuratie: Een V-vormige laminair voorverwarmde H2/lucht-vlam wordt gestabiliseerd op een keramische staaf (diameter 1 mm) in een stroming met een uniforme snelheid ($u_{bulk} = 2.28$ m/s).
• Diagnostiek: Er wordt gebruikgemaakt van OH-PLIF (Planar Laser-Induced Fluorescence van hydroxylradicalen) om het vlamfront te visualiseren.
• Vlamgedrag: De vlam vertoont twee distincte regimes:
  • Stabiele tak (upstream): Een glad, quasi-stabiel vlamfront met een kleine hellingshoek ($\theta_s$).
  • Instabiele tak (downstream): Een gebied waar TDI's optreden, gekenmerkt door sterke rimpeling (wrinkling) en cellulaire structuren, wat leidt tot een steilere hellingshoek ($\theta_u$).
• Variatie: Metingen worden uitgevoerd voor verschillende equivalentieratio's ($\phi$) tussen 0,35 en 0,40.

2. Data-analyse en Berekening van $I_0$:
De methode benut de abrupte overgang tussen de stabiele en instabiele tak:

• Snelheidsverhouding: Uit de geometrische verhouding van de hellingshoeken ($\theta_u$ en $\theta_s$) kan de verhouding van de verbrandingssnelheden worden afgeleid: $S_u/S_s = \sin(\theta_u) / \sin(\theta_s)$. Dit wordt geïnterpreteerd als de genormaliseerde vlamverbruiksnelheid ($S_c/S_L$).
• Oppervlaktevergroting: Drie verschillende beeldverwerkingsalgoritmen ("f-canny", "f-otsu", "f-otsuCanny") worden gebruikt om de toename van het vlamoppervlak door rimpeling te kwantificeren ($A/A_0$), waarbij $A$ het gerimpelde oppervlak is en $A_0$ het gladde referentieoppervlak.
• Berekening Rekfactor: De rekfactor wordt berekend als:
  $$I_0 = \frac{S_u}{S_s} \cdot \frac{A_0}{A}$$
• Numerieke Validatie: De experimentele bevindingen worden vergeleken met 2D-numerieke simulaties (met de solver CONVERGE) die dezelfde configuratie nabootsen. Hierbij wordt de volledige thermochemische staat gebruikt om het vlamfront (via een progressievariabele $Y_c$) te definiëren, wat onzekerheden in de experimentele OH-metingen omzeilt.

Belangrijkste Bijdragen

1. Nieuwe Meetmethode: Een directe methode om $I_0$ te bepalen op basis van OH-PLIF-beeldvorming van een staaf-geankerde V-vlam, waarbij de overgang van stabiel naar instabiel wordt gebruikt als natuurlijke "schakelaar".
2. Minimalisatie van Configuratie-effecten: In tegenstelling tot bolvormige vlammen of straalbranders, minimaliseert deze V-vlam-configuratie door de divergerende takken interacties tussen vlamfronten en opgelegde kromming, waardoor TDI's vrijer kunnen ontwikkelen.
3. Validatie door Simulatie: De combinatie van experimenten en 2D-simulaties bevestigt de kwalitatieve overgang en de geldigheid van de gebruikte analyse-methodiek.
4. Kwantificering van Onzekerheid: Het gebruik van drie verschillende beeldverwerkingsstrategieën voor het bepalen van het vlamoppervlak biedt inzicht in de robuustheid en de onzekerheden van de oppervlaktebepaling.

Resultaten

• Overgangspositie: De positie waar de TDI's beginnen (de overgang van glad naar gerimpeld) verschuift naar de brander toe naarmate de equivalentieratio ($\phi$) toeneemt. Dit correleert met de afname van de vlamdikte.
• Hellingshoeken: De hellingshoek van de instabiele tak ($\theta_u$) is altijd steiler dan die van de stabiele tak ($\theta_s$), wat aangeeft dat de effectieve vlamfrontsnelheid toeneemt door TDI's.
• Snelheids- en Oppervlakteverhoudingen:
  • De verhouding $S_u/S_s$ neemt monotoon af bij toenemende $\phi$ (van ca. 1,5 bij $\phi=0.35$ naar lager bij $\phi=0.4$).
  • De oppervlakteverhouding $A/A_0$ blijft over het onderzochte bereik vrijwel constant (ongeveer 1,2 - 1,3), hoewel numerieke data een lichte toename suggereert.
• Rekfactor ($I_0$):
  • De berekende $I_0$-waarden nemen monotoon af bij toenemende equivalentieratio.
  • Waarden variëren van ongeveer 1,1–1,3 bij $\phi = 0,35$ tot 0,8–0,9 bij $\phi = 0,4$.
  • Deze trend komt overeen met theoretische voorspellingen en eerdere numerieke modellen (zowel 2D als 3D).
  • De hogere waarden bij magere mengsels bevestigen dat TDI's sterker zijn bij lagere $\phi$ en dat de vlamfronten lokaal reactiever zijn door superadiabatische verbranding.
• Discrepanties: Er zijn kleine verschillen tussen experiment, simulatie en modellen, voornamelijk toe te schrijven aan de definitie van het vlamfront, drempelwaarden en het ontbreken van 3D-effecten in de 2D-simulaties en het meetvlak.

Betekenis en Conclusie

Dit onderzoek biedt een waardevolle basis voor de ontwikkeling van theorieën en de validatie van numerieke modellen voor vlammen die gevoelig zijn voor thermodiffusieve instabiliteiten. De voorgestelde methode maakt het mogelijk om de rekfactor $I_0$ direct en experimenteel te bepalen zonder complexe metingen van lokale stromingsvelden of brandstofverbruik.

De bevindingen bevestigen dat de toename van de globale vlamverbruiksnelheid in magere waterstofvlammen zowel wordt veroorzaakt door een toename van het vlamoppervlak als door een verhoging van de lokale vlamfrontsnelheid (door $I_0 > 1$). De studie benadrukt ook het belang van het rekening houden met 3D-effecten en lokale rek in toekomstig onderzoek, aangezien deze in de huidige 2D-benadering slechts gedeeltelijk worden vastgelegd. De methode is veelbelovend voor het verbeteren van ontwerpen voor veilige en efficiënte waterstofverbrandingssystemen.