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.
temperatura T*, К

Solución analítica
temperatura T, К

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.
temperatura T*, К

Solución analítica
temperatura T, К

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.
temperatura T*, К

Solución analítica
temperatura T, К

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.
temperatura T*, К

Solución analítica
temperatura T, К

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

autofem.com

 

Volver al índice