A HHO formulation for variable density incompressible flows where the density is purely advected

Questo lavoro propone una formulazione Hybrid High-Order (HHO) per flussi incomprimibili a densità variabile che garantisce la conservazione esatta del volume e l'avvezione pura della densità, sfruttando discretizzazioni spaziali ibride e schemi temporali ESDIRK per simulare con precisione e robustezza miscele di fluidi immiscibili, come dimostrato attraverso analisi di convergenza e lo studio dell'instabilità di Rayleigh-Taylor.

L'essenziale

Immagina di dover mescolare due liquidi diversi, come olio e acqua, in un contenitore. Se li agiti, si muovono, si mescolano e cambiano forma, ma non si fondono mai davvero: rimangono due cose distinte che scivolano l'una sull'altra. Simulare questo comportamento al computer è una sfida enorme per gli scienziati, perché i liquidi sono complessi e i computer tendono a "confondersi" quando le cose si muovono velocemente o quando le densità cambiano drasticamente.

In questo articolo, Lorenzo Botti e Francesco Carlo Massa presentano un nuovo "modo di fare i calcoli" (chiamato HHO) per risolvere questi problemi. Ecco una spiegazione semplice, usando analogie quotidiane.

1. Il Problema: Il "Gioco del Telefono Senza Fili"

Immagina di voler tracciare il movimento di una goccia d'inchiostro in un fiume. Se il tuo metodo di calcolo è impreciso, la goccia potrebbe sembrare che si allarghi da sola, che sparisca o che si muova dove non dovrebbe. Nel mondo reale, se l'acqua è incomprimibile (non si può schiacciare), il volume deve rimanere esattamente lo stesso. Se il computer sbaglia anche di poco, l'acqua potrebbe "sparire" o "crearsi dal nulla", rendendo il risultato sbagliato.

2. La Soluzione: Il Metodo "Ibrido" (HHO)

Gli autori hanno creato un metodo chiamato HHO (Hybrid High-Order). Pensa a questo metodo come a un team di ispettori che controllano un edificio stanza per stanza.

• Ispettori interni (Cellule): Ogni stanza (o "cella" della griglia del computer) ha i suoi ispettori che controllano cosa succede al centro della stanza.
• Ispettori di confine (Scheletro): Ma c'è di più! Ci sono anche ispettori che stanno esattamente sulle porte e sui muri che separano le stanze.

Questa è la parte "Ibrida". Invece di guardare solo il centro della stanza, il metodo controlla anche i confini. Questo permette di scambiare informazioni tra le stanze in modo molto preciso, come se gli ispettori si passassero dei bigliettini sulle porte per assicurarsi che nessuno stia rubando o creando acqua dal nulla.

3. La Magia: "Pure Advection" (Il Treno dei Densità)

Uno dei punti di forza di questo metodo è che tratta la densità (quanto è "pesante" il fluido in un punto) come un passeggero su un treno.

• Se il treno (il fluido) si muove, il passeggero (la densità) si sposta con lui.
• Il passeggero non cambia peso, non sparisce e non ne nasce uno nuovo dal nulla.
• Il metodo garantisce che il volume totale dell'acqua rimanga esattamente lo stesso, fino all'ultima virgola del computer. È come se avessi un conto in banca che non può avere errori di arrotondamento: se hai 100 euro, ne avrai sempre 100, anche dopo mille transazioni.

4. La Robustezza: Il "Paracadute" contro gli Errori

