Modelling of pressure drop in periodic square-bar packed beds

Diese Studie untersucht numerisch und experimentell den Druckabfall in periodischen Packungen aus quadratischen Stäben, wobei sie zeigt, dass sich die Strömungsgeometrie je nach Rotationswinkel in kanal- oder gitterartige Konfigurationen unterteilt und die Verwendung eines winkelabhängigen äquivalenten Moduldurchmessers eine erfolgreiche Vorhersage des Reibungsfaktors mittels der Ergun-Korrelation ermöglicht.

Das Wesentliche

Das große Rätsel: Wie fließt Wasser durch einen Steinhaufen?

Stell dir vor, du hast einen riesigen Haufen aus Steinen (oder in diesem Fall aus quadratischen Stäben), durch den Wasser fließen muss. Ingenieure müssen genau wissen, wie viel Widerstand das Wasser spürt, um zum Beispiel Filter zu bauen oder chemische Reaktoren zu optimieren.

Das Problem ist: Die meisten bisherigen Studien haben sich nur mit runden Kugeln beschäftigt. Das ist wie wenn man lernt, wie Autos fahren, indem man nur Rennwagen betrachtet. Aber in der echten Welt sind die "Steine" oft eckig, unregelmäßig oder wie hier: quadratische Stäbe, die wie ein Puzzle gestapelt sind.

Das Experiment: Ein sich drehendes Turm-Modell

Die Forscher haben sich ein spezielles Experiment ausgedacht, das wie ein Turm aus drehbaren Tellern aussieht:

• Jeder "Teller" (ein Modul) hat darin quadratische Stäbe.
• Man kann jeden Teller im Verhältnis zum Teller darunter ein bisschen drehen.
• Durch diese Drehung entstehen völlig neue Wege für das Wasser. Mal sind die Wege gerade wie ein Autobahn-Tunnel (wenn die Stäbe fast übereinander liegen), mal sind sie wie ein Labyrinth oder ein Gitter, in dem das Wasser hin und her springen muss.

Sie haben das Wasser durch diesen Turm gepumpt und dabei die Geschwindigkeit variiert – von ganz langsam (wie ein schlendernder Spaziergänger) bis schnell (wie ein Sprinter).

Die wichtigsten Entdeckungen (in Alltagssprache)

1. Die zwei Welten: Autobahn vs. Labyrinth
Je nachdem, wie stark man die Teller dreht, ändert sich das Verhalten des Wassers komplett:

• Wenig Drehung (bis ca. 10°): Die Stäbe liegen fast genau übereinander. Es entstehen lange, gerade Röhren. Das Wasser fließt hier wie auf einer Autobahn – schnell und ohne viel Widerstand. Die Forscher nennen das "kanalartig".
• Viel Drehung (ab ca. 15°): Die Stäbe verschieben sich. Die Wege werden zu einem Labyrinth. Das Wasser muss ständig abbiegen, wird gebremst und wirbelt herum. Das nennt man "gitterartig". Hier ist der Widerstand viel höher.

2. Der "schlimmste" Winkel hängt von der Geschwindigkeit ab
Das ist das Überraschendste: Welcher Winkel den größten Widerstand erzeugt, hängt davon ab, wie schnell das Wasser fließt!

• Langsames Wasser (wie Honig): Der größte Widerstand entsteht bei 25°. Warum? Weil sich die Öffnungen an dieser Stelle so verengen, dass das Wasser wie durch einen sehr engen Schlund gequetscht wird.
• Schnelles Wasser: Bei hoher Geschwindigkeit ist 60° der schlimmste Winkel. Hier passiert etwas anderes: Das Wasser schießt in die Öffnungen, prallt ab und wirbelt wild (wie ein Auto, das zu schnell um eine enge Kurve fährt und ins Schleudern gerät). Der Widerstand kommt hier nicht vom "Quetschen", sondern vom "Wirbeln".

3. Die neue Formel für Ingenieure
Bisher haben Ingenieure Formeln benutzt, die für runde Kugeln gemacht waren. Die haben bei diesen eckigen Stäben oft daneben gelegen.
Die Forscher haben eine neue Methode entwickelt: Sie betrachten nicht den einzelnen Stab, sondern das ganze Modul als eine Art "Super-Teilchen". Wenn man diese neue Größe in die alten Formeln einsetzt, passen die Vorhersagen plötzlich viel besser. Es ist, als würde man statt der Größe eines einzelnen Ziegelsteins die Größe der ganzen Mauer messen, um den Luftzug vorherzusagen.

