🔬 condensed matter

A Smoluchowski equation for a sheared suspension of frictionally interacting rods

Dit artikel ontwikkelt een Smoluchowski-vergelijking en spannings-tensor voor dichte, afgeschoven suspensies van wrijvingsinteragerende staafjes door de Onsager-variatiemethode van Doi toe te passen om zowel vaste als gesmeerde wrijving te incorporeren, terwijl ook een gegeneraliseerd model voor het gemiddelde aantal interacties tussen staafjes wordt voorgesteld.

Oorspronkelijke auteurs: Genevieve Quiñones, Peter D. Olmsted

Gepubliceerd 2026-07-20

📖 8 min leestijd🧠 Diepgaand

CC BY 4.0

Oorspronkelijke auteurs: Genevieve Quiñones, Peter D. Olmsted

Oorspronkelijk artikel gelicentieerd onder CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ Dit is een AI-gegenereerde uitleg van het onderstaande artikel. Het is niet geschreven of goedgekeurd door de auteurs. Raadpleeg het oorspronkelijke artikel voor technische nauwkeurigheid. Lees de volledige disclaimer

Stel je een wereld voor waarin de dingen die we drinken, verven of printen niet slechts een eenvoudige vloeistof zijn, maar een chaotische dans van kleine, stijve stokjes die in een vloeistof zweven. Dit is het domein van suspensies, een tak van de natuurkunde die bestudeert hoe deeltjes bewegen wanneer ze gemengd zijn in een vloeistof. Al meer dan een eeuw begrijpen wetenschappers hoe deze stokjes zich gedragen als ze ver uit elkaar staan, drijvend als bladeren in een zachte stroom. Ze weten ook hoe ze zich gedragen als ze druk op elkaar staan, maar nog steeds langs elkaar glijden op een dun laagje gladde vloeistof, zoals schaatsers op een bevroren vijver. Maar er is een derde, rommeliger scenario: wat gebeurt er wanneer de menigte zo dicht wordt dat het "ijs" breekt? Wanneer de stokjes tegen elkaar botsen, hun gladde laagje verliezen en tegen elkaar beginnen te schuren met wrijving? Dit is de vraag die de moderne rheologie (de studie van stroming) aandrijft. Het begrijpen hiervan is cruciaal omdat het verklaart waarom sommige materialen plotseling veranderen van een dunne vloeistof in een solide blok wanneer je ze te snel roert, een fenomeen dat industriële leidingen kan verstoppen of 3D-printers kan laten vastlopen.

In dit artikel pakken Genevieve Quiñones en Peter D. Olmsted de rommelige wiskunde van deze wrijvingsvolle botsingen aan. Ze bouwen een nieuwe set regels, de Smoluchowski-vergelijking genoemd, om te voorspellen hoe een dichte menigte wrijvende staafjes zich gedraagt wanneer deze wordt afgeschuurd (geroerd). Denk eraan als het upgraden van een verkeerssimulatie van een model waarbij auto's alleen op gladde, wrijvingsloze wegen rijden naar een model waarbij ze daadwerkelijk kunnen slippen, schuren en vastlopen. De auteurs suggereren dat wanneer deze staafjes dicht op elkaar gepakt zitten, ze niet alleen glijden; ze blijven soms ook plakken en schrapen, wat een "wrijvingskoppel" (frictional torque) creëert dat hun rotatie verandert. Hun werk suggereert dat deze wrijving werkt als een rem, die de wilde, tollende bewegingen die de staafjes normaal gesproken maken onderdrukt en hen dwingt om netter met de stroming mee te lijnen. Ze stellen ook een nieuwe manier voor om te tellen hoe vaak een staafje tegen zijn buren botst, waarbij ze ontdekken dat hoe ordelijker de menigte is, hoe minder botsingen er plaatsvinden. Hoewel ze nog geen fysieke machine hebben gebouwd om dit te bewijzen, laten hun computersimulaties zien dat het toevoegen van wrijving de danspassen van de staafjes verandert, waardoor ze niet uit de hand draaien, wat potentieel verklaart waarom sommige dikke vloeistoffen plotseling nog dikker worden wanneer je er harder op drukt.

