Gestión Energética II

Guı́a de ejercicios Certamen N°4

Tercer trimestre 2021

1 Calentamiento de agua mediante colectores solares

En la casa de su amigo poseen un estanque en donde se acumula agua calentada a través de un
colector solar. El estanque es capaz de almacenar m kg de agua, la cual tiene un calor especifico cp . Se
sabe que el colector calienta agua gracias a una potencia neta recibida debido a un nivel de irradiancia
It .
Según el manual del colector, las especificaciones técnicas de éste son:

• Coeficiente de transmitividad del vidrio en la cubierta τ .

• Coeficiente de absortividad de la placa α.

• Área efectiva del colector Ac .

• Resistencia térmica total del colector solar Rt .

A su amigo le interesa conocer cómo varı́a la temperatura mientras pasa el tiempo y le pide a usted
que realice los cálculos, pero él asegura que el colector no pierde energı́a en el tiempo y le obliga a
proceder considerando este supuesto. Además, le dice que la irradiancia varı́a con respecto al tiempo
de forma sinusoidal. Al ser medida obtuvo que en un ciclo de 24 horas la máxima irradiancia se registra
a las 12 del dı́a, mientras que a las 12 de la noche no registra irradiancia. El colector solar fue puesto
en marcha exactamente a las 6 de la mañana (t=0), y a esa hora fue medida la temperatura del fluido
la cual era de Tm . Las especificaciones del colector además dicen que la eficiencia de transferencia de
calor de la placa al fluido es igual a np .
Para el siguiente análisis es correcto asumir que los parámetros (a excepción de la irradiancia) no
varı́an con el tiempo.
1. Dados los supuestos del problema, muestre que el perfil de temperatura T (t) es:
   
np τ αAc H 12 πt
T (t) = − cos + t + C1
mcp π 12

Donde H es la amplitud de la función sinusoidal.

2. Dado lo obtenido en el punto anterior y las condiciones de borde muestre que:

12np τ αAc H
C1 = Tm +
πmcp

3. ¿Son válidos los supuestos que le obligó a usar su amigo? ¿Qué consecuencia tiene el supuesto de
no pérdida de energı́a al querer determinar la temperatura máxima?
Solución
Supuestos
• El colector no cuenta con pérdidas a través de ningún mecanismo de transferencia.
• Radiación solar varı́a en el tiempo de forma sinusoidal.
• La radiación incide perpendicular al colector.
• Propiedades del agua y colector constantes en el tiempo.
• Las paredes del estanque son adiabáticas.
1. Como se tiene un estanque con masa de agua, se puede realizar un balance de energı́a entre
potencia efectiva que traspasa el colector al fluido y la ganancia neta de calor del fluido que se
almacena en el estanque, lo cual se puede expresar de la siguiente manera:
dTf
ṁcp = np τ αAc IT − np Ac Ul (Tf − Ta )
dt
Pero si consideramos las perdidas despreciables:
dTf
ṁcp = np τ αAc IT
dt
Se debe obtener una expresión para la irradiancia en función del tiempo, para lo cual se sabe:
1) se comporta como una función seno Asen(ax + b)
2) Se mide en un ciclo de 24 horas, con saltos de 1 hora.
3) A las 6 A.M (t=0), la irradancia se encuentra a la mitad de su valor máximo.
De acuerdo a lo señalado, la variación de la radiación puede expresarse como:
2πt
IT (t) = H(sen( ) + 1)
24
Reemplazando y despejando:
np τ αAc H πt
dTf = (sen( ) + 1)dt
ṁcp 12
Integrando:
np τ αAc H 12 πt
Tf (t) = (− cos( ) + t) + C1
ṁcp π 12
2. Utilizando la condición de borde:
12np τ αAc H
Tf (0) = Tm = − + C1
π ṁcp
12np τ αAc H
Tm + = C1
π ṁcp
3. El supuesto que considera el amigo es erróneo, puesto que no es posible que no existan perdidas
de calor hacia el ambiente desde el colector solar. Al querer calcular la temperatura máxima se
generan problemas debido a que:
lı́m Tf (t) = ∞
t→∞
Lo que significarı́a que la temperatura puede aumentar de manera infinita en el colector,y que el
colector puede almacenar energı́a ı́limitada lo que claramente no es verdadero.

Taller C4 - GE II Página 2
2 Análisis de un colector solar térmico.
Un colector de placa plana de transmisividad del vidrio τ = 0, 87, absortividad de la placa del
colector α = 0, 91. A través de distintos estudios se pudo determinar que el porcentaje del calor que es
traspasado desde la placa colectora al fluidoes 70 %. Se sabe también que su coeficiente de pérdidas
W
globales toma un valor de UL = 4.9 .
m2 · ℃ W 
Considere el siguiente modelo para la irradiancia en m 2 , en función de la hora del dı́a t (las 00:00
corresponde a t = 0, las 01:00 a t = 1, las 13:00 a t = 13 . . . ):

IT (t) = −3, 1468 · t2 + 75, 62 · t (1)

