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