Investigation of Differential Diffusion and Strain Coupling in Large Eddy Simulations of Hydrogen-Air Flames

Deze studie toont aan dat Large Eddy Simulations met een vlamlet-gebaseerde thermochemie, zonder gebruik van een database voor vervormde vlamletten, de complexe interacties tussen differentieel diffusie, rek en kromming in turbulente waterstof-luchtvlammen succesvol kan voorspellen, wat leidt tot betere overeenstemming met experimentele data en modelvereenvoudiging.

De kern

Waar gaat dit over? De "Hydrogen-Flits" en de Kunst van het Branden

Stel je voor dat je een heel efficiënte, schone motor wilt bouwen die op waterstof werkt. Waterstof is fantastisch omdat het geen CO2 uitstoot, maar het is ook een onruststoker. Het brandt heel snel, is heel licht en gedraagt zich anders dan benzine of aardgas.

De onderzoekers van deze studie (uit Italië, Nederland en het VK) wilden begrijpen hoe je waterstof-vlammen het beste kunt simuleren op een computer. Ze keken specifiek naar een vlam die wordt vastgehouden door een "bluff-body" (een kegelvormig obstakel in de stroom), vergelijkbaar met hoe een vlam in een brander wordt vastgehouden.

Het Grote Probleem: De "Snelvoetige" Waterstof

In een normale vlam (bijvoorbeeld op aardgas) bewegen alle deeltjes (zuurstof en brandstof) ongeveer even snel. Maar bij waterstof is dat anders.

• De Analogie: Stel je een renwedstrijd voor. De waterstof-atomen zijn de sprinters die razendsnel vooruit rennen. De zuurstof-atomen zijn de wandelaars die langzamer gaan.
• Het Effect: Omdat de waterstof sneller is, "loopt hij vooruit" op de zuurstof. Dit noemen ze differentiële diffusie. Hierdoor wordt de mix van brandstof en lucht lokaal verrijkt of verarmd, net als een groepje vrienden die uit elkaar lopen tijdens een wandeling.

De onderzoekers wilden weten: Kunnen we dit simuleren zonder dat we een supercomputer nodig hebben die de hele wereld in detail berekent?

De Oplossing: Een Slimme "Receptenboek"-Benadering

Om vlammen te simuleren, gebruiken wetenschappers vaak een "flamelet"-model.

• De Analogie: Denk aan een kookboek. In plaats van elke keer opnieuw te berekenen hoe een pan met soep kookt, kijken ze naar een standaardrecept (een "flamelet") voor een perfecte vlam. Ze gebruiken dit recept om te voorspellen wat er in de grote pan (de motor) gebeurt.

Het probleem is dat de meeste recepten in dit kookboek zijn gemaakt voor vlammen die niet worden uitgerekt of gekromd. Maar in een echte motor wordt de vlam wel uitgerekt (door wind/stroming) en gekromd (door turbulentie).

De onderzoekers gebruikten een nieuwere, slimmere versie van dit receptenboek. Ze voegden een correctie toe voor die "snelvoetige" waterstof-atomen. Ze wilden testen of hun model het gedrag van de vlam goed kon voorspellen, zelfs als de vlam wordt uitgerekt door de stroming.

Wat Vonden Ze? (De Resultaten)

