Rayleigh-Taylor Unstable Flames: Thin and Thick

Dit artikel onderzoekt of Rayleigh-Taylor-instabiele vlammen met behulp van turbulente verbrandingstheorie kunnen worden beschreven, en concludeert dat hoewel deze vlammen door hun eigen gegenereerde turbulentie kunnen verdikken, hun structuur fundamenteel verschilt van die van klassieke turbulente vlammen, wat belangrijke implicaties heeft voor subgrid-modellering in toepassingen zoals Type Ia-supernova's.

De kern

Vuur in een storm: Waarom vlammen soms dikker worden en soms dunner

Stel je een vlam voor die niet rustig op een kaars staat, maar zich bevindt tussen twee lagen: een zware, koude laag (brandstof) en een lichte, warme laag (rook). Als de zware laag boven de lichte laag komt te liggen, gebeurt er iets raars: de zware laag wil naar beneden, de lichte naar boven. Ze proberen van plek te wisselen. Dit noemen wetenschappers de Rayleigh-Taylor-instabiliteit. Het is alsof je een olievat op water zet; de olie zakt en het water stijgt, maar dan in een chaotische dans van bellen en spitsen.

In dit artikel onderzoekt de auteur, E.P. Hicks, wat er gebeurt als er branding (verbranding) plaatsvindt op het grensvlak tussen deze twee lagen. Denk aan een vlam die probeert te branden terwijl de zware lucht erboven haar naar beneden duwt en de lichte lucht eronder haar omhoog trekt.

Hier is de kern van het verhaal, vertaald naar alledaagse taal:

1. Het mysterie: Dik of Dun?

In de wereld van de brandkunde weten we dat als je een vlam in een storm zet (turbulentie), de vlam vaak dikker wordt. De wind blaast de vlam uit, waardoor hij meer ruimte inneemt.

Maar bij deze specifieke vlammen (die door de zwaartekracht en de instabiliteit worden aangedreven) dachten wetenschappers jarenlang dat het tegenovergestelde gebeurde. Ze dachten: "De vlam wordt zo hard uitgerekt door de zwaartekracht dat hij juist dunner wordt, als een elastiekje dat uit elkaar wordt getrokken."

De auteur heeft echter gekeken naar een groot aantal simulaties en ontdekt dat het antwoord niet zo simpel is. Het hangt af van twee dingen:

1. De 'dikte' van de vlam zelf: Hoe snel de vlam brandt en hoe 'ruw' de chemische reactie is.
2. De 'stroperigheid' van de vloeistof: Hoe moeilijk het is voor de vloeistof om te bewegen (in de natuurkunde noemen we dit het Prandtl-getal).

Het verrassende resultaat: Bij bepaalde omstandigheden (vooral als de vloeistof heel 'stroperig' is en de vlam van nature al wat dikker is), wordt de vlam inderdaad dikker door de eigen turbulentie die hij zelf maakt. Maar! De manier waarop hij dik wordt, is heel anders dan bij een normale vlam in de wind.

2. De 'Zwevende Sandwich' (De structuur van de vlam)

Dit is het meest creatieve deel van het verhaal.

• Normale vlam in de wind: Stel je voor dat je een stukje papier in de wind houdt. De wind blaast de voorkant van het papier uit, waardoor het papier dikker en onrustig wordt. De vlam wordt dikker van voren.
• Deze speciale vlam: Deze vlam gedraagt zich als een zwevende sandwich met een dunne korst en een dikke vulling.
  • De voorkant (de kant die de brandstof raakt) wordt extreem dun uitgerekt door de zwaartekracht. Het is als een stukje zijde dat strak wordt getrokken.
  • De achterkant (de kant waar de rook uitkomt) wordt juist dik en rommelig door de turbulentie die de vlam zelf genereert. Het is alsof de achterkant van de vlam in een modderbad valt en uit elkaar begint te vallen.

De vlam wordt dus niet uniform dikker. Hij wordt een rare mix: een superdunne, strakke voorkant en een dikke, chaotische achterkant.

3. Waarom is dit belangrijk? (De 'Regelboeken' zijn verouderd)

