Wir freuen uns über ein neues CastNet-Release mit zahlreichen neuen Funktionen. Ein Schwerpunkt: OpenFOAM wird für Windows-Anwender deutlich einfacher nutzbar. Wir liefern jetzt native Windows-Installer für die Keysight/OpenCFD-Versionen (v2512, v2606) und die Foundation-Versionen (v12, v13) – integriert in eine einsatzbereite Basisumgebung mit ParaView und nützlichen OpenFOAM-Add-ons. BlueCFD wird weiterhin unterstützt.

Auch Strukturmechanik ist abgedeckt – über den nichtlinearen FE-Löser CalculiX. Als Highlight liefern wir zusätzlich preCICE mit Windows-Installer aus und ermöglichen so FSI-Kopplung zwischen OpenFOAM und CalculiX. Damit steht eine komplette Suite zur Lösung von CAE-Aufgaben in einem durchgängigen Workflow bereit.

Das Video zeigt Highlights am Beispiel eines Wärmeübertragers: unseren hybriden Vernetzungsansatz (Hex-Core mit einfach erstellbaren Extrudes über mehrere Regionszonen) und die Multi-Region-Unterstützung für CHT. Weitere Neuerungen: deutlich erweitertes RunGui-Monitoring, RBM-Unterstützung für Mesh-Morphing und Overset-Meshes sowie neue Meshing-Funktionen.

Fast jedes Zerstäubungsproblem hat eine Lücke in den Zeitskalen: Der primäre Strahlzerfall geschieht in wenigen Millisekunden und lässt sich nur mit adaptiver Verfeinerung, winzigen Zeitschritten und zig Millionen Zellen auflösen. Entscheidend für das Produkt ist aber meist, was weiter stromabwärts und über viel längere Zeiträume passiert – etwa beim Lackieren und Beschichten, bei der Metallzerstäubung (Tropfen kühlen ab und erstarren zu Pulver) oder bei Wandfilmen und Sprühchemie wie AdBlue/SCR.

Da sich der Strahlzerfall nicht über die gesamte Prozessdauer auflösen lässt, trennen wir beide Bereiche: Kurz hinter unserem VoF-zu-Lagrange-Übergang platzieren wir eine Sampling-Ebene, die jeden durchtretenden Tropfen als einzelnes Ereignis erfasst – Größe, Geschwindigkeit, Richtung und lokaler Massenstrom – über sieben bis acht Pulse, sodass die natürliche kHz-Pulsation erhalten bleibt.

Im Video vergleichen wir verschiedene Methoden, die Partikeldaten in ein rein Lagrange'sches Setup wieder einzuspeisen: voll aufgelöste transiente Injektion, mit oder ohne adaptive Netzverfeinerung (AMR) der Tropfen, sowie die Umwandlung in einen OpenFOAM-Cone-Injector. So lassen sich die durch Vereinfachungen entstehenden Fehler gezielt eingrenzen.

Hier das dritte Beispiel unseres Membran-FSI-Lösers, direkt in das Rigid-Body-Motion-Framework (RBM) von OpenFOAM integriert – diesmal eine VoF-Anwendung (isoAdvector). Angelehnt an OpenFOAMs floating-object-Tutorial modellieren wir eine Box aus Baffle-Wänden, von denen zwei elastisch sind. Die Box ist in Wasser getaucht und innen mit Luft gefüllt: In den ersten Sekunden stellt sich über eine relaxierte FSI-Kopplung das Gleichgewicht zwischen Auftrieb und Wandverformung ein, danach wird sie Wellen ausgesetzt.

Dank der eng integrierten Kopplung bleibt die Gesamtlaufzeit mit rund 2 Stunden auf 16 Kernen (ca. 600.000 Zellen) kurz. Ein zentraler Vorteil: Die Zeitschrittweite lässt sich dynamisch an Kriterien wie die alpha-Courant-Zahl, die Bewegungsgeschwindigkeit der Box oder die Verformungsrate der Membranen anpassen.

Zahlreiche Methoden zur Stabilisierung sind direkt im Löser integriert und machen ihn besonders effizient für Membrananwendungen: verschiedene Relaxationsschemata (Aitken, IQN-ILS), zeitliche Filter und räumliche Glättung des Drucks, Vorspann-Techniken sowie diverse Dämpfungsansätze wie Rayleigh-Dämpfung.

Wir haben den klassischen Mok-Benchmark gerechnet: einen Strömungskanal, der von einer elastischen Membran begrenzt wird, die sich unter dem Unterdruck einer zeitlich veränderlichen Einström-Geschwindigkeit verformt – mit ausgezeichneter Übereinstimmung zur veröffentlichten Referenz.

Unser FSI-Löser wurde für elastische Membranen entwickelt, die an starre Körper gekoppelt sind (über OpenFOAMs rigidBodyMotion-Bibliothek). Wie Mok zeigt, kann dieselbe Membran aber ebenso gut die äußere Berandung des Strömungsgebiets selbst sein.

