Validare la mesh CFD: i criteri che fanno la differenza

Oggi parliamo di una tematica in cui, nella mia vita da analista, ho visto troppo spesso combinare pasticci: la validazione della mesh. Spesso si tende a pensare che, una volta cliccato su Generate Mesh e visto comparire un fitto reticolo di celle, il lavoro sia finito. Al massimo si chiama un collega più esperto che guarda la mesh e si affida all’occhiometria per testarne la validità. In realtà, è proprio lì che iniziano i guai (o i successi). Una mesh esteticamente gradevole non è necessariamente una mesh che produce risultati corretti.

In questo articolo esploreremo le metriche fondamentali per capire se la vostra discretizzazione è un capolavoro di ingegneria o una bomba a orologeria numerica, con un occhio di riguardo agli strumenti che Simcenter STAR-CCM+ ci mette a disposizione.

Prima di tuffarci nei numeri, facciamo un passo indietro. La Fluidodinamica Computazionale non risolve le equazioni di Navier-Stokes in modo continuo nello spazio (semplicemente, non si può), ma ha bisogno di discretizzarle, ossia di trasformare un dominio infinito in un numero finito di piccoli volumi di controllo, in ognuno dei quali vengono risolte le equazioni NS.

Il problema? Le equazioni sono estremamente schizzinose. Se una cella è troppo schiacciata, troppo storta o troppo diversa dalla sua vicina, l’errore numerico cresce, la convergenza rallenta e, nel peggiore dei casi, il residuo della pressione decide di decollare verso l’infinito (ed oltre).

Validare la mesh significa assicurarsi che la griglia che regge la nostra fisica non introduca errori artificiali tali da invalidare la soluzione. Non è solo una questione di precisione, è una questione di affidabilità dei risultati.

Le prime metriche da controllare sono puramente geometriche. Non ci dicono nulla sulla fisica, ma ci dicono molto sulla salute matematica dei nostri volumi. In STAR-CCM+, queste metriche sono calcolate automaticamente e accessibili tramite lo strumento “Check Mesh”.

La Face Validity misura se il centro di una cella si posiziona correttamente rispetto alle sue facce. In una cella di buona qualità, le normali delle facce puntano verso l’esterno, lontano dal centroide della cella. In una cella con validità delle facce scadente, una o più normali delle facce puntano verso l’interno, verso il centroide della cella.

Una validità delle facce pari a 1.0 significa che tutte le normali delle facce puntano correttamente lontano dal centroide della cella, mentre valori inferiori a 0.5 indicano celle non-convesse o peggio, celle che si auto-intersecano (con volume negativo).

Se avete celle con validità inferiore a 1.0, il solutore potrebbe faticare; se sono sotto lo 0.8, aspettatevi problemi seri di convergenza.

Star-CCM+ utilizza un indice di qualità ibrido che combina diversi fattori, come la deviazione dei vettori normali tra celle adiacenti.

Una Cell Quality bassa, vicina allo zero, indica solitamente celle molto distorte. Un buon analista punta a tenere la stragrande maggioranza delle celle sopra lo 0.1, idealmente molto più in alto.

In parole semplici, lo Skewness Angle (θ) misura quanto una cella è “storta” rispetto alla sua vicina. È il parametro che ci dice se l’informazione (pressione, velocità) può viaggiare correttamente da un centro cella all’altro senza “perdersi” o esplodere matematicamente.

Immagina due celle adiacenti. Per calcolare come variano le quantità tra loro, il software guarda due vettori: il vettore normale alla faccia che le separa e il vettore che unisce i due centri delle celle; l’angolo di skewness θ è l’angolo tra questi due vettori.  Significa che la mesh è perfettamente ortogonale. Più l’angolo cresce, più la mesh è distorta.

L’Aspect Ratio misura quanto una cella sia “allungata”. Immaginate un cubo perfetto: il suo rapporto è 1:1. Immaginate ora uno spaghetto: il rapporto è altissimo.

Nelle zone di core del flusso, vogliamo un Aspect Ratio vicino a 1. Tuttavia, negli strati limite, è normale e spesso desiderabile avere Aspect Ratio elevati (anche 10:1 o più) per catturare i gradienti normali alla parete senza sprecare milioni di celle nella direzione parallela al flusso.

Questa è forse la metrica più sottovalutata. Rappresenta il rapporto tra il volume di una cella e quello delle sue vicine.

I solutori odiano i salti bruschi: se passate da una cella millimetrica a una centimetrica istantaneamente, l’errore di troncamento nei gradienti esploderà.

Una volta che la geometria è “sana”, dobbiamo chiederci: la mesh è adatta alla fisica che vogliamo simulare? Qui le metriche non sono più solo geometriche, ma diventano soluzione-dipendenti.

Al Wall y+ dedicheremo non uno ma ben due articoli differenti di questa serie. Qui ricapitoliamo gli effetti chiari ma vi invito a leggerli appena li pubblicheremo.
Se state simulando turbolenza (e quasi certamente lo state facendo), il y+ è il vostro migliore amico.