Het verhaal van de schurende stokjes

Stel je voor dat je op een drukke dansvloer staat. Als de ruimte ruim is, kan iedereen vrij draaien en tollen. Dit is als een dilute suspensie (dunne suspensie) van staafjes, waarbij de deeltjes ver uit elkaar liggen en slechts incidenteel tegen elkaar botsen. Natuurkundigen weten al heel lang hoe ze deze dans kunnen voorspellen; ze gebruiken vergelijkingen (zoals de vergelijking van Jeffery) die beschrijven hoe een enkele staaf roteert in een stromende vloeistof.

Stel je nu voor dat de dansvloer vol raakt. De dansers zijn niet langer alleen aan het draaien; ze botsen tegen elkaar aan. In een semi-dilute (half-geconcentreerde) menigte glijden ze misschien nog steeds langs elkaar heen dankzij een dun laagje vloeistof (zoals zweet of olie) dat als smeermiddel fungeert. Maar in een concentrated (geconcentreerde) menigte kan dat smeerlaagje breken. Plotseling zijn de dansers niet meer aan het glijden; ze zijn hun schoenen tegen elkaar aan het schuren. Dit is vaste wrijving (solid friction).

Al een lange tijd worstelen wetenschappers met het opschrijven van de wiskunde voor dit schuren. Eerdere pogingen gingen ervan uit dat zelfs wanneer de staafjes schuren, ze nog steeds dezelfde "danspassen" volgen (de vergelijking van Jeffery) alsof ze in een ijle menigte zijn, maar dan met een extra "duw" toegevoegd. Maar toen de auteurs van dit artikel probeerden het gemiddelde effect van deze duwen te berekenen, zagen ze dat de wiskunde wegviel naar nul. Het was alsof de wrijving onzichtbaar was. Dit maakte voor hen geen zin, want we weten uit de praktijk dat wrijving zaken vertraagt en verandert hoe ze bewegen.

De nieuwe aanpak: De "energierekening"

Quiñones en Olmsted besloten een andere strategie te proberen. In plaats van te raden hoe de staafjes bewegen en dan wrijving toe te voegen, begonnen ze bij de energie. Ze gebruikten een methode genaamd de Rayleighian, wat in essentie een manier is om de "energierekening" van het hele systeem te berekenen.

Denk aan de Rayleighian als een totale kostencalculator. Deze telt twee dingen op:

De Vrije Energie: Hoe graag de staafjes zichzelf op een bepaalde manier willen ordenen (zoals magneten die op één lijn liggen).
De Dissipatie: Hoeveel energie verloren gaat aan warmte door wrijving en weerstand.

De auteurs realiseerden zich dat om de ware beweging van de staafjes te vinden, je de totale "energierekening" moet minimaliseren. Toen zij de energie die verloren gaat aan vaste wrijving (het schuren) en grensvlak-gesmeerde wrijving (het glijden met een dunne film) in deze rekening opnamen, gebeurde er iets magisch. De wiskunde produceerde vanzelf een nieuwe "koppel" (een draaiende kracht) die niet nul was.

Ze ontdekten dat wrijving werkt als een rotatierem. Wanneer de staafjes tegen elkaar schuren, wordt het moeilijker voor hen om te draaien. Dit verandert de Smoluchowski-vergelijking, de hoofdberekening die voorspelt hoe de oriëntatie van alle staafjes in de loop van de tijd verandert.

Wat ze ontdekten: De dans verandert

De auteurs voerden computersimulaties uit om te zien wat er gebeurt als je deze wrijvende staafjes roert. Dit is wat zij ontdekten:

