la micrometeorologia e la dispersione degli inquinanti ... - ARPA Lazio
la micrometeorologia e la dispersione degli inquinanti ... - ARPA Lazio la micrometeorologia e la dispersione degli inquinanti ... - ARPA Lazio
CAPITOLO 3 TEORIA DI BASE DELLA DISPERSIONE DEGLI INQUINANTI IN ATMOSFERA Fin qui si è concentrata tutta l’attenzione nello studio qualitativo e quantitativo del PBL, individuando un modello matematico capace di descriverne realisticamente le caratteristiche principali, in particolare la turbolenza. A questo punto diviene possibile cominciare lo studio della dispersione di sostanze (solide, aerosol e gassose) emesse entro il PBL, in modo da poterne stimare la concentrazione nello spazio e nel tempo. Nel seguito, per semplicità, tali sostanze emesse verranno indicate genericamente con il termine inquinanti, anche se spesso esse sono di origine totalmente naturale, come accade, per esempio, per le emissioni gassose e solide provenienti dai vulcani e per le spore ed i virus. Ci sono due approcci metodologici profondamente differenti per descrivere la dispersione degli inquinanti nel PBL. Il primo metodo prende il nome di approccio Euleriano. Secondo tale approccio, la dinamica della dispersione delle sostanze chimiche inquinanti è descritta avendo come riferimento spaziale un sistema fisso di coordinate cartesiane ortogonali solidale con la superficie terrestre. Questo è il modo consueto con cui vengono trattati in Fisica i fenomeni di trasferimento di calore e di massa e, nel caso della Micrometeorologia, anche il modo consueto usato per descrivere matematicamente i campi meteorologici. Come sarà chiaro nel seguito, il metodo euleriano cerca di formulare le principali variabili statistiche legate alla concentrazione dei diversi inquinanti sulla base delle proprietà statistiche delle velocità euleriane del fluido, cioè delle velocità misurate in un punto fisso del PBL. Una formulazione di questo tipo è molto utile non solo perché la statistica euleriana è facilmente misurabile con i normali strumenti usati in meteorologia, ma anche perché le espressioni matematiche che ne derivano sono immediatamente applicabili in situazioni in cui si ha la presenza di reazioni chimiche. Sfortunatamente, l’approccio euleriano comporta seri ostacoli. Il primo ostacolo risiede nel fatto che il modello matematico cui tale approccio dà origine non consente soluzioni analitiche sufficientemente generali ed immediatamente utilizzabili. Il secondo ostacolo, invece, è connesso con il problema della chiusura che può rendere difficile la rappresentazione realistica della dispersione degli inquinanti in un PBL fortemente convettivo. Oltre a ciò va rilevato che un modello di tipo euleriano normalmente si presenta come un’equazione differenziale alle derivate parziali che richiede quindi una soluzione numerica e, come noto, gli algoritmi numerici impiegati possono condizionare pesantemente i risultati finali. Il secondo approccio è l’approccio Lagrangiano in cui i cambiamenti di concentrazione sono descritti relativamente al moto del fluido. Ma come sarà più chiaro nel seguito, la peculiarità principale di tale approccio non è tanto il differente sistema di riferimento, quanto piuttosto l’atteggiamento che sta alla base di tale descrizione. Mentre nell’approccio euleriano si ipotizza di essere effettivamente in grado di descrivere con esattezza la fisica che sta alla base del trasporto e della dispersione degli inquinanti in aria, nell’approccio Lagrangiano si parte dal presupposto che la nostra ignoranza sulla fisica del fenomeno sia talmente grande da essere costretti ad usare la statistica come mezzo descrittivo. In pratica, le tecniche lagrangiane cercano di descrivere la statistica delle concentrazioni in termini di proprietà statistiche degli spostamenti di gruppo di particelle rilasciate nel fluido. La matematica di tale approccio, anche se non sembra a prima vista, è più trattabile di quella dei metodi Euleriani, non essendoci espliciti problemi di chiusura, ma l’applicabilità pratica diretta delle equazioni che ne derivano è limitata dal fatto che queste equazioni 203
TEORIA DI BASE DELLA DISPERSIONE DEGLI INQUINANTI IN ATMOSFERA non consentono una descrizione diretta delle reazioni chimiche non lineari spesso molto importanti nella normale dinamica del PBL. Come si vedrà nel seguito, i due approcci comportano differenti tipi di relazioni matematiche per descrivere la concentrazione degli inquinanti, espressioni matematiche differenti che però, in ultima analisi, possono essere messe in relazione le une con le altre. Ciascun approccio è una valida descrizione della dispersione turbolenta: la scelta di quale approccio adottare in una data situazione dovrà essere deciso caso per caso, nelle situazioni reali. Il questo capitolo vengono presentati, in termini abbastanza generali, i due differenti approcci, riferendoci alla celebre lezione di Lamb (1980) e a quanto presentato in Seinfeld e Pandis (1998). Nei capitoli successivi introdurremo poi differenti famiglie di modelli che possono essere applicati sulla pratica e che comunque derivano da uno dei due approcci. 3.1 I FONDAMENTI DELLA TEORIA EULERIANA. 3.1.1 La formulazione istantanea di base. L’ipotesi alla base dell’approccio euleriano è la sicurezza che le leggi della fluidodinamica siano in grado di descrivere in modo completo e perfetto (e quindi deterministico) il fenomeno della dispersione degli inquinanti nel PBL.E’questo, quindi, un approccio completamente deterministico ed in ultima analisi, si basa sulla scrittura dell’equazione di conservazione delle sostanze inquinanti presenti nel PBL. A differenza di quanto abbiamo fatto nel Cap.2, in questo caso non è più possibile considerare la miscela gassosa che costituisce il PBL come un’entità unica (che abbiamo chiamato genericamente aria) , ma al contrario siamo costretti: • a descrivere individualmente le diverse sostanze inquinanti a cui siamo interessati; • a tener conto del fatto che alcune di esse non sono gas inerti, ma sostanze più o meno reattive chimicamente; • a tener conto che alcune possono presentare un decadimento radioattivo; • a tener conto che alcune sono particelle solide o aerosol. Si consideri inizialmente un solo inquinante di tipo gassoso. Se si concentra l’attenzione su un piccolo volume d’aria (una particella, dunque) come mostrato in Fig.3.1 e si trascura la diffusività molecolare, è possibile cercare di applicare ad esso la conservazione della massa dell’inquinante considerato. In particolare, si consideri il flusso di inquinante che transita nel volume lungo la direzione x in un intervallo di tempo δt.Tale flusso è pari a: dove u 1 è la componente u della velocità dell’aria entrante nel volume attraverso la faccia di sinistra e u 2 è la velocità uscente dalla faccia di destra. Dato che δx è piccolo, è possibile realizzare la seguente espansione in serie di Taylor: Da queste considerazioni risulta che il flusso in transito è pari a: 204
- Page 153 and 154: LA MICROMETEOROLOGIA E LA CAPACITA
- Page 155 and 156: LA MICROMETEOROLOGIA E LA CAPACITA
- Page 157 and 158: LA MICROMETEOROLOGIA E LA CAPACITA
- Page 159 and 160: LA MICROMETEOROLOGIA E LA CAPACITA
- Page 161 and 162: LA MICROMETEOROLOGIA E LA CAPACITA
- Page 163 and 164: LA MICROMETEOROLOGIA E LA CAPACITA
- Page 165 and 166: LA MICROMETEOROLOGIA E LA CAPACITA
- Page 167 and 168: LA MICROMETEOROLOGIA E LA CAPACITA
- Page 169 and 170: LA MICROMETEOROLOGIA E LA CAPACITA
- Page 171 and 172: LA MICROMETEOROLOGIA E LA CAPACITA
- Page 173 and 174: LA MICROMETEOROLOGIA E LA CAPACITA
- Page 175 and 176: LA MICROMETEOROLOGIA E LA CAPACITA
- Page 177 and 178: LA MICROMETEOROLOGIA E LA CAPACITA
- Page 179 and 180: LA MICROMETEOROLOGIA E LA CAPACITA
- Page 181 and 182: LA MICROMETEOROLOGIA E LA CAPACITA
- Page 183 and 184: LA MICROMETEOROLOGIA E LA CAPACITA
- Page 185 and 186: LA MICROMETEOROLOGIA E LA CAPACITA
- Page 187 and 188: LA MICROMETEOROLOGIA E LA CAPACITA
- Page 189 and 190: LA MICROMETEOROLOGIA E LA CAPACITA
- Page 191 and 192: LA MICROMETEOROLOGIA E LA CAPACITA
- Page 193 and 194: LA MICROMETEOROLOGIA E LA CAPACITA
- Page 195 and 196: LA MICROMETEOROLOGIA E LA CAPACITA
- Page 197 and 198: LA MICROMETEOROLOGIA E LA CAPACITA
- Page 199 and 200: LA MICROMETEOROLOGIA E LA CAPACITA
- Page 201 and 202: LA MICROMETEOROLOGIA E LA CAPACITA
- Page 203: BIANCA
- Page 207 and 208: TEORIA DI BASE DELLA DISPERSIONE DE
- Page 209 and 210: TEORIA DI BASE DELLA DISPERSIONE DE
- Page 211 and 212: TEORIA DI BASE DELLA DISPERSIONE DE
- Page 213 and 214: TEORIA DI BASE DELLA DISPERSIONE DE
- Page 215 and 216: TEORIA DI BASE DELLA DISPERSIONE DE
- Page 217 and 218: TEORIA DI BASE DELLA DISPERSIONE DE
- Page 219 and 220: TEORIA DI BASE DELLA DISPERSIONE DE
- Page 221 and 222: TEORIA DI BASE DELLA DISPERSIONE DE
- Page 223 and 224: TEORIA DI BASE DELLA DISPERSIONE DE
- Page 225 and 226: TEORIA DI BASE DELLA DISPERSIONE DE
- Page 227 and 228: MODELLI DI TIPO STAZIONARIO 226 omo
- Page 229 and 230: MODELLI DI TIPO STAZIONARIO il camm
- Page 231 and 232: MODELLI DI TIPO STAZIONARIO f z è
- Page 233 and 234: MODELLI DI TIPO STAZIONARIO di tipo
- Page 235 and 236: MODELLI DI TIPO STAZIONARIO nio di
- Page 237 and 238: MODELLI DI TIPO STAZIONARIO Proprio
- Page 239 and 240: MODELLI DI TIPO STAZIONARIO Questo
- Page 241 and 242: MODELLI DI TIPO STAZIONARIO ⇒ Rel
- Page 243 and 244: MODELLI DI TIPO STAZIONARIO Tab.4.4
- Page 245 and 246: MODELLI DI TIPO STAZIONARIO rise).
- Page 247 and 248: MODELLI DI TIPO STAZIONARIO dove U
- Page 249 and 250: MODELLI DI TIPO STAZIONARIO si ha:
- Page 251 and 252: MODELLI DI TIPO STAZIONARIO Terzo c
- Page 253 and 254: MODELLI DI TIPO STAZIONARIO più di
TEORIA DI BASE DELLA DISPERSIONE DEGLI INQUINANTI IN ATMOSFERA<br />
non consentono una descrizione diretta delle reazioni chimiche non lineari spesso<br />
molto importanti nel<strong>la</strong> normale dinamica del PBL.<br />
Come si vedrà nel seguito, i due approcci comportano differenti tipi di re<strong>la</strong>zioni<br />
matematiche per descrivere <strong>la</strong> concentrazione <strong>degli</strong> <strong>inquinanti</strong>, espressioni<br />
matematiche differenti che però, in ultima analisi, possono essere messe<br />
in re<strong>la</strong>zione le une con le altre. Ciascun approccio è una valida descrizione<br />
del<strong>la</strong> <strong>dispersione</strong> turbolenta: <strong>la</strong> scelta di quale approccio adottare in una data<br />
situazione dovrà essere deciso caso per caso, nelle situazioni reali. Il questo<br />
capitolo vengono presentati, in termini abbastanza generali, i due differenti<br />
approcci, riferendoci al<strong>la</strong> celebre lezione di Lamb (1980) e a quanto presentato<br />
in Seinfeld e Pandis (1998). Nei capitoli successivi introdurremo poi differenti<br />
famiglie di modelli che possono essere applicati sul<strong>la</strong> pratica e che<br />
comunque derivano da uno dei due approcci.<br />
3.1 I FONDAMENTI DELLA TEORIA EULERIANA.<br />
3.1.1 La formu<strong>la</strong>zione istantanea di base.<br />
L’ipotesi al<strong>la</strong> base dell’approccio euleriano è <strong>la</strong> sicurezza che le leggi del<strong>la</strong> fluidodinamica<br />
siano in grado di descrivere in modo completo e perfetto (e quindi<br />
deterministico) il fenomeno del<strong>la</strong> <strong>dispersione</strong> <strong>degli</strong> <strong>inquinanti</strong> nel PBL.E’questo,<br />
quindi, un approccio completamente deterministico ed in ultima analisi, si basa<br />
sul<strong>la</strong> scrittura dell’equazione di conservazione delle sostanze <strong>inquinanti</strong> presenti nel PBL.<br />
A differenza di quanto abbiamo fatto nel Cap.2, in questo caso non è più possibile<br />
considerare <strong>la</strong> misce<strong>la</strong> gassosa che costituisce il PBL come un’entità unica<br />
(che abbiamo chiamato genericamente aria) , ma al contrario siamo costretti:<br />
• a descrivere individualmente le diverse sostanze <strong>inquinanti</strong> a cui siamo interessati;<br />
• a tener conto del fatto che alcune di esse non sono gas inerti, ma sostanze più<br />
o meno reattive chimicamente;<br />
• a tener conto che alcune possono presentare un decadimento radioattivo;<br />
• a tener conto che alcune sono particelle solide o aerosol.<br />
Si consideri inizialmente un solo inquinante di tipo gassoso. Se si concentra l’attenzione<br />
su un piccolo volume d’aria (una particel<strong>la</strong>, dunque) come mostrato in<br />
Fig.3.1 e si trascura <strong>la</strong> diffusività moleco<strong>la</strong>re, è possibile cercare di applicare ad<br />
esso <strong>la</strong> conservazione del<strong>la</strong> massa dell’inquinante considerato. In partico<strong>la</strong>re, si<br />
consideri il flusso di inquinante che transita nel volume lungo <strong>la</strong> direzione x in<br />
un intervallo di tempo δt.Tale flusso è pari a:<br />
dove u 1 è <strong>la</strong> componente u del<strong>la</strong> velocità dell’aria entrante nel volume attraverso<br />
<strong>la</strong> faccia di sinistra e u 2 è <strong>la</strong> velocità uscente dal<strong>la</strong> faccia di destra. Dato che δx<br />
è piccolo, è possibile realizzare <strong>la</strong> seguente espansione in serie di Taylor:<br />
Da queste considerazioni risulta che il flusso in transito è pari a:<br />
204