4. Wann wird es turbulent?
Normalerweise denkt man, Wasser wird erst bei sehr hoher Geschwindigkeit turbulent (unruhig). Aber bei diesen speziellen Gittern beginnt das "Unwetter" schon viel früher (bei einer bestimmten Geschwindigkeit), weil die Ecken und Kanten das Wasser schneller durcheinanderbringen als runde Kugeln.

Warum ist das wichtig?

Diese Studie ist wie ein Bauplan für die Zukunft. Wenn Ingenieure wissen, wie sich das Wasser in solchen komplexen, eckigen Strukturen verhält, können sie:

• Bessere Filter für die Umwelttechnik bauen.
• Effizientere Reaktoren für die Chemieindustrie entwerfen.
• Energie speichern, indem sie den Wärmefluss besser verstehen.

Kurz gesagt: Die Forscher haben bewiesen, dass die Form und Anordnung der "Steine" im Filter genauso wichtig ist wie ihre Größe. Und sie haben gezeigt, dass man mit ein paar cleveren Drehungen am Design den Fluss des Wassers genau steuern kann – entweder wie auf einer Autobahn oder wie in einem wilden Labyrinth.

Technische Zusammenfassung

Titel: Modellierung des Druckabfalls in periodischen Packungen aus quadratischen Stäben

1. Problemstellung und Motivation

Das Verständnis von Strömungen durch poröse Medien mit komplexen Geometrien ist entscheidend für die Optimierung von Festbettreaktoren in der chemischen und verfahrenstechnischen Industrie. Die meisten bestehenden Studien und Modelle konzentrieren sich jedoch auf kugelförmige Packungen. Dies führt zu einem Mangel an genauen Modellen für unregelmäßige Zwischenraumgeometrien (interstitielle Geometrien), wie sie bei nicht-kugelförmigen Partikeln (z. B. Zylinder, Polyeder) auftreten.

Das vorliegende Ziel dieser Studie ist es, den Einfluss der Geometrievariation auf den induzierten Strömungswiderstand (Drag) und die Permeabilität in einem dichten Festbett aus quadratischen Stäben zu quantifizieren. Ein zentrales Problem ist, dass Standard-Korrelationen wie die Ergun-Gleichung ursprünglich nur für Kugeln entwickelt wurden und ihre Übertragbarkeit auf komplexe, nicht-sphärische Anordnungen bei konstanten makroskopischen Eigenschaften (z. B. Porosität) oft unzureichend ist.

2. Methodik

Die Studie basiert auf numerischen Simulationen mittels Direct Numerical Simulation (DNS) mit dem Open-Source-CFD-Code OpenFOAM.

• Geometrie: Untersucht wurde eine modulare Festbettanordnung, bestehend aus scheibenförmigen Modulen (Dicke $B = 10$ mm), die quadratische Stäbe enthalten. Die Module sind vertikal gestapelt und können relativ zueinander um einen Winkel $\alpha$ gedreht werden ($0^\circ$ bis $90^\circ$ in Schritten von $5^\circ$). Dies ermöglicht die Erzeugung verschiedener Zwischenraumgeometrien bei konstanter Porosität ($\phi = 0.322$) und Querschnittsfläche.
• Strömungsbedingungen: Die Simulationen decken einen Reynolds-Zahlen-Bereich von $Re_p = 0.1$ bis $200$ ab (basierend auf der Stabgröße $B$ und der intrinsischen Durchschnittsgeschwindigkeit). Dies umfasst sowohl den viskosen (Kriechströmung) als auch den trägheitsdominierten Bereich.
• Validierung: Die numerischen Ergebnisse wurden für ausgewählte Fälle ($\alpha = 30^\circ$, $Re_p = 100, 200$) mit experimentellen Daten aus Particle Image Velocimetry (PIV) Messungen (Velten et al., 2026) validiert.
• Numerisches Vorgehen: Es wurden sowohl stationäre als auch instationäre Solver verwendet. Eine detaillierte Gitterunabhängigkeitsstudie bestätigte, dass ein Gitter mit ca. 1 Million Zellen pro Modul für die meisten Fälle ausreichend ist. Periodische Randbedingungen wurden angewendet, um unendlich wiederholte Strukturen abzubilden.

3. Wichtige Beiträge und Erkenntnisse

A. Klassifizierung der Geometrien:
Die Studie identifiziert zwei fundamentale Strömungsregime, die durch den Rotationswinkel $\alpha$ bestimmt werden:

1. Kanal-ähnlich (Channel-like): Für kleine Winkel ($\alpha \le 10^\circ$). Die Stäbe sind weitgehend ausgerichtet, was zu geradlinigen, wenig behinderten Strömungspfaden führt.
2. Gitter-ähnlich (Lattice-like): Für größere Winkel ($\alpha \ge 15^\circ$). Die Zwischenräume sind mehrfach verbunden, was zu komplexeren, verzweigten Strömungspfaden führt.
   • Klassifikationskriterien: Der Übergang lässt sich quantitativ durch das Verhältnis von Bettdurchmesser zu äquivalentem Partikeldurchmesser ($D/d_{eq} \approx 3$) und die hydraulische Tortuosität ($\tau \approx 1.1$) definieren.

B. Reibungsfaktor und Widerstandscharakteristik:

• Der Reibungsfaktor $f_p$ hängt stark vom Winkel und dem Reynolds-Zahl-Bereich ab.
• Viskoses Regime: Das Maximum des Reibungsfaktors liegt bei $\alpha = 25^\circ$. Dies wird auf eine starke Verengung der Strömungspfade und lokale Scherspannungen zurückgeführt.
• Trägheitsregime: Das Maximum verschiebt sich zu $\alpha = 60^\circ$. Hier dominieren Formwiderstände durch Strömungsablösung und Wirbelbildung in Kontraktions-Expansions-Sequenzen.
• Bei $\alpha = 90^\circ$ ist der Widerstand trotz der "Gitter"-Klassifizierung relativ gering, da die orthogonale Anordnung minimale Behinderung bietet.

C. Permeabilitätsmodellierung und charakteristische Längenskalen:

• Die Wahl der charakteristischen Längenskala ist entscheidend für die Korrelation mit der Ergun-Gleichung.
• Der modulbasierte äquivalente Durchmesser ($d_{eq}$), der die winkelabhängige benetzte Oberfläche berücksichtigt, führt dazu, dass die Reibungsfaktoren für gitterähnliche Geometrien auf die Ergun-Korrelation kollabieren.
• Im Gegensatz dazu versagt der konstante Durchmesser eines einzelnen Stabs ($d_{sb}$), da er die geometrischen Änderungen der benetzten Oberfläche nicht erfasst und den Reibungsfaktor systematisch überschätzt.
• Die Blake-Kozeny-Beziehung basierend auf $d_{eq}$ liefert für gitterähnliche Geometrien die beste Vorhersage der Permeabilität (mittlere absolute prozentuale Abweichung MAPE = 12,3 %).

D. Übergang zur Trägheitsströmung:

• Der Übergang vom Darcy-Regime (viskos) zum Forchheimer-Regime (träge) erfolgt bei den untersuchten Geometrien typischerweise bei $Re_p \approx 7,5$ (definiert als Punkt, an dem der träge Anteil 5 % des Gesamtwiderstands ausmacht).
• Bei $\alpha = 55^\circ$ tritt dieser Übergang früher auf ($Re_p \approx 5$), was auf eine empfindliche Topologie der Zwischenräume und frühe lokale Rezirkulation hindeutet.
• Kanal-ähnliche Geometrien ($\alpha \le 10^\circ$) zeigen eine verzögerte Transition, da die ausgerichteten Pfade die Strömungsablösung unterdrücken.

4. Signifikanz und Fazit

Diese Studie liefert einen wichtigen Beitrag zum Verständnis von Strömungen in nicht-sphärischen porösen Medien, die in der Praxis häufiger vorkommen als ideale Kugelpackungen.

• Validiertes Framework: Die Kombination aus DNS und PIV-Validierung bietet ein zuverlässiges Werkzeug zur Analyse komplexer periodischer Strukturen.
• Geometrie-Wirkungs-Beziehung: Es wurde gezeigt, dass makroskopische Parameter wie Porosität allein nicht ausreichen, um den Strömungswiderstand vorherzusagen; die spezifische Anordnung (Winkel) und die daraus resultierende Topologie der Zwischenräume sind entscheidend.
• Verbesserte Modelle: Die Einführung des winkelabhängigen äquivalenten Durchmessers ($d_{eq}$) ermöglicht eine deutlich genauere Anwendung klassischer Korrelationen (wie Ergun) auf komplexe Geometrien als bisherige Ansätze mit konstanten Partikeldurchmessern.
• Zukünftige Anwendungen: Die Ergebnisse sind relevant für das Design von Festbettreaktoren, Wärmespeichersystemen und Filtern, insbesondere wenn nicht-kugelförmige Materialien verwendet werden. Die Studie unterstreicht die Notwendigkeit, geometrieangepasste Modelle für den Forchheimer-Koeffizienten und die Permeabilität zu entwickeln, die über den laminaren Bereich hinausgehen.