Revisiting creeping viscoelastic cross-slot flow: Global linear stability and structural sensitivity analyses

Technische Samenvatting: Een herziening van de kruis-slot visco-elastische stroming bij lage Reynoldsgetallen

Probleemstelling
Het artikel behandelt de langdurig onopgeloste fysieke oorsprong van de elastische instabiliteit in een tweedimensionale kruis-slot stroming (cross-slot flow) van Oldroyd-B vloeistoffen bij een nagenoeg nul Reynoldsgetal. Deze stroming, die bijna een halve eeuw geleden voor het eerst experimenteel werd waargenomen en twee decennia geleden numeriek werd gereproduceerd, vertoont een symmetriebrekende transitie van een stationaire symmetrische staat naar een stationaire asymmetrische staat wanneer het Weissenberggetal ($Wi$) een kritische drempel overschrijdt. Ondanks uitgebreide fenomenologische studies blijft de exacte mechanica die deze instabiliteit aandrijft onderwerp van debat. Eerdere hypothesen hebben de aanvang toegeschreven aan diverse factoren, waaronder compressiestroming stroomopwaarts van het stagnatiepunt, extensiestroming stroomafwaarts, geometrische singulariteiten bij scherpe hoeken, sterke afschuiving nabij de hoeken, of de interactie tussen elasticiteit en afschuifverdunning (shear-thinning). De auteurs merken op dat een theoretische voorspelling vanuit de eerste principes van de onset van de instabiliteit en de modusstructuur tot nu toe ontbrak, deels vanwege de numerieke moeilijkheden geassocieerd met het simuleren van visco-elastische stromingen (het "high Weissenberg number problem") in complexe geometrieën die geen analytische basistoestanden bezitten.

Methodologie
De auteurs hanteren een uitgebreide, veelzijdige aanpak waarbij vier analytische en numerieke technieken worden gecombineerd voor een diffusieve Oldroyd-B vloeistof in een standaard 2D kruis-slot kanaal bij een verwaarloosbaar Reynoldsgetal:

1. Directe Numerieke Simulaties (DNS): Gebruikmakend van de eindelementenmethode (FEM) via FreeFem++, lossen de auteurs de incompressibele Stokes-vergelijkingen op die gekoppeld zijn aan de diffusieve Oldroyd-B constitutieve wet. Zij maken gebruik van de discrete elastic-viscous split-stress gradient (DEVSS-G) formulering en het streamline upwind Petrov–Galerkin (SUPG) schema voor numerieke stabilisatie.
2. Globale Lineaire Stabiliteitsanalyse (LSA): De auteurs lineariseren de regerende vergelijkingen rond de stationaire basistoestand om complexe eigenwaarden en eigenmodi op te lossen. Dit maakt de voorspelling van het kritische $Wi$ en de temporele groeicijfers van perturbaties mogelijk.
3. Structurele Sensitiviteitsanalyse: Door zowel directe als adiabatische (adjoint) eigenmodi te berekenen, definiëren de auteurs structurele sensitiviteitsvelden (zowel op basis van snelheid als op basis van spanning) om de "wavemaker"-regio te identificeren: de ruimtelijke kern waar de instabiliteit het meest gevoelig is voor lokale feedback.
4. Energie-budgetanalyse: De auteurs deconstrueren de evolutie van de perturbatie-energievergelijking om de specifieke fysieke mechanismen te identificeren die bijdragen aan energieproductie of -dissipatie. Zij richten zich op de koppelingsterm tussen de basistoestand en de perturbaties, waarbij zij specifiek de polymere vermogensdichtheid analyseren om de energetische bron van de instabiliteit te lokaliseren.

De studie onderzoekt ook twee geometrische varianten van de standaard kruis-slot: één met afgeronde hoeken (om geometrische singulariteiten te verwijderen) en één met een gecentreerde cilinder (om het vrije stagnatiepunt te elimineren), om de robuustheid van de geïdentificeerde mechanismen te testen.

Belangrijkste Resultaten

