CARTAS DE HEISLER Y DIFERENCIAS FINITAS

CARTAS DE HEISLER Y DIFERENCIAS FINITAS

Universidad de Córdoba – Sede Berástegui.
Luis carlos Altamiranda Sanchez
Carlos Albeiro Salazar Vergara
Gustavo Andres Ramirez Oviedo
Rafael David Rangel Terraza
Roberto Jose Yances Acevedo

OBJETIVOS

 Conocer el alcance y las limitaciones que ofrecen las cartas de Heisler

en la solución de problemas de transferencia de calor por conducción en
estado transitorio.

RESUMEN

En ésta práctica se observó el proceso de transferencia de calor por

convección natural, utilizando una probeta de arena con una temperatura inicial
uniforme que es expuesta a un medio con una temperatura diferente,
desarrollándose un gradiente de temperatura a lo largo de la transferencia de
calor que evolucionó con el tiempo. Para este experimento se utilizó una
probeta de arena, tres termómetros digitales, una mufla, un cronómetro y un pie
de rey.
Con los datos obtenidos se calculó el número de Fourier, el número
adimensional de temperatura y se graficó cada una de estas entre sí y entre la
temperatura con el fin de ver su comportamiento en comparación con las
cartas de Heisler, el cual dio muy parecidas en tan sólo un tramo de la gráfica
de éstas.
Palabras Claves: Cartas de Heisler, Diferencias Finitas, Propiedades Térmicas
de la Arena, Convección Natural.

ABSTRACT
In practice it was observed the process of heat transfer by natural convection,
using a sand specimen with a uniform initial temperature is exposed to a
medium with a different temperature to develop a temperature gradient along
the heat transfer evolved over time. For this experiment we used a sand
specimen, three digital thermometers, a furnace, a stopwatch and a caliper.

With the data we calculated the number of Fourier, the dimensionless number is
plotted temperature and each of these among themselves and between the
temperature in order to see their behavior as compared with the Heisler charts,
which gave very similar only a portion of the graph of these.

Keywords: Letters of Heisler, Finite Differences, Thermal Properties of the

Sand, Natural Convection.

LABORATORIO DE TRANSFERENCIA DE CALORPágina 1

 CARTAS DE HEISLER Y DIFERENCIAS FINITAS

INTRODUCCIÓN

Donde h es el coeficiente de
La convección es una de las tres convección (ó coeficiente de
formas de transferencia de calor y película), As es el área del cuerpo
se caracteriza porque se produce en contacto con el fluido, Ts es la
por intermedio de un fluido (aire, temperatura en la superficie del
agua) que transporta el calor entre
cuerpo y es la temperatura del
zonas con diferentes temperaturas.
fluido lejos del cuerpo. (Wikipedia,
La convección se produce
2006).
únicamente por medio de materiales
fluidos. Éstos, al calentarse,
En 1947, M. P. Heisler presentó los
aumentan de volumen y, por lo
diagramas de temperatura
tanto, disminuyen su densidad y
transitoria para una pared plana
ascienden desplazando el fluido que
grande, un cilindro largo y una
se encuentra en la parte superior y
esfera y se llama diagrama de
que está a menor temperatura. Lo
Heisler. H. Grober las complementó
que se llama convección en sí, es el
en 1961 con los diagramas de
transporte de calor por medio de las
transferencia de calor en régimen
corrientes ascendente y
transitorio. Existen tres diagramas
descendente del fluido.
asociados con cada figura
geométrica: la primera es para
La transferencia de calor implica el
determinar la temperatura To en el
transporte de calor en un volumen y
centro de la configuración en un
la mezcla de elementos
instante dado t. la segunda es para
macroscópicos de porciones
determinar la temperatura en otros
calientes y frías de un gas o un
lugares en el mismo instante, en
líquido. Se incluye también el
términos de To. La tercera es para
intercambio de energía entre una
determinar la cantidad total de
superficie sólida y un fluido o por
transferencia de calor hasta el
medio de una bomba, un ventilador
instante t. éstas gráficas son validas
u otro dispositivo mecánico
para τ > 0.2.
(convección mecánica o asistida).
Nótese el caso 1/Bi = k/hL = 0
En la transferencia de calor libre o
corresponde a h → ∞, lo cual
natural en la cual un fluido es más
corresponde al caso de temperatura
caliente o más frío y en contacto
superficial T∞ especificada. Es
con una superficie sólida, causa una
decir, el caso en el que las
circulación debido a las diferencias
superficies se llevan súbitamente a
de densidades que resultan del
la temperatura T∞ en t = 0, y
gradiente de temperaturas en el
mantenerse en T∞ en todo
fluido.
momento puede manejarse al hacer
que h tienda al infinito.
La transferencia de calor por
convección se expresa con la Ley
La temperatura del cuerpo cambia
del Enfriamiento de Newton:
de la temperatura inicial Ti a la de
los alrededores T∞ al final del

