LABORATORIO DE OPERACIONES UNITARIAS: CALOR Y MASA

PRÁCTICA 1
TRANSFERENCIA DE CALOR EN RÉGIMEN TRANSITORIO: PROCESO
CON DEPENDENCIA SÓLO DEL TIEMPO Y PROCESOS CON
DEPENDENCIA DEL TIEMPO Y EL ESPACIO

1. OBJETIVO
Estudiar el mecanismo de transferencia de calor en régimen transitorio por convección
con flujo externo sobre superficies cilíndricas

1.1 OBJETIVOS ESPECÍFICOS

 Determinar experimentalmente los coeficientes de la correlación de Hilpert para flujo
externo en cilindros metálicos
 Hallar los coeficientes de convección promedio de transferencia de calor.
 Encontrar la correlación experimental del número de Nusselt en función del número
Reynolds y Prandlt.
 Hallar el perfil de temperatura, experimental y teórica, en el centro de un cilindro de
resistencia interna considerable en función del tiempo.

1.2 ASPECTOS TEÓRICOS

En la convección forzada, el movimiento relativo entre el fluido y la superficie se mantiene

por medios externos como un ventilador o una bomba, y no por fuerzas de empuje

debidas a los gradientes de temperatura en el fluido.

En la práctica existen numerosas situaciones en las que se necesita analizar o describir

la transferencia de calor asociada con el flujo de un fluido alrededor de la superficie

exterior de un sólido, particularmente una esfera o un cilindro.

El número de Nusselt se puede determinar experimentalmente, mediante la medición de

la transferencia de calor en condiciones controladas y correlacionar los datos en términos

de parámetros adimensionales apropiados.

1
Zulamita Zapata Benabithe
Oscar Vasco Echeverri
 LABORATORIO DE OPERACIONES UNITARIAS: CALOR Y MASA

̅̅̅̅ = 𝑓 (𝑅𝑒, 𝑃𝑟)

𝑁𝑢 (1.1)

El número Nusselt promedio para un flujo externo alrededor de un cilindro circular se

expresa de la forma:
ℎ𝐷̅
̅̅̅̅
𝑁𝑢𝐷 = 𝑘 (1.2)

Los resultados asociados con un fluido dado, y por ello un número de Prandtl fijo, caen

cerca de una línea recta, que se puede representar mediante la forma:

̅̅̅̅𝐷 = 𝐶𝑅𝑒𝐷𝑚 𝑃𝑟 𝑛
𝑁𝑢 (1.3)

Los valores de C, m y n a menudo son independientes de la naturaleza del fluido. Las

propiedades necesarias para utilizar la ecuación 1.3, se evalúan a una temperatura media
de la capa límite o temperatura de película.

En ingeniería, el principal interés es hallar un coeficiente convectivo promedio (del aire)

para toda la superficie. Las diferentes correlaciones existentes, entre ellas la de Hilpert,

proporcionan una forma sencilla de obtener dicho coeficiente.

1
̅̅̅̅
𝑁𝑢𝐷 = 𝐶𝑅𝑒𝐷𝑚 𝑃𝑟 ⁄3 𝑃𝑟 ≥ 0.7 (1.4)

Donde C y m son constantes que dependen del número de Reynolds.

Conducción transitoria
En el proceso de transferencia de calor transitoria (calentamiento o enfriamiento) tiene
lugar en el tiempo intermedio antes de alcanzarse el equilibrio (transferencia de calor
estacionaria), debe modificarse el análisis para tener en cuenta la variación de la energía
interna del cuerpo, y las condiciones de contorno deben ajustarse para que encajen con
la situación física que se plantea en el problema de transferencia de calor n estacionaria.

2
Zulamita Zapata Benabithe
Oscar Vasco Echeverri
 LABORATORIO DE OPERACIONES UNITARIAS: CALOR Y MASA

Muchos de los problemas más prácticos de conducción de calor no implican el flujo de

