Large eddy simulation of turbulent swirl-stabilized flames using the front propagation formulation: impact of the resolved flame thickness

Dit onderzoek breidt het frontpropagatiemodel uit voor LES-simulaties van turbulente draalgestabiliseerde vlammen en toont aan dat een nauwkeurige modellering van de opgeloste vlamdikte essentieel is om de vorming van ingevangen vlamzakken en secundaire temperatuurpieken correct te voorspellen.

De kern

Stel je voor dat je een enorme, chaotische danszaal hebt. In het midden draait een groep mensen (de vlam) razendsnel rond, terwijl er overal om hen heen andere mensen (de lucht) door de zaal stormen. De uitdaging voor wetenschappers is om te voorspellen hoe deze dans eruitziet, hoe heet het wordt en of de dansers niet per ongeluk uit elkaar vallen.

Dit papier gaat over het simuleren van zo'n danszaal: een turbo-brander in een vliegtuigmotor of gasturbine. De brander gebruikt een draaiende luchtstroom (swirl) om de vlam stabiel te houden. Maar omdat de dans zo snel en chaotisch is, is het bijna onmogelijk om elke beweging op de computer te tekenen. De computers zijn niet snel genoeg.

Hier is wat de onderzoekers hebben gedaan, vertaald in alledaags taal:

1. Het Probleem: De "Wazige" Foto

Stel je voor dat je een foto maakt van een snel bewegende vlam, maar je camera is niet scherp genoeg. In plaats van een scherpe lijn zie je een dikke, wazige vlek.
In de computerwereld noemen we dit een opgelost rooster. Omdat de computer niet oneindig veel details kan verwerken, wordt de vlam op de computer vaak "te dik" weergegeven.

• Het gevaar: Als de vlam te dik is, denkt de computer dat de vlam zich op een verkeerde manier verspreidt. Het is alsof je denkt dat een vlam zich als een dikke, trage wolk verspreidt, terwijl hij in werkelijkheid als een snelle, scherpe lijn door de lucht snijdt. Dit leidt tot fouten: de computer denkt dat de vlam minder brandt dan hij eigenlijk doet.

2. De Oplossing: De "Slimme Scharnier" (FPF)

De onderzoekers hebben een nieuwe methode gebruikt, genaamd Front Propagation Formulation (FPF).

• De analogie: Stel je voor dat je de vlam niet tekent als een dikke muur, maar als een slimme, flexibele scharnierlijn. Zelfs als de computer "wazig" is, zorgt deze methode ervoor dat de lijn zich gedraagt alsof hij scherp is.
• Ze hebben deze methode aangepast voor twee extra dingen:
  1. Hitteverlies: In de echte wereld koelt de vlam iets af door de koude wanden (net als een warme soep die afkoelt in een metalen kom). De computer moest dit ook kunnen rekenen.
  2. De "Dikke Vlam" Correctie: Ze hebben een trucje bedacht om te voorkomen dat de vlam op de computer onnodig dik wordt. Ze zorgen ervoor dat de chemische reactie de vlam weer "dichtknijpt", net zoals een elastiekje dat terugveert.

3. De Experimenten: De "TECFLAM" Danszaal

Ze hebben hun methode getest op een beroemde brander (de TECFLAM), die al door veel anderen is bestudeerd.

• Ze draaiden de computer aan en keken of de resultaten overeenkwamen met echte metingen in het lab.
• Het resultaat: Het klopte! De computer voorspelde precies waar de vlam zat, hoe snel de lucht bewoog en hoe heet het was. Zelfs de kleine, lastige details kwamen overeen.

4. De Grote Ontdekking: De "Vlamzakjes"

Dit is het meest interessante deel. De onderzoekers ontdekten iets moois in de buitenste rand van de draaiende luchtstroom:

• Wat er gebeurt: Grote, draaiende wervels (zoals kleine tornadoatjes) in de buitenste laag van de brander grijpen de vlam vast. Ze rekken de vlam uit en rollen hem op, net als een deken die wordt opgerold. Hierdoor ontstaan er losse zakjes met hete lucht die vastzitten in de draaiende lucht.
• Waarom dit belangrijk is: Deze zakjes zorgen voor een tweede piek in de temperatuur. Als je dit niet ziet, mis je een belangrijk stukje van het verhaal.

