A numerical study on the coefficient of restitution of wet collisions

In dieser Studie wird mithilfe von SPH-Simulationen der Rückprallkoeffizient bei nassen Kollisionen untersucht, wobei eine Skalierungsgesetz identifiziert wird, das dessen Abhängigkeit von der Stokes-Zahl und einer dimensionslosen Filmdicke in zwei unterschiedlichen Regimen beschreibt.

Das Wesentliche

Das große Bild: Wenn Bälle nass werden

Stell dir vor, du wirfst einen Gummiball auf den Boden. Er prallt ab, fliegt wieder hoch und verliert dabei ein bisschen Energie. Das ist eine „trockene" Kollision. Aber was passiert, wenn du den Ball erst in eine Pfütze wirfst, bevor er den Boden berührt?

Die Forscher aus diesem Papier haben genau das untersucht: Wie stark prallt ein harter Ball ab, wenn er durch einen dünnen Film Wasser (oder Öl) fällt, bevor er auf den Boden knallt?

Das ist wichtig für viele Dinge im echten Leben: Von Schlammlawinen, die über Straßen rutschen, bis hin zu Maschinen, die Medikamente in Tabletten formen. Wenn man weiß, wie viel Energie dabei verloren geht, kann man diese Prozesse besser planen.

Die Methode: Ein digitaler Sandkasten

Da es sehr schwer und teuer ist, Tausende von Experimenten im Labor durchzuführen, haben die Wissenschaftler einen digitalen Sandkasten benutzt. Sie nannten es „Smoothed Particle Hydrodynamics" (SPH).

Stell dir das so vor:

• Statt den Ball als einen festen, glatten Stein zu sehen, haben sie ihn im Computer aus Millionen winziger, unsichtbarer Sandkörner aufgebaut.
• Das Wasser ist auch nicht eine flüssige Masse, sondern eine Ansammlung von kleinen Tropfen.
• Der Computer berechnet dann, wie sich diese Sandkörner und Tropfen gegenseitig stoßen, drücken und abprallen.

Das Tolle an ihrer Methode ist, dass sie den Ball im Inneren nicht komplett aus Sandkörnern gebaut haben, sondern nur die Hülle (wie eine leere Schale). Das spart dem Computer enorm viel Rechenzeit, genau wie wenn man einen Kuchen backt, aber nur die Kruste modelliert, weil das Innere für das Ergebnis egal ist.

Die Entdeckung: Es gibt zwei verschiedene Welten