calor estacionario entre dos receptáculos térmicos, sino en vez de eso, la repuesta
transitoria a los cambios de temperatura. La calefacción de los edificios, el proceso de
cocción de los alimentos, el vaciado de hierro de fundición y el tratamiento térmico de los
metales (recocido, temple, revenido, etc.), todos ellos son casos en donde la tasa se
cambio de la temperatura con el tiempo son importantes.

Formulación Concentrada: resistencia interna despreciable

El análisis denominado formulación concentrada, constituye un análisis simplificado al
problema de conducción transitoria. Si un cuerpo de forma arbitraria que inicialmente se
encuentra a una temperatura inicial T i y que es sumergido en un recipiente de grandes
dimensiones que se encuentra a una temperatura T ∞. Se satisface que Ti >T∞ de manera
que el cuerpo estará sometido a un proceso de enfriamiento.
La hipótesis más importante que se realiza en el análisis de Resistencia interna
despreciable, es que se considera que la distribución de temperatura dentro del cuerpo
es uniforme, de tal manera que la temperatura depende exclusivamente del tiempo, esta
suposición en principio es aceptable bajo las condiciones de: cuerpo de dimensiones
pequeñas, alta conductividad térmica del cuerpo y bajo coeficiente de transferencia de
calor por convección.

De la primera Ley de la Termodinámica se puede llegar a la Ecuación 3.1, la cual permite

determinar la distribución de temperaturas.

𝑑𝑇 ℎ𝐴𝑠
= (𝑇 − 𝑇∞) (3.1)
𝑑𝑡 𝜌 ∙ 𝑉 ∙ 𝐶𝑝

La Ecuación 3.1 es de primer orden y se necesita suponer una condición inicial, T=T i en
t=0. Incorporando las variables adimensionales del número de Biot, Fourier y una variable
que involucre las temperaturas (Ec. 3.2).

𝜃 ∗ = 𝑒 (−𝐵𝑖𝑚∙ 𝐹𝑜𝑚) (3.2)

3
Zulamita Zapata Benabithe
Oscar Vasco Echeverri
 LABORATORIO DE OPERACIONES UNITARIAS: CALOR Y MASA

Donde:
θ*: Es la relación de las temperaturas, θ*=(T-T∞)/(Ti-T∞)
Bim: Número de Biot modificado
Fom: Número de Fourier modificado

Interpretación Física del Número de Biot y del Número de Fourier

El número de Biot compara los valores relativos de la resistencia a la transferencia de
calor por convección en la superficie y la resistencia interna a la conducción. Un valor
muy bajo significa que la resistencia a la conducción es despreciable en comparación con
la resistencia convectiva de la superficie. Esto implica que la temperatura será
prácticamente uniforme a lo largo del sólido.

El Número de Biot modificado posee la interpretación física mostrada en la Ec. 3.3.

𝐿𝑐⁄
[𝑅𝑒𝑠𝑖𝑠𝑡𝑒𝑛𝑐𝑖𝑎 𝑝𝑜𝑟 𝑐𝑜𝑛𝑑𝑢𝑐𝑐𝑖ó𝑛] 𝑘
𝐵𝑖𝑚 = = (3.3)
[𝑅𝑒𝑠𝑖𝑠𝑡𝑒𝑛𝑐𝑖𝑎 𝑝𝑜𝑟 𝑐𝑜𝑛𝑣𝑒𝑐𝑐𝑖ó𝑛] 1⁄
ℎ

Donde la longitud característica (Lc), corresponde a la relación entre el volumen y el

área superficial del cuerpo.

Bim > 0.1 Implica que la resistencia por conducción domina

Bim < 0.1 Implica que la resistencia por conducción es despreciable

El número de Fourier (Ec. 3.4) es una medida del calor conducido a través de un cuerpo
con relación al calor almacenado en él. Por tanto, un valor grande del número de Fourier
indica una propagación más rápida de calor a través del cuerpo. Además compara una
longitud característica del cuerpo con un valor aproximado de la longitud hasta la que
penetra la onda de temperatura en un tiempo dado t.

