A Volume of Fluid Immersed Boundary Method for Industrial Polymer Mixing

Dit artikel presenteert een nieuwe block-gekoppelde Volume-of-Fluid-geïmmanseerde grenslaag (BC-VOF-IB) solver, geïmplementeerd in OpenFOAM, die numerieke instabiliteiten veroorzaakt door hoge viscositeitscontrasten overwint om vrije-oppervlakte polymeermenging in gedeeltelijk gevulde industriële extruders nauwkeurig te simuleren.

De kern

Stel je voor dat je probeert een gigantische ketel met dikke, plakkerige honing met lucht te mengen in een draaiende machine. Dit is in wezen wat er gebeurt bij industriële polymerenmengen, waarbij bedrijven zoals Pirelli gesmolten kunststof moeten mengen met additieven om banden, medische hulpmiddelen of auto-onderdelen te maken. Het doel is om alles perfect te mengen zodat het eindproduct sterk en uniform is.

Het simuleren van dit proces op een computer is echter een nachtmerrie voor wiskundigen en ingenieurs. Hier is waarom, en hoe dit artikel het oplost, met behulp van eenvoudige analogieën:

Het Probleem: De Strijd tussen "Dikke Honing en Dunne Lucht"

In deze machines heb je twee zeer verschillende vloeistoffen:

1. Polymerensmelt: Extreem dik, plakkerig en traag bewegend (zoals koude honing).
2. Lucht: Zeer dun en snel bewegend.

Wanneer je probeert te simuleren hoe deze twee met elkaar interageren in een machine met draaiende schroeven, raken standaard computerprogramma's in de war. Het is alsof je probeert de beweging van een slak en een raceauto op hetzelfde spoor te berekenen met dezelfde set regels. De computer probeert heel kleine, heel kleine stappen te nemen om te voorkomen dat de "slak" (de dikke kunststof) te snel beweegt, wat de simulatie ongelooflijk traag maakt – soms duurt het dagen om een paar seconden real-time mengen te voltooien.

Bovendien hebben de machines complexe, draaiende onderdelen (schroeven) die zich binnen een vaste container bewegen. Traditioneel moet je, om dit te simuleren, een digitaal rooster (een raster van kleine vakjes) bouwen dat perfect om de draaiende schroeven heen past. Terwijl de schroeven draaien, moet dit rooster zich voortdurend opnieuw vormen, wat net zo is als proberen een trui te breien terwijl de persoon die hem draagt een marathon loopt. Het is rommelig, moeilijk en vatbaar voor fouten.

De Oplossing: Een Nieuw "Slim Rooster" en een "Team-aanpak"

De auteurs van dit artikel hebben een nieuwe manier ontwikkeld om deze simulaties uit te voeren met behulp van de software OpenFOAM. Ze combineerden twee krachtige technieken:

1. De Immersed Boundary Method (De "Spookmuur"-truc)
In plaats van het rooster opnieuw te vormen om de draaiende schroeven te passen, hielden ze het rooster vast en stijf (zoals een vast blok ijs). Ze vertelden de computer vervolgens: "Hé, er zit een draaiende schroef in dit blok ijs."

• De Analogie: Stel je een zwembad voor met een vast raster van tegels op de bodem. In plaats van de tegels te verplaatsen om een zwemmer te passen, zeg je gewoon tegen het water: "Ga niet door de zwemmer heen." De computer gebruikt wiskunde om een "spookmuur" rond de schroef te creëren, waardoor de vloeistof eromheen moet stromen zonder dat het rooster ooit opnieuw hoeft te worden gebouwd. Dit maakt het verwerken van complexe, bewegende vormen veel eenvoudiger.

2. De Volume of Fluid (VOF) Method (De "Verf Volgen"-truc)
Om te zien waar de dikke kunststof eindigt en de lucht begint, gebruiken ze een "verf" die de cellen vult.

• De Analogie: Stel je het computerraster voor als een 3D-schakenbord. Sommige vakken zijn 100% kunststof, sommige 100% lucht en sommige een mengsel. De computer volgt hoeveel "kunststofverf" er in elk vakje zit om het oppervlak van de vloeistof te zien.

