la micrometeorologia e la dispersione degli inquinanti ... - ARPA Lazio

MODELLO EULERIANO
5.2.5 Considerazioni conclusive
La prima considerazione da fare è che l’ordine con cui si realizza lo schema di
frazionamento del tempo non è rigido e si può organizzare in maniere differenti.
Spesso si utilizza una sequenza (per esempio quella presentata nei paragrafi
precedenti) ad un certo passo temporale e nel passo successivo la si inverte e così
nei passi successivi. Questo metodo consente di minimizzare le instabilità ancora
presenti nella soluzione numerica.
Un elemento importante nello schema numerico è costituito dalle condizioni
iniziali, cioè lo stabilire quale concentrazione caratterizza ciascuna cella del
dominio di calcolo all’inizio della simulazione. Come si è visto la determinazione
delle condizioni iniziali dipende dalle informazioni che risultano disponibili:
• se si dispone di dati sperimentali rilevati in alcuni punti del dominio di calcolo,
si può interpolare questi valori in tutto il dominio e utilizzare questa interpolazione
come condizione iniziale;
• se non si hanno informazioni, si può porre a zero il campo di concentrazione
iniziale e si può ripetere la simulazione dell’istante temporale iniziale finché
non si raggiunge una situazione di equilibrio.
Delle condizioni al contorno già si è parlato in precedenza.
Da ultimo vanno considerati i coefficienti di diffusività turbolenta K xx ,K yy e K zz .
Sono molte le relazioni proposte per la loro stima. Qui vengono presentate solo
quelle che tengono conto esplicitamente delle differenti condizioni di turbolenza
del PBL.
Per quanto riguarda K zz sono disponibili forme funzionali differenti a seconda
dello stato di turbolenza del PBL. In particolare:
• nelle situazioni convettive
dove la funzione f(z/zi) risulta pari a:
quando 0.6 > z/z i ≥ 0.1,
quando 1.1 > z/z i ≥ 0.6 da
quando z/zi > 1.1 da
308
• nelle situazioni adiabatiche la relazione per K zz è la seguente: