Continuum modeling of Soft Glassy Materials under shear

Dieser Artikel stellt ein räumlich aufgelöstes, nicht-lokales Kontinuumsmodell auf Fluiditätsbasis vor, das das scherverursachte Fließen und die damit verbundenen komplexen Phänomene wie Spannungsüberschwingen und Scherbänderung in weichen glasartigen Materialien quantitativ beschreibt.

Das Wesentliche

Wenn Festes flüssig wird: Eine Reise durch den „weichen Glaskörper"

Stellen Sie sich vor, Sie haben ein Glas voll mit winzigen, weichen Kügelchen – vielleicht wie winzige Wasserbällchen oder kleine Gummibärchen, die so dicht gepackt sind, dass sie sich gar nicht bewegen können. Wenn Sie das Glas ruhig stehen lassen, verhält sich diese Masse wie ein Feststoff. Sie können sie mit dem Finger berühren, und sie gibt nicht nach.

Aber wenn Sie nun anfangen, das Glas zu schütteln oder die Masse zu rühren (das nennt man in der Physik „Scherung"), passiert etwas Magisches: Die Masse wird plötzlich flüssig. Dieser Übergang von fest zu flüssig ist das Herzstück dieser Forschungsarbeit.

Die Wissenschaftler in diesem Papier haben sich gefragt: Wie genau passiert das? Und wie schnell?

1. Das Problem: Der „Stress-Ausbruch"

Wenn Sie so eine Masse (wie Zahnpasta, Ketchup oder spezielle Gele) plötzlich schnell rühren, passiert Folgendes:

• Zuerst widersteht sie stark. Der Widerstand (die Spannung) steigt schnell an.
• Dann erreicht sie einen Höhepunkt (einen „Stress-Overshoot"). Es ist, als würde man versuchen, einen schweren Koffer über einen rutschigen Boden zu ziehen: Man muss erst richtig stark ziehen, bis er sich löst.
• Sobald er sich löst, wird es plötzlich viel leichter, ihn zu ziehen. Die Masse fließt.

Das Interessante ist: Dieser Übergang passiert nicht immer überall gleichzeitig. Oft bricht die Masse an einer Stelle auf (wie ein Riss in Eis), und diese flüssige Zone wächst langsam durch den Rest der Masse, bis alles fließt. Das nennt man „Scherbanden" (Schubzonen).

2. Die Lösung: Ein neues Modell („Fluidity")

Die Forscher haben ein neues mathematisches Modell entwickelt, um dieses Verhalten vorherzusagen. Statt jedes einzelne Kügelchen zu berechnen (was unmöglich wäre, da es Milliarden davon sind), betrachten sie die Masse als eine Art „flüssiges Feststoff-Gemisch".

Sie führen eine neue Größe ein, die sie „Fluidity" (Fließfähigkeit) nennen.

• Stellen Sie sich die Fluidity wie eine Art „Unruhe" oder „Energie" vor.
• In den festen Bereichen ist die Unruhe null.
• In den flüssigen Bereichen ist die Unruhe hoch.

Das Geniale an ihrem Modell ist, dass sie berücksichtigen, wie sich diese Unruhe von einem Ort zum anderen ausbreitet. Wenn ein Kügelchen sich bewegt, gibt es einen kleinen „Stoß" an seine Nachbarn weiter. Das ist wie eine Welle im Stadion: Wenn eine Person aufsteht (sich bewegt), bringt sie die nächste Person dazu, aufzustehen. Das Modell beschreibt genau, wie schnell sich diese Welle durch den ganzen Raum ausbreitet.

3. Was das Modell erklärt

Mit diesem Modell konnten die Forscher zwei wichtige Dinge erklären, die man im Labor beobachtet:

• Wie schnell rührt man, bestimmt den Ausbruch: Je schneller man rührt, desto höher ist der Widerstand, bevor die Masse fließt. Das Modell sagt genau voraus, wie dieser Peak aussieht.
• Wie lange dauert es, bis alles fließt? Wenn man langsam rührt, dauert es lange, bis die flüssige Zone die ganze Masse durchdringt. Das Modell zeigt, dass diese Zeit nicht linear abnimmt, sondern einem ganz bestimmten mathematischen Gesetz folgt.

4. Der „Gleit-Effekt" an den Wänden

Ein weiterer wichtiger Punkt ist, was an den Wänden passiert. Wenn die Masse gegen eine glatte Wand gedrückt wird, kann sie dort leicht rutschen (wie auf einer Eisbahn). Wenn die Wand rau ist, bleibt sie haften.
Die Forscher haben ihr Modell so erweitert, dass es diesen Gleit-Effekt berücksichtigt. Sie haben entdeckt, dass dieser Effekt die Art und Weise verändert, wie die Masse fließt, besonders bei sehr langsamen Bewegungen. Es ist, als würde man versuchen, einen Teppich auf einem glatten Parkett zu verschieben – er rutscht leicht, aber auf einem rauen Teppich bleibt er stehen.

5. Warum ist das wichtig?

Dieses Modell ist wie eine Landkarte für Ingenieure.

• In der Industrie: Wenn man Zahnpasta, Farben, Klebstoffe oder sogar 3D-Druck-Materialien herstellt, muss man genau wissen, wie sie sich verhalten, wenn man sie durch Rohre pumpt oder auf Oberflächen aufträgt.
• Die Vorhersage: Mit diesem Modell kann man berechnen, ob ein Material beim Pumpen plötzlich „platzen" wird (zu schnell fließen) oder ob es stecken bleibt. Man kann also den Prozess optimieren, bevor man überhaupt ein Experiment macht.

Zusammenfassung in einem Satz

Die Forscher haben eine Art „Wettervorhersage" für weiche, zähe Materialien entwickelt, die erklärt, wie diese Materialien unter Druck von fest zu flüssig wechseln, indem sie beschreiben, wie sich kleine Bewegungen durch das Material ausbreiten – ähnlich wie eine Welle, die durch eine Menschenmenge läuft.

Technische Zusammenfassung

1. Problemstellung und Hintergrund

Weiche glasartige Materialien (Soft Glassy Materials, SGMs) wie Carbopol-Mikrogel, Emulsionen oder dichte kolloidale Suspensionen zeigen im Ruhezustand feste Eigenschaften, können jedoch unter ausreichend großer Scherung in einen flüssigen Zustand übergehen. Dieser Übergang ist durch komplexe Phänomene gekennzeichnet:

• Stress-Overshoot: Unter konstanter Scherrate baut sich die Spannung zunächst linear auf (elastisches Verhalten), erreicht ein Maximum und relaxiert dann auf einen stationären Wert.
• Räumliche Heterogenität: Der Fließübergang ist oft nicht homogen. Es entstehen Scherbänder (Shear Bands), bei denen fluide Bereiche mit erstarrten, feststoffähnlichen Bereichen koexistieren.
• Modellierungsdefizite: Bestehende Modelle auf makroskopischer Ebene (Kontinuumsmechanik) erfassen oft nicht die nicht-lokalen Effekte, die für die Bildung und Dynamik dieser Scherbänder sowie für die Abhängigkeit des Stress-Overshoots von der Scherrate verantwortlich sind. Zudem fehlt oft eine konsistente Beschreibung von Randbedingungen, insbesondere bei weichen Partikeln, die an festen Wänden gleiten (Elasto-Hydrodynamik).

2. Methodik: Das nicht-lokale Fluiditäts-Modell

Die Autoren erweitern einen bereits existierenden Kontinuumsansatz, der auf dem Konzept der Fluidität ($f$) basiert. Die Fluidität wird als Ordnungsparameter definiert, der die lokale Rate plastischer Ereignisse repräsentiert.

Die Modellierung stützt sich auf drei Hauptkomponenten:

1. Dynamikgleichung der Fluidität: Die zeitliche Entwicklung der Fluidität wird durch eine Gleichung beschrieben, die als Gradientenfluss eines freien Energie-Funktionals $F[f]$ formuliert ist:
   $$\frac{\partial f}{\partial \tilde{t}} = f \left( \xi^2 \Delta f + m f - f^{3/2} \right)$$
   Hierbei ist $\xi$ eine Kooperativitäts-Längenskala, die nicht-lokale Effekte (wie die Ausbreitung plastischer Ereignisse in die Nachbarschaft) beschreibt. Der Term $m$ hängt von der Spannung $\Sigma$ ab und bestimmt, ob das Material fließt ($m^2 > 0$) oder fest ist.
2. Maxwell-Gleichung für die Spannung: Die Spannung $\Sigma$ entwickelt sich gemäß einer Maxwell-Beziehung, die elastische und plastische Dehnung kombiniert:
   $$\frac{d\Sigma}{d\tilde{t}} = \frac{1}{\tau} \left[ \dot{\Gamma} - \langle f \rangle \Sigma \right]$$
   wobei $\langle f \rangle$ der räumliche Mittelwert der Fluidität ist.
3. Randbedingungen: Es werden spezifische Randbedingungen für die Fluidität an den Wänden definiert (z. B. feste Fluidität $f_w = m^2$ an der bewegten Wand), um die Keimbildung des flüssigen Phasenzustands zu steuern.

Erweiterung um Elasto-Hydrodynamik (EHD):
In einem zweiten Schritt wird das Modell um EHD-Effekte erweitert, die bei weichen, verformbaren Partikeln (wie Mikrogelen) auftreten. Dies geschieht durch Hinzufügen eines zusätzlichen Dehnungsbeitrags $\Gamma_{EHD}$ zur Maxwell-Gleichung, der proportional zum Quadrat der Spannung skaliert ($\dot{\Gamma}_{EHD} \propto \Sigma^2$). Dies modelliert das Gleiten der Partikel an den Wänden durch dünne Lösungsmittelfilme.

3. Wichtige Beiträge und Ergebnisse

A. Quantitative Erfassung des Stress-Overshoots

Das Modell sagt erfolgreich voraus, wie das maximale Spannungsüberschuss ($\Sigma_M$) von der Scherrate $\dot{\Gamma}$ abhängt. Die Analyse zeigt zwei Skalierungsregime:

• Diffusives Regime (kleine Scherraten): Der Anstieg des Scherbands wird durch Diffusion dominiert.
• Asymptotisches Regime (große Scherraten): Die Dynamik folgt einem anderen Potenzgesetz.
  Die resultierende Skalierungsgleichung lautet:
  $$\Sigma_M - 1 \sim B \left( \frac{\dot{\Gamma}}{\xi \tau^{1/2}} \right)^{\beta(n)} + C \left( \frac{\dot{\Gamma}}{\xi \tau} \right)^{\alpha(n)}$$
  mit Exponenten $\beta(n) = 2n/3$ und $\alpha(n) = 4n/(9-n)$, wobei $n$ der Scherverdünnungs-Exponent im Herschel-Bulkley-Gesetz ist. Diese Vorhersagen stimmen hervorragend mit experimentellen Daten an Carbopol-Mikrogelen überein.

B. Dynamik der Scherbänder und Fluidisierungszeit

Das Modell beschreibt die zeitliche Entwicklung des Scherbands ($\ell_b$), das von der Wand aus in das Material hineinwächst.

• Es wird eine Fluidisierungszeit $T_f$ definiert, bis der gesamte Probenraum fließt.
• Für konstante Scherrate ergibt sich die Skalierung $T_f \sim \dot{\Gamma}^{-9/4}$.
• Für konstante Spannung ergibt sich eine andere Skalierung. Das Verhältnis der Exponenten hängt direkt vom Herschel-Bulkley-Exponenten $n$ ab, was experimentell bestätigt wurde.

C. Einfluss von EHD-Wechselwirkungen und Randbedingungen

Durch die Einbeziehung der EHD-Interaktionen (Modellierung des Wandgleitens):

• Wird die stationäre Fließkurve bei niedrigen Scherraten korrekt wiedergegeben (Abweichung vom reinen Herschel-Bulkley-Verhalten).
• Ändert sich die Skalierung des Stress-Overshoots im EHD-dominierten Bereich zu $\Sigma_M \sim \dot{\Gamma}^{1/2}$.
• Zeigt sich, dass die Art der Randbedingungen (Festhalten der Fluidität vs. Festhalten des Gradienten) entscheidet, ob das Material duktil (durch Scherbänder) oder spröde (durch abrupten Spannungsabfall) versagt.

4. Bedeutung und Fazit

Der Artikel liefert einen pädagogischen und rigorosen Rahmen zur Beschreibung des Fließübergangs weicher glasartiger Materialien.

• Theoretische Leistung: Die Autoren zeigen, dass ein nicht-lokales Kontinuumsmodell, das auf der Fluidität basiert, in der Lage ist, komplexe, zeitabhängige Phänomene (Overshoot, Scherbänder) quantitativ vorherzusagen, ohne auf mikroskopische Simulationen zurückgreifen zu müssen.
• Universalität: Die gefundenen Skalierungsgesetze sind robust und gelten für verschiedene Materialien, solange die makroskopischen Parameter (wie $n$ und $\xi$) bekannt sind.
• Praktische Relevanz: Das Modell erklärt, warum Randbedingungen (Rauigkeit, Gleiten) einen entscheidenden Einfluss auf das makroskopische Verhalten haben. Dies ist für Anwendungen wie additive Fertigung, Beschichtungsverfahren und Haftbewegung von großer Bedeutung.
• Offene Fragen: Als Ausblick wird gefordert, das Modell auf zeitabhängige Protokolle (z. B. Rampen) zu erweitern, um Hysterese-Effekte zu untersuchen, und die mikroskopischen Ursachen der Randbedingungen weiter zu entschlüsseln.

Zusammenfassend stellt dieser Ansatz eine Brücke zwischen mikroskopischen Mechanismen (plastische Ereignisse) und makroskopischen rheologischen Beobachtungen dar und bietet ein leistungsfähiges Werkzeug zur Vorhersage des Verhaltens von SGMs unter komplexen mechanischen Belastungen.