Entropy-Compatible Reconstruction for High-Weissenberg Viscoelastic Flow

Sintesi Tecnica: Ricostruzione Compatibile con l'Entropia per Flussi Viscoelastici ad Alto Numero di Weissenberg

Definizione del Problema
Il problema dell'alto numero di Weissenberg rimane un ostacolo centrale nella simulazione dei flussi viscoelastici. Nei modelli a tensore di conformazione, lo stato polimerico è rappresentato da un tensore simmetrico definito positivo (SPD) $A$. Sebbene la positività sia essenziale per il significato fisico dell'energia libera elastica ($\log \det A$) e del logaritmo di matrice, l'articolo sostiene che la positività da sola è insufficiente per la stabilità numerica ad alti numeri di Weissenberg.

I metodi esistenti, come la log-conformazione ($A = \exp \Psi$) e le ricostruzioni radice quadrata ($A = BB^T$), riescono a imporre il vincolo SPD. Tuttavia, gli esperimenti numerici mostrano che questi metodi possono comunque produrre overshoot di stress non fisici e una crescita non fisica dell'energia elastica. Il problema principale identificato è la disconnessione tra il vincolo di cono (positività) e il vincolo di bilancio accoppiato (compatibilità con l'energia libera). Nello specifico, se il tensore utilizzato nel lavoro di stress dell'equazione del momento differisce dal tensore utilizzato nel lavoro di stretching dell'equazione di conformazione o nelle variabili di entropia, si genera un difetto di lavoro a segno indefinito. Inoltre, le ricostruzioni strettamente positive possono introdurre tre meccanismi non fisici specifici:

1. Bias entropico di tipo Jensen: le mappe non lineari (come l'esponenziale) applicate a perturbazioni a media nulla nello spazio logaritmico generano un bias positivo della media nel tensore fisico e nell'entropia.
2. Amplificazione esponenziale: in stati altamente stirati, piccole perturbazioni logarithmiche vengono amplificate dall'esponenziale di matrice, scalando con il valore proprio massimo $\lambda_{\max}(A)$.
3. Difetti di lavoro a segno indefinito: l'uso di tensori ricostruiti discordanti per i termini di stress e di entropia impedisce la cancellazione esatta del lavoro polimerico, portando a un'iniezione o rimozione incontrollata di energia.

Metodologia
L'articolo propone una ricostruzione logaritmica corretta e compatibile con l'entropia. Piuttosto che sostituire i metodi esistenti di preservazione della positività, questo approccio agisce come un filtro di ammissibilità applicato a una ricostruzione logaritmica ad alto ordine grezza prima che entri nel bilancio accoppiato stress-lavoro.

La metodologia prevede:

• Definizione del Percorso Logaritmico: Dato un predittore fisico $\hat{A}$ e una ricostruzione ad alto ordine grezza $\tilde{A}$, viene definito un percorso nello spazio logaritmico: $A(\theta) = \exp(\hat{\Psi} + \theta(\tilde{\Psi} - \hat{\Psi}))$ per $\theta \in [0, 1]$.
• Vincolo del Budget di Entropia: Il metodo cerca il valore massimo di $\theta$ tale che l'energia libera discreta totale (entropia) non superi l'entropia del predittore fisico più un budget locale $\tau_h$. Il budget è scalato per corrispondere all'errore di troncamento atteso dello schema ad alto ordine (es. $O(h^{2k+2})$).
• Algoritmo di Bisezione: l'ottimo $\theta$ viene trovato tramite una ricerca locale per bisezione cella per cella. Se la ricostruzione grezza soddisfa il budget, $\theta=1$ (nessuna modifica). Altrimenti, il metodo torna indietro lungo il percorso fino al punto massimo ammissibile.
• Imposizione della Consistenza: Il tensore selezionato è utilizzato coerentemente in tutti i termini di accoppiamento: lavoro di stress del momento, stretching di conformazione, variabili di entropia e quadratura.

Contributi Chiave e Risultati Teorici
L'articolo fornisce un'analisi rigorosa di questo principio di ricostruzione, dimostrando diverse proprietà chiave:

• Esistenza e Convessità: il funzionale dell'entropia lungo il percorso logaritmico è convesso, garantendo l'esistenza di un parametro ammissibile massimo $\theta^*$.
• Minimo Smorzamento: tra tutti gli stati ammissibili sul percorso, il $\theta^*$ selezionato minimizza la distanza dalla ricostruzione grezza, rendendo la correzione "a minimo smorzamento".
• Controllo Spettrale: la correzione preserva l'indifferenza al frame e mantiene i limiti spettrali derivati dai logaritmi degli estremi.
• Inattività Asintotica: se la ricostruzione grezza ad alto ordine è già entropicamente ammissibile (ovvero, il suo difetto è entro la scala dell'errore di discretizzazione), la correzione diventa inattiva ($\theta^* = 1$), preservando l'accuratezza formale di alto ordine.
• Vantaggio di Risoluzione in Regimi di Alto Stiramento: nei regimi di alto stiramento, il metodo mantiene un'elevata accuratezza nelle metriche log-relative e di entropia, mentre le ricostruzioni dirette nello spazio fisico soffrono di problemi di scaling delle derivate proporzionali a $\lambda_{\max}(A)$.
• Stima dell'Energia: gli autori dimostrano una stima dell'energia libera discreta per il sistema accoppiato, mostrando che il budget di ricostruzione $\tau_h$ agisce come una sorgente di entropia controllata e limitata piuttosto che come viscosità artificiale.

Risultati Numerici
L'articolo valida la teoria attraverso una serie di diagnostiche:

• Test Scalari: è stato dimostrato che le perturbazioni logaritmiche a media nulla creano un bias positivo dell'entropia (bias di Jensen) e che le ricostruzioni lineari, pur essendo a basso bias, non preservano la positività.
• Test di Amplificazione: ha confermato che i difetti di forza-stress nelle ricostruzioni logaritmiche grezze crescono linearmente con lo stiramento di fondo $\lambda_{\max}(A)$.
• Verifica del Difetto di Lavoro: ha mostrato che l'uso di tensori discordanti per stress ed entropia porta a difetti energetici a segno indefinito.
• Loop Dinamici ad Alto Numero di Weissenberg: nelle simulazioni con $Wi=80$ e alto stiramento, le ricostruzioni logaritmiche standard accumulavano errori di entropia e di forza-stress nel tempo, mentre la ricostruzione corretta controllava questi difetti, riducendo gli errori di forza di ordini di grandezza.
• Benchmark di Risoluzione: il metodo corretto ha mantenuto la compatibilità entropica e controllato i difetti di forza attraverso i raffinamenti della mesh e i diversi livelli di stiramento, mentre i metodi grezzi mostravano errori crescenti con l'aumento dello stiramento.

Significato e Rivendicazioni
L'articolo posiziona questo lavoro come colui che colma una lacuna critica tra la preservazione della positività (gestita dai metodi log-conformation/radice quadrata) e la compatibilità termodinamica (gestita dai metodi dissipativi dell'energia libera). Gli autori affermano che:

1. La positività è una condizione necessaria ma insufficiente per la robustezza ad alto numero di Weissenberg.
2. La correzione proposta è uno strato locale, computabile e controllato spettralmente che può essere aggiunto agli schemi esistenti ad alto ordine senza alterarne la struttura fondamentale.
3. Il metodo previene l'iniezione non fisica di energia garantendo che il tensore che entra nell'equazione del momento sia compatibile con le variabili di entropia.
4. Offre un vantaggio di risoluzione condizionale nei regimi di alto stiramento operando nello spazio delle variabili logaritmiche, che è regolare, pur rispettando il budget di entropia.

Gli autori concludono che, per i calcoli ad alto numero di Weissenberg, monitorare il valore minimo dell'autovalore non è sufficiente; è necessario anche tracciare l'eccesso di entropia della ricostruzione e i difetti di forza-stress per garantire la fedeltà fisica. La ricostruzione logaritmica corretta fornisce un meccanismo per imporre questa compatibilità rimanendo asintoticamente inattiva per flussi sufficientemente risolti.