Het "kayaking" stopt: In een wrijvingsloze menigte voeren de staafjes vaak een beweging uit die "kayaking" wordt genoemd, waarbij ze kantelen en tollen in een complexe, uit het vlak lopende baan. De auteurs ontdekten dat wanneer wrijving wordt ingeschakeld, deze wilde draaiende beweging wordt onderdrukt. De staafjes stoppen met "kayaken" en hebben in plaats daarvan de neiging om in lijn te blijven met de stroming, zoals een school vissen die in een rechte lijn zwemt.
De link met "jamming": Ze suggereren dat deze wrijvingsrem een sleutelcomponent is van Discontinuous Shear Thickening (DST). Dit is het vreemde fenomeen waarbij een vloeistof (zoals maïszetmeel en water) plotseling verandelt in een vaste stof wanneer je er hard op slaat. Het model van de auteurs suggereert dat naarmate je harder roert, de wrijving tussen de staafjes toeneemt, de "remmen" blokkeren en het materiaal vastloopt (jamming).
Het tellen van de botsingen: Ze kwamen ook met een nieuwe formule om te schatten hoe vaak een staafje zijn buren aanraakt. Ze ontdekten dat het aantal contacten afhangt van hoe geordend de menigte is. Als de staafjes perfect op één lijn liggen (hoge orde), raken ze minder buren aan. Als ze rommelig en ongeordend zijn, raken ze meer buren aan. Dit is een generalisatie van een eerder idee dat alleen werkte voor rommelige, ongeordende menigten.

De beperkingen van het model

Het is belangrijk om te vermelden wat dit artikel niet beweert. De auteurs benadrukken dat hun model ervan uitgaat dat de staafjes altijd langs elkaar glijden. Ze kunnen nog niet het moment modelleren waarop de staafjes zo erg vast komen te zitten dat ze een enorme, onbeweeglijke cluster vormen (statische wrijving). In de echte wereld zijn deze vastgelopen clusters waarschijnlijk de belangrijkste reden dat materialen volledig "jammen". De auteurs geven toe dat hun model dit "vastzittende" deel mist, omdat statische wrijving geen warmte genereert (dissipatie), en dus niet voorkomt in hun energierekening.

Bovendien zijn hun resultaten gebaseerd op simulaties en wiskundige afleidingen, niet op nieuwe fysieke experimenten. Ze suggereren dat hun bevindingen verklaren waarom wrijving ertoe doet, maar ze hebben het nog niet in een laboratorium gemeten. Ze gaan er ook vanuit dat de suspensie perfect gemengd en uniform is, waarbij ze de mogelijkheid negeren dat wrijving ervoor kan zorgen dat de staafjes op specifieke plekken samenklonteren.

De kern van het verhaal

In simpele termen hebben Quiñones en Olmsted een nieuwe wiskundige gereedschapskist gebouwd die eindelijk rekening houdt met het "schuren" tussen deeltjes in een dikke vloeistof. Ze lieten zien dat als je wrijving negeert, je de remmen mist die de staafjes ervan weerhouden wild rond te tollen. Door deze wrijving op te nemen, vonden ze een zelfconsistente manier om te voorspellen hoe deze drukke staafjes bewegen en hoeveel spanning ze leggen op het vat dat hen vasthoudt. Hoewel ze het mysterie nog niet hebben opgelost van waarom vloeistoffen volledig vastlopen (dat vereist het begrijpen van de "vastzittende" staat), hebben ze een cruciaal puzzelstuk geleverd: het wrijvingskoppel dat de dans vertraagt voordat de blokkade optreedt.

Technische Samenvatting: Een Smoluchowski-vergelijking voor een Geschoren Suspensie van Frictie-interagerende Staafjes

