A Solid-Based Approach for Modeling Simple Yield-Stress Fluids: Rheological Transitions, Overshoot and Relaxation

In deze studie wordt een nieuw constitutief model voor eenvoudige vloeistoffen met vloeigrens voorgesteld dat, gebaseerd op een visco-elastisch solide element, zowel stationaire als transienten rheologische gedragingen zoals spanningsoverschrijding tijdens opstartende schuifstroom nauwkeurig voorspelt en aantoont dat dit overschot het gevolg is van een homogeen mechanisme waarbij normaalspanningsverschillen de plastische respons versnellen.

De kern

Een Nieuwe Manier om "Vloeibare Vaste Stoffen" te Begrijpen: Een Simpele Uitleg

Stel je voor dat je een pot honing hebt die zich soms als een stevige gelatine gedraagt, maar als je er hard op duwt, plotseling als water vloeit. Of denk aan tandpasta: hij blijft op je tandenborstel zitten (vast), maar stroomt eruit als je knijpt (vloeibaar). In de wetenschap noemen we dit vloeistoffen met een vloeigrens (yield-stress fluids). Ze zijn overal: van modderige rivieren tot de inkt die wordt gebruikt om 3D-objecten te printen of batterijen te maken.

Deze nieuwe studie van onderzoekers uit Zuid-Korea probeert een wiskundige formule te vinden die precies voorspelt hoe deze rare materialen zich gedragen. Hier is de kern van hun ontdekking, vertaald naar alledaags taal:

1. Het Probleem: De "Verwarde" Vloeistof

Vroeger dachten wetenschappers dat ze deze materialen konden beschrijven met simpele regels voor vloeistoffen of vaste stoffen. Maar dat werkt niet goed.

• Het Overshoot-probleem: Als je zo'n materiaal plotseling begint te roeren, stijgt de weerstand eerst heel snel (een piek) en daalt daarna pas naar een normaal niveau. Alsof je een deur opent die eerst vastzit, dan met een klap open gaat, en daarna soepel draait. Oude modellen konden dit "pieken" niet verklaren.
• Het Relaxatie-probleem: Als je stopt met roeren, verdwijnt de kracht niet helemaal tot nul. Er blijft een beetje spanning over, alsof het materiaal nog steeds "onthoudt" dat het bewogen is.

2. De Oplossing: Een Nieuw Bouwpakket

De onderzoekers hebben een nieuw model bedacht dat deze materialen ziet als een mix van een elastisch net en een vloeibaar solvent.

Stel je het materiaal voor als een dichte menigte mensen op een drukke markt (de microgels):

• Het Elastische Net (De Gel): De mensen houden elkaar vast. Als je ze een beetje duwt, buigen ze mee en komen ze terug als je stopt (zoals een veer). Dit is het "vaste" deel.
• De Vloeibare Tussenruimte (Het Oplosmiddel): Tussen de mensen zit water. Als de mensen zich verplaatsen, moet dit water ook meebewegen, wat weerstand biedt. Dit is het "vloeibare" deel.

Het nieuwe model combineert deze twee delen op een slimme manier:

1. Een veer (voor de elasticiteit).
2. Een demper (zoals een schokdemper in een auto, voor de plastische, permanente vervorming).
3. Een extra demper die parallel loopt, voor het water dat tussen de deeltjes zit.

3. De Grote Doorbraak: Waarom de "Pieken" Ontstaan

Het meest interessante deel van dit onderzoek is hoe ze de "overshoot" (die plotselinge piek in weerstand) verklaren.

Vroeger dachten mensen dat deze piek kwam door complexe veranderingen in de structuur van het materiaal (zoals mensen die in groepjes gaan staan). Maar dit onderzoek toont aan dat het geen ingewikkelde structuurveranderingen nodig heeft.

De Analogie:
Stel je voor dat je een groep mensen in een krappe gang duwt.

• In het begin duw je ze recht vooruit. Ze buigen mee (elastisch).
• Maar door de druk beginnen ze ook een beetje naar de zijkanten te duwen (dit noemen ze "normale spanning").
• Deze zijwaartse duwkracht maakt dat de mensen sneller uit elkaar geduwd worden dan je verwachtte.
• Hierdoor daalt de weerstand plotseling, wat de "piek" veroorzaakt.

Het model laat zien dat dit een homogeen effect is: het gebeurt overal tegelijk, puur door de wiskundige manier waarop de krachten in drie dimensies werken. Je hoeft niet aan te nemen dat het materiaal "moe" wordt of dat de deeltjes zich anders gaan ordenen. Het is puur een gevolg van de fysica van het duwen in 3D.