LABORATORIO DE TRANSFERENCIA DE CALORPágina 2

 CARTAS DE HEISLER Y DIFERENCIAS FINITAS

proceso transitorio de conducción  Cronómetro.

de calor. Por tanto, la cantidad  Rectángulo de fibra de vidrio.
máxima de calor que un cuerpo  Pie de Rey.
puede ganar o perder es
sencillamente el cambio en el
contenido de energía del cuerpo.

Entonces, el uso de los diagramas

de Heisler queda limitado a las
condiciones especificadas, el cual el
cuerpo está inicialmente a una
temperatura uniforme, la
temperatura del medio que lo
circunda y el coeficiente de
transferencia de calor por
convección son constantes y
uniformes y no hay generación de
energía en dicho cuerpo. (Cengel,
1996).
Figura 1: Montaje de la práctica
“cartas de Heisler”.

MATERIALES Y MÉTODOS PROCEDIMIENTO

La práctica se realizó en las

instalaciones del laboratorio de
operaciones unitarias de la 1. Se calentó previamente la
Universidad Pontificia Bolivariana, probeta de arena hasta una
sede Medellín, con una temperatura temperatura aproximada de
promedio de 21ºC. para este 100 a 120ºC en una mufla
experimento se empleó el método durante 2 horas con un
de transferencia de calor por setpoint.
convección natural, éste se 2. Cuando las temperaturas en
fundamenta en la perdida de calor todos los nodos y puntos de
por acción de un gradiente de la probeta sean iguales, se
temperatura provocada una realizaron las mediciones de
temperatura ambiental mucho las temperaturas.
menor al cual es expuesto 3. Se instalaron tres
produciendo un equilibrio entre el termómetros digitales en tres
cuerpo y el exterior. (Cengel, 1996). termocuplas a diferentes
distancias, medidas por un
En el montaje de la práctica (Figura pie de rey.
1), se emplearon los siguientes 4. luego se midió la temperatura
materiales: a la cual estaba cada punto
de la probeta indicados por lo
 Mufla. termómetros, en un lapso de
tiempo de cada 5 minutos
 Probeta de Arena.
hasta el enfriamiento.
 Termómetro Digital.
5. Se anotaron los datos y se
 Selector de Temperatura con tabularon.
6 termocuplas.

LABORATORIO DE TRANSFERENCIA DE CALORPágina 3

 CARTAS DE HEISLER Y DIFERENCIAS FINITAS

CALCULOS Y RESULTADOS

Propiedades Térmicas y Físicas

de la Arena
Densida
d 1400 Kg/m^3
Cp 795 J/Kg*K
K 0,455 W/m* ºC
α 0,000004085 M^2/s
h 10 W/m^2*K
Diámetro 0,1375 M
Too 21,4 ºC
Ti 120 ºC
r0 0,06875 M

Tabla 1: Propiedades termofísicas

de la arena.

D =0.062m

Lc=0.0155m

mRadio
1 0,00064
Radio 2 0,0254
Radio 3 0,0516

Tabla 2: Radios de la Probeta.

LABORATORIO DE TRANSFERENCIA DE CALORPágina 4

 CARTAS DE HEISLER Y DIFERENCIAS FINITAS

Tiempo(s) T1(ºC) T2(ºC) T3(ºC) 0 0)1 0)2 0)3