Probleemstelling
De rheologie van dichte suspensies van rigide staafvormige deeltjes onder afschering (shear) wordt niet volledig gevangen door bestaande theorieën, met name in regimes die discontinue afscheringsverdikking (DST), afscheringsjamming en vorticiteit-kanteling vertonen. Hoewel klassieke modellen (bijv. Jeffery, Doi, Advani-Tucker) succesvol zijn in het beschrijven van oriëntatiedynamica in ijle en semi-ijle regimes via hydrodynamische interacties en uitgesloten volume-effecten, falen zij in het verklaren van solide contacten tussen deeltjes. Recente experimentele en computationele bewijzen suggeren dat wrijvingscontacten tussen staafjes de primaire drijfveer zijn achter deze onverklaarde rheologische gedragingen. Echter, eerdere pogingen om wrijving in de dynamica van staafjes te integreren (bijv. Sandstrom & Tucker, Bounoua et al.) leverden nulresultaten op voor het gemiddelde wrijvingskoppel, grotendeels omdat zij aannamen dat de beweging van de staafjes de vergelijking van Jeffery volgde en simpelweg een koppelterm toevoegden, wat leidde tot annuleringen bij het middelen.

Methodologie
De auteurs ontwikkelen een constitutief kader voor dichte, geschorende suspensies van frictie-interagerende staafjes door de variatiemethode van Onsager toe te passen, zoals geformuleerd door Doi. De kern van de methodologie bestaat uit het construeren van een Rayleighiaan ( $R$ ), gedefinieerd als de som van de verandering in vrije energie ( $\dot{A}$ ) en een dissipatiefunctie ( $\Phi$ ) die afhankelijk is van de hoeksnelheid ( $\omega$ ) van de staafjes.

Dissipatiefuncties: De auteurs leiden dissipatiefuncties af voor twee verschillende wrijvingsmechanismen:
- Grensvlak-gesmeerde (BL) Wrijving: Schaalt lineair met de relatieve tangentiële contactsnelheid ( $v_c$ ), gemodelleerd als $F \propto v_c$ .
- Solide Amontons-Coulomb Wrijving: Schaalt met de grootte van de normale kracht ( $|F_N|$ ) en is enkel afhankelijk van de richting van $v_c$ , gemodelleerd als $F \propto |F_N|$ .
Mean Field Benadering: Om de veel-deeltjes-natuur van contacten aan te pakken, gebruiken de auteurs een mean field benadering. Zij leiden een uitdrukking af voor het gemiddelde aantal contacten per staafje, $\langle c \rangle$ , als functie van de volumefractie ( $\phi$ ), aspectratio ( $L/D$ ) en de nematische ordeparameter ( $S$ ). Dit generaliseert de willekeurige contactvergelijking van Philipse voor wanordelijke pakkingen om oriëntatie-orde te bevatten.
Variatieve Afleiding: In tegen tegenstelling tot eerdere benaderingen die een vooraf bepaalde bewegingsvergelijking aannamen, minimaliseren de auteurs de totale Rayleighiaan met betrekking tot $\omega$ . Dit levert een zelfconsistentie-vergelijking van beweging op waarbij het wrijvingskoppel van nature een functie is van de hoeksnelheid zelf.
Smoluchowski-vergelijking en Spanning: De minimalisering van de Rayleighiaan produceert een nieuwe Smoluchowski-vergelijking voor de oriëntatieverdelingsfunctie $\psi(\hat{u})$ en een bijbehorende spanningstensor $\sigma$ . Het model incorporeert een schakelfunctie $\Theta(p)$ om over te gaan tussen de gesmeerde en de solide wrijvingsregimes op basis van de gereduceerde druk.

Belangrijkste Bijdragen

