🔬 condensed matter

Criticality of the viscous to inertial transition near jamming in non-Brownian suspensions

Mithilfe von Diskret-Element-Methode-Simulationen zeigt diese Studie, dass der Übergang von viskosem zu inertiellem Scherverdicken in dichten reibungsfreien nicht-Brownschen Suspensionen nahe dem Verklemmungspunkt durch eine divergierende mikrostrukturelle Längenskala gesteuert wird, die dazu führt, dass die charakteristische Übergangsscherrate verschwindet, wenn sich der Volumenanteil dem Verklemmungslimit nähert.

Ursprüngliche Autoren: Nishanth Murugan, Donald Koch, Sarah Hormozi

Veröffentlicht 2026-07-16

📖 5 Min. Lesezeit🧠 Tiefgang

CC BY 4.0

Ursprüngliche Autoren: Nishanth Murugan, Donald Koch, Sarah Hormozi

Originalarbeit lizenziert unter CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ Dies ist eine KI-generierte Erklärung des untenstehenden Papers. Sie wurde nicht von den Autoren verfasst oder gebilligt. Für technische Genauigkeit konsultieren Sie das Originalpaper. Vollständigen Haftungsausschluss lesen

Stellen Sie sich eine Welt vor, in der sich die Regeln des Fließens ändern, je nachdem, wie fest man sie drückt. Dies ist das faszinierende Reich der Fluiddynamik, das sich speziell auf „Suspensionen“ konzentriert – Mischungen, bei denen winzige feste Partikel in einer Flüssigkeit schweben, wie Sand in Wasser oder Schokoladenstückchen in Milch. Normalerweise denken wir bei Fluiden entweder an dünnflüssig (wie Wasser) oder dickflüssig (wie Honig), aber diese partikelfüllten Mischungen können tückisch sein. Manchmal, wenn man sie langsam rührt, verhalten sie sich wie eine normale Flüssigkeit. Aber wenn man sie schnell rührt, können sie plötzlich steif werden, fast wie ein Festkörper. Dieses Phänomen wird als „Scherverdickung“ bezeichnet. Wissenschaftler sind besessen davon, genau zu verstehen, wann und warum dieser Wechsel stattfindet, da dies alles beeinflusst – von der Verarbeitung von Zement und Keramik bis hin zur Vorhersage gefährlicher Naturereignisse wie Schlammströme und Lawinen. Das große Rätsel war bisher: Was ist das exakte „Tempolimit“, bei dem die Flüssigkeit entscheidet, aufzuhören, flüssig zu sein und anzfangen, steif zu werden? Ist es eine feste Zahl oder ändert sie sich je nach der Menge der Partikel?

In dieser Studie nutzten die Forscher Nishanth Murugan, Donald Koch und Sarah Hormozi von der Cornell University leistungsstarke Computersimulationen, um ein Teil dieses Puzzles zu lösen. Sie konzentrierten sich auf eine spezifische Art von Mischung: dichte Suspensionen aus Partikeln, die nicht aneinander haften (reibungslos) und zu groß sind, um durch Wärme herumgewirbelt zu werden (nicht-Brownsches Verhalten). Denken Sie an eine überfüllte Tanzfläche, auf der die Tänzer glatte Anzüge tragen und nicht gegeneinander stoßen können. Das Team wollte den präzisen Moment finden, in dem die Tanzfläche von einem geschmeidigen Gleiten zu einem chaotischen, drängelnden Gedränge wechselt.

Ihre Simulationen zeigten, dass es kein einzelnes, universelles Tempolimit für diesen Übergang gibt. Stattdessen hängt das „Tempolimit“ vollständig davon ab, wie dicht gepackt die Tanzfläche ist. Je näher die Partikel dem Punkt kommen, an dem sie vollständig verklemmt sind und sich gar nicht mehr bewegen können (ein Zustand, der als „Jamming“ oder Verkeilung bezeichnet wird), desto niedriger sinkt die Geschwindigkeit, die erforderlich ist, um das versteifende Verhalten auszulösen – bis hin zu fast Null. Mit anderen Worten: Je gedrängter das System wird, desto leichter lässt es sich in einen Festkörper verwandeln.