3. Het Block-Coupled Schema (De "Teamhuddle")
Dit is de belangrijkste doorbraak. In standaard simulaties lost de computer de snelheid van de vloeistof in de X-, Y- en Z-richting één voor één op, alsof drie mensen om de beurt praten. Wanneer de vloeistof superdik is (zoals polymeren), zorgt deze "om de beurt"-aanpak ervoor dat de simulatie crasht of tot een slakken tempo vertraagt, omdat de dikke vloeistof alle richtingen strak aan elkaar koppelt.

De auteurs veranderden dit in een Block-Coupled-aanpak.

• De Analogie: In plaats van dat drie mensen om de beurt praten, huddelen ze allemaal samen en lossen ze het probleem tegelijkertijd op. Door de beweging in alle richtingen te behandelen als één groot, onderling verbonden team, kan de computer het enorme verschil tussen de dikke kunststof en de dunne lucht verwerken zonder vast te lopen.

De Resultaten: Van Uren naar Minuten

Het team testte hun nieuwe methode op twee scenario's:

1. Een Kanaal in Hondenbot-vorm: Een testgeval waarbij dikke kunststof wordt geïnjecteerd in een smal, kronkelend kanaal.

   • Oude manier: Het standaard computerprogramma crashte of duurde 7 uur om een paar seconden te simuleren omdat het gedwongen werd om heel kleine stappen te nemen.
   • Nieuwe manier: Hun nieuwe "Teamhuddle"-methode voltooide dezelfde taak in slechts 16 minuten en crashte niet, zelfs niet toen de kunststof extreem dik werd.

2. Echte Industriële Machines: Ze simuleerden echte single-schroef en twin-schroef extruders (de machines die worden gebruikt om kunststofkorrels te maken).

   • Ze toonden succesvol aan hoe de kunststof de machine vult, hoe de druk opbouwt en hoe de lucht wordt weggeduwd.
   • Ze bewezen dat hun "Spookmuur"-methode net zo goed werkt als de oude, moeilijke methode van het opnieuw vormen van het rooster, maar veel sneller en makkelijker in te stellen is.

Wat Komt Er Vervolgens?

Het artikel concludeert dat dit een grote stap voorwaarts is voor de industrie. Het overbrugt de kloof tussen academische wiskunde en echte fabrieksbehoeften. De auteurs merken echter op dat hun huidige model ervan uitgaat dat de temperatuur gelijk blijft (isotherm). In werkelijkheid genereert het mengen van kunststof warmte, wat verandert hoe dik de kunststof is. Het toevoegen van temperatuureffecten en meer complexe "rekbaar" kunststofgedrag zijn de volgende stappen voor toekomstig onderzoek.

Kortom: Ze bouwden een snellere, stabielere manier om te kijken hoe dikke kunststof met lucht wordt gemengd in draaiende machines op een computer, waardoor een proces dat vroeger uren duurde, nu minuten duurt, zonder dat je elke keer dat een schroef draait de digitale wereld hoeft te herbouwen.

Technische samenvatting

Technische Samenvatting: Een Volume-of-Fluid Immersed Boundary-methode voor industriële polymerenmenging

Probleemstelling
De numerieke simulatie van mengprocessen voor polymeren vormt aanzienlijke uitdagingen vanwege de complexe wisselwerking tussen niet-Newtoniaanse reologie, meerfasige dynamiek met een vrij oppervlak en bewegende complexe geometrieën. Industriële mixers, zoals enkelschroefextruders (SSE) en tweeschroefextruders (TSE), werken vaak onder gedeeltelijk gevulde omstandigheden waarbij hoogviskeuze polymeersmelten samengaan met lucht. Het simuleren van deze scenario's vereist het vastleggen van het interface tussen fasen, terwijl er tegelijkertijd met roterende vaste grensvlakken wordt omgegaan.

Standaard Computational Fluid Dynamics (CFD)-benaderingen staan in deze context voor twee hoofdhobbels:

1. Geometrische Complexiteit: Het genereren van lichaam-aangepaste (conforme) meshes voor roterende schroeven met complexe geometrieën is rekenkundig duur en vereist vaak dynamische mesh-technieken die kunnen leiden tot meshvervorming.
2. Numerieke Instabiliteit: Het extreme viscositeitscontrast tussen polymeersmelten en lucht (vaak hoger dan $10^5$ tot $10^8$) veroorzaakt ernstige numerieke instabiliteiten in standaard gesegregeerde Volume-of-Fluid (VOF)-oplossers. Specifiek legt de expliciete behandeling van de getransponeerde snelheidsgradiënt in de viskeuze diffusie-term zware tijdstapbeperkingen op ($\Delta t \propto h^2$), wat simulaties van hoogviskeuze stromingen rekenkundig onbetaalbaar of instabiel maakt.

Methodologie
Dit werk stelt een nieuw numeriek raamwerk voor dat een Volume-of-Fluid (VOF) interface-vastleggende aanpak integreert met een niet-conforme Immersed Boundary (IB)-methode, geïmplementeerd binnen de OpenFOAM Finite Volume-bibliotheek.

• Immersed Boundary (IB)-Aanpak: De methode maakt gebruik van een vast achtergrondrooster waarbij vaste grensvlakken (bijv. roterende schroeven) worden weergegeven door getrianguleerde STL-oppervlakken. Een Weighted Least Squares (WLS)-interpolatie-operator wordt ingezet om randvoorwaarden (Dirichlet voor snelheid, Neumann voor volumefractie) af te dwingen op cellen die door het immersed boundary worden gesneden. Dit elimineert de noodzaak voor het genereren van conforme meshes en maakt het hanteren van bewegende contactlijnen langs complexe, roterende geometrieën mogelijk.
• VOF-Formulering: Het interface wordt gevolgd door het oplossen van een transportvergelijking voor de volumefractie van de fase ($\alpha$). De methode maakt gebruik van het MULES-algoritme (Multidimensional Universal Limiter with Explicit Solution) om scherpte en begrenzing van het interface te behouden, inclusief een compressieve snelheidsterm om numerieke diffusie tegen te gaan.
• Block-gekoppelde (BC)-Schematiek: Om de stabiliteitsproblemen die voortvloeien uit hoge viscositeitscontrasten aan te pakken, introduceren de auteurs een block-gekoppeld schema voor de impulsvergelijking. In tegenstelling tot standaard gesegregeerde oplossers die snelheidscomponenten sequentieel behandelen, behandelt de BC-aanpak de volledige viskeuze diffusie-term $\nabla \cdot (\mu(\nabla \mathbf{u} + \nabla^\top \mathbf{u}))$ impliciet. Dit houdt in dat de drie snelheidscomponenten worden gekoppeld binnen één enkel lineair systeem. De discretisatie van de tangentiële gradiëntterm maakt gebruik van een Finite Area (FA)-methode in combinatie met WLS-interpolatie op uitgebreide stencilen (waarbij puntburen worden betrokken in plaats van alleen vlakburen), aangepast van toepassing op elasticiteit [11].
• Oplossingsalgoritme: De solver maakt gebruik van een gewijzigd PIMPLE-algoritme (combinatie van SIMPLE en PISO) om de druk-snelheidskoppeling te behandelen. De IB-beperkingen worden geïntegreerd in het algebraïsche systeem, waarbij de matrixcoëfficiënten en brontermen worden gewijzigd om randvoorwaarden op de immersed oppervlakken te voldoen, terwijl massabehoud wordt gehandhaafd.

Belangrijkste Bijdragen