𝛼∙𝑡
𝐹𝑜 = (3.4)
𝑠2

4
Zulamita Zapata Benabithe
Oscar Vasco Echeverri
 LABORATORIO DE OPERACIONES UNITARIAS: CALOR Y MASA

En donde α, es la difusividad térmica del material, la cual es una propiedad térmica del
sólido, que se define como la mostrada en la Ecuación 3.5.
[𝐸𝑛𝑒𝑟𝑔í𝑎 𝑞𝑢𝑒 𝑐𝑜𝑛𝑑𝑢𝑐𝑒] 𝑘
𝛼= = (3.5)
[𝐸𝑛𝑒𝑟𝑔í𝑎 𝑞𝑢𝑒 𝑎𝑙𝑚𝑎𝑐𝑒𝑛𝑎] 𝜌 ∙ 𝐶𝑝

Determinación del calor total

La temperatura del sólido en cualquier tiempo t puede calcularse a partir de la Ecuación
3.6, a partir del sólido de resistencia interna despreciable (Bi < 0.1).

𝑞(𝑡) = ℎ𝑐 ∙ 𝐴 ∙ (𝑇 − 𝑇𝑎 ) (3.6)

Si se sustituye la temperatura instantánea T de la Ec. 3.1 en la Ec. 3.6 se llega a la

Ecuación 3.7.
ℎ ∙𝐴
−( 𝑐 ∙ 𝜌∙𝑉)∙𝑡
𝐶𝑝
𝑞(𝑡) = ℎ𝑐 ∙ 𝐴 ∙ (𝑇𝑜 − 𝑇∞ ) ∙ 𝑒 (3.7)

Para determinar la cantidad total de calor transferida por el sólido desde el tiempo t = 0
hasta t = t, se puede integrar la Ec. 3.7 y se obtiene la Ec. 3.8.

ℎ ∙𝐴
−( 𝑐 ∙𝜌∙𝑉)∙𝑡
𝐶𝑝
𝑞 = 𝐶𝑝 ∙ 𝜌 ∙ 𝑉 ∙ (𝑇𝑜 − 𝑇∞ ) ∙ [1 − 𝑒 ] (3.8)

Donde To es la temperatura en el instante t = 0.

La Ec. 3.8 expresa el número de calor que puede ser transferido por el cuerpo, que se
presenta cuando transcurre un tiempo lo suficientemente largo, tiempo tiende a infinito,
en donde se alcanza el equilibrio térmico y por tanto el cuerpo adquiere la temperatura
del medio circundante.

En general se acepta que la teoría de la resistencia interna despreciable es válida siempre

y cuando el número del Bi sea inferior a 0.1, criterio para la aplicación del análisis de
resistencia interna despreciable.
5
Zulamita Zapata Benabithe
Oscar Vasco Echeverri
 LABORATORIO DE OPERACIONES UNITARIAS: CALOR Y MASA

Cálculo de la longitud característica para las geometrías primarias, Lc.

La longitud característica, Lc, es una dimensión característica del cuerpo, obtenida a

partir de la Ecuación 3.9.

[𝑉𝑜𝑙𝑢𝑚𝑒𝑛 𝑑𝑒𝑙 𝑐𝑢𝑒𝑟𝑝𝑜] 𝑉

𝐿𝑐 = = (3.9)
[Á𝑟𝑒𝑎 𝑐𝑜𝑛𝑣𝑒𝑐𝑡𝑖𝑣𝑎 𝑑𝑒𝑙 𝑐𝑢𝑒𝑟𝑝𝑜] 𝐴

La longitud característica para el cilindro es el radio (R/2) con volumen πR2L, área
superficial de 2πRL y si la longitud característica se considera infinita es R.

En algunas ocasiones la longitud característica puede simplificarse debido a la