5. De Les: Dikte Matters!

Hier komt de kern van hun ontdekking:

• Ze hebben twee simulaties gedaan: één met de "slimme scharnier" (waarbij de vlam scherp blijft) en één waarbij ze de scherpte hebben uitgeschakeld (de vlam werd dik).
• Bij de scherpe vlam: De computer zag de wervels de vlam grijpen, uitrekken en de zakjes vormen. Dit klopte met de werkelijkheid.
• Bij de dikke vlam: De vlam was te dik en "zweefde" te hoog. De wervels konden de vlam niet meer grijpen. De vlam raakte los van de buitenste luchtstroom. Resultaat: De computer zag de zakjes niet en voorspelde geen tweede temperatuurpiek.

Conclusie in één zin:
Om te voorspellen hoe een vlam in een turbulente motor zich gedraagt, moet je de computer niet toestaan dat de vlam "dik en wazig" wordt; je moet zorgen dat hij scherp blijft, anders mis je de spannende details (zoals de losse zakjes hete lucht) die bepalen hoe goed de motor werkt.

Kortom: Scherp is beter dan dik, zelfs in de wereld van computermodellen voor vlammen.

Technische samenvatting

Probleemstelling

Gestabiliseerde voorverbranding met werveling (swirl) wordt veel gebruikt in moderne gasturbines en vliegtuigmotoren om hoge efficiëntie en lage NOx-emissies te bereiken. Het simuleren van deze processen met Large Eddy Simulation (LES) vormt echter een grote uitdaging vanwege de complexe interacties tussen turbulentie en de vlam.

Een specifiek probleem in de praktijk is dat de rekengrids vaak vergelijkbaar met of groter zijn dan de laminare vlamdikte. Dit leidt tot een ondertekende vlamoplossing (under-resolved flame), waarbij de gefilterde reactiezone niet correct wordt opgelost. Dit resulteert in:

1. Onnauwkeurige voorspellingen van het totale reactietempo.
2. Spurious (kunstmatige) voortplanting van de vlamfronten.
3. Een discrepantie tussen de werkelijke opgeloste vlamdikte en de gebruikte filtergrootte, vooral als chemische verdunningseffecten (steepening) niet correct worden gemodelleerd.

Bestaande methoden, zoals de Thickened Flame LES (TFLES) of G-equation modellen, lossen dit deels op, maar de toepassing van de Front Propagation Formulation (FPF) op complexe, swirl-gestabiliseerde verbranding met brandstofstratificatie en warmteverlies was nog niet onderzocht.

Methodologie

De auteurs hebben de FPF-methode uitgebreid en toegepast op de LES van turbulente, voorverbrande swirl-vlammen in de TECFLAM-brander. De kernpunten van de methodologie zijn:

1. Uitgebreide FPF-methode:

   • De chemische reactiebronterm wordt gemodelleerd als het product van de onbrandbare dichtheid, de sub-filter vlamconsumptiesnelheid en een geregulariseerde Dirac-deltafunctie.
   • Deze deltafunctie ($\psi$) is ontworpen om de chemische verdunningseffecten (steepening) van de vlam te reproduceren, waardoor de voortplantingsfouten bij onopgeloste vlammen worden geminimaliseerd.
   • De vorm van $\psi$ wordt bepaald door modelparameters ($\alpha$ en $\gamma$) die afhankelijk zijn van de verhouding tussen filtergrootte en laminare vlamdikte.

2. Aanpassingen voor realistische omstandigheden:

   • Niet-adiabatische effecten: Er is een bibliotheek voor laminare vlamsnelheid en -dikte ontwikkeld op basis van 1D-simulaties met warmteverlies, gekoppeld aan de totale enthalpie.
   • Verbeterde sub-filter vlamsnelheid: De schatting van de sub-filter vlamverwringing (wrinkling) is aangepast om inconsistenties op te lossen die ontstaan door warmteontwikkelingseffecten op lokale turbulentie binnen de vlamzone. In plaats van turbulentieparameters van de onbrandbare stroming te gebruiken, worden deze geëxtrapoleerd vanuit de onbrandbare zijde naar de vlamzone.