Die Autoren entdeckten, dass dies geschieht, weil eine verborgene, wachsende Struktur innerhalb der Mischung entsteht. Wenn die Partikel enger zusammenrücken, beginnen sie, riesige, temporäre Cluster zu bilden, die wie eine einzige, massive Einheit agieren. Die Größe dieser Cluster wächst immer weiter an, während sich das System dem Jamming-Zustand nähert. Die Forscher fanden heraus, dass der Übergang von „flüssig“ zu „steif“ tatsächlich durch die Trägheit (den Impuls) dieser riesigen Cluster ausgelöst wird, nicht nur durch die einzelnen Partikel. Sie bewiesen dies, indem sie zeigten, dass die Größe dieser Cluster und die Geschwindigkeit, mit der die Flüssigkeit steif wird, demselben mathematischen Muster folgen.

Interessanterweise argumentiert die Arbeit gegen die Idee, dass dieser Übergang durch die Größe eines einzelnen Partikels oder eine feste Kraft gesteuert wird. In Systemen, in denen Partikel Reibung haben (wie grober Sand), findet der Übergang bei einer anderen, konstanten Geschwindigkeit statt. Aber für diese glatten, rutschigen Partikel ist der Übergang „kritisch“, was bedeutet, dass er tief mit der Dichte des Systems verwoben ist. Die Studie legt nahe, dass man sich dem Jamming-Punkt nähert, wächst eine „Längenskala“ (ein Maß dafür, wie groß diese Cluster werden) ins Unendliche, was das System extrem empfindlich für selbst den kleinsten Stoß macht.

Die Forscher verwendeten eine Methode namens „Discrete Element Method“ (DEM), was wie ein supergenaues Videospiel ist, bei dem sie die Kräfte zwischen Millionen von Partikeln berechnen. Sie simulierten sowohl 2D- als auch 3D-Versionen dieser überfüllten Systeme. Ihre Ergebnisse zeigten, dass die „Stokes-Zahl“ (eine schicke Art, das Gleichgewicht zwischen der Klebrigkeit der Flüssigkeit und dem Impuls der Partikel zu messen) gegen Null geht, wenn der Volumenanteil sich dem Jamming-Limit nähert. Sie fanden heraus, dass die Größe der Cluster gemäß einem spezifischen Potenzgesetz wächst, mit einem Exponenten von etwa 0,79 in 2D und 0,84 in 3D.

Indem sie testeten, wie sich die Ergebnisse veränderten, wenn sie die Simulationsbox kleiner oder größer machten, bestätigten sie, dass diese riesigen Cluster wirklich existieren und der Schlüssel zum Übergang sind. Wenn die Simulationsbox kleiner als der Cluster war, änderten sich die Ergebnisse, was bewies, dass der Cluster versuchte, „durch die Wände zu greifen“. Dies bestätigte, dass der Übergang ein nicht-lokales Ereignis ist, was bedeutet, dass das, was in einem Teil der Flüssigkeit passiert, von dem abhängt, was weit entfernt geschieht, verbunden durch diese massiven, unsichtbaren Strukturen.

Kurz gesagt legt die Arbeit nahe, dass für glatte, überfüllte Partikel der Wechsel von flüssig zu fest nicht über die einzelnen Tänzer entscheidet, sondern über die Bildung einer massiven, unsichtbaren Menge, die wächst, bis sie den ganzen Raum übernimmt. Diese Erkenntnis hilft Wissenschaftlern zu verstehen, dass in der Nähe des Jamming-Punktes die Regeln des Fließens kritisch werden und von diesen wachsenden Strukturen abhängen, was den Weg für bessere Theorien darüber ebnet, wie sich diese komplexen Materialien verhalten.

Technische Zusammenfassung: Kritikalität des Übergangs von viskosem zu inertiellem Verhalten nahe dem Verkeilen in nicht-Brownschen Suspensionen

