header blog

Blog

Storie di innovazione e facili tutorial per conoscere il mondo dell'opensource

Tutorial: Dimensionamento della Mesh per Simulazioni CFD in OpenFOAM

Introduzione

Nelle simulazioni di fluidodinamica computazionale (CFD), il dimensionamento della mesh è cruciale per ottenere convergenza e risultati affidabili. Le mesh per CFD richiedono una risoluzione adeguata principalmente vicino ai muri (wall-condition) per catturare accuratamente gli effetti dello strato limite. In OpenFOAM, come in molti altri software, uno degli approcci per dimensionare la mesh vicino alle pareti solide si basa sul concetto di Y+. Questo tutorial spiegherà come utilizzare Y+ per dimensionare correttamente la mesh in una simulazione CFD bilanciando tempi di calcolo e affidabilità.

 

resolution mesh orange concept
Forse abbiamo esagerato con l'approssimazione della geometria reale?

Cos'è Y+?

Y+ è un parametro adimensionale che rappresenta la distanza dalla parete alla prima cella della mesh, in unità di lunghezze di scala dello strato limite viscoso. Prima di generare la mesh è quindi necessario un'idea preliminare della dimensione delle celle vicino alle pareti basata sull'Y+ che vogliamo ottenere per le condizioni al contorno che andremo a impostare. Valori tipici di Y+ per diversi modelli di turbolenza sono:

  • Per modelli con trattamento integrale dello strato limite (ad es. k-omega SST): Y+ ~ 1.
  • Per modelli che usano funzioni di parete (ad es. k-epsilon standard): Y+ compreso tra 30 e 300.

Y+ è definito come:

y+=uτyν

dove:

  • uτ = velocità di attrito (m/s), calcolabile da uτ=τwρ
  • τw = sforzo di parete (Pa)
  • y = distanza normale dalla parete alla prima cella della mesh (m)
  • ν = viscosità cinematica del fluido (m²/s), ν=μρ

Per calcolare y+, è spesso necessario stimare uτ che dipende dallo sforzo di parete τw.

Esempio: Se hai un flusso d'aria (densità ρ=1.225 kg/m³, viscosità μ=1.81×105 Pa·s) con uno sforzo di parete τw=0.5 Pa  e vuoi un valore di y+=30, la distanza y dalla parete alla prima cella della mesh sarà:

uτ=0.51.2250.64 m/s

ν=1.81×1051.2251.48×105 m²/s

y=30×1.48×1050.646.94×104 m

Quindi, la distanza dalla parete alla prima cella della mesh dovrebbe essere di circa 0.694 mm per raggiungere un y+ di 30.


Approssimare la velocità di attrito dalla velocità media del flusso.

Per calcolare Y+ utilizzando la velocità del fluido, è possibile adottare un approccio leggermente diverso. Y+ è definito in termini di velocità di attrito (uτ), ma uτ può essere strettamente correlato alla velocità media del flusso in alcune situazioni, specialmente in flussi turbolenti completamente sviluppati lungo superfici piane. Tuttavia, si deve tenere presente che questa correlazione può non essere precisa in tutti i casi, specialmente in geometrie complesse o in flussi non turbolenti.

La formula per Y+ rimane la stessa:

y+=uτyν

Per collegare la velocità di attrito alla velocità del fluido, si può utilizzare una relazione empirica o teorica basata sullo specifico caso di flusso. Ad esempio, in un canale o in un flusso turbolento lungo una piastra piana, la velocità di attrito può essere stimata come una frazione della velocità del flusso basata sul numero di Reynolds.

Una relazione approssimativa in tali casi è:

uτURe

Dove:

  • U = velocità media del flusso
  • Re = numero di Reynolds basato sulla lunghezza caratteristica

Questa è una stima molto approssimativa e dovrebbe essere usata con cautela. Per calcoli più accurati, è consigliato utilizzare i risultati di simulazioni precedenti o dati sperimentali per ottenere una stima più precisa di uτ.


Calcolo del numero di Reynolds

Il numero di Reynolds è un parametro non dimensionale utilizzato nella meccanica dei fluidi per prevedere i pattern di flusso in diversi tipi di problemi fluidodinamici. È definito come:

Re=ρULμ

dove:

  • ρ = densità del fluido (kg/m³)
  • U = velocità caratteristica del flusso (m/s)
  • L = lunghezza caratteristica (m), ad esempio il diametro di un tubo o la lunghezza di un'ala
  • μ = viscosità dinamica del fluido (Pa·s)

Esempio: Se hai un flusso d'aria (densità ρ=1.225 kg/m³, viscosità μ=1.81×105 Pa·s) sopra un'ala con una lunghezza caratteristica di 1 m e una velocità del flusso di 30 m/s, il numero di Reynolds sarà:

Re=1.225×30×11.81×1052.03×106

 


Esempio Conclusivo

Supponiamo di voler avere un valore di Y+ = 30 in un flusso d'aria con le seguenti proprietà:

  • Velocità media del flusso, U=15 m/s.
  • Viscosità cinematica dell'aria, ν=1.5×105 m²/s.
  • Numero di Reynolds basato su una lunghezza caratteristica, Re=106.

Prima calcoliamo uτ:

uτ15106=0.015 m/s

Ora usiamo la formula rielaborata di Y+ per trovare y:

y=30×1.5×1050.015=3×105 m=0.03 mm

Quindi, per ottenere un valore di Y+ di 30, la distanza dalla parete alla prima cella della mesh dovrebbe essere circa 0.03 mm.

Ricorda che questi calcoli sono basati su stime e approssimazioni. La precisione può variare a seconda della complessità del flusso e della geometria specifica del caso in esame. In una simulazione CFD reale, potrebbero essere necessarie ulteriori iterazioni per ottimizzare la dimensione della mesh e raggiungere il valore desiderato di Y+.


 Suggerimenti

Ora che hai capito la teoria alla base del dimensionamento della mesh per una simulazione di openFOAM ti suggerisco di utilizzare il pratico calcolatore automatico sviluppato da CFDFEASERVICE per l'utilizzo con snappyHexMesh di openFOAM.

Vuoi saperne di più su questo argomento?

 

pisa profile Andrea          

Vuoi approfondire le tematiche presentate nelle nostre simulazioni o nei casi studio? Contattaci liberamente per una discussione dettagliata o per qualsiasi domanda tu abbia.

Siamo pronti a offrirti il nostro supporto e la nostra esperienza.

Clicca qui per inviarci un messaggio e ti risponderemo al più presto.

L.THEME
We use cookies

Utilizziamo i cookie sul nostro sito Web. Alcuni di essi sono essenziali per il funzionamento del sito, mentre altri ci aiutano a migliorare questo sito e l'esperienza dell'utente (cookie di tracciamento). Puoi decidere tu stesso se consentire o meno i cookie. Ti preghiamo di notare che se li rifiuti, potresti non essere in grado di utilizzare tutte le funzionalità del sito.