Flujo de calor en una esfera
Consideremos una esfera hueca con el radio interior r1, radio r 2 externos, que tiene coeficiente constante de conductividad térmica λ. La superficie interna de la esfera se mantiene a temperatura T1. Un intercambio de calor con el medio ambiente que tiene la temperatura T2 se lleva a cabo en la superficie externa. La intensidad de la transferencia de calor convectivo se caracteriza por el coeficiente de transferencia de calor β.
El modelo de elementos finitos con las cargas térmicas aplicadas |
Solución analítica del problema tiene la forma:
Para el cálculo numérico considerar 1/8 parte de la esfera hueca (ver figura). En los bordes laterales especificamos condiciones de simetría (el flujo de calor a través de los bordes laterales es Igual a 0).
Usemos los siguientes datos: radio interno de la esfera r1 = 150 mm radio externo de la esfera r2 = 250 mm.Coeficiente de conductividad térmica λ del material de la esfera es igual a 47 W / m.K.
La temperatura T1 en la superficie interna de la esfera es 373.15 K (or 100 oC). La temperatura del entorno ambiental T2 es igual to 298.15 K (or 25 oC), coeficiente de transferencia de calor β es igual to 100 W / (m2 . oC) .
Después de llevar a cabo el cálculo con la ayuda de AutoFEM, se obtienen los siguientes resultados:
Tabla 1. Parámetros de malla de elementos finitos
Tipo de elementos finitos. |
Número de Nodos |
Número de elementos finitos. |
Tetraedro lineal |
4871 |
25915 |
Tabla 2. Resultado " temperatura " de r= (3 r1+r2) / 4 = 0.175 m
Solución numérica. |
Solución analítica |
Error δ = 100%* |T* - T| / |T| |
3.66000671E+002 |
3.66138033E+002 |
0.04 |
Tabla 3.Resultado " temperatura "de r = (r1+r2) / 2 = 0.200 m
Solución numérica. |
Solución analítica |
Error δ = 100%* |T* - T| / |T| |
3.60805206E+002 |
3.60879058E+002 |
0.02 |
Tabla 4. Resultado " temperatura " de r = (r1+ 3 r2) / 4 = 0.225 m
Solución numérica. |
Solución analítica |
Error δ = 100%* |T* - T| / |T| |
3.56699341E+002 |
3.56788743E+002 |
0.025 |
Tabla 5. Resultado " temperatura " de r = r2 = 0.250 m
Solución numérica. |
Solución analítica |
Error δ = 100%* |T* - T| / |T| |
3.53451569E+002 |
3.53516492E+002 |
0.02 |
Los resultados de pruebas numéricas dependen de la malla de elementos finitos y pueden diferir ligeramente de los que figuran en la tabla.
Lea más acerca de AutoFEM Térmico Analysis