Wetenschappers gebruiken vaak 'regels' (diagrammen) om te voorspellen hoe vlammen zich gedragen in complexe situaties, zoals in vliegtuigmotoren of bij explosies in sterren. Deze regels zijn gebaseerd op hoe vlammen reageren op wind.

De auteur laat zien dat deze oude regels niet werken voor deze vlammen.

• Als je de oude regels gebruikt, denk je dat de vlam zich gedraagt als een normale vlam in de wind.
• Maar in werkelijkheid is het een heel ander dier. Als je deze oude regels gebruikt om vliegtuigmotoren of sterrenexplosies te simuleren, krijg je de verkeerde uitkomsten. Je zou bijvoorbeeld denken dat de vlam sneller brandt, terwijl hij in werkelijkheid juist vertraagt omdat hij dikker wordt.

4. De nieuwe kaart (Het regime-diagram)

Omdat de oude regels niet werken, heeft de auteur een nieuwe kaart getekend.
Stel je voor dat je een oude landkaart hebt waarop een berg staat, maar de berg is eigenlijk een vulkaan. Je kunt die kaart niet gebruiken om te voorspellen of er lava uitkomt.

De nieuwe kaart van Hicks houdt rekening met de 'stroperigheid' en de 'dikte' van de vlam. Hiermee kunnen ingenieurs nu beter voorspellen:

• Vliegtuigmotoren: Om brandstof efficiënter te verbranden.
• Koudemiddelen: Om te weten hoe gevaarlijk nieuwe, milieuvriendelijke gassen (zoals ammoniak) branden als ze lekken.
• Sterren (Supernova's): Om te begrijpen hoe sterren exploderen. In een sterrenexplosie (Type Ia supernova) brandt de vlam door de ster heen. Als de vlam dikker wordt, kan dit leiden tot een enorme explosie. De auteur suggereert dat deze vlammen misschien wel kunnen overgaan in een 'distributie-verbranding' (een soort super-dikke vlam), wat de explosie nog krachtiger maakt.

Samenvatting in één zin

Deze vlammen zijn geen gewone vlammen die door de wind worden uitgedreven; ze zijn een bizarre strijd tussen zwaartekracht (die ze dun trekt) en eigen turbulentie (die ze dik maakt), en we hebben een nieuwe manier nodig om ze te begrijpen, omdat de oude regels van de brandkunde hier niet werken.

De les voor de leek: Als je iets probeert te begrijpen dat zowel zwaartekracht als branden combineert, moet je niet aannemen dat het zich gedraagt zoals een gewone vlam. Soms is het dun van voren en dik van achteren, en dat verandert alles.

Technische samenvatting

Hier is een gedetailleerde technische samenvatting van het artikel "Rayleigh-Taylor Unstable Flames: Thin and Thick" van E.P. Hicks, geschreven in het Nederlands.

Titel: Rayleigh-Taylor Instabiele Vlammen: Dun en Dik

Auteur: E.P. Hicks (CIERA, Northwestern University & Epsilon Delta Labs)
Doelwit: Journal of Fluid Mechanics

---

1. Het Probleem

Rayleigh-Taylor (RT) instabiele vlammen zijn dunne brandende interfaces die zich bevinden tussen een zware brandstoflaag en een lichte aschlaag. Deze fenomenen spelen een cruciale rol in diverse complexe systemen, waaronder:

• Aviatie: Ontwikkeling van lichte, efficiënte gasturbines (waarbij centrifugale krachten de RT-instabiliteit gebruiken om de verbrandingssnelheid te verhogen).
• Alternatieve Brandstoffen: Veiligheid van opslag en transport van brandbare mengsels zoals ammoniak en nieuwe koelmiddelen (R-32, R-1234yf), waarbij opwaartse stroming en RT-instabiliteit de vlamverspreiding beïnvloeden.
• Astrofysica: Type Ia supernova's, waarbij de overgang van een subsonische vlam (deflagratie) naar een supersonische vlam (detonatie) mogelijk wordt gestuurd door RT-instabiele vlammen in witte dwergen.

De Kernvraag:
Bestaande subgrid-modellen voor deze systemen baseren zich vaak op theorieën voor turbulente verbranding. Echter, RT-instabiele vlammen genereren hun eigen turbulentie (via barocliene vorticiteit) in plaats van dat ze worden blootgesteld aan externe turbulentie. Het is onduidelijk of de structuur van RT-vlammen kan worden beschreven met bestaande modellen voor turbulente verbranding, of dat ze fundamenteel anders gedragen.

2. Methodologie

De auteur voerde een uitgebreide parameterstudie uit met behulp van Directe Numerieke Simulaties (DNS) en bijna-DNS van de Boussinesq-vergelijkingen met een gemodelleerde reactieterm.

• Numerieke Opzet:
  • 3D-simulaties in een lange, vierkante schacht met periodieke zijwanden.
  • De vlam plant zich omhoog tegen de zwaartekracht in.
  • Gebruik van de open-source code Nek5000 (spectrale elementenmethode).
• Controleparameters:
  Om de lengte van de turbulente cascade binnen de vlamfront te controleren, werden vier parameters gemanipuleerd:
  1. Zwaartekracht ($G$): Bepaalt de sterkte van de RT-instabiliteit.
  2. Domeinbreedte ($L$): Bepaalt de integraalschaal van de turbulentie.
  3. Prandtl-getal ($Pr$): Bepaalt de verhouding tussen viscositeit en thermische diffusiviteit (cruciaal voor Type Ia supernova's waar $Pr \approx 10^{-5}$, versus aardse toepassingen waar $Pr \approx 1$).
  4. Reactiemodel ($p$): Een nieuwe parameter die de "steilheid" van de vlam bepaalt.
• De $p$-vlammen:
  De auteur ontwikkelde een familie van oplossingen voor een een-staps reactie waarbij de laminare vlamsnelheid constant blijft ($s_0=1$), maar de fysieke vlamdikte varieert door de parameter $p$.
  • Hoge $p$: Dunne vlammen.
  • Lage $p$: Dikke vlammen.
    Dit maakte het mogelijk om systematisch het effect van de reactiedikte op de dynamica te bestuderen zonder de snelheid te veranderen.
• Metingen:
  Directe metingen van de interne vlamstructuur (isovolume-slices tussen temperatuurcontouren $T=0.1$ en $T=0.9$) om de vlamdikte en oppervlakte-evolutie te analyseren.

3. Belangrijkste Bijdragen

1. Ontdekking van "Dikke + Langzame" Vlammen: Voor het eerst wordt aangetoond dat RT-instabiele vlammen onder bepaalde omstandigheden (lage $p$ en lage $Pr$) dikker worden dan hun laminare tegenhangers en langzamer branden. Dit staat haaks op eerdere bevindingen dat RT-vlammen altijd dunner en sneller zijn.
2. Interne Structuur-analyse: Het is het eerste werk dat de interne structuur van RT-vlammen direct meet. Er wordt een uniek patroon ontdekt: de vlam is dun aan de voorkant (door RT-rek) maar dik aan de achterkant (door zelf gegenereerde turbulentie).
3. Nieuwe Regime-diagram: Bestaande diagrammen voor turbulente verbranding (zoals die van Peters of Skiba) blijken ongeschikt om RT-vlammen te classificeren. De auteur introduceert een nieuw regime-diagram gebaseerd op een herformuleerde Karlovitz-getal ($Ka_x$) dat rekening houdt met $Pr$ en de fysieke vlamdikte ($p$).
4. Kritische Evaluatie van Modellen: Het werk waarschuwt voor het gebruik van kunstmatig verdikte vlammodellen (zoals het populaire KPP-model) in simulaties, aangezien dit de fysica van de zelfgegenereerde turbulentie en de vlamversnelling kan verstoren.

4. Resultaten

• Overgang van Dun naar Dik:
  Er is een duidelijke overgang van "dun + snel" naar "dik + langzaam" bij lage waarden van $p$ (dikke laminare vlam) en lage $Pr$. Deze overgang wordt gedreven door de interactie met de zelfgegenereerde turbulentie.
• De Rol van Turbulentie:
  De verdikking wordt veroorzaakt door kleine wervels die de vlam sneller mengen dan deze verbrandt. Dit wordt bevestigd door het feit dat de data goed collapseert wanneer geplot tegen een aangepast Karlovitz-getal ($Ka_x$) dat $Pr$ en de vlamdikte incorporeert.
  • Formule: $Ka_x = (u'/s_0)^{3/2} (L/\delta_F)^{-1/2} (\delta_F/\delta)^{1/2} Pr^{-1/2}$.
• Unieke Interne Structuur:
  • Turbulente vlammen: Verdikken van voren (preheat-zone) terwijl de reactiezone laminair blijft.
  • RT-vlammen: De voorkant blijft dun (gerekt door RT-instabiliteit), maar de achterkant (asch-zijde) wordt dik en gedesynchroniseerd door turbulentie. De temperatuurcontouren aan de achterkant evolueren onafhankelijk van de voorkant.
• Gedistribueerde Verbranding (Distributed Burning):
  De studie concludeert dat een volledige overgang naar een toestand van "gedistribueerde verbranding" (waar turbulentie de verbranding volledig domineert en de RT-instabiliteit overwint) onwaarschijnlijk is. De RT-instabiliteit is nodig om de turbulentie te genereren; als de vlam te dik wordt, neemt de turbulentieproductie af en wordt de vlam weer dun. Er is dus een cyclisch evenwicht. Echter, bij zeer hoge $Ka_x$ en $G$ kan een vlam ontstaan met een zeer dunne voorkant en een zeer dikke, turbulente "borstel" achteraan, wat lijkt op gedistribueerde verbranding.
• Ongeschiktheid van Bestaande Diagrammen:
  Bestaande regime-diagrammen (theoretisch en experimenteel) falen omdat ze geen rekening houden met $Pr$ en omdat ze aannemen dat verdikking van voren komt. RT-vlammen vallen niet in de bestaande regio's zoals "Thin Reaction Zones" op een betekenisvolle manier.

5. Betekenis en Implicaties

• Voor Modellering:

  • Flamelet-modellen: Deze zijn mogelijk niet geschikt voor RT-vlammen omdat de interne structuur (dunne voorkant, dikke achterkant) fundamenteel verschilt van laminare vlammen, zelfs als de vlam gemiddeld dun is.
  • Subgrid-modellen: Simulaties die kunstmatig verdikte vlammen gebruiken (zoals het KPP-model) kunnen foutieve resultaten opleveren, vooral bij lage $Pr$ (zoals in supernova's). Ze kunnen vlammen te traag maken en de turbulentieproductie onderschatten.
  • Aanbeveling: Gebruik realistische reactiemodellen en pas correcties toe voor turbulentie-verdikking alleen als $Ka_x$ hoog is.

• Voor Type Ia Supernova's:
  De bevindingen suggereren dat naarmate de vlam naar de buitenste lagen van de ster reist (waar de dichtheid daalt en de vlam dikker wordt), de vlam dikker kan worden door zelfgegenereerde turbulentie. Dit zou de voorwaarden kunnen scheppen voor een Deflagratie-naar-Detonatie Overgang (DDT) via mechanismen zoals het Zel'dovich-gradiëntmechanisme of turbulente zelfversnelling, zelfs zonder vooraf bestaande turbulentie.

• Voor Aardse Toepassingen:
  Voor de veiligheid van alternatieve brandstoffen en koelmiddelen is het essentieel om te beseffen dat RT-instabiliteit de vlamversnelling kan verhogen, maar dat de vlamstructuur anders is dan bij standaard turbulente verbranding. Dit vereist specifieke subgrid-modellen voor RANS/LES-simulaties.

Conclusie:
RT-instabiele vlammen worden bepaald door een drievoudige competitie tussen de RT-instabiliteit (die de vlam verdunt), zelfgegenereerde turbulentie (die de vlam verdikt) en verbranding. Ze gedragen zich niet als standaard turbulente vlammen en vereisen nieuwe theoretische kaders en regime-diagrammen om correct gemodelleerd te worden.