Con esto en mente se le solicita:

a) A través del modelo de Hottel-Whillier, encontrar la eficiencia instantánea ηi a las 14:00

(t = 14). Se sabe que la temperatura de entrada al colector a esa hora es de Te = 25◦ C y la
temperatura ambiente es Tamb = 15◦ C.

b) Graficar la eficiencia instantánea, tomando como eje ordenado la eficiencia instantánea ηi y en

el eje de las abscisas la variable TeI−T amb
T (t)
, indicando claramente el valor de la pendiente y el valor
de intercepto con el eje ordenado.

Taller C4 - GE II Página 3
Solución
Primeramente, se debe encontrar la irradiancia que existe en esa zona a esta hora. Para esto se reem-
plazara en el modelo de irradiancia.

I(14) = −3, 1468 · 142 + 75.62 · 14

I(14) = 441, 9 ≈ 442
Con esto en mente se puede conocer el valor de la eficiencia instantanea segun el modelo de Hottel
Whiller

(Te − Tamb )
ηi = ηpf · τ α − ηpf · UL ·
I(t)
Como ya se conocen todos los datos se procede a encontrar la eficiencia instantanea.

(25 − 15)
ηi (t = 14) = 0, 7 · 0, 87 · 0, 91 − 0, 7 · 4, 788 ·
442
ηi (t = 14) = 0, 47

Pregunta B
Dado que se conocen los parámetros que caracterizan a este colector, se puede plantear el modelo de
Hottel Whiller de la siguiente manera.
(Te − Tamb )
ηi = 0, 5542 − 3, 35
I(t)
Se puede apreciar que esto es una ecuación de la recta de la forma.
y = 0, 5542 − 3, 35 · x

Figura 1: Gráfico de eficiencia instantánea del colector

Taller C4 - GE II Página 4
3 Intercambiador de Calor concéntricos a contraflujo
Se utiliza un intercambiador de calor de tubo concéntrico a contraflujo para enfriar el aceite lubricante
de un gran motor de turbina a gas industrial. El caudal de agua utilizado para la refrigeración que
pasa a través del tubo interior (Di= 25 mm) es de 0,2 kg/s, por otro lado, el flujo de aceite caliente a
través del anillo exterior (Do = 45 mm) es 0,1 kg/s. El coeficiente de convección asociado con el flujo
de aceite es ho =40W/m2 K.

El aceite y el agua entran a temperaturas de 100 y 30° C, respectivamente, tal como se muestran
en el esquema. Para el caso del agua, dada la curva en la figura, se puede asumir que las propiedades
óptimas se encuentran cuando el fluido está a 35°C.

Figura 2: Esquema Intercambiador de Calor a Contraflujo

En base a los datos mencionados, responda.

1. ¿Cuál es la temperatura de salida del agua?

2. ¿Cuál es la longitud de tubo requerida para una temperatura de salida de aceite de 60 C?

Hint: Para calcular el número de Nusselt local de flujo turbulento en un tubo circular liso se utiliza
la correlación de Dittus- Boelter, que se muestra a continuación. Donde n es igual a 0,4 para calentar
(Ts ¿Tm) y 0.3 para enfriamiento (Ts ¡Tm).
 
0, 6 ≤ P r ≤ 160
4/5
N uD = 0, 023ReD P rn  ReD ≥ 10.000 
L/D ≥ 10

Taller C4 - GE II Página 5
 Figura 3: Propiedades termo fı́sicas del aceite saturado

Figura 4: Propiedades termo fı́sicas del agua saturada

Taller C4 - GE II Página 6
3.1 Solución:
Supuestos:
• Pérdida de calor insignificante con el entorno.

• Efectos de la energı́a cinética y potencial insignificantes.

• Propiedades constantes.

• Resistencia térmica de la pared del tubo y factores de incrustación insignificantes.

• Condiciones completamente desarrolladas para el flujo de agua.

En primer lugar se calculan las propiedades del agua a Tc = 35°C = 308 K, de lo que se obtiene:
• Cp = 4.178 J/kg·K

• µ = 725 x 10−6 N·/m2

• k = 0,625 W/m·K

• Pr = 4,85
Por otra parte, para obtener las propiedades del aceite se toma la temperatura promedio de las
paredes del tubo, la que es T¯h =80°C=353 K, en donde el Cp del aceite es igual a 2.131 J/kg·K.

Con los datos presentados, se puede calcular la tasa de transferencia de calor mediante un balance
de energı́a del fluido caliente (aceite).

q = ṁh Cp,h (Th,i − Th,o ) = 0, 1kg/s · 2.131J/kg · K(100 − 60)C = 8.524W

Tomando en consideración que el calor que entrega el fluido caliente es el mismo que absorbe el
frı́o, se hace un balance de energı́a a este último, en donde la temperatura de salida del agua es:
q 8.524W
Tc,o = + Tc,i = + 30C = 40, 2C
ṁc Cp,c 0, 2kg/s · 4.178J/kg · K
La longitud requerida para que cumplir con la temperatura solicitada se puede obtener a partir de
la ecuación de velocidad de convección:
q = U A∆Tlm
Dada la geometrı́a, el área está dada por A=πDi L.
Como ya son conocidas todas las temperaturas, se puede calcular la diferencia de temperatura
media logarı́tmica.

