Discontinuous change of viscosity in a sheared granular gas with velocity-dependent restitution

Dit artikel toont aan dat een granulaire gas met een snelheidsafhankelijke terugkaatsingscoëfficiënt een discontinuïteit in de viscositeit vertoont bij een kritische schuifsnelheid, wat voortkomt uit zuiver kinetische effecten in plaats van wrijving of verstopping.

De kern

De Viscositeits-Magie: Waarom Korrelachtig Gas Soms "Plotseling" Harder Wordt

Stel je voor dat je een bak vol met kleine, harde balletjes hebt. Normaal gesproken gedragen deze balletjes zich als een vloeibare soep: hoe harder je roert (de "schuifsnelheid"), hoe makkelijker ze bewegen. Maar wat als deze balletjes een geheime eigenschap hebben die ze laat zien afhankelijk van hoe hard ze tegen elkaar botsen?

Dat is precies wat deze studie van de Japanse onderzoekers ontdekt. Ze keken naar een "korrelgas" (een verzameling van kleine deeltjes die niet aan elkaar plakken) en ontdekten een vreemd fenomeen: de vloeistof kan plotseling van "zeer soepel" naar "zeer stroperig" springen, zonder dat er iets anders verandert dan de snelheid waarmee je roert.

Hier is de uitleg in gewone taal, met een paar creatieve vergelijkingen:

1. De Regels van het Spel: De "Botsings-Rem"

In de echte wereld verliezen balletjes bij een botsing altijd een beetje energie (ze worden trager). Dit noemen we een restitutiecoëfficiënt.

• Normaal: Hoe harder ze botsen, hoe meer energie ze verliezen (ze remmen harder).
• In dit experiment: De onderzoekers verzonnen een speciaal type balletje dat een omgekeerde logica volgt.
  • Bij langzame botsingen gedragen ze zich als slapende leeuwen: ze verliezen heel veel energie (ze remmen hard).
  • Bij snelle botsingen gedragen ze zich als wilde panthers: ze verliezen heel weinig energie (ze stuiteren bijna perfect terug).

2. Het Experiment: De Dansvloer

Stel je een dansvloer voor waar duizenden balletjes rondrennen. Je begint langzaam te draaien (de dansvloer roteren).

• Langzaam roteren: De balletjes botsen zachtjes. Omdat ze bij zachte botsingen veel energie verliepen, gedragen ze zich als een dichte, stroperige massa. Ze bewegen moeizaam.
• Snel roteren: De balletjes botsen hard. Nu gedragen ze zich als snelle, energieke panthers die weinig energie verliezen. Ze bewegen soepel en snel.

3. Het Verrassende Moment: De "S-vormige" Sprong

Wat de onderzoekers vonden, is dat er een gevaarlijk middengebied is.
Stel je voor dat je de snelheid van de dansvloer langzaam opvoert.

1. Aan het begin is het stroperig (langzaam).
2. Je draait sneller, en plotseling springt de vloeistof naar een heel andere staat.
3. In dit middengebied kan de vloeistof drie verschillende soorten stroperigheid hebben bij exact dezelfde snelheid.

Dit is als een auto die op een bepaalde snelheid ineens niet meer kan accelereren, maar plotseling van versnelling springt naar een heel ander toerental. De grafiek van de stroperigheid ziet eruit als een S (of een trapje). Als je de snelheid iets verhoogt, gebeurt er niets... tot het punt waarop de vloeistof plotseling van "stroperig" naar "soepel" (of andersom) springt.

4. Waarom is dit speciaal? (De Vergelijking)

In de wetenschap kennen we dit gedrag al van dikke suspensies (zoals modder of verf met veel deeltjes). Daar gebeurt het omdat de deeltjes tegen elkaar blijven hangen en vastlopen (verstoppen).

• Het oude verhaal: "De deeltjes blijven aan elkaar plakken door wrijving, daarom wordt het hard."
• Dit nieuwe verhaal: "Er is geen wrijving en geen plakken. De balletjes zijn perfect glad. Het enige dat telt, is hoe hard ze tegen elkaar bonken."

Het is alsof je een dansvloer hebt waar de dansers niet vastlopen, maar waar hun stijl van dansen plotseling verandert afhankelijk van hoe snel ze bewegen. De "sprong" komt puur voort uit de manier waarop ze energie verliezen bij een botsing, niet omdat ze vastlopen.