Problemstellung
Die Rheologie dichter nicht-Brownscher Suspensionen ist entscheidend für das Verständnis natürlicher Phänomene (z. B. Erdrutsche, Schlammlawinen) und industrieller Prozesse (z. B. Zementleim, Keramik). Ein zentrales Merkmal dieser Systeme ist der Übergang von einem viskosen Regime, das durch eine Newtonsche Rheologie ( $\Sigma \propto \dot{\gamma}$ ) gekennzeichnet ist, zu einem inertiellem Regime, das durch eine Bagnoldsche Rheologie ( $\Sigma \propto \dot{\gamma}^2$ ) gekennzeichnet ist. Dieser Übergang wird traditionell der zunehmenden Dominanz der Partikelträgheit gegenüber den hydrodynamischen Kräften zugeschrieben.

Der physikalische Mechanismus, der diesen Übergang auslöst, und die charakteristische Scherrate, bei der er auftritt, bleiben jedoch Gegenstand von Debatten. Die experimentelle und numerische Literatur präsentiert widersprüchliche Ergebnisse bezüglich der kritischen Stokes-Zahl ( $St_{v \to i}$ ), eines dimensionslosen Maßes für das Verhältnis zwischen inertiellem und viskosem Stress. In reibungsfreien Suspensionen (z. B. Polystyrolkugeln mit elektrostatischer Abstoßung) wurde beobachtet, dass $St_{v \to i}$ mit dem Volumenanteil ( $\phi$ ) variiert und gegen Null geht, wenn das System sich dem Verkeilungsvolumenanteil ( $\phi_m$ ) nähert. Im Gegensatz dazu scheint $St_{v \to i}$ in reibungsbehafteten Suspensionen (z. B. raue PMMA-Partikel) unabhängig von $\phi$ zu sein und bleibt endlich. Das Fehlen eines Konsenses über den physikalischen Mechanismus, insbesondere die Rolle der Mikrostruktur in reibungsfreien Systemen, motiviert diese Studie.

Methodik
Die Autoren verwenden Discrete-Element-Method (DEM)-Simulationen, um dichte, nicht-Brownsche, reibungsfreie Suspensionen über ein breites Spektrum an Scherraten und Volumenanteilen zu untersuchen.

Kraftmodell: Die Dynamik wird durch die Newtonschen Bewegungsgleichungen bestimmt, die zwei primäre Kräfte berücksichtigen: hydrodynamische Schmierkräfte ( $F^L_{ij}$ ), die aus der relativen Bewegung im Zwischenraum der Partikel resultieren, und repulsive Kontaktkräfte ( $F^R_{ij}$ ).
Reibungsfreie Näherung: Um die in Experimenten vorkommende elektrostatische Abstoßungsbarriere (die direkten Kontakt verhindert) zu modellieren, wird die repulsive Kraft als Hookesche Feder ( $F^R_{ij} = K\delta n$ ) modelliert, die aktiviert wird, wenn der Oberflächenabstand $h$ unter einen minimalen Schwellenwert $h_{min}$ fällt.
Simulationsaufbau: Die Simulationen werden sowohl in 2D- als auch in 3D-Domänen unter Verwendung von Lees-Edwards-Randbedingungen durchgeführt, um eine Scherströmung zu simulieren. Das System verwendet eine polydisperse Partikelverteilung (16 %), um Kristallisation zu verhindern. Die Schmierkräfte werden unter der Annahme eines minimalen Abstands von $h_{min}$ berechnet, wenn Partikel versuchen, sich näher zu kommen, was effektiv die Normalkonstruktion einschränkt.
Parameter: Die Studie untersucht den Parameterraum, der durch den Volumenanteil ( $\phi$ ), die Stokes-Zahl ( $St = 4\rho a^2 \dot{\gamma} / \eta_f$ ) und das dimensionslose Steifigkeitsverhältnis ( $\Gamma$ ) definiert ist. Die Simulationen werden in einem Regime durchgeführt, in dem $\Gamma \gg 1$ gilt, wodurch sichergestellt ist, dass die repulsive Barriere gegenüber den Schmierkräften bei kleinen Abständen dominiert.

Wesentliche Beiträge und Ergebnisse