0 118,9 95,6 81 0,78194726 0,98884381 0,7525355 0,60446247
300 117,9 94,4 77
0,76098715 0,97870183 0,74036511 0,56389452
600 116,5 92,6 74
0,74002705 0,96450304 0,72210953 0,53346856
900 114,6 90,4 71
1200 112,4 88,2 68 0,71602434 0,94523327 0,69979716 0,5030426
1500 110,1 85,8 66 0,69100744 0,92292089 0,67748479 0,47261663
1800 107,6 83,5 63 0,668357 0,89959432 0,65314402 0,45233266
2100 104,9 81,3 61 0,64198783 0,87423935 0,62981744 0,42190669
2400 102,3 79 59 0,61866126 0,84685598 0,60750507 0,40162272
2700 99,6 76,7 57 0,59533469 0,82048682 0,5841785 0,38133874
3000 96,8 74,5 56 0,57167005 0,79310345 0,56085193 0,36105477
0,55138607 0,76470588 0,53853955 0,35091278
Tabla 3: Mediciones realizadas
durante la práctica.

Tabla 5: Cálculo de los Números a

dimensionales de Temperatura

Temperatur
a Central
(ºC) Lc (m) Fo
98,5 0,034375 0
1,0371173
96,4333333   6
2,0742347
94,3666667   1
3,1113520
92   7
4,1484694
89,5333333   2 Para la realización de los cálculos
5,1855867 se tuvieron en cuenta los siguientes
87,3   8
6,2227041 modelos matemáticos:
84,7   3
7,2598214
82,4   9
8,2969388
80,1   4
77,7666667   9,3340562
10,371173
75,7666667   6

Tabla 4: Cálculos y resultados de la

práctica.

LABORATORIO DE TRANSFERENCIA DE CALORPágina 5

 CARTAS DE HEISLER Y DIFERENCIAS FINITAS

T 1 +T 2 +T 3 configuraciones geométricas,
Tcentral= pero la solución comprende
2
series infinitas, las cuales son
difíciles de resolver. Sin
Diámetro embargo, los términos en la
Lc=
4 solución convergen con rapidez
al crecer el tiempo y para
α∗t Fo>0.2, conservar el primer
Fo=
Lc 2 término y despreciar todos los
restantes en la serie conlleva a
un error no mayor del 2%. Por
( Tcentral−T ∞ )
θ0 = consiguiente resulta muy
( T i−T ∞ ) apropiado expresar la solución
utilizando la expresión a un
( T 1 −T ∞ ) término como solución de la
θ/θ01= ecuación diferencial
( T i −T ∞ ) adimensional que rige el
fenómeno de transferencia de
calor en la probeta, este sería:
( T 2 −T ∞ )
θ/θ02=
( T i −T ∞ ) ( T (r , t )−T ∞ ) −λ
12
θ( r , t ) = = A1e J 0 ( λ1 r /r 0 )
( T i −T ∞ )
( T 3 −T ∞ ) Donde las constantes A1 y λ1
θ/θ03= dependen de Bi, cuyos valores
( T i −T ∞ ) se encuentran ya tabulados. La
función J0 es la función de
Los cálculos se realizaron con Excel Bessel de primera clase y de
y evaluados para cada tiempo. orden cero cuyo valor está en
tablas y r el radio del cilindro.
Estas relaciones suelen
simplificarse para dar las
PREGUNTAS relaciones en el centro de la
figura de la siguiente manera:

1. ¿Cuál es la solución
( T 0−T ∞ ) −λ
12
θ0 = = A1 e
aproximada de la ecuación ( T i−T ∞ )
diferencial adimensional que
rige el fenómeno de Entonces las condiciones para
transferencia de calor en la que se puedan usar las cartas
probeta? Detalle los términos de Heisler es que Bi>0.1 y
que la componen. ¿Cuáles Fo>0.2.(Cengel, 1996).
son las condiciones de Bi y
Fo para utilizar las cartas?

RTA: El problema de conducción 2. ¿Se pueden utilizar las cartas

de calor unidimensional en de Heisler para analizar
régimen transitorio se puede probetas con resistencia
resolver son más exactitud para externa despreciable?
cualquiera de las tres Explique.

LABORATORIO DE TRANSFERENCIA DE CALORPágina 6

 CARTAS DE HEISLER Y DIFERENCIAS FINITAS

RTA: El número de Biot es la 0,9

razón de la resistencia interna de 0,8
un cuerpo a la conducción de 0,7
calor con respecto a su 0,6
resistencia externa a la
0,5
convección de calor, así:
0,4
y = -0,0235x + 0,8082
Re sistenciaConducción 0,3 R2 = 0,9995
Bi=
Re sistenciaConvección 0,2
0,1
Entonces, si trabajamos con una 0
resistencia externa despreciable, 1 2 3 4 5 6 7 8 9 10 11
lo cual me indica que es una
valor muy pequeño, esto nos
daría como resultado un número La extrapolación es el método más
de Biot muchísimo mayor a 0.1, habitual de pronóstico. Se basa en
lo que nos dice que se pueden suponer que el curso de los
usar las cartas de Heisler ya que acontecimientos continuará en la
al ser Bi>0.41 se cumple la misma dirección. Despejando x de
condición para el usos de éstos la ecuación me daría un punto (Fo)
diagramas, a pesar de ser alejado que se intercepta con 0 de
menos preciso ya que a medida 0.001 y debe de seguir la misma
de que Bi sea mayor los datos tendencia con un decrecimiento.
son más imprecisos que cuando Con esto el valor sería:
Bi = 0 ó Bi<0.1. (Cengel, 1996).
Y =−0. 0235 X +0 . 8082
3. Grafique 0 vs Fo (Para Y −0 . 8082
Fo>0.2). Extrapole hasta 0= X=
−0 .0235
0.001. Compare la curva 0 .001−0 .8082
obtenida con la de la carta de X=
Heisler. −0 .0235
X =34 .35

0,9
0,8
0,7
0,6
0,5
0,4
0,3
0,2
0,1
0
1 2 3 4 5 6 7 8 9 10 11 12

RTA:

Como se puede ver la curva

mostrada en la imagen es parecida
a las curvas de las cartas de Heisler

LABORATORIO DE TRANSFERENCIA DE CALORPágina 7

 CARTAS DE HEISLER Y DIFERENCIAS FINITAS

lo cual nos indica que la probeta temperatura y es por esto que

experimenta el enfriamiento por experimenta un decrecimiento a
convección natural ya que la curva medida que pasa el tiempo.
empieza a decrecer a medida que
pasa el tiempo. 5. ¿Detalle las causas de error
en los resultados y como
4. Compare la distribución de afectan estas en la práctica?
temperaturas en el cilindro
(0) con las de la carta de RTA: existen varias causas por
Heisler. Como es la variación las cuales se presentaron
de (0) con respecto al errores al momento de calcular,
tiempo para Fo>0.2. Explique estas son:
que se debe.
 Se declararon como
1,2 constantes variables
que en realidad
1 dependen de la
0,8
temperatura como h,
Serie2
k, Cp, entre otros.
0,6
Serie3
 Los datos obtenidos
Serie4
están en un lapso
0,4
corto de tiempo lo cual
0,2 no nos permite ver el
comportamiento en la
0 graficas del proceso,
1 2 3 4 5 6 7 8 9 10 11
con esto, la
comparación de estos
RTA: Se observó que las curvas con los de las cartas
presentan la tendencia que de Heisler se dificulta
demarcan las curvas de las ya que sólo se tiene
cartas de Heisler, aunque las un tramo de la
mostradas en la anterior trayectoria.
ilustración sólo es una trama de  Errores de operarios
la curva completa ya que los se evidenciaron en la
datos obtenidos no fueron los práctica ya que hubo
suficientes para observar el mucha corriente de
comportamiento como tal del aire lo cual causa
proceso de transferencia de alteración en el
calor por convección natural. proceso, ya que este
debe darse de forma
Esta tendencia que presenta natural.
esta gráfica es debido a que
inicialmente la probeta estaba a
una Ti equivalente
aproximadamente a 120ºC, ésta
es sujeta a convección de
manera natural por toda la
superficie de la probeta hacia un
medio ambiente a la temperatura 6. Graficar la curva de
T∞, la probeta alcanza ésta enfriamiento Ti vs t para cada

