Ensemble algorithm for parametrized flow problems with energy stable open boundary conditions

Sintesi Tecnica: Algoritmo d'Insieme per Problemi di Flusso Parametrizzati con Condizioni al Contorno Aperte Energeticamente Stabili

Definizione del Problema
Il documento affronta la sfida computazionale di risolvere insiemi (ensemble) di equazioni di Navier-Stokes incomprimibili soggetti a variazioni delle condizioni iniziali, delle forze di corpo e dei coefficienti di viscosità ($\nu_j$). In applicazioni quali l'analisi di sensibilità, la modellazione a ordine ridotto e la previsione d'insieme, il numero di membri dell'insieme ($J$) può essere elevato. I metodi semi-impliciti standard richiedono l'assemblaggio e la risoluzione di un sistema lineare distinto per ogni membro $j$ ad ogni passo temporale, a causa della dipendenza della matrice di sistema da $\nu_j$ e dai termini non lineari. Ciò comporta costi computazionali proibitivi. Sebbene gli approcci completamente espliciti permettano un singolo assemblaggio del sistema, essi impongono condizioni sul passo temporale molto restrittive, specialmente per flussi con numeri di Reynolds elevati. I precedenti metodi d'insieme per viscosità costante o variazioni specifiche della viscosità spesso richiedevano assunzioni restrittive sui coefficienti di viscosità (ad esempio, differenze relative limitate) per garantire la stabilità.

Metodologia
Gli autori propongono un nuovo metodo di calcolo d'insieme che disaccoppia la matrice del sistema lineare dall'indice del membro dell'insieme $j$, consentendo la soluzione di un singolo sistema lineare con molteplici termini noti ad ogni passo temporale. Il nucleo della metodologia prevede:

1. Strategia di Scissione: Il termine convettivo non lineare viene scisso in una componente media dell'insieme e una componente fluttuante. La componente media è trattata semi-implicitamente, mentre la componente fluttuante è trattata esplicitamente.
2. Trattamento della Viscosità: A differenza di lavori precedenti (ad esempio, [8]) che richiedevano vincoli specifici sulla variazione dei coefficienti di viscosità per garantire la stabilità, questo articolo introduce un trattamento differente del termine diffusivo. Definendo una viscosità massima $\nu_\infty = \max_j \nu_j$ e scindendo il termine di diffusione, gli autori derivano uno schema in cui la matrice del sistema lineare è indipendente da $j$ senza imporre restrizioni sull'entità delle variazioni di viscosità.
3. Condizioni al Contorno Aperte (OBC): Il metodo è esteso a domini con confini aperti ($\Gamma_N$). Il confine è partizionato in regioni di outflow (uscita) e backflow (riflusso). Gli autori utilizzano una condizione al contorno energeticamente stabile che coinvolge una scala di lunghezza caratteristica $L$ e una funzione di tipo Heaviside per gestire la natura convettiva dell'outflow.
4. Discretizzazione: Il documento presenta sia schemi di avanzamento temporale del primo ordine (Algoritmi 1–3) che del secondo ordine (Algoritmi 4–6). La discretizzazione spaziale utilizza spazi di elementi finiti conformi (specificamente coppie Taylor-Hood $P_2/P_1$) con forme trilinari skew-simmetrizzate per garantire la stabilità energetica.
5. Analisi della Stabilità: Prove teoriche stabiliscono la stabilità energetica per gli schemi del primo ordine sotto specifiche condizioni CFL che coinvolgono le fluttuazioni di velocità. Per gli schemi del secondo ordine, sebbene non venga fornita una prova rigorosa di stabilità senza restrizioni sulla viscosità (a causa della complessità del termine $\nabla \partial_{tt} u$), esperimenti numerici suggeriscono la stabilità sotto le stesse condizioni degli schemi del primo ordine.

Contributi Chiave

• Stabilità della Viscosità Senza Restrizioni: Il principale contributo teorico è uno schema d'insieme per coefficienti di viscosità variabili che garantisce la stabilità senza l'assunzione restrittiva $|\nu_j - \nu|/\nu \leq \sqrt{\mu}$ presente nella letteratura precedente.
• Estensione alle Condizioni al Contorno Aperte: Gli autori estendono i metodi d'insieme a problemi con confini aperti, derivando formulazioni deboli e provando la stabilità energetica per schemi del primo ordine sotto specifici vincoli del passo temporale.
• Efficienza Algoritmica: Gli algoritmi proposti risultano in un singolo matrice dei coefficienti per ogni passo temporale, risolvibile efficientemente mediante solver per sistemi con molteplici termini noti, riducendo significativamente il costo computazionale rispetto alla risoluzione di $J$ sistemi indipendenti.

Risultati
Esperimenti numerici sono stati condotti utilizzando il pacchetto FreeFem++:

1. Studio di Convergenza: Un test di convergenza su un dominio con soluzioni esatte ha confermato i tassi di convergenza del secondo ordine attesi sia per la velocità che per la pressione per i proposti schemi del secondo ordine. L'accuratezza è stata comparabile a quella delle simulazioni sequenziali indipendenti.
2. Flusso Attorno a un Cilindro: Gli algoritmi sono stati testati su un problema benchmark riguardante il flusso attorno a un cilindro con tre diversi valori di viscosità. I risultati (drag, portanza e caduta di pressione) hanno concordato con i valori di riferimento della letteratura che utilizzano condizioni di Dirichlet e di "do-nothing" in uscita. La condizione al contorno aperta ha permesso ai vortici di uscire naturalmente dal dominio, evitando i problemi di riflessione artificiale associati ai profili di outflow parabolici prescritti.
3. Flusso in un Canale con Contrazione: È stato simulato uno scenario di flusso complesso con due uscite e viscosità variabili. I risultati dell'insieme hanno concordato qualitativamente con le simulazioni DNS indipendenti, nonostante le simulazioni d'insieme siano state eseguite su una mesh più grossolana stabilizzata da un filtro non lineare adattivo.

Significatività e Rivendicazioni
L'articolo sostiene di offrire un quadro robusto ed efficiente per le simulazioni di flussi parametrizzati. Rimuovendo le assunzioni restrittive sulle variazioni di viscosità richieste dai precedenti metodi d'insieme, l'algoritmo proposto amplia l'applicabilità delle simulazioni d'insieme a problemi con parametri fisici altamente variabili. L'estensione alle condizioni al contorno aperte con stabilità energetica dimostrabile è un passo significativo verso simulazioni realistiche di flussi non limitati. Gli autori concludono che, sebbene gli schemi d'insieme funzionino bene qualitativamente e quantitativamente rispetto alle simulazioni indipendenti, sono necessarie ulteriori ricerche per quantificare pienamente i vantaggi computazionali, in particolare riguardo ai confronti del tempo di CPU su mesh adattivamente raffinate per flussi con numeri di Reynolds elevati.