Transferencia de calor

Asignación 5. Revisión completa del ejercicio adicional 1er parcial

Cajigas Velásquez Grecia
Viernes 11 de noviembre 2022
Objetivo: Reafirmar los conceptos aprendidos sobre transferencia de calor en estado transitorio.
Se va a realizar un experimento con el fin de determinar el coeficiente de transferencia de calor sobre las superficies de
tomates que se colocan en agua fría a 7°C. Los tomates (k = 0.59 W/m °C, alfa = 0.141 x 10 -6 m2/s, densidad = 999 kg/m3,
cp = 3.99 kJ/kg °C) con una temperatura inicial uniforme de 30°C tienen forma esférica con un diámetro de 8 cm. Después
de un periodo de 2 horas, las temperaturas en el centro y en la superficie de los tomates son de 10.0°C y 7.1°C,
respectivamente.
a) Aplicando el método analítico de aproximación de un término (no el de las gráficas de Heisler), calcula el
coeficiente de transferencia de calor.

𝛼𝑡 (0.141𝑥10−6 𝑚 2 /𝑠)(7200 𝑠)
𝜏= = = 0.6345
𝑟2 (0.04 𝑚)2
𝜏 = 0.6345 > 0.2 Aplicar aproximación con u valor de 𝜆𝑛 𝑦 𝐴𝑛
𝑇 − 𝑇∞
𝜃=
𝑇𝑖 − 𝑇∞
2 𝑠𝑒𝑛 𝜆1
𝜃 = 𝐴1 𝑒 −𝜆1 𝜏
𝜆1
𝜃 𝑠𝑒𝑛 𝜆1
= = 0.0333
𝜃0 𝜆1
Graficando las funciones 0.0333 𝜆1 = 𝑠𝑒𝑛 𝜆1 obtenemos el valor para 𝜆1 ≈ 3.03, y eso nos da un Biot de 30, por lo
tanto:
𝐵𝑖 𝑘
ℎ=
𝑟
h= 442.5 w/𝑚 2 °𝐶
 b) Asumiendo asimismo la aproximación de un término, deriva las expresiones para el flux de calor en la superficie
de los tomates. Considera ambas opciones para evaluar el flux: (1) el gradiente de temperatura en el tomate y (2)
la convección al fluido. Evalúa para t = 0, 20, 40 y 60 minutos. Grafica el flux (t) y estima en forma gráfica la
energía cedida por un tomate en 60 min. Confirma esta estimación con los cálculos analíticos.
𝑞𝑐𝑜𝑛𝑣 = ℎ (𝑇(𝑥, 𝑡) − 𝑇∞ ) Flux por convección
(𝑇∞ −𝑇(𝑥,𝑡))
𝑞𝑐𝑜𝑛𝑑 = − 𝑘 𝑟
Flux por convección

𝑇(𝑥, 𝑡) = 𝑇∞ + 𝜃(𝑇𝑖 − 𝑇∞ ) Hay una diferencia de temperaturas negativa, porque

por la esencia de una delta, es estado final menos
estado inicial, y tiene sentido porque nuestro sistema
se enfría, el producto final es positivo por la constante
k negativa

El área en la superficie del tomate es constante 4𝜋𝑟 2 = 4𝜋(0.04 𝑚 )2 = 0.02 𝑚 2

t(s) 𝜏 𝜃 T(°C) 𝑞𝑐𝑜𝑛𝑣 𝑞𝑐𝑜𝑛𝑣 ∗ á𝑟𝑒𝑎 Δ𝑡
0 0 1 30 10177.5 203.55 0
1200 0.10575 0.0130866 7.30099 133.188 2.6637 1200
2400 0.2115 4.9337x10^-3 7.1134 50.1795 1.003 1200
3600 0.31725 1.86003x10^-3 7.0427 18.8947 0.37789 1200

t(s) 𝜏 𝜃 T(°C) 𝑞𝑐𝑜𝑛𝑑 𝑞𝑐𝑜𝑛𝑑 ∗ á𝑟𝑒𝑎 Δ𝑡

