An LES model with finite-rate phase change and subgrid spray based on a thermodynamically consistent four-equation multiphase model

Dit artikel presenteert een nieuw LES-model voor meervasige simulaties dat gebruikmaakt van een thermodynamisch consistent vier-vergelijkingenmodel, een verbeterde $\Sigma$-sprayformulering en een eindige faseovergang om complexe spuitprocessen nauwkeurig en efficiënt te voorspellen.

De kern

Stel je voor dat je een supercomputer probeert te gebruiken om een perfecte kop koffie te zetten, maar de computer moet niet alleen de temperatuur berekenen, maar ook hoe de melk zich mengt, hoe de stoom opstijgt en hoe de kleine druppeltjes vloeistof door de lucht vliegen. Dat is precies wat deze wetenschappers proberen te doen, maar dan voor de extreem ingewikkelde brandstof die in een raketmotor of een dieselmotor wordt gespoten.

Hier is de uitleg van hun onderzoek in begrijpelijke taal:

Het probleem: De "Digitale Mist"

Als wetenschappers een motor simuleren op een computer, proberen ze de vloeistof (brandstof) te laten verdampen tot gas. Maar in de echte wereld gebeurt dit razendsnel en op een chaotische manier.

Het probleem is dat computers niet "slim" genoeg zijn om elk minuscuul druppeltje te zien. Als je probeert elk klein druppeltje apart te tekenen, raakt de computer oververhit en loopt hij vast. Het is alsof je probeert een foto te maken van een waterval waarbij je elk individueel watermolecuul apart wilt fotograferen: dat lukt simpelweg niet.

De oplossing: De "Slimme Sluier"

De onderzoekers hebben een nieuwe manier gevonden om dit op te lossen zonder de computer te laten crashen. Ze gebruiken een soort "slimme sluier" (het $\Sigma$-model).

In plaats van elk druppeltje te tellen, berekenen ze de totale oppervlakte van alle onzichtbare druppeltjes samen. Denk aan een mistbank: de computer ziet niet elk waterdruppeltje, maar hij weet wel hoe dik de mist is en hoeveel oppervlakte die mist heeft. Hoe meer oppervlakte de "mist" heeft, hoe sneller de brandstof kan verdampen.

De "Thermodynamische Rem"

Een ander probleem bij computersimulaties is dat ze soms "valsspelen". Ze kunnen berekenen dat brandstof plotseling verandert in gas, zelfs als dat natuurkundig gezien onmogelijk is (bijvoorbeeld omdat het te koud is). Het is alsof je een film ziet waarin ijs smelt in een vriezer; dat klopt niet.

De onderzoekers hebben een nieuwe "thermodynamische rem" ingebouwd. Ze hebben een wiskundige regel toegevoegd die zegt: "Ho even! Je mag alleen maar zoveel brandstof laten verdampen als de natuurwetten op dit moment toelaten." Ze gebruiken hiervoor de Gibbs vrije energie—een soort natuurkundige boekhouder die controleert of de balans tussen vloeistof en gas nog wel klopt.

Hoe werkt het in de praktijk? (De vergelijking)

Stel je een drukke discotheek voor:

1. De Brandstof: De mensen in de club.
2. Verdamping: De mensen die de club verlaten via de uitgang.
3. Het $\Sigma$-model (De oppervlakte): In plaats van elke persoon te tellen, kijkt de computer naar hoe breed de deuren openstaan. Hoe breder de deuren (meer oppervlakte), hoe meer mensen er tegelijk naar buiten kunnen.
4. De Finite-Rate Phase Change (De snelheid): De mensen rennen niet allemaal tegelijk naar buiten; ze hebben een bepaalde snelheid.
5. De Thermodynamische Rem: De uitsmijter bij de deur. Zelfs als de deuren wagenwijd openstaan, laat de uitsmijter niemand naar buiten als de club volgens de regels al "leeg" is.

Waarom is dit belangrijk?

Door deze slimme wiskundige trucjes kunnen we met veel minder rekenkracht toch heel nauwkeurige voorspellingen doen. Dit helpt ingenieurs om:

• Betere motoren te bouwen die minder brandstof verbruiken.
• Veiligere raketten te ontwerpen die efficiënter brandstof verbranden.
• Minder vervuiling te veroorzaken, omdat we precies weten hoe de brandstof verbrandt.

Kortom: Ze hebben een manier gevonden om de chaos van vloeistof en gas te vangen in een slimme, snelle en natuurkundig correcte digitale formule.

Technische samenvatting

Technische Samenvatting: LES-model met eindige faseverandering en subgrid-spray op basis van een thermodynamisch consistent vier-vergelijkingen multiphase model

1. Probleemstelling (Het Probleem)

Het modelleren van complexe, meervoudige fase- en multicomponent-systemen (zoals brandstofinjectie in motoren of raketvoortstuwing) op industriële schaal is computationeel extreem zwaar. De grootste uitdagingen zijn:

