Strumenti e criteri per tenere sotto controllo Wall y+

Se lavori nel mondo della fluidodinamica computazionale, sai che ci sono poche certezze nella vita: una di queste è che il y+ (ti ricordo che si legge “wuai plas”) non è mai quello che vorresti al primo colpo.

Il Wall y+ è quel parametro adimensionale che ci dice se la nostra mesh vicino alle pareti è giusta per il modello di turbolenza che abbiamo scelto; è il ponte tra la realtà fisica dello strato limite e l’astrazione numerica del software. Ma come si controlla davvero in Simcenter STAR-CCM+? Non basta lanciare il calcolo e sperare. Serve una strategia di valutazione rigorosa tra mappe colorate, istogrammi e soglie critiche.

Prima di tutto, un passo indietro per su alcuni concetti di base. Lo strato limite è dove succede tutta la magia (e il dolore): attrito, distacco della vena, scambio termico.

Per decidere come trattare la zona vicino alla parete, abbiamo visto che esistono in STAR-CCM+ due approcci:

• High y+: La prima cella si trova nella turbulent region: qui vogliamo un y+ > 30 (idealmente tra 30 e 300).
• Low y+: Siamo interessati ad una maggior definizione nel substrato viscoso: qui vogliamo un y+ ∼ 1.

In realtà esiste una terza via, un modello del software che permette di discernere in quale zona siamo tra le due e applicare le wall function corrette, chiamato All-y+, disponibile in STAR-CCM+. In questo caso, risulta fondamentale essere al di fuori della buffer zone con y+ compreso tra 10 e 30.

Sbagliare il y+ significa sbagliare il calcolo del drag, del calore o, peggio, ignorare un distacco del flusso che cambierebbe completamente i tuoi risultati. Fortunatamente, in STAR-CCM+, abbiamo tutti gli strumenti per non sbagliare, o almeno per capire dove abbiamo sbagliato prima di consegnare i risultati.

La prima cosa che facciamo tutti è creare una Scalar Scene. In questa fase, la Scalar Scene interpreta il ruolo della Vedetta. Il suo compito è quella di fornire una ricognizione aerea del problema, permettendoci di evidenziare a colpo d’occhio le macro-aree critiche. Vedere la mappa ti permette di capire lo stato generale della mesh vicino alla parete. Magari hai un ottimo y+ sul corpo principale, ma nei pressi di un raccordo o di un bordo d’uscita, la velocità accelera e il tuo y+ schizza alle stelle.

Se stai puntando a un approccio Low-Re (dove y+ < 1), imposta il massimo della scala a 5 o 10. Tutto ciò che è rosso “bruciato” è fuori target. Se invece usi le Wall Functions logaritmiche, imposta un range tra 30 e 300. Le zone blu (sotto 30) sono quelle dove il modello matematico sta soffrendo perché la mesh è troppo fitta per le assunzioni della funzione di parete.

Vediamo un caso pratico: immaginiamo di calcolare la resistenza di un sottomarino immerso, con approccio High y+; abbiamo dunque provato a dimensionare la nostra mesh con un target Wall y+ compreso tra 30 e 300 ovunque.

Questa prima immagine è una classica color map: molto smooth, con tanti livelli e colori sfumati, transizioni gradevoli… e non ci dice niente!

Questa seconda immagine invece è molto meno sfumata, con zone più nette: sicuramente meno gradevole, ma permette di evidenziare subito dove ci sono i problemi evidenziando le celle con y+ > 30.

Pro Tip: puoi sovrapporre un Data Focus alle scalar scene per mettere in evidenza zone in determinati range di field function mentre visualizzi grandezze differenti. È molto utile per capire ad esempio la velocità nelle zone di basso y+, per andare ancora più nel dettaglio della fisica dietro questi fenomeni.

Le mappe sono belle, ma i revisori tecnici vogliono i numeri. Quanta percentuale della tua superficie è fuori target? Per rispondere, dobbiamo passare a strumenti più tecnici come gli Histogram Plot.

Per una analisi asettica che rappresenti in maniera oggettiva lo stato della mia mesh abbiamo bisogno dell’Intelligence: l’Histogram Plot.

L’Intelligence non si emoziona. Prende tutti i dati di tutte le facce della parete e li trasforma in un rapporto statistico freddo e oggettivo. Se la Vedetta ci dice dove sono i problemi, l’Istogramma ci dice quanto sono gravi.

L’istogramma ti dà una visione d’insieme della qualità della mesh. Se vedi un picco stretto intorno al valore di target, hai fatto un lavoro da manuale. Se vedi una distribuzione spalmata che va da 0.5 a 1000, la tua mesh non è uniforme e i risultati potrebbero essere instabili.

Pro Tip: Ricorda di pesare l’istogramma per Area. In questo modo, le facce più grandi contano di più nell’analisi statistica rispetto a piccoli frammenti di mesh negli angoli. A sinistra è mostrato un istogramma pesato per area e a destra uno per frequenza: è facile notare come la percezione di quello di destra sia di un Wall y+ più spalmato di quanto non sia in realtà.

A volte la mappa è confusa e l’istogramma ti dice che c’è qualcosa che non va, ma non capisci dove. Qui entra in gioco il nostro cecchino (il threshold), uno dei tool più sottovalutati di STAR-CCM+; il suo compito è isolare ed eliminare (o meglio, identificare per la correzione) solo le celle che si comportano come ribelli.

Vedrai illuminarsi come per magia i pezzi di mesh cattivi. È un modo chirurgico per identificare dove devi intervenire con un Volumetric Control o dove devi ridurre il Base Size della mesh superficiale.

Questo metodo è fondamentale per i report di validazione: “Il 95% della superficie ha un y+ > 30, le zone in cui la prima cella ricade nella buffer zone (visibili nel Threshold) sono localizzate solo in corrispondenza dei punti di ristagno e non influenzano il drag globale”.

Valutare il y+ non è un’attività da fare alla fine, ma un processo iterativo. Spesso richiede di tornare indietro al Mesher, regolare il Prism Layer Thickness o il numero di strati, e ricalcolare.

In STAR-CCM+, grazie alla pipeline automatizzata, puoi cambiare un parametro nel mesher e vedere l’istogramma aggiornarsi in tempo reale dopo qualche iterazione di calcolo. Non accontentarti di una media globale; usa mappe, istogrammi e soglie per garantire che la tua fisica sia solida quanto la tua geometria.

La buona notizia però è che per geometrie non eccessivamente complicate, e con la giusta dose di esperienza e conoscenza, ci sono tanti metodi di calcolo per stimare lo spessore di prima cella a priori, garantendo di ridurre al minimo (se sei fortunato, a zero) il numero di iterazioni necessarie per raggiungere un buon dimensionamento di strato limite.