Impianti Eolici -...
Transcript of Impianti Eolici -...
PO
LI
di M
I
te
cn
ic
o
la
no
Corso di
Impianti Eolici
Alessandro CroceDipartimento di Scienze e Tecnologie Aerospaziali
Politecnico di MilanoMilano
Anno Accademico 2017-18
Cors
o d
i Im
pia
nti
Eolici–
a.a
. 2
01
7-18
POLI-Wind Research Lab
Principi di Aerodinamica
Aerodinamica: studio dell’interazione tra solidi e gas in moto relativo
Definizioni:
Approccio Lagrangiano ►linee di corrente (pathlines)
Approccio Euleriano ►linee di flusso (streamlines)
Trasformazione Galileana
Moto stazionario
Fluido comprimibile VS incomprimibile
Cors
o d
i Im
pia
nti
Eolici–
a.a
. 2
01
7-18
POLI-Wind Research Lab
Risultante azioni aerodinamiche su un corpo avente superficie esterna S:
poiché t (sforzo) ha un contributo di pressione e uno di attrito:
e quindi:
Si ha (integrando):
Principi di Aerodinamica
𝑭 = න𝑆
𝝉 𝑑𝐴
𝝉 = 𝝉(𝑝, 𝒏, 𝒗, 𝜇)
𝝉 = 𝝉𝒑 + 𝝉𝝁
𝑭 = 𝑭(𝑝∞, 𝑽∞, 𝜇∞, 𝑆)
Cors
o d
i Im
pia
nti
Eolici–
a.a
. 2
01
7-18
POLI-Wind Research Lab
Ogni equazione fisica del tipo
è esprimibile nella forma
Ove sono numeri adimensionali funzioni delle variabili attraverso
le equazioni del tipo:
Teorema di Buckingham (π)
Teorema dell’analisi dimensionale, esso asserisce che ogni equazione fisica,
dipendente da n variabili fisiche, esprimibili in termini di k quantità fisiche
fondamentali indipendenti, è rappresentabile come funzione di (n-k) variabili
adimensionali costruite con le variabili originali.
Cors
o d
i Im
pia
nti
Eolici–
a.a
. 2
01
7-18
POLI-Wind Research Lab
Quantità primitive Quantità aerodinamiche
Coefficiente di viscosità
Lunghezza
Tempo
Massa
F: Forza viscosa
Teorema di Buckingham (π)
Cors
o d
i Im
pia
nti
Eolici–
a.a
. 2
01
7-18
POLI-Wind Research Lab
Considerando una generica forza aerodinamica F
e applicando il modello di gas ideale:
e ricordando che la velocità del suono è
si ha che (semplificando la notazione e passando alla forma scalare):
Da cui si può ipotizzare la seguente relazione generica:
Gli esponenti di tale espressione si ottengono applicando il teorema di Buckingham
Teorema di Buckingham (π)
𝑭 = 𝑭(𝑝∞, 𝑽∞, 𝜇∞, 𝑆)
𝐹 = 𝐹(𝜌, 𝑎, 𝑉, 𝛼, 𝜇, 𝑆)
𝐹 = 𝑘(𝛼) 𝜌𝑒𝜌 𝑎𝑒𝑎 𝑉𝑒𝑉 𝜇𝑒𝜇 𝑆𝑒𝑆
𝑝 = 𝜌𝑅𝑇 𝑎 = 𝛾𝑅𝑇
Cors
o d
i Im
pia
nti
Eolici–
a.a
. 2
01
7-18
POLI-Wind Research Lab
Dalle quantità primitive:
Da cui:
E risolvendo rispetto 𝑒𝜇 ed 𝑒𝑎
Teorema di Buckingham (π)
𝐹 = 𝑘(𝛼) 𝜌𝑒𝜌 𝑎𝑒𝑎 𝑉𝑒𝑉 𝜇𝑒𝜇 𝑆𝑒𝑆
[𝑀𝐿𝑇−2] = [𝑀𝐿−3]𝑒𝜌 [𝐿𝑇−1]𝑒𝑎 [𝐿𝑇−1]𝑒𝑉 [𝑀𝐿−1𝑇−1]𝑒𝜇 [𝐿2]𝑒𝑆
𝑀 → 1 = 𝑒𝜌+𝑒𝜇
𝐿 → 1 = −3𝑒𝜌+𝑒𝑎 + 𝑒𝑉 − 𝑒𝜇 + 2𝑒𝑆
𝑇 → −2 =-𝑒𝑎 − 𝑒𝑉 − 𝑒𝜇
𝐹 = 𝑘(𝛼) 𝜌𝑉2𝑆𝜌𝑉 𝑆
𝜇
−𝑒𝜇𝑉
𝑎
−𝑒𝑎
Cors
o d
i Im
pia
nti
Eolici–
a.a
. 2
01
7-18
POLI-Wind Research Lab
Ovvero definendo i DUE (come previsto dal teorema 𝜋) parametri adimensionali:
numero di Reynolds
numero di Mach
essendo 𝐶𝐹 il coefficiente adimensionale di forza aerodinamica
Teorema di Buckingham (π)
𝐹 = 𝑘(𝛼) 𝜌𝑉2𝑆𝜌𝑉 𝑆
𝜇
−𝑒𝜇𝑉
𝑎
−𝑒𝑎
𝐹 =1
2𝜌𝑉2𝑆𝐶𝐹 𝛼, 𝑅𝑒,𝑀
𝑀 ≔𝑉
𝑎
𝑅𝑒 ≔𝜌𝑉𝑙
𝜇
Cors
o d
i Im
pia
nti
Eolici–
a.a
. 2
01
7-18
POLI-Wind Research Lab
α: angolo di incidenza
c: corda1: linea di portanza nulla2: bordo d'attacco3: cerchio osculatore del bordo d'attacco4: curvatura5: spessore6: dorso7: bordo d'uscita8: linea d'inarcamento media9: ventre.
Profili aerodinamici
Cors
o d
i Im
pia
nti
Eolici–
a.a
. 2
01
7-18
POLI-Wind Research Lab
Nomenclatura (ENG)
Source: J.F. Manwell, J.G.McGowan,A.L.Rogers, Wind Energy Explained, Theory, Design and Application, Wiley, 2002
Profili aerodinamici
Cors
o d
i Im
pia
nti
Eolici–
a.a
. 2
01
7-18
POLI-Wind Research Lab
Lift, Drag, Moment
Source: J.F. Manwell, J.G.McGowan,A.L.Rogers, Wind Energy Explained, Theory, Design and Application, Wiley, 2002
Profili aerodinamici
Cors
o d
i Im
pia
nti
Eolici–
a.a
. 2
01
7-18
POLI-Wind Research Lab
Profili aerodinamici
Cors
o d
i Im
pia
nti
Eolici–
a.a
. 2
01
7-18
POLI-Wind Research Lab
Caratteristiche aerodinamiche
Cors
o d
i Im
pia
nti
Eolici–
a.a
. 2
01
7-18
POLI-Wind Research Lab
Profili per turbine eoliche. Esempio: DU Wind Airfoils
Source: http://gcep.stanford.edu/pdfs/energy_workshops_04_04/wind_van_rooij.pdf
These airfoils are used by various wind turbine manufacturers world wide (GAMESA GE‐Wind, REpower, Dewind, Suzlon, Gamesa, LM Glasfiber, NOI Rotortechnik, Fuhrlander, Pfleiderer, EUROS, NEG Micon, Umoeblades, Ecotecnia ) in over 10 different rotor blades for turbines with rotor diameters ranging from 29 m to over 100 m, corresponding to machines with maximum power ranging from 350 kW to 3.5 MW
Profili aerodinamici
Cors
o d
i Im
pia
nti
Eolici–
a.a
. 2
01
7-18
POLI-Wind Research Lab
Source: experimental data + Viterna extrapolationAngle of attack [deg]
Profili aerodinamici
Cors
o d
i Im
pia
nti
Eolici–
a.a
. 2
01
7-18
POLI-Wind Research Lab
Source: experimental data + Viterna extrapolationAngle of attack [deg]
Profili aerodinamici
Cors
o d
i Im
pia
nti
Eolici–
a.a
. 2
01
7-18
POLI-Wind Research Lab
Influenza del numero di Reynolds
Cors
o d
i Im
pia
nti
Eolici–
a.a
. 2
01
7-18
POLI-Wind Research Lab
Galleria del vento1.4MW Civil-Aeronautical Wind Tunnel(CAWT):
• 13.8x3.8m, 14m/s, civil section:
- turbulence index < 2%
- turbulence generators = 25%
- 13m turntable
• 4x3.8m, 55m/s, aero section:
- turbulence <0.1%
- open-closed test section
Controllable motors, complete flow measurement and data acquisition systems
Cors
o d
i Im
pia
nti
Eolici–
a.a
. 2
01
7-18
POLI-Wind Research Lab
Galleria del vento
Turbulence (boundary layer) generators
Turn-table
13 m
• Low speed testing in the presence of vertical wind profile
• Multiple wind turbine testing (wake-machine interaction)
• High speed testing
• Aerodynamic characterization (Cp-TSR-β & CF-TSR-β curves)
Cors
o d
i Im
pia
nti
Eolici–
a.a
. 2
01
7-18
POLI-Wind Research Lab
Galleria del vento
http://www.dnw.aero/mod_popup.htm?/uploadedImages/airfoil_in_test_sec_krgupdatervar1000F.jpg
http://www.scielo.org.ar/scielo.php?pid=S0327-07932005000300003&script=sci_arttext
Cors
o d
i Im
pia
nti
Eolici–
a.a
. 2
01
7-18
POLI-Wind Research Lab
V
a
V
V
a
Fonte: http://www.diam.unige.it/~irro/lecture.html
Campo di moto attorno ad un profilo
Cors
o d
i Im
pia
nti
Eolici–
a.a
. 2
01
7-18
POLI-Wind Research Lab
Un forte gradiente di recupero di pressione conduce alla transizione e causa separazione
Pressioni al bordo d’uscita comportano maggior recupero di pressione
Minimo valore di Cp comporta la massima supervelocità locale ed eventualmente formazione di onda d’urto, determina l’entità del recupero di pressione
Gradiente di pressione in aspirazione comporta un flusso laminare, bassa resistenza
V
Distribuzione e risultante degli sforzi agenti sul profilo
Cors
o d
i Im
pia
nti
Eolici–
a.a
. 2
01
7-18
POLI-Wind Research Lab
Freccia uscente dal profilo: pressione minore della pressione atmosferica p∞
Freccia entrante nel profilo: pressione maggiore della pressione atmosferica p∞
Versore normale
Contorno del profilo
Coefficiente di forza
Corda del profilo
Vettore sforzo
Pressione in condizioni statiche
Densità
Velocità
Distribuzione e risultante degli sforzi agenti sul profilo
Cors
o d
i Im
pia
nti
Eolici–
a.a
. 2
01
7-18
POLI-Wind Research Lab
Freccia uscente dal profilo: pressione minore della pressione atmosferica p∞
Freccia entrante nel profilo: pressione maggiore della pressione atmosferica p∞
Distribuzione e risultante degli sforzi agenti sul profilo al variare di a
Cors
o d
i Im
pia
nti
Eolici–
a.a
. 2
01
7-18
POLI-Wind Research Lab
www.av8n.com/irro/profilo_e.html
Distribuzione di velocità e pressione
Cors
o d
i Im
pia
nti
Eolici–
a.a
. 2
01
7-18
POLI-Wind Research Lab
1. La variazione di velocità di un fluido è causata da un gradiente di pressione, NON il contrario.
2. Un gradiente di pressione normale alle linee di flusso è generato dalla deflessione delle linee di flusso stesse.
Bernoulli
Cors
o d
i Im
pia
nti
Eolici–
a.a
. 2
01
7-18
POLI-Wind Research Lab
Equazione di continuità(conservazione massa):
Equazione di conservazione dell’energia:
Bernoulli
Cors
o d
i Im
pia
nti
Eolici–
a.a
. 2
01
7-18
POLI-Wind Research Lab
Bernoulli
– Flusso stazionario
– Incomprimibile
– Non viscoso
– No scambio calore
– No variazione quota
Quindi scambio pressione-velocità… non chiaro però causa-effetto
Bernoulli
Cors
o d
i Im
pia
nti
Eolici–
a.a
. 2
01
7-18
POLI-Wind Research Lab
Bernoulli – 2
Bernoulli
► Bernoulli deriva da equazione
di Newton, quindi l’accelerazione
del fluido deriva da un gradiente
di pressione.
Cors
o d
i Im
pia
nti
Eolici–
a.a
. 2
01
7-18
POLI-Wind Research Lab
Accelerazione normale
Bernoulli
Se le linee di flusso vengono
deflesse, nasce un gradiente di
pressione attraverso le linee di
flusso stesse, con la pressione
che aumenta allontanandosi dal
centro di curvatura.►
Cors
o d
i Im
pia
nti
Eolici–
a.a
. 2
01
7-18
POLI-Wind Research Lab
Accelerazione normale
Source: H. Babinsky, How do wings work? http://iopscience.iop.org/0031-9120/38/6/001/pdf/pe3_6_001.pdf
Bernoulli
pA=pC=pinf
pD>pC & pA>pB
pD>p
B
Cors
o d
i Im
pia
nti
Eolici–
a.a
. 2
01
7-18
POLI-Wind Research Lab
Strato limite e flusso viscoso
Cors
o d
i Im
pia
nti
Eolici–
a.a
. 2
01
7-18
POLI-Wind Research Lab
Separazione
Distribuzione e risultante degli sforzi agenti sul profilo
Cors
o d
i Im
pia
nti
Eolici–
a.a
. 2
01
7-18
POLI-Wind Research Lab
Transizione da laminare a turbolento
Cors
o d
i Im
pia
nti
Eolici–
a.a
. 2
01
7-18
POLI-Wind Research Lab
Transizione da laminare a turbolento
Cors
o d
i Im
pia
nti
Eolici–
a.a
. 2
01
7-18
POLI-Wind Research Lab
Transizione da laminare a turbolento
Cors
o d
i Im
pia
nti
Eolici–
a.a
. 2
01
7-18
POLI-Wind Research Lab
Flusso laminare (parte bassa nellafoto) e flusso turbolento (parte altanella foto) sullo scafo di un sommergibile
Transizione da flusso laminare a flussoturbolento su un profilo.
Transizione
Transizione da laminare a turbolento
Cors
o d
i Im
pia
nti
Eolici–
a.a
. 2
01
7-18
POLI-Wind Research Lab
Transizione
Source: http://www.aerospaceweb.org/question/aerodynamics/q0215.shtml
Cors
o d
i Im
pia
nti
Eolici–
a.a
. 2
01
7-18
POLI-Wind Research Lab
Incidenza di portanza nulla
Incidenza aerodinamica
Al di sotto dello stallo
Valori tipici
Valore teoricoTeoria dei profili sottili
Tratto lineare tipico
Caratteristiche aerodinamiche
Cors
o d
i Im
pia
nti
Eolici–
a.a
. 2
01
7-18
POLI-Wind Research Lab
Valori tipici
Non dipende significativamente dalla forma e non dipende dalla velocità
Dipende in misura significativa dalla forma
Caratteristiche aerodinamiche
Cors
o d
i Im
pia
nti
Eolici–
a.a
. 2
01
7-18
POLI-Wind Research Lab
Caratteristiche aerodinamicheProfilo spesso, stallo di bordo d’uscita
Profilo sottile, stallo di bordo d’attacco
Profilo spesso
Profilo sottile
Cors
o d
i Im
pia
nti
Eolici–
a.a
. 2
01
7-18
POLI-Wind Research Lab
DF = Da + Dp
Resistenza di profilo
Resistenza di attrito
Bassa per laminare
Alta per turbolento
Fortemente dipendente da Re
Resistenza di pressione o forma
Alta per laminare
Bassa per turbolento
Poco dipendente da Re
Caratteristiche aerodinamiche
Cors
o d
i Im
pia
nti
Eolici–
a.a
. 2
01
7-18
POLI-Wind Research Lab
influenza dello spessoreCaratteristiche aerodinamiche
Cors
o d
i Im
pia
nti
Eolici–
a.a
. 2
01
7-18
POLI-Wind Research Lab
Visualizzazione con fumo della scia vorticosa generata dalla portanza. Modello di ala in galleria del vento.
Risulta ben visibile alle estremità l’avvolgimento dei filetti vorticosi, che comporta una intensa concentrazione di vorticità.
Ala di apertura finita
Cors
o d
i Im
pia
nti
Eolici–
a.a
. 2
01
7-18
POLI-Wind Research Lab
Visualizzazione delle linee di corrente del campo di motoattorno ad un’ala di apertura finita
Nell'ala d’apertura finita, la differenzadi pressione tra dorso e ventre siannulla alle estremità alari,comportando un gradiente dipressione in direzione dell’aperturaalare stessa.
Come conseguenza si sviluppa unacomponente trasversale di velocità, disegno opposto fra dorso e ventre, che,a sua volta, origina uno stratovorticoso che si diparte dal bordod’uscita dell’ala.I filetti vorticosi di rotazione concordetendono ad avvolgersi, dando origine,in corrispondenza dell’estremità alare,a due zone a forma cilindrica di forteintensità vorticosa.
Resistenza indotta
Cors
o d
i Im
pia
nti
Eolici–
a.a
. 2
01
7-18
POLI-Wind Research Lab
Schema di strato vorticoso convortici concentrati alle estremità
Resistenza indotta
Cors
o d
i Im
pia
nti
Eolici–
a.a
. 2
01
7-18
POLI-Wind Research Lab
La presenza della scia vorticosaaltera il campo di moto nonsolamente a valle ma anche a montedell’ala che si trova ad essereinvestita da un flusso d’aria nonuniforme contestualmente allagenerazione della portanza.
Resistenza indotta
Cors
o d
i Im
pia
nti
Eolici–
a.a
. 2
01
7-18
POLI-Wind Research Lab
dL
dL
dD
Resistenza indotta
Cors
o d
i Im
pia
nti
Eolici–
a.a
. 2
01
7-18
POLI-Wind Research Lab
Scia di un rotore
http://www.oneplusplus.com/AeroDynamic2/WindTech.html
A rare, unretouched aerial photo of an offshore windfarm in Denmark clearly shows how turbulence generated by large turbine rotors continues to build with each successive row of turbines.