Shear Banding in Simulations of Polymer Melts

Deze studie vergelijkt numerieke simulaties van polymerschuifstroming met een Rolie-Poly-model dat Convected Constraint Release (CCR) incorporeert, en concludeert dat er een semi-kwantitatieve overeenkomst is waarbij schuifbanding optreedt wanneer het evenwichtsverstrengelingsgetal een kritieke waarde overschrijdt die afhangt van de krommingstijfheid.

De kern

Scheuren in een soep van plastic: Waarom sommige polymeren in banen stromen

Stel je voor dat je een grote pan met heel dik, plakkerig honing of gesmolten plastic (een polymeer) hebt. Als je deze soep rustig roert, stroomt hij soepel en gelijkmatig. Maar wat gebeurt er als je heel snel en krachtig roert? Soms gedraagt het zich raar: in plaats van dat alles even snel beweegt, splitst de soep zich op in verschillende "banen". Er ontstaan snelle stromen naast trage stromen. Dit fenomeen noemen wetenschappers schuifbanding (shear banding).

In dit onderzoek kijken drie wetenschappers (Lucas, Gary en Peter) naar waarom dit gebeurt in gesmolten plastic, en of ze dit kunnen voorspellen met een wiskundig model. Ze doen dit door te kijken naar computer-simulaties, want in het echte leven is dit heel lastig te meten zonder dat de meting zelf de resultaten verstoort.

Hier is de uitleg in simpele taal, met wat creatieve vergelijkingen:

1. De Verwikkeldheid: Een bos van spaghetti

Stel je de moleculen in het plastic voor als lange, slappe slierten spaghetti. Als je een potje spaghetti hebt, liggen ze allemaal door elkaar.

• Verwikkeldheid (Entanglements): Als de slierten heel lang zijn, raken ze enorm verstrikt. Ze kunnen niet zomaar langs elkaar heen glijden. Ze zitten vast in een soort "buis" van andere slierten.
• Het getal Z: Dit is een maat voor hoe verstrikt de spaghetti is. Een klein getal betekent een paar verstrengelingen; een groot getal betekent een enorme, ondoordringbare knoop.

2. De "Borstel" en het "Loslaten" (CCR)

Wanneer je de pan snel roert (shear flow), proberen de spaghetti-slierten zich uit te strekken en uit te lijnen met de stroming.

• Het probleem: Als ze te snel worden uitgerekt, beginnen ze de verstrengelingen te "schuiven" of los te laten. Dit noemen ze CCR (Convected Constraint Release).
• De analogie: Stel je voor dat je door een drukke menigte loopt. Als je rustig loopt, moet je wachten tot mensen opzij gaan (verstrengelingen). Maar als je hard rent, duw je mensen letterlijk opzij en maken ze ruimte voor je. Je komt sneller vooruit, maar de "buis" waarin je liep, is nu minder strak.
• De parameter β: Dit is een knop die bepaalt hoe makkelijk de spaghetti's loslaten als je hard roert. Sommige soorten plastic (zoals stijve staafjes) laten makkelijker los dan andere (zoals zachte slierten).

3. De Voorspelling: Wanneer breekt het?

De wetenschappers gebruiken een wiskundig model (het DO-model) om te voorspellen wanneer deze spaghetti-soep in banen zal splitsen.

