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

Sintesi Tecnica: Una Nuova Prospettiva sul Flusso Cross-Slot Viscoelastico in Regime di Creeping Flow

Definizione del Problema
Il documento affronta l'origine fisica, da lungo tempo non risolta, dell'instabilità elastica nel flusso cross-slot bidimensionale di fluidi Oldroyd-B in regime di creeping flow. Osservato per la prima volta sperimentalmente quasi mezzo secolo fa e riprodotto numericamente due decenni fa, questo flusso presenta una transizione di rottura della simmetria da uno stato stazionario simmetrico a uno stato stazionario asimmetrico quando il numero di Weissenberg ($Wi$) supera una soglia critica. Nonostante estesi studi fenomenologici, il meccanismo preciso che guida questa instabilità rimane oggetto di dibattito. Le ipotesi precedenti hanno attribuito l'insorgenza a vari fattori, tra cui il flusso compressivo a monte del punto di stagnazione, il flusso estensionale a valle, le singolarità geometriche agli angoli acuti, il forte shear vicino agli angoli o l'interazione tra elasticità e shear-thinning. Gli autori osservano che è mancata una previsione teorica basata sui primi principi dell'insorgenza dell'instabilità e della sua struttura modale, anche a causa delle difficoltà numeriche associate alla simulazione di flussi viscoelastici (il "problema dell'alto numero di Weissenberg") in geometrie complesse prive di stati base analitici.

Metodologia
Gli autori impiegano un approccio multidimensionale completo combinando quattro tecniche analitiche e numeriche per un fluido Oldroyd-B diffusivo in un canale cross-slot 2D standard con numero di Reynolds nullo:

1. Simulazioni Numeriche Dirette (DNS): Utilizzando il metodo degli elementi finiti (FEM) tramite FreeFem++, gli autori risolvono le equazioni di Stokes incomprimibili accoppiate con la legge costitutiva di Oldroyd-B diffusiva. Utilizzano la formulazione del gradiente di stress diviso elastico-viscoso discreto (DEVSS-G) e lo schema streamline upwind Petrov–Galerkin (SUPG) per la stabilizzazione numerica.
2. Analisi di Stabilità Lineare Globale (LSA): Gli autori linearizzano le equazioni governanti attorno allo stato base stazionario per risolvere gli autovalori complessi e gli automodi. Ciò consente di prevedere il $Wi$ critico e i tassi di crescita temporale delle perturbazioni.
3. Analisi della Sensibilità Strutturale: Calcolando sia gli automodi diretti che quelli adimensionali (adjoint), gli autori definiscono i campi di sensibilità strutturale (sia basati sulla velocità che sullo stress) per identificare la regione del "wavemaker" — il nucleo spaziale dove l'instabilità è più sensibile al feedback localizzato.
4. Analisi del Bilancio Energetico: Gli autori decompongono l'equazione dell'evoluzione dell'energia della perturbazione per identificare i meccanismi fisici specifici che contribuiscono alla produzione o alla dissipazione di energia. Si concentrano sul termine di accoppiamento tra lo stato base e le perturbazioni, analizzando specificamente la densità di potenza polimerica per localizzare la sorgente energetica dell'instabilità.

Lo studio investiga anche due varianti geometriche del cross-slot di base: una con angoli arrotondati (per rimuovere le singolarità geometriche) e una con un cilindro centrato (per eliminare il punto di stagnazione libero), al fine di testare la robustezza dei meccanismi identificati.

Risultati Chiave

• Biforcazione e Validazione: La DNS riproduce con successo la canonica biforcazione pitchfork supercritica, con il flusso che transita da simmetrico ad asimmetrico a un numero di Weissenberg critico $Wi_{cr} \approx 0.36$. La LSA predice coerentemente questa soglia e i tassi di crescita temporale, che sono in linea con i dati DNS e i risultati della Decomposizione Modale Dinamica (DMD).
• Struttura del Modo di Instabilità: L'automodo diretto principale consiste in una perturbazione chirale di velocità-stress. La perturbazione di velocità forma una struttura vorticosa a forma di diamante centrata nel punto di stagnazione. Fondamentalmente, la perturbazione dello stress normale polimerico ($\hat{\tau}_{yy}$) presenta una parità antisimmetrica, che inclina e ruota il filamento birifrangente generato dal flusso estensionale. Questa rotazione devia le separatrici del flusso, portando a una distribuzione del flusso diseguale nei rami di uscita.
• Nucleo Spaziale (Wavemaker): L'analisi della sensibilità strutturale rivela che l'instabilità non è guidata dagli angoli o dal punto di stagnazione stesso. Invece, il wavemaker basato sullo stress è localizzato in strette creste ad alto tasso di estensione all'interno del flusso estensionale a valle del punto di stagnazione. Queste creste esibiscono un pattern a "spina di pesce" (herringbone) di perturbazione della conformazione polimerica, caratterizzato da simultaneo stiramento e compressione fianco a fianco.
• Sorgente Energetica: L'analisi del bilancio energetico identifica l'accoppiamento tra lo stato base e la perturbazione come la fonte di energia dominante. Nello specifico, la densità di potenza polimerica della perturbazione ($P$) è positiva e massima lungo le stesse creste verticali identificate dall'analisi del wavemaker. Ciò indica un trasferimento localizzato di energia elastica immagazzinata dalle catene polimeriche deformate verso l'energia cinetica della perturbazione, guidando l'instabilità.
• Varianti Geometriche: In entrambe le varianti (angoli arrotondati e cilindro centrato), l'instabilità persiste. Il picco di densità di potenza polimerica rimane localizzato lungo le creste verticali nelle regioni di flusso estensionale. Ciò dimosta che né gli angoli acuti né un punto di stagnazione centrale libero sono l'ingrediente destabilizzante essenziale.

Significatività e Rivendicazioni
Il documento sostiene di aver risolto la controversia riguardante l'origine fisica dell'instabilità elastica del cross-slot. Gli autori concludono che l'instabilità non deriva dalle singolarità geometriche (angoli) o dal punto di stagnazione in sé, ma piuttosto dal rilascio di energia elastica nelle regioni dominate dall'estensione, caratteristiche del flusso cross-slot.

Nello specifico, il meccanismo prevede:

1. La generazione di un filamento birifrangente con stress normali e gradienti sostanziali nel flusso estensionale a valle del punto di stagnazione.
2. Una perturbazione chirale che inclina e ruota questo filamento.
3. Il rilascio dell'energia elastica immagazzinata in strette creste ad alto tasso di estensione (la regione del wavemaker), che alimenta la crescita dell'instabilità.

Gli autori affermano che i loro risultati riconciliano le precedenti interpretazioni contrastanti: mentre gli angoli acuti non sono necessari (come dimostrato dagli esperimenti OSCER) e il punto di stagnazione non è l'unico driver (come dimostrato dagli esperimenti di inserimento di un cilindro), l'ingrediente essenziale è l'alto tasso di estensione del flusso e il relativo stiramento dei polimeri. Lo studio convalida inoltre l'utilità dell'integrazione di analisi di stabilità globale, sensibilità strutturale e bilancio energetico per comprendere le instabilità in flussi viscoelastici complessi dove non sono disponibili stati base analitici.