Flux de chaleur dans une sphère

Examiner une sphère creuse avec le rayon intérieur r1, rayon extérieur r2, ayant un coefficient de conductivité thermique constant λ. Surface interne de la sphère est maintenue à la température T1. Un échange de chaleur avec l'environnement ayant la température T2 e déroule sur la surface externe. L'intensité du transfert de chaleur par convection est caractérisé par le coefficient de transfert de chaleur β.

Le modèle éléments finis appliquées aux charges thermiques

Solution analytique du problème est de la forme:

Pour le calcul numérique considérer seulement 1/8 de la sphère creuse (voir figure). Sur les bords latéraux on précise des conditions de symétrie (le flux de chaleur à travers les bords latéraux est égale à 0). Nous laisser utiliser les données suivantes: rayon intérieur de la sphère r1 = 150 mm, rayon extérieur de la sphère r2 = 250 mm. Coefficient de conductivité thermique λ du matériau de la sphère est égale à 47 W / m.K.
La température T1  à la surface interne de la sphère est 373.15 K (ou 100 oC). La température de l'environnement ambiant T2 est égal à  298.15 K (ou 25 oC), le coefficient de transfert de chaleur β est égale à 100 W / (m2 . oC) .

Après avoir effectué le calcul à l'aide de AutoFEM, les résultats suivants sont obtenus:

Tableau 1. Paramètres de maillage éléments finis

Type d'élément fini

Nombre de nœuds

Nombre d'éléments finis

tétraèdre linéaire

1660

7365

Tableau 2. Résultat "Temperature" at r= (3 r1+r2) / 4 = 0.175 m

Solution numérique
Temperature T*, К

Solution analytique
Temperature T, К

Erreur δ = 100%* |T* - T| / |T|

3.66257111E+002

3.66138033E+002

0.03

Tableau 3. Résultat «Température» à r = (r1+r2) / 2 = 0.200 m

Solution numérique
Temperature T*, К

Solution analytique
Temperature T, К

Erreur δ = 100%* |T* - T| / |T|

3.60996582E+002

3.60879058E+002

0.03

Table 4. Résultat «Température» à r = (r1+ 3 r2) / 4 = 0.225 m

Solution numérique
Temperature T*, К

Solution analytique
Temperature T, К

Erreur δ = 100%* |T* - T| / |T|

3.56892792E+002

3.56788743E+002

0.03

Table 5. Résultat «Température» à r = r2 = 0.250 m

Solution numérique
Temperature T*, К

Solution analytique
Temperature T, К

Erreur δ = 100%* |T* - T| / |T|

3.53626129E+002

3.53516492E+002

0.03

 

 

*Les résultats des tests numériques dépendent du maillage éléments finis et peuvent différer légèrement de celles indiquées dans le tableau.

 

En savoir plus sur AutoFEM Thermal Analysis

 

Retour au sommaire

autofem.com