1. Ontwikkeling van BC-VOF-IB-oplossers: Het artikel introduceert de UCInterIbFoam-solver, die de niet-conforme IB-methode combineert met de block-gekoppelde VOF-aanpak. Dit is de eerste implementatie van een dergelijke gekoppelde methode binnen het OpenFOAM-10-raamwerk voor tweefasige stromingen.
2. Stabiliteitsverbetering via Block-koppeling: Het werk toont aan dat de volledig impliciete behandeling van viskeuze diffusie de stabiliteitsbeperkingen voor de tijdstap aanzienlijk verlicht. Dit maakt simulaties van hoogviskeuze polymeerstromingen met viscositeitsverhoudingen tot $10^8$ mogelijk zonder dat de simulatie crasht door schendingen van het Courant-getal.
3. Hanteren van Bewegende Contactlijnen: Het raamwerk adresseert succesvol de geometrische en algoritmische moeilijkheden van bewegende contactlijnen waar het tweefasige interface immersed grensvlakken snijdt, een kritieke vereiste voor gedeeltelijk gevulde mengapparatuur.

Numerieke Resultaten
De voorgestelde methoden zijn gevalideerd en getest op twee categorieën problemen:

• Benchmarks (Dog-bone Injectie- molding):

  • Er werd een vergelijking uitgevoerd tussen gesegregeerde (interFoam), block-gekoppelde conforme (UCInterFoam) en block-gekoppelde niet-conforme (UCInterIbFoam) oplossers voor een injectiegeval van een hoogviskeuze Power-law-vloeistof.
  • De standaard gesegregeerde solver faalde om binnen redelijke tijdframes stabiele resultaten te produceren vanwege ernstige tijdstapbeperkingen ($\Delta t \approx 10^{-6}$ s).
  • De block-gekoppelde oplossers behaalden stabiele oplossingen met tijdstappen tot $10^{-3}$ s, waardoor de rekentijd werd gereduceerd van 7 uur naar 16 minuten (conform) en 58 minuten (niet-conform).
  • De niet-conforme IB-oplossing toonde kwalitatieve overeenstemming met de conforme oplossing, met kleine afwijkingen in viscositeitspieken nabij geometrische restricties, toegeschreven aan de niet-conforme aard van het rooster.

• Industriële Toepassingen (SSE en TSE):

  • Simulaties werden uitgevoerd op vereenvoudigde enkelschroefextruders (SSE) en realistische tweeschroefextruders (TSE) onder zowel volledig gevulde als gedeeltelijk gevulde (starve-fed) omstandigheden.
  • De oplossers slaagden erin de evolutie van het vrije oppervlak, vulgraden en drukvelden vast te leggen.
  • De resultaten kwamen overeen met fysieke verwachtingen: volledig gevulde scenario's ($Q_{in} > Q_0$) vertoonden negatieve drukgradiënten van inlaat tot uitlaat, terwijl gedeeltelijk gevulde scenario's ($Q_{in} < Q_0$) drukstijgingen nabij de uitlaat lieten zien.
  • De IB-methode bleek effectief in het hanteren van de complexe, verstrengelde geometrie van de TSE zonder dat dynamische meshregeneratie nodig was.

Betekenis en Beweringen
Het artikel beweert dat het voorgestelde BC-VOF-IB-raamwerk een betekenisvolle stap voorwaarts is om de kloof tussen academisch CFD-onderzoek en industriële eisen voor polymerenverwerking te overbruggen. Door de numerieke instabiliteiten geassocieerd met hoge viscositeitscontrasten te overwinnen en de noodzaak voor complexe lichaam-aangepaste meshgeneratie voor roterende geometrieën te elimineren, maakt de methode fysiek consistente voorspellingen van snelheids- en drukvelden mogelijk in industriële, gedeeltelijk gevulde mengapparatuur.

De auteurs erkennen dat het huidige werk beperkt is tot isotherme processen en gegeneraliseerde Newtoniaanse reologie. Zij stellen dat toekomstige ontwikkelingen thermische effecten (energievergelijking) en visco-elastische constitutieve modellen moeten omvatten om de fysica van industriële polymerenmenging volledig vast te leggen. Het huidige raamwerk legt echter een rigoureuze basis voor deze uitbreidingen en biedt een robuust hulpmiddel voor het simuleren van stromingen met een vrij oppervlak in complexe, bewegende geometrieën.