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.
5.1 AVVIO IN CONFIGURAZIONE SAVONIUS 80<br />
Figura 5.5: Volumi ideali per il controllo volumetrico della mesh<br />
Considerando le simulazioni numeriche e i relativi risultati<br />
ottenuti da Di Paolo (2007), è stato <strong>ad</strong>ottato per la configurazione<br />
d’avvio il modello <strong>di</strong> turbolenza Standard k-ε.<br />
Per quanto riguarda la scelta del time-step, dal punto <strong>di</strong> vista Time-step<br />
numerico vale l’in<strong>di</strong>cazione riportata nelle sezioni precedenti che<br />
prevede un valore limite per il numero <strong>di</strong> Courant nel campo<br />
pari a circa 10. Tuttavia nel caso in esame, la modellizzazione<br />
6-DOF introduce un ulteriore vincolo sulla scelta del time-step:<br />
al fine <strong>di</strong> produrre soluzioni in<strong>di</strong>pendenti dalla <strong>di</strong>mensione del<br />
time-step quest’ultimo deve risultare sufficientemente piccolo e<br />
tale da garantire che in un singolo passo temporale l’interfaccia<br />
ruoti per <strong>una</strong> <strong>di</strong>stanza inferiore della <strong>di</strong>mensione caratteristica <strong>di</strong><br />
cella locale. Nel nostro caso la <strong>di</strong>mensione caratteristica <strong>di</strong> cella<br />
(L) in corrispondenza dell’interfaccia è <strong>di</strong> L = 0, 007 m mentre la<br />
<strong>di</strong>stanza dall’<strong>asse</strong> <strong>di</strong> rotazione dell’interfaccia è <strong>di</strong> Ri = 0, 23 m.<br />
Al fine <strong>di</strong> garantire la proprietà sopra citata, il time step limite,<br />
ovvero quello per il quale durante un singolo passo temporale la<br />
mesh si muove <strong>di</strong> <strong>una</strong> <strong>di</strong>stanza pari esattamente alla <strong>di</strong>mensione<br />
locale <strong>di</strong> cella, può essere calcolato come<br />
∆t = ∆θ<br />
ω<br />
= L<br />
Riω<br />
Poiché la velocità angolare varierà da zero fino <strong>ad</strong> un valore<br />
incognito, scegliamo un time-step iniziale <strong>di</strong> ∆t = 0, 0005 s che ci<br />
assicura la con<strong>di</strong>zione <strong>di</strong> sicurezza sopra citata fino a velocità an-