4. Wat Betekent Dit voor de Wereld?

Deze formule is als een nieuwe "GPS" voor ingenieurs.

• Batterijen: Bij het maken van lithium-batterijen wordt een modderige pasta (elektrode-suspensie) gebruikt. Als je deze te snel of te traag verwerkt, kan de structuur kapot gaan. Met dit model kunnen fabrikanten precies voorspellen hoe de pasta zich gedraagt tijdens het storten of printen.
• 3D-Printen: Bij het printen met inkt moet de inkt eerst vast zijn om zijn vorm te houden, maar vloeibaar zijn om door de printkop te gaan. Dit model helpt om de perfecte "vloeigrens" te vinden.
• Algemene Toepassing: Het model werkt niet alleen voor simpele bewegingen, maar ook voor complexe, driedimensionale situaties (zoals wanneer je een materiaal in alle richtingen verdraait).

Conclusie

Kortom, deze onderzoekers hebben een nieuwe wiskundige taal ontwikkeld die zegt: "Deze rare vloeistoffen zijn eigenlijk vaste stoffen die kunnen vloeien, en hun vreemde gedrag (zoals pieken en niet helemaal ontspannen) komt door de manier waarop ze in drie dimensies op elkaar drukken."

Het is een stap in de richting van het beter begrijpen en ontwerpen van materialen die ons dagelijks leven makkelijker maken, van onze tandpasta tot de batterij in je telefoon.

Technische samenvatting

Titel: Een op vaste stof gebaseerde aanpak voor het modelleren van eenvoudige vloeistoffen met vloeigrens: Reologische overgangen, overshoot en relaxatie

1. Probleemstelling

Vloeistoffen met een vloeigrens (yield-stress fluids), zoals Carbopol-dispersies en elektrode-suspensies voor batterijen, vertonen complex reologisch gedrag dat niet volledig kan worden beschreven door alleen stationaire eigenschappen (zoals de Herschel-Bulkley-relatie). Bestaande constitutieve modellen hebben moeite om zowel stationaire als transiënte gedragingen tegelijkertijd nauwkeurig te voorspellen. Specifiek zijn er vier uitdagingen die vaak niet door één enkel model worden opgelost:

1. Stress overshoot: Het verschijnsel waarbij de schuifspanning tijdelijk een piek bereikt tijdens het opstarten van een schuifstroom, afhankelijk van de stroomcondities.
2. Herschel-Bulkley-gedrag: De relatie tussen spanning en schuifsnelheid in stationaire toestand.
3. Niet-nul relaxatiespanning: Na het stoppen van de stroom (stress relaxation) daalt de spanning niet volledig naar nul, maar bereikt een eindige, niet-nul waarde (karakteristiek voor vaste stoffen).
4. Creep-overgang: De overgang van vast-achtig gedrag (krimp saturatie) naar vloeistof-achtig gedrag (continue krimp) afhankelijk van de aangebrachte spanning.

Veel bestaande modellen attribueren stress overshoot aan thixotropie (tijdsafhankelijke microstructuurverandering) of isotrope hardening. Echter, experimenten met eenvoudige vloeistoffen met verwaarloosbare thixotropie (zoals Carbopol) tonen overshoot aan, wat suggereert dat een ander mechanisme verantwoordelijk is. Daarnaast voorspellen sommige modellen geen relaxatie naar een niet-nul spanning, wat in strijd is met recente experimentele waarnemingen.

2. Methodologie

De auteurs stellen een nieuwe constitutieve vergelijking voor die gebaseerd is op een visco-elastisch vastestofmodel in plaats van een visco-elastisch vloeistofmodel. De aanpak omvat de volgende stappen:

• Fysisch model: Het gedrag wordt gemodelleerd als een verzadigd systeem van microgel-deeltjes (gel-fase) met daartussen een interstitieel oplosmiddel.
  • De gel-fase wordt beschreven door een Zener-model (een veer in serie met een Kelvin-Voigt-element). De veer vertegenwoordigt de elastische vervorming binnen de deeltjes, en het Kelvin-Voigt-element (veer + demper) vertegenwoordigt de plastische herschikking van de aggregaten.
  • De oplosmiddelfase wordt gemodelleerd als een lineaire demper parallel aan het Zener-element, wat zorgt voor extra viskeuze dissipatie.
• Mechanische decompositie:
  • Voor het 1D-model (simple shear) wordt een additieve decompositie van de rek gebruikt ($\gamma = \gamma_e + \gamma_p$).
  • Voor het 3D-model (algemene vervorming) wordt de Kröner-Lee decompositie ($F = F^e \cdot F^p$) toegepast om de principes van materiële frame-invariantie te waarborgen bij eindige vervormingen.
