A Volume of Fluid Immersed Boundary Method for Industrial Polymer Mixing

Dieser Beitrag stellt einen neuartigen blockgekoppelten Volume-of-Fluid-Immersed-Boundary (BC-VOF-IB)-Löser vor, der in OpenFOAM implementiert ist, numerische Instabilitäten infolge hoher Viskositätskontraste überwindet und eine präzise Simulation der Freioberflächen-Mischpolymerisation in teilweise gefüllten industriellen Extrudern ermöglicht.

Das Wesentliche

Stellen Sie sich vor, Sie versuchen, einen riesigen Kessel aus dickflüssigem, klebrigem Honig mit Luft in einer rotierenden Maschine zu mischen. Genau dies geschieht bei der industriellen Polymermischung, wo Unternehmen wie Pirelli geschmolzenen Kunststoff mit Additiven vermischen müssen, um Reifen, medizinische Geräte oder Autoteile herzustellen. Das Ziel ist es, alles perfekt zu mischen, damit das Endprodukt stabil und gleichmäßig ist.

Die Simulation dieses Prozesses am Computer ist jedoch ein Albtraum für Mathematiker und Ingenieure. Hier ist der Grund und wie diese Arbeit das Problem löst, unter Verwendung einfacher Analogien:

Das Problem: Der Kampf „Dickflüssiger Honig gegen Dünne Luft"

In diesen Maschinen befinden sich zwei sehr unterschiedliche Fluide:

1. Polymer-Schmelze: Extrem dickflüssig, klebrig und langsam bewegend (wie kalter Honig).
2. Luft: Sehr dünn und schnell bewegend.

Wenn man versucht zu simulieren, wie diese beiden in einer Maschine mit rotierenden Schneckenelementen interagieren, geraten Standard-Computerprogramme in Verwirrung. Es ist wie der Versuch, die Bewegung einer Schnecke und eines Rennwagens auf derselben Strecke unter Verwendung desselben Regelwerks zu berechnen. Der Computer versucht, winzige, winzige Schritte zu machen, um zu verhindern, dass die „Schnecke" (der dicke Kunststoff) zu schnell wird, was die Simulation unglaublich langsam macht – manchmal dauert es Tage, um ein paar Sekunden Echtzeit-Mischung zu simulieren.

Darüber hinaus verfügen die Maschinen über komplexe, rotierende Teile (Schnecken), die sich innerhalb eines festen Behälters bewegen. Traditionell muss man für die Simulation ein digitales Gitter (ein Raster aus winzigen Kästchen) erstellen, das sich perfekt um die rotierenden Schnecken legt. Während sich die Schnecken drehen, muss dieses Gitter sich ständig neu formen, was so ist, als würde man einen Pullover stricken, während die Person, die ihn trägt, einen Marathon läuft. Es ist chaotisch, schwierig und fehleranfällig.

Die Lösung: Ein neues „Intelligentes Gitter" und ein „Team-Ansatz"

Die Autoren dieser Arbeit entwickelten eine neue Methode, um diese Simulationen mit einer Software namens OpenFOAM durchzuführen. Sie kombinierten zwei leistungsstarke Techniken:

1. Die Immersed-Boundary-Methode (Der „Geisterwand"-Trick)
Anstatt das Gitter an die rotierenden Schnecken anzupassen, behielten sie das Gitter fest und starr (wie einen festen Eisblock). Dann sagten sie dem Computer: „Hey, in diesem Eisblock befindet sich eine rotierende Schnecke."

• Die Analogie: Stellen Sie sich einen Swimmingpool vor, dessen Boden mit einem festen Gitter aus Fliesen bedeckt ist. Anstatt die Fliesen zu bewegen, um sie an einen Schwimmer anzupassen, sagen Sie dem Wasser einfach: „Geh nicht durch den Schwimmer." Der Computer nutzt Mathematik, um eine „Geisterwand" um die Schnecke zu erzeugen, die das Fluid zwingt, daran vorbeizuströmen, ohne dass das Gitter jemals neu aufgebaut werden muss. Dies macht den Umgang mit komplexen, sich bewegenden Formen viel einfacher.

2. Die Volume-of-Fluid-Methode (VOF) (Der „Farbverfolgungs"-Trick)
Um zu sehen, wo der dicke Kunststoff aufhört und die Luft beginnt, verwenden sie eine „Farbe", die die Zellen füllt.

• Die Analogie: Stellen Sie sich das Computergitter als ein 3D-Schachbrett vor. Einige Felder sind zu 100 % Kunststoff, einige zu 100 % Luft und einige sind eine Mischung. Der Computer verfolgt, wie viel „Kunststofffarbe" sich in jedem Feld befindet, um die Oberfläche der Flüssigkeit zu erkennen.

3. Das Block-gekoppelte Schema (Das „Team-Huddle")
Dies ist der wichtigste Durchbruch. Bei Standard-Simulationen berechnet der Computer die Geschwindigkeit des Fluids in den X-, Y- und Z-Richtungen nacheinander, wie drei Personen, die sich abwechselnd äußern. Wenn das Fluid extrem dick ist (wie Polymer), führt dieser „Abwechseln"-Ansatz dazu, dass die Simulation abstürzt oder sich ins Unendliche verlangsamt, da das dicke Fluid alle Richtungen eng miteinander koppelt.

Die Autoren änderten dies zu einem block-gekoppelten Ansatz.

• Die Analogie: Anstatt dass drei Personen sich abwechseln, huddeln sie sich zusammen und lösen das Problem gleichzeitig. Indem sie die Bewegung in allen Richtungen als ein einziges großes, miteinander verbundenes Team behandeln, kann der Computer den massiven Unterschied zwischen dem dicken Kunststoff und der dünnen Luft bewältigen, ohne stecken zu bleiben.

Die Ergebnisse: Von Stunden zu Minuten

Das Team testete ihre neue Methode an zwei Szenarien:

1. Ein hundeförmiger Kanal: Ein Testfall, bei dem dicker Kunststoff in einen schmalen, gewundenen Kanal injiziert wird.

   • Alter Weg: Das Standard-Computerprogramm stürzte ab oder benötigte 7 Stunden, um ein paar Sekunden zu simulieren, da es gezwungen war, winzige Schritte zu machen.
   • Neuer Weg: Ihre neue „Team-Huddle"-Methode erledigte denselben Job in nur 16 Minuten und stürzte nicht ab, selbst wenn der Kunststoff extrem dick wurde.

2. Echte Industriemaschinen: Sie simulierten reale Einschnecken- und Doppelschnecken-Extruder (die Maschinen zur Herstellung von Kunststoffpellets).

   • Sie zeigten erfolgreich, wie der Kunststoff die Maschine füllt, wie sich der Druck aufbaut und wie die Luft verdrängt wird.
   • Sie bewiesen, dass ihre „Geisterwand"-Methode genauso gut funktioniert wie die alte, schwierige Methode des Neuformens des Gitters, jedoch viel schneller und einfacher einzurichten ist.

Was kommt als Nächstes?

Die Arbeit kommt zu dem Schluss, dass dies ein großer Schritt nach vorn für die Industrie ist. Sie schließt die Lücke zwischen akademischer Mathematik und echten Fabrikbedürfnissen. Die Autoren weisen jedoch darauf hin, dass ihr aktuelles Modell davon ausgeht, dass die Temperatur konstant bleibt (isotherm). In der Realität erzeugt das Mischen von Kunststoff Wärme, was die Dicke des Kunststoffs verändert. Das Hinzufügen von Temperatureffekten und komplexeren „dehnbaren" Kunststoffverhalten sind die nächsten Schritte für zukünftige Forschung.

Kurz gesagt: Sie entwickelten eine schnellere, stabilere Möglichkeit, das Mischen von dickem Kunststoff mit Luft in rotierenden Maschinen am Computer zu beobachten und verwandelten einen Prozess, der früher Stunden dauerte, in einen, der Minuten dauert, ohne dass jedes Mal die digitale Welt neu aufgebaut werden muss, wenn sich eine Schnecke dreht.

Technische Zusammenfassung

Technisches Fazit: Ein Volume-of-Fluid-Immersed-Boundary-Verfahren für das industrielle Polymermischen

Problemstellung
Die numerische Simulation von Polymermischprozessen stellt aufgrund des komplexen Zusammenspiels nicht-newtonscher Rheologie, mehrphasiger Freiflächendynamik und sich bewegender komplexer Geometrien erhebliche Herausforderungen dar. Industrielle Mischer, wie Einschneckenextruder (SSE) und Doppelschneckenextruder (TSE), arbeiten häufig unter Bedingungen der Teilladung, bei denen hochviskose Polymerschmelzen mit Luft koexistieren. Die Simulation dieser Szenarien erfordert die Erfassung der Phasengrenzfläche bei gleichzeitiger Handhabung rotierender fester Grenzen.

Standardansätze der Strömungsmechanik (CFD) stoßen in diesem Kontext auf zwei Haupthürden:

1. Geometrische Komplexität: Die Erzeugung körperangepasster (konformer) Netze für rotierende Schnecken mit komplexen Geometrien ist rechenintensiv und erfordert oft dynamische Netztechniken, die zu Netzverzerrungen führen können.
2. Numerische Instabilität: Der extreme Viskositätskontrast zwischen Polymerschmelzen und Luft (oft $10^5$ bis $10^8$ übersteigend) verursacht schwere numerische Instabilitäten in standardmäßigen entkoppelten Volume-of-Fluid (VOF)-Lösern. Insbesondere die explizite Behandlung des transponierten Geschwindigkeitsgradienten im viskosen Diffusionsterm setzt strenge Zeitschrittbeschränkungen ($\Delta t \propto h^2$) durch, was Simulationen hochviskoser Strömungen rechnerisch prohibitiv oder instabil macht.

Methodik
Diese Arbeit schlägt ein neuartiges numerisches Rahmenwerk vor, das einen Volume-of-Fluid (VOF)-Ansatz zur Erfassung von Grenzflächen mit einer nicht-konformen Immersed-Boundary (IB)-Methode integriert, implementiert innerhalb der OpenFOAM-Finite-Volumen-Bibliothek.

• Immersed-Boundary (IB)-Ansatz: Die Methode nutzt ein festes Hintergrundgitter, wobei feste Grenzen (z. B. rotierende Schnecken) durch triangulierte STL-Oberflächen repräsentiert werden. Ein Weighted-Least-Squares (WLS)-Interpolationsoperator wird eingesetzt, um Randbedingungen (Dirichlet für die Geschwindigkeit, Neumann für den Volumenanteil) an Zellen durchzusetzen, die von der Immersed Boundary geschnitten werden. Dies eliminiert die Notwendigkeit der Erzeugung konformer Netze und ermöglicht die Handhabung sich bewegender Kontaktlinien entlang komplexer, rotierender Geometrien.
• VOF-Formulierung: Die Grenzfläche wird durch Lösen einer Transportgleichung für den Phasenvolumenanteil ($\alpha$) verfolgt. Die Methode verwendet den MULES-Algorithmus (Multidimensional Universal Limiter with Explicit Solution), um Schärfe und Beschränktheit der Grenzfläche zu erhalten, einschließlich eines kompressiven Geschwindigkeitsterms zur Kompensation numerischer Diffusion.
• Block-gekoppeltes (BC) Schema: Um die Stabilitätsprobleme, die aus hohen Viskositätskontrasten resultieren, zu adressieren, führen die Autoren ein block-gekoppeltes Schema für die Impulsgleichung ein. Im Gegensatz zu standardmäßigen entkoppelten Lösern, die Geschwindigkeitskomponenten sequenziell behandeln, behandelt der BC-Ansatz den vollen viskosen Diffusionsterm $\nabla \cdot (\mu(\nabla \mathbf{u} + \nabla^\top \mathbf{u}))$ implizit. Dies beinhaltet die Kopplung der drei Geschwindigkeitskomponenten innerhalb eines einzigen linearen Systems. Die Diskretisierung des tangentialen Gradiententerms nutzt eine Finite-Area (FA)-Methode in Kombination mit WLS-Interpolation auf erweiterten Stencils (unter Einbeziehung von Punkt-Nachbarn statt nur Flächen-Nachbarn), angepasst aus Anwendungen der Elastizitätstheorie [11].
• Lösungsalgorithmus: Der Löser verwendet einen modifizierten PIMPLE-Algorithmus (Kombination aus SIMPLE und PISO) zur Handhabung der Druck-Geschwindigkeits-Kopplung. Die IB-Bedingungen werden in das algebraische System integriert, wobei die Matrixkoeffizienten und Quellterme modifiziert werden, um Randbedingungen auf den Immersed-Oberflächen zu erfüllen und gleichzeitig die Massenerhaltung zu gewährleisten.

Hauptbeiträge

1. Entwicklung von BC-VOF-IB-Lösern: Die Arbeit stellt den UCInterIbFoam-Löser vor, der die nicht-konforme IB-Methode mit dem block-gekoppelten VOF-Ansatz kombiniert. Dies ist die erste Implementierung einer solchen gekoppelten Methode innerhalb des OpenFOAM-10-Rahmenwerks für Zweiphasenströmungen.
2. Stabilitätsverbesserung durch Block-Kopplung: Die Arbeit zeigt, dass die vollständig implizite Behandlung der viskosen Diffusion die Stabilitätsbeschränkungen für Zeitschritte erheblich lockert. Dies ermöglicht die Simulation hochviskoser Polymerströmungen mit Viskositätsverhältnissen bis zu $10^8$, ohne dass die Simulation aufgrund von Courant-Zahl-Verletzungen abstürzt.
3. Handhabung sich bewegender Kontaktlinien: Das Rahmenwerk adressiert erfolgreich die geometrischen und algorithmischen Schwierigkeiten sich bewegender Kontaktlinien, an denen die Zweiphasengrenzfläche Immersed Boundaries schneidet, eine kritische Anforderung für teilweise gefüllte Mischvorrichtungen.

Numerische Ergebnisse
Die vorgeschlagenen Methoden wurden an zwei Kategorien von Problemen validiert und getestet:

• Benchmark-Fall (Dog-bone-Spritzgießen):

  • Ein Vergleich zwischen entkoppelten (interFoam), block-gekoppelten konformen (UCInterFoam) und block-gekoppelten nicht-konformen (UCInterIbFoam) Lösern wurde an einem Injektionsfall für hochviskose Power-law-Fluide durchgeführt.
  • Der standardmäßige entkoppelte Löser konnte innerhalb angemessener Zeiträume keine stabilen Ergebnisse liefern, bedingt durch strenge Zeitschrittbeschränkungen ($\Delta t \approx 10^{-6}$ s).
  • Die block-gekoppelten Löser erreichten stabile Lösungen mit Zeitschritten bis zu $10^{-3}$ s und reduzierten die Rechenzeit von 7 Stunden auf 16 Minuten (konform) bzw. 58 Minuten (nicht-konform).
  • Die nicht-konforme IB-Lösung zeigte eine qualitative Übereinstimmung mit der konformen Lösung, mit geringfügigen Abweichungen bei Viskositätsspitzen in der Nähe geometrischer Einschränkungen, die der nicht-konformen Natur des Gitters zugeschrieben werden.

• Industrielle Anwendungen (SSE und TSE):

  • Simulationen wurden an vereinfachten Einschneckenextrudern (SSE) und realistischen Doppelschneckenextrudern (TSE) unter sowohl vollständig gefüllten als auch teilweise gefüllten (starve-fed) Bedingungen durchgeführt.
  • Die Löser erfassten erfolgreich die Entwicklung der Freifläche, Füllverhältnisse und Druckfelder.
  • Die Ergebnisse stimmten mit physikalischen Erwartungen überein: Vollständig gefüllte Szenarien ($Q_{in} > Q_0$) zeigten negative Druckgradienten vom Einlass zum Auslass, während teilweise gefüllte Szenarien ($Q_{in} < Q_0$) Druckanstiege in der Nähe des Auslasses zeigten.
  • Die IB-Methode erwies sich als effektiv bei der Handhabung der komplexen, ineinandergreifenden Geometrie des TSE, ohne dass eine dynamische Netzregenerierung erforderlich war.

Bedeutung und Behauptungen
Die Arbeit behauptet, dass das vorgeschlagene BC-VOF-IB-Rahmenwerk einen bedeutenden Schritt darstellt, um die Lücke zwischen akademischer CFD-Forschung und industriellen Anforderungen an die Polymerverarbeitung zu schließen. Durch die Überwindung der mit hohen Viskositätskontrasten verbundenen numerischen Instabilitäten und die Eliminierung der Notwendigkeit komplexer körperangepasster Netzgenerierung für rotierende Geometrien ermöglicht die Methode physikalisch konsistente Vorhersagen von Geschwindigkeits- und Druckfeldern in industriell relevanten, teilweise gefüllten Mischvorrichtungen.

Die Autoren räumen ein, dass die aktuelle Arbeit auf isotherme Prozesse und generalisierte newtonsche Rheologie beschränkt ist. Sie stellen fest, dass zukünftige Entwicklungen thermische Effekte (Energiegleichung) und viskoelastische konstitutive Modelle einschließen müssen, um die Physik des industriellen Polymermischens vollständig zu erfassen. Das aktuelle Rahmenwerk etabliert jedoch eine rigorose Grundlage für diese Erweiterungen und bietet ein robustes Werkzeug zur Simulation von Freiflächenströmungen in komplexen, sich bewegenden Geometrien.