Métodos Numéricos en La Conducción de Calor

Métodos numéricos en la
conducción de calor
Andrés Felipe Rodriguez Valencia M.Sc
Universidad Autónoma de Occidente

Facultad de Ingeniería
Contenido
 Métodos numéricos en la conducción de calor
 ¿Por qué los métodos numéricos?
 Razones para la búsqueda de métodos alternativos de resolución
 Formulación en diferencias finitas de ecuaciones diferenciales
 Condiciones de frontera
 Conducción bidimensional de calor en estado estacionario
Métodos numéricos en la conducción de calor
• Muchos problemas que se encuentran en la practica de ingeniería comprenden configuraciones geométricas
complicadas, con condiciones de frontera complejas o propiedades variables, y no se pueden resolver
analíticamente. En estos casos se pueden obtener soluciones suficientemente exactas por medio de computadoras
utilizando un método numérico.
•En este capítulo se usará principalmente el enfoque de balance de energía, que se basa en los conocidos balances
de energía en volúmenes de control y no en pesadas formulaciones matemáticas.
¿PORQUÉ LOS MÉTODOS NUMÉRICOS?
Para el caso de la esfera en la figura podemos encontrar
una solución analítica para determinar la temperatura y la
velocidad de transferencia de calor, aplicando las
ecuaciones diferenciales:
1   2 T  1   T  1   T   T
 kr   
 kr 
   ksen    gen  c
r r 
2
r  r sen      r sen   
2 2 2 2
z  t
Solución analítica
RAZONES PARA LA BUSQUEDA DE MÉTODOS
ALTERNATIVOS DE RESOLUCIÓN
1. LIMITACIONES:
Los métodos analíticos de soluciones se
limitan a problemas simplificados en
configuraciones geométricas simples.
(cilindros, paredes planas y esferas)
2- UNA MEJOR ELABORACIÓN DE MODELOS:
La solución numérica aproximada de un problema del mundo real puede ser mas precisa que la
solución exacta (analítica) de un modelo simplificado en exceso (con el fin de hacer que el modelo
matemático sea suficientemente simple como para justificar una solución analítica)
3-FLEXIBILIDAD:
Los problemas de ingeniería a menudo requieren estudios paramétricos extensos con el fin de
entender la influencia de algunas variables sobre la solución y así elegir el conjunto correcto de
variables y dar respuesta a algunas preguntas de “¿Qué sucede si..?. Se trata de un proceso
iterativo que es muy tedioso y demorado. Los métodos numéricos y las computadoras resultan
ideales para estos cálculos.
• 4-COMPLICACIONES:
Algunos problemas se pueden resolver
analíticamente, pero el procedimiento de solución
es tan complejo y las expresiones resultantes de la
solución son tan complicadas que no vale la pena
todo el esfuerzo.
• La gran mayoría de los problemas de conducción
de calor llevan a ecuaciones diferenciales
parciales cuyas soluciones requieren
conocimientos profundos en transformada de
Laplace y de Fourier, funciones de Bessel y
Legendre, series de Taylor, etc.
5-NATURALEZA HUMANA:
• En la actualidad prácticamente todas las oficinas de ingeniería están equipadas con
computadoras de alto poder, con refinados paquetes de software, con diagramas a todo
color, con estilos de presentación gráficos o tabular y con resultados tan exactos como
los analíticos.
Formulación en diferencias finitas de ecuaciones
diferenciales
La pared se divide en M secciones de espesor
igual en la dirección x (Δx=L/M), separados por
los planos que pasan por los M+1 puntos
0,1,2,…,m-1,m,,+1,…M llamados nodos o
puntos nodales.
La ecuación de conducción de calor comprende la
segunda derivada de la temperatura con respecto a
Esquema de los nodos y las
temperaturas nodales usados en el las variables espaciales, tales como d2T/dx2 y la
desarrollo de la formulación en formulación en diferencias finitas se basa en el
diferencias finitas de la
transferencia de calor en una pared reemplazo de las segundas derivadas por
plana diferencias apropiadas
diferenciales
diferenciales
Unidimensional
Ecuación diferencial
La ecuación diferencial es válida en todo punto

medio, en tanto que la ecuación en diferencias
finitas sólo es válida en puntos discretos (los
nodos)
diferenciales
Bidimensional
Malla de diferencias finitas para

la conducción bidimensional en
coordenadas rectangulares.
Conducción unidimensional de calor en estado
estacionario (Balance de energía)
• El método del balance de energía se basa en
la subdivisión del medio en un número
suficiente de elementos de volumen, y a
continuación aplicar un balance de energía
en cada elemento.
Balance de Energía
Puntos nodales y elementos de Energía generada

volumen para la formulación en
diferencias finitas de la conducción
unidimensional en una pared plana.
Conducción unidimensional de calor en estado
estacionario (Balance de energía)
En la formulación en diferencias
Donde
finitas se supone que la temperatura
varía linealmente entre los nodos
Condiciones de frontera
• Hasta ahora hemos visto que para la obtención de la ecuación de diferencias finitas en un
nodo se requiere de dos nodos mas a cada lado, pero que pasa si nos toca un nodo en la
frontera?
Debemos aplicar un balance de energía en los elementos de volumen de los nodos frontera
• Las condiciones de frontera mas comunes en la

practica son las de temperatura especifica, de flujo
especifico de calor, convección y radiación.
Balance de Energía
Flujo de calor constante en la superficie
Superficie aislada
Superficie sometida a convección

Esquema para la formulación en
diferencias finites del nodo
frontera de la izquierda de una Superficie sometida a radiación
pared plana.
Superficie sometida a convección y radiación
Superficie sometida a convección y radiación
Superficie sometida a convección, radiación y flujo

de calor constante
TRATAMIENTO DE LOS NODOS EN UNA FRONTERA AISLADA COMO
NODOS INTERIORES: EL CONCEPTO DE IMAGEN ESPECULAR
• Una manera de obtener la formulación en diferencias
finitas de un nodo sobre una frontera aislada es tratar
el aislamiento como flujo de calor cero, es decir no le
entra calor por la izquierda. (como se vio
anteriormente).
• Otra forma mas practica es tratar el nodo sobre una
frontera aislada como uno interior. Esto se realiza al
reemplazar el aislamiento sobre la frontera por un
ESPEJO y considerar la reflexión del medio como su
extensión. (ver figura). De esta manera, el siguiente
nodo al nodo frontera aparece a ambos lados de este
último debido a la simetría, al convertirlo en un nodo
interior.
Conducción bidimensional de calor en estado
estacionario
• Considere una región rectangular en la cual la conducción
de calor es significativa en x y y.
• Divida la región en una malla rectangular de distancias
iguales en Δx y Δy. (ver figura) y considere una profundidad
unitaria Δz.
• El objetivo es determinar las temperaturas en los nodos y
resulta conveniente enumerarlos y describir su posición por
los números, en lugar de las coordenadas reales. Ejemplo
(m,n) donde m=0,1,2,…M en la dirección x y n=0,1,2…N
en la dirección y
Red nodal para la formulación en diferencias
finitas de la conducción bidimensional, en
coordenadas rectangulares.
estacionario
• Considere un elemento de volumen Δχ x Δy x 1,
con centro en un nodo interior, con generación de
calor y conductividad térmica constante, como
muestra la figura.
Balance de Energía
Dividiendo entre Δx y Δy tenemos:

estacionario
• En el análisis con diferencias finitas por lo común se usa, por sencillez una MALLA
CUADRADA (excepto cuando las magnitudes de los gradientes de temperatura en las
direcciones x y y son muy diferentes) y, por tanto Δx y Δy se consideran iguales. Entonces
Δx=Δy=l y la relación antes dada se simplifica a
Nodos Frontera
• La formulación en diferencias finitas de un nodo frontera se obtiene al escribir un balance de

energía sobre su elemento de volumen.
Métodos Numéricos en la conducción de calor:
Ejemplos
Si se considera estado estacionario y se toma k= 45W/m*°C, encuentre las temperaturas
en los puntos indicados (ΔX= ΔY=2cm).
Ejemplos
Solución
Por simetría
Balance de energía en los nodos

Ejemplos
Los gases calientes de la combustión de un horno fluyen por una chimenea de concreto (k= 1.4w/m°C) de
sección transversal rectangular. La sección de flujo de la chimenea tiene 20 cm X 40 cm y el espesor de la
pared es de 10 cm. La temperatura promedio de los gases calientes en la chimenea es de 280°C y el
coeficiente de transferencia de calor por convección dentro de esta última es de hi= 75 W/m2°C). la
chimenea pierde calor por convección desde su superficie exterior hacia el aire ambiente que está a T0=
15°C, con un coeficiente de transferencia de calor por convección de h0= 18W/m2°C, y hacia el cielo por
radiación. La emisividad de la superficie exterior de la pared es de e= 0.9 y se estima que la temperatura
del cielo es de 250 K. con Dx=Dy= 10 cm calcule:
• Formula en diferencias finitas para la transferencia de calor bidimensional en estado estacionario.
• Determine las temperaturas en todos los puntos nodales.
• Evalúe la razón de transferencia de calor para una sección de 1 m de largo de chimenea.
Métodos Numéricos en la conducción de calor: Ejemplos
a)
Métodos Numéricos en la conducción de calor: Ejemplos
b)
c)
REFERENCIAS
[1] CENGEL, Yunus A. Transferencia de Calor. McGraw-Hill Higher Education, 2003. 792
p.

Métodos Numéricos en La Conducción de Calor

Cargado por

Copyright:

Formatos disponibles

Métodos Numéricos en La Conducción de Calor

Cargado por

Información del documento

Descripción original:

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Métodos Numéricos en La Conducción de Calor

Cargado por

Copyright:

Formatos disponibles

Métodos numéricos en la

Andrés Felipe Rodriguez Valencia M.Sc

Universidad Autónoma de Occidente

La ecuación diferencial es válida en todo punto

Malla de diferencias finitas para

Puntos nodales y elementos de Energía generada

• Las condiciones de frontera mas comunes en la

Flujo de calor constante en la superficie

Superficie sometida a convección

Superficie sometida a convección y radiación

Superficie sometida a convección, radiación y flujo

Dividiendo entre Δx y Δy tenemos:

• La formulación en diferencias finitas de un nodo frontera se obtiene al escribir un balance de

Balance de energía en los nodos

También podría gustarte