Quando i fluidi si muovono molto velocemente (come in un uragano o in un'esplosione), i calcoli tendono a diventare instabili.

• Pressure-Robustness (Robustezza alla pressione): Immagina di spingere un carrello. Se spingi in modo sbagliato, il carrello potrebbe oscillare. Questo metodo è così intelligente che, anche se la "pressione" (la forza che spinge) è calcolata con un po' di imprecisione, il movimento del fluido (la velocità) rimane perfetto. È come avere un carrello con le ruote che si auto-aggiustano: non importa quanto spingi male, il carrello va dritto.
• Stabilità: Il metodo usa delle "stabilizzazioni" (come dei paracadute) per evitare che i numeri impazziscano quando il fluido scorre veloce.

5. Il Tempo: Il "Cronometro Intelligente" (ESDIRK)

Per simulare come le cose cambiano nel tempo, usano un metodo chiamato ESDIRK.
Immagina di dover attraversare un fiume a passi. Un metodo normale farebbe un passo grande e spera di non cadere. Questo metodo, invece, fa molti piccoli passi intermedi prima di decidere dove atterrare.

• Controlla il terreno prima di fare il passo vero.
• Se il terreno è scivoloso, si adatta immediatamente.
• Questo permette di usare passi temporali più grandi senza perdere precisione, risparmiando tempo di calcolo.

6. Risparmiare Memoria: La "Condensazione Statica"

I computer hanno una memoria limitata. Calcolare tutto per ogni punto della griglia richiederebbe troppa memoria.

• L'analogia: Immagina di dover risolvere un puzzle gigante. Invece di tenere in mano tutti i pezzi, questo metodo risolve prima i pezzi interni di ogni sezione e poi butta via i pezzi interni, tenendo solo i pezzi che toccano i bordi.
• In questo modo, il puzzle finale da risolvere è molto più piccolo, ma il risultato è lo stesso. Questo permette di fare calcoli molto complessi anche su computer non superpotenti.

7. Il Test Finale: L'Instabilità di Rayleigh-Taylor

Per dimostrare che funziona, hanno simulato un classico esperimento: mettere un fluido pesante sopra uno leggero (come olio sopra acqua, ma al contrario).

• Cosa succede: Il fluido pesante affonda creando delle "piume" o spirali intricate.
• Il risultato: Il loro metodo riesce a vedere queste spirali con una precisione incredibile, anche quando i fluidi sono molto diversi tra loro (alta densità) o quando si muovono molto velocemente.
• Hanno mostrato che usando un metodo "ad alta precisione" (polinomi di grado alto) su una griglia grossa, ottengono lo stesso risultato di un metodo "a bassa precisione" su una griglia piccolissima, ma molto più velocemente. È come se avessero trovato una scorciatoia magica per ottenere risultati da supercomputer usando risorse più modeste.

In Sintesi

Questo articolo presenta un nuovo modo di simulare fluidi che:

1. Non sbaglia mai il volume (l'acqua non sparisce).
2. È molto preciso anche con griglie non perfette.
3. È veloce perché usa trucchi matematici per risparmiare memoria.
4. Resiste agli errori quando i fluidi si muovono velocemente.

È come passare da una mappa disegnata a mano, piena di errori, a un GPS satellitare di ultima generazione che ti dice esattamente dove sono le auto, anche nel traffico più caotico.

Sommario tecnico

Ecco un riassunto tecnico dettagliato del documento in lingua italiana, strutturato secondo le sezioni richieste.

Titolo: Formulazione HHO per flussi incomprimibili a densità variabile con densità puramente advettiva

Autori: Lorenzo Botti e Francesco Carlo Massa (Università di Bergamo, Italia)

---

1. Il Problema

Il lavoro si concentra sulla simulazione numerica di flussi di fluidi incomprimibili a densità variabile, tipicamente associati a miscele di fluidi immiscibili. Le sfide principali affrontate includono:

• Conservazione esatta del volume: È fondamentale garantire che la densità venga trasportata puramente per avvezione (senza diffusione numerica spuria) per mantenere la massa e il volume corretti.
• Robustezza: La necessità di gestire regimi dominati dalla convezione e forze di corpo irrotazionali senza che gli errori nella pressione influenzino la velocità (proprietà di pressure-robustness).
• Stabilità: La difficoltà di mantenere la stabilità e la positività della densità (evitando valori negativi) in presenza di alti rapporti di densità (numeri di Atwood elevati) e forti gradienti interfacciali.
• Efficienza computazionale: La gestione di problemi ad alta accuratezza spaziale e temporale (ordine elevato) con un footprint di memoria ridotto.

2. Metodologia

Gli autori propongono una formulazione ibrida ad alto ordine (HHO - Hybrid High-Order) combinata con schemi di integrazione temporale ESDIRK (Explicit Singly Diagonal Implicit Runge-Kutta).

• Discretizzazione Spaziale (HHO):

  • Utilizza spazi ibridi per velocità, densità e pressione, definiti sia sugli elementi della mesh che sulle loro facce (scheletro).
  • Costruzione di un campo di velocità H(div): Grazie alla scelta di spazi polinomiali locali e alla pressione ibrida, viene ricostruita un'approssimazione della velocità conforme a $H(\text{div})$ e a divergenza nulla punto per punto.
  • Trattamento dell'avvezione: I termini convettivi per la quantità di moto e la massa sono trattati con formulazioni skew-simmetriche (per la stabilità) e stabilizzazioni upwind.
  • Condizioni al contorno: Imposizione debole delle condizioni al contorno di Dirichlet e Neumann.
  • Proprietà chiave: La conservazione del volume è imposta cella per cella fino alla precisione della macchina, garantendo che l'equazione di conservazione della massa sia una derivata temporale materiale della densità.

• Discretizzazione Temporale (ESDIRK):

  • Utilizzo di schemi ESDIRK di ordine elevato, specificamente concepiti per Equazioni Differenziali-Algebriche (DAE) di indice 2, necessarie per gestire le condizioni al contorno non passive e la natura del problema.
  • Il metodo è implicito e completamente ibrido.

• Algebra e Soluzione:

  • Condensazione Statica: Le incognite interne agli elementi vengono eliminate localmente, lasciando accoppiate globalmente solo le incognite sulle facce. Questo riduce drasticamente la dimensione e la densità della matrice Jacobiana globale.
  • Precondizionamento: Per la risoluzione dei sistemi lineari all'interno del metodo di Newton, viene utilizzata una strategia di precondizionamento p-multigrid, che sfrutta la struttura degli spazi polinomiali per migliorare l'efficienza degli solver iterativi.

3. Contributi Chiave

1. Formulazione Ibrida per Densità Variabile: Estensione dei metodi HHO (precedentemente limitati a densità costante) a flussi a densità variabile, garantendo la purezza dell'avvezione della densità.
2. Pressure-Robustness: La formulazione garantisce che le forze di corpo irrotazionali influenzino solo il campo di pressione, rendendo gli errori di velocità e densità indipendenti dall'accuratezza della rappresentazione della pressione.
3. Stabilità e Conservazione:
   • Dimostrazione teorica e numerica della stabilità dissipativa della formulazione.
   • Conservazione esatta del volume e della massa a livello discreto.
   • Capacità di preservare i limiti della densità (bound-preserving) a ordini bassi ($k=0$), essenziale per alti rapporti di densità.
4. Efficienza: La combinazione di condensazione statica e precondizionamento p-multigrid permette di ottenere soluzioni ad alta accuratezza con un costo computazionale e un uso di memoria significativamente ridotti rispetto alle formulazioni tradizionali.

4. Risultati Numerici

I risultati sono validati attraverso diversi casi di studio:

• Test di Kovasznay (Flusso Stazionario):

  • Dimostrazione della purezza dell'avvezione della densità (la densità rimane costante e uniforme).
  • Convergenza spaziale ottimale: ordine $k+2$ per la velocità e $k+1$ per la pressione e il gradiente di velocità in norma $L^2$.
  • Confronto tra le varianti HHO-$\pi^k$ e HHO-$\pi^{k+1}$: la variante $\pi^{k+1}$ mostra maggiore accuratezza per $k=0$, mentre per $k>0$ entrambe raggiungono i tassi teorici.

• Soluzione Manufactured di Guermond e Quartapelle (Flusso Transitorio):

  • Validazione della convergenza temporale con schemi ESDIRK a 3 e 6 stadi.
  • Conferma della pressure-robustness: l'errore di velocità è indipendente dall'errore di pressione.
  • Ottenimento di ordini di convergenza temporali elevati (fino al 4° ordine per la densità con lo schema a 6 stadi).

• Instabilità di Rayleigh-Taylor (RTI):

  • Basso numero di Atwood ($At=0.5$): Simulazioni a $Re=1000$ e $Re=5000$. Confronto tra mesh grossolane con ordine alto ($k=6$) e mesh fini con ordine basso ($k=1$). I risultati mostrano un ottimo accordo, evidenziando come l'ordine alto su mesh grossolane possa essere più efficiente (riduzione di 14 volte della dimensione della matrice Jacobiana) rispetto all'ordine basso su mesh fini.
  • Alto numero di Atwood ($At=0.75$): Simulazione con rapporto di densità 7:1. L'uso di $k=0$ garantisce la preservazione dei limiti della densità (nessuna densità negativa), permettendo di gestire l'instabilità senza instabilità numeriche, sebbene con una risoluzione dell'interfaccia inferiore rispetto agli ordini superiori.

5. Significato e Conclusioni

Questo lavoro rappresenta un passo significativo nello sviluppo di strumenti computazionali ad alta accuratezza per la fluidodinamica multifase.

• Impatto Scientifico: Dimostra che i metodi ibridi (HHO) possono essere estesi con successo a problemi complessi di flussi a densità variabile, offrendo proprietà di conservazione e robustezza superiori rispetto ad altri metodi (come DG o FV tradizionali).
• Impatto Pratico: La capacità di simulare flussi con alti rapporti di densità mantenendo la stabilità e l'efficienza (grazie alla condensazione statica e al p-multigrid) rende questa formulazione promettente per applicazioni industriali e di ricerca avanzata, come la dinamica dei fluidi multifase.
• Prospettive Future: Gli autori indicano come passo successivo l'estensione di questa formulazione al sistema Cahn-Hilliard-Navier-Stokes (CHNS) per modellare flussi multifase più fisicamente rilevanti.

In sintesi, la formulazione proposta combina l'accuratezza spaziale e temporale di alto ordine con la robustezza necessaria per i flussi reali, offrendo un compromesso ottimale tra precisione fisica ed efficienza computazionale.