Steady Incremental Viscosity Splitting Method for solving the stationary Navier-Stokes equation

Riepilogo Tecnico: Metodo di Splitting Incrementale di Viscosità Stazionario per le Equazioni di Navier-Stokes Stazionarie

Enunciato del Problema
Il lavoro affronta la soluzione numerica delle equazioni di Navier-Stokes stazionarie e incomprimibili con stabilizzazione grad-div. Le equazioni governative sono:
$$-\nu\Delta \mathbf{u} + \mathbf{u} \cdot \nabla \mathbf{u} + \nabla p - \gamma\nabla(\nabla \cdot \mathbf{u}) = \mathbf{f}, \quad \nabla \cdot \mathbf{u} = 0$$
in un dominio $\Omega$, con condizioni al contorno di Dirichlet omogenee. Gli autori notano che la risoluzione dei sistemi lineari a punto di sella non simmetrici risultanti (derivanti dalla linearizzazione standard di Picard) è computazionalmente costosa, in particolare per problemi 3D dove i solutori diretti falliscono a causa di vincoli di memoria e i solutori iterativi faticano con la struttura a punto di sella. Sebbene i metodi di splitting siano ben consolidati per flussi non stazionari, la loro applicazione a problemi genuinamente stazionari è meno sviluppata.

Metodologia
Gli autori propongono un metodo di Splitting Incrementale di Viscosità Stazionario (SIVS). Questo approccio adatta lo schema di Splitting Incrementale di Viscosità (IVS), precedentemente utilizzato per flussi non stazionari, al contesto stazionario. Il metodo opera come uno schema iterativo a punto fisso:

1. Passo 1 (Previsione della Velocità): Data l'iterata precedente $(\mathbf{u}^{k-1}, p^{k-1})$, risolvere un'equazione differenziale alle derivate parziali ellittica per una velocità intermedia $\tilde{\mathbf{u}}^k$:
   $$a(\tilde{\mathbf{u}}^k, \mathbf{v}) + c^*(\mathbf{u}^{k-1}, \tilde{\mathbf{u}}^k, \mathbf{v}) + b(p^{k-1}, \mathbf{v}) = \langle \mathbf{f}, \mathbf{v} \rangle$$
   Qui, la forma bilineare $a(\cdot, \cdot)$ include il termine viscoso e il termine di stabilizzazione grad-div ($\gamma$).

2. Passo 2 (Correzione della Pressione e Proiezione): Risolvere per la velocità aggiornata $\mathbf{u}^k$ e la pressione $p^k$ correggendo la soluzione intermedia:
   $$a(\mathbf{u}^k - \tilde{\mathbf{u}}^k, \mathbf{v}) + b(p^k - p^{k-1}, \mathbf{v}) = 0$$
   $$b(q, \mathbf{u}^k) = 0$$
   Questo passo garantisce che la velocità finale sia a divergenza nulla e soddisfi le condizioni al contorno.

Interpretazione Algebrica
Il lavoro interpreta il SIVS come un metodo di splitting algebrico. Eseguendo una decomposizione LU incompleta (ILU) sulla matrice del sistema, gli autori mostrano che il SIVS corrisponde a una specifica fattorizzazione a blocchi in cui la matrice del complemento di Schur $S = -B\tilde{A}^{-1}B^T$ (con $\tilde{A}$ che rappresenta l'operatore viscoso lineare) è spettralmente equivalente alla matrice di massa della pressione. Crucialmente, a differenza delle iterazioni standard di Picard che richiedono l'inversione di un diverso complemento di Schur ad ogni passo, lo schema SIVS richiede la risoluzione di un sistema con una matrice SPD fissa per la variabile pressione attraverso tutte le iterazioni non lineari. Ciò consente il pre-calcolo e il riutilizzo dei precondizionatori, riducendo significativamente il costo computazionale.

Contributi Chiave e Risultati Teorici

• Schema Innovativo: Introduzione del metodo SIVS per le equazioni di Navier-Stokes stazionarie con stabilizzazione grad-div, garantendo che il campo di velocità discreto sia a divergenza nulla e soddisfi le condizioni al contorno di Dirichlet complete.
• Analisi di Convergenza: Gli autori dimostrano la limitatezza uniforme delle iterazioni e la convergenza geometrica dello schema sotto una condizione di piccolezza sui dati ($\Lambda < 1/\sqrt{2}$).
• Efficienza Algebrica: Il metodo è dimostrato essere un approccio di splitting algebrico che disaccoppia le risoluzioni di velocità e pressione in sistemi lineari più semplici rispetto al sistema accoppiato standard a punto di sella.

Risultati Numerici
Gli autori validano il metodo utilizzando FreeFem++ con elementi Taylor-Hood ($P_2, P_1$) su tre casi di test:

1. Test di Convergenza: Utilizzando una soluzione prodotta artificialmente, il metodo ha raggiunto i tassi di convergenza ottimali attesi per velocità, pressione ed errori di divergenza. Il numero di iterazioni non lineari è rimasto uniforme attraverso i raffinamenti della mesh.
2. Cavità 2D Azionata dal Coperchio: Lo schema è stato testato a numeri di Reynolds fino a $Re=5000$. I risultati per i profili di velocità e le linee di flusso hanno corrisposto accuratamente ai dati di riferimento (Ghia et al.). I conteggi delle iterazioni erano comparabili al metodo Incrementale Picard-Yosida (IPY), ma gli autori prevedono un tempo CPU inferiore a causa dei sistemi lineari più semplici.
3. Flusso 3D Intorno a un Cilindro: Un problema di riferimento è stato risolto con $Re \approx 450$ (basato su $\nu=10^{-3}$). Lo schema è convergente in 11 iterazioni non lineari, producendo coefficienti di resistenza, portanza e caduta di pressione entro gli intervalli di riferimento forniti in letteratura.
4. Influenza dei Parametri: Test sul parametro grad-div $\gamma$ hanno indicato che valori più grandi di $\gamma$ accelerano la convergenza allo stato stazionario, sebbene possano deteriorare il condizionamento del sistema lineare nel Passo 1.

Significato
Il lavoro afferma che il metodo SIVS offre un'alternativa computazionalmente efficiente alle linearizzazioni standard di Picard per flussi incomprimibili stazionari. Sfruttando la stabilizzazione grad-div e lo splitting algebrico, il metodo produce sistemi lineari più semplici che sono più facili da precondizionare e risolvere iterativamente. Ciò rende l'approccio particolarmente adatto per problemi 3D su larga scala dove i solutori diretti sono non fattibili. Gli autori concludono che il metodo mantiene bassi conteggi di iterazioni riproducendo accuratamente i benchmark standard, suggerendo la sua fattibilità per applicazioni complesse di multiphysics.