5. De Conclusie: Een Nieuw Mechanisme

De onderzoekers tonen aan dat je niet nodig hebt dat deeltjes aan elkaar plakken of vastlopen om deze plotselinge veranderingen te krijgen. Het is genoeg dat de manier waarop ze energie verliezen, verandert afhankelijk van hun snelheid.

Samengevat in één zin:
Deze studie toont aan dat als je deeltjes hebt die bij langzame botsingen "slapen" en bij snelle botsingen "wakker worden", je een vloeistof kunt maken die plotseling van staat verandert, puur door de snelheid te veranderen, zonder dat er ook maar één deeltje vastloopt.

Het is een beetje alsof je een magische soep hebt die bij langzaam roeren als honing is, maar bij snel roeren ineens als water stroomt, en in het midden een raar punt heeft waar je niet weet of je een lepel of een hamer nodig hebt.

Technische samenvatting

Probleemstelling

Het artikel onderzoekt de reologie (de stroming en vervorming) van een verdund granulaire gas, bestaande uit harde bollen met een snelheidsafhankelijke restitutiecoëfficiënt. Granulaire systemen vertonen rijk niet-evenwichtsgedrag door inelastische botsingen. Hoewel modellen met een snelheidsafhankelijke restitutiecoëfficiënt ($e$) al bestudeerd zijn voor homogene koeltoestanden (homogeneous cooling state), is het gedrag onder stabiele afschuiving (steady shear) weinig onderzocht.

De kernvraag is of de competitie tussen verschillende dissipatiemechanismen (afhankelijk van de botsingssnelheid) kan leiden tot niet-triviale constitutieve relaties, zoals een discontinu verloop van de viscositeit. Specifiek wordt onderzocht wat er gebeurt als de restitutiecoëfficiënt schakelt tussen twee waarden ($e_1$ en $e_2$) bij een drempelsnelheid ($v_c$). Interessant is dat het artikel een scenario onderzoekt waarbij $e_1 < e_2$ (dissipatie is sterker bij lage snelheden dan bij hoge snelheden), wat tegenstrijdig lijkt met gedrag van droge granulaire materialen, maar wel relevant kan zijn voor geladen of adhezieve deeltjes.

Methodologie

De auteurs hanteren een combinatie van kinetische theorie en event-driven moleculaire dynamica (EDMD) simulaties:

1. Fysisch Model:

   • Een verdund, wrijvingsloos granulair gas van monodisperse harde bollen (massa $m$, diameter $d$, volumefractie $\phi = 0.01$).
   • De normale restitutiecoëfficiënt $e(v_n)$ hangt af van de relatieve snelheid loodrecht op de contactlijn ($v_n$):
     $$e(v_n) = e_1 - (e_1 - e_2)\Theta(v_n - v_c)$$
     Hierbij is $v_c$ de drempelsnelheid en $\Theta$ de Heaviside-stapfunctie.

2. Kinetische Theorie:

   • De dynamica wordt beschreven door de Boltzmann-vergelijking voor een uniform afschuivende stroming met afschuifsnelheid $\dot{\gamma}$.
   • De spanningstensor wordt afgeleid uit de tweede momenten van de snelheidsverdelingsfunctie.
   • Omdat een exacte oplossing van de verdelingsfunctie niet beschikbaar is, wordt Grad's benadering (een momentenmethode) gebruikt om de verdelingsfunctie te benaderen als een Maxwelliaanse verdeling gecorrigeerd met anisotropie-termen.
   • Hieruit worden de vergelijkingen voor de evolutie van de granulaire temperatuur ($T$), temperatuur-anisotropie ($\Delta T$) en de afschuifspanning ($P_{xy}$) afgeleid.

3. Stabiele Toestand en Viscositeit:

   • In de stabiele toestand worden de tijdsafgeleiden nul gesteld, waardoor een stelsel van gekoppelde algebraïsche vergelijkingen ontstaat.
   • De afschuifviscositeit $\eta$ wordt gedefinieerd als $\eta = -P_{xy} / \dot{\gamma}$.
   • De theorie voorspelt dat alle macroscopische grootheden in de stabiele toestand uitsluitend kunnen worden uitgedrukt als functies van de granulaire temperatuur $T$.

