Universidad Tecnológica de Panamá
Facultad de Ingeniería Eléctrica
Ingeniería Electromecánica
Conducción Radial: Dependencia con la potencia de calentamiento y del flujo de
refrigeración
Nombre: Amaya, Joshua
Cédula: 8-903-85
Nombre: De León, Ovidio
Cédula: 8-904-976
Nombre: Duarte, Moisés
Cédula: 3-736-1999
Nombre: Thomas, Gabriela
Cédula: 8-888-965
Descriptores: conductividad térmica, resistencia térmica, termopar.
Conductividad térmica: Capacidad de una sustancia de transferir la energía en forma de calor de sus moléculas a otras
adyacentes o a sustancias con las que está en contacto.
Resistencia térmica: Representa la capacidad del material de oponerse al flujo del calor. En el caso de materiales
homogéneos es la razón entre el espesor y la conductividad térmica del material; en materiales no homogéneos la
resistencia es el inverso de la conductancia térmica.
Termopar: sensor para medir la temperatura. Se compone de dos metales diferentes, unidos en un extremo. Cuando la
unión de los dos metales se calienta o enfría, se produce una tensión que es proporcional a la temperatura.
INTRODUCCIÓN
La conducción (transferencia de calor por difusión) se refiere al transporte de energía en un medio debido a un gradiente de
temperatura, y el mecanismo físico es el de la actividad aleatoria atómica o molécula. La capacidad para conducir calor es
una propiedad de cada material, y esta se conoce como la conductividad térmica; en esta práctica se realizan dos tipos de
ensayos para determinar esta propiedad.
Si suponemos la conducción estacionaria de calor a través de un tubo de agua caliente. El calor se pierde de forma continua
hacia el exterior a través de la pared del tubo, la transferencia de calor a través de este se efectúa de forma normal a la
superficie y no se tiene alguna transferencia significativa en otras direcciones, por lo tanto, la transferencia de calor puede
considerarse unidimensional, puesto que la temperatura del tubo depende de una sola dirección (radial). Y la temperatura
puede expresarse como T=T(r).
Para la realización de esta experiencia de laboratorio se utilizó el equipo TXC-CR de EDIBON. El objetivo de esta práctica
es demostrar que la distribución de calor se realiza de forma radial y observar cómo se afecta la distribución de temperatura
al variar la potencia de entrada y el caudal del flujo de refrigeración. Presentaremos un ensayo por conducción radial el cual
está constituido por un disco, provisto de una unidad de calentamiento eléctrico central y una unidad de enfriamiento
periférico; desde el centro a la periferia se cuentan con seis sensores de temperatura colocados a 10mm de distancia cada uno.
OBJETIVO GENERAL:
Estudiar la conducción de calor radial observando su dependencia de la potencia de entrada.
OBJETIVOS ESPECÍFICOS:
1. Demostrar que la distribución de calor tiene lugar en la superficie plana.
2. Analizar la distribución de calor para diferentes potencias de entrada.
FÓRMULAS DE TRABAJO
Para calcular el flujo de calor en una superficie que se transmite de forma radial se utiliza la siguiente ecuación:
𝑄=2𝜋(𝑇𝑖−𝑇𝑗) / ln(𝑟𝑗/𝑟𝑖)
Dónde:
Universidad Tecnológica de Panamá
Facultad de Ingeniería Eléctrica
Ingeniería Electromecánica
• 𝑘 es la conductividad térmica del material.
• 𝑙 es la longitud del cilindro, en este caso es de 3 mm.
• 𝑇𝑗 y 𝑇𝑖 son las temperaturas en el punto j e i respectivamente.
• 𝑟𝑗 y 𝑟𝑖 son los radios en donde se encuentran 𝑇𝑗 y 𝑇𝑖 ubicados.
EQUIPAMIENTO
• Sistema de adquisición y control SCADA.
• Equipo TXC-CR
RECOLECCIÓN DE DATOS
Parte A: Dependencia con la potencia de calentamiento.
Radio(mm) Temp.(C) Ln(r)
13 46.8 2.56
21 36.7 3.04
29 36.1 3.36
37 33.2 3.61
45 30.4 3.8
53 28.3 3.97
RESULTADOS
Hacer un gráfico T (r) vs r y ajustar a una ecuación cuadrática.
Temp.(C)
50
45
40
35
30 y = 0.0087x2 - 0.9828x + 56.576
R² = 0.9339
25
20
15
10
5
0
0 10 20 30 40 50 60
Gráfico logarítmico lineal de la temperatura (lineal) en función del radio (logarítmico).
Universidad Tecnológica de Panamá
Facultad de Ingeniería Eléctrica
Ingeniería Electromecánica
Temp vs Ln(r)
4.5
3.5
3
y = -0.0782x + 6.1452
2.5
R² = 0.9498
2
1.5
0.5
0
0 10 20 30 40 50
Calcular la curva teórica; asumiendo que se conocen las siguientes condiciones de frontera:
R1 = T1 y R6= T6
Q = 42W
K = 155.59 W/ (m ºC).
T(r) = C1 ln(r)+C2
→ Sabiendo que 𝑇(𝑟1 ) = 𝑇1 y 𝑇(𝑟6 ) = 𝑇6 , entonces:
𝑇(𝑟1 ) = 𝑇1 = 𝐶1 ln(𝑟1 ) + 𝐶2
𝑇(𝑟6 ) = 𝑇6 = 𝐶1 ln(𝑟6 ) + 𝐶2
→ Luego:
𝑇(𝑟6 ) − 𝑇(𝑟1 ) = 𝑇6 − 𝑇1 = 𝐶1 ln(𝑟6 ) + 𝐶2 − 𝐶1 ln(𝑟1 ) − 𝐶2
𝑟6
𝑇6 − 𝑇1 = 𝐶1 ln ( )
𝑟1
𝑇6 − 𝑇1
𝐶1 = 𝑟
ln ( 6 )
𝑟1
𝑇6 −𝑇1
→ Resultando: 𝑇(𝑟) = ( 𝑟 ) ln(𝑟) + 𝐶2
ln( 6 )
𝑟1
Para hallar C2:
𝑇6 − 𝑇1
𝑇(𝑟1 ) = 𝑇1 = ( 𝑟 ) ln(𝑟1 ) + 𝐶2
ln ( 6 )
𝑟1
𝑇6 − 𝑇1
𝐶2 = 𝑇1 − ( 𝑟 ) ln(𝑟1 )
ln ( 6 )
𝑟1
→ La ecuación de temperatura queda así:
Universidad Tecnológica de Panamá
Facultad de Ingeniería Eléctrica
Ingeniería Electromecánica
𝑇6 − 𝑇1 𝑇6 − 𝑇1
𝑇(𝑟) = ( 𝑟6 ) ln(𝑟) + 𝑇1 − ( 𝑟 ) ln(𝑟1 )
ln ( ) ln ( 6 )
𝑟1 𝑟1
→ Ahora podemos determinar las temperaturas de esta experiencia de laboratorio, mediante la función anterior:
28.3 − 46.8 28.3 − 46.8
(𝑟1 ) = ( ) ln(13) + 46.8 − ( ) ln(13) = 46.8º𝐶
53 53
ln ( ) ln ( )
13 13
28.3 − 46.8 28.3 − 46.8
𝑇(𝑟2 ) = ( ) ln(121) + 46.8 − ( ) ln(13) = 40.48º𝐶
53 53
ln ( ) ln ( )
13 13
28.3 − 46.8 28.3 − 46.8
𝑇(𝑟3 ) = ( ) ln(29) + 46.8 − ( ) ln(13) = 36.24º𝐶
53 53
ln ( ) ln ( )
13 13
28.3 − 46.8 28.3 − 46.8
𝑇(𝑟4 ) = ( ) ln(37) + 46.8 − ( ) ln(13) = 33.03º𝐶
53 53
ln ( ) ln ( )
13 13
28.3 − 46.8 28.3 − 46.8
𝑇(𝑟5 ) = ( ) ln(45) + 46.8 − ( ) ln(13) = 30.45º𝐶
53 53
ln ( ) ln ( )
13 13
28.3 − 46.8 28.3 − 46.8
𝑇(𝑟6 ) = ( ) ln(53) + 46.8 − ( ) ln(13) = 28.3º𝐶
53 53
ln ( ) ln ( )
13 13
𝑄 = 𝑚̇𝑐𝑝 ∆𝑇
𝑔 𝐽
𝑄 = (42.55 ) (4.1813 ) (28.3 − 27) = 231.28 𝑊
𝑠 𝑔𝐾
Calcular la razón de calor utilizándola ecuaciones teóricas y la ecuación cuadrática. Y compara los resultados.
→ Razón de calor experimental:
Qentrada (W) QCalculada (W) Error ( % )
42 44.06 4.67
→ Razón de calor teórica: calculamos el calor promedio.
2𝜋𝑘𝐿(𝑇1 − 𝑇2 )
𝑄̇ = 𝑟
ln ( 2 )
𝑟1
2𝜋(3𝑥10−3 )(155.59)(46.8 − 36.7)
𝑄12̇ = = 61.77 𝑊
21
ln ( )
13
2𝜋(3𝑥10−3 )(155.59)(36.7 − 36.1)
̇
𝑄23 = = 5.45 𝑊
29
ln ( )
21
Universidad Tecnológica de Panamá
Facultad de Ingeniería Eléctrica
Ingeniería Electromecánica
2𝜋(3𝑥10−3 )(155.59)(36.1 − 33.2)
̇ =
𝑄34 = 34.91 𝑊
37
ln ( )
29
2𝜋(3𝑥10−3 )(155.59)(33.2 − 30.4)
̇ =
𝑄45 = 41.95 𝑊
45
ln ( )
37
2𝜋(3𝑥10−3 )(155.59)(30.4 − 28.3)
𝑄56̇ = = 37.64 𝑊
53
ln ( )
45
𝑆𝑢𝑚𝑎 𝑄
𝑄̇𝑃𝑅𝑂𝑀 = ( 𝑁 )
61.77 𝑊+ 5.45 𝑊+ 34.91 𝑊+41.95 𝑊 +37.64 𝑊
𝑄̇𝑃𝑅𝑂𝑀 = ( )
5
𝑸̇𝑷𝑹𝑶𝑴 = 𝟑𝟔. 𝟑𝟒𝟒 𝑾
Durante la experiencia de laboratorio pudimos notar que el termopar 2 estaba en malas condiciones, por lo que la
𝑸̇𝑷𝑹𝑶𝑴 calculada nos da un valor fuera del rango. Si se calcula nuevamente la 𝑸̇𝑷𝑹𝑶𝑴 sin tomar en cuenta la medición
del termopar 2, nos da un valor más cercano al rango.
61.77 𝑊+ 34.91 𝑊+41.95 𝑊 +37.64 𝑊
𝑄̇𝑃𝑅𝑂𝑀 = ( ) = 44.06 W
4
Calcular la k del material y comparar con el material real.
𝑟
ln ( 2 )
𝑟1
𝑘 = 𝑄̇
2𝜋𝐿(𝑇1 − 𝑇2 )
21
ln (
) 𝑊
𝑘12 = (42𝑊) 13 = 106
−3
2𝜋(3𝑥10 𝑚)((46.8 − 36.7) 𝑚 º𝐶
29
ln(21) 𝑊
𝑘23 = (42𝑊) = 1198.6
2𝜋(3𝑥10−3 𝑚)((36.7−36.1) 𝑚 º𝐶
37
ln (
) 𝑊
𝑘34 = (42𝑊) 29 = 187.9
2𝜋(3𝑥10−3 𝑚)((36.1 − 33.2) 𝑚 º𝐶
45
ln ( ) 𝑊
𝑘45 = (42𝑊) 37 = 154
−3
2𝜋(3𝑥10 𝑚)((33.2 − 30.4) 𝑚 º𝐶
53
ln ( ) 𝑊
𝑘56 = (42𝑊) 45 = 173.6
−3
2𝜋(3𝑥10 𝑚)((30.4 − 28.3) 𝑚 º𝐶
𝑆𝑢𝑚𝑎 𝑘
𝑘𝑃𝑅𝑂𝑀 = ( )
𝑁
→ Para este cálculo decidimos obviar la k del termopar 2, ya que se pudo demostrar en el cálculo de Qprom que
esta medición es errada, por las malas condiciones del termopar.
106+187.9 +154 +173.6
𝑘𝑃𝑅𝑂𝑀 = ( )
4
𝑾
𝒌𝑷𝑹𝑶𝑴 = 𝟏𝟓𝟓. 𝟑𝟕𝟓
𝒎 º𝑪
Universidad Tecnológica de Panamá
Facultad de Ingeniería Eléctrica
Ingeniería Electromecánica
Podemos decir que la 𝒌𝑷𝑹𝑶𝑴 calculada se encuentra en el rango
de los metales.
CONCLUSIONES
En este experimento, observando que el valor calculado del flujo de calor tomando como referencia dos
termopares no resultó constante (o muy parecidos entre sí) para cada par de puntos, el valor de 𝑄𝑐𝑎𝑙𝑐𝑢𝑙𝑎𝑑𝑜 es el
promedio de estos valores.
Al realizar la práctica se tiene como un factor importante la medida de la temperatura, la cual en algunos puntos
tuvo error, ya que los termopares que son los elementos para medirla arrojaban valores que por inspección no cabían
en el rango de medición. Además, la variación del caudal afecta en la medida de disipación de calor y este no era
muy constante durante la toma de los datos. Los cálculos presentan errores muy altos debido a que las condiciones
del medio cambiaban continuamente, ya que el tiempo de estabilización del sistema no era el apropiado.
Muy generalmente se puede decir que mientras más cercanos a la fuente de calor estén los termopares, más
cercano al valor real será el valor de Q.
Para los cálculos de Q y k experimentales obviamos los datos tomados por el termopar 2 ya que parecían
inconsistentes con los demás datos, lo que se puede corroborar con la gráfica, el punto del termopar 2 muestra un
cambio drástico en la función.
Gabriela Thomas
REFERENCIAS BIBLIOGRÁFICAS
Cengel, Y., 2007. Transferencia de calor y masa. México DF: McGraw-Hill. EDIBON, 2013. TXC/CR manual de
prácticas.