O il vostro peggiore incubo.

È una coordinata adimensionale che misura quanto è fitta la mesh vicino alle pareti.

• High-y+ (Wall Functions): 30 < y+ < 100 (ideale per risparmiare risorse).
• Low-y+ (All y+ Wall Treatment): y+ ≈ 1 (necessario per catturare il sotto-strato viscoso).

Sbagliare il y+ significa sbagliare il calcolo della resistenza e del calore scambiato. STAR-CCM+ permette di visualizzare il y+ come una scena scalare sulla superficie dopo poche iterazioni: fatelo sempre!

Nessun articolo sulla validazione sarebbe completo senza menzionare lo studio di indipendenza. Il concetto è semplice: se raddoppio il numero di celle e il mio risultato (es. la portanza di un’ala) cambia del 10%, allora la mia mesh originale non era abbastanza densa.

Il gold standard è l’utilizzo del Grid Convergence Index (GCI), basato sull’estrapolazione di Richardson.

Dove Fs è un fattore di sicurezza, ε è l’errore relativo tra due mesh, r è il rapporto di refinement e ρ è l’ordine di accuratezza dello schema numerico.

Ci siamo quasi. Le metriche geometriche sono verdi, il y+ è centrato, ma prima di lasciare il solutore a macinare dati per tutta la notte (se vi va bene, perché rischiate che diverga dopo 5 iterazioni e ve ne accorgete solo al mattino dopo!), c’è un ultimo step fondamentale. Dobbiamo assicurarci che la mesh sia “topologicamente sana”. Saltare questi tre controlli finali significa esporsi a crash improvvisi o a una sottile ma devastante perdita di accuratezza.

• La caccia ai Volumi Negativi (Elementi Invertiti): È il controllo più banale ma più critico. Se durante la deformazione della mesh (magari per seguire un raggio di raccordo stretto o una feature millimetrica) il mesher genera elementi con volume zero o negativo, avete una bomba a orologeria nel dominio. Il solutore, non appena proverà a dividere per quel volume nelle equazioni discretizzate, si troverà davanti a una divisione per zero e interromperà il calcolo con un errore fatale. Elementi negativi = crash istantaneo.
• La trappola delle Interfacce Non Conformi: Quando collegate domini differenti — ad esempio la zona stazionaria di un vano motore con il dominio rotante di una ventola, o semplicemente zone con topologie di mesh diverse (es. da esaedrici a poliedrici) — STAR-CCM+ crea delle interfacce. Verificate sempre che non ci siano “buchi” o accoppiamenti parziali. Se le superfici non combaciano geometricamente o se la tolleranza di accoppiamento è errata, il flusso rimarrà letteralmente bloccato o, peggio, si creeranno ricircoli artificiali che invalideranno la fisica del problema.
• L’allineamento al flusso (Non-orthogonality numerica): Questa, lasciatemelo dire, è finezza da analisti esperti. Nelle zone dove il fluido viaggia ad alta velocità con una direzione preferenziale (come nel canale di un radiatore o lungo la carrozzeria), la direzione delle celle conta quanto la loro dimensione. Se le facce degli elementi sono fortemente disallineate rispetto alle linee di flusso reali, gli schemi di discretizzazione (anche quelli del secondo ordine o upwind avanzati) introdurranno una massiccia diffusione numerica. Orientare la mesh (usando ad esempio celle direzionali o estrusioni) in modo che le facce siano il più possibile parallele o perpendicolari ai vettori di velocità riduce l’errore artificiale, regalandovi una soluzione incredibilmente più nitida e vicina alla realtà.

Ingegneri e analisti passano ore a discutere sulle finezze dei modelli di turbolenza, sui setup dei solutori segregati o accoppiati e sulle condizioni al contorno più esotiche. Ma la verità è che potete impostare la fisica più raffinata del mondo, ma se la griglia che la sostiene è debole, state solo accelerando verso un risultato sbagliato.

Ricordate: una mesh pessima con un solutore del secondo ordine produrrà sempre risultati peggiori di una mesh eccellente con un solutore del primo ordine. Investite tempo nella griglia, e la fisica vi ringrazierà.

Validare la mesh non è un esercizio accademico per appesantire i report tecnici; è il primo paracadute che avete per assicurarvi che i millisecondi spesi dalla CPU (e i soldi delle licenze) stiano calcolando la realtà e non l’errore numerico del software.

Prima di lanciare il prossimo calcolo pesante e abbandonare la workstation per la notte, aprite il pannello di diagnostica di STAR-CCM+, analizzate i grafici di distribuzione della qualità e fatevi una domanda onesta: questa mesh merita davvero il tempo del mio solutore? Investite quel 10% di tempo in più sulla griglia oggi. La convergenza dei vostri residui, ma soprattutto la credibilità dei vostri risultati domani, vi ringrazierà.

Buon meshing.