3. Numeriek Opzet:

   • De simulaties zijn uitgevoerd voor de TECFLAM-brander (Re = 10.000, swirl-getal 0.75, $\phi$ = 0.833).
   • Er zijn twee scenario's vergeleken om het effect van de opgeloste vlamdikte te onderzoeken:
     • Geval A (Correct): $\gamma = 2.5$ (chemische verdunningseffecten zijn actief).
     • Geval B (Onjuist): $\gamma = 1.0$ (chemische verdunningseffecten zijn afwezig, wat leidt tot een overgeschatte vlamdikte).

Belangrijkste Bijdragen

• Eerste toepassing van FPF op complexe swirl-verbranding: De methode is succesvol uitgebreid naar configuraties met brandstofstratificatie, warmteverlies en sterke wervel-vlam interacties.
• Oplossing voor turbulentie-consistentie: Een nieuwe aanpak voor het schatten van sub-filter verwringing die rekening houdt met de verandering in turbulentie door warmteontwikkeling binnen de vlam.
• Inzicht in de rol van vlamdikte: Het artikel demonstreert kwantitatief hoe de nauwkeurigheid van de opgeloste vlamdikte direct invloed heeft op de voorspelling van complexe vlamdynamica, zoals het vormen van "flame pockets".

Resultaten

De resultaten tonen een uitstekende overeenkomst met experimentele metingen voor snelheid, mengfractie en temperatuur. De belangrijkste bevindingen zijn:

1. Validatie: De uitgebreide FPF-methode reproduceert de gemiddelde en fluctuerende snelheidsprofielen, inclusief het centrale recirculatiegebied en de shear-lagen, zeer nauwkeurig.
2. Invloed van $\gamma$ op vlamdikte:
   • Bij $\gamma = 2.5$ (met verdunningseffecten) blijft de opgeloste vlamdikte relatief dun en neemt de vlamoppervlakte toe door verwringing.
   • Bij $\gamma = 1.0$ (zonder verdunningseffecten) wordt de vlamdikte aanzienlijk overschat (ongeveer drie keer zo dik nabij de uitlaat). Hierdoor wordt de totale vlamoppervlakte onderschat.
3. Secundaire temperatuurpieken en Flame Pockets:
   • Experimenten tonen secundaire temperatuurpieken in de buitenste shear-laag. Deze worden alleen correct voorspeld in de simulatie met $\gamma = 2.5$.
   • Mechanisme: Grote wervelstructuren in de buitenste shear-laag rekken het vlamfront uit, wat leidt tot het vormen en vastlopen van geïsoleerde vlamzakken ("flame pockets" of "peninsula's").
   • Bij $\gamma = 1.0$ is de vlam te dik en "ontkoppelt" deze zich van de buitenste shear-laag. Hierdoor worden deze wervel-vlam interacties niet vastgelegd, worden er geen vlamzakken gevormd, en verdwijnt de secundaire temperatuurpiek.

Betekenis en Conclusie

De studie benadrukt dat in LES van turbulente swirl-vlammen de opgeloste vlamdikte een kritieke parameter is. Als chemische verdunningseffecten niet correct worden gemodelleerd (wat leidt tot een te dikke vlam), wordt het verbruik van brandstof onderschat en de vlamhoek te klein. Dit veroorzaakt een fysiek onjuiste ontkoppeling van de vlam van de shear-laag, waardoor complexe fenomenen zoals het vastlopen van vlamzakken en de bijbehorende temperatuurprofielen niet kunnen worden voorspeld.

De conclusie is dat voor betrouwbare LES van complexe verbrandingssystemen niet alleen de filtergrootte, maar vooral de nauwkeurige weergave van de opgeloste vlamdikte (via correcte modellering van chemische steepening) essentieel is om de juiste vlamdynamica en warmteoverdracht te simuleren. De uitgebreide FPF-methode biedt hiervoor een robuust kader.