• Constitutieve componenten:
  • Vloeibaarheidsmodel: Er worden twee niet-lineaire viscositeitsmodellen gebruikt: het Eyring-model en het Carreau-Yasuda-model. Deze zijn afhankelijk van de spanningsinvariantie en beschrijven de schuifverdunning.
  • Vloeiregel: De plastische vervormingssnelheid wordt bepaald door een vloeiregel die de spanning in het Kelvin-Voigt-element relateert aan de plastische snelheid.
  • Back-stress evolutie: De evolutie van de back-stress (de spanning in de veer van het Kelvin-Voigt-element) volgt de Armstrong-Frederick-theorie, wat kinematische harding beschrijft zonder isotrope harding.

3. Belangrijkste Bijdragen

• Unificatie van vastestof- en vloeistofkarakteristieken: Het model behandelt yield-stress fluids fundamenteel als vaste stoffen die kunnen vloeien, wat verklaart waarom ze een niet-nul relaxatiespanning vertonen.
• Voorspelling van overshoot zonder thixotropie: Het model toont aan dat stress overshoot tijdens opstart-schuiven kan worden voorspeld puur door de tensoriële structuur van het 3D-model en de interactie tussen elastische belasting en relaxatie, zonder de noodzaak van thixotrope evolutie of ruimtelijke heterogeniteit.
• 3D Uitbreiding: De auteurs hebben het model succesvol uitgebreid van een scalair 1D-model naar een tensoriële 3D-formulering die geldig is voor complexe stromingen en voldoet aan de principes van continuümmechanica.

4. Resultaten

De resultaten werden gevalideerd via analytische afleidingen en numerieke integratie voor drie testscenario's:

• Opstart-schuiven (Start-up Shear):
  • Het 1D-model voorspelde een monotoon toenemende spanning naar een stationaire waarde, maar geen overshoot.
  • Het 3D-model voorspelde wel een duidelijke stress overshoot. De oorzaak bleek te liggen in de tensoriële aard van de spanning: normale spanningsverschillen verhogen de spanningsinvariantie, wat de plastische reactie versnelt en tijdelijk de elastische belasting overtreft. Dit mechanisme is homogeen en vereist geen microstructuur-heterogeniteit.
• Stress Relaxatie:
  • Na het stoppen van de schuifstroom daalt de spanning en bereikt een niet-nul, volledig gerelaxeerde spanning.
  • In het 3D-model is deze relaxatiespanning afhankelijk van de voorafgaande schuifsnelheid (hoger voorafgaand schuiven leidt tot lagere relaxatiespanning), wat overeenkomt met experimentele data. Het 1D-model voorspelde een relaxatiespanning die onafhankelijk was van de snelheid.
• Creep (Krimp):
  • Het model reproduceerde de overgang van vast-achtig gedrag (rek saturatie) naar vloeistof-achtig gedrag (continue rek) wanneer de aangebrachte spanning de vloeigrens overschrijdt.
  • De drempelwaarde voor deze overgang correspondeerde goed met de vloeigrens ($\tau_y$) afgeleid uit de Herschel-Bulkley-fit.
• Stationaire Stroom:
  • Bij lage schuifsnelheden volgt het model de Herschel-Bulkley-relatie. Bij hoge snelheden vertoont het een tweede Newtoniaans plateau, veroorzaakt door de parallelle demper (oplosmiddel).

5. Betekenis en Conclusie

De studie biedt een robuust theoretisch raamwerk voor het modelleren van yield-stress fluids zoals Carbopol en batterij-suspensies. De belangrijkste implicaties zijn:

1. Fundamenteel inzicht: Stress overshoot is een inherent gevolg van de visco-elastische vaste-stofkarakteristieken en de tensoriële koppeling in 3D, niet noodzakelijk een teken van thixotropie.
2. Praktische toepasbaarheid: Het model kan worden gebruikt in numerieke simulaties (zoals FEM) voor industriële processen waar complexe vervormingen optreden, zoals bij het gieten van batterij-elektroden of direct ink writing.
3. Verbeterde voorspellingskracht: In tegenstelling tot eerdere modellen, slaagt dit model erin om de volledige reeks van reologisch gedrag (overshoot, relaxatie naar niet-nul spanning, en creep-overgang) in één samenhangend kader te beschrijven.

De auteurs concluderen dat hun 3D constitutieve vergelijking een solide basis vormt voor kwantitatieve analyses van complexe stromingen van yield-stress fluids, met de mogelijkheid om thixotropie in de toekomst als een extra component toe te voegen voor nog complexere materialen.