1. Het Model Werkt:
   Ze vergeleken hun computer-simulatie met echte experimenten (met camera's die de vlam en de snelheid van de lucht zagen). Het nieuwe model gaf een heel nauwkeurig beeld van hoe de vlam eruitzag en hoe snel de lucht bewoog.

2. De "Sprinters" Veranderen Alles:
   Door rekening te houden met het feit dat waterstof sneller is dan zuurstof, zag het model iets belangrijks:

   • De vlam werd korter en stabieler.
   • De reactie (het branden) vond dichter bij het obstakel plaats.
   • Dit komt omdat de "sprinters" (waterstof) de brandstof lokaal verrijken waar de stroming snel is. Dit zorgt voor een hevigere reactie, waardoor de vlam zich dichter bij de bron vastzet.

3. De "Grote Stroom" is Belangrijker dan de "Kleine Bochten":
   Vaak denken mensen dat de kromming van de vlam (de bochten) het belangrijkste is. Maar in dit specifieke geval (met de kegel) bleek dat de strekkracht (de stroming die de vlam uitrekt) de hoofdrolspeler was.

   • De Les: Zelfs als je alleen simuleert met "rechte" recepten (zonder rek), kun je het gedrag van de vlam goed voorspellen, zolang je maar goed berekent hoe de stroming de vlam uitrekt. Je hebt dus geen duizenden complexe, vooraf berekende "recepten voor uitgerekte vlammen" nodig. Dat bespaart enorm veel rekenkracht!

Waarom is dit Belangrijk?

Dit onderzoek is een grote stap voorwaarts voor de toekomst van schone energie.

• Efficiëntie: Het betekent dat ingenieurs minder rekenkracht nodig hebben om nieuwe, veilige waterstof-branders te ontwerpen.
• Veiligheid: Omdat waterstof snel kan "terugslaan" (flashback) in leidingen, helpt dit model om vlammen beter te begrijpen en te beheersen.
• Klimaat: Als we waterstof-branders beter kunnen ontwerpen, kunnen we meer schone energie gebruiken zonder schadelijke stikstofoxiden (NOx) te produceren.

Samenvatting in één zin:

De onderzoekers hebben bewezen dat je met een slimme, simpele computer-simulatie (die rekening houdt met de snelle waterstof-deeltjes) precies kunt voorspellen hoe een waterstof-vlam zich gedraagt in een sterke windstroom, zonder dat je een onmogelijk complexe berekening nodig hebt.

Technische samenvatting

Probleemstelling

Waterstofverbranding biedt grote kansen voor de decarbonisatie, maar brengt unieke uitdagingen met zich mee vanwege de fysische en chemische eigenschappen van waterstof, zoals hoge moleculaire diffusiviteit, lage activeringsenergie en een hoge verbrandingssnelheid. Een cruciaal fenomeen bij waterstofvlammen is differentiële diffusie (preferentiële diffusie), waarbij lichtere moleculen (zoals waterstof) sneller diffunderen dan zwaardere moleculen. Dit leidt tot lokale verrijking of verarming van het mengsel, wat de vlamdynamica sterk beïnvloedt.

In Large Eddy Simulaties (LES) worden vaak vlamlet-gebaseerde modellen gebruikt. Traditionele vlamlet-modellen gaan echter uit van eenheidse Lewis-getallen (gelijke diffusie voor alle soorten), wat de effecten van differentiële diffusie en de koppeling met rek (strain) en kromming negeert. Voor waterstofvlammen is deze koppeling echter essentieel, zelfs bij matige rek-niveaus. De centrale vraag is of een vlamlet-database gebaseerd op ongetrokken (unstretched) vlammen voldoende is om de complexe interactie tussen differentiële diffusie, tangentiële rek en kromming in turbulente waterstofvlammen correct te voorspellen, zonder de noodzaak van een uitgebreide database van getrokken vlammen.

Methodologie

De auteurs hebben een LES-simulatie uitgevoerd van een voorverwarmde, voorverbrande waterstof-luchtvlam die wordt gestabiliseerd door een bluff-body (een conische obstakel). De simulatie is gebaseerd op het ongeconfineerde testgeval van de NTNU (Noorse Universiteit voor Wetenschap en Technologie).

• Model: Er is gebruik gemaakt van een vlamlet-gebaseerde thermochemie gekoppeld aan een presumed Filtered Density Function (FDF) sluiting.
• Differentiële Diffusie: Een recent ontwikkelde methode (Ferrante et al., 2024b) is toegepast om differentiële diffusie-effecten te modelleren. Dit gebeurt door correcties toe te passen op de transportvergelijkingen voor de mengfractie ($z$) en de vooruitgangvariabele ($c$), gebaseerd op een lineaire combinatie van species-massafracties en hun Lewis-getallen.
• Vlamlet-database: De database bestaat uit ongetrokken (unstretched) 1D laminair voorverbrande vlammen bij verschillende equivalentieratio's. Er is geen database gebruikt met expliciete rek-effecten.
• Vergelijking: Twee simulaties zijn uitgevoerd:
  1. Een case met differentiële diffusie ($Le_k \neq 1$).
  2. Een case met eenheidse Lewis-getallen ($Le_k = 1$).
• Validatie: De resultaten zijn gevalideerd tegen experimentele data (PIV voor snelheid en OH*-chemiluminescentie voor vlamstructuur) van het NTNU-experiment.

Belangrijkste Bijdragen

1. Validatie van het model: Het model is succesvol gevalideerd tegen experimentele data voor zowel het snelheidsveld als de vlamstructuur.
2. Koppeling Rek en Diffusie: Het onderzoek toont aan dat een model gebaseerd op ongetrokken vlammen, wanneer gecombineerd met een correctie voor differentiële diffusie, in staat is om de macroscopische effecten van rek en kromming op de mengfractieverdeling en reactiesnelheid te vangen.
3. Dominantie van Rek: In deze specifieke bluff-body configuratie wordt aangetoond dat tangentiële rek de dominante factor is voor veranderingen in de mengfractie, meer dan krommingseffecten (in tegenstelling tot eerdere studies met slot-branders waar thermische-diffusieve instabiliteiten een grotere rol speelden).
4. Vereenvoudiging van Modellering: Het bewijs dat een database van ongetrokken vlammen toereikend is voor het modelleren van waterstofvlammen met hoge rek, zolang de rek grotendeels wordt opgelost (resolved) in de LES. Dit biedt aanzienlijke vereenvoudigingen voor de modellering van turbulente waterstofverbranding.

Resultaten

• Vlamstructuur: De simulatie met differentiële diffusie ($Le_k \neq 1$) voorspelde een kortere vlam en een verschuiving van het reactiepiekgebied dichter bij de bluff-body, wat beter overeenkwam met experimentele waarnemingen dan de case met eenheidse Lewis-getallen.
• Mengfractie en Temperatuur:
  • Het model voorspelde correct een verrijking van het mengsel (toename van de mengfractie) in gebieden met positieve tangentiële rek en een afname in gebieden met negatieve kromming.
  • Dit leidde tot superadiabatische temperaturen (pieken tot ca. 2000 K) in het verankeringpunt van de vlam, veroorzaakt door de lokale verrijking.
• Reactiesnelheid: Er werd een toename van de reactiesnelheid waargenomen in gebieden met hoge positieve rek. Dit gedrag, kenmerkend voor vlammen met een negatief Markstein-getal, werd correct weergegeven door het model, ondanks het gebruik van ongetrokken vlamlet-data.
• Interactie Rek/Kromming: De analyse van de waarschijnlijkheidsdichtheidsfuncties (PDF) toonde aan dat de rek in deze configuratie de thermische-diffusieve instabiliteiten (die sterke kromming veroorzaken) onderdrukt. De vlam "jaagt" dus het gebied met hoge rek na, wat leidt tot stabilisatie dichter bij de bluff-body.

Betekenis en Conclusie

De studie concludeert dat de macroscopische effecten van differentiële diffusie en rek op voorverbrande waterstofvlammen, zelfs bij significante rek-niveaus, succesvol kunnen worden voorspeld met een vlamlet-gebaseerde aanpak die gebruikmaakt van een database van ongetrokken vlammen.

Dit is een belangrijke doorbraak omdat het de noodzaak van complexe en rekenintensieve databases met getrokken vlammen (strained flamelets) wegneemt. Zolang de meeste rek en kromming in de LES worden opgelost en er een geschikte correctie wordt toegepast voor differentiële diffusie, kunnen deze modellen de juiste thermochemische toestanden reproduceren. Dit biedt waardevolle inzichten voor het ontwerp van nieuwe verbranders en vereenvoudigt de modellering van turbulente waterstofverbranding aanzienlijk.