• Bifurcatie en Validatie: De DNS reproduceert succesvol de canonieke superkritische pitchfork-bifurcatie, waarbij de stroming overgaat van symmetrisch naar asymmetrisch bij een kritisch Weissenberggetal $Wi_{cr} \approx 0,36$. De LSA voorspelt deze drempel en de temporele groeicijfers consistent, welke overeenkomen met de DNS-data en Dynamic Mode Decomposition (DMD) resultaten.
• Instabiliteitsmodusstructuur: De leidende directe eigenmode bestaat uit een chirale snelheid-spanning-perturbatie. De snelheidspersubatie vormt een ruitvormige vortexstructuur gecentreerd rond het stagnatiepunt. Cruciaal is dat de normale polymere spanning-perturbatie ($\hat{\tau}_{yy}$) een antisymmetrische pariteit vertoont, die de dubbelbrekende streng (birefringent strand) gegenereerd door de extensiestroming doet kantelen en roteren. Deze rotatie buigt de stromingsseparatrices af, wat leidt tot een ongelijke fluxverdeling in de uitlaatvertakkingen.
• Ruimtelijke Kern (Wavemaker): Structurele sensitiviteitsanalyse onthult dat de instabiliteit niet wordt gedreven door de hoeken of het stagnatiepunt zelf. In plaats daarvan is de spanning-gebaseerde wavemaker gelokaliseerd in smalle, hoog-extensiesnelheids-richels binnen de extensiestroming stroomafwaarts van het stagnatiepunt. Deze richels vertonen een "herringbone"-patroon van polymeerconfiguratie-perturbatie, gekenmerkt door simultane zij-aan-zij rek en compressie.
• Energetische Bron: Energie-budgetanalyse identificeert de koppeling tussen de basistoestand en de perturbatie als de dominante energiebron. Specifiek is de verstoring van de polymere vermogensdichtheid ($P$) positief en maximaal langs dezelfde verticale richels die geïdentificeerd zijn door de wavemaker-analyse. Dit duidt op een gelokaliseerde overdracht van opgeslagen elastische energie van de vervormde polymeerketens terug naar de kinetische energie van de verstoring, wat de instabiliteit aandrijft.
• Geometrische Varianten: In zowel de variant met afgeronde hoeken als de variant met een gecentreerde cilinder blijft de instabiliteit bestaan. De piek in de polymere vermogensdichtheid blijft gelokaliseerd langs verticale richels in de extensiestromingsgebieden. Dit demonstreert dat noch scherpe hoeken, noch een vrij centraal stagnatiepunt de essentiële destabiliserende ingrediënten zijn.

Betekenis en Claims
Het artikel claimt de controverse over de fysieke oorsprong van de kruis-slot elastische instabiliteit te hebben opgelost. De auteurs concluderen dat de instabiliteit niet voortkomt uit geometrische singulariteiten (hoeken) of het stagnatiepunt per se, maar uit de vrijgave van elastische energie in extensie-gedomineerde regio's die kenmerkend zijn voor de kruis-slot stroming.

Specifiek omvat het mechanisme:

1. De generatie van een dubbelbrekende streng met aanzienlijke normale spanningen en gradiënten in de extensiestroming stroomafarts van het stagnatiepunt.
2. Een chirale perturbatie die deze streng doet kantelen en roteren.
3. De vrijgave van opgeslagen elastische energie in smalle, hoog-extensiesnelheids-richels (de wavemaker-regio), die de groei van de instabiliteit voedt.

De auteurs stellen dat hun bevindingen eerdere conflicterende interpretaties verzoenen: hoewel scherpe hoeken niet noodzakelijk zijn (zoals aangetoond door OSCER-experimenten) en het stagnatiepunt niet de enige drijfveer is (zoals aangetoond door experimenten met cilinder-insertie), is het essentiële ingrediënt de hoog-extensiesnelheidsstroming en de daarmee gepaard gaande rek van de polymeren. De studie valideert verder het nut van het integreren van globale stabiliteit, structurele sensitiviteit en energie-budgetanalyses voor het begrijpen van instabiliteiten in complexe visco-elastische stromingen waar analytische basistoestanden niet beschikbaar zijn.