Roberto Sozzi (ARPA Lazio) Teodoro Georgiadis (CNR-IBIMET ...
Roberto Sozzi (ARPA Lazio) Teodoro Georgiadis (CNR-IBIMET ... Roberto Sozzi (ARPA Lazio) Teodoro Georgiadis (CNR-IBIMET ...
5.ANALISI SPETTRALE DELLA TURBOLENZA DEL PBL.—————————————————————⎯⎯——————————————pari a due volte l’energia cinetica turbolenta E che, nell’ambito dell’analisi di Fourier, porta a:3∑ u uiuiuiR ( 0) = 2E= ∫ ∫ ∫∑ B ( k ) ⋅ dk = ∫ ΨE( k)i=1∞ 3∞1i2 ⋅ dk[5.15]2 3− ∞ i=1 0dove u i è la generica componente del vettore vento, k 2 =k 1 2 +k 2 2 +k 3 2 e Ψ E è detto densità spettraletridimensionale dell’energia cinetica.Ovviamente non è facile misurare uno spettro tridimensionale e quindi l’attenzione viene focalizzataprevalentemente su ciò che realisticamente è possibile misurare nel PBL. In effetti solo gli spettrimonodimensionali (longitudinale e trasversali), che dipendono dal numero d’onda in una particolaredirezione, possono essere agevolmente misurati e su di essi concentriamo la nostra attenzione, anche senon sono completamente appropriati per descrivere la struttura tridimensionale della turbolenza del PBL(Sorbjan, 1989). In particolare, se si considera la direzione del vento medio (k 1 ), lo spettrolongitudinale e laterale possono essere definiti come:HH11221π∞∫−∞( k ) R ( r,0,0)1uu − r= ⋅ ejk 1dr1π∞∫−∞( k ) R ( r,0,0)1vv − r= ⋅ ejk 1dr[5.16a][5.16b]Se si è in condizioni isotropiche, k 1 può essere sostituito da k (da non confondere nelle formule con lacostante di von Karman).Consideriamo ora con maggior dettaglio la densità spettrale dell’energia cinetica Ψ E presentataschematicamente nella Fig.5.1. Essa rappresenta la distribuzione dimensionale dei vortici presenti entroil PBL, dato che Ψ E rappresenta il contributo all’energia cinetica totale portato dalle armonichecaratterizzate da un numero d’onda compreso tra k e k+dk. Per comprendere meglio il suosignificato è opportuno richiamare quanto già presentato alla fine del Cap.1. Kolmogorov ha mostratocome lo spettro dell’energia cinetica possa essere suddiviso qualitativamente in tre zone distinte(Trombetti e Tagliazucca, 1994):1) la zona dei vortici di maggior dimensione a carattere permanente che acquisiscono l’energiacinetica dall’instabilità dei moti a grande scala derivante dallo shear e dal galleggiamento.Pertanto, le caratteristiche della densità spettrale dipendono dalle condizioni di formazione (velocitàdel vento, stabilità, rugosità superficiale ed estensione verticale del PBL) ed in questa zone laturbolenza è permanente, cioè ∂Ψ E (k,t)/∂t è piccolo.2) Zona dei vortici di grandi dimensioni contenenti energia; in questa zona l’energia cineticapassa dai vortici a maggior dimensione a quelli a dimensione inferiore senza perdite dissipative. Ladensità spettrale ad ogni numero d’onda dipende quindi dall’energia che transita in questa zona. Latipica lunghezza di scala nello Strato Superficiale è dell’ordine dell’estensione verticale del PBL.3) Zona di isotropia locale (o zona di equilibrio universale). In questa zona i vortici di dimensioneinferiore ad una lunghezza d’onda critica λ c hanno subito un decadimento sufficiente ad avercompletamente perso le caratteristiche dei vortici più grandi da cui hanno avuto origine. Così, i———————————————————————⎯⎯————————————- 179 -
5.ANALISI SPETTRALE DELLA TURBOLENZA DEL PBL.—————————————————————⎯⎯——————————————vortici di dimensione inferiore a λ c presentano proprietà statistiche indipendenti dai vortici didimensioni superiori e sono caratterizzati da isotropia, omogeneità e stazionarietà locali. Inquesta zona, quindi, le caratteristiche statistiche della turbolenza saranno universali, cioèindipendenti dalle condizioni esterne e determinate unicamente dalla viscosità dinamica ν e daltasso di dissipazione dell’energia cinetica turbolenta ε (Prima ipotesi di Kolmogorov). Questazona non è completamente omogenea, ma può essere suddivisa in due sottozone distinte:3a) Inertial subrange: in questa zona il rapporto tra le forze di inerzia e le forze viscose (cioè ilNumero di Reynolds) è così elevato che le forze viscose, e quindi ν, possono esseretrascurate. Il comportamento della turbolenza dipenderà esclusivamente da ε (Secondaipotesi di Kolmogorov) e dall’analisi Dimensionale si ha che:2 3 −53( k) = C kΨE k ε [5.17a]in cui C k è la costante universale adimensionale di Kolmogorov. Questa relazione vale finchéla lunghezza caratteristica del vortice risulta inferiore a µ, microscala di Kolmogorov, cherappresenta la tipica dimensione dei vortici che dissipano energia in calore (dell’ordine deimillimetri) e che è definita come:( ν ) 3 ε1 4µ = [5.17b]3b) viscous subrange: intervallo spettrale in cui l’energia dei vortici di dimensione piccolissimaviene irrimediabilmente trasformata in moti molecolari e rimossa dal sistema sotto forma dicalore.Vortici dipendenti dallacondizione di formazioneVortici indipendentidalla condizione di formazioneVortici a caratterepermanenteVortici altamenteenergeticiZona di Equilibrio UniversaleRange InerzialeDissipazione viscosaFig. 5.1: rappresentazione schematica della densità spettrale dell’energia cinetica (Trombetti eTagliazucca, 1994).———————————————————————⎯⎯————————————- 180 -
- Page 141 and 142: 3. ANALISI ENERGETICA DEL PLANETARY
- Page 143 and 144: 3. ANALISI ENERGETICA DEL PLANETARY
- Page 145 and 146: 3. ANALISI ENERGETICA DEL PLANETARY
- Page 147 and 148: 3. ANALISI ENERGETICA DEL PLANETARY
- Page 149 and 150: 3. ANALISI ENERGETICA DEL PLANETARY
- Page 151 and 152: 4. TEORIA DELLA SIMILARITÀ——
- Page 153 and 154: 4. TEORIA DELLA SIMILARITÀ——
- Page 155 and 156: 4. TEORIA DELLA SIMILARITÀ——
- Page 157 and 158: 4. TEORIA DELLA SIMILARITÀ——
- Page 159 and 160: 4. TEORIA DELLA SIMILARITÀ——
- Page 161 and 162: 4. TEORIA DELLA SIMILARITÀ——
- Page 163 and 164: 4. TEORIA DELLA SIMILARITÀ——
- Page 165 and 166: 4. TEORIA DELLA SIMILARITÀ——
- Page 167 and 168: 4. TEORIA DELLA SIMILARITÀ——
- Page 169 and 170: 4. TEORIA DELLA SIMILARITÀ——
- Page 171 and 172: 4. TEORIA DELLA SIMILARITÀ——
- Page 173 and 174: 4. TEORIA DELLA SIMILARITÀ——
- Page 175 and 176: 4. TEORIA DELLA SIMILARITÀ——
- Page 177 and 178: 4. TEORIA DELLA SIMILARITÀ——
- Page 179 and 180: 4. TEORIA DELLA SIMILARITÀ——
- Page 181 and 182: 4. TEORIA DELLA SIMILARITÀ——
- Page 183 and 184: 4. TEORIA DELLA SIMILARITÀ——
- Page 185 and 186: 4. TEORIA DELLA SIMILARITÀ——
- Page 187 and 188: 5.ANALISI SPETTRALE DELLA TURBOLENZ
- Page 189 and 190: 5.ANALISI SPETTRALE DELLA TURBOLENZ
- Page 191: 5.ANALISI SPETTRALE DELLA TURBOLENZ
- Page 195 and 196: 5.ANALISI SPETTRALE DELLA TURBOLENZ
- Page 197 and 198: 5.ANALISI SPETTRALE DELLA TURBOLENZ
- Page 199 and 200: 5.ANALISI SPETTRALE DELLA TURBOLENZ
- Page 201 and 202: 5.ANALISI SPETTRALE DELLA TURBOLENZ
- Page 203 and 204: 5.ANALISI SPETTRALE DELLA TURBOLENZ
- Page 205 and 206: 5.ANALISI SPETTRALE DELLA TURBOLENZ
- Page 207 and 208: 5.ANALISI SPETTRALE DELLA TURBOLENZ
- Page 209 and 210: 5.ANALISI SPETTRALE DELLA TURBOLENZ
- Page 211 and 212: 5.ANALISI SPETTRALE DELLA TURBOLENZ
- Page 213 and 214: 6. LA STRUTTURA DEL PBL IN CONDIZIO
- Page 215 and 216: 6. LA STRUTTURA DEL PBL IN CONDIZIO
- Page 217 and 218: 6. LA STRUTTURA DEL PBL IN CONDIZIO
- Page 219 and 220: 6. LA STRUTTURA DEL PBL IN CONDIZIO
- Page 221 and 222: 6. LA STRUTTURA DEL PBL IN CONDIZIO
- Page 223 and 224: 6. LA STRUTTURA DEL PBL IN CONDIZIO
- Page 225 and 226: 6. LA STRUTTURA DEL PBL IN CONDIZIO
- Page 227 and 228: 6. LA STRUTTURA DEL PBL IN CONDIZIO
- Page 229 and 230: 6. LA STRUTTURA DEL PBL IN CONDIZIO
- Page 231 and 232: 6. LA STRUTTURA DEL PBL IN CONDIZIO
- Page 233 and 234: 6. LA STRUTTURA DEL PBL IN CONDIZIO
- Page 235 and 236: 6. LA STRUTTURA DEL PBL IN CONDIZIO
- Page 237 and 238: 6. LA STRUTTURA DEL PBL IN CONDIZIO
- Page 239 and 240: 6. LA STRUTTURA DEL PBL IN CONDIZIO
- Page 241 and 242: 6. LA STRUTTURA DEL PBL IN CONDIZIO
5.ANALISI SPETTRALE DELLA TURBOLENZA DEL PBL.—————————————————————⎯⎯——————————————pari a due volte l’energia cinetica turbolenta E che, nell’ambito dell’analisi di Fourier, porta a:3∑ u uiuiuiR ( 0) = 2E= ∫ ∫ ∫∑ B ( k ) ⋅ dk = ∫ ΨE( k)i=1∞ 3∞1i2 ⋅ dk[5.15]2 3− ∞ i=1 0dove u i è la generica componente del vettore vento, k 2 =k 1 2 +k 2 2 +k 3 2 e Ψ E è detto densità spettraletridimensionale dell’energia cinetica.Ovviamente non è facile misurare uno spettro tridimensionale e quindi l’attenzione viene focalizzataprevalentemente su ciò che realisticamente è possibile misurare nel PBL. In effetti solo gli spettrimonodimensionali (longitudinale e trasversali), che dipendono dal numero d’onda in una particolaredirezione, possono essere agevolmente misurati e su di essi concentriamo la nostra attenzione, anche senon sono completamente appropriati per descrivere la struttura tridimensionale della turbolenza del PBL(Sorbjan, 1989). In particolare, se si considera la direzione del vento medio (k 1 ), lo spettrolongitudinale e laterale possono essere definiti come:HH11221π∞∫−∞( k ) R ( r,0,0)1uu − r= ⋅ ejk 1dr1π∞∫−∞( k ) R ( r,0,0)1vv − r= ⋅ ejk 1dr[5.16a][5.16b]Se si è in condizioni isotropiche, k 1 può essere sostituito da k (da non confondere nelle formule con lacostante di von Karman).Consideriamo ora con maggior dettaglio la densità spettrale dell’energia cinetica Ψ E presentataschematicamente nella Fig.5.1. Essa rappresenta la distribuzione dimensionale dei vortici presenti entroil PBL, dato che Ψ E rappresenta il contributo all’energia cinetica totale portato dalle armonichecaratterizzate da un numero d’onda compreso tra k e k+dk. Per comprendere meglio il suosignificato è opportuno richiamare quanto già presentato alla fine del Cap.1. Kolmogorov ha mostratocome lo spettro dell’energia cinetica possa essere suddiviso qualitativamente in tre zone distinte(Trombetti e Tagliazucca, 1994):1) la zona dei vortici di maggior dimensione a carattere permanente che acquisiscono l’energiacinetica dall’instabilità dei moti a grande scala derivante dallo shear e dal galleggiamento.Pertanto, le caratteristiche della densità spettrale dipendono dalle condizioni di formazione (velocitàdel vento, stabilità, rugosità superficiale ed estensione verticale del PBL) ed in questa zone laturbolenza è permanente, cioè ∂Ψ E (k,t)/∂t è piccolo.2) Zona dei vortici di grandi dimensioni contenenti energia; in questa zona l’energia cineticapassa dai vortici a maggior dimensione a quelli a dimensione inferiore senza perdite dissipative. Ladensità spettrale ad ogni numero d’onda dipende quindi dall’energia che transita in questa zona. Latipica lunghezza di scala nello Strato Superficiale è dell’ordine dell’estensione verticale del PBL.3) Zona di isotropia locale (o zona di equilibrio universale). In questa zona i vortici di dimensioneinferiore ad una lunghezza d’onda critica λ c hanno subito un decadimento sufficiente ad avercompletamente perso le caratteristiche dei vortici più grandi da cui hanno avuto origine. Così, i———————————————————————⎯⎯————————————- 179 -