consideración de una longitud infinita. Por ejemplo, para un cilindro infinito la longitud
característica es el radio.

Uso de los diagramas de Heisler

Cuando el número de Biot modificado supera el valor de 0.1, la resistencia interna ya no
es despreciable, de manera que hay que recurrir a los diagramas de Heisler o utilizar una
solución exacta. Para la aplicación de los diagramas de Heisler deben determinarse los
números de Biot y los números de Fourier. Además de las longitudes características para
las diferentes geometrías infinitas. Los cálculos para los diagramas de Heisler se han
realizado truncando las series infinitas de los problemas y reduciéndolas a unos pocos
términos, esto restringe el campo de aplicación de los diagramas para valores de Fo ≥
0.2

Heisler se dio a la tarea de graficar dichas soluciones con la evaluación del primer término,
el cual determina la Ecuación 3.10.

𝑇0 − 𝑇∞ 𝜃0
𝜃∗ = = (3.10)
𝑇𝑖 − 𝑇∞ 𝜃𝑖

6
Zulamita Zapata Benabithe
Oscar Vasco Echeverri
 LABORATORIO DE OPERACIONES UNITARIAS: CALOR Y MASA

Donde T0 es la temperatura del sólido en cualquier instante de tiempo, Ti es la

temperatura inicial del sólido, las dos medidas en la superficie (se supone que ésta es
uniforme en todo el sólido) y T ∞ es la temperatura del medio en el que se encuentra el
sólido.

Luego de haber empleado las cartas de Heisler, se tiene que la temperatura en el centro
del sólido después de un determinado tiempo t está dado por la Ec. 3.11.

𝜃0
𝑇0 = ∙ (𝑇𝑖 − 𝑇∞ ) + 𝑇∞ (3.11)
𝜃𝑖

NOTA: Tener en cuenta que para la placa en los diagramas de Heisler se lee con 1/Biot,
para un cilindro 1/(2*Biot), y para una esfera 1/(3*Biot).

1.3 EQUIPO Y MATERIAL REQUERIDO

 Mufla
 Probetas de metálicas de 2"
 2 termopares
 2 indicadores de temperatura
 Anemómetro
 Túnel de viento
 Cronómetro
 Soporte cerámico
 Baño de Temperatura
 Salchicha tipo llanera

1.4 PREINFORME

Se debe tener clara la teoría sobre convección externa sobre una superficie cilíndrica y
la determinación experimental del coeficiente de convección en función del tiempo. Se
debe llevar a la práctica el modelo de cálculo.

7
Zulamita Zapata Benabithe
Oscar Vasco Echeverri
 LABORATORIO DE OPERACIONES UNITARIAS: CALOR Y MASA

1.5 PROCEDIMIENTO EXPERIMENTAL

1.5.1 Convección forzada flujo externo

 Calentar las probetas de igual diámetro a una temperatura de 150°C.
 Medir las condiciones de velocidad y temperatura del aire a las cuales se van a hacer
el enfriamiento de las probetas.
 Sacar una de las probetas de la estufa y ubicarla sobre un soporte de cerámico justo
en el extremo del túnel.
 Iniciar el registro de la temperatura cuando la temperatura haya llegado a un valor
estable e inicie el descenso en convección natural.
 Encender el ventilador del túnel de viento a una velocidad establecida y continuar con
el registro de temperatura cada 30 s.
 Realizar 4 enfriamientos a diferentes velocidades para las diferentes probetas.
 Medir para cada probeta masa y dimensiones.
1.5.1 Convección forzada flujo natural
 Tomar la salchicha e introducirle un termopar teniendo en cuenta que sea por el
centro.
 Sumergir la salchicha con el termopar al baño térmico y registrar la variación de la
temperatura en el tiempo hasta alcanzar un valor estable en el centro.
 Sacar la salchicha del baño, secarla y dejarla enfriar por convección natural. Registrar
