SISTEMAS TERMODINÁMICOS
ABIERTOS
Apellidos y nombre: Velázquez Martí, Borja ([email protected])1
Departamento/Centro: 1 Departamento de Ingeniería Rural y Agroalimentaria
Universitat Politècnica de València
Índice general
1. Resumen de las ideas clave 2
2. Introducción 2
3. Objetivos 2
4. Ecuaciones generales de balance de masas y balance de energía en STA 2
5. Sistemas estacionarios 3
6. Estado no estacionario 5
7. Cierre 8
8. Ejercicios propuestos 8
9. Bibliografía 10
Página 1 de 10
1 Resumen de las ideas clave
En este artículo se muestra cómo se aplican el balance de masas y el balance de energía en un
sistema termodinámico abierto (STA) con el objetivo de predecir las condiciones con las que los
fluidos salen del sistema (volumen de control) a partir de las condiciones de entrada, conociendo
la energía intercambiada; o bien, determinar las condiciones a las que deben entrar los fluidos en
el sistema para lograr unas condiciones de salida prefijadas. Se analizan primeramente sistemas
en estado estacionario, es decir, sistemas en los que los flujos de masas y energía que entran y
salen del sistema son equivalentes, por tanto la energía contenida en el volumen de control es
constante. Posteriormente, se analizan sistemas no estacionarios, sistemas en los que la energía
interna del volumen de control se modifica con el tiempo.
2 Introducción
El control climático de un recinto, tal como un invernadero, una sala de curación de embutidos,
una granja o un edificio de oficinas, implica aplicar balances de masas y energía en un sistema
termodinámico abierto. Por otra parte, los mismos principios aplicados en el control climático se
pueden aplicar a procesos industriales, tales como control térmico de reactores, salas de secado
o dispositivos de acondicionamiento de fruta.
3 Objetivos
Una vez que el alumno se lea con detenimiento este documento, será capaz de:
1. Aplicar las ecuaciones generales del balance de masas y balance de energía en sistemas
termodinámicos abiertos.
2. Analizar sistemas termodinámicos estacionarios con mezclas de fluidos
3. Analizar sistemas termodinámicos estacionarios sin mezclas de fluidos
4. Analizar sistemas termodinámicos no estacionarios
4 Ecuaciones generales de balance de masas y balance de
energía en STA
a) Balance de masas: La variación de la masa contenida en el volumen de control es igual a la
diferencia entre el flujo que entra y sale del mismo.
dmvc
= ṁe − ṁs (1)
dt
b) Balance de energía: La variación de energía contenida en el volumen de control por unidad
de tiempo ( dEdtvc ) es igual al calor que entra menos el trabajo que sale, más la energía asociada
al flujo másico de entrada, menos la energía asociada al flujo de salida.
La energía asociada a los flujos másicos es el producto de cada flujo por su entalpía, más la
energía cinética y potencial del mismo.
dEvc 1 1
= Q̇ − Ẇ + ṁe · (he + · c2e + g · ze ) − ṁs · (hs + · c2s + g · zs ) (2)
dt 2 2
Página 2 de 10
Simbología
mvc es la masa contenida en el volumen de control.
ṁe y ṁs es el flujo másico de entrada y salida respectivamente (kg/s).
t es el tiempo (s).
Evc es la energía total contenida en el volumen de control (J).
Q̇ es el flujo de calor que intercambia el volumen de control con el exterior (W).
Ẇ es el trabajo introducido o extraído del volumen de control por unidad de tiempo (W).
he y js es la entalpía específica del fluido de entrada y salida del volumen de control (J/kg).
ce y cs es la velocidad del flujo de entrada y salida del volumen de contrl respectivamente (m/s).
g es la gravedad (m/s2 ).
ze y zs es la cota de entrada y salida de los flujos en el volumen de control (m).
5 Sistemas estacionarios
Se dice que un sistema está en estado estacionario cuando la masa y la energía interna contenidas
en el volumen de control permanecen constantes. Esto significa que se cumplen dos condiciones:
dmvc
=0
dt
dEvc
=0
dt
- Todo el flujo de masa que entra en el sistema es igual al que sale
ṁe = ṁs
- Toda la energía que entra en el sistema es igual a la que sale
1
0 = Q̇ − Ẇ + ṁ · (he − he + · (c2e − c2s ) + g · (ze − zs ))
2
En numerosos sistemas las diferencias de energía cinética y potencial son despreciaples en
comparación de las transferencias de energía en forma de calor y trabajo. Es por ello que
pueden eliminarse de la ecuación resultando:
0 = Q̇ − Ẇ + ṁ · (he − hs )
Procedamos a mostrar como se aplican estos balances con un ejemplo.
Página 3 de 10
Ejemplo 1
Se dispone de un depósito mezclador que funciona a 10 bares de presión y posee tres
entradas y una salida. Por la primera entrada circula un flujo de agua de 2 l/min a 100ºC;
por el segundo circula un flujo de 3 l/min a 45ºC; y por la tercera entrada circula un flujo
de agua de 5 l/min a 220ºC. ¿Cual es el caudal y la temperatura del flujo de salida? ¿Cual
es el estado de ese flujo: líquido, vapor o mezcla de ambos? En caso de ser mezcla, calcula
el título de la mezcla.
Resolución
Hay que advertir que la mezcla de caudales sólo es posible si todos ellos están a la misma
presión. En este caso están todos a 10 bares. La temperatura de saturación del agua a
esta presión es 179,9 ºC. Esto significa que los flujos 1 y 2 se encuentran en estado líquido
subenfriado. El flujo 3 se encuentra en estado vapor sobrecalentado.
Inicialmente buscamos las entalpías de los fluidos que entran en el mezclador (volumen de
control). En el caso de líquidos las entalpías se pueden calcular mediante la expresión:
h = Cp · T
Donde para el agua Cp = 4, 18 kJ/kg ºC, tal que
h1 = 4, 18 · 45 = 188, 1 kJ/kg
h2 = 4, 18 · 100 = 418 kJ/kg
Para la determinación de la entalpía del flujo 3, que está en estado vapor, nos ayudamos
de las tablas de las propiedades termodinámicas del agua. Interpolando se comprueba que
h3 = 2874, 15 kJ/kg
Página 4 de 10
El flujo se salida se calcula como
ṁ4 = ṁ1 + ṁ2 + ṁ3 = 3 + 2 + 5 = 10 kg/min
Y la entalpía del flujo de salida (ṁ4 ) se calcula como
ṁ1 · h1 + ṁ2 · h2 + ṁ3 · h3
h4 = = 1599, 2 kJ/kg
ṁ4
Se observa que la entalpía del flujo de salida se encuentra entre las entalpías de líquido
saturado y la de vapor saturado a 10 bares.
hliqsat = 762, 81 kJ/kg
hvapsat = 2778, 11 kJ/kg
Eso significa que la temperatura del flujo de salida es exactamente 179,9 ºC, siendo una
mezcla de líquido y vapor.
h4 = x · hvapsat + (1 − x) · hliqsat
h4 − hliqsat 1599, 2 − 762, 81
x= = = 0, 415
hvapsat − hliqsat 2778, 11 − 762, 81
Esto significa que el 41,5 % del flujo saliente es vapor y un 59,5 % es líquido.
6 Estado no estacionario
En un estado no estacionario bien la masa contenida en el sistema varía con el tiempo, o la
energía acumulada en el sistema varía con el tiempo. También se puede dar el caso de que
ambos varíen con el tiempo.
Dado que la variación energía almacenada del sistema a presión constante se puede expresar
como la ecuación (3), donde m es la masa del sistema y Cp es el calor específico (J/kg ºC),
el análisis de estos sistemas generalmente requiere resolver una ecuación diferencial de primer
orden, bien de variables separables o una ecuación diferencial lineal.
dE dT
= m · Cp · (3)
dt dt
Vamos a ilustrar este caso con un ejemplo.
Recordatorio resolución ecuaciones diferenciales lineales de primer orden
La ecuación general de una ecuación diferencial lineal de primer orden es:
y 0 + a(x) · y = b(x)
Página 5 de 10
Este tipo de ecuación siempre tiene solución a través de la fórmula:
R
Z R R
−
y(x) = e a(x)·dx
· b(x) · e a(x)·dx
+ C · e− a(x)·dx
Si la ecuación diferencial lineal es de coeficientes constantes se tiene
y0 + a · y = b
Por tanto la resolución es sencilla
R
Z R R
−
y(x) = e a·dx
· b·e a·dx
+ C · e− a·dx
b a·x
y(x) = e−a·x · · e + C · e−a·x
a
b
y(x) = + C · e−a·x
a
La C es una constante de integración que se calcula a partir de las condiciones de contorno.
Ejemplo 2
Se posee un calentador de agua de 500 l en el que existe una resistencia de 2000 W. La
temperatura inicial del depósito es 80 ºC. El flujo de entrada y salida de agua del mismo
permanece constante, siendo de 1 litro/min. El agua que entra en el depósito es de 15 ºC.
¿Cual se la temperatura del agua de salida al cabo de 3 min? ¿Cual es la temperatura a la
que se estabiliza en sistema?
Resolución:
Se aplica la ecuación del balance de energía en STA, donde Ẇ = 0 y se consideran
despreciables las variaciones de energía cinética y potencial.
dEvc
= Q̇ − Ẇ + ṁ · (he − hs )
dt
La entalpía de los líquidos puede calcularse como h = Cp · T . En caso del agua Cp = 4, 18
kJ/kg ºC. Por otra parte, téngase en cuenta que la temperatura de salida (T ) es igual a
la temperatura de mezcla en el depósito. Entonces se tiene que:
T
m · Cp · = Q̇ + ṁ · Cp(Te − T )
dt
Ordenando términos se tiene
dT ṁ Q̇ ṁ
+ ·T = + · Te
dt m m · Cp m
que es una ecuación diferencial lineal de coeficientes constantes: y 0 + a · y = b , donde
Página 6 de 10
ṁ Q̇ ṁ
y = T, a = ,yb= + · Te
m m · Cp m
Este tipo de ecuación siempre tiene solución a través de la fórmula:
b
y(x) = + C · e−a·x
a
Sustituyendo:
ṁ
Q̇ − ·t
T = + Te + C · e m
ṁ · Cp
Para calcular C nos fijamos en las condiciones de contorno. En el tiempo inicial, t = 0 ⇒
T = T1
Q̇
T1 = + Te + C
ṁ · Cp
Q̇
C = (T1 − Te ) −
ṁ · Cp
" # ṁ
Q̇ Q̇ − ·t
T = + Te + (T1 − Te ) − ·e m
ṁ · Cp ṁ · Cp
Esta ecuación también se puede escribir del siguiente modo:
ṁ ṁ ṁ
Q̇ − ·t − ·t − ·t
T = · 1−e m
+ Te · 1 − e m
+ T1 · e m
ṁ · Cp
Fíjense que la temperatura en el volumen de control tiende a un valor fijo con el tiem-
po suficiente, ‘por lo que pasaría de un estado no estacionario (transitorio) a un estado
estacionario.
Q̇
t→∞⇒T = + Te
ṁ · Cp
Sustituyendo las variables por su valor se tiene que
t = 180 s
m = 1000 l = 1000 kg
ṁ = 1kg/min= 0, 0166 kg/s
Cp = 4, 18 kJ/kg ºC
Página 7 de 10
0, 0166
2 2 − ·180
T = + 15 + (80 − 15) − · e 1000 = 50ºC
0, 0166 · 4, 18 0, 0166 · 4, 18
A los 3 min la temperatura de salida del agua del depósito será 50ºC
2
T = + 15 = 43, 71ºC
0, 0166 · 4, 18
Con el tiempo suficiente la temperatura del agua de salida se estabilizará en 43,71 ºC
Si se representa gráficamente la ecuación del transitorio se puede estimar fácilmente la
duración del transitorio.
7 Cierre
El análisis de los sistemas termodinámicos abiertos requiere la aplicación de las ecuaciones de
balance de materia y energía.
Figura 1: Análisis del sistema
En caso que el sistema esté en estado estacionario, la determinación de las variables suele
obtenerse mediante la resolución de un sistema de ecuaciones lineales. En caso que el sistema
esté en estado no estacionario (transitorio), la determinación de las variables suele obtenerse
mediante la resolución de un sistema de ecuaciones diferenciales.
Tras el periodo transitorio la temperatura del sistema se estabiliza alcanzando un estado esta-
cionario.
8 Ejercicios propuestos
Ejercicio 1
Un mezclador de agua que funciona en estado estacionario, tiene 2 entradas y una salida. En la
entrada 1 el vapor entra a presión p1 =7 bar, con temperatura T1 = 200ºC, y un flujo másico
de 0,10 kg/s. En la entrada 2 entra agua líquida a p2 =7 bar y T2 =40ºC, un flujo másico de
0,04 kg/s Determínese la salida del flujo, temperatura y el título a la salida del mezclador.
Sol: ṁs = 0, 14 kg/s, Ts = 165ºC, x = 0, 33
Página 8 de 10
Ejercicio 2
Un mezclador de agua que funciona en estado estacionario, tiene 2 entradas y una salida y una
resistencia de 150 kW. En la entrada 1 el vapor entra a presión p1 =7 bar, con temperatura
T1 =200ºC, y un flujo másico de 0,10 kg/s. En la entrada 2 entra agua líquida a p2 =7 bar y
T2 =40ºC, un flujo másico de 0,04 kg/s Determínese temperatura a la salida.
Sol: ṁs = 0, 14 kg/s, Ts = 344, 21ºC, x = 1
Ejercicio 3
En un intercambiador de un condensador entra vapor a 1 bar con título 0,45 y el condensado
líquido sale a 1 bar y 45ºC. El agua de refrigeración entra como líquido saturado al condensador
a 20ºC y sale a 35ºC sin cambio de presión. Determínese la relación de flujos másicos entre el
agua de refrigeración y el vapor a condensar y la energía transferida entre fluidos por unidad
de masa. (Los flujos no se mezclan, sólo intercambian calor).
Figura 2: Esquema de un intercambiador de calor
Sol: ṁ2 = −45, 72 · ṁ1 kg/s, Q̇ = 2276, 96 kJ/kg
Ejercicio 4
Un tanque tiene inicialmente 45 kg de agua a 45ºC. Tiene una entrada y una salida de flujos
másicos iguales de 270 kg/h. El agua líquida de entrada está a 45ºC. Un serpentín refrigerante
inmerso en el agua captura 7,6 kW de energía. Por otra parte el tanque tiene un agitador que
proporciona 0,6 kW. Determínese como varía la temperatura en función del tiempo.
Sol: T = 23, 63 + 21, 38 · e−0,0016·t
Ejercicio 5
Dado el calentador de la figura, siendo el rendimiento de la instalación del 80 %, determina:
a) El flujo y la temperatura de salida y potencia de la caldera
b) El consumo específico de combustible si se utilizan pélets con poder calorífico de 25 MJ/kg
c) Volumen del silo de almacenamiento de pélets para una autonomía de 7 días si poseen una
densidad aparente de 500 kg/m3
Sol: ms = 20 l/s, Ts = 50 ºC,Q̇c = 4762 W
Página 9 de 10
Figura 3: Sistema de calefacción de agua
9 Bibliografía
Albright, L. D. (1990). Environment control for animals and plants. Ed. American Society of
Agricultural Engineers.
ATECYR (2010). Fundamentos de climatización: para instaladores e ingenieros recién titulados.
Ed. Asociación Técnica Española de Climatización y Refrigeración (Atecyr)
ATECYR (2016). Fundamentos de refrigeración. Ed. Asociación Técnica Española de Climati-
zación y Refrigeración (Atecyr)(Fernández Seara, J. et al.)
Bakker, J. C., Bot, G. P. A., Challa, H., van de Braak, N. J. (Eds.). (1995). Greenhouse climate
control: an integrated approach. Wageningen Academic Publishers.
Carbó Ballester, J. (2021). Reglamento de Instalaciones Térmicas en los Edificios. Ed. Mar-
combo.
DTIE 3.01, psicrometría (Pinazo Ojer, José Manuel | García Lastra, Arcadio | ATECYR)
DTIE 7.05, cálculo de cargas térmicas (Pinazo Ojer, José Manuel | Soto Francés, Víctor Manuel
| García Lastra, Arcadio | ATECYR)
Fenton, Donald L. (2010). Fundamentals of refrigeration I-P. Ed. American Society of Heating,
Refrigerating and Air-Conditioning Engineers (ASHRAE).
Pinazo, JM (1996) Cálculo de instalaciones frigoríficas. Editorial UPV
Real Decreto 1027/2007, de 20 de julio, por el que se aprueba el Reglamento de Instalaciones
Térmicas en los Edificios.
Real Decreto 552/2019, de 27 de septiembre, por el que se aprueban el Reglamento de seguridad
para instalaciones frigoríficas y sus instrucciones técnicas complementarias.
Stoecker, Wilbert F. (1998). Industrial refrigeration handbook. Ed. McGraw-Hill Education
Velázquez Martí, B. (2017). Aprovechamiento de la biomasa para uso energético. Ed. UPV-
Reverté. 885 pp.
Página 10 de 10