LABORATORIO DE TRANSFERENCIA DE CALORPágina 8

 CARTAS DE HEISLER Y DIFERENCIAS FINITAS

punto de medición i. (Todas De acuerdo a los datos

las curvas en una misma obtenidos se puede notar que no
gráfica). está descrito completamente el
comportamiento de la
140 temperatura y demás variables
120
en función del tiempo, lo
podemos observar en las
100 graficas, las cuales sólo nos
80
muestra una trama de la curva la
Serie2
cual debería experimentar un
60 Serie3
leve decrecimiento al pasar del
40
Serie4 tiempo; esto es debido, a que no
se tomaron los suficientes datos
20
de tiempo y temperatura al
0 momento de la realización de la
1 2 3 4 5 6 7 8 9 10 11 práctica, es decir, no se terminó
el proceso de transferencia de
calor por convección natural.
7. Graficar la curva de Esto no nos dice que los datos
temperaturas T vs r para 0s, tomados sean incorrectos, sino
600s, 1200s, 3600s, 5400s y que sólo podemos observar la
7200s. (Todas las curvas en primera parte de la curva, esto lo
una misma grafica). podemos ver al momento de
hacer la comparación con las
140 cartas de Heisler para la figura
120
cilíndrica ya que presenta un
decrecimiento leve en la parte
100 superior y en la parte final el
80 decrecimiento es muchísimo
mayor. Lo que se observa en
60
nuestras gráficas es el leve
40 decrecimiento mas no lo
siguiente que es la caída brusca
20
de la temperatura con respecto
0 al tiempo, lo mismo sucede al
1 2 3 4 5 6 7 8 9 10 11 comparar el número
adimensional de temperatura, el
número de Fourier y demás con
respecto al tiempo.

El proceso se hizo lento debido a

que la probeta fue calentada por
encima de los 100ºC, con una
temperatura de 120ºC, lo que
hizo que el comportamiento
observado en las graficas sea
menor y sólo se vea esa leve
disminución en la temperatura y
ANÁLISIS DE RESULTADOS demás variables con respecto al
tiempo.

LABORATORIO DE TRANSFERENCIA DE CALORPágina 9

 CARTAS DE HEISLER Y DIFERENCIAS FINITAS

También se incurrieron en
errores y es al declarar variables
dependientes de la temperatura BIBLIOGRAFÍAS
como lo son h, k, Cp, Densidad,
etc., como constantes, esto hace  Yunus A. Cengel.
que los cálculos sean menos Termodinámica, 2ª Edición,
precisos y la trayectoria Editorial Mc Graw Hill, 1996.
demarcada por la gráfica no sea 581 Páginas.
la esperada. Reiterando que los  [Internet]. URL:
valores obtenidos sólo nos http://es.wikipedia.org/wiki/Co
permite ver un intervalo corto de nvecci%C3%B3n. Accedido:
tiempo y temperatura y esto [26/11/2009]
hace que no se observe el
comportamiento como tal, esto lo
podemos demostrar al momento
de hacer la extrapolación ya que
al buscar un valor para Fo con
un valor de 0 de 0.001,
notamos que decae
bruscamente la curva indicando
que la probeta y el proceso de
transferencia de calor que
ocurrió en ella es el indicado y
antes estudiado.

CONCLUSIONES

 Las cartas de Heisler

presentan limitaciones y
es en su poca exactitud
en los valores que brinda.
 Estos diagramas nos
permiten evaluar
procesos largos y con una
gran cantidad de datos
como los obtenidos
durante la práctica de
transferencia de calor por
convección natural.
 Para poder analizar por
medio de las cartas de
Heisler un proceso de
transferencia de calor por
convección natural se
necesita que el proceso
se de completamente ya
que de esta manera el
comportamiento en las
graficas es el adecuado.

LABORATORIO DE TRANSFERENCIA DE CALORPágina 10