Flusso di calore in una parete sferica
Consideriamo un guscio sferico con il raggio interno r1, raggio esterno r2, avente costante il coefficiente conducibilità termica λ, che dipende dal materiale. La superficie interna del guscio è mantenuta alla temperatura T1, mentre la superficie esterna è mantenuta alla temperatura T2. Tra la superficie interna del guscio e quelle esterna vi è trasferimento di calore secondo convezione caratterizzato dal coefficiente di trasferimento di calore β.
Il modello ad elementi finiti con applicati carichi termici |
Soluzione analitica del problema ha la forma:
Per il calcolo numerico consideriamo 1/ 8 parte del guscio sferico (si veda figura). Sui bordi laterali sono definiti le condizioni al contorno (il flusso di calore attraverso i margini laterali è uguale a 0).
Usiamo i seguenti dati: r1 = 150 mm, r2 = 250 mm, λ = 47 W/m.K, T1 = 373.15 K (o 100 oC),T2 = 298.15 K (o 25 oC), β = 100 W / (m2 . oC) .
Dopo aver effettuato il calcolo con l'aiuto di AutoFEM, sono stati ottenuti questi risultati:
Tabella 1. I parametri della maglia ad elementi finiti
Tipo di elemento finito |
Numero di nodi |
Numero di elementi finiti |
tetraedro lineare |
1660 |
7365 |
Tabella 2. Risultato "Temperatura" at r= (3 r1+r2) / 4 = 0.175 m
Soluzione numerica |
Soluzione analitica |
Errore δ = 100%* |T* - T| / |T| |
3.66257111E+002 |
3.66138033E+002 |
0.03 |
Tabella 3.Risultato "Temperatura" at r = (r1+r2) / 2 = 0.200 m
Soluzione numerica |
Soluzione analitica |
Errore δ = 100%* |T* - T| / |T| |
3.60996582E+002 |
3.60879058E+002 |
0.03 |
Tabella 4. Risultato "Temperatura" at r = (r1+ 3 r2) / 4 = 0.225 m
Soluzione numerica |
Soluzione analitica |
Errore δ = 100%* |T* - T| / |T| |
3.56892792E+002 |
3.56788743E+002 |
0.03 |
Tabella 5. Risultato "Temperatura" at r = r2 = 0.250 m
Soluzione numerica |
Soluzione analitica |
Errore δ = 100%* |T* - T| / |T| |
3.53626129E+002 |
3.53516492E+002 |
0.03 |
* I risultati dei test numerici dipendono dalla maglia ad elementi finiti e possono differire leggermente da quelli indicati nella tabella.
Approfondire riguardo ad AutoFEM Analisi Termica