Warum wir den Löser vollständig in OpenFOAM integriert halten: Stabilisierungs-Werkzeuge sitzen genau dort, wo sie gebraucht werden – im Löser selbst; das ist entscheidend für Membranen, die keinen Druck aufnehmen können und beim Vorzeichenwechsel sofort zu falten beginnen. Es gibt keinen externen Kopplungs-Overhead, und die implizite FSI (IQN-ILS) mit Quasi-Newton-Beschleunigung bewältigt starke Added-Mass-Fälle zuverlässig.

Gekoppelte FSI ohne OpenFOAM zu verlassen: Für geeignete Fälle integriert DHCAE die Strukturmechanik der Membran direkt in den Löser – mit OpenFOAMs eigener Infrastruktur. Für gekoppelte Multiphysik sind preCICE und solids4Foam etablierte, leistungsfähige Optionen; in bestimmten Szenarien zahlt es sich jedoch aus, die Strukturmechanik direkt einzubetten: wenn der strukturelle Aufwand gegenüber der Strömung klein ist, eine zuverlässige Elementformulierung vorliegt und spezielle Stabilisierungs- oder Beschleunigungstechniken für robuste Konvergenz zählen.

Genau das haben wir gebaut: Membranelemente im Rigid-Body-Framework (RBM) von OpenFOAM. Bewegte Körper können starre Teile und elastische Membranbereiche kombinieren. Eingebaut für robuste FSI sind u. a. ein Tension-Field-Wrinkling-Modell, statische Vorspannung (pOffset), räumliche und zeitliche Glättung der Membran-Drucklast, FSI-gesteuerte adaptive Zeitschrittweite, Anlauframpen für einen stoßfreien Start sowie Rayleigh-Dämpfung.

Die komplette FSI läuft als äußere Schleife um PIMPLE, mit Aitken- oder IQN-ILS-Relaxation – geeignet für transiente und quasi-statische Anwendungen. Einphasen- und VoF-Varianten sind verfügbar. Validiert gegen die Hencky-Membran und den Mok-Benchmark. Typische Anwendungen: Wind- und Strömungslasten auf Segeln und Zugmembranen, aufblasbare Strukturen, hydraulische Membranventile, Fallschirme sowie mobile Hochwasserschutz-Barrieren.

Wir nutzen adaptive Netzverfeinerung (AMR) in OpenFOAM intensiv innerhalb unseres hybriden VoF-zu-Lagrange-Übergangslösers, um Tropfen und Ligamente während der Zerstäubung aufzulösen, bevor sie in Lagrange'sche Parcels umgewandelt werden. Danach wird der Spray auf einem bewusst gröberen Netz propagiert, um den Rechenaufwand beherrschbar zu halten.

Daraus entsteht ein Zellgrößen-Dilemma im Sprühbereich: Das Netz muss grob genug sein, damit die Punktpartikel-Annahme gilt (Zellen größer als die Tropfen), aber fein genug, damit die LES-aufgelöste Gasströmung die Nachläufe und Querströmungen zwischen großen Tropfen tatsächlich trägt – denn diese treiben Einmischung und Impulsaustausch der kleinen Partikel.

Deshalb haben wir unser AMR so erweitert, dass es auch innerhalb des Lagrange-Sprays verfeinert – gesteuert durch ein Partikel-Einflussfeld statt durch den VoF-Indikator. Die Verfeinerung reicht in den Sprühkegel hinein, bleibt aber deutlich oberhalb der Tropfenskala und erhält so die Gültigkeit des Punktpartikel-Ansatzes – ohne die Kosten einer Verfeinerung des gesamten Gebiets.

Industrielle Schlauchfilter sind die Arbeitspferde der Staubabscheidung – doch das Beladungsmuster im Inneren eines laufenden Filterhauses ist schwer messbar und notorisch ungleichmäßig. Wohin wandert der Staub tatsächlich? Welche Schläuche tragen die Last? Wie verteilt sich der Luftstrom um, während sich das Medium füllt?

Um das zu untersuchen, haben wir einen vollständigen Beladungszyklus eines Filterhauses mit 16 Schläuchen mit Lagrange'scher Partikelverfolgung simuliert, gekoppelt an ein Filtermedium, dessen lokaler Widerstand mit der Staublast steigt. Die iterative Kopplung zwischen Partikelablagerung und Strömungsfeld übernimmt unser Filter-Löser auf Basis von OpenFOAM v2512.

Einige Beobachtungen: Der Einström-Jet treibt Staub bevorzugt auf die zentralen Schläuche – und dort zuerst auf die unteren Bereiche. Die Schläuche an der Rückwand sammeln über die Zeit die höchste Gesamtlast. Der Luftstrom verlagert sich fortlaufend zu den weniger beladenen Schläuchen. Der Nutzen: die Strömungsverteilung optimieren, um ungeplante Stillstände zu vermeiden, die installierte Filterfläche effizient nutzen und Reinigungszyklen an die tatsächliche Last statt an feste Intervalle koppeln.

Credits: Die Filterhaus-Geometrie ist ein Demonstrationsmodell von wahyu suryanto auf GrabCAD – danke für die freie Bereitstellung.