Flamelet Connection to Turbulence Kinetic Energy Dissipation Rate

Dit artikel stelt voor om het turbulentie-kinetische-energie-dissipatievermogen ($\epsilon$) te gebruiken als koppelingsvariabele tussen opgeloste schalen in RANS- of LES-simulaties en een subgrid-vlamletmodel, waarbij de invloed van vorticiteit op de nauwkeurigheid van de berekende verbrandingsparameters wordt onderzocht.

De kern

De Vlam in de Turbulentie: Hoe een Simpel Getal de Branding van Vuur in Motoren Kan Voorspellen

Stel je voor dat je naar een enorme, woelige rivier kijkt. Dit is de lucht die door een straalmotor of een industriële brander stroomt. Op het eerste gezicht zie je alleen grote, draaiende stromingen (zoals grote golven). Maar als je door een microscoop kijkt, zie je dat deze grote stromingen uit miljoenen kleine, razendsnelle draaikolletjes bestaan. Deze kleine draaikolletjes zijn waar het vuur eigenlijk ontstaat en brandt.

Het probleem voor ingenieurs is dat computers te traag zijn om al die miljoenen kleine draaikolletjes in één keer te berekenen. Het zou te veel rekenkracht kosten. Daarom gebruiken ze een slimme truc: ze berekenen alleen de grote stromingen en gebruiken een "model" om te raden wat er in die kleine, onzichtbare draaikolletjes gebeurt.

Dit artikel van Sirignano en zijn team (van de Universiteit van Californië) introduceert een nieuwe, slimmere manier om dat model te laten werken.

Het Oude Probleem: De "Vreemde Variabele"

Vroeger probeerden ingenieurs de kleine vlammetjes te beschrijven door een kunstmatige variabele te gebruiken die ze een "progressievariabele" noemden. Dat is alsof je probeert te voorspellen hoe snel een auto rijdt door te kijken naar hoe ver de wielen gedraaid zijn, in plaats van naar de snelheidsmeter. Het werkt, maar het is niet direct gekoppeld aan de echte fysica van de stroming. Het is alsof je een kaart tekent zonder de wegen erop te hebben getekend; je moet er een beetje op gokken.

De Nieuwe Oplossing: De "Energie-Verbrandingsmeter" ($\epsilon$)

De auteurs zeggen: "Waarom gokken we? Laten we kijken naar de energie die in de stroming wordt verbruikt."

In de natuurkunde is er een getal dat $\epsilon$ (epsilon) heet. Dit is de dissipatie van turbulente kinetische energie. In gewone taal: dit is een maat voor hoe snel de energie van de grote, woelige stromingen wordt omgezet in warmte en wrijving door de aller-kleinste draaikolletjes.

Stel je voor dat de grote stroming een enorme, zware hamer is. De kleine draaikolletjes zijn de duizenden kleine spijkers waar die hamer op slaat. $\epsilon$ vertelt ons hoe hard die hamer op de spijkers slaat.

Het nieuwe idee in dit artikel is simpel maar krachtig:

1. De computer berekent de grote stroming en meet daar de "hamerkracht" ($\epsilon$).
2. Dit getal wordt direct doorgegeven aan het model van de kleine vlammetjes.
3. Het model gebruikt deze kracht om te berekenen hoe snel de lucht draait (vorticiteit) en hoe hard de vlam wordt uitgerekt (rek).

De Creatieve Analogie: De Dansende Vlam

Laten we een analogie gebruiken om dit te begrijpen:

Stel je een danser voor (de vlam) die op een dansvloer staat.

• De oude methode: De danser kreeg een instructie: "Dans zo snel als je kunt, gebaseerd op een getal dat we zelf hebben verzonnen." De danser wist niet precies hoe de vloer bewoog.
• De nieuwe methode: De danser krijgt een sensor op zijn voet die direct de trillingen van de vloer meet. Als de vloer heftig trilt (hoge $\epsilon$), weet de danser direct: "Oké, ik moet harder dansen en mezelf strakker houden."

Maar er is nog iets belangrijks in dit artikel: De draaiing (Vorticiteit).

In de oude modellen werd de draaiing van de lucht vaak genegeerd. Het was alsof je een danser liet dansen op een vloer die alleen op en neer ging, maar nooit ronddraaide.
De auteurs tonen aan dat de kleine draaikolletjes in een motor vaak draaien (zoals een kleine tornado). Deze draaiing zorgt voor een centrifugale kracht (net als wanneer je in een draaimolen zit en naar buiten wordt geduwd).

• Zonder draaiing: De vlam is als een platte lap stof die wordt uitgerekt.
• Met draaiing: De vlam is als een lap stof die in een centrifuge wordt gegooid. De draaiing duwt de lucht naar buiten, wat de vlam anders beïnvloedt dan alleen uitrekken.

Het artikel laat zien dat als je deze draaiing meeneemt in je berekeningen (het "Rotational Flamelet Model"), je veel nauwkeurigere resultaten krijgt. De vlam kan bijvoorbeeld tegen een hogere snelheid blijven branden voordat hij dooft, omdat de draaiing de brandstof en zuurstof op een specifieke manier samendrukt.