• Computationele kosten: Volledige niet-evenwichtsmodellen (zoals zeven-vergelijkingen modellen) die de mechanische en thermische gradiënten tussen fasen expliciet berekenen, zijn te duur voor Large Eddy Simulations (LES) op grote schaal.
• Fysische nauwkeurigheid: Veel modellen gebruiken het Homogeneous Equilibrium Model (HEM), dat ervan uitgaat dat faseverandering onmiddellijk gebeurt. Dit is onrealistisch voor subgrid-structuren (zoals kleine druppels of vloeistofligamenten) die een eindige oppervlakte hebben en dus een eindige snelheid van faseverandering vertonen.
• Subgrid-fysica: In LES-simulaties zijn veel spray-kenmerken (zoals de interface-oppervlakte en kleine druppels) kleiner dan de roosterresolutie, wat leidt tot onnauwkeurigheden in massatransport en spray-verspreiding.

2. Methodologie (De Aanpak)

De auteurs stellen een uitgebreid numeriek kader voor dat gebruikmaakt van een vier-vergelijkingen multiphase model. Dit model gaat uit van een strikt subgrid-evenwicht voor druk, temperatuur en snelheid, wat enorme rekenvoordelen biedt. Om de nauwkeurigheid te behouden, voegen ze de volgende componenten toe:

• Thermodynamisch gebonden eindige faseverandering: In plaats van een constante relaxatietijd te gebruiken, introduceren ze een model gebaseerd op de Hertz-Knudsen benadering. Cruciaal is dat dit model "thermodynamisch begrensd" is: de faseverandering wordt gestuurd door de relaxatie van de chemische potentiaal (Gibbs vrije energie), waardoor het model nooit fysiek onmogelijke toestanden bereikt (zoals overmatige verdamping voorbij het evenwicht).
• Phase-Confined $\Sigma$-spray model: Om de subgrid-oppervlakte ($\Sigma$) te modelleren, gebruiken ze een aangepaste Eulerian $\Sigma$-vergelijking. Een innovatie hier is de "phase-confined" vorm: de diffusie van de oppervlakte-dichtheid wordt beperkt tot de gasfase. Dit voorkomt "unphysical leakage" (het onrealistisch lekken van vloeistofkenmerken naar de gasfase door numerieke diffusie).
• Spray-regime indicatoren: Het model kan automatisch onderscheid maken tussen verschillende fysieke regimes:
  • Dense vs. Dilute: Gebruikmakend van de vloeistofvolumefractie en de Hinze-schaal.
  • Flashing vs. Evaporation: Gebruikmakend van het Jakob-getal om te bepalen of de faseverandering wordt gedomineerd door snelle verdamping (flashing) of geleidelijke verdamping.

3. Belangrijkste Bijdragen (Key Contributions)

1. Efficiëntie door evenwicht: Het bewijzen dat een vier-vergelijkingen model (met mechanisch/thermisch evenwicht) zeer nauwkeurige resultaten kan leveren als de faseverandering en subgrid-spray correct worden gemodelleerd.
2. Nieuwe $\Sigma$-transportvergelijking: Een verbeterde transportvergelijking voor subgrid-oppervlakte die numeriek robuust is en lekken over interfaces voorkomt.
3. Consistentie in faseverandering: Een nieuw mechanisme waarbij de snelheid van faseverandering (Hertz-Knudsen) wordt gekoppeld aan de thermodynamische limiet van het chemische evenwicht.
4. Hybride brontermen: Integratie van brontermen voor turbulente breuk (breakup), coalescentie en flashing binnen één consistent LES-kader.

4. Resultaten (De Resultaten)

Het model is gevalideerd met de Engine Combustion Network (ECN) Spray A case (injectie van n-dodecaan in stikstof).

• Niet-verdampende condities: De simulaties vertoonden uitstekende overeenkomst met experimentele data voor de projectie van de massadichtheid (PMD), de oppervlakteverdeling (PSA) en de Sauter Mean Diameter (SMD). Dit bevestigt dat de spray-dynamica en de subgrid-oppervlaktemodellering correct werken.
• Verdampende condities: De simulaties voorspelden de penetratielengte van zowel de vloeistof als de gaswolk zeer nauwkeurig. De overeenkomst met experimentele DBI- en Mie-scattering-metingen was zeer hoog, wat aantoont dat de gekoppelde faseveranderingsmodellen de complexe interactie tussen spray-verspreiding en massatransport correct vangen.

5. Betekenis (Significance)

Dit werk biedt een krachtig en efficiënt instrument voor ingenieurs die complexe multiphase-stromingen willen simuleren. Het overbrugt de kloof tussen de lage rekenkosten van evenwichtsmodellen en de hoge fysieke nauwkeurigheid van niet-evenwichtsmodellen. De methodologie is direct toepasbaar op het optimaliseren van brandstofinjectiesystemen in motoren, raketmotoren en brandblussystemen, waarbij zowel de spray-structuur als de thermodynamische faseovergangen cruciaal zijn.