(Th,i − Tc,o ) − (Th,o − Tc,i ) 59, 8 − 30

∆Tlm = = = 43, 2C
ln[(Th,i − Tc,o )/(Th,o − Tc,i )] ln(59, 8/30)
Recordando que el coeficiente de transferencia de calor puede ser expresado como:
1
= Req
UA
para este caso solo se considera la convección, por lo que el coeficiente queda de la siguiente manera:
1
U=
(1/hi ) − (1ho )

Taller C4 - GE II Página 7
 Para estimar el coeficiente convectivo del fluido frı́o (hi ) se calcula el número de Reynolds, en pos
de caracterizar el flujo y seleccionar la correlación adecuada.

4ṁc 4 · 0, 2kg/s
ReD = = = 14.050
πDi µ π(0, 025m)(725 · 106 N · s/m2 )
Dado el valor de Reynolds obtenido, el flujo se define como turbulento y el coeficiente de convección
puede estimarse utilizando la correlación Dittus-Boelter, con un n = 0,4 puesto que Ts ¿Tm.
4/5
N uD = 0, 023ReD P r0,4 = 0, 023(14.050)4/5 (4, 85)0,4 = 90

k 90 · 0, 625W/m · K
hi = N uD = = 2.250W/m2 · K
Di 0, 025m
como el coeficiente de convección para el fluido caliente es ho =40W/m2 K, el coeficiente general
es:
1
U= = 39, 3W/m2 · K
(1/2250W/m2 · K) + (1/40W/m2 · K)
Por último, de la ecuación de la tasa de convección se deduce que la longitud requerida del inter-
cambiador es:

q 8.5224W
L= = 2
= 63, 9m
U πDi ∆Tlm 39, 3W/m · Kπ(0, 025m)(43, 2C)

Taller C4 - GE II Página 8
4 Intercambiador dde tubos y coraza
Un intercambiador de tubos y coraza (cuatro pasos de los tubos y dos pasos por la coraza) se usa para
calentar 10.000 kg/h de agua presurizada desde 35 °C a 120 °C, usando 5000 kg/h de agua caliente
que entra al intercambiador a 300 °C.
Si el coeficiente global de transferencia de calor es 1500 W/m2 · K, responda usando el método de la
efectividad de NUT:

a) Determine la tasa de transferencia de calor entre los fluidos.

b) Obtenga el valor del ratio de capacidad c.

c) Determine cual es la efectividad del IdC.

d) ¿Cuál es el área necesaria que debe tener el IdC, para que cumpla con el aumento de 35 a 120 °C?

Considere: cp,c 4.195 J/kg · K y cp,h 4660 J/kg · K

Escriba claramente los supuestos utilizados.

Figura 5: Perfiles de temperatura.

Figura 6: Efectividad versus NUT para el IdC.

Taller C4 - GE II Página 9
4.1 Solución:
Supuestos:

• Pérdida de calor insignificante a los alrededores

• Energı́a potencial y cinética insignificante
• Condiciones estacionarias

Pregunta a)
De las ecuaciones del balance de energı́a se puede obtener:

Q̇ = ṁc cp,c · (Tc,o − Tc,i ) (2)

10.000 kg J
Q̇ = h
s · 4195 (120 − 35)◦ C = 9, 905 · 105 W (3)
3600 h kg · K

Pregunta b)
El ratio de capacidades queda expresado por:

Cmin ṁh · Cp,h

= (4)
Cmax ṁc · Cp,c

Cmin 5000 kg J
h · 4660 kg·K
= ≈ 0, 56 (5)
Cmax 10000 kg J
h · 4195 kg·K

Pregunta c)
Para obtener el valor de la efectividad, se utiliza la siguiente fórmula:

q Cmax · (Tc,o − Tc,i )

ε= = (6)
qmax Cmin · (Th,i − Th,o )
1 (120 − 35)◦ C
ε= × = 0, 57 (7)
0, 56 (300 − 35)◦ C

Pregunta d)
Con los valores obtenidos, es posible entrar con ε y Cmin /Cmax a la Figura 6, donde se obtiene un
valor de N T U ≈ 1 1 .
Con esto se puede obtener el valor del área de intercambio 2 :

N T U · Cmin 1 × 5000 kg J
h · 4660 kg·K
AS = = J s
= 4, 314[m2 ] (8)
U 1500 s·m2 ·K 3600 h
Por lo que se necesitarı́a un intercambiador de al menos 4.314[m2 ] para poder elevar la temperatura
del agua presurizada de 35 a 120 [°C]
1
Considerar correcto entre 1 y 1,1
2
Si se usa NTU = 1,1 se obtiene 4,75[m2 ]

Taller C4 - GE II Página 10