cfdmfFTFoam: A front-tracking solver for multiphase flows on general unstructured grids in OpenFOAM

Dit artikel introduceert cfdmfFTFoam, een nieuwe OpenFOAM-oplosser die de Front-Tracking-methode implementeert op algemene ongestructureerde roosters om de beperkingen van bestaande open-source software voor meerfasestromingen te overwinnen en verdere onderzoek en algoritmeverbetering te faciliteren.

De kern

🌊 De Digitale Kookpan: Een nieuwe manier om vloeistoffen te simuleren

Stel je voor dat je een grote pan met water en olie op het vuur zet. De olie vormt druppels, die samensmelten, weer uit elkaar gaan en rondzwemmen. Als je wilt weten hoe dit precies werkt, moet je een computerprogramma gebruiken dat dit na-ijlt. Dit heet een "stroom-simulatie".

Maar er is een probleem: computers zijn niet zo slim als ze lijken. Als je twee vloeistoffen laat samensmelten (zoals een zeepbel die in de lucht vliegt en uit elkaar valt), raken de oude methoden vaak de draad kwijt. Ze "vergeten" precies waar de rand van de druppel zit, of ze laten de druppel op een rare manier uit elkaar vallen alsof het een digitale glitch is.

De auteurs van dit paper hebben een nieuw, krachtig programma gemaakt genaamd cfdmfFTFoam. Laten we kijken hoe dit werkt, zonder de saaie wiskunde.

1. Het Probleem: De "Vage" Rand

In de oude methoden (zoals VOF, wat vaak wordt gebruikt in OpenFOAM) wordt de rand tussen twee vloeistoffen als een "wazige foto" behandeld. Het is alsof je probeert de rand van een ijsklontje te tekenen door alleen de kleur van de lucht eromheen te veranderen. Dat werkt goed voor simpele dingen, maar als de druppel heel erg uitgerekt wordt, wordt de foto wazig en gaat de simulatie fout.

2. De Oplossing: Een Onzichtbaar Net (Front-Tracking)

De nieuwe methode, Front-Tracking, doet iets heel anders.

• De Oude Manier: Kijkt naar de hele pan en vraagt: "Is hier water of olie?"
• De Nieuwe Manier (cfdmfFTFoam): Plakt een onzichtbaar, digitaal net direct op de rand van de druppel.

Stel je voor dat je een balletje hebt en je bedekt het met een strakke, elastische sok van duizenden kleine driehoekjes. Die sok is het "front". De computer houdt die sok precies in de gaten. Waar de sok is, is de rand van de druppel. Dat is veel scherper en nauwkeuriger.

3. Het Grote Nadeel van de Oude Netten

Het probleem met die "sok" is dat computers die vaak alleen op simpele, rechthoekige roosters (zoals een schaakbord) kunnen gebruiken. Maar in de echte wereld zijn buizen en machines vaak krom en onregelmatig. Je kunt een schaakbord niet perfect om een ronde buis leggen zonder gaten of rare stukken.

Tot nu toe was er geen gratis software die dit "sok-net" kon gebruiken op onregelmatige roosters (zoals een net dat perfect om een kromme buis past).

4. De Innovatie: cfdmfFTFoam

De auteur, Ehsan Amani, heeft een brug gebouwd tussen twee werelden:

1. De "Sok" (Front-Tracking): De nauwkeurige methode om de rand te volgen.
2. OpenFOAM: Een beroemde, gratis software die al heel goed is in het simuleren van stromingen op onregelmatige roosters.

Hij heeft de code van de "sok" (die oorspronkelijk voor simpele roosters was gemaakt) aangepast zodat hij perfect werkt binnen OpenFOAM.

Wat doet dit programma nu precies?
Het programma doet een dans van zes stappen, telkens als de simulatie een fractie van een seconde verder gaat:

1. De Maatlat (Volume Correctie): Soms rekent de computer de grootte van de druppel net iets verkeerd uit (alsof de druppel een beetje krimpt of groeit zonder reden). Het programma meet dit en duwt de "sok" een heel klein beetje uit of in om de grootte perfect te houden.
2. De Knopen (Remeshing): Als de druppel heel erg uitrekt, worden de driehoekjes in de "sok" te groot of te klein. Het programma knipt ze dan netjes bij en plakt nieuwe stukjes erbij, zodat de sok strak blijft zitten.
3. De Spanning (Oppervlaktespanning): Het berekent hoe hard de "huid" van de druppel trekt (net zoals een zeepbel).
4. De Boodschapper (Communicatie): De "sok" moet de krachten vertellen aan het water eromheen. Het programma fungeert als een tolk die de boodschap van de druppel naar het water overbrengt.
5. De Kaart (Indicator): Het maakt een kaart van wie er waar zit (water hier, olie daar).
6. De Dans (Advectie): De hele "sok" wordt een stukje verplaatst volgens de stroming van het water.

5. Waarom is dit geweldig?

De auteurs hebben het getest met verschillende proefballen:

• De Deformatie-test: Ze lieten een druppel in een wervelende stroming dansen. De oude methoden lieten de druppel vaak "kapot" gaan of uit elkaar vallen. De nieuwe methode hield de druppel heel en nauwkeurig.
• De Zeepbel: Ze lieten een bel stijgen in een vloeistof. De nieuwe software voorspelde de vorm en snelheid veel beter dan de oude methoden, zelfs als de bel heel erg vervormde.
• De Taylor-bubbel: Een lange, rechte bel die door een buis stijgt. Dit is lastig omdat de bel heel lang kan zijn. De software hield de vorm perfect vast.

Het belangrijkste voordeel:
Het werkt op elk type rooster. Of je nu een simpele kubus simuleert of een ingewikkelde motor met kromme leidingen, dit programma past zich er perfect aan.

Conclusie

Voorheen was het heel moeilijk om complexe vloeistofproblemen (zoals hoe olie en water zich mengen in een ingewikkelde machine) nauwkeurig te simuleren zonder duizenden euro's aan software te betalen of zonder dat de computer vastliep.

Met cfdmfFTFoam hebben de onderzoekers een gratis, open-source tool gemaakt dat:

• De "sok-methode" gebruikt voor super-nauwkeurige randen.
• Werkt op complexe, onregelmatige vormen.
• Zorgt dat druppels niet "verdwijnen" of "kapotgaan" in de simulatie.

Het is alsof ze een nieuwe, superkrachtige bril hebben ontworpen voor ingenieurs, zodat ze precies kunnen zien wat er gebeurt als vloeistoffen met elkaar spelen, zelfs in de meest ingewikkelde situaties.

Technische samenvatting

Probleemstelling

De numerieke simulatie van meerfasestromingen (multiphase flows) vereist vaak methoden die de interface tussen vloeistoffen volledig oplossen (Fully Resolved Simulation - FRS). Hoewel er open-source software beschikbaar is voor interface-vangmethoden (zoals VOF en Level-Set), is er een aanzienlijk tekort aan open-source oplossingen voor de Front-Tracking Methode (FTM) op algemene ongestructureerde roosters.

• Bestaande pure FTM-codes (zoals PARIS, FronTier++, TrioCFD) zijn beperkt tot gestructureerde Cartesische roosters.
• Bestaande FTM-implementaties in OpenFOAM (zoals lentFoam) zijn hybride (Level-Set/FTM) en gebruiken geen zuiver Lagrangiaans oppervlak, wat de nauwkeurigheid van kromming en oppervlaktespanning beperkt.
• De complexiteit van het coderen van communicatie tussen een Lagrangiaans front en een Euleriaans rooster, vooral in parallelle omgevingen met ongestructureerde cellen, vormt een grote barrière.

Methodologie

De auteurs hebben een nieuwe open-source solver ontwikkeld genaamd cfdmfFTFoam, geïntegreerd in het OpenFOAM v9 CFD-framework. De kern van de methode is de integratie van de bestaande Fortran-code Ftc3D (ontwikkeld door Tryggvason en Dabiri) in de C++-omgeving van OpenFOAM.

Belangrijkste technische componenten:

1. Hybride Grid-Architectuur: De solver gebruikt een Euleriaans volumetrisch rooster (algemeen ongestructureerd, inclusief tetraëders, hexaeders, polyhedra) voor de Navier-Stokes vergelijkingen, en een Lagrangiaans oppervlakrooster (driehoekige mesh) om de interface te volgen.
2. Communicatie (Front-to-Field): Om data (zoals oppervlaktespanning) van het Lagrangiaanse front naar het Euleriaanse rooster te transfereren, wordt de Reproducing Kernel Particle Method (RKPM) geïmplementeerd. Dit lost het probleem op van het verdelen van krachten op een ongestructureerd rooster terwijl momenten (kracht en koppel) behouden blijven.
3. Sub-algoritmen: De solver bevat een uitgebreide bibliotheek (frontTracking) met specifieke FTM-algoritmen:
   • Volume Correctie (VC2): Een algoritme om massabehoud te garanderen door de posities van front-vertices te corrigeren op basis van een kubieke vergelijking.
   • Remeshing: Dynamisch verfijnen (refining), grover maken (coarsening) en gladmaken (smoothing/undulation removal) van het frontrooster om de kwaliteit van de driehoeken te handhaven.
   • Oppervlaktespanning: Berekening via een directe element-gebaseerde methode.
   • Indicator Functie: Constructie van de fase-indicator via een Poisson-vergelijking of Closest Point Transform (CPT).
4. Parallelisatie: De code is ontworpen voor domeindecompositie in parallelle runs, waarbij synchronisatie tussen de processoren en de front-vertices wordt beheerd.

Kernbijdragen

• Eerste Pure FTM voor OpenFOAM: Dit is de eerste open-access, pure Front-Tracking solver die werkt op algemene ongestructureerde roosters binnen OpenFOAM.
• Integratie en Uitbreiding: Het succesvol integreren van legacy Fortran-code (Ftc3D) in een moderne C++-framework, aangevuld met nieuwe algoritmen voor volume-correctie en communicatie op ongestructureerde grids.
• Open Source Beschikbaarheid: De code is beschikbaar via GitHub en de CPC Library, wat de toegankelijkheid voor de onderzoeksgemeenschap vergroot.
• Flexibiliteit: Ondersteuning voor verschillende roostertypen (hexaedrisch, prismatisch, polyhedrisch, tetraëdrisch) en schaalbaarheid op meerdere processors.

Resultaten en Validatie

De solver is gevalideerd tegen diverse standaard benchmarks:

1. 3D Deformatie Benchmark: De solver toonde een advectiefout van slechts 1%, vergeleken met 36-41% voor traditionele VOF-methoden (MULES en isoAdvector). Dit demonstreert het vermogen van FTM om scherpe interfaces te behouden zonder kunstmatige breuk.
2. Stilstaande Druppel: De berekende drukverschilafwijking was slechts 0,046% en de parasitaire stromingen waren zeer laag, wat superieur is aan VOF-resultaten (12,6% fout).
3. Druppel in Poiseuille-stroming: De solver leverde nauwkeurige resultaten voor verschillende roostertypen (inclusief tetraëders) en toonde goede parallelle schaalbaarheid (speedup), hoewel iets lager dan de standaard interFoam solver.
4. Taylor-bubbel en Opstijgende Bellen: De solver kon complexe fenomenen zoals het rijzen van een Taylor-bubbel en een bel met hoge dichtheidsverhoudingen (1000:1) simuleren. Hoewel er bij zeer hoge dichtheidsverhoudingen nog enige onderpredictie van de vervorming was (geassocieerd met discretisatie-inconsistenties in de momentum-vergelijking), bevestigden de resultaten de robuustheid van de methode.

Betekenis en Toekomstperspectief

De paper introduceert een krachtig nieuw instrument voor de CFD-gemeenschap dat de kloof overbrugt tussen de hoge nauwkeurigheid van de Front-Tracking Methode en de flexibiliteit van ongestructureerde roosters in OpenFOAM.

• Wetenschappelijke Impact: Het stelt onderzoekers in staat om complexe meerfasestromingen (zoals coalescentie, breuk, en contactlijndynamica) met hogere nauwkeurigheid te simuleren dan met traditionele interface-vangmethoden mogelijk is.
• Toekomst: De basis is gelegd voor verdere ontwikkeling, waaronder het toevoegen van contactlijndynamica, topologiewijzigingen (splitsen/samenvoegen), en interfaciale warmte- en massatransport. De code dient als een open platform voor het verbeteren van FTM-algoritmen.

Samenvattend biedt cfdmfFTFoam een robuust, nauwkeurig en open-source alternatief voor bestaande methoden, specifiek ontworpen voor complexe geometrieën en parallelle simulaties van meerfasestromingen.