• De theorie: Het model zegt: "Als de spaghetti's erg verstrikt zijn (hoog Z-getal) én ze niet te snel loslaten (kleine β), dan gaat het mis."
• Het resultaat: De soep wordt onstabiel. Het kan niet meer gelijkmatig stromen. Het splitst zich op in een snelle strook (waar de spaghetti's loslaten en makkelijk glijden) en een trage strook (waar ze nog verstrikt zitten).
• De grens: Er is een kritiek punt. Als je plastic niet verstrikt genoeg is, stroomt het gewoon. Als het wel verstrikt genoeg is, krijg je die rare banen.

4. De Simulatie: De Digitale Keuken

De auteurs draaiden enorme computer-simulaties om dit te testen. Ze bouwden virtuele spaghetti's en keken wat er gebeurde als ze deze virtueel roerden.

• Wat zagen ze? Ja, ze zagen de banen! De computer-simulaties lieten precies zien wat het model voorspelde: bij bepaalde snelheden en bij bepaalde soorten plastic splitste de stroming op.
• De overeenkomst: De simulaties kwamen heel goed overeen met de wiskundige voorspellingen. Het model was dus een goede "recept" voor het gedrag van deze plastic soep.

5. Waarom zien we dit niet altijd in het echt?

Je vraagt je misschien af: "Als dit zo duidelijk is in de computer, waarom zien we het niet in de fabriek of in een flesje verf?"

• De "Stijfheid" van het plastic: Het onderzoek toont aan dat het type plastic belangrijk is. Sommige plasticsoorten (zoals die met een specifieke chemische structuur) hebben een heel hoge "loslaat-knop" (β). Hierdoor laten ze zo snel los dat ze nooit in de gevaarlijke zone komen waar banding ontstaat.
• De "Grenzen" van de theorie: Het model werkt goed, maar niet perfect. In de echte wereld zijn er extra factoren, zoals wrijving tegen de wanden van de bak of breuken aan de randen, die de banen kunnen verstoren of verbergen.
• De geschiedenis: Soms hangt het gedrag af van hoe je begon. Als je de plastic soep eerst heel snel roert en dan vertraagt, gedraagt hij zich anders dan als je langzaam begint en versnelt. Het is alsof de spaghetti's een "herinnering" hebben aan hoe je ze hebt behandeld.

Conclusie: Wat leren we hieruit?

Dit onderzoek is als het vinden van de perfecte recept voor het bakken van een taart, maar dan voor vloeibaar plastic.

1. Het model werkt: We kunnen nu redelijk goed voorspellen welk plastic in banen zal stromen en welk niet, op basis van hoe verstrikt het is en hoe stijf de moleculen zijn.
2. De oorzaak: Het fenomeen wordt veroorzaakt door een strijd tussen de verstrengeling (die de soep dik maakt) en het loslaten van die verstrengeling door de beweging.
3. Toekomst: Door dit te begrijpen, kunnen ingenieurs beter voorspellen hoe plastic zal stromen in productieprocessen. Misschien kunnen ze zelfs materialen ontwerpen die niet in banen stromen, zodat ze makkelijker te verwerken zijn, of juist materialen die wel in banen stromen voor specifieke toepassingen.

Kortom: Het is een succesvol experiment om te begrijpen waarom sommige plasticsoepen "kapot" gaan in de stroming, en het helpt ons om de complexe dans van moleculen beter te doorgronden.

Technische samenvatting

Titel: Shear Banding in Simulaties van Polymeersmelten

Auteurs: Lucas L. Nelson, Gary S. Grest en Peter D. Olmsted.

---

1. Het Probleem

Polymeervloeistoffen vertonen vaak complexe stromingsgedragingen wanneer ze onder hoge schuifspanning staan. Een van de meest gecontesteerde fenomenen is shear banding (schuifbandvorming), waarbij de vloeistof in stationaire toestand inhomogene gebieden met verschillende schuifsnelheden vormt.

• Experimentele onduidelijkheid: In experimenten met polymeersmelten is het bewijs voor shear banding vaak onbeslist of afwezig. Waarnemingen worden vaak verstoord door randeffecten (zoals "edge fracture") en wandglijding, wat het moeilijk maakt om een zuivere stationaire toestand te meten.
• Theoretische uitdaging: Bestaande modellen, zoals het Doi-Edwards (DE) model, voorspellen een niet-monotone relatie tussen schuifspanning en schuifsnelheid, wat leidt tot bandvorming. Echter, om dit te corrigeren voor realistischere polymeren, werd het mechanisme van Convected Constraint Release (CCR) geïntroduceerd. CCR versnelt de ontvlechting van polymeren door stroming, wat de spanning verhoogt en de niet-monotonie (en dus bandvorming) kan onderdrukken. De sterkte van dit effect wordt bepaald door een parameter $\beta$.
• Kernvraag: Onder welke omstandigheden (karakteristieke parameters zoals het aantal vervlechtingen $Z$ en de CCR-parameter $\beta$) treedt shear banding op in polymeersmelten, en waarom wordt dit in experimenten soms niet waargenomen?

2. Methodologie

De auteurs combineren theoretische modellering met uitgebreide moleculair-dynamica (MD) simulaties om deze vraag te beantwoorden.

• Theoretisch Model (DO-model):

  • Ze gebruiken het Dolata-Olmsted (DO) model, een thermodynamisch consistente uitbreiding van het Rolie-Poly model.
  • Dit model koppelt de evolutie van de buisconformatie (tube conformation) aan het verminderde aantal primitieve pad-knikken ($\nu$), wat de mate van ontvlechting (disentanglement) door CCR beschrijft.
  • Het model bevat vijf dimensieloze parameters: het aantal vervlechtingen $Z$, de CCR-parameter $\beta$, de relaxatie-anisotropie $\alpha$, de maximale rek $\lambda_{max}$ en de viscositeitsverhouding $\varepsilon$.
  • De auteurs berekenen een "fasediagram" dat aangeeft wanneer de constitutieve kromme (spanning vs. schuifsnelheid) monotoon of niet-monotoon is. Niet-monotonie is een voorwaarde voor shear banding.

• Simulaties:

  • Er zijn non-equilibrium moleculair-dynamica (NEMD) simulaties uitgevoerd met Kremer-Grest (KG) bol-veer polymeren.
  • Twee verschillende methoden werden gebruikt:
    1. SLLOD: Dwingt een uniforme snelheidsprofiel op (onderdrukt bandvorming, nuttig voor het meten van $\beta$).
    2. Deformeerde cel (DPD thermostaat): Staat inhomogene stroming toe en kan shear banding vormen.
  • Variatie in parameters: De stijfheid van de ketens werd gevarieerd via een hoekpotentiaal ($k_\theta = 0, 1.5, 2.0$) om het effect op $\beta$ te bestuderen. Ketenglengten varieerden van 200 tot 2000 monomeren.
  • Data-analyse: Het aantal vervlechtingen ($Z_k$) werd gemeten met de Z1+ code om de ontvlechting $\nu$ te kwantificeren als functie van de Weissenberg-getallen.

3. Belangrijkste Bijdragen

1. Kwantificering van $\beta$: De auteurs hebben een methode ontwikkeld om de CCR-parameter $\beta$ empirisch te bepalen uit simulaties door de voorspelde ontvlechting (via het DO-model) te fiten aan de gemeten ontvlechting in de SLLOD-simulaties.
2. Relatie Stijfheid en $\beta$: Ze tonen aan dat $\beta$ sterk afhankelijk is van de ketenstijfheid (gekwantificeerd als de verhouding Kuhn-lengte tot verpakkingslengte, $b_K/p$). Stijvere polymeren hebben een hogere $\beta$.
3. Validatie van het Fasediagram: Ze vergelijken de voorspelde grenzen voor bandvorming uit het DO-model direct met de waarnemingen in de simulaties voor verschillende waarden van $Z$ en $\beta$.

4. Resultaten

• Overeenkomst Model en Simulatie: Er is een semi-kwantitatieve overeenkomst tussen het DO-model en de simulaties. Het model voorspelt correct dat shear banding optreedt wanneer het aantal vervlechtingen $Z$ een kritieke waarde $Z_c$ overschrijdt.
• Invloed van $\beta$:
  • Een hogere $\beta$ (sterkere CCR) zou volgens het model de kritieke waarde $Z_c$ verhogen, waardoor bandvorming moeilijker optreedt.
  • De simulaties tonen echter een tegenstrijdigheid: voor de onderzochte stijfheden lijkt $Z_c$ te afnemen bij toenemende $\beta$, wat tegengesteld is aan de theorie. De auteurs speculeren dat dit komt door de beperkingen van het buis-model bij hoge schuifsnelheden of door methodologische problemen bij het meten van $\beta$ in sterk uitgelijnde stromingen.
• Stijfheid en Bandvorming:
  • Flexibele polymeren (lage $k_\theta$, lage $\beta$) vertonen shear banding bij lagere $Z$-waarden.
  • Stijvere polymeren (hoge $k_\theta$, hoge $\beta$) vertonen minder of geen bandvorming in de simulaties, wat overeenkomt met het idee dat sterke CCR de instabiliteit onderdrukt.
• Normale Spanningsverhouding: De verhouding $-N_2/N_1$ (gerelateerd aan $\alpha$) is moeilijk nauwkeurig te bepalen in experimenten en simulaties bij lage schuifsnelheden. De data suggereert dat $\alpha$ voor smelten mogelijk hoger is dan voor oplossingen, wat de voorspellingen van het model beïnvloedt.
• Geschiedenisafhankelijkheid: De vorming van banden is sterk afhankelijk van de stromingsgeschiedenis (metastabiliteit). Een systeem dat vanuit rust wordt gestart, kan banden vormen, terwijl hetzelfde systeem dat vanuit een hoge schuifsnelheid wordt afgekoeld, uniform kan blijven.

5. Betekenis en Conclusie

• Interpretatie van Experimenten: De resultaten suggereren dat het ontbreken van shear banding in experimenten met polymeersmelten mogelijk te wijten is aan de chemische specificiteit van de polymeren. Polymeren met complexe interne vrijheidsgraden (zoals dihedral-rotaties in echte polymeren) kunnen een zeer hoge $\beta$ hebben.
• Kritieke Vervlechting: Voor polymeren met een hoge $\beta$ (zoals polyetheen) zou de kritieke waarde $Z_c$ voor bandvorming extreem hoog kunnen zijn (bijv. $Z > 20-300$). Aangezien veel experimentele smelten $Z$-waarden onder deze drempel hebben, wordt bandvorming daar niet waargenomen.
• Voorspelling: Flexibele polymeren met een lage $\beta$ (zoals PDMS) zouden volgens dit model veel eerder shear banding moeten vertonen dan stijvere polymeren.
• Beperkingen: Het DO-model (en buis-modellen in het algemeen) faalt bij zeer hoge schuifsnelheden waar de vervlechtingen volledig verdwijnen en de "buis" zelf niet meer gedefinieerd is. Dit verklaart de afwijkingen bij hoge schuifsnelheden in de simulaties.

Kortom, dit werk biedt een brug tussen microscopische simulaties en macroscopische rheologische modellen, en biedt een verklaring voor waarom shear banding in sommige polymeersystemen wel en in andere niet wordt waargenomen, gebaseerd op de chemische stijfheid en de sterkte van de convectieve ontvlechting.