Analisi numerica di una turbina eolica ad asse verticale - Atomino FVG
Analisi numerica di una turbina eolica ad asse verticale - Atomino FVG
Analisi numerica di una turbina eolica ad asse verticale - Atomino FVG
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
4.2 GENERALITÀ SULLE SIMULAZIONI NUMERICHE 50<br />
viscoso (y + < 5). Per tali modelli apposite funzioni tengono conto<br />
della transizione dalla zona dominata dagli sforzi viscosi alla<br />
zona dominata dagli effetti turbolenti: l’utilizzo <strong>di</strong> questi modelli<br />
permette pertanto <strong>di</strong> risolvere completamente la struttura dello<br />
strato limite, ma al tempo stesso è richiesto un elevato numero<br />
<strong>di</strong> celle. I modelli <strong>di</strong> tipo High-Reynolds number sono invece solitamente<br />
associati all’utilizzo <strong>di</strong> mesh meno fitte in prossimità<br />
delle pareti con il primo strato <strong>di</strong> celle a parete localizzate all’interno<br />
della regione logaritmica (y + > 30). Le informazioni sulle<br />
con<strong>di</strong>zioni <strong>di</strong> flusso per tali celle sono fornite dalle wall function.<br />
Esse sono rappresentate da un set <strong>di</strong> relazioni matematiche utilizzate<br />
per ottenere le con<strong>di</strong>zioni al contorno per le equazioni <strong>di</strong><br />
conservazione: il modello <strong>di</strong> turbolenza risulta pertanto valido<br />
solo all’esterno del sottostrato viscoso e quest’ultima zona non<br />
viene risolta esplicitamente.<br />
In tutte le simulazioni è stato usato il trattamento a parete<br />
all-y + che è un modello ibrido che cerca <strong>di</strong> emulare l’high-y + per<br />
mesh r<strong>ad</strong>e e il low-y + per mesh più fitte. Inoltre questo modello<br />
è stato formulato per dare soluzioni sod<strong>di</strong>sfacenti per mesh <strong>di</strong><br />
risoluzione interme<strong>di</strong>a, cioè quando il centro della cella c<strong>ad</strong>e nel<br />
buffer layer dello strato limite.<br />
Le con<strong>di</strong>zioni al contorno <strong>asse</strong>gnate per tutte le serie <strong>di</strong> simu- Con<strong>di</strong>zioni al<br />
lazioni sono le seguenti<br />
contorno e iniziali<br />
• Il valore <strong>di</strong> velocità per ogni simulazione viene fissato nella<br />
sezione <strong>di</strong> ingresso (velocity inlet) come con<strong>di</strong>zione al<br />
contorno e in tutto il dominio come con<strong>di</strong>zione iniziale.<br />
• La pressione viene fissata nella sezione <strong>di</strong> uscita (pressure<br />
outlet) come con<strong>di</strong>zione al contorno e in tutto il dominio come<br />
con<strong>di</strong>zione iniziale, per tutte le simulazione la pressione<br />
nella sezione <strong>di</strong> uscita è pari a 1 atm (101325 Pa).<br />
• Per gli elementi del dominio che rappresentano profili, assi<br />
e vele si <strong>asse</strong>gna la con<strong>di</strong>zione <strong>di</strong> wall. Questa con<strong>di</strong>zione<br />
equivale <strong>ad</strong> imporre le con<strong>di</strong>zioni <strong>di</strong> impermeabilità ed<br />
<strong>ad</strong>erenza del fluido alla parete (no-slip).<br />
• Nelle simulazioni 3D si assume che le sezioni superiore,<br />
inferiore e laterali siano piani <strong>di</strong> simmetria. Star-CCM+ attribuisce<br />
flusso nullo a tutte le grandezze <strong>di</strong> trasporto attraverso<br />
un piano <strong>di</strong> simmetria, non esiste cioè flusso convettivo<br />
(la componente della velocità normale al piano è nulla) e<br />
neanche flusso <strong>di</strong>ffusivo (i gra<strong>di</strong>enti normali al piano, <strong>di</strong><br />
tutte le variabili sono nulli). Dato che non esiste sforzo<br />
<strong>di</strong> taglio sul piano <strong>di</strong> simmetria, quest’ultimo può essere<br />
interpretato come <strong>una</strong> parete senza attrito (slip wall).