Pillole Di CFD #2 - Modelli di turbolenza - Come orientarsi per la scelta corretta?
Aggiornamento: 13 giu
I flussi di un fluido possono essere categorizzati in due classi:
Flussi laminari: quando non si ha rimescolamento e dissipazione di energia e velocità;
Flussi turbolenti: quando si è in presenza di una certa oscillazione di una grandezza che può portare ad un rimescolamento ed a una dissipazione.
Questo aspetto è governato da una grandezza molto importante, ossia il famoso numero di Reynolds.
Figura 1: numero di Reynolds
Questo numero, che è un parametro adimensionale, rappresenta il rapporto tra le forze di inerzia del fluido e quelle viscose.
Se le forze di inerzia sono superiori a quelle viscose, allora si ha un elevato rimescolamento, mentre invece se a vincere sono le forze viscose, allora prevalgono gli effetti dissipatori.
I parametri che mette in relazione sono:
la densità rho;
la velocità del fluido V;
la dimensione caratteristica L;
la viscosità dinamica ni.
A seconda del suo valore e della tipologia di flusso, possiamo trovarci in una situazione di flusso laminare oppure di flusso turbolento.
In tabella sono rappresentati i possibili casi con i relativi valori di riferimento:
Tabella 1: Valori dl numero di Reynolds per capire se il flusso è laminare o turbolento
Come vediamo, a seconda che il flusso avvenga internamente o esternamente, si hanno diversi valori di Reynolds.
Dobbiamo solo capire cosa si intende per turbolenza.
COSA SI INTENDE PER TURBOLENZA?
Molto spesso sentiamo dire questa frase :" la turbolenza è la presenza di vorticità all'interno di un fluido".
Questa definizione, detta in questo modo. è fortemente sbagliata. Vortici e turbolenza sono due aspetti collegati ma staccati.
La turbolenza, nella realtà, altro non è che la fluttuazione di una certa grandezza all'interno del fluido.
Figura 2: effetti della turbolenza nell'andamento della velocità
Nel grafico che vediamo sopra, il profilo di velocità è dato dalla somma della velocità media Vm e della fluttuazione della velocità V'.
Questa fluttuazione genera un rimescolamento delle particelle fluide che, rimescolandosi, possono dar luogo ad una vorticosità che tende a smorzare la grandezza di interesse.
Quindi possiamo dire che la turbolenza genera dei vortici ma non il contrario, ossia che i vortici generano turbolenza, anche perchè si possono avere situazioni di flusso stazionario con situazioni di vortice.
QUALI SONO GLI APPROCCI COMPUTAZIONALI:
Nel corso del tempo sono stati sviluppati diversi approcci al modo di calcolare la turbolenza. Tra questi i più importanti sono:
Direct Numerical Simulation (DNS);
Large Eddy Simulation (LES)
Reynolds averaged Navier-Stokes simulation (RANS)
Ognuno di questi approcci ha i suoi pro e contro.
DIRECT NUMERICAL SIMULATION
L'approccio DNS implementa le equazioni della turbolenza direttamente nelle equazioni di Navier-Stokes. Non vi è nessuna idealizzazione. Per tale motivo questo metodo è solitamente relegato al campo della ricerca e poche volte viene utilizzato nelle applicazioni industriali.
LARGE EDDY SIMULATION
L'approccio LES risolve direttamente i vortici più grandi mentre modella quelli più piccoli.
Come metodo risulta più leggero del DNS ma molto spesso ancora troppo costoso per applicazioni pratiche
REYNOLDS AVERAGED NAVIER STOKES
Simile al LES, ma la risoluzione del campo si ferma molto prima e si da più spazio alla parte modellistica.
Nell'immagine qui sotto è molto più chiaro quello che succede. Alla base vi è l'ipotesi che i vortici (le fluttuazioni), da che sono grandi, a cascata, diventino piccoli.
A seconda dell'approccio che si utilizza, viene scelto a che punto effettuare una modellazione matematica del vortice invece che ricercare la sua forma esatta nel campo.
L'approccio DNS risolve tutto il campo nella sua interezza, mentre invece il LES ed il RANS arrivano fino ad un certo punto e poi approssimano.
Figura 3: confronto tra le diverse mentalità di approccio al problema della turbolenza.
Mentre per la parte di campo che viene risolta non vi è nulla da dire (sono le equazioni di Navier-Stokes ad essere risolte), dobbiamo comprendere come viene affrontata la modellazione delle fluttuazioni.
Poiché a livello industriale il RANS è il più utilizzato, in questo articolo parleremo solamente dei possibili modelli di turbolenza che vengono accoppiati con il RANS.
MODELLI DI TURBOLENZA
Per il RANS il problema viene chiuso modellando il tensore degli sforzi di Reynolds tramite un approccio ad una o due equazioni. Si ottengono quindi:
Modelli Reynolds-Stress Models (RSM): In questo caso le componenti del tensore sono risolte tramite equazioni del trasporto. Offre notevoli vantaggi in flussi 3D con linee di flusso che si avvitano su loro stesse. I modelli sono complessi e il costo computazionale elevato.
Modelli Eddy Viscosity: si utilizzano due equazioni basate sulla turbolenza e sulla viscosità. Costo computazionale basso. Utile per quei casi in cui gli sforzi a parete non sono elevati (non si ha distacco della vena fluidica).
Tra i modelli Eddy Viscosity, i più utilizzati (che trovate maggiormente nei principali software) sono:
Modello k-eps: modello Eddy Viscosity a due equazioni, ossia l'energia cinetica turbolenta e la dissipazione. La viscosità turbolenta è modellata come prodotto della velocità di turbolenza e della grandezza della turbolenza. E' un ottimo compromesso tra il costo computazionale e l'accuratezza. ATTENZIONE: se siamo in presenza di flussi che ruotano, distacco di vena fluidica o situazioni simili, questo metodo sbaglia completamente i risultati.
Modello k-omega: modello Eddy Viscosity a due equazioni, ossia l'energia cinetica turbolenta e frequenza. La viscosità turbolenta è modellata come prodotto della velocità di turbolenza e della grandezza della turbolenza. Perfetto per molti problemi ingegneristici, riesce a modellare meglio i flussi nelle condizioni al contorno, in quanto utilizza delle funzioni di parete (questo argomento sarà trattato in un altro articolo).
Modello Shear Stress Transport (SST): E' un modello via di mezzo tra la facilità di convergenza del modello k-eps e l'accuratezza del modello k-omega. Dal primo prende la capacità di riuscire a predire bene il comportamento all'intenro del fluido mentre dal secondo la capacità di approssimare il comportamento a parete. E' difficile da far convergere.
In immagine è spiegato come il modello SST ragiona a livello di zone.
Figura 4: schematizzazione del modello SST, con approccio k-omega alla aprete e k-eps nel fluido
Vediamo adesso i diversi risultati per il modello K-eps e per il modello K-omega. In figura è rappresentata l'energia cinetica di turbolenza. Si ricorda che il flusso, in entrambi i casi, presenta le stesse condizioni al contorno.
Figura 5: energia cinetica turbolenta nel caso di modello k-eps
Figura 6: energia cinetica turbolenta nel caso di modello k-omega
Quello che si nota è come il modello k-eps, per quanto sia in grado di convergere molto velocemente, non è in grado di simulare al meglio il distacco del fluido, cosa che il modello SST riesce a fare perfettamente (anche se con problematiche nella convergenza).
Si può quindi dire che l'approccio più conveniente è il metodo RANS con modello k-eps o k-omega a seconda delle applicazioni (che variano a seconda se è o meno presente un possibile distacco di fluido).
Ovviamente negli esempi qui sopra, la mesh andrebbe adattata al meglio per il software e per riuscire a coprire al meglio le condizioni di parete. questo è dovuto alla capacità del modello di riuscire ad approssimare lo strato limite che si genera intorno ad un corpo.
Come si modella questo strato? Ci sono tecniche particolari?
Queste domande saranno affrontate in altri articoli delle #PillolediCFD
Ing. Francesco Grispo
---------------------------------------------------------------------------------------------------------------------------
FGCAEANALYST è al servizio delle aziende. Ti seguiamo passo passo nell'analisi dei tuoi componenti mettendoci la stessa attenzione che metteremmo se fosse nostro.
Contattaci per una consulenza gratuita cliccando al seguente link.
Se invece vuoi imparare ad applicare correttamente i modelli di turbolenza, contattami per dei corsi di CFD al seguente link.