el cambio de temperatura en función del tiempo durante en el enfriamiento.
 Medir las dimensiones y peso de la salchicha. Algunos datos para la salchicha se
pueden consultar en Matteson y Sommerfeld (1984).
 Suponer el valor del coeficiente de convección del aire.

1.6 ELABORACIÓN DEL INFORME

El informe se debe entregar cuando se realice la siguiente práctica de laboratorio, de

forma grupal. Debe contener los datos y resultados tabulados, análisis de resultados,
conclusiones, recomendaciones y bibliografía.
8
Zulamita Zapata Benabithe
Oscar Vasco Echeverri
 LABORATORIO DE OPERACIONES UNITARIAS: CALOR Y MASA

 Tabular y graficar los datos de las curvas de enfriamiento (t vs tiempo).

 Determinar la temperatura promedio de película y hallar Pr,  y k del aire a la
temperatura medida.
 Calcular los números Biot, Reynolds y Nusselt por medio de la correlación de Hilpert
y experimentalmente mediante el balance de energía.
 Graficar el Nusselt experimental y por correlaciones versus el número de Reynolds
para cada probeta.
 A partir de la gráfica anterior obtener una regresión lineal y determinar los coeficientes
para construir una nueva correlación (C y m). Comparar los valores obtenidos con los
de la correlación de Hilpert.
 Mediante solución exacta y cartas de Heisler predecir el valor de la temperatura en
el centro de la salchicha y comparar con los valores experimentales.
 Determinar el valor de la cantidad total de calor transferida por el sólido para ambos
procedimientos.

1.7 TOMA DE DATOS EXPERIMENTALES

Tabular los valores de temperatura en función del tiempo para un valor de velocidad
determinado. Medir las dimensiones de las probetas y determinar las propiedades
necesarias (Cp, k).

Tabla 1. Toma de datos en el laboratorio

Biot < 0.1 Biot < 0.1 Biot < 0.1 Biot < 0.1
Biot > 0.1
Probeta 1 Probeta 2 Probeta 3 Probeta 4
Tiempo Temp. Tiempo Temp. Tiempo Temp. Tiempo Temp. Tiempo Temp.
(seg) (°C) (seg) (°C) (seg) (°C) (seg) (°C) (seg) (°C)

9
Zulamita Zapata Benabithe
Oscar Vasco Echeverri
 LABORATORIO DE OPERACIONES UNITARIAS: CALOR Y MASA

Biot < 0.1 Biot < 0.1 Biot < 0.1 Biot < 0.1
Biot > 0.1
Probeta 1 Probeta 2 Probeta 3 Probeta 4
Tiempo Temp. Tiempo Temp. Tiempo Temp. Tiempo Temp. Tiempo Temp.
(seg) (°C) (seg) (°C) (seg) (°C) (seg) (°C) (seg) (°C)

1.8 BIBLIOGRAFÍA
 Vasco, Oscar y Ocampo, Carlos. Manual de laboratorio, transferencia de calor.
Editorial Universidad Pontificia Bolivariana. ISBN: 978-958-764-097-7. 2013.
 Vasco, Oscar. Laboratorio de transferencia de masa. Editorial Universidad Pontificia
Bolivariana. ISBN: 978-958-764-184-4. 2014.
 Frank Incropera, David DeWitt. Fundamentos de Transferencia de Calor. Prentice
Hall. México 1999 (4° Ed). ISBN 0-471-30460-3. 888 pp.
 Frank Kreith y Mark Bohn. Principios de Transferencia de Calor. Thomson Learning.
México 2001 (6° ed). 700 pp.
 Lienhard IV, John H. and Lienhard V, John H. A heat transfer textbook. Third edition.
Cambridge,MA. 2005.
 J.P. Holman Transferencia de Calor 8.a Edición (1.a Español) McGraw Hill.
 Matteson, M. y Sommerfeld, J. Thermal conductivity of a Hotdog. ChE division. 1984.

10
Zulamita Zapata Benabithe
Oscar Vasco Echeverri