Niet-nul Wrijvingskoppel: Het artikel demonstreert dat door de dynamica af te leiden uit een dissipatiefunctie in plaats van een koppel toe te voegen aan de vergelijking van Jeffery, een niet-nul gemiddeld wrijvingskoppel wordt verkregen. Dit lost het "nulresultaat" op dat werd gevonden in eerdere literatuur (bijv. Bounoua et al.), welke voortkwam uit het middelen van oneven machten van contactpositievariabelen.
Gegeneraliseerd Contactgetal: De auteurs stellen een vorm voor van het gemiddelde aantal contacten $\langle c \rangle$ die afhankelijk is van de nematische ordeparameter $S$ , specifiek $\langle c \rangle \approx \phi (L/D) (1 - S^2)$ . Dit koppelt de microstructurele aard van contacten direct aan de macroscopische oriëntatie-orde.
Zelfconsistente Dynamica: De resulterende Smoluchowski-vergelijking is zelfconsistent; de effectieve rotatieweerstandscoëfficiënt hangt af van de contactsnelheid, die op zijn beurt afhankelijk is van de hoeksnelheid. Deze niet-lineariteit is inherent aan het Coulomb-wrijvingsmodel.
Spanningstensor Formule: Een nieuwe spanningstensor is afgeleid die bijdragen bevat van zowel hydrodynamische weerstand als wrijvingsdissipatie, waarbij het wrijvingscomponent schaalt met het gemiddelde wrijvingskoppel.

Resultaten en Numerieke Studie
De auteurs voeren een numerieke studie uit van de afgeleide Smoluchowski-vergelijking (benaderd via een getrunceerde reeksexpansie en Langevin-dynamica) voor een eenvoudige afschering (simple shear flow).

Gesmeerd Limiet: In de afwezigheid van vaste wrijving reduceert het model tot de Doi-Edwards theorie met een gerenormaliseerde rotatieweerstandconstante, wat de rotatie van de staafjes effectief vertraagt.
Vaste Wrijvingslimiet: Wanneer vaste wrijving actief is ( $\mu_k = 1$ $μ_{k} = 1$ ), verandert de dynamica aanzienlijk. Het model voorspelt een onderdrukking van de director-beweging in het stroom-gradiëntvlak.
- Flow Alignment: In tegenstelling tot het wrijvingsloze geval waar de director "kayaking" (uit-het-vlak tuimelen) of periodieke banen kan vertonen bij intermediaire Peclet-getallen, stabiliseert de inclusie van wrijving de director in een flow-aligned staat.
- Initiële Condities: Het wrijvingsmodel laat zien dat zelfs uitgaande van een isotrope verdeling, het systeem de neiging heeft om met de stroming mee te lijnen en daar te blijven, met slechts sporadische excursies naar kayaking-banen.
- Vorticiteit-kanteling: Bij lage Peclet-getallen drijft de vorticiteit-component van de director nabij nul, in contrast met het grillige dwalen gezien in het wrijvingsloze geval.

Significantie
Het artikel beweert een theoretische beschrijving te bieden van dichte staaf-suspensies die expliciet rekening houdt met wrijvingscontacten, een factor die eerder werd genegeerd of onjuist werd behandeld in constitutieve modellen. Door het variatieprincipe van Onsager te gebruiken, vestigen de auteurs een zelfconsistent kader waarin wrijvingsdissipatie van nature de noodzakelijke koppel genereert om de oriëntatiedynamica te veranderen. Het werk suggereert dat wrijvingscontacten een cruciaal mechanisme zijn voor de onderdrukking van complexe oriëntatiebewegingen (zoals kayaking) en mogelijk verantwoordelijk zijn voor het ontstaan van afscheringsverdikking en jamming in dichte staaf-suspensies. De auteurs merken op dat hoewel hun model een ruimtelijke homogeniteit en glijdende wrijving veronderstelt (exclusief statische kleverige clusters), het een rigoureuze fundering biedt voor het begrijpen van de rheologische anomalieën waargenomen in geconcentreerde staaf-suspensies.

Het verhaal van de schurende stokjes

De nieuwe aanpak: De "energierekening"

Wat ze ontdekten: De dans verandert

De beperkingen van het model

De kern van het verhaal

Meer zoals dit