Flusso di calore in una parete cilindrica
Si consideri una parete cilindrica di lunghezza L (tubo) con raggio interno r1, raggio esterno r2, avente costante il coefficiente di conducibilità termica λ. La superficie interna del tubo è mantenuta a temperatura T1 da fonti di calore distribuite uniformemente qν. Il calore generato sulla parete interna viene dissipato nell'ambiente circostante attraverso la superficie esterna del tubo (si veda figura).
|
La soluzione generale di questo problema ha la forma:
T = C1 ln (r) + C2 - qν r2/ 4λ.
Le costanti C1 e C2 sono determinate da condizioni prescritte sulle superfici interne (r = r1) ed esterna (r = r2) del tubo:
T |r=r1 = T1 ,
λ (dT / dr) |r=r2 = q .
Quindi,
![]()
.
Usiamo i seguenti dati: r1 = 100 mm, r2 = 250 mm, L = 1000 mm, λ = 43 W/m.K.
L'energia Q delle fonti di calore all'interno del tubo è pari a 4500 W. Poiché le fonti di calore qν sono distribuite uniformemente su tutto il volume del tubo,
.
Il flusso di calore specifico sulla superficie esterna del tubo q = -15000 W/m2. La Temperatura T1 sulla superficie interna del tubo è pari a 373.15 K (or 100 oC), mentre quella sulla superficie esterna T2 è pari alla temperatura ambiente 298.15 K (o 25 oC).
|
|
Il modello ad elementi finiti con applicati carichi termici |
Dopo aver effettuato il calcolo con l'aiuto di AutoFEM, sono stati ottenuti questi risultati:
Tabella 1. I parametri della maglia degli elementi finiti
Tipo di elemento finito |
Numero di nodi |
Numero di elementi finiti |
tetraedro lineare |
3309 |
14699 |
Tabella 2. Risultato "Temperatura" at r = r2 = 0.250 m *
Soluzione numerica |
Soluzione analitica |
Errore δ = 100* | T* - T | / | T |, % |
3.03744202E+002 |
3.03081513E+002 |
0.21 |
* I risultati dei test numerici dipendono dalla maglia ad elementi finiti e possono differire leggermente da quelli indicati nella tabella.
c