Transferencia Lab 2
Transferencia Lab 2
Transferencia Lab 2
1. INTRODUCCIÓN
1.1. ANTECEDENTES
La característica del régimen transitorio que son procesos de transferencia de calor en donde
la temperatura no solo depende de parámetros espaciales sino que también depende del
tiempo; los modelos referentes de este régimen se encuentran fuertemente relacionados con
el mecanismo de convección.
1.2. OBJETIVOS
En la realidad toda transferencia de calor pasa por una etapa transitoria, es decir que la
temperatura varía en función del espesor y del tiempo transcurrido (𝑇 = 𝑇(𝑥, 𝜃). Sin embargo,
en algunos casos, esta etapa es corta en comparación al régimen permanente de tal forma
que se la puede despreciar.
1.3.2. CONVECCIÓN
La ecuación que nos permite conocer la cantidad de calor que se transfiere mediante este
mecanismo se debe a la ley de enfriamiento de Newton, que indica: “El flujo de calor por
unidad de área es directamente proporcional a la diferencia de temperaturas que existe entre
la pared sólida y el fluido global, donde el factor de proporcionalidad se conoce como el
coeficiente pelicular de convección”
𝑄 = 𝐴 ℎ (𝑇𝑊 − 𝑇∞ ) (1.2)
Dónde:
𝑞 [𝑊⁄ 2 ] ∶ 𝐹𝑙𝑢𝑗𝑜 𝑑𝑒 𝑐𝑎𝑙𝑜𝑟 𝑝𝑜𝑟 𝑢𝑛𝑖𝑑𝑎𝑑 𝑑𝑒 𝑎𝑟𝑒𝑎.
𝑚
Este caso particular se obtiene cuando la temperatura del sistema sujeto a una respuesta
térmica transitoria es casi uniforme, por lo que se pueden ignorar las pequeñas diferencias de
temperatura dentro del mismo; las variaciones en la energía interna del sistema se pueden
expresar en función de las variaciones de la temperatura uniforme, aproximación que se
conoce también como modelo de la capacidad térmica global. La suposición de que la
temperatura es uniforme se justifica, puesto que la resistencia térmica a la conducción dentro
del sólido es pequeña comparada con la resistencia exterior convectiva. Se cumple para un
valor de Bi < 0,1 lo que garantiza el que la temperatura en el interior del sólido, (placa, cilindro,
esfera), no diferirá de la de la superficie en más de un 5%, lo que equivale a suponer que el
fluido debe tener un bajo coeficiente de convección (gases), por lo que el sólido tiende a la
temperatura del fluido en un tiempo relativamente grande, permaneciendo uniforme en cada
instante la temperatura del sólido, incluida su superficie.
La fuerza ascensional por unidad de volumen es igual al peso del aire frio menos el peso del
aire caliente.
𝐹𝑎 = 𝛾𝑎𝑓 − 𝛾𝑎𝑐 = 𝑔(𝜌𝑜 − 𝜌) (1.3)
Del principio de dilatación de los gases se tiene:
𝜌𝑜 − 𝜌 = 𝜌𝑜 𝛽 ∆𝑇 (1.4)
Donde:
𝛽: 𝐸𝑠 𝑒𝑙 𝑐𝑜𝑒𝑓𝑖𝑐𝑖𝑒𝑛𝑡𝑒 𝑑𝑒 𝑒𝑥𝑝𝑎𝑛𝑠𝑖𝑜𝑛 𝑣𝑜𝑙𝑢𝑚𝑒𝑡𝑟𝑖𝑐𝑎.
𝐹𝑎 = 𝑔 (𝜌𝑜 − 𝜌) = 𝑔 𝜌𝑜 𝛽 ∆𝑇 (1.5)
𝑁𝑢 = 𝑓(𝐺𝑟; Pr)
Numero de Grashof:
𝜌2 𝑔 𝐿3 𝛽 ∆𝑇
𝐺𝑟 = (1.6)
𝜇2
Si en vez de partir de un fluido estacionario que se agita por el mero efecto de las diferencias
de temperatura, forzamos el movimiento relativo con otros medios (con un ventilador en
marcha, por ejemplo), tenemos lo que se conoce como convección forzada. Si la velocidad
relativa de partida es lo bastante grande, la que sería provocada por los cambios de
temperatura se hace irrelevante. La «constante» de convección del modelo del enfriamiento
de Newton se hace muy insensible a la temperatura. Además de esto, como la velocidad
relativa puede ser muy grande, la eficacia de la convección forzada puede ser mucho mayor
que la de la convección natural. Esto es algo que se tiene en cuenta, por ejemplo, al diseñar
sistemas de refrigeración para dispositivos electrónicos: cuando la refrigeración es por
convección natural —refrigeración pasiva—, la superficie necesaria para disipar el calor puede
ser muchísimo mayor que la necesaria con un sistema de refrigeración por convección forzada
(refrigeración activa).
La refrigeración por convección forzada puede funcionar cuando la convección natural ni
siquiera es una opción, bien por cuestiones geométricas (quizá sería necesario un disipador
enorme), bien por cuestiones físicas (como sucede en las cargas de pago de la Estación
Espacial Internacional, donde la convección natural no es una opción por el entorno de
microgravedad).
UNIVERSIDAD TECNICA DE ORURO
FACULTAD NACIONAL DE INGENIERIA
CARRERA DE INGENIERIA MECANICA
ELECTROMECANICA
LABORATORIO DE MAQUINAS TERMICAS
TRANSFERENCIA MEC 2251
DE CALOR REGIMEN TRANSITORIO
Y CONVECCION
LABORATORIO 2 Pág. 6 de 34
En el caso del coeficiente de transferencia de calor por convección (h), se utiliza el parámetro
adimensional que se denomina número de Nusselt (Nu)
Proporciona una medida de la transferencia de calor por convección que ocurre en una
superficie.
NUMERO DE PRANDTL
La mejor manera de describir el espesor de las capas límites de velocidad y térmica es por
medio del parámetro adimensional conocido como número de Prandtl.
Tabla 1.1. Rangos Típicos de los Números de Prandtl para fluidos comunes.
NUMERO DE REYNOLDS.
NUMERO DE GRASHOF
𝑅𝑎 = 𝐺𝑟 ∗ 𝑃𝑟 (1.13)
𝑃𝑒 = 𝑅𝑒 ∗ 𝑃𝑟 (1.14)
Flujo Couette, flujo laminar interno completamente desarrollado para condición de pared flujo
de calor constante, solución de Blasius para flujo laminar externo.
Las soluciones exactas son importantes porque las suposiciones que se hacen se pueden
especificar con precisión y su validez puede comprobarse en forma experimental. Presupone
mecanismos físicos conocidos para ser descriptos en lenguaje matemático. Sirven como base
de comparación y como una comprobación para los métodos simples, pero aproximados.
Esto evita la descripción detallada del flujo en la capa límite. Si se usa una ecuación más
simple (ecuaciones polinómicas) para describir las distribuciones de velocidad y temperatura
en la capa límite. El problema se analiza sobre la base macroscópica, aplicando la ecuación
de movimiento y la ecuación de la energía para el agregado de partículas de fluido contenidas
dentro de la capa límite. Da soluciones aproximadas cuya exactitud en general depende del
perfil considerado.
UNIVERSIDAD TECNICA DE ORURO
FACULTAD NACIONAL DE INGENIERIA
CARRERA DE INGENIERIA MECANICA
ELECTROMECANICA
LABORATORIO DE MAQUINAS TERMICAS
TRANSFERENCIA MEC 2251
DE CALOR REGIMEN TRANSITORIO
Y CONVECCION
LABORATORIO 2 Pág. 8 de 34
𝜕𝑇
Necesidad de conocer: 𝜕𝑦 y conocer 𝑇 = 𝑓(𝑦)
Sabiendo que:
𝜕𝑇
−𝑘𝑓𝑙𝑢𝑖𝑑𝑜 ( )
𝜕𝑦 𝑦=0
ℎ= (1.15)
(𝑇𝑤 − 𝑇∞ )
1.3.5. CUESTIONARIO
1. Recopilar y conformar dos tablas, bien referenciadas (mínimo cuatro referencias, con nota
de pie de página), de valores representativos para los coeficientes de convección natural
y forzada, las cuales serán usadas referentemente para adoptar los valores de los
coeficientes de convección en el experimento.
UNIVERSIDAD TECNICA DE ORURO
FACULTAD NACIONAL DE INGENIERIA
CARRERA DE INGENIERIA MECANICA
ELECTROMECANICA
LABORATORIO DE MAQUINAS TERMICAS
TRANSFERENCIA MEC 2251
DE CALOR REGIMEN TRANSITORIO
Y CONVECCION
LABORATORIO 2 Pág. 9 de 34
TABLA 1
CORRELACIONES CONVECCION LIBRE FLUJO EXTERNO
Placa plana Vertical Inclinada
Ta
N° Correlación Condiciones de Aplicación Propiedades Nombre
Tª media de
1 ̅̅̅̅
𝑁𝑢𝐿 = 𝐶 𝑅𝑎𝐿𝑛 Medio, Ts cte película. Mc Adams
1
6
̅̅̅̅ 0,387𝑅𝑎𝐿 2 Tª media de Churchill y
𝑁𝑢𝐿 = {0,825 + 0,492 9⁄ }
2 [1+( ⁄ ) 16 ]4⁄9
𝑃𝑟 Medio, Ts cte, 10−1 < 𝑅𝑎 𝐿 ≤ 1012 película. Chu
1⁄4
̅̅̅̅𝐿 = 0,68 + 0,670 𝑅𝑎𝐿 Tª media de Churchill y
𝑁𝑢 0,492 9
3 [1+( ⁄𝑃𝑟 ) ⁄16]4⁄9 Laminar, medio, Ts cte, 10−1 < 𝑅𝑎 𝐿 ≤ 109 película. Chu
Placa plana Inclinada
Correlaciones (1), (2) y (3) sustituyendo RaL 0° < θ < 60°, Superficie superior de placa fría o Tª media de
4 por RaL cosθ superficie inferior de placa caliente. película.
Placa plana Horizontal
Medio, Ts cte, Superficie superior de placa
caliente o superficie inferior de placa fría, Tª media de
̅̅̅̅𝐿 = 𝐶 𝑅𝑎𝐿𝑛
5 𝑁𝑢 Longitud característica: L = As /P película. McAdams
Medio, Ts cte, Superficie superior de placa
caliente o superficie inferior de placa fría, Tª media de
⁄
𝑁𝑢𝐿 = 0,27 𝑅𝑎𝐿1 4
6 ̅̅̅̅ Longitud característica: L = As /P película. McAdams
Cilindro Largo Horizontal
Tª media de
7 ̅̅̅̅
𝑁𝑢𝐷 = 𝐶 𝑅𝑎𝐷𝑛 Medio, Ts cte película. Morgan
1
6
̅̅̅̅ 0,387𝑅𝑎𝐷 2 Tª media de Churchill y
𝑁𝑢𝐷 = {0,6 + 9 }
8 [1+(0,559⁄ ) ⁄16 ]8⁄27
𝑃𝑟 Medio, Ts cte, 10−5 < 𝑅𝑎 𝐷 ≤ 1012 película. Chu
Flujo externo esfera
1⁄4
̅̅̅̅𝐷 = 2 + 0,589 𝑅𝑎𝐷 Tª media de
𝑁𝑢 9
9 [1+(0,469⁄𝑃𝑟 ) ⁄16 ]4⁄9 Medio, Ts cte, 𝑅𝑎 𝐷 ≤ 1011 , 𝑃𝑟 ≥ 0,7 película. Churchill
Recinto Rectangular
T. media Globe y
⁄
𝑁𝑢𝐿 = 0,069 𝑅𝑎𝐿1 3 𝑃𝑟 0,0704
10 ̅̅̅̅ Medio, Ts cte, θ = 0º, 3 ∗ 105 < 𝑅𝑎 𝐿 ≤ 7 ∗ 109 superficies Dropkin
0,29 Medio, Ts cte, θ = 90º, 103 < Pr RaL/0,2+Pr T. media
̅̅̅̅𝐿 = 0,18 ( 𝑅𝑎𝐿 𝑃𝑟 )
𝑁𝑢
11 0,2+𝑃𝑟 10−3 < 𝑅𝑎 𝐿 ≤ 105 , 1 < H/L < 2 superficies Catton
𝑅𝑎𝐿 𝑃𝑟 0,29 𝐻 −1⁄4 Medio, Ts cte, θ = 90º, 103 < RaL < 1010 Pr < T. media
̅̅̅̅
12 𝑁𝑢𝐿 = 0,22 (0,2+𝑃𝑟) ( )
𝐿 105 , 2 < H/L < 10 superficies Catton
⁄ 𝐻 −0,3 Medio, Ts cte, θ = 90º, 104 < RaL < 107 1 < Pr < T. media MacGregor
𝑁𝑢𝐿 = 0,42 𝑅𝑎𝐿1 4 𝑃𝑟 0,012 ( )
13 ̅̅̅̅ 2x104 , 10 < H/L < 40 superficies y Emery
𝐿
Medio, Ts cte, θ = 90º, 106 < RaL < 109 1 < Pr < T. media MacGregor
⁄
𝑁𝑢𝐿 = 0,046 𝑅𝑎𝐿1 3
14 ̅̅̅̅ 20, 1 < H/L < 40 superficies y Emery
1708 ∗
̅̅̅̅
𝑁𝑢𝐿 = 1 + 1,444 [1 − ] [1 −
𝑅𝑎𝐿 𝑐𝑜𝑠𝜃
∗
1708 (𝑠𝑒𝑛1,8𝜃)1,6 𝑅𝑎𝐿 𝑐𝑜𝑠𝜃 1⁄3 Medio, Ts cte, 0 < θ ≤ θcr, H/L ≥ 12 H/L 1 3 6 12 T. media
15 ] + [( ) − 1] >12 θcr 25º 53º 60º 67º 70º superficies Hollands
𝑅𝑎𝐿 𝑐𝑜𝑠𝜃 5830
̅̅̅̅𝐿 = 𝑁𝑢
𝑁𝑢 ̅̅̅̅𝐿 (𝜃 = 0°)
𝜃 ⁄𝜃𝑐𝑟
̅̅̅̅
𝑁𝑢𝐿 (𝜃 = 0°) T. media
[ ] (𝑠𝑒𝑛𝜃𝑐𝑟 )𝜃⁄4𝜃𝑐𝑟
16 ̅̅̅̅𝐿 (𝜃 = 90°)
𝑁𝑢 Medio, Ts cte, 0 < θ ≤ θcr, H/L < 12 superficies Catton
UNIVERSIDAD TECNICA DE ORURO
FACULTAD NACIONAL DE INGENIERIA
CARRERA DE INGENIERIA MECANICA
ELECTROMECANICA
LABORATORIO DE MAQUINAS TERMICAS
TRANSFERENCIA MEC 2251
DE CALOR REGIMEN TRANSITORIO
Y CONVECCION
LABORATORIO 2 Pág. 10 de 34
T. media Ayaswamy y
17 ̅̅̅̅ ̅̅̅̅𝐿 (𝜃 = 90°)(𝑠𝑒𝑛𝜃)1/4
𝑁𝑢𝐿 =𝑁𝑢 Medio, Ts cte, θcr < θ < 90 superficies Catton
T. media
̅̅̅̅𝐿 = 1 + [𝑁𝑢
18 𝑁𝑢 ̅̅̅̅𝐿 (𝜃 = 90°) − 1]𝑠𝑒𝑛𝜃 Medio, Ts cte, 90 < θ ≤ 180 superficies Arnold
Recinto Cilindrico Concentrico
𝐷𝑒
𝑄 = 2𝜋𝑘𝑒𝑓 𝐿(𝑇𝑖 − 𝑇𝑒 )/ln( )
𝐷𝑖
Medio, Ts cte
𝑘 𝑅𝑎∗ 𝑃𝑟 1/4
= 0,386 [ ] 102 ≤ 𝑅𝑎 ∗ ≤ 107
𝑘𝑒𝑓 0,861 + 𝑃𝑟 δ = De – Di
𝐷𝑒 4
[ln( )]
𝑅𝑎 ∗ = 3 𝐷𝑖 𝑅𝑎
𝛿 [𝐷𝑒 −3/5 − 𝐷𝑖 −3/5 ]5 𝛿 T. media
19 superficies
Recinto esferico concéntrico
Correlacion 5
Tipo de flujo RaL C N
Laminar 104-107 0.54 1/4
Turbulento 107-1011 0.15 1/3
Correlacion 7
Ra D C n
10-10-10-2 0.675 0.058
10-2-102 1.020 0.148
102-104 0.850 0.188
104-107 0.480 1/4
107-1012 0.125 1/3
UNIVERSIDAD TECNICA DE ORURO
FACULTAD NACIONAL DE INGENIERIA
CARRERA DE INGENIERIA MECANICA
ELECTROMECANICA
LABORATORIO DE MAQUINAS TERMICAS
TRANSFERENCIA MEC 2251
DE CALOR REGIMEN TRANSITORIO
Y CONVECCION
LABORATORIO 2 Pág. 11 de 34
UNIVERSIDAD TECNICA DE ORURO
FACULTAD NACIONAL DE INGENIERIA
CARRERA DE INGENIERIA MECANICA
ELECTROMECANICA
LABORATORIO DE MAQUINAS TERMICAS
TRANSFERENCIA MEC 2251
DE CALOR REGIMEN TRANSITORIO
Y CONVECCION
LABORATORIO 2 Pág. 12 de 34
UNIVERSIDAD TECNICA DE ORURO
FACULTAD NACIONAL DE INGENIERIA
CARRERA DE INGENIERIA MECANICA
ELECTROMECANICA
LABORATORIO DE MAQUINAS TERMICAS
TRANSFERENCIA MEC 2251
DE CALOR REGIMEN TRANSITORIO
Y CONVECCION
LABORATORIO 2 Pág. 13 de 34
UNIVERSIDAD TECNICA DE ORURO
FACULTAD NACIONAL DE INGENIERIA
CARRERA DE INGENIERIA MECANICA
ELECTROMECANICA
LABORATORIO DE MAQUINAS TERMICAS
TRANSFERENCIA MEC 2251
DE CALOR REGIMEN TRANSITORIO
Y CONVECCION
LABORATORIO 2 Pág. 14 de 34
UNIVERSIDAD TECNICA DE ORURO
FACULTAD NACIONAL DE INGENIERIA
CARRERA DE INGENIERIA MECANICA
ELECTROMECANICA
LABORATORIO DE MAQUINAS TERMICAS
TRANSFERENCIA MEC 2251
DE CALOR REGIMEN TRANSITORIO
Y CONVECCION
LABORATORIO 2 Pág. 15 de 34
Referencia Bibliográfica
- Frank Incropera, David P. Dewitt, Fundamentos de Transferencia de Calor, Pearson,
Mexico, 1999
2. METODOLOGÍA
El trabajo de “Régimen Transitorio y Convección” fue realizado el día 20 de septiembre a las
14:30 p.m. hasta las 16:00 p.m.
El trabajo fue realizado por un grupo de 17 universitarios en los ambientes del laboratorio de
térmicas con el asesoramiento del Ing. Edgar Peñaranda M.
Modelo: ------------------------------
Marca : ------------------------------
Industria: ------------------------------
Color: Verde
Accesorios Manguera Fig. 2.2
Color: Amarillo
Fig. 2.3
Modelo: ------------------------------
Marca : ------------------------------
Industria: ------------------------------
Modelo: ------------------------------
Marca : ------------------------------
Industria: ------------------------------
Tipo: Paralelepípedo
Material: Aluminio
a = 32.8 [mm]
Dimensiones: b = 32.86 [mm]
c= 30.3 [mm]
Se ingresó al laboratorio llenando las fichas con los materiales y equipos a utilizar para el
préstamo
Al mismo tiempo se introduce el cuerpo solido (cubo) dentro del recipiente con agua, para
así lograr el equilibrio térmico del sólido y asegurar 88° C en toda la superficie exterior
Después de 10 minutos de ebullición del agua, este tiempo para asegurar la temperatura
deseada del cuerpo, retiramos el sólido para exponerlo al aire libre y obtener los datos
objeto de nuestro laboratorio.
De inmediato se toman los puntos de medición de temperatura en lados opuestos del cubo
y cada 5 minutos se toman mediciones de la temperatura reconociendo el proceso de
enfriamiento que tiene el sistema.
Para este cometido se utilizó un termómetro infrarrojo con el cual se debe tener mucho
cuidado con su manipulación, como por ejemplo el posicionamiento le termómetro de
manera perpendicular al punto de medición.
Por parte del ingeniero se realizó una explicación breve acerca del proceso analítico en el
laboratorio.
Para finalizar se tomaron medidas de los objetos medidos, en este caso el cubo para los
posteriores cálculos
Condiciones Ambientales:
Temperatura: T=18 ± 1 ºC
2.5. CALCULOS
Proponer, aplicar y generar resultados con el método analítico de régimen transitorio que
mejor corresponda a la prueba del enfriamiento del laboratorio y luego construir en una sola
gráfica el historial de temperaturas real y analíticas. Se construirán dos propuestas analíticas:
una con t inicial de 88.3 ºC y la otra con la t inicial medida
𝑇𝐴 +𝑇𝐵
N t[min] TA[°C] TB[°C] T[°C] ( )
2
1 0 62,4 60,2 𝟔𝟏, 𝟑
2 5 55,4 𝟓𝟓. 𝟒
3 10 49,1 𝟒𝟗. 𝟏
4 15 41,8 𝟒𝟏. 𝟖
5 20 36,5 36,2 𝟑𝟔. 𝟒
6 25 33,6 34,2 𝟑𝟑. 𝟗
7 30 31,1 31,4 𝟑𝟏. 𝟑
8 35 27,4 28,4 𝟐𝟕. 𝟗
ℎ∗𝐿
𝑁𝑢 =
𝐾
De donde:
𝑵𝒖 ∗ 𝑲𝒂𝒊𝒓𝒆
𝒉=
𝑳
Calculo del Numero de Nusselt
De la ecuación (1.7), C, n son constantes que dependen de la geometría del cuerpo, que
también nos permite escoger la expresión para calcular Ø.
Para una placa vertical (cubo) se tienen los siguientes datos: (Por Tablas)
1
−
1 1 4
𝑪 = 0.8 𝒏= ∅ = [1 + (1 + )]
4 √𝑃𝑟
𝑚2
Viscosidad cinemática del aire: 𝒗 = 1.702 ∗ 10−5
𝑠
Con estos valores se puede calcular el Número de Grashof (Ecuación 1.6), y el valor de Ø:
𝑔 ∗ 𝐿3 ∗ 𝛽 ∗ ∆𝑇
𝐺𝑟 =
𝑣2
𝑚
Donde g: Gravedad (9.76 𝑠2 𝑝𝑎𝑟𝑎 𝑂𝑟𝑢𝑟𝑜)
L: Longitud de una cara del cubo (se considerara la longitud mayor)
1
𝛽: Coeficiente de expansión volumétrica (𝑇 ) (Tf en K)
𝑓
∆𝑇: Diferencia de temperaturas entre la pared del cubo y el ambiente
v: Viscosidad cinemática del aire
UNIVERSIDAD TECNICA DE ORURO
FACULTAD NACIONAL DE INGENIERIA
CARRERA DE INGENIERIA MECANICA
ELECTROMECANICA
LABORATORIO DE MAQUINAS TERMICAS
TRANSFERENCIA MEC 2251
DE CALOR REGIMEN TRANSITORIO
Y CONVECCION
LABORATORIO 2 Pág. 24 de 34
Antes de calcular el Número de Nusselt se comprobara que el flujo es Laminar con la Ec (1.13)
𝑅𝑎 = 𝐺𝑟 ∗ Pr = 165562.83 ∗ 0.7255 = 𝟏. 𝟐 ∗ 𝟏𝟎𝟓
𝑵𝒖 = 𝟏𝟏. 𝟏𝟓𝟖
𝑉𝑜𝑙𝑢𝑚𝑒𝑛
𝑳𝑪 : Longitud característica del cuerpo (𝐴𝑟𝑒𝑎 (𝑡𝑜𝑡𝑎𝑙))
𝑊
K: Constante de conductividad del cuerpo (Aluminio (puro) = 237𝑚 ℃
UNIVERSIDAD TECNICA DE ORURO
FACULTAD NACIONAL DE INGENIERIA
CARRERA DE INGENIERIA MECANICA
ELECTROMECANICA
LABORATORIO DE MAQUINAS TERMICAS
TRANSFERENCIA MEC 2251
DE CALOR REGIMEN TRANSITORIO
Y CONVECCION
LABORATORIO 2 Pág. 25 de 34
𝑊
9.04 ∗ 5.47 ∗ 10−3 𝑚
𝑩𝒊 = 𝑚2 ℃ = 𝟐. 𝟎𝟖𝟔 ∗ 𝟏𝟎−𝟒
𝑊
237 𝑚 ℃
Como Bi < 0.1, la temperatura es solo función del tiempo y no del espacio.
𝑇(𝜃) − 𝑇∞
= 𝑒 −𝐵𝑖𝐹0
𝑇0 − 𝑇∞
𝑚2
En donde 𝛼: Difusividad térmica (Aluminio = 97.1 ∗ 10−6 )
𝑠
𝜃: Tiempo
𝐿𝐶 : Longitud característica del cuerpo
−𝟒 1
𝑇(𝜃) = (61.3 − 18)𝑒 −𝟐.𝟎𝟖𝟔∗𝟏𝟎 ∗𝟑.𝟐𝟒𝟓[ ] 𝜽
𝑠 + 18
Curva Analítica
Curva experimental
Fig. 2.12 Curva Experimental y curva Analítica de la Temperatura en función del tiempo
UNIVERSIDAD TECNICA DE ORURO
FACULTAD NACIONAL DE INGENIERIA
CARRERA DE INGENIERIA MECANICA
ELECTROMECANICA
LABORATORIO DE MAQUINAS TERMICAS
TRANSFERENCIA MEC 2251
DE CALOR REGIMEN TRANSITORIO
Y CONVECCION
LABORATORIO 2 Pág. 27 de 34
𝑇𝐴 +𝑇𝐵
N t[min] TA[°C] TB[°C] T[°C] ( )
2
1 0 88.3 88.3 𝟖𝟖, 𝟑
2 5 55,4 𝟓𝟓. 𝟒
3 10 49,1 𝟒𝟗. 𝟏
4 15 41,8 𝟒𝟏. 𝟖
5 20 36,5 36,2 𝟑𝟔. 𝟒
6 25 33,6 34,2 𝟑𝟑. 𝟗
7 30 31,1 31,4 𝟑𝟏. 𝟑
8 35 27,4 28,4 𝟐𝟕. 𝟗
ℎ∗𝐿
𝑁𝑢 =
𝐾
De donde:
𝑵𝒖 ∗ 𝑲𝒂𝒊𝒓𝒆
𝒉=
𝑳
Calculo del Numero de Nusselt
De la ecuación (1.7), C, n son constantes que dependen de la geometría del cuerpo, que
también nos permite escoger la expresión para calcular Ø.
Para una placa vertical (cubo) se tienen los siguientes datos: (Por Tablas)
1
−
1 1 4
𝑪 = 0.8 𝒏= ∅ = [1 + (1 + )]
4 √𝑃𝑟
𝑚2
Viscosidad cinemática del aire: 𝒗 = 1.798 ∗ 10−5 𝑠
Con estos valores se puede calcular el Número de Grashof (Ecuación 1.6), y el valor de Ø:
𝑔 ∗ 𝐿3 ∗ 𝛽 ∗ ∆𝑇
𝐺𝑟 =
𝑣2
𝑚
Donde g: Gravedad (9.76 2 𝑝𝑎𝑟𝑎 𝑂𝑟𝑢𝑟𝑜)
𝑠
L: Longitud de una cara del cubo (se considerara la longitud mayor)
1
𝛽: Coeficiente de expansión volumétrica ( ) (Tf en K)
𝑇𝑓
∆𝑇: Diferencia de temperaturas entre la pared del cubo y el ambiente
v: Viscosidad cinemática del aire
Antes de calcular el Número de Nusselt se comprobara que el flujo es Laminar con la Ec (1.13)
𝑅𝑎 = 𝐺𝑟 ∗ Pr = 230787.016 ∗ 0.7228 = 𝟏. 𝟔𝟕 ∗ 𝟏𝟎𝟓
𝑵𝒖 = 𝟏𝟐. 𝟏𝟏
𝑉𝑜𝑙𝑢𝑚𝑒𝑛
𝑳𝑪 : Longitud característica del cuerpo (𝐴𝑟𝑒𝑎 (𝑡𝑜𝑡𝑎𝑙))
𝑊
K: Constante de conductividad del cuerpo (Aluminio (puro) = 237𝑚 ℃
𝑊
10.8 ∗ 5.47 ∗ 10−3 𝑚
𝑩𝒊 = 𝑚2 ℃ = 𝟐. 𝟒𝟗𝟑 ∗ 𝟏𝟎−𝟒
𝑊
237 𝑚 ℃
Como Bi < 0.1, la temperatura es solo función del tiempo y no del espacio.
𝑇(𝜃) − 𝑇∞
= 𝑒 −𝐵𝑖𝐹0
𝑇0 − 𝑇∞
𝛼∗𝜃
𝐹𝑜𝑢𝑟𝑖𝑒𝑟 𝐹𝑜 =
𝐿𝑐 2
𝑚2
En donde 𝛼: Difusividad térmica (Aluminio = 97.1 ∗ 10−6 𝑠
)
𝜃: Tiempo
𝐿𝐶 : Longitud característica del cuerpo
UNIVERSIDAD TECNICA DE ORURO
FACULTAD NACIONAL DE INGENIERIA
CARRERA DE INGENIERIA MECANICA
ELECTROMECANICA
LABORATORIO DE MAQUINAS TERMICAS
TRANSFERENCIA MEC 2251
DE CALOR REGIMEN TRANSITORIO
Y CONVECCION
LABORATORIO 2 Pág. 30 de 34
−𝟒 1
𝑇(𝜃) = (88.3 − 18)𝑒 −𝟐.𝟒𝟗𝟑∗𝟏𝟎 ∗𝟑.𝟐𝟒𝟓[ ] 𝜽
𝑠 + 18
Curva Analítica
Curva experimental
Fig. 2.12 Curva Experimental y curva Analítica de la Temperatura en función del tiempo