Térmico flujo de calor en un disco isotrópico

Consideremos ahora el problema de calentamiento de una placa circular con las fuentes de calor distribuidos sobre la superficie. Vamos a prescribir el flujo de calor desde ambos lados. Alrededor de la placa a lo largo de los bordes vamos a prescribir la Temperatura constante.

Como ejemplo consideramos una placa delgada con la conductividad térmica K=75 W/(m • oK) [Hierro Dúctil (SN)]. El radio de la placa es de R = 100 mm, grosor D = 5 mm. El valor de la densidad del flujo de calor es F=60 W/m2.  temperatura de los bordes alrededor de la placa esT=293.15 oK. (ver la figura).

Los datos de entrada y solución analítica se pueden ver aquí.

Vamos a localizar los sensores de la temperatura como se muestra en la figura y hacer un complot para ellos en r =0; 20; 40; 60 mm. En los puntos dados compararemos la solución numérica. obtenido utilizando el análisis AutoFEM con la solución analítica.

El modelo de elementos finitos con condiciones de frontera aplicados

Después de llevar a cabo el cálculo 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

857

1648

Tabla 2. Resultado " temperatura "

Radios r, mm

Solución numérica.
temperatura T*, K

Solución analítica
temperatura T, К

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

0

293.9491

293.9500

3.06E-04

20

293.9165

293.9110

1.87E-03

40

293.8206

293.8220

4.76E-04

60

293.6603

293.6620

5.79E-04

 

Conclusiones:

 

El error relativo de la solución numérica. en comparación con la solución analítica no exceda de 0,002% para los elementos lineales.

El problema se resolvió casi exactamente con el mínimo coste computacional. Es importante tener en cuenta que en un punto r = 20 mm la solución con elementos lineales resulta ser más preciso. Esto fue realmente posible debido a las propiedades de convergencia de la solución en el método de elementos finitos (convergencia en un sentido de la norma integral y menor número de puntos / argumentos juega un papel aquí). Sin embargo, en general, es imposible predecir la apariencia y la ubicación de estos puntos.

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