0 0 1 30 339.25 6.785 0
1200 0.10575 0.0130866 7.30099 4.439 0.08878 1200
2400 0.2115 4.9337x10^-3 7.1134 1.6726 0.03345 1200
3600 0.31725 1.86003x10^-3 7.0427 0.629 0.01258 1200

q (Energía/área)
Flux por conveción
12000

10000

8000

6000

4000

2000

0
0 1200 2400 3600

Flux conv Tiempo (s)

Gráfica del Flux por convección

con respecto al tiempo
 q (Energía/área) Flux por conducción
400

350

300

250

200

150

100

50

0
0 1200 2400 3600
Tiempo (s)
Flux cond

Gráfica del Flux por conducción con

respecto al tiempo

Q conv = 4683.6 J
Q cond= 147 J

Y de forma analítica tenemos la siguiente expresión: 𝑄 = á𝑟𝑒𝑎 ∫ 𝑞(𝑡) ⅆ𝑡

𝑄 = á𝑟𝑒𝑎∫ ℎ (𝑇(𝑥, 𝑡) − 𝑇∞ ) ⅆ𝑡 Por convección

(𝑇∞ −𝑇(𝑥,𝑡))
𝑄 = á𝑟𝑒𝑎∫ − 𝑘 𝑟
ⅆ𝑡 Por conducción
𝛼𝑡
Localizando t de 𝜏 = 𝑟2 que está en el término de 𝜃 adimensional, del que depende T (x, t) y sacando todo lo que es
constante de la integral tenemos:
𝐴1 𝑠𝑒𝑛 𝜆1 2 𝛼𝑡
𝑄 = á𝑟𝑒𝑎 ℎ (𝑇𝑖 − 𝑇∞ ) ∫ 𝑒 −𝜆1 𝑟2 ⅆ𝑡
𝜆1

(𝑇𝑖 −𝑇∞ ) 𝐴1 𝑠𝑒𝑛 𝜆1 2 𝛼𝑡

𝑄 = á𝑟𝑒𝑎 ℎ ∫ 𝑒 −𝜆1 𝑟2 ⅆ𝑡
𝑟 𝜆1

Y resolviendo la integral:
𝐴1 𝑠𝑒𝑛 𝜆1 −𝑟2 2 𝛼𝑡
𝑄𝑐𝑜𝑛𝑣 = á𝑟𝑒𝑎 ℎ (𝑇𝑖 − 𝑇∞ ) 2 𝑒 −𝜆1 𝑟2
𝜆1 𝜆1 𝛼

𝛼𝑡
(𝑇𝑖 −𝑇∞ ) 𝐴1 𝑠𝑒𝑛 𝜆1 −𝑟2 −𝜆1 2 2
𝑄𝑐𝑜𝑛ⅆ = á𝑟𝑒𝑎 ℎ 𝑟 𝜆1 2 𝑒 𝑟
𝜆1 𝛼

Q conv = 4683.6 J
Q cond= 147 J
 c) Finalmente evalúa con otra perspectiva la energía absorbida por un tomate en 1 hora (a partir de la energía
acumulada en cada capa de un tomate: considera la capacidad calorífica y el aumento de temperatura dado por el
perfil T(r) en el tomate). Realiza el cálculo numérico y el analítico.
Discusión: El tomate no absorbe energía, es un intercambio de calor y el tomate se enfría, por lo tanto, la pierde.
𝑄 = 𝜌𝑉𝐶𝜌 (𝑇(𝑟, 𝑡) − 𝑇∞ )

4
𝑄 = 𝜌 𝜋𝑟 3 𝐶𝜌 (𝑇(𝑟, 𝑡) − 𝑇∞ )
3

De forma numérica:
Para t= 1 hora= 3600s 𝜏 𝑒𝑠 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒, 𝜏 = 0.31725
Considerando un Δ𝑟 = 0.02𝑚

Radio (m) 𝜃 T (°C) Q

0 0.10786 9.4878 0
0.02 1.8823x10^-3 7.04329 5.78x10^-3
0.04 1.86003x10^-3 7.04278 0.04571

Q=0.05149 kJ = 51.492 J