4. Validatie:

   • De theoretische voorspellingen worden gevalideerd met EDMD-simulaties voor harde bollen met de specifieke snelheidsafhankelijke restitutie.

Belangrijkste Bijdragen en Resultaten

1. S-vormige Constitutieve Kromme:

   • Wanneer de restitutiecoëfficiënt bij lage snelheden lager is dan bij hoge snelheden ($e_1 < e_2$), vertoont de viscositeit een S-vormige afhankelijkheid van de afschuifsnelheid ($\dot{\gamma}$).
   • Dit betekent dat er binnen een bepaald bereik van afschuifsnelheden drie mogelijke waarden voor de viscositeit bestaan voor één enkele opgelegde snelheid.

2. Discontinue Overgang:

   • Als de afschuifsnelheid geleidelijk wordt verhoogd, volgt het systeem de onderste tak van de S-kromme totdat de helling divergeert. Op dit kritieke punt ondergaat het systeem een discontinue sprong naar de bovenste tak (een plotselinge toename van de viscositeit).
   • Dit fenomeen lijkt op discontinue afschuivingverdikking (Discontinuous Shear Thickening - DST) die vaak wordt waargenomen in dichte suspensies.

3. Mechanisme: Kinetisch vs. Frictie:

   • Het cruciale onderscheid met bestaande modellen (zoals het Wyart-Cates-model voor suspensies) is dat deze overgang puur kinetisch is.
   • In tegenstelling tot suspensies, waar DST wordt veroorzaakt door de overgang van wrijvingsloze naar wrijvingscontacten en "jamming" (verstopping), bestaat dit systeem uitsluitend uit een verdund gas zonder blijvende contacten.
   • De overgang wordt gedreven door de competitie tussen twee Bagnold-type regimes:
     • Bij lage $\dot{\gamma}$ domineren botsingen met $v_n < v_c$ (effectieve restitutie $e_1$).
     • Bij hoge $\dot{\gamma}$ domineren botsingen met $v_n > v_c$ (effectieve restitutie $e_2$).
   • De S-vorm ontstaat door de niet-lineaire koppeling tussen de temperatuur en de afschuifsnelheid in het overgangsgebied waar beide regimes mengen.

4. Faseovergang en Voorwaarden:

   • De S-vormige structuur treedt alleen op als $e_1 < e_2$ en het verschil tussen $e_1$ en $e_2$ groot genoeg is.
   • Als $e_1 > e_2$ (het meer gebruikelijke geval voor droge korrels), is er geen shear thinning en blijft de viscositeit een eenduidige functie van de afschuifsnelheid.
   • De auteurs presenteren een fase-diagram in het $(e_1, e_2)$-vlak dat de grens aangeeft waar de discontinuïteit optreedt.

5. Validatie:

   • De theoretische krommen komen uitstekend overeen met de EDMD-simulaties, zowel in de asymptotische limieten (waar Bagnold-schaling geldt) als in het overgangsgebied.

Significantie en Conclusie

De studie toont aan dat discontinue reologische responsen niet noodzakelijk afhankelijk zijn van complexe interacties zoals wrijving, jamming of hoge dichtheden. Het is voldoende om te beschikken over een verdund gas met concurrerende dissipatiemechanismen die afhankelijk zijn van de botsingssnelheid.

Dit heeft belangrijke implicaties voor:

• Fundamentele fysica: Het biedt een minimaal model om te begrijpen hoe niet-lineaire dissipatie kan leiden tot complexe stromingsgedrag zonder de noodzaak van structurele verstopping.
• Toepassingen: Het verklaart mogelijk het gedrag van systemen waar extra dissipatie optreedt bij lage snelheden, zoals geladen korrels, adhezieve deeltjes, of systemen met hydrodynamische interacties.
• Toekomstig onderzoek: De auteurs wijzen erop dat het interessant is om te onderzoeken of dit mechanisme standhoudt in dichter systemen waar collisionele overdracht en ruimtelijke correlaties een rol spelen, en hoe robuust het is tegenover realistischere restitutiewetten en rotatievrijheidsgraden.

Kortom, het artikel onthult een nieuw type kinetische overgang die de link legt tussen microscopische botsingsdynamica en macroscopische discontinuïteiten in de viscositeit.