Roberto Sozzi (ARPA Lazio) Teodoro Georgiadis (CNR-IBIMET ...
Roberto Sozzi (ARPA Lazio) Teodoro Georgiadis (CNR-IBIMET ... Roberto Sozzi (ARPA Lazio) Teodoro Georgiadis (CNR-IBIMET ...
8. MODELLI NUMERICI DEL PBL—————————————————————⎯⎯——————————————8.1.3.4 L’equazione che contiene il termine di sorgenteSe si eliminano dalle due equazioni che rappresentano la conservazione della quantità di moto i terminiconvettivi ed avvettivi, si ottiene il seguente sistema che evidenzia il termine di sorgente dovutoall'accelerazione di Coriolis:∂u = fv∂t∂v= − fu∂t[8.33a][8.33b]La semplicità di questo sistema di equazioni ci indurrebbe ad adottare schemi risolutivi estremamentesemplici, che però si mostra essere completamente instabili. Seguendo Pielke (1984), uno schemastabile e semplice per risolvere tale sistema di equazioni è invece il seguente:uvn+1in+i− u∆t1− vi∆tnin= fvni= − fun+1i[8.34a][8.34b]In pratica, all’istante n+1 si opera nel modo seguente:• si calcola il nuovo campo di u impiegando i valori di u e v ottenuti al passo temporale n;• si calcola il nuovo campo di v, impiegando il campo di v ottenuto al passo temporale n ed il campodi u ottenuto a n+1.8.1.3.5 L’equazione del bilancio di massaL’importanza dell’equazione (8.5) sta nel fatto che è l’unica equazione da cui sia possibile ottenere lacomponente verticale media del vento. A tal proposito, consideriamo il fatto che dalla (8.5) è immediatoottenere:∂w⎡∂u∂v⎤= −⎢+∂⎥[8.35]z ⎣∂x∂y⎦Si consideri il solito dominio di calcolo in cui un generico nodo interno sia univocamente individuato dagliindici i (asse x), j (asse y) e k (asse z). Se impieghiamo uno schema forward per l’approssimazionedella derivata rispetto a z e uno schema centrato per le altre due derivate, alla fine otteniamo:wi,j,k + 1− w∆zi,j,k⎡u= −⎢⎢⎣i+1, j,k− u2∆xi−1,j,kv+i,j+1, k− v2∆yi,j−1,k⎤⎥⎥⎦[8.36]Se i nodi al suolo sono caratterizzati da k = 0, si ha per la no-slip condition che w 0 = 0 qualunquevalore assumano gli indici i e j. Alla luce di tutto ciò, se sono noti i campi u e v e i valori assunti da wai nodi di indice verticale k, il valore assunto da w ai nodi con indici k+1 sarà facilmente ottenutoimpiegando direttamente la (8.36).—————————————————————————————————————- 287 -
8. MODELLI NUMERICI DEL PBL—————————————————————⎯⎯——————————————8.1.3.6 L’operazione di splittingSe consideriamo l’equazione (8.4a) notiamo che in essa sono combinati tutti i tipi di equazionidifferenziali fin qui considerate singolarmente. Il problema che ci poniamo è se sia possibile in qualchemodo sfruttare quanto abbiamo fin qui sviluppato per risolvere numericamente questa equazione e lealtre che (ad eccezione dell’equazione del bilancio di massa) risultano formalmente identiche. Una dellerisposte possibili è l’operazione di time splitting. Per semplicità, consideriamo una generica equazionedifferenziale alle derivate parziali del tipo:∂u= Ψ∂t( u)[8.37]in cui Ψ sia un generico operatore differenziale, definito come la somma di operatori che agisconoadditivamente sulla variabile u nel modo seguente:( u) = Ψ ( u) + Ψ ( u ) + ⋅⋅ ⋅ + Ψ ( u)Ψ1 2m[8.38]Immaginiamo, inoltre, che per i vari operatori Ψ 1 , Ψ 2 , ..,Ψ m siano disponibili schemi di approssimazionestabili ed efficienti. In questo caso, si può dimostrare che la risoluzione della (8.37) può essere ottenutanel modo seguente.Si supponga di conoscere il campo di u all’istante n. Basandoci su tale campo si risolva(numericamente) l’equazione:∂u= Ψ 1∂t( u)[8.39a]da cui si ottiene un nuovo campo di u (che indichiamo come u (1) ) che però non sarà il campo volutoall’istante n+1.Usando u (1) come campo iniziale, se risolviamo l’equazione:∂u= Ψ∂t( 1())2 u[8.39b]otterremo un nuovo campo intermedio u (2) che terrà conto sia dell’azione dell’operatore Ψ 1 chedell’operatore Ψ 2 . Continuando questa metodologia fino all’esaurimento degli m operatori in cui è statopossibile scomporre l’operatore differenziale originario, alla fine si otterrà il campo u (m) che èeffettivamente il campo cercato per l’istante temporale n+1 (cioè u n+1 = u (m) ).Questo metodo rappresenta la soluzione ai nostri problemi. In effetti, se consideriamo ancora una voltala (8.4a), vediamo come la sequenza di operazioni richiesta dallo splitting sia la seguente:∂u∂u= −u∂t∂x∂u∂u= −v∂t∂y∂u∂u= −w∂t∂z[8.40a][8.40b][8.40c]—————————————————————————————————————- 288 -
- Page 249 and 250: 7. IL PBL IN SITUAZIONI SUPERFICIAL
- Page 251 and 252: 7. IL PBL IN SITUAZIONI SUPERFICIAL
- Page 253 and 254: 7. IL PBL IN SITUAZIONI SUPERFICIAL
- Page 255 and 256: 7. IL PBL IN SITUAZIONI SUPERFICIAL
- Page 257 and 258: 7. IL PBL IN SITUAZIONI SUPERFICIAL
- Page 259 and 260: 7. IL PBL IN SITUAZIONI SUPERFICIAL
- Page 261 and 262: 7. IL PBL IN SITUAZIONI SUPERFICIAL
- Page 263 and 264: 7. IL PBL IN SITUAZIONI SUPERFICIAL
- Page 265 and 266: 7. IL PBL IN SITUAZIONI SUPERFICIAL
- Page 267 and 268: 7. IL PBL IN SITUAZIONI SUPERFICIAL
- Page 269 and 270: 7. IL PBL IN SITUAZIONI SUPERFICIAL
- Page 271 and 272: 7. IL PBL IN SITUAZIONI SUPERFICIAL
- Page 273 and 274: 7. IL PBL IN SITUAZIONI SUPERFICIAL
- Page 275 and 276: 7. IL PBL IN SITUAZIONI SUPERFICIAL
- Page 277 and 278: 7. IL PBL IN SITUAZIONI SUPERFICIAL
- Page 279 and 280: 7. IL PBL IN SITUAZIONI SUPERFICIAL
- Page 281 and 282: 7. IL PBL IN SITUAZIONI SUPERFICIAL
- Page 283 and 284: 8. MODELLI NUMERICI DEL PBL——
- Page 285 and 286: 8. MODELLI NUMERICI DEL PBL——
- Page 287 and 288: 8. MODELLI NUMERICI DEL PBL——
- Page 289 and 290: 8. MODELLI NUMERICI DEL PBL——
- Page 291 and 292: 8. MODELLI NUMERICI DEL PBL——
- Page 293 and 294: 8. MODELLI NUMERICI DEL PBL——
- Page 295 and 296: 8. MODELLI NUMERICI DEL PBL——
- Page 297 and 298: 8. MODELLI NUMERICI DEL PBL——
- Page 299: 8. MODELLI NUMERICI DEL PBL——
- Page 303 and 304: 8. MODELLI NUMERICI DEL PBL——
- Page 305 and 306: 8. MODELLI NUMERICI DEL PBL——
- Page 307 and 308: 8. MODELLI NUMERICI DEL PBL——
- Page 309 and 310: 8. MODELLI NUMERICI DEL PBL——
- Page 311 and 312: 8. MODELLI NUMERICI DEL PBL——
- Page 313 and 314: 8. MODELLI NUMERICI DEL PBL——
- Page 315 and 316: 8. MODELLI NUMERICI DEL PBL——
- Page 317 and 318: 8. MODELLI NUMERICI DEL PBL——
- Page 319 and 320: 8. MODELLI NUMERICI DEL PBL——
- Page 321 and 322: 8. MODELLI NUMERICI DEL PBL——
- Page 323 and 324: 8. MODELLI NUMERICI DEL PBL——
- Page 325 and 326: 8. MODELLI NUMERICI DEL PBL——
- Page 327 and 328: 9 - TECNICHE PER L’OSSERVAZIONE D
- Page 329 and 330: 9 - TECNICHE PER L’OSSERVAZIONE D
- Page 331 and 332: 9 - TECNICHE PER L’OSSERVAZIONE D
- Page 333 and 334: 9 - TECNICHE PER L’OSSERVAZIONE D
- Page 335 and 336: 9 - TECNICHE PER L’OSSERVAZIONE D
- Page 337 and 338: 9 - TECNICHE PER L’OSSERVAZIONE D
- Page 339 and 340: 9 - TECNICHE PER L’OSSERVAZIONE D
- Page 341 and 342: 9 - TECNICHE PER L’OSSERVAZIONE D
- Page 343 and 344: 9 - TECNICHE PER L’OSSERVAZIONE D
- Page 345 and 346: 9 - TECNICHE PER L’OSSERVAZIONE D
- Page 347 and 348: 9 - TECNICHE PER L’OSSERVAZIONE D
- Page 349 and 350: 9 - TECNICHE PER L’OSSERVAZIONE D
8. MODELLI NUMERICI DEL PBL—————————————————————⎯⎯——————————————8.1.3.4 L’equazione che contiene il termine di sorgenteSe si eliminano dalle due equazioni che rappresentano la conservazione della quantità di moto i terminiconvettivi ed avvettivi, si ottiene il seguente sistema che evidenzia il termine di sorgente dovutoall'accelerazione di Coriolis:∂u = fv∂t∂v= − fu∂t[8.33a][8.33b]La semplicità di questo sistema di equazioni ci indurrebbe ad adottare schemi risolutivi estremamentesemplici, che però si mostra essere completamente instabili. Seguendo Pielke (1984), uno schemastabile e semplice per risolvere tale sistema di equazioni è invece il seguente:uvn+1in+i− u∆t1− vi∆tnin= fvni= − fun+1i[8.34a][8.34b]In pratica, all’istante n+1 si opera nel modo seguente:• si calcola il nuovo campo di u impiegando i valori di u e v ottenuti al passo temporale n;• si calcola il nuovo campo di v, impiegando il campo di v ottenuto al passo temporale n ed il campodi u ottenuto a n+1.8.1.3.5 L’equazione del bilancio di massaL’importanza dell’equazione (8.5) sta nel fatto che è l’unica equazione da cui sia possibile ottenere lacomponente verticale media del vento. A tal proposito, consideriamo il fatto che dalla (8.5) è immediatoottenere:∂w⎡∂u∂v⎤= −⎢+∂⎥[8.35]z ⎣∂x∂y⎦Si consideri il solito dominio di calcolo in cui un generico nodo interno sia univocamente individuato dagliindici i (asse x), j (asse y) e k (asse z). Se impieghiamo uno schema forward per l’approssimazionedella derivata rispetto a z e uno schema centrato per le altre due derivate, alla fine otteniamo:wi,j,k + 1− w∆zi,j,k⎡u= −⎢⎢⎣i+1, j,k− u2∆xi−1,j,kv+i,j+1, k− v2∆yi,j−1,k⎤⎥⎥⎦[8.36]Se i nodi al suolo sono caratterizzati da k = 0, si ha per la no-slip condition che w 0 = 0 qualunquevalore assumano gli indici i e j. Alla luce di tutto ciò, se sono noti i campi u e v e i valori assunti da wai nodi di indice verticale k, il valore assunto da w ai nodi con indici k+1 sarà facilmente ottenutoimpiegando direttamente la (8.36).—————————————————————————————————————- 287 -