Früher dachten die Forscher, es gäbe nur eine Regel, die bestimmt, wie hoch der Ball zurückprallt. Diese Regel hängt hauptsächlich von der Geschwindigkeit und der Zähigkeit der Flüssigkeit ab (in der Wissenschaft „Stokes-Zahl" genannt).

Aber die neuen Simulationen haben gezeigt: Das stimmt nicht immer. Es gibt eigentlich zwei verschiedene Szenarien, je nachdem, wie dick der Wasserfilm im Verhältnis zur Größe des Balls ist.

Man kann sich das wie zwei verschiedene Arten von Spielen vorstellen:

Szenario 1: Der dicke Film (Der „Schleim"-Effekt)

Stell dir vor, du wirfst einen kleinen Murmel durch eine sehr dicke Schicht Honig.

• Was passiert? Der Honig hat keine Zeit, sich schnell zu bewegen. Der Ball muss sich durch den Widerstand „durchwühlen".
• Das Ergebnis: Die Energie geht fast nur durch Reibung verloren. Die Geschwindigkeit des Balls spielt hier eine große Rolle. Je schneller er kommt, desto mehr Energie wird in Wärme umgewandelt.
• Die Regel: Hier hängt das Abprallen stark von der Geschwindigkeit ab.

Szenario 2: Der dünne Film (Der „Luftkissen"-Effekt)

Stell dir vor, du wirfst einen riesigen Basketball durch eine hauchdünne Schicht Wasser auf dem Boden.

• Was passiert? Der Ball ist so groß und schwer, dass er das Wasser nicht einfach beiseiteschieben kann. Das Wasser wird extrem schnell zusammengedrückt und baut einen enormen Druck auf (wie ein Luftkissen).
• Das Ergebnis: Es entstehen wirbelnde Strömungen und kleine Wirbel im Wasser, die Energie „fressen". Die Geschwindigkeit des Balls spielt plötzlich eine viel geringere Rolle.
• Die Regel: Hier hängt das Abprallen viel mehr davon ab, wie dünn der Film im Verhältnis zum Ball ist.

Die große Erkenntnis: Ein neuer Maßstab

Die Forscher haben eine neue Formel gefunden, die beide Welten beschreibt. Sie sagen:

„Um zu wissen, wie hoch ein Ball abprallt, musst du nicht nur wissen, wie schnell er kommt, sondern auch, wie dick der Wasserfilm im Verhältnis zum Ball ist."

Sie haben zwei verschiedene mathematische „Landkarten" erstellt:

1. Eine für kleine Bälle in dicken Filmen (hier zählt die Geschwindigkeit).
2. Eine für große Bälle in dünnen Filmen (hier zählt das Verhältnis von Film zu Ball).

Warum ist das cool?

Bisher haben Ingenieure oft nur eine einfache Regel benutzt, die bei kleinen Bällen oder dicken Filmen nicht funktioniert hat. Mit dieser neuen Erkenntnis können sie jetzt viel genauere Vorhersagen treffen.

Ein einfaches Bild zum Schluss:
Stell dir vor, du versuchst, einen Regenschirm in einer Pfütze zu öffnen.

• Wenn die Pfütze tief ist (Szenario 1), musst du gegen den Wasserwiderstand ankämpfen.
• Wenn die Pfütze nur eine dünne Schicht ist (Szenario 2), wird das Wasser unter dem Schirm extrem komprimiert und drückt ihn wieder nach oben, aber auf eine ganz andere Art.

Die Forscher haben herausgefunden, wann welche Art von Widerstand dominiert, und eine einfache Formel dafür geschrieben. Das hilft uns, alles besser zu verstehen – vom Schlamm in einer Lawine bis zur Herstellung von Tabletten in der Pharmaindustrie.

Technische Zusammenfassung

Titel: Numerische Studie zum Rückstoßkoeffizienten bei nassen Kollisionen

1. Problemstellung und Motivation

Die Abschätzung der Energiedissipation bei Kollisionen von Partikeln in Flüssigkeiten ist für die Modellierung zahlreicher granulärer Strömungsphänomene von entscheidender Bedeutung (z. B. Schlammlawinen, industrielle Granulierung, Sedimenttransport). Während der Rückstoßkoeffizient (COR, Coefficient of Restitution) bei trockenen Kollisionen oft als materialkonstant und nur schwach von der Aufprallgeschwindigkeit abhängig betrachtet wird, verändert das Vorhandensein eines interstitiellen Flüssigkeitsfilms das Bild grundlegend.

Bisherige makroskopische Modelle gehen oft davon aus, dass nasse Kollisionen allein durch die Stokes-Zahl ($St$) bestimmt werden. Experimentelle Studien (z. B. Gollwitzer et al.) zeigten jedoch, dass diese Annahme versagt, wenn sich das Verhältnis von Flüssigkeitsfilmdicke zu Partikeldurchmesser ($\delta/D$) ändert oder wenn Trägheitseffekte in der Flüssigkeitsschicht signifikant werden. Es besteht eine Lücke im Verständnis der Skalierungsgesetze, die den COR in Abhängigkeit von sowohl der Stokes-Zahl als auch der geometrischen Filmdicke beschreiben.

2. Methodik

Die Autoren nutzen Inkompressible Smoothed Particle Hydrodynamics (ISPH)-Simulationen, um den senkrechten Aufprall einer starren Kugel auf einen dünnen Flüssigkeitsfilm im Regime moderater bis hoher Weber-Zahlen zu untersuchen.

• Numerischer Ansatz:
  • Verwendung eines mesh-freien, updated-Lagrangian-Frameworks.
  • Lösung der Druck-Poisson-Gleichung mittels BiCGSTAB-Verfahren zur Sicherstellung inkompressibler Strömungen.
  • Kopplung mit der Starrkörperdynamik: Die Kugel wird als Cluster von Lagrange-Partikeln modelliert, wobei nur der Bereich in der Nähe der Flüssigkeit diskretisiert wird (eine sphärische Schale), um den Rechenaufwand zu minimieren, während die physikalischen Eigenschaften (Masse, Trägheitsmoment) der Vollkugel beibehalten werden.
  • Die Starrkörperbewegung (Translation und Rotation) wird über Impulsbilanzen und die Euler-Gleichungen gelöst.
• Kollisionsmodellierung:
  • Der Aufprall auf den festen Untergrund wird durch einen trockenen Rückstoßkoeffizienten ($e_r \approx 0.98$) modelliert.
  • Die Benetzungseffekte und die Energiedissipation entstehen durch die Wechselwirkung mit der Flüssigkeitsschicht (viskose Kräfte, Druckaufbau, Formwiderstand).
  • Der COR wird als Verhältnis der Geschwindigkeit beim Verlassen des Films zur Geschwindigkeit beim Eintritt definiert.
• Validierung:
  • Die Simulationen wurden gegen experimentelle Daten von Gollwitzer et al. (Aufprall von Glaskugeln auf Silikonöl) validiert.
  • Eine Konvergenzstudie zeigte, dass bei einer Auflösung von ca. 178.000 Flüssigkeitspartikeln eine Abweichung von unter 1 % gegenüber den Experimenten erreicht wird.

3. Wichtige Beiträge und Ergebnisse

A. Identifikation neuer dimensionsloser Parameter
Die Studie zeigt, dass die Stokes-Zahl allein nicht ausreicht, um den COR zu charakterisieren. Als zweiten entscheidenden Parameter wurde die dimensionslose Filmdicke eingeführt:
$$\gamma = \frac{\delta}{D}$$
wobei $\delta$ die Filmdicke und $D$ der Partikeldurchmesser ist.

B. Entdeckung zweier Skalierungsregime
Durch systematische Variation von Aufprallgeschwindigkeit, Partikelgröße, Viskosität und Filmdicke wurden zwei distincte Regime identifiziert, die durch unterschiedliche Potenzgesetz-Exponenten gekennzeichnet sind. Die Beziehung wird als $\text{COR} \sim St^a \gamma^b$ modelliert:

1. Regime R1 (Niedrige $\gamma$, größere Partikel):

   • Hier zeigt der COR eine starke Abhängigkeit sowohl von der Stokes-Zahl als auch von der dimensionslosen Filmdicke.
   • Skalierungsgesetz: $\text{COR} \propto St^{0.17} \gamma^{-0.16}$.
   • Physikalisch dominiert hier der viskose Widerstand und der lineare Druckaufbau.

2. Regime R2 (Hohe $\gamma$, kleinere Partikel):

   • Die Abhängigkeit von der Stokes-Zahl ist hier vernachlässigbar schwach ($St^{0.01}$), während die geometrische Dicke ($\gamma$) den dominierenden Einfluss hat ($\gamma^{-0.22}$).
   • Skalierungsgesetz: $\text{COR} \propto St^{0.01} \gamma^{-0.22}$.
   • In diesem Regime spielen sekundäre Strömungsphänomene (Wirbelbildung) eine größere Rolle bei der Energiedissipation.

C. Physikalische Erklärung der Regime
Die Analyse der kinetischen Energie im Fluid während des Aufpralls liefert eine Erklärung für die zwei Regime:

• Im Regime R1 dissipiert die kinetische Energie des Fluids monoton nach dem Kontakt mit dem Boden. Dies deutet auf einen primären Strömungsmechanismus hin, der stark von der Viskosität und damit von $St$ abhängt.
• Im Regime R2 tritt ein sekundärer Peak in der kinetischen Energie auf, was auf die Bildung von Wirbeln und komplexeren sekundären Strömungen hindeutet. Diese Mechanismen machen den COR weniger empfindlich gegenüber der Stokes-Zahl und stärker abhängig von der Geometrie ($\gamma$).

4. Bedeutung und Fazit

Die Studie liefert einen wichtigen Fortschritt in der Modellierung nasser Partikelkollisionen:

• Skalierungsgesetz: Sie etabliert ein neues, zweiparametriges Skalierungsgesetz, das die Grenzen der bisherigen, rein stokes-basierten Modelle überwindet.
• Anwendbarkeit: Die entwickelten Modelle und Algorithmen sind direkt in diskrete Elementmethoden (DEM) und andere statistische Simulationsansätze integrierbar, um die Energiedissipation in granularen Strömungen mit Flüssigkeitsbrücken präziser vorherzusagen.
• Methodischer Beitrag: Die vorgestellte effiziente Kopplung von ISPH mit Starrkörperdynamik (unter Verwendung einer sphärischen Schale statt einer Vollkugel) ermöglicht hochauflösende Simulationen bei deutlich reduzierten Rechenkosten.

Zusammenfassend demonstriert die Arbeit, dass das Verhalten nasser Kollisionen nicht durch einen einzigen Parameter bestimmt wird, sondern durch das Zusammenspiel von Trägheit (Stokes-Zahl) und Geometrie (dimensionslose Filmdicke), wobei der Übergang zwischen den Regimen durch den Wechsel der dominierenden Dissipationsmechanismen (viskose Reibung vs. sekundäre Strömungen) erklärt werden kann.