Umidità e flussi evaporativi
INTRODUZIONE
In questa esperienza abbiamo verificato la predizione teorica dell'andamento del flusso evaporativo in funzione del raffreddamento e del contenuto ponderale d'acqua e abbiamo verificato 1'andamento temporale della temperatura durante la fase di evaporazione.
Il flusso evaporativo (Fev) è una grandezza che dipende linearmente dalla variazione di peso del campione (DP), dal tempo (Dt) e dalla superficie interessata dal fenomeno di evaporazione (S), come è chiaramente esposto nell'equazione seguente:
(1)
Il contenuto ponderale d'acqua (W) è invece il rapporto tra la massa d'acqua presente nel campione in esame (ma) e la sua massa “secca” (ms). In simboli:
,
(2), (3)
dove, naturalmente, m = ma + ms è la massa totale del campione umido.
Considerando l'insieme degli scambi energetici tra un corpo e l'ambiente esterno, si arriva a una semplice equazione di bilancio:
(4)
dove
msup = densità superficiale di massa
c = calore specifico del campione umido
=
variazione della temperatura in funzione del tempo
e = emissività del campione
s
= costante di Stephan-Boltzman =
T = temperatura superficiale del campione
a = coefficiente di assorbimento del campione
eamb = emissività dell'ambiente
Tamb = temperatura dell'ambiente
h = coefficiente di scambio convettivo
=
coefficiente di conduzione =
,
dove k = conducibilità termica del corpo
Tint = temperatura del campione a profondità Dx
P = costante solare
Lev
= calore latente di evaporazione dell'acqua ≈ 2,4 ·
106
Per quello che riguarda in specifico la nostra esperienza in laboratorio, è stato possibile operare alcune semplificazioni, che hanno dato all'equazione una forma più facilmente indagabile.
Si può infatti
considerare P = 0, si può assumere che Tamb = Tint ,
cioè che il corpo abbia la stessa temperatura dell'ambiente esterno,
e che
Così facendo si ottiene l'equazione su cui in realtà abbiamo lavorato:
(5)
STRUMENTI UTILIZZATI
TERMOCAMERA
BILANCIA DI PRECISIONE
FITOTRONE: questo strumento, nato in origine per la conservazione di organismi vegetali, consente di mantenere all'interno della sua camera una temperatura e un grado di umidità pressoché costanti, mediante un opportuno pannello di regolazione.
SONDA PER MISURARE L'UMIDITA': essa è stata posta all'interno del fitotrone per permettere una precisa lettura del livello di umidità nella camera, poiché, nonostante il valore sia impostato con precisione, esso subisce delle leggere variazioni, dovute principalmente al termostato dell'apparecchio.
CRONOMETRO
PARAFILM “M” LABORATORY FILM: questa pellicola trasparente, essendo completamente impermeabile e facilmente deformabile, ci ha permesso di isolare i campioni da analizzare, sigillandone tutte le facce non interessanti, e lasciando che il fenomeno di evaporazione avvenisse solo sulla superficie scelta.
I campioni su cui abbiamo lavorato sono quattro carote di mattone da muratura. Pur essendo completamente uguali per dimensioni (5 cm il diametro e 2 cm l'altezza) e per materiale, possiedono caratteristiche diverse a seconda del diverso tipo di trattamento che hanno subito, e il grado di traspirazione che li caratterizza è radicalmente differente.
Vengono distinti in A, Ca, B, Cb.
Oltre ai quattro campioni analizzati abbiamo utilizzato anche due ulteriori carote di mattone, uguali in tutto e per tutto alle altre, ma completamente secche. Poste all'interno del fitotrone ci hanno fornito due standard di riferimento per monitorare direttamente le condizioni climatiche reali.
PROCEDURA SPERIMENTALE
PARTE 1: Flusso e raffreddamento
La prima fase dell'esperienza consisteva nel verificare l'andamento lineare del flusso evapoorativo in funzione del raffreddamento (inteso come differenza tra la temperatura all'inizio del processo e quella alla fine).
Abbiamo cominciato impostando i corretti valori di temperatura e umidità sul fitotrone, regolando la prima su un valore stabile di 24,5°C, e variando di volta in volta la seconda da un minimo di 30% a un massimo di 80%.
Per mantenere costanti i valori impostati si poteva utilizzare un dispositivo di ventilazione interno al fitotrone, la cui regolazione però ci ha fornito qualche problema. Infatti un livello di ventilazione alto, se da un lato permetteva di mantenere costanti temperatura e grado di umidità, dall'altro ostacolava la lettura del peso dei campioni sulla bilancia facendola oscillare, e soprattutto aumentava sensibilmente il processo di evaporazione dei campioni, falsando notevolmente i risultati.
Siamo riusciti a raggiungere la condizione di lavoro migliore regolando la ventilazione a metà circa dell'intensità, ponendo degli schermi di cartone tra la bilancia e la grata di fuoriuscita del getto di aria, in modo che questo non investisse direttamente i campioni, e mantenendo gli sportelli laterali della bilancia chiusi durante le pesate, lasciando aperto solo quello superiore per garantire comunque un certo grado di ricircolo di aria.
Dopo aver posto la termocamera sopra al fitotrone, in modo che dalla finestra sulla faccia superiore riuscisse ad inquadrare i campioni all'interno, dopo aver messo in bolla e tarato la bilancia, e dopo aver impostato sulla camera il corretto valore di emissività (0,81 per il mattone), abbiamo incominciato col prendere le misure.
Abbiamo misurato il peso dei campioni, secchi e senza pellicola isolante, e, dopo averli lasciati in acqua per qualche giorno, li abbiamo posti all'interno del fitotrone, dando inizio al processo di evaporazione. A turno ciascuno dei quattro mattoncini è stato posto sulla bilancia per la prima pesata e, nel momento in cui il valore del peso si è stabilizzato, è stato dato il via al cronometro. Dopo circa tre minuti si sono ripetute le stesse misure di peso, riportando il tempo trascorso, e ottenendo così due valori per la massa (P1, P2) e due valori per il tempo (t1, t2). Quindi si è atteso che i campioni raggiungessero la temperatura di equilibrio, e quando ciò è avvenuto si sono ripetute per verifica le stesse misure di prima, che ci hanno permesso in sede di analisi dei dati di considerare come valore effettivamente misurato la media matematica dei due singoli valori, riducendo così l'errore sperimentale.
Questo procedimento e' stato ripetuto per ciascuno dei diversi valori di umidità impostati sul fitotrone, e l'intera gamma di misure è stata svolta nuovamente il giorno dopo, quando si è passati da un contenuto ponderale d'acqua di circa il 20% a uno del 10%.
ANALISI DEI DATI
Riportiamo di seguito, dopo una tabella riassuntiva dei valori misurati e dei calcoli effettuati, i grafici relativi ai risultati ottenuti. E' da notare che il primo grafico si riferisce alle misure effettuate con quattro diversi valori di umidità, mentre nel secondo grafico abbiamo deciso di infittire i dati prendendo altre misure con altri due distinti valori.
A ciascun valore di umidità corrisponde una diversa colorazione dei dati ottenuti, così da rendere facilmente riconoscibili i relativi punti sui grafici.
FLUSSO IN FUNZIONE DEL RAFFREDDAMENTO
giorno 1
UMIDITA' |
CAMPIONE |
tempo 1 |
massa 1 |
tempo 2 |
massa 2 |
FLUSSO |
temp.secco |
temp.eq. |
raffreddamento |
FLUSSO MEDIO |
|
|
|
|
|
|
|
|
|
|
|
80 % |
A |
55 |
78,0539 |
258 |
78,0327 |
5,33E-05 |
23,09 |
20,36 |
2,73 |
4,6188E-05 |
|
B |
37 |
76,3086 |
238 |
76,2878 |
5,28E-05 |
23,09 |
20,33 |
2,76 |
4,74752E-05 |
|
Ca |
88 |
76,7308 |
276 |
76,7112 |
5,32E-05 |
23,09 |
20,39 |
2,7 |
4,7689E-05 |
|
Cb |
0 |
78,056 |
295 |
78,0204 |
6,16E-05 |
23,09 |
20,44 |
2,65 |
5,0694E-05 |
|
|
|
|
|
|
|
|
|
|
|
|
A |
1786 |
77,8998 |
2017 |
77,8821 |
3,91E-05 |
|
|
|
|
|
B |
1730 |
76,1534 |
2104 |
76,1225 |
4,22E-05 |
|
|
|
|
|
Ca |
1750 |
76,578 |
2033 |
76,5546 |
4,22E-05 |
|
|
|
|
|
Cb |
1700 |
77,8962 |
2078 |
77,8667 |
3,98E-05 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
60 % |
A |
24 |
77,782 |
201 |
77,7568 |
7,26E-05 |
22,39 |
18,46 |
3,93 |
7,42802E-05 |
|
B |
0 |
76,0377 |
222 |
76,0053 |
7,45E-05 |
22,39 |
18,56 |
3,83 |
7,28146E-05 |
|
Ca |
60 |
76,4486 |
272 |
76,4155 |
7,97E-05 |
22,39 |
18,44 |
3,95 |
7,79161E-05 |
|
Cb |
43 |
77,7725 |
250 |
77,7434 |
7,17E-05 |
22,39 |
18,63 |
3,76 |
7,19358E-05 |
|
|
|
|
|
|
|
|
|
|
|
|
A |
1626 |
77,5576 |
1919 |
77,514 |
7,59E-05 |
|
|
|
|
|
B |
1447 |
75,8315 |
1837 |
75,7771 |
7,12E-05 |
|
|
|
|
|
Ca |
1467 |
76,2427 |
1895 |
76,1788 |
7,62E-05 |
|
|
|
|
|
Cb |
1515 |
77,5615 |
1870 |
77,5113 |
7,21E-05 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
45 % |
A |
0 |
77,4439 |
295 |
77,3896 |
9,39E-05 |
22,82 |
16,12 |
6,7 |
8,97677E-05 |
|
B |
30 |
75,6915 |
322 |
75,6361 |
9,68E-05 |
22,82 |
16,03 |
6,79 |
9,42048E-05 |
|
Ca |
60 |
76,0928 |
350 |
76,0366 |
9,89E-05 |
22,82 |
15,83 |
6,99 |
9,5729E-05 |
|
Cb |
92 |
77,4163 |
375 |
77,3623 |
9,74E-05 |
22,82 |
16,24 |
6,58 |
9,54934E-05 |
|
|
|
|
|
|
|
|
|
|
|
|
A |
2994 |
76,8832 |
3398 |
76,8154 |
8,56E-05 |
|
|
|
|
|
B |
3020 |
75,1276 |
3426 |
75,0547 |
9,16E-05 |
|
|
|
|
|
Ca |
3045 |
75,498 |
3482 |
75,4187 |
9,26E-05 |
|
|
|
|
|
Cb |
3067 |
76,8524 |
3510 |
76,7711 |
9,36E-05 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
30 % |
A |
0 |
76,74 |
317 |
76,6676 |
0,000117 |
21,68 |
14,8 |
6,88 |
0,000116488 |
|
B |
44 |
74,9778 |
357 |
74,907 |
0,000115 |
21,68 |
14,59 |
7,09 |
0,000112729 |
|
Ca |
74 |
75,3381 |
394 |
75,2621 |
0,000121 |
21,68 |
14,56 |
7,12 |
0,000117089 |
|
Cb |
98 |
76,6966 |
427 |
76,6184 |
0,000121 |
21,68 |
14,72 |
6,96 |
0,000113413 |
|
|
|
|
|
|
|
|
|
|
|
|
A |
1432 |
76,4101 |
1779 |
76,3309 |
0,000116 |
|
|
|
|
|
B |
1460 |
74,6609 |
1804 |
74,5867 |
0,00011 |
|
|
|
|
|
Ca |
1484 |
75,0085 |
1833 |
74,9312 |
0,000113 |
|
|
|
|
|
Cb |
1507 |
76,3767 |
1870 |
76,3016 |
0,000106 |
|
|
|
|
giorno 2
UMIDITA' |
CAMPIONE |
tempo 1 |
massa 1 |
tempo 2 |
massa 2 |
FLUSSO |
temp.secco |
temp.eq. |
raffreddamento |
FLUSSO MEDIO |
|
|
|
|
|
|
|
|
|
|
|
30 % |
A |
50 |
71,1345 |
483 |
71,0715 |
7,42E-05 |
22,72 |
18,37 |
4,35 |
6,70833E-05 |
|
B |
78 |
70,4334 |
511 |
70,403 |
3,58E-05 |
22,72 |
20,22 |
2,5 |
3,34347E-05 |
|
Ca |
24 |
70,299 |
536 |
70,2528 |
4,6E-05 |
22,72 |
20,17 |
2,55 |
4,26782E-05 |
|
Cb |
0 |
71,3768 |
559 |
71,3164 |
5,51E-05 |
22,72 |
19,65 |
3,07 |
5,15641E-05 |
|
|
|
|
|
|
|
|
|
|
|
|
A |
1345 |
70,956 |
1677 |
70,917 |
5,99E-05 |
|
|
|
|
|
B |
1378 |
70,346 |
1705 |
70,3261 |
3,1E-05 |
|
|
|
|
|
Ca |
1408 |
70,1798 |
1735 |
70,1546 |
3,93E-05 |
|
|
|
|
|
Cb |
1441 |
71,2276 |
1762 |
71,1974 |
4,8E-05 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
60 % |
A |
0 |
70,8796 |
256 |
70,8634 |
3,23E-05 |
22,9 |
19,98 |
2,92 |
3,41202E-05 |
|
B |
22 |
70,3092 |
283 |
70,2996 |
1,88E-05 |
22,9 |
21,19 |
1,71 |
1,94112E-05 |
|
Ca |
46 |
70,132 |
328 |
70,1183 |
2,48E-05 |
22,9 |
20,95 |
1,95 |
2,50131E-05 |
|
Cb |
67 |
71,172 |
361 |
71,1547 |
3E-05 |
22,9 |
20,71 |
2,19 |
2,9732E-05 |
|
|
|
|
|
|
|
|
|
|
|
|
A |
1201 |
70,7964 |
1499 |
70,7754 |
3,6E-05 |
|
|
|
|
|
B |
1232 |
70,2641 |
1522 |
70,2527 |
2,01E-05 |
|
|
|
|
|
Ca |
1263 |
70,0719 |
1546 |
70,0579 |
2,52E-05 |
|
|
|
|
Legenda:
A,B = carote di mattone non trattate
Ca,Cb = carote di mattone trattate
CONCLUSIONI
Dai grafici ottenuti si riconosce subito l'andamento lineare, come ragionevolmente previsto, anche se, come sottolinea il relativo valore dell'indice di correlazione lineare, l'interpolazione dei dati non è particolarmente buona. Infatti al primo set di misure corrisponde un valore di R2 di 0,8967 che non garantisce una buona correlazione lineare; mentre, relativamente al secondo grafico, ottenuto infittendo le misure, si ottiene un valore di 0,9153 che, rimanendo comunque al di sotto delle aspettative, costituisce comunque un miglioramento.
Anche i valori di intercetta con l'asse delle ordinate costituiscono una prova dell'imprecisione dei risultati. Infatti teoricamente a un raffreddamento uguale a zero dovrebbe corrispondere un flusso uguale a zero, mentre i valori da noi ottenuti, pur essendo molto bassi (2E-5 per il primo grafico e 4E-6 per il secondo) non sono trascurabili rispetto agli ordini di grandezza del flusso.
PARTE 2: Flusso e Contenuto ponderale d'acqua
L'obiettivo di questa seconda parte dell'esperienza è stato il mostrare la correlazione esistente tra flusso evaporativo e contenuto ponderale d'acqua.
Abbiamo per prima cosa collegato la bilancia ad un PC, in modo che ogni 15 secondi venisse prelevato il peso del campione in esame e automaticamente registrato su un foglio di lavoro EXCEL.
Abbiamo posto i due standard secchi all'interno della bilancia, sempre per avere un riscontro diretto delle condizioni climatiche all'interno, e abbiamo fatto si che la termocamera sulla sommità del fitotrone inquadrasse direttamente il piatto della bilancia.
Dopo aver impostato i corretti valori di temperatura, sensibilità, emissività e media temporale sulla termocamera, abbiamo dato il via alle misure, cominciando a prelevare immagini e pesi nello stesso istante. Abbiamo programmato la termocamera in modo che registrasse la prima ventina di immagini ogni 15 secondi (perché la prima fase dell'evaporazione è la più importante), poi una decina ogni 45 secondi, poi ogni minuto e infine ogni 2 minuti.
ANALISI DEI DATI
E' stato possibile ottenere il grafico relativo all'andamento del flusso evaporativo in funzione del contenuto ponderale d'acqua poiché eravamo a conoscenza sia dell'andamento del peso del campione in funzione del tempo, ottenuto direttamente in uscita dalla bilancia (grafici 3,4,5,6), sia del flusso in funzione del tempo, calcolato come differenza di peso diviso per il tempo e per l'area interessata, come riferito nella formula “1” (grafici 7,8,9,10).
Riportiamo ora i grafici relativi:
Come si può notare dai grafici l'andamento della massa del campione in funzione del tempo è pressoché lineare, se si eccettuano degli scostamenti dovuti probabilmente alle non perfettamente stabili condizioni interne al fitotrone.
Inoltre è interessante osservare che, nonostante qualche piccola oscillazione, il flusso evaporativo è abbastanza costante nel tempo (ricordiamo infatti che le variazioni di flusso interessano valori dell'ordine di 10-5 Kg/m2s). Per ovviare a questi scostamenti abbiamo pensato di mediare i valori di flusso attribuendo ad ogni singolo dato la media con i successivi nove; questo metodo permette di regolarizzare gli andamenti troppo rumorosi.
Riportiamo ora i grafici relativi all'andamento del flusso evaporativo in funzione del contenuto ponderale d'acqua.
Dall'osservazione dei grafici non si riesce a riconoscere nessuna relazione tra le due grandezze. Probabilmente un'analisi su tempi più lunghi e su un intervallo di contenuti ponderali più ampio permetterebbe un'indagine più approfondita. Infatti dalla rapida osservazione di un grafico mostratoci dal dottor Ludwig, ottenuto da una misura su una scala di tempi più ampia, abbiamo potuto riconoscere a grandi linee gli andamenti nelle zone corrispondenti ai nostri intervalli di misura, riscontrando un minimo accordo.
Inoltre abbiamo notato che localmente il flusso subisce variazioni significative dovute ad una forte instabilità del fitotrone nel mantenere costanti l'umidità, la temperatura e la ventilazione.
Parte terza: temperatura in funzione del tempo
Come ultima misura abbiamo studiato l'andamento della temperatura dei campioni durante la fase di evaporazione.
Abbiamo effettuato queste misure in contemporanea con quelle di peso, prelevando con la termocamera sia la temperatura del campione in esame sia quella del campione “secco” di riferimento. La temperatura di quest'ultimo è servita per tener conto nell'analisi delle oscillazioni delle condizioni ambientali all'interno del fitotrone. Infatti abbiamo notato che le variazioni anomale di temperatura nei due tracciati erano simili e contemporanee; questa è stata la prova che tali oscillazioni erano dovute al termostato del fitotrone. Abbiamo preso come temperatura di riferimento del campione secco la minima registrata e l'abbiamo sottratta alle temperature registrate istante per istante. I valori così ottenuti, indici delle oscillazioni, sono stati successivamente sottratti alle temperature dei campioni in esame. Questo è stato il miglor sistema per ovviare alla non perfetta stabilità delle condizioni interne del fitotrone, altrimenti non direttamente risolvibili.
La termocamera ci ha fornito delle immagini a falsi colori; successivamente con vari software abbiamo convertito le immagini in tonalità di grigio ed infine in dati di temperatura.
Ecco di seguito i grafici ottenuti.
Per mancanza di tempo non abbiamo approfondito il nostro studio con interpolazioni numeriche, però risulta abbastanza evidente dai grafici che l'andamento della temperatura è un esponenziale negativo. Questo è in accordo col modello del bilancio energetico nell'ipotesi che il raffreddamento per irraggiamento sia costante, per il quale la soluzione dell'equazione (4) è del tipo:
VERIFICA DEI VALORI DI FLUSSO OTTENUTI
Come ultima parte dell'esperienza abbiamo verificato la validità dei valori di flusso misurati, utilizzando la formula seguente che fornisce il flusso evaporativo in funzione della temperatura.
Questa non è altro che l'equazione di bilancio (formula 4) scritta in funzione del flusso nell'ipotesi che il campione abbia raggiunto la temperatura di equilibrio (dT/dt = 0), e che non vi siano fenomeni di natura convettiva (h=0).
Riportiamo di seguito la tabella recante, per ciascun campione e ciascuna condizione di lavoro, i due diversi valori di flusso
Campioni |
T fitotrone |
Umidità |
T equilibrio |
Flusso pesato |
Flusso calcolato |
---|---|---|---|---|---|
A |
23,09 |
80 % |
20,36 |
4,619E-05 |
1,22551E-05 |
B |
|
|
20,33 |
4,748E-05 |
1,23898E-05 |
Ca |
|
|
20,39 |
4,769E-05 |
1,21203E-05 |
Cb |
|
|
20,44 |
5,069E-05 |
1,18958E-05 |
A |
22,39 |
60 % |
18,46 |
7,428E-05 |
1,76437E-05 |
B |
|
|
18,56 |
7,281E-05 |
1,71947E-05 |
Ca |
|
|
18,44 |
7,792E-05 |
1,77335E-05 |
Cb |
|
|
18,63 |
7,194E-05 |
1,68804E-05 |
A |
22,82 |
45 % |
16,12 |
8,977E-05 |
3,00814E-05 |
B |
|
|
16,03 |
9,420E-05 |
3,04855E-05 |
Ca |
|
|
15,83 |
9,573E-05 |
3,13835E-05 |
Cb |
|
|
16,24 |
9,549E-05 |
2,95426E-05 |
A |
21,68 |
30 % |
14,8 |
1,165E-04 |
3,08897E-05 |
B |
|
|
14,59 |
1,127E-04 |
3,18326E-05 |
Ca |
|
|
14,56 |
1,171E-04 |
3,19673E-05 |
Cb |
|
|
14,72 |
1,134E-04 |
3,12489E-05 |
Come si può notare i valori di flusso riportati nelle due colonne, pur non essendo completamente in accordo fra loro, mantengono comunque lo stesso ordine di grandezza. La non perfetta accordanza tra i risultati è sicuramente imputabile alle tutt'altro che stabili condizioni ambientali all'interno del fitotrone.