Con questo articolo si puo' calcolare il decadimento della temperatura nel tempo t, quando un corpo ad una certa temperatura iniziale, si raffredda in un ambiente a temperatura minore. Questa teoria e' applicabile a tutti quei dispositivi che si riscaldano e che per un certo tempo non possono essere toccati e le applicazioni vanno anche fino ai circuiti elettronici di potenza.
Si sa che la distribuzione delle temperature nella trasmissione per conduzione del calore, sia in regime stazionario e non, viene regolata dalla equazione di Fourier a tre dimensioni data da:
con T temperatura assoluta, k conducibilita' ed ,
con ρ densita' e c calore specifico a pressione costante del mezzo.
α e' detta diffusivita'.
Il termine qv e' relativo alla potenza di calore assorbita o persa dal sistema (temine di sorgente o pozzo).
Inoltre vale anche in maniera piu' semplificata la relazione di Fourier nella trasmissione di calore unidimensionale, che e' data, per unita' di superficie trasversale alla propagazione del calore, da:
con k conducibilita' termica.
Ma nell'aria e nel vuoto abbiamo anche due modalita' di trasmissione del calore: convezione ed irraggiamento. Infatti il calore si trasmette sia per i moti delle particelle dell'aria a temperature diverse, sia perche' in corpo caldo puo' assimilarsi per la legge di Stephen-Boltzmann ad una sorgente elettromagnetica. La somma di queste due ultime modalita' di trasmissione del calore e' detta adduzione.
Andiamo a esaminare il caso di raffreddamento di un corpo. Se ipotizziamo una lastra metallica quadrata di lato l e spessore d ad una certa temperatura T, si puo' arrivare a determinare la legge di raffreddamento nel tempo, considerando congiuntamente i fenomeni di conduzione attraverso la sezione della lastra, la convezione vicino alla lastra e l'irraggiamento provocato dalla temperatura della lastra stessa. Esaminiamo il caso di convezione. Infatti s'ipotizza che vi sia a contatto della superficie della piastra uno strato laminare di aria dove il calore si trasmette per convezione,la cui profondita' e' nell'ordine dei millimetri. La quantita' di calore per unita' di tempo che viene scambiata per convezione e' data da:
con hc costante di convezione e Tc e Ta temperature del corpo lastra e del mezzo al di la' della lastra, ovvero l'aria.
D'altro canto, nella lastra la quantita' elementare di calore che fluisce tra due sezioni vicine distanti in temperatura e' data da:
sulla superficie vale la condizione limite che le quantita' trasmesse per conduzione e per convezione debbano essere uguali per cui si arriva all'equazione differenziale seguente in T(t)
Considerando l'irraggiamento, l'equazione di trasmissione del calore e' della stessa forma di quella di convezione per cui, applicando anche per questo caso la condizione al contorno sulla superficie della lastra, uguagliando flusso di calore per conduzione e per irraggiamento si giunge all'equazione differenziale:
con hi costante d'irraggiamento.
Sommando i due casi di convezione ed irraggiamento, l'equazione differenziale diventa:
Manipolando ancora un po', si ottiene:
con ha = hc + hi
valendo l'ipotesi tacitamente ammessa, di linearizzazione per piccole variazioni di temperatura. Il risultato dell'equazione e' la seguente funzione di temperatura nel tempo, che da' la legge di raffreddamento da una temperatura Tc ad una Ta.
La quantita' e' la costante di tempo del sistema considerato, tempo in cui la temperatura scende da quella massima Tc a una inferiore di
Per il calcolo dei coefficienti hc e hi occorre riferirsi alle teorie della trasmissione di calore per conduzione ed irraggiamento.
Nell'ambito della convezione si ha che il coefficiente hc e' legato al numero di Nusselt, secondo la relazione:
dove Nu e' il numero di Nusselt, Lc la dimensione della zona di convezione, ovvero la larghezza dello strato laminare nell'aria nel nostro caso e k la conducibilita' dell'aria.
Il numero di Nusselt e' legato all'andamento del numero Ra di Rayleigh, che e' il prodotto tra il numero di Grashoff Gr e il numero di Prandtl Pr. Si ha che
e' il numero di Grashoff
e' il numero di Prandtl
e quindi
dove, riepilogando
g e' l'accelerazione della gravita'
β il coefficiente di dilatazione dei gas perfetti, pari a con T temperatura assoluta.
Lc dimensione caratteristica della convezione
ν viscosita' dinamica del fluido
α viscosita' dinamica del fluido
Se
si ha che Nu = 0.59Ra0.25 ovvero si e' in convezione naturale con regime laminare;
se si e' in convezione naturale turbolenta con
Nu = 0.13Ra0.33
Il nostro caso sara' il primo dei due.
Per l'irraggiamento hc deriva dalla legge della potenza di calore per unita' di superficie irradiata da un corso grigio che vale
Confrontando con l'equazione della quantita' di calore trasmessa per irraggiamento vista precedentemente, si ricava che
In questa formula σ0 e' il coefficiente di Stephen - Boltzmann, che vale , mentre
a e' il coefficiente di assorbimento del corpo irradiante, assimilato ad un corpo grigio.