Waarom is dit belangrijk voor de wereld?

Dit klinkt als pure theorie, maar het heeft grote gevolgen:

1. Efficiëntere Motoren: Als we precies weten hoe vlammen werken in de kleinste details, kunnen we motoren bouwen die minder brandstof verbruiken en minder schadelijke uitstoot hebben.
2. Veiligheid: Het helpt bij het voorspellen wanneer een vlam dooft (extinctie) of wanneer hij te heet wordt, wat cruciaal is voor het ontwerp van veilige branders en raketten.
3. Geen "Gokwerk" meer: Door $\epsilon$ te gebruiken, hoeven ingenieurs geen kunstmatige variabelen meer te verzinnen. Ze koppelen de grote wereld (wat de computer ziet) direct aan de kleine wereld (waar het vuur brandt) via een natuurkundig getal dat al bestaat.

Samenvatting in één zin

In plaats van te gokken over hoe kleine vlammetjes zich gedragen in een wirwar van lucht, gebruiken deze onderzoekers de "energie-uitputting" van de grote stroming als een directe, betrouwbare schakel om de vlammen te voorspellen, waarbij ze ook rekening houden met de draaiende beweging van de lucht die de vlam een extra duw geeft.

Het is alsof ze eindelijk de sleutel hebben gevonden om de taal van de grote storm te vertalen naar de taal van de kleine vlam, zonder dat er een vertaler (een kunstmatige variabele) tussen hoeft te zitten.

Technische samenvatting

Probleemstelling

In de simulatie van turbulente verbranding (met RANS of LES) kunnen de kleinste schalen van turbulentie en chemische reacties niet direct worden opgelost. Deze sub-grid schalen worden traditioneel gemodelleerd met Flamelet-modellen. De huidige benaderingen, zoals het Klassieke Flamelet-model (CFM) of modellen gebaseerd op een Progress Variable (FPV), hebben echter fundamentele tekortkomingen:

1. Afhankelijkheid van aannames: Ze vertrouwen vaak op een vooraf verondersteld profiel voor de scalair dissipatiesnelheid ($\chi$) en gebruiken een "presumed PDF" (Probability Density Function) om turbulentiefluctuaties te koppelen aan de verbranding.
2. Verlies van fysieke koppeling: Het gebruik van een progressieve variabele (FPV) als koppeling tussen de opgeloste schaal en het sub-grid model leidt tot een verlies van de mechanische koppeling. De intrinsieke relatie tussen de rek snelheid (strain rate) en de vlamstructuur gaat verloren, omdat de progressieve variabele onafhankelijk evolueert van de lokale rek snelheid.
3. Verwaarlozing van werveling (Vorticity): De meeste modellen negeren het effect van werveling op de kleinste schalen (Kolmogorov-schaal), hoewel dit een centrifugaal effect heeft dat de moleculaire transportprocessen en de verbrandingssnelheid beïnvloedt.
4. Noodzaak van tracking variabelen: Er is een behoefte aan een methode die de sub-grid fysica direct koppelt aan de opgeloste schaal zonder kunstmatige tracking-variabelen te creëren.

Het doel van dit onderzoek is een nieuwe koppelingsprocedure te ontwikkelen die de turbulente kinetische energie dissipatiesnelheid ($\epsilon$) gebruikt als de primaire variabele om de mechanische randvoorwaarden (rek snelheid en werveling) voor het sub-grid flamelet-model te bepalen.

Methodologie

De auteurs stellen een theoretisch raamwerk voor dat de Rotational Flamelet Model (RFM) koppelt aan de opgeloste schaal via $\epsilon$.

1. Theoretische Basis:

   • De analyse richt zich op de kleinste turbulentieschalen (Kolmogorov-schaal), waar viskeuze dissipatie domineert en diffusie-gedreven verbranding het snelst verloopt.
   • Er wordt een transformatie uitgevoerd naar een roterend referentiekader (rotatie met de helft van de wervelingssterkte $\omega$) om de schuifrek te elimineren en drie hoofd-rek snelheden te definiëren.
   • De auteurs leiden relaties af tussen de kinetische energie dissipatiesnelheid ($\epsilon$), de viskeuze dissipatie ($\Phi$), de wervelingssterkte ($\omega^2$, enstrophy) en de rek snelheid ($S^*$).

2. Koppelingsmechanisme:

   • In plaats van $\chi$ (scalair dissipatie) als invoer te gebruiken, wordt $\epsilon$ (berekenbaar in RANS/LES) gebruikt om de mechanische constraints te bepalen.
   • Uit de vergelijkingen voor viskeuze dissipatie en kinetische energie wordt afgeleid dat $\epsilon$ direct gerelateerd is aan zowel de rek snelheid $S^*$ als de wervelingssterkte $\omega$.
   • De volgende relaties worden afgeleid (waarbij $C_{vd}$ en $C_{ke}$ dimensieloze coëfficiënten zijn):
     • $S^* \propto \sqrt{\epsilon/\nu}$
     • $\omega \propto \sqrt{\epsilon/\nu}$
   • Een cruciale voorwaarde voor het bestaan van een tegenstroom (counterflow) flamelet is dat de Laplacian van de druk negatief moet zijn. Dit vereist een balans waarbij de wervelingssterkte niet de rek snelheid overheerst ($C_{vd}/2 < C_{ke} < C_{vd}$).

3. Numerieke Implementatie:

   • Er worden drie modellen vergeleken:
     1. CFM: Klassiek flamelet-model zonder momentumvergelijking en zonder werveling (gebaseerd op FlameMaster).
     2. RFM (zonder werveling): Rotational Flamelet Model met $S_1=S_2$ en $\omega=0$.
     3. RFM (met werveling): Volledig RFM met $\omega \neq 0$.
   • Twee brandstof-oxidator combinaties worden getest: $H_2/N_2 - O_2$ en JP-5 (jetbrandstof) - lucht.
   • De resultaten worden geanalyseerd als functie van $S^*$ en $\sqrt{\epsilon/\nu}$.

Belangrijkste Bijdragen

1. $\epsilon$ als Koppelingsvariabele: Het paper introduceert $\epsilon$ niet als een "progress variable", maar als een koppelingsvariabele die de opgeloste schaal direct verbindt met de sub-grid mechanische constraints. Dit elimineert de noodzaak voor kunstmatige tracking-variabelen of progressieve variabelen.
2. Inclusie van Werveling: Het is een van de eerste werken dat systematisch de invloed van werveling op de kleinste schalen integreert in een flamelet-model. Dit leidt tot een meerwaarde in de fysieke nauwkeurigheid.
3. Meerwaardige Koppeling: De auteurs tonen aan dat er geen eenduidige relatie bestaat tussen de rek snelheid ($S^*$) en de scalair dissipatiesnelheid ($\chi$) wanneer werveling wordt meegenomen. Voor een gegeven $S^*$ kunnen verschillende $\chi$-waarden optreden afhankelijk van de werveling, wat leidt tot verschillende vlamstructuren.
4. Schalingswetten: Er wordt een schalingsrelatie afgeleid waarbij de scalair gradiënten schalen met $\epsilon^{1/4}$ en de scalair dissipatiesnelheid met $\epsilon^{1/2}$.

Resultaten

De numerieke resultaten voor zowel waterstof als JP-5 brandstof tonen de volgende inzichten:

• Temperatuur en Verbrandingssnelheid: Wanneer de resultaten worden geparametriseerd op basis van de scalair dissipatie ($\chi$), lijken de effecten van werveling minimaal (de S-vormige krommen overlappen grotendeels). Echter, wanneer geparametriseerd op basis van de rek snelheid ($S^*$) of $\epsilon$, blijkt dat werveling de extinctie-grens significant verhoogt (de vlam kan hogere rek snelheden weerstaan).
• Invloed op Verbranding: De aanwezigheid van werveling in het RFM leidt tot een significante afname van de geïntegreerde verbrandingssnelheid (heat release rate) in vergelijking met het model zonder werveling of het klassieke model. Dit komt door de centrifugale krachten die de menging en transport beïnvloeden.
• Scalair Dissipatie Profielen: De profielen van $\chi(Z)$ zijn niet langer uniek voor een gegeven $S^*$ wanneer werveling aanwezig is. Dit betekent dat het gebruik van een vooraf bepaald $\chi(Z)$ profiel (zoals in CFM) misleidend is en de fysieke realiteit van turbulente verbranding niet correct weergeeft.
• JP-5 vs. Waterstof: Hoewel JP-5 een lagere verbrandingssnelheid heeft dan waterstof, is het percentage toename van de ontvlambaarheidsgrens door werveling aanzienlijk groter voor JP-5 dan voor waterstof.

Betekenis en Conclusie

De studie concludeert dat het gebruik van de turbulente kinetische energie dissipatiesnelheid ($\epsilon$) als koppelingsvariabele een robuuste en fysiek onderbouwde methode biedt om sub-grid flamelet-modellen te koppelen aan RANS- of LES-simulaties.

• Noodzaak van Werveling: Een flamelet-model dat werveling negeert, kan de opgeloste schaal parameters niet correct vertalen naar de sub-grid fysica. Dit leidt tot onnauwkeurige voorspellingen van warmteafgifte en vlamstabiliteit.
• Voordeel van RFM: Het Rotational Flamelet Model, gekoppeld via $\epsilon$, biedt een vollediger beeld van de verbrandingsfysica door de centrifugale effecten van werveling mee te nemen.
• Toekomstige Toepassing: Hoewel de huidige studie zich richt op de theorie en steady-state resultaten, biedt deze methode een fundament voor toekomstige toepassing in complexe, tijdsafhankelijke turbulente verbrandingssimulaties zonder de complexiteit van extra transportvergelijkingen voor progressieve variabelen.

Kortom, dit werk verlegt de focus van het modelleren van "gemiddelde" scalair dissipatie naar het direct modelleren van de mechanische constraints (rek en werveling) die worden gedreven door de energie dissipatie van de turbulentie, wat leidt tot een nauwkeurigere en fysiek consistenter beschrijving van turbulente niet-gemengde verbranding.