Kritikalität des Übergangs: Die Simulationen bestätigen, dass für reibungsfreie Suspensionen die charakteristische Stokes-Zahl für den Übergang zu inertiellem Scherverdickung ( $St_{v \to i}$ ) nicht konstant ist. Stattdessen sinkt sie, wenn der Volumenanteil steigt, und nähert sich Null an, wenn $\phi \to \phi_m$ . Die Beziehung folgt einem Potenzgesetz: $St_{v \to i} \sim (\phi_m - \phi)^\beta$ .
Skalierungsrahmen: Um dieses Verhalten zu erklären, schlagen die Autoren einen Skalierungsrahmen vor, der auf der Hypothese basiert, dass eine mikrostrukturelle Längenskala beim Annähern des Systems an die Verkeilung divergiert. Sie führen eine modifizierte Stokes-Zahl ein, $St_{Cluster} = B St (\phi_m - \phi)^{-\beta}$ , welche die inertiellem Effekte relativ zu einer wachsenden Clustergröße $L_c \sim (\phi_m - \phi)^{-\beta/2}$ skaliert.
Datenkollaps: Unter Verwendung dieses Skalierungsrahmens kollabieren die rheologischen Daten (relative Scherviskosität $\eta_r$ $η_{r}$ vs. $St$) für verschiedene Volumenanteile auf einer einzigen Masterkurve in sowohl 2D- als auch 3D-Simulationen.
- In 2D ist der Exponent $\beta \approx 0,79$ .
- In 3D ist der Exponent $\beta \approx 0,84$ .
- Diese Werte stimmen mit experimentellen Beobachtungen überein, die darauf hindeuten, dass $St_{v \to i} \sim (\phi_m - \phi)^{-1}$ .
Systemgrößenabhängigkeit: Um die Existenz der vorgeschlagenen divergierenden Längenskala zu validieren, führten die Autoren Tests zur Systemgrößenabhängigkeit durch. Sie beobachteten, dass die Viskosität ansteigt, wenn die Größe der Simulationsbox kleiner als die charakteristische Clustergröße ist, da der Cluster durch die periodischen Randbedingungen mit sich selbst interagiert. Durch Anpassung der Viskositätsdaten an ein Potenzgesetz in Abhängigkeit von der Systemgröße konnten sie erfolgreich eine Cluster-Längenskala $L_c$ extrahieren, die mit einem Exponenten konsistent mit der rheologischen Skalierung ( $L_c \sim (\phi_m - \phi)^{-\beta/2}$ ) gegen $\phi \to \phi_m$ divergiert.

Bedeutung und Behauptungen
Das Paper behauptet, dass der Übergang von viskoser zu inertiellem Scherverdickung in reibungsfreien Suspensionen ein kritisches Phänomen ist, das durch eine divergierende mikrostrukturelle Längenskala (Cluster) getrieben wird, während sich das System der Verkeilung nähert.

Mechanismus: Die Kritikalität des Übergangs wird der Wachstumsrate dieser Cluster zugeschrieben, welche die Skala verändert, bei der die inertiellem Effekte dominant werden. Da die Clustergröße divergiert, verschwindet die charakteristische Scherrate, die erforderlich ist, um die inertiellem Effekte auszulösen.
Nicht-lokale Rheologie: Die Existenz einer divergierenden Längenskala impliziert, dass die Rheologie dieser Suspensionen nicht-lokal ist; die Spannungsübertragung hängt von der Bildung und Zerstörung dieser großräumigen Cluster ab.
Kontrast zu reibungsbehafteten Systemen: Die Autoren merken an, dass dieses kritische Verhalten spezifisch für reibungsfreie Systeme ist. In reibungsbehafteten Suspensionen, in denen die Partikelkontakte durch Oberflächenrauheit statt durch langreichweitige Abstoßung bestimmt werden, scheint der Übergang unabhängig vom Volumenanteil zu sein, was darauf hindeutet, dass der mikrostrukturelle Mechanismus sich grundlegend unterscheidet.
Zukünftige Ausrichtung: Die Arbeit motiviert die Entwicklung nicht-lokaler Theorien, die die Spannungsübertragung durch diese wachsenden Cluster berücksichtigen. Sie hebt zudem die Notwendigkeit hervor, weitere Studien zu den Skalierungsgesetzen für Clustergrößen in reibungsbehafteten Suspensionen durchzuführen, um eine vereinheitlichte Theorie über verschiedene Strömungsregime hinweg zu etablieren.

Mehr davon