1.

Propiedades de una
sustancia pura

Profesor:
Joaquín Zueco Jordán
Área de Máquinas y Motores Térmicos
 Principio de estado
Objetivo de la Termodinámica es relacionar las variables
termodinámicas de un sistema , en equilibrio termodinámico

El nº de propiedades independientes es (1+ nº de modos de

trabajo cuasiestático relevantes)

Sistema simple compresible posee 1 forma de trabajo, p dV

Estado termodinámico definido conocido 2 props. independientes
M.C. P, V, T  Variables termodinámicas
Sistema simple fundamentales

P,V,T f (P,V,T) = 0

La ecuación de estado nos permite hallar una P =P(V,T)

V=V(P,T)
variable fundamental conocidas las otras dos. T=T(P,V)
 Comportamiento de los fluidos, superficie P-v-T
Sustancia que se
Sólido y líquido contrae al solidificar

Líquido Gas

Sólido
o p
j k
Líquido
Vapor
Presión y vapor

Sólido y vapor
Temperatura

Volumen específico
Diagrama PVT para el
agua (caso especial) Sustancia que se
expande al solidificar
Sólido y líquido

Gas
Líquido

Vapor
Presión

Sólido
Temperatura
Proyecciones p-T, p-v y T-v
Diagrama de fases
Diagrama de fases

P T P T

P AGUA
Curva de fusión

Líquido
Punto crítico
Sólido
Curva vaporización

A B
Gas
Punto triple Vapor
C D
Curva de sublimación

T
Diagrama de fases (ejemplos)

i) A presión atmosférica el CO2 sólo

puede existir como sólido o vapor.
Cuando se suministra calor al CO2
sólido, en vasija abierta a 1 bar, se
transforma directamente en vapor (se
sublima), sin pasar por la fase líquida,
de ahí el nombre de hielo seco.

ii) El CO2 líquido sólo puede existir a una

presión mayor de 5,11 atm. Las botellas de
acero en las cuales se almacena
corrientemente contienen líquido y vapor
saturado. La presión en ellas es la presión
de vapor del CO2 a la Tª de la botella.

 Si la temperatura es de 20 ºC, la presión

de vapor es de 56 atm.
Diagrama de fases (ejemplos)

i) La Tª de ebullición aumenta con la

presión. Esta Tª baja a medida que
nos alejamos de la superficie
terrestre ya que la presión disminuye.
ii) Sólido, líquido y vapor pueden
permanecer en equilibrio en el punto
triple, para el cual la presión de vapor es
de 4,5 mm Hg y la Tª es 0,016 ºC.

iii) Si se reduce lo suficiente la Tª de un

vapor cuya presión < P punto triple, el
vapor pasa directamente a sólido. Estas
son las condiciones bajo las cuales se
forma la escarcha y nieve.
El proceso inverso (sublimación) ocurre
también a la nieve (no hay charcos).
 Comportamiento de los fluidos
Introducimos un gas en un cilindro y medimos
P, V en distintos estados de igual temperatura
P
Tª=cte
Tª=cte
Tª=cte
Tª=cte
Tª=cte Tª=cte
Gas
Líquido C
Pc Proceso isotérmico
Un gas no se puede
B2 T3
A2 licuar comprimiéndolo
Vapor
húmedo
Tc a T=cte > Tc
B1 A1 Vapor T2
sobrecalentado
T1
Líquido Vapor saturado
saturado
V
 Comportamiento de los fluidos
Introducimos un gas en un cilindro y medimos
T, V en distintos estados de igual presión
T Vapor sobrecalentado
pc Líquido y vapor
Gas p2 Líquido
Pa
C p1
Tc
Líquido
B2 A2
Vapor Vapor
húmedo sobrecalentado
B1 A1

Líquido Vapor saturado

saturado
Proceso isobárico
V
Fuente: http://www2.udec.cl/~dfiguero/curso/capacidadfases/capacidadfases.htm
T P3 Diagrama T-s
Pc= 221,2 bar
Vapor saturado seco
374,15 ºC C P2 hg sg vg ug
Tc Vapor
sobrecalentado
Líquido saturado
Vapor P1 hf sf vf uf
Líquido húmedo

hfg =TS (sg - sf) Entalpía:

sf sg h=u+pv
s

Calor latente de vaporización

hfg = ug - uf + PS (vg - vf) Ts Tª de saturación
Ps Presión de saturación
Ejemplo: Entropía del propano en estado líquido
a 25ºC y 1 atm hfg = TS sfg
 Temperatura y presión de saturación

Altura P. atm. (kPa) Tsat(ºC) Psat(kPa)

0 101.33 100.0
600
1000 89.55 96.3 AGUA
2000 79.50 93.2 400
5000 54.05 83.0 Tsat= f(Psat)
200
10000 26.50 66.2 0
20000 5.53 34.5 0 50 100 150 200 Tsat(ºC)

• Al psat tiempos de cocción más cortos y ahorro energético

• Al altura disminuye la psat (y Tsat) tiempos de cocción mayores

• Por cada 1000 metros de altura la Tsat disminuye unos 3ºC

 (a presión atmosférica)
Fuente: http://www2.udec.cl/~dfiguero/curso/capacidadfases/capacidadfases.htm
 Diagrama de Mollier, h-s
h

Vapor
Líquido sobrecalentado
comprimido
C

Vapor húmedo

s
Fuente: Termodinámica, 7ª ed. Cengel & Boles
Fuente: Termodinámica, 7ª ed. Cengel & Boles
Fuente: Termodinámica, 7ª ed. Cengel & Boles
 Problema 1.1
200 g de agua en estado de líquido saturado es completamente
vaporizada a 100kPa. Determinar:
- Cambio de volumen
- Cantidad de energía aportada
 Título o calidad de un vapor húmedo
Título de un vapor Grado de humedad
mg mg mf
x =  =  y = 
m mf+ mg m
x kg (1 – x) kg
1 kg
vapor húmedo
= vapor saturado + líquido saturado
seco
P V = Vliq + Vvap
C V Vf Vg
Pc v =   
m m m
Vapor
húmedo
mf mg
T v =  vf +  vg
x m m

vf v vg v
 Título o calidad de un vapor húmedo

mf mg
v =  vf +  vg
m m
v =  x vf + x vg
mg mf
x =   x 
m m

P
C Otras propiedades
Pc termodinámicas
Vapor
húmedo u = ( 1 – x ) uf + x ug
T h = ( 1 – x ) hf + x hg
x
s = ( 1 – x ) sf + x sg
vf v vg v
 Problema 1.2
Un recipiente de 80 litros contiene 4 kg de refrigerante R134a
a 160 kPa. Determinar:
- Tª
- Calidad o título de vapor
- Entalpía
- Volumen ocupado por la fase vapor
Fuente: Termodinámica, 7ª ed. Cengel & Boles
 Aproximaciones a líquido saturado de
un líquido comprimido
Cuando no se disponen datos (a presiones bajas)

Volumen específico: v(P,T)  vf(T)

Energía interna: u(P,T)  uf(T)

Entalpía: h(P,T)  hf(T) + vf(T) [P-Psat(T)]

Ejemplo:
Aprox. líquido saturado a 20ºC (20ºC, 10MPa)
vf (20ºC) = 0,001002 m3/kg 0,0009973 m3/kg
uf (20ºC) = 83,91 kJ/kg 83,31 kJ/kg
h(P,T) = 93,93 kJ/kg 93,28 kJ/kg
 Problema 1.3
Determinar la energía interna de agua en estado de líquido
comprimido a 80ºC y 5 MPa.
 Calor específico
Cantidad de calor necesaria para elevar un grado la
temperatura de la unidad de masa de una sustancia

1 dQ dq Unidades:
c =    c =c(T,P) • J/kg K
m dT dT
• J/kg ºC

cv  c. e. a volumen cte. cp  c.e. a presión cte.

dq dq
cv =  cP = 
dT v dT p
u h
cv =  cP = 
V cte Qv T v P cte Qp T p
cv(v,T) cp(P,T)
 Sustancia incompresible, v≈constante
u u
v ≈ cte  dv ≈ 0 du ≈  dT   dv
 v v T
u = u (T, v) u
du ≈  dT du ≈ cv dT
 v
h=u+Pv sustancia incompresible
h = h (T, P) dh ≈ du + Pdv + v dP
v≈ cte  dv ≈ 0
En un proceso dh ≈ du + v dP
a presión cte dh ≈ du h
P= cte  dP = 0 cv ≈ 
du ≈ cv dT T p
h
cp ≈ cv ≈  cp(T) ≈ cv (T) ≈ c(T)
T p
 Sustancia incompresible, v=constante
2 2
h2 – h1 =  c(T) dT + v (P2 – P1) u2 – u1 =  c(T) dT
1 1

h2 – h1  c (T2-T1) +v (P2 – P1) u2 – u1  c (T2 – T1)

c (kJ/kgK)
Hierro Los calores específicos varían
muy poco con la presión pero
0,4 Cobre sí mucho con la Tª

0,2
Plomo Para gases cp > cv ya que a
P=cte se permite que el
sistema se expanda y ello
100 200 300 400 500 conlleva más gasto de energía
Temperatura (K)
Fuente: http://www.vaxasoftware.com
 Problema 1.4
Determinar el estado termodinámico del agua:
2. Comportamiento
p-v-T en gases

Profesor:
Joaquín Zueco Jordán
Área de Máquinas y Motores Térmicos
 _ Ecuación de estado en gases ideales
P v _
 lim P v J
T P0  T
= R=8,3143 
mol k
_
Pv
 = R
R _T V PV = nRT
v = 
n
P

m m
n =  PV =  R PV = mR T Pv = R T
M M
R .
R =  V V . .
M v =  
.
m m P V = mR T
 Ecuaciones de estado en gases reales
Factor de compresibilidad Z
_ PV = Z n RT
Pv Pv lim
Z =  =  1
RT R T P0  = 1 Pv = Z R T

PV = Z mR T
. .
PV = Z mR T
 Otras ecuaciones de estado en gases reales

Ecuación de Van der Waals

a
P +  v – b = RT
v2
a, b => constantes

- Ecuación de Beattie-Bridgeman

- Ecuación de Benedict-Webb-Rubin
Constantes de van der Waals de algunos gases

http://www.vaxasoftware.com/doc_edu/qui/vanderwaals_esp.pdf
Modelo de gas ideal – Ecuaciones del Virial

Pv PV=nRT
Z =  
R Pv=RT
PV=mRT
. .
PV=mRT

Z = 1 + B(T) P + C(T) P2 + D(T) P3 + ...

Z = 1 + B´(T)/V + C´(T)/V2 + D´(T)/V3 + ...

Cuando P 0 o V  las fuerzas de interacción entre las

moléculas de un gas se debilitan y Z=1 en el límite
 Ley de los estados correspondientes
Conocidas dos variables reducidas Variables reducidas
está determinada la tercera. P v T
f (Pr , vr , Tr ) = 0 Pr =  vr =  Tr = 
Pc vc Tc
1.2 2.5
5.0 3.0
1.1 3.5
Z

1.0 2.0
Factor de compresibilidad

0.9 1.8
0.8 1.6
1.4 Tr
0.7 0.8 1.3
0.6 0.9 1.2
0.5 0.95
1.1
0.4
0.3 Tr=1.0 Diagrama de Nelson-Obert
0.2
0.1
0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 10
Presión reducida Pr
 Ley de los estados correspondientes
Inexactitud en la cercanía del punto crítico

1.2 2.5
5.0 3.0
1.1 3.5
Z

1.0 2.0
Factor de compresibilidad

0.9 1.8
0.8 1.6
1.4
0.7 0.8 1.3
0.6 0.9 1.2
0.5 0.95 Gráfica no muy exacta
1.1
0.4 Zcrítico= 0,27
0.3 En realidad (0,23 - 0,33)
0.2 Tr=1.0
0.1
0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 10
Presión reducida Pr
Fórmula, Peso molecular, R y propiedades críticas

Fuente: Termodinámica, 7ª ed. Cengel & Boles

 Problema 2.1

Determinar el volumen específico de agua a 100 bar y 400ºC

usando:
- Modelo de gas ideal
- Gráficas del factor generalizado de compresibilidad
 Problema 2.2

¿Cuál es el error de considerar CO2 como un gas ideal en los

siguientes casos?
- 30 bar y 10ºC
- 30 bar y 100ºC
 Problema 2.3

CO2 es comprimido siguiendo un proceso politrópico (pv1,2 =cte)

desde el estado inicial de 10 bar y 200ºC. Sabiendo que la presión
final es 30 bar. Determinar la Tª final.
- Modelo de gas ideal
- Ecuación de van der Waals
 Problema 2.4

Calcular el volumen que ocupa un mol de oxígeno a 100 atm y 298 K

considerando que se comporta como un gas de van der Waals.
a= 1,36 L2 atm mol-2; b= 0,0318 L mol-1

La ec. de van der Waals es cúbica con respecto al volumen. Para

resolver este problema utilizamos un método numérico:
Método de Newton Raphson
v = 0,00075625 m3/kg
 Problema 2.5

Una muestra de H2 se encuentra a 25,34 atm y 34,585 K.

a) ¿En qué condiciones de temperatura y presión debe estar una
muestra de cloro para estar en estados correspondientes con el
hidrógeno?
b) ¿Cuál es el factor de compresibilidad de cada uno de estos
gases?
c) ¿Cuál es el volumen molar de cada uno de estos gases?
d) ¿Cuál es el volumen reducido de cada uno de estos gases?
 Experimento de Joule (1843) Modelo de gas ideal
termómetro
Estado inicial (P,V,T)
A BB
gas
idealgas ideal vacio
V
V V Estado final (P/2 , 2V ,T)
agua
Aplicando el
primer principio: du
u =u(T) cv(T) =  u2 - u1 =  cv(T) dT
Q12 = U12 + W12
U =U(T) dT
U12=0 U1= U2

h=u+Pv dh
cP (T)=  h2 - h1 =  cp(T) dT
u(T) + RT = h(T) dT
 ¿Es el vapor de agua un gas ideal?
Error = (vtabla-videal)/vtabla
T 5,0 2,4 0,5 0 0
Error <1%
17,6 0,8 0,1 0
7,4
49,5 2,2 0,1 0
100 bar GAS
25,7 0,2 0
50 bar IDEAL
7,6 0
10 bar
1,6 0
1 bar 0,1
0,1 bar

v
Adecuado: Acondicionamiento de aire, presión es muy baja
No adecuado: Centrales térmicas
Vapor sobrecalentado: T=cte, p h Gas ideal: T=cte, h=cte

s
 Signo de los calores específicos
Proceso adiabático: c = 0 1 δQ δq
c =   
Proceso isotermo: c = ±∞ m dT dT
Calores específicos positivos
P
Proceso 1-4: +/+ = +
2 4
Proceso 1-5: -/- = +
1
Calores específicos negativos
3 isoterma
5 Proceso 1-2: -/+ = -
adiabática Proceso 1-3: +/- = -
v
Un calor específico negativo significa
que la sustancia cede calor aunque la Tª
aumente, o bien, que absorbe calor
aunque la Tª disminuya
cp cv


Fuente: Termodinámica, 7ª ed. Cengel & Boles

Calor específico cp para algunos gases según la Tª

 

Fuente: Termodinámica, 7ª ed. Cengel & Boles

Ecuación polinómica para obtener el calor específico cp
para algunos gases según la Tª

cp = [a + b T + c T2 + d T3 ]

(T en K y cp en kJ/kmol K)
a b c d

Fuente: Termodinámica, 7ª ed. Cengel & Boles

 T
Tabla de gases h(T) = href +  cp(T) dT
ideales Tref

T
Estado de referencia: 0K, h=0 h(T) =  cp(T) dT
0

u(T) = h(T) –R T
T h Pr u vr so
(K)
kJ/kg kJ/kg
0 0 0
200 199,97 142,56
300 300,19 214,07
400 400,98 286,16
500 503,02 359,49
600 607,02 434,78
700 713,27 512,33
800 821,95 592,30
900 932,93 674,58
1000 1046,04 758,94
 Modelo de un gas perfecto

cp cp(T)
2
Es un gas ideal con
calores específicos
1 constantes (valores
medios)
cp,m
cp,m= (cp,1+cp,2)/2

T1 T2 T

u2 – u1= cv,m (T2 – T1)

h2 – h1 = cP,m (T2 – T1)

 Ley de Mayer
diferenciando
du = cv(T) dT h=u+Pv dh du
  R
dT d
dh = cp(T) dT Pv=RT
(gases ideales) cpcv + R

cp 
 =  (  exponente adiabático
cv 
En general varía con la Tª,
disminuyendo al aumentar ésta

Para gases perfectos R

J R = cp - cv cp = 
R = cp - cv   
mol k cp
=  R
J cv
R = cp - cv  cv = 
kg k  
 Valores de los calores específicos y de  para el aire
T (K) cp (kJ/kgK) cv (kJ/kgK) 

cp > cv ya que a p=cte el gas se dilata realizando cierto trabajo para

vencer la presión exterior, y entonces se necesita absorber una
cantidad de energía equivalente a ese trabajo
Valores de los calores específicos para gases de combustión
Problema 2.6
Obtener la variación de energía interna del H2 cuando sufre un
cambio de Tª desde 200K a 800K.
3. El Primer Principio
de la Termodinámica
en sistemas cerrados
Profesor:
Joaquín Zueco Jordán
Área de Máquinas y Motores Térmicos
 Trabajo, energía cinética y energía potencial

Energía cinética Ec = 1/2 m c2 +Q -W

Energía potencial Ep = m g z SISTEMA

+W
-Q
Trabajo W =  F ds

Energía interna Energía almacenada en un sistema

Modos de transferencia de calor

CONDUCCIÓN Energía total
CONVECCIÓN E = U + Ec + Ep
RADIACIÓN
 Trabajo en el eje ds
d
Weje = F ds = M/r (r d) = M d
r
. Weje
Weje=   = 2n M 
dt
. Par: M = F r F
Weje =  Weje dt =  2n M dt

Trabajo eléctrico
I
Wele = V Qc = V I dt

. V
Wele
Wele=  V I
dt
 Enunciados generales del primer principio
Experimental: Procesos adiabáticos en sistemas cerrados

T21,p21 I T21,p21 Weje= 50 kJ

V Batería

Wele= 50 kJ
Proceso A Proceso B

Joule a
principios de
XIX ΔE12 = E2- E1 = -W12,adiab
 Equivalente mecánico del calor
Joule demostró que el calor era una forma de energía y cálculo la equivalencia
entre la caloría (unidad de calor) y el julio (unidad de trabajo).

Si el bloque de masa M
desciende una altura h, la
energía potencial disminuye en
Mgh, y ésta es la energía que
se utiliza para calentar el agua
(despreciando otras pérdidas).

Joule encontró que la disminución de energía potencial es proporcional al

incremento de Tª del agua. La constante de proporcionalidad (el calor
específico de agua) es igual a 4.186 J/(g ºC). Por tanto, 4.186 J de energía
mecánica aumentan la Tª de 1g de agua en 1º C.
Equivalente mecánico del calor
 • La energía no se crea ni se destruye solo se transforma

Energía que + Incremento de Energía que

=
entra - energía almacenada + sale
• Otras expresiones

ΔE12 = Q12 - W12 dE = Q - W

. .
dE/dt = Q12 - W12
Régimen estacionario:
. .
Q12 y W12 (no dependen del tiempo)

. .
dE/dt = 0 Q12 = W12
 • Función de estado
P 1
a W1a2  W1b2 Depende del camino seguido
Q1a2  Q1b2 No son función de estado
b
2 Depende del estado inicial y
ΔU1a2 = ΔU1b2 final no del camino seguido
v

Una función de estado es una variable o función (energía interna,

entalpía, etc) que SOLO depende del estado del sistema (presión, Tª,
etc) y no depende de como el sistema llegó a ese estado.

Las funciones de estado pueden verse como propiedades del sistema,

mientras que las funciones que no son de estado representan
procesos en los que las funciones de estado varían.
 Proceso cuasiestático
A B C D

A-B-C => Proceso cuasiestático (no implica proceso sea reversible)

La compresión pasa por una serie de estados de equilibrio termodinámico
ya que todos los parámetros del sistema varían lentamente. Sucesión de
procesos infinitesimales en los que no se realiza trabajo disipativo.
P
D C-D => Proceso irreversible
Se produce una onda de presión, luego
C la presión no es la misma en todas
B partes del sistema y por lo tanto no hay
A estados de equilibrio termodinámico
v
Trabajo de un sistema (de frontera) Pe
Al aumentar el volumen el sistema realiza
un trabajo contra las fuerzas de la Presión A dx
externa Pe
ztá
P
W =Pe A dx W=Pe dV
dV = A dx

• proceso cuasiestático

 PdV
2
El sistema pasa por una sucesiva
serie de estados de equilibrio W12 =
Pe=P  W=P dV 1
  PdV
2

W12 =
1
P P 2
1

2 1
W>0 W<0
V1 dV V2 V V2 dV V1 V
p 1

B Trabajo depende del

A 2
proceso
Trabajo para
proceso A

V
 Trabajo de un sistema
• proceso irreversible
P Proceso no
Los estados intermedios, no están en cuasiestático
equilibrio termodinámico, debido al 1
desplazamiento rápido del pistón

 PdV
2
V1 V2 V
W12 =
1
 Trabajo de rozamiento (proceso irreversible)
F

+Q -W
Wr
gas
Wr SISTEMA
+W
(solo puede ser recibido) -Q
Calentamiento con ventilador

En realidad, la velocidad finita del pistón va a producir turbulencias

en el sistema, con el correspondiente rozamiento interno (trabajo
disipativo)

 PdV- W
2

W12=WPdV - Wr W12 = r
1
 Trabajo de rozamiento (proceso irreversible)

Q12 + Wr = U2 –U1 + WPdV


2

Q12 + Wr = U2 – U1 + PdV
1

Forma diferencial Q + Wr = dU + P dV

En procesos adiabáticos
dU + P dV = 0
y reversibles
Problema 3.1

Proceso isóbaro en un dispositivo cilindro-pistón hasta

alcanzar el estado de vapor saturado. Determinar:
- Tª del estado final
- Calor cedido al medio ambiente
P Análisis energético de ciclos

• en un proceso cíclico Q=W

Q calor entregado al sistema
v W trabajo neto

F.C.
Motor térmico
Q1
Q1 calor entregado del F.C. al sistema
W
sistema Q2 calor rechazado por el sistema al F.F.
W trabajo neto
Q2 W = Q1 - Q2
W Q1 – Q2 Q2
F.F.  =  =  = 1 - 
Q1 Q1 Q1
P
Ciclos inversos

Máquina frigorífica:
v
Q2 Q2
C.O.P.MF  = 
F.C. W Q1 - Q2

Q1 Bomba de calor
Q1 Q1
sistema
W C.O.P.BC  = 
W Q1 - Q2
Q2
C.O.P. Coeficiente operación
F.F.
 FC T1 FC Q1 T1
Q1

W W

FF FF
Q2 T2<T1 Q2 T2<T1

MOTOR TÉRMICO MÁQUINA FRIGORÍFICA

/BOMBA DE CALOR
 Procesos politrópicos

P vn = K n => (  ,+ ) n = exponente politrópico

P vn = K P1 vn1 = P2 vn2 (cualquier sustancia)

v2


2

 dv 1
w12 = Pdv =  w12=  ( P2v2 – P1v1 )
v n 1–n
1 v1

(GAS IDEAL)
n-1 n-1
Pvn = K P1 vn1 = P2 vn2 
T2 P2 n v1
  
Pv P1v1 P2v2 T1 P1 v2
  R  
T T1 T2
R
w12=  ( T2 – T1 )
1–n
 Procesos isócoros e isobáricos
Isócoros Isóbaros Si n =0
P 2  P v= K
Si n =  P
 P= K
 P v  = K
 P0 v = K 1 2
v=K w12
1

v v


2
w12 = P (v2 –v1)
w12 = Pdv = 0 w12 = 0
1 GAS IDEAL
GAS IDEAL P=K
v1 v2
v=K P1v1 P2v2  
P1 P2   T1 T2
P1v1 P2v2   T1 T2
  T1 T2
T1 T2 w12 = P (v2 –v1) = R ( T2 – T1 )
 Procesos adiabáticos GAS IDEAL
En un proceso cuasiestático: Dividiendo (2) por (1)
δq = du + Pdv cp dT vdP
 
δq = 0 Adiabático cv dT =-Pdv (1) cv dT -Pdv
du = cv dT dP dv
Gas perfecto    0
dh = cp dT P v
cp dT = vdP (2)
dh = du + Pdv + vdP ln P +  ln v = C
δq = 0 Pv = K n=
P
2 -1
 -1

Pv = K T2 P2 v1
Si n =    
T1 P1 v2
1
w12 R
w12=  ( T2 – T1 )
1-
v
 Procesos isotérmicos GAS IDEAL

P
2 Si n = 1 Pv= K

1 P1 v1 = P2 v2

w12 P1v1 P2v2 T 1 = T2

 
T1 T2 isotérmico
v

P1v1 = P2v2

v2


2
dv v2 v2 P1

w12  Pdv  ln  R ln  R ln 
1 v1 v v1 v1 P2
Procesos politrópicos, resumen

P n=
n=1

n=0

n=- n= 
v

Isoentrópicos n = 

Procesos Isotérmicos n = 1 T = cte

politrópicos
Isócoros n =  V = cte

Isóbaros n=0 P = cte

 Comparación: Trabajo de expansión y compresión
1
Procesos politrópicos  w12=  ( P2v2 – P1v1 )
1–n

P P
Trabajo Ahorro en
extra consumo
n=1 n=1

n=  n= 

v1 v2 v v2 v1 v
Expansión Compresión

Interesa que ambos procesos se hagan sin variación de T

Problema 3.2

Una habitación es calentada mediante una resistencia

eléctrica, permaneciendo la Tª del aire inalterada. Sabiendo
que las pérdidas de calor al exterior son 6500kJ//h.
Determinar la potencia del calentador eléctrico.
 Problema 3.3

El aire (a 15ºC) de una habitación adiabática es movido por un

ventilador eléctrico que consume 150 W. Aire puede ser
considerado un gas perfecto. Obtener la Tª transcurridas 10 horas.
Problema 3.4

15 kg de aire es calentado de 25 ºC a 77ºC a presión constante

mediante un calentador eléctrico. Las pérdidas de calor al
exterior son 60 kJ.
Determinar la cantidad de energía eléctrica aportada en kWh.
Problema 3.5

Un dispositivo cilindro-pistón contiene 0,8 kg de nitrógeno

(gas perfecto) y sufre un proceso politrópico (n=1,3) hasta que
su volumen se hace la mitad. Determinar:
- Trabajo, en kJ
- Calor, en kJ
Problema 3.6
Dos tanques contienen vapor de agua en diferentes
condiciones. Se mezclan en proceso no adiabático hasta
alcanzar la presión de 300 kPa. Determinar:
- Tª de la mezcla
- Calor cedido al medio ambiente
4. El Primer Principio
de la Termodinámica
en sistemas abiertos

Profesor:
Joaquín Zueco Jordán
Área de Máquinas y Motores Térmicos
 Aplicación del primer principio a sistemas abiertos

•Conservación de la masa
La masa, como la energía, es una propiedad que se conserva, y
no se crea ni se destruye.
(la masa cambia al cambiar la energía), excepto
E = m c2
en reacciones nucleares, éste cambio es muy pequeño

(Einstein)
ment - msal = mVC

Masa que _ Masa que = Cambio de la

entra V.C. sale V.C. masa en el V.C.

.
Caudal másico: m, (kg/s) . .
ment - msal = dmVC/dt
 Energía nuclear
E = m c2 Energía
generada
(masa
pérdida por
(Einstein) nucleón)

Generación Generación
de energía de energía
(fusión) (fisión)
Número
másico
Núcleos Núcleos Núcleos
ligeros intermedios pesados

Núcleos
intermedios

Núcleos Núcleos
pesados intermedios
Fisión Núcleos
ligeros
Fusión
236 U  144 Ba  90 Kr  2 neutrones
92 56 36
 • Ecuación de continuidad
Flujo unidimensional: Las propiedades en la frontera son
uniformes en la sección transversal
.
1 2 m = A x =  A c t m =c A
c
Ecuaciones de Bernuillí y de
continuidad . c1A1 c2A2
m =  = 
x t
v1 v2
.
c = x/t m = 1A1 c1 = 2A2 c2


.
c c V = cndA
A


.
Perfil unidimensional m = ρ cndA
Perfil real
A
 Interacciones de trabajo en un V.C.
WPdV

Weje Entrada de masa

Salida de gas
masa Wele

. Trabajo de flujo
.
W = P (A c) = P m v

PA c (VC)
. . .
Wflujo = me (P v)e – ms (P v)s
 Aplicación del primer principio a sistemas abiertos

c1 q •Conservación de la energía
1 2
z1 c2 EVC = Q - W + Eent - Esal
w VC
z2

Cambio
Energía de Energía que _ Energía que = energía
calor y trabajo + entra VC sale VC del VC
c12 c22
q – w + u1+ gz1  P1v1 – u2 – gz2 –  – P2v2 = eVC
2 2
Entalpía especifica => h= u + P v
 Aplicación del primer principio a sistemas abiertos

•Conservación de la energía

c12 c22
q – w + h1+ gz1 – h2 – gz2 –  = eVC
 2

V.C. con más de una entrada y una salida

dEVC . . . c2 . c2
 = Q – W + m (h + ---- + g z)  -  m (h + ----+ g z)
dt Entrada 2 Salida 2
 Aplicación del primer principio a sistemas abiertos
Caso de régimen estacionario

- Ninguna propiedad dentro del V.C. cambia con el tiempo

- Ninguna propiedad en la frontera del V.C. cambia con el tiempo

- Los flujos de calor y trabajo no cambian con el tiempo

Como el volumen del V.C. es “cte” el trabajo de frontera WPdV = 0

. .
Conservación de la masa ment = msal

. =m
. . .
m1 2 V1 = c 1 A1  V2 = c 2 A2
 Ecuación estacionaria de la energía
. . . c2 . c2
 = Q – W +  m (h + ---- + g z)  -  m (h + ---- + g z)
Entrada 2 Salida 2
Sistema c22-c21
estacionario con 1 q = h2-h1 +  + g (z2-z1) + w
entrada y 1 salida 2
q  calor específico
w  trabajo específico q , h , g z , c2/2 ,
c2/2  energía cinética por unidad de masaw
gz  energía potencial por unidad de masa
J/kg
2 c2
1 -c2
Q = m [h2-h1 +  + g(z2-z1)] + W
2
Q  calor Q
W  trabajo W Julios
m  masa m  kg
 . . c22-c21 . . .
Q = m[ h2-h1 +  + g(z2-z1)] + W Q=mq
2
. .
Q  flujo de calor Q Watios
. . . .
W potencia W W=mw
. .
m gasto másico m  kg/s

EC y EP suelen ser despreciables:

c de 45 m/s supone EC=1 kJ/kg
z de 102 m supone EP=1 kJ/kg

• Ecuación de Bernuillí En un líquido  = cte  v1 = v2= v

P2-P1 c22-c21 (explicación en un tema posterior)

0 =  +  + g (z2-z1)
 2
Válida para flujo estacionario, no viscoso e incompresible
 Aplicaciones típicas del primer principio a
sistemas abiertos

Turbinas
Compresores
Bombas y ventiladores
Toberas y difusores

Válvulas y tubos aislados

Intercambiadores de calor

Cerrados y abiertos
 Representación Turbina axial
simbólica

1 Ecuación de la
energía aplicada a
turbinas
w
T

c22-c21
q = h2-h1 + + g(z2-z1) + w
2
c1 ≈ c2
q = 0 => Adiabático w = h -h
1 2
z2-z1=> Se desprecia
Problema 4.1

3 kg/s de vapor de agua a 80 bar y 500ºC se expanden en una

turbina adiabática hasta alcanzar los 20 kPa. Sabiendo que la
potencia suministrada es de 2,5 MW. Determinar la Tª de
salida del vapor.
Problema 4.2

Argón (gas perfecto) se expanden en una turbina adiabática

(ver datos en figura). Sabiendo que la potencia suministrada
es de 250 kW. Determinar la Tª de salida del gas.
Representación
simbólica Compresores
c22-c21
2 q = h2-h1 + + g(z2-z1) + w
Ecuación de la 2
w
energía aplicada a

c1 ≈ c2
compresores y
bombas

C q = h -h + w
z2-z1=> Despreciable 2 1

1 q
Proceso adiabático w = h1-h2

de aire rotativo axial

Problema 4.3

0,5 kg/s de CO2 (gas ideal) es comprimido en un compresor

adiabático. Determinar:
- Caudal volumétrico en la entrada
- Trabajo de compresión en kW
 Bombas
cte c22-c21
0 = h2-h1 +  + g(z2-z1) + w
2
2
1 B c1 ≈ c2
w w = h1-h2
z2-z1=> Despreciable

P2-P1
Proceso isoentrópico w  = 


axial axial y radial alternativa

 Procesos de derrame: W =0
Los fluidos experimentan grandes variaciones de velocidad

1 Tobera Difusor 2
2 1
c1 c2 c1 c2

dA<0 dA>0

c22-c21
q = h2-h1 + + g(z2-z1) + w
2
q = 0 =>Proceso adiabático c21 c22
z2-z1=> Se desprecia h1 + = h2 +
2 2
w =0
Para todo tipo de materia (sólido, líquido o gas), el aumento de presión
(dp>0), origina siempre una disminución de volumen (dV<0 o d)

.
m = A c = cte  d+ dA/A dc/c = 0

TOBERA DIFUSOR
dA<0, ddc>0 dA>0, ddc<0

dp<0 dp>0
c2 > c1 c2 < c1
P2 < P1 P2 > P1
Una tobera es un dispositivo Por el contrario, un difusor
diseñado para transformar transforma energía cinética en
entalpía en energía cinética entalpía
La tobera es la parte final del motor de reacción. En ella los gases a muy
alta presión se expanden saliendo muy acelerados y por el principio de
acción-reacción de Newton, causan un empuje hacia delante.
 Válvula y tuberías
1 2
c22-c21
Válvula q = h2-h1 + + g(z2-z1) + w
laminación 2

c1 ≈ c2 Proceso
q = 0 => Adiabático isoentálpico
z2 = z1 h2=h1
w =0

Tubos
c22-c21
q q = h2-h1 + + g(z2-z1) + w
1 2 2
c1 c2 c1 = c2
z2 = z1 q = h2-h1
w=0
 Medición del título de un vapor húmedo
T2
h
P1 válvula
P2
Tubo de P1
muestra P2
1 2 T2

x
calorímetro Salida de vapor
s recalentado
P2 Proceso de laminación
Estado 2 T2 h2 h 2 = h1
P1 h1 – h1f
Estado 1 x = 
h1 Mirando en las tablas h1g – h1f
Problema 4.4
Obtener la presión de salida, estado termodinámico y
energía interna, de R134a que recorre una válvula.
 Intercambiador cerrado
1 P1 = P2
2 P3 = P4
fluido 1
Ecuación de la energía aplicada a
intercambiadores de calor  1 > T2
fluido 2  4 > T3
4 3

c22-c21 c1 = c2
q12 = h2-h1 + + g(z2-z1) + w12 z2 = z1
2 w12 = 0
c24-c23 c3 = c4
q34 = h4-h3 + + g(z4-z3) + w34 z3 = z4
2
w34 = 0
Calor cedido = Calor absorbido
q12 = q34
. .
h2-h1  = h4-h3 mfluido1 = mfluido2
Problema 4.5

8 kg/s de R134a es recorrido por un intercambiador cerrado.

Obtener el caudal másico de agua fría.
 Intercambiador abierto (procesos de mezcla)
1 P1 = P2 =P3
fluido 1
.
Ecuación de la energía aplicada a
intercambiadores de calor fluido 3 Q=0
fluido 2 3
2

T Son despreciables las

p variaciones de EC y EP
. . .
1 0 = m1 h1 + m2 h2  m3 h3
3 .
m  kg/s

V
Aplicación principal:
Calentadores de agua de alimentación Ejemplo
 2 Generador de vapor
vapor de agua
Ecuación de la energía
aplicada a calderas o
generadores de vapor

q Generador
de vapor

agua liquida
1

c22-c21
q = h2-h1 + + g(z2-z1) + w
2

c1 = c2
z2-z1=> Se desprecia q = h2-h1
w=0
 Análisis energético de un compresor IDEAL: Comparación

P a) Sistema cerrado:

C Proceso A-B: WAB = P1 (VB - VA)= P1 V1

D P2
1
Proceso B-C: WBC= ----- ( P2 V2 – P1 V1 )
1 –n
A B P1 Proceso C-D: WCD = P2 (VD - VC) = - P2 V2

V2 V1 V
n
WTOTAL = WAB + WBC + WCD = ------ ( P2 V2 – P1 V1 )
1-n

b) Sistema abierto: Se obtendrá la misma expresión

 Aplicación del primer principio a sistemas abiertos
Caso de régimen no estacionario c 1
1 q
En cada punto del VC
P=P(t) z1
tendremos en cuenta la VC 2
T=T(t)
variación de masa y de c2
c=c(t)
energía w
z z2
•Balance de materia :
Masa que Masa que Variación de dm .- .
entra V.C. sale V.C. = masa en V.C.  = m -  m
dt Entrada Salida
•Balance de energía:
Energía que - Energía que
entra V.C. sale V.C. = Variación de energía en V.C.

dE c 2 c2
. . . .
 = Q –W + m (h + ---- + g z)  -  m (h + ---- + g z)
dt Entrada 2 Salida 2
 Ejemplos de procesos no estacionarios

P,T
P1,T1 Turbina
P1,T1
Q
P2,T2

Carga de un depósito Descarga de un depósito

Arranque de turbinas, compresores y calderas

Problema 4.6
Un depósito rígido y aislado de 0,5 m3 de volumen, que
inicialmente se encuentra lleno de He a P= 1 bar y Tª de 600 K,
está conectado a través de una válvula cerrada a una tubería
que contiene He a P= 50 bar Tª= 300 K. El He se puede
considerar como gas perfecto.
Se abre la válvula para permitir que el He entre lentamente en
el depósito hasta que se equilibran las presiones. Calcular la
temperatura final del contenido del depósito.
 La energía del VC es solamente energía interna (el depósito está
quieto); es adiabático y rígido, y no hay más que un conducto de
entrada, en el que no se considera Ec ni Ep.

Eliminando dt e integrando entre el estado inicial y el final queda:

La energía interna y la entalpía son proporcionales a la Tª; y el gas es

perfecto,
 5. El Segundo
Principio de la
Termodinámica

Profesor:
Joaquín Zueco Jordán
Área de Máquinas y Motores Térmicos
 Dirección de los procesos

Q
Ti>To T To

Tiempo

Los procesos inversos no son

Pi P
Po To posibles espontáneamente

Tiempo

Masa
Masa
zi Masa
zi
 Otros ejemplos de la dirección de los procesos

No todos los procesos que verifican

el 1er principio son viables
Q

1ª Ley 2ª Ley
Masa

La transferencia de calor Q a una

hélice no provocará que gire
Es necesario el 2º principio:
- Saber si un proceso es espontáneo
I=0
- Saber la dirección de un proceso
Q V - Saber si un proceso puede ocurrir

La transferencia de calor Q a un Oportunidad de producir trabajo

alambre no generará electricidad
 Más utilidades del Segundo Principio

Si hay posibilidad de producir trabajo:

- ¿Cuál es la máxima cantidad de trabajo?
- ¿Cuáles son los factores que hacen imposible obtenerlo?

-Establecer las condiciones de equilibrio

-Definición de una escala de temperatura independiente

de la sustancia termométrica

- Desarrollo de expresiones para evaluar u y h en función

de otras propiedades más fácilmente medibles
P • Según el primer principio, en un
proceso cíclico
Introducción al segundo Q =principio
W
Q calor entregado al sistema
v W trabajo neto

• Según el segundo principio

F.C.
Motor térmico Q>W
Q1
Q1 calor entregado del F.C. al sistema
W
sistema Q2 calor rechazado por el sistema al F.F.
W trabajo neto
Q2 W = Q1 - Q2
W Q1 – Q2 Q2
F.F.  =  =  = 1 -  < 1
Q1 Q1 Q1
P Ciclos inversos
Máquina frigorífica:
Q2 Q2
C.O.P.MF  = 
v W Q1 - Q2

F.C. Bomba de calor

Q1 Q1
Q1 C.O.P.BC  =  > 1
W Q1 - Q2
W C.O.P. Coeficiente operación
sistema

Q2

F.F.
 Enunciados del segundo principio
Kelvin Plank Clausius
Es imposible con un motor Es imposible construir una
térmico, producir un trabajo neto, máquina, que funcionando con
en un ciclo completo, un ciclo, no produzca otro efecto,
intercambiando calor que transferir calor desde un
solamente, con un cuerpo a una cuerpo a otro de mayor
temperatura fija. temperatura.

F.C.
T= CTE
T1
Q1 T1 >T2
Q
W
sistema
F.F.
Q2= 0 T2
 Físicos que enunciaron el segundo principio

Kelvin Plank Clausius

 Procesos irreversibles
Aquellos que una vez ha sucedido, es imposible devolver al
sistema y al entorno a sus estados iniciales

Transferencia de calor
Expansión libre de un fluido
Reacción química espontánea
Mezcla espontánea de sustancias con diferente
composición o estado
Rozamiento (fricciones mecánicas/ deslizamiento por la
viscosidad de un fluido)
Internas
Todos los procesos reales son irreversibles
Externas
 Procesos reversibles
Aquellos que una vez ha sucedido, es posible devolver al
sistema y al entorno a sus estados iniciales

Estos procesos no ocurren, sólo se definen ya que algunos

procesos reales son casi reversibles
Saber límite IDEAL

T1
ztá Tobera-difusor subsónicos
Q12 T1 – T2  0
P
Péndulo ideal
T2
Wexp = -Wcomp
 Máquinas de movimiento perpetuo
Son máquinas que no funcionan, al no verificar una de las 2 Leyes

MMP de 1ª especie MMP de 2ª especie

Wnet,sal Qent

CALDERA
Calentador eléctrico CALDERA

Bomba
T T
Wnet,sal
GENERADOR
Bomba
CONDENSADOR Sin condensador

Qsal

(viola la 1ª Ley) (viola la 2ª Ley)

Ciclo de Carnot Isotérmicos
P A
Adiabáticos
B
(Gas ideal)

D C
v
Isotérmicos AB CD
vB
Q1 = QAB = WAB= m R T1 ln 
vA
vD
Q2 = QCD = WCD= - m R T2 ln 
vC
Adiabáticos BC DA
-1 -1
T1 vB = T2 vC vB vC
ln  = ln 
-1 -1 vA vD
T1 vA = T2 vD
 Rendimiento del Ciclo de Carnot
Un motor térmico logrará un rendimiento máximo si funciona
con un ciclo reversible entre dos niveles de temperatura.

vC
T2 ln  T2
Q2 vD
=1 -  = 1 -  c= 1 - 
T1
Q1 vB
T1 ln 
vA

El ciclo de Carnot es totalmente reversible porque:

a) No tiene irreversibilidades internas: los cuatro procesos se supone
que son internamente reversibles
b) No tiene irreversibilidades externas: los procesos isotermos de
intercambio de calor se realizan a la misma temperatura que sus focos

El ciclo de Carnot es un ciclo ideal, irrealizable, pero que se usa

como comparación con otros ciclos.
 Ciclo de Carnot inverso
T P
P2 A
T1 A B B
P1
T2 C
D D C
v
v
T2
•Maquina frigorífica C.O.P.MF 
T1-T2

T1
•Bomba de calor (C.O.P.)BC 
T1-T2
 T1 F.C. Teorema de Carnot
WR WI
R I I  R
Q2 T2
1 -   1 - 
T2 F.F. Q1 T1

Corolarios
T1 F.C.
= f( T1 ,T2)
WR WR
R1 R2
 f - fluido operante
- tipo de máquina
T2 F.F.
R1 = R2
Problema 5.1

Un motor térmico extrae 200 Julios de un foco caliente a una

temperatura de 373 K, y realiza 48 Julios de trabajo, cediendo
152 Julios al entorno, que se encuentra a 273 K. ¿Cuánto
trabajo se pierde debido a las irreversibilidades del proceso?
Problema 5.2

Un kg de agua realiza un ciclo de Carnot, limitado entre las

isóbaras a 1 bar y 15 bar. El inicio de la expansión isotérmica
se realiza con un título de vapor del 25% terminando en el
estado de vapor saturado. Al inicio de la compresión
isotérmica el título de vapor es del 84,9 %.
Determinar:
- Trabajo intercambiado
- Calor intercambiado
- Eficiencia del ciclo
 Ciclo de Carnot, x1 = 0,25, x3 = 0,849
- Evaluar Q, W y eficiencia térmica

T (ºC) 1kg H2O

QAB = UAB + WAB
P2=15 bar
1 2


198,3 B
P1=1 bar
99,63
D ~z wx WAB = PdV
4 3 tzât
A

v
Proceso 1-2: v1 = v1f + x1 (v1g – v1f) = 33,815  10-3 m3/kg
h1 = h1f + x1 (h1g – h1f) = 1331,6 kJ/kg

W12 = m P2 (v2 – v1) = 1  15  (131,8 – 33,815)  102 = 147 kJ

Q12 = m (u2 –u1) + W12 = m (h2 –h1) = 1(2792,2 – 1331,6) = 1460,5 kJ

 Proceso 2-3 Q23 = 0

T (ºC) W23 = m (u2 – u3) = 1  (2594,5 – 2190,7) = 403,8 kJ

u3 = u3f + x3 (u3g – u3f) = 2190,7 kJ/kg
P2
198,3 1 2 Proceso 3-4
P1
99,63 W34 = m P3 (v4 – v3) = - 86,68 kJ
4 3
v3 = v3f + x3 (v3g – v3f) = 1,438 m3/kg
v
Ciclo de Qced TF Q34 373
    Q34 = 1156,6 kJ/kg
Carnot Qapo TC Q12 471

Q34 = m (h4 –h3) h4 = 1177,9 kJ/kg

x4 = (h4 – h4f ) / (h4g – h4f) = 0,337
v4 = v4f + x4 (v4g – v4f ) = 0,5716 m3/kg
 Proceso 4-1 Q41 = 0

T (ºC) W41 = m (u4 – u1) = 1  (1121,3 – 1281,0) = -159,7 kJ

u1 = u1f + x1 (u1g – u1f) = 1281,0 kJ/kg
P2
198,3 1 2 u4 = u4f + x4 (u4g – u4f) = 1121,3 kJ/kg
P1
99,63
4 3
Wneto=W12+W23+W34+W41= 304,4 kJ
v

Wneto
Eficiencia  =  
Qapo

TC
Eficiencia  = 1 -  
TF
Problema 5.3

0,5 kg de aire realiza un ciclo de Carnot (que posee un

rendimiento del 50%), partiendo del estado a 7 bar y 0,12 m3 al
inicio de la expansión isotérmica, donde se aporta la cantidad
de calor de 40 kJ.
Determinar:
- Trabajo intercambiado
- Calor intercambiado
 Ciclo de Carnot, P1= 7 bar, V1= 0,12 m3
= 0,5, Q12= 40 kJ
P 1 - Evaluar Q, W y temperaturas de los focos
TC
2 QAB = UAB + WAB


B
4 C?H ~z
TF 3 WAB = PdV
wx t|Üx
v A
(gas ideal)
T1 = TC = P1 V1/(mR) = 585,4 K

TF/TC = 1 -  = 0,5 TF = 292,7 K

Proceso 1-2: 0
Q12 = U12 + W12

W12 = Q12 = m R TC ln(V2/V1) = 40 kJ V2 = 0,193 m3

 Proceso 2-3: Q23 = 0
P 1 W23 = m (u2 – u3) = 0,5 (423,7 – 208,8) = 107,5 kJ
TC Proceso 3-4: 0
2 Q34 = U34 + W34

Qced TF Q34 292,7

4 TF 3    
Qapo TC Q12 585,4
v
Ciclo de Carnot
Q34 = 20 kJ Q34 = W34 = - 20 kJ
Proceso 4-1: Q41 = 0
W41 = m (u4 – u1) = 0,5 (208,8 - 423,7) = - 107,5 kJ

Se puede comprobar: Wneto 40+107,5-20-107,5

 =   
Qapo 40
6. Entropía

Profesor:
Joaquín Zueco Jordán
Área de Máquinas y Motores Térmicos
 Desigualdad de Clausius
 de un ciclo irreversible   de un ciclo reversible = Reversible
< Irreversible
2
Q2 T2 Q1 Q2 Qi
1 -   1 -         
Q1 T1 T1 T2 i=1 Ti
Tomamos los calores con su respectivo signo

P δQi En un ciclo infinitesimal de Carnot

δQi
   
Ti
i

En todo el ciclo

Ciclo descompuesto en v
infinitos ciclos de Carnot
 dQi
  
Ti
Entropía, una nueva propiedad del sistema 2


C
δQi
   Procesos B
Ti reversibles
A
(sistema cerrado)
1

   
2 2
2
δQ
1
δQ δQ δQ
 +  =0
 = 
T T 1
T 1
T
1 2
A C A B

 
2 En un proceso:
δQ 1
δQ


 +  = 0 2
T T δQ
1 2 S12 = S2 – S1 = 
B C T
1
 Entropía, una nueva propiedad del sistema 2

rev


δQi
 
2 1
  δQ δQ
Ti  +  <0 irr
1
T 2
T
irr rev
1
S21 = S1 – S2 Con proceso irreversible

 
2
2
δQ S2 – S1  δQ
 < S2 – S1 = S12 
T 1
T
1
irr irr/rev

2 2

  
2 SG =0 proceso reversible
δQ δQ δWr
S2 – S1 = 
T + S G =  +
T

T
SG >0 proceso real
1 1 1 SG <0 proceso imposible
 δQ kJ
Entropía S dS   
T K
Entropía específica s
S δq kJ
s   ds   
m T kg K

T Proceso reversible D T Ciclo reversible

2 i
1 T a
g W
Q12 r
a
m Q2
S1 dS S2 S a S
s


2
δQ =T dS Q12= T dS Q1 = Q2 + W
T-S
1
Problema 6.1

Un depósito fijo de 0,5 m3 está lleno de R134a (título de vapor

0,4) y 200 kPa. Es suministrado calor (mediante una fuente
calorífica que se encuentra a 35ºC) al depósito hasta que
aumenta la presión a 400 kPa. Determinar:
- Variación de entropía del sistema
- Variación de entropía de la fuente calorífica
- Variación de entropía total
Problema 6.2

Un dispositivo cilindro-pistón contiene inicialmente vapor de

agua saturado a 200ºC. Sucede un proceso isotermo y
reversible de aporte de calor, hasta alcanzar la presión de 800
kPa. Determinar:
- Calor aportado
-Trabajo desarrollado
 Demostración del c para un ciclo de vapor
T (Vapor) Ciclo de Carnot
P2
Q1
T1 A B
P1
T2 D C
Q2

Q2 - T2 (SD – SC) T2 (SB – SA) T2

=1 - 
Q1
= 1 -    c= 1 - 
T (S – S )
1 B A
T1 (SB – SA) T1
 Ejemplo de ciclo de máximo rendimiento
T
-No es el ciclo de Carnot
A Q1
T1 B
(ciclo regenerativo)
QDA
QBC - Existen muchos ciclos
de máximo rendimiento
T2
D A´´ C B´´
Q2 T2
= 1 - 
T1
D´ A´ C´ B´ S

Q1 y Q2 son los mismos que en Carnot

QBC = QDA (caso ideal o reversible, en la práctica no es posible)
 Diagrama T-s
T

s(p,T)  sf(T)

P=cte

h=cte
v=cte

x=0,2 x=0,9
s
 Diagrama h-s
h

s
 Trabajo reversible e irreversible en un SC
Proceso internamente irreversible:
Q
Q = dU + Wreal dS =  SG
T

dU = T (dS - SG) - Wreal

Wreal = T dS - T SG - dU

Proceso internamente reversible: Wrev= T dS - dU

Wreal,ent  Wrev,ent
Wrev - Wreal = T SG  0
Wreal,sal  Wrev,sal
 Ecuaciones TdS
(combinación del primer y segundo principio)
Proceso internamente reversible:

Q = dU + W
T dS = dU + W
Q
dS = 
T W = p dV

T dS = dU + P dV
dH = dU + P dV + V dP
T dS = dH - V dP (válidas para todo tipo de procesos,
al ser la entropía una propiedad )
Ejemplo: Cambio de fase
T dS = dH - V dP dh = T ds

Ejemplo: Sustancia incompresible

T dS = dU + p dV du = T ds

T 2

Área

s
Tercer principio de la termodinamica
La entropía de una sustancia pura, en equilibrio
termodinámico, tiende a cero, a medida que la
temperatura absoluta tiende a cero.
lim
S=0 Hernst, 1906
T0
Es posible acercarse indefinidamente al cero absoluto, pero nunca
se puede llegar a él. En el cero absoluto el sistema tiene la mínima
energía posible (cinética más potencial).

Esto tiene una importancia considerable en el análisis de

reacciones químicas y en la física cuántica. Los materiales
presentan propiedades extrañas cuando se enfrían a Tªs muy bajas.
Algunos pierden por completo su resistencia eléctrica.
 Cambio de entropía en gases ideales (visto en Primer Curso)

du P • du = cv(T)dT
• T ds = du +Pdv ds =   dv
T T
• dh = cp(T)dT
dh v
• T ds = dh – v dP ds =    dP
T T •Pv=RT

dT dv dT v2
ds = cv(T)  R  s12 =  cv(T)  R ln 
T v T v1

dT dP dT P2
ds = cp(T)   R  s12 =  cp(T)   R ln 
T P T P1

T dT P2
so(T) =  cp(T)  s12 = so(T 2) – so(T 1) - R ln 
0 T P1
Problema 6.3

Aire entra en un compresor. Determinar la variación de

entropía.
 Proceso isoentrópico de gases ideales
P2 P2
0= so(T 2) – so(T 1) - R ln  so(T 2) – so(T 1) = R ln 
P1 P1

P2 exp [so(T2)/R] Pr = exp[so(T)/R]

 
P1 exp [so(T1)/R]
Presión relativa

Volumen relativo

P2 Pr2 V2 Vr2
   
P1 Pr1 V1 Vr1
Cambio de entropía en gases perfectos (visto en Primer Curso)

du P
•T ds = du +Pdv ds =  dv
T T
•du = cv dT
P R
•P v = R T  
T v

dT dv
ds = cv  R 
T v

T2 v2
s12= cv ln  R ln 
T1 v1
Cambio de entropía en gases perfectos

dh v
dh = Tds + vdP ds =  - dP
T T
•dh = cp dT
v R dT dP
•Pv = R T   ds = cp   R 
T P T P

T2 P2
s12= cp ln   R ln 
T1 P1

Análisis de los siguientes procesos reversibles:

Proceso isotérmico Proceso isócoro

Proceso adiabático Proceso isóbaro
 Cambio de entropía en un proceso isócoro

Isócoro v = cte n = 
v5
T
v4
T2 v2
v3 s12= cv ln  R ln 
T1 v1
2 v2
v1
1 T2
s12 = cv ln 
T1
Q12

s
 Cambio de entropía en un proceso isóbaro

Isóbaro P = cte n=0

P5 T2 P2
P4 s12= cp ln   R ln 
T1 P1
P3
2
P2 T2
1 s12 = cp ln 
T1
P1
Q12

s
 Cambio de entropía en un proceso isotérmico

Isotérmico T=cte n=1

P1
R T ln 
q 12 P2 P1
s12 =    R ln 
T T T P2
T2 v2
T4 s12= cv ln  R ln 
T1 v1
1 2 T3
T2 P2
T2 s12= cp ln   R ln 
T1 P1
q12 = w12 T1
(válido también para
s1 s2 gases ideales)
s
 Cambio de entropía en un proceso adiabático

Adiabático Q12 = 0 n =

s12 = s2 – s1= 0
T s 2 = s1
2 P2
T2 v2
0 = cv ln  R ln 
T1 v1

P1
1 T2 P2
0 = cp ln   R ln 
T1 P1
s1 = s2
s
 Cambio de entropía en un proceso adiabático

T2 P2
T cp ln  R ln 
T1 P1
2 P2
T2 P2 R/cp
 
T1 P1

T2 P2 -1/ v1 -1
P1   
1 T1 P1 v2
s1 = s2 
s Pv =K

Se demuestra que un proceso isoentrópico de un gas perfecto

es un proceso politrópico con n=
Balances de
entropía

Profesor:
Joaquín Zueco Jordán
Área de Máquinas y Motores Térmicos
 Aplicación del segundo principio a sistemas cerrados


WPdV
2
Q
S2 – S1 =  SG
1
T
Variación Generación Weje
de entropía Transferencia de entropía
de entropía Debido a irreversibilidades gas
internas, no es una propiedad
Transferencia de entropía que
acompaña a la transferencia de calor

SG > 0 proceso real

SG = 0 proceso reversible

S2 – S1 > 0
S2 – S1 = 0
S2 – S1 < 0
 Aplicación del segundo principio a sistemas abiertos
Acumulación de = Entropía que _ Entropía que + Generación de
entropía en el V.C. entra al V.C. sale del V.C. entropía


. . . dQi dS
SG = m2 s2 – m1 s1 -  +   0
Ti dt

. . . dS
Régimen estacionario m1 = m2 = m  = 0
dt
.
Proceso adiabático
 dQ
i 
Ti
Variación de entropía de una sustancia incompresible

• u = f(T) du P du
ds =  dv ds = 
T T T
• dv = 0

du = c dT T2
s2 - s1 = c ln 
T1

Si el calor específico depende de Tª

c(T)
s2 - s1 =   d
T
Ejemplo: Entropía del agua en estado líquido a 30ºC y 1 atm

s(p,T)  sf(T=30ºC) = 0.4369 kJ/kg K

T2
s2 - s1 = c ln 
T1
T2 (273+30)
s2 = s1 + c ln   4.18 ln  kJ/kg K
T1 (273+25)
 Trabajo en proceso de expansión adiabática de
una sustancia incompresible
Turbina hidráulica
w12 = (h2 - h1) + (c22 – c12)/2 + g (z2 – z1)

T P1
h2 - h1 = (u2-u1) + v (P2- P1)

P2 T2
SG = s12 = c ln   
T1
1
2 T 2 > T1 Real T2  T1+0,5

s
Ideal
w12 = v (P2- P1) + (c22 – c12)/2 + g (z2 – z1)
 Trabajo en procesos estacionarios y reversibles

0 = q – w – (h2 – h1) – ec – ez

w = q – dh – dec – dez


dQi
SG = (s2 – s1) – 
Ti q = Tds – T SG

Despejando el trabajo: w = Tds – T SG– dh – dec – dez

Sabiendo: Tds = dh – v dP

w = dh – v dP – T SG– dh – dec – dez

c12 – c22 2
wr w =    v dP - wr
 1
Trabajo en procesos reversibles
c12 – c22 2
P
w =    v dP - wr 2
 1

wreversible= –  v dP 1

v
Procesos de derrame, W=0
c12 – c22 2
0 =    v dP - wr
 1

Son procesos casi reversibles

Toberas: P y T c12 – c22 2
0 =    v dP
Difusores: P y T  1
Ecuación de Bernoullí w = - v dP – T SG – dec – dez
En un líquido:  = cte  v1 = v2= v
Proceso de derrame: w=0
Proceso reversible, T SG = Wr = 0 (fluidos ideales)

P2-P1 c22-c21
0 =  +  + g (z2-z1)
 2
Tubo horizontal,
c2
 + (P + g z) = cte z2=z1
2
Presión dinámica + Presión estática = cte c2 P2
Problema 6.4

En una tobera adiabática entra vapor de agua a 10m/s, 600ºC

y 40 bar. A la salida la velocidad es de 250 m/s.
Determinar cual de estas presiones puede obtenerse a la
salida:
- 60 bar
- 30 bar
 Proceso en tobera (adiabática), vapor H2O

- Evaluar las tendencias posibles de T2 y P2

1
c2 > c1 2 c21 c22
h1 + = h2 +
2 2

c2= 250 m/s

h2= 3674,4 - 31,2 = 3643,2 kJ/kg
c1=10 m/s
T1= 600ºC h1= 3674,4 kJ/kg
s1= 7,3688 kJ/kgK
P1= 40 bar
ii) Para P2= 60 bar i) Para P2= 30 bar
s2≈ 7,13 kJ/kgK < s1 s2= 7,4620 kJ/kgK > s1
T2≈ 590 ºC T2= 582,7 ºC
(proceso imposible) (proceso posible)
Expansión en tobera (adiabática)

c2 h2 ; T2 y P2
h
h1 1
2
h2 2

Estados posibles
s2-s1 > 0
s
Problema 6.5

En un difusor adiabático entra vapor de agua a 250 m/s,

600ºC y 40 bar. A la salida la velocidad es de 10 m/s.
Determinar cual de estas presiones puede obtenerse a la
salida:
- 60 bar
- 30 bar
- 42 bar
 Proceso en difusor (adiabático), vapor H2O
- Evaluar las tendencias posibles de T2 y P2
2
1 c2 < c1
c21 c22
c2=10 m/s h1 + = h2 +
2 2
h2= 3674.4 + 31.2 = 3705.6 kJ/kg
c1=250 m/s
T1= 600ºC h1= 3674.4 kJ/kg
P1= 40 bar s1= 7.3688 kJ/kgK

ii) Para P2= 60 bar i) Para P2= 30 bar

s2≈ 7.20 kJ/kgK < s1 s2= 7.5345kJ/kgK > s1
T2≈ 620 ºC T2= 610.22 ºC
(proceso imposible) (proceso posible)
 Proceso en difusor (adiabático), vapor H2O

Doble
iii) Para P2= 42 bar interpolación
s2= 7.3848 kJ/kgK > s1
T2= 614.15 ºC
(proceso posible)
Proceso en difusor (adiabático)

c2 h2 ; T2 y P2?
h
h2 2 2
2
h1 1

Estados posibles
s2-s1 > 0
s
Trabajo en procesos politrópicos, estacionarios y reversibles
Procesos politrópicos  P vn = K n => (  ,+ )
P vn = K P1 vn1 = P2 vn2 n = exponente politrópico

v2

 
2
dP -n
w12 = - v dP = -  w12=  ( P2v2 – P1v1 )
P1/n n–1
1 v 1

Gas ideal
n-1 n-1

Pv P1v1 P2v2 T2 P2 n v1
  R     
T T1 T2 T1 P1 v2

-nR
w12=  ( T2 – T1 )
n-1
 Comparación: Trabajo de expansión y compresión
-n
Procesos politrópicos  w12=  ( P2v2 – P1v1 )
n-1

P P
Trabajo Ahorro en
extra P2 consumo
P1
n=1 P1 n=1
P2
n=  n= 

v v
Expansión Compresión

Interesa que ambos procesos se hagan sin variación de Tª,

esto en la práctica es verdaderamente complicado
 Estados permitidos según el 2º Principio
T Estados P2 s2 – s1 = f(Q) + SG
prohibidos
2s
2 P1
i) Q=0 (adiabático) o
Q>0 (aporte de calor)
1
s

T Estados no P2
prohibidos
2s
2 2 P1
ii) Q<0 (pérdida de calor)

1
s Válido para sistemas
cerrados y abiertos
 Rendimiento interno o isentrópico
Expansión adiabática P1
h 1
1 w h1 – h2
w ST =  
T wS h1 – h2s wS P2
w
2 cp (T1 – T2) 2
ST  2s
cp (T1 – T2s)
s
* Gas perfecto h = cp T

DIFUSOR c12 – c22

 h2 – h1
2
EcS (c2s2 – c12 )/2 h2s – h1
SDif =   
Ec 2 2
(c2 – c1 )/2 h2 – h1
 Rendimiento interno o isentrópico
TOBERA

Ec (c22 – c12)/2 h1 – h2

STob =   
EcS (c2s – c1 )/2 h1 – h2s
2 2

(considerando sustancia incompresible)

c12 – c22
  h2 – h1 = c (T2-T1) + v (P2 – P1)
2
T2 (c22 – c12)/2
(isoentrópico) SG = 0 = c ln  STob = 
T1 v (P1 – P2)
 Rendimiento interno o isentrópico
Compresión adiabática P2
h 2
2 2s
|wS| h2s – h1
w SC =  
C | w | h2 – h1 w
wS P1

1 cp (T2s – T1)
SC  1
cp (T2 – T1)

* Gas perfecto h = cp T s

BOMBA (considerando sustancia incompresible)

h2 – h1 = c (T2-T1) + v (P2 – P1)

v (P2 – P1)
T2 SB = 
SG = 0 = c ln  (isoentrópico) h2 – h1
T1
Problema 6.6

Obtener el trabajo y la entropía generada en los siguientes

casos:
- Compresión en una bomba adiabática
- Compresión en un compresor adiabático
Nota: Se trata de agua y el rendimiento isoentrópico es de 0.85
en ambos casos

x=0 x=1
Problema 6.7

Aire (gas ideal) recorre una tobera. Determinar:

- Velocidad de salida
- Eficiencia isoentrópica
- Entropía generada
Problema 6.8

Una turbina adiabática de vapor de agua (rendimiento

isoentrópico del 85%) posee dos extracciones (10% y 90% del
total, respectivamente) a dos presiones diferentes. Determinar
el trabajo suministrado.
 Ejemplo del uso del diagrama T-S Gas ideal
T
WR = (h1 – h2s) = (h1´ – h2s) P1 (adiabático)

WI = (h1 – h2) = (h1´ – h2)

P2

E
x 1 P=cte dh = TdS
p 1´
a
n WR
s
i 2s 2 WI
ó
n WP WP = WR - WI
Trabajo perdido por
3 4 5 irreversibilidades
s
 Ejemplo del uso del diagrama T-S Gas ideal
T
WR = (h2s – h1) = ( h2s– h1´) P2 (adiabático)
2
WI = (h2 – h1) = (h2 –h1´) 2s
P1

C 1´ P=cte dh = TdS
o
m 1
p WP
r WR
e
s WI
i
ó WP = WI -WR
n
Trabajo adicional por
3 4 5 irreversibilidades
s
 Ejemplo del uso del diagrama T-S Cualquier sustancia
T
(isotermo y reversible)
-w12 = -q12 + (h1- h2 ) P2 (no adiabático)

h2 P1
C
o 2 1 P=cte dh = TdS
m
p A q12 = T (s2-s1)
r
e
s
i w12 h1-h2 (h1>h2)
ó
n
3 4 5
s
Problema 6.9

Un intercambiador de calor abierto se produce la mezcla de

nitrógeno en estado gaseoso según se muestra en la figura.
Determinar si el proceso verifica los principios de la
termodinámica.
 - ¿Puede ocurrir según la Termodinámica?

1 bar
. 2 355 K
1 bar 1 Q=0 . . .
Nitrógeno: m2 = m3 = 0.5 m1
.
294 K W=0 3 1 bar
233 K

Son despreciables las variaciones de EC y EP

- 1er Principio
. . . . .
0 = Q  W - m1 h1  m2 h2 m3 h3

h- href = cp (T- Tref); Tref= 0K , href= 0

T1 = 294 = 0.5 T2 + 0.5 T3 = 177.5 + 116.5 (se verifica)

 . 0


- 2o Principio dQi . . . .
0=  - m1 s1 + m2 s2 + m3 s3 + SG
Ti

. .
SG = SG/m1 = s1 – 0.5 s2 – 0.5 s3
. .
SG = SG/m1 = 0.5 (s1 – s2) + 0.5 (s1 – s3)

(s1 – s2) = cp ln(T1/T2) – R ln(P1/P2)

(s1 – s3) = cp ln(T1/T3) – R ln(P1/P3)

SG = 0.022 kJ/kg K > 0 (verifica)

Problema 6.10

Vapor de agua a 7 bar y 240ºC se encuentra inicialmente en

un depósito de 10 m3 de volumen. Mediante un orificio
minúsculo se permite la descarga muy lenta del vapor al
ambiente hasta alcanzar la presión de 1.5 bar.
Determinar:
- Temperatura del vapor en el estado final
- Cantidad de masa evacuada al ambiente
 P1= 7 bar, T1 = 240ºC y P2 = 1.5 bar
- Evaluar la Tª final y la cantidad de masa que sale
T (ºC)
P1 Tabla A-4:
Vapor 1
P1,T1 3 s1 = 7.0641kJ/kg K
V=10m P2
T2
v1 = 0.3292 m3/kg
2
s

Descarga muy lenta proceso reversible T2(P2) = 111.4ºC

x2 = (s2 – s2f) / (s2g – s2f) = 0.9725

v2 = v2f + x2 (v2g – v2f) = 1.127 m3/kg

m = m1 - m2 = V/v1 - V/v2 = 21.51 kg

7. Sistemas de refrigeración
y bomba de calor

Profesor:
Joaquín Zueco Jordán
Área de Máquinas y Motores Térmicos
 INTRODUCCIÓN
Máquina frigorífica y Bomba de Calor
Coeficientes de operación
QF QC
COPMF  COPBC 
W W

Relación entre ambos coeficientes

COPBC  COPMF  1

Máquina frigorífica Bomba de calor

Máquina frigorífica y Bomba de Calor

INVIERNO VERANO
Climatización geotérmica
CICLO INVERTIDO DE Es el ciclo inverso del de potencia de
CARNOT Carnot, que tenía:
• Máxima eficiencia térmica
PROCESOS: • Ciclo totalmente reversible

• 4-1: Evaporación isoterma

• 1-2: Compresión adiabática
• 2-3: Condensación isoterma
• 3-4: Expansión adiabática

Q F
m TF (s a  s b ) TF
Coeficiente de operación: COPMF   

W W  (TC  TF )(s a  s b ) TC  TF
(Máquina frigorífica) COMP  TURB
m
 m

 CICLO INVERSO DE CARNOT
EN LA PRÁCTICA

• Es el ciclo de refrigeración más eficiente que opera entre

dos niveles de temperatura específicos
• Es un modelo inadecuado para los ciclos de refrigeración
+ Los procesos 2-3 y 4-1 pueden ser aproximados en evaporadores y
condensadores reales (procesos isotérmicos de transferencia de
calor)
INCONVENIENTES:
– El proceso 1-2 en la práctica suele evitarse, ya que la presencia de
gotas de líquido puede averiar el compresor
– La expansión en la turbina, proceso 3-4, produce poco trabajo
comparado con el necesario en el proceso de compresión. Además,
la presencia de agua liquida estropea la turbina.
 CICLO IDEAL DE REFRIGERACIÓN POR
COMPRESIÓN DE VAPOR

Diferencias respecto al ciclo invertido de

Carnot:

(1) Existencia de un salto de temperaturas

entre condensador y foco caliente, y
entre evaporador y foco frío

(2) Compresión seca (isoentrópica)

(3) La turbina se sustituye por una válvula

de expansión con el consiguiente
ahorro de inversión inicial y
mantenimiento
 dE vc   Ce2  Cs2 
Balance de energía:   h e  h s  
 Q vc  Wvc  m  gz e  z s 
dt  2 

Coeficientes de operación de
la MF y BC
Q F h h Q C h h
COPMF   1 4 COPBC   2 3
W neto h2  h1 W neto h2  h1

Calor intercambiado en EVAPORADOR

= Área diagrama T-s

Q F  m TEVAP ( s1  s4 )
Problema 7.1
Ciclo refrigeración ideal con R134a, conocidas las presiones y QL. Determinar:
(a) caudal de refrigerante, (b) potencia del compresor,(c) caudal de agua en
condensador si ésta sufre un recalentamiento de 8ºC. Cp(agua) = 4,18 kJ/kgºC
 CICLO REAL POR COMPRESIÓN DE VAPOR

(1) Irreversibilidades durante el proceso WC s h2 s  h1

de compresión, contabilizadas mediante la C  
WC  h2  h1
eficiencia adiabática del compresor
(2) Caídas de presión en evaporador,
condensador y tuberías de conexión
(3) Existencia de subenfriamiento a la
salida del condensador y recalentamiento
a la salida del evaporador
Problema 7.2
Ciclo refrigeración no-isentrópico con R134a, conocidas las presiones y los
calores desprendidos por equipos y desde el exterior. Determinar (a) Tª salida del
compresor, (b) Qpersonas, (c) COP, (d) Caudal volumétrico mínimo en 1.
CLASIFICACIÓN: Por las presiones de trabajo:
- Una mezcla zeotrópica (R4XX) están formadas por dos o más componentes
(refrigerantes puros) de diferente volatilidad. Cuando estas mezclas se evaporan o
se condensan en un sistema de refrigeración, su composición y su temperatura de
saturación cambian (es menor que lo de los componentes individuales).

Al hervir en el evaporador, se evapora un mayor porcentaje del componente más

volátil, y cambia el punto de ebullición del líquido remanente.
- Una mezcla azeotrópica (R5XX) es una mezcla líquida de dos o más componentes
que posee un único punto de ebullición constante y fijo, y que al pasar al estado
gaseoso se comporta como un compuesto puro, o sea como si fuese un solo
componente.
- Si ocurre una fuga de refrigerante liquido o gaseoso, la composición seguirá siendo
la misma, por lo que se efectúa la recarga sin inconvenientes. No así con las mezclas
zeotrópicas.

Quizás la ventaja más importante de un azeótropo es su inusitada capacidad para

mezclar ingredientes inflamables y no inflamables para producir una mezcla no
inflamable estable.

Otra ventaja de los azeótropos es que son muy fáciles de destilar y recuperar.

Nos permiten personalizar (propiedades físicas) la mezcla para adaptarla a cualquier

aplicación o presupuesto. Las posibilidades son casi ilimitadas.
Derivados halogenados de hidrocarburos

- CFC (ChloroFluoroCarbons): R11 (CCl3F) y R12 (CCl2F2)

- HCFC (HydroChloroFluoroCarbons): R22
Aplicación Refrigerante Sustituto
CFC & HCFC HFC
Alta Temperatura
R22 R407C
(AA)
Media T
R12 R134a
(Conservación)
Baja T R502 R404A
 EVOLUCIÓN DE ACUERDOS INTERNACIONALES
Protocolo de Kyoto en 1997, se hizo la clasificación de los
refrigerantes HFC como fluidos de efecto invernadero

El sector de la Refrigeración Comercial es uno de los más implicados

en el calentamiento global del planeta debido a su alto índice de
fugas (afecta directamente) de refrigerante (HFCs) y al significativo
consumo eléctrico (afecta indirectamente).
Reglamento (UE) 517/2014
Limita el uso de refrigerantes en algunas aplicaciones que posean
un GWP elevado >150 (ver en la tabla siguiente que R134a, R407c
y R404a ya no se emplearían).
 REFRIGERANTES ALTERNATIVOS
Debido al alto valor del GWP de los HFC, algunos países han decidido eliminar y
cambiar a NH3, hidrocarburos (HC), Hidro-fluor-olefina (HFO) o CO2, que se
considera que tienen un prometedor futuro con cero ODP y bajo GWP.

Sustituye al R134a en refrigeración de automóviles.

Un nuevo refrigerante HFO, R1234ze. Tiene un GWP ultra-bajo<1 (en

comparación con R134a, que tiene un potencial de 1300, y del hidrocarburo R290,
que tiene un potencial de 3. Su ODP es cero. Además, su duración en la atmósfera
es de sólo 18 días, que es muy inferior a los 13 años para R134a.
PROPIEDADES DE LOS REFRIGERANTES
Criterios de selección
PROPIEDADES DE LOS REFRIGERANTES
Criterios de selección
PROPIEDADES DE LOS REFRIGERANTES
Criterios de selección (otros)
• Bajo punto de congelación.
• No tóxico, no corrosivo, no inflamable
• Estable químicamente
• Bajo coste
• Compatible con lubricantes y aceites comerciales
• Fácilmente detectable por indicaciones de fugas
PROCESO DE EBULLICION: El agua se comporta de
forma similar a un buen refrigerante.

Tª.ebullición Presión
(º C) (bar)
151,1 5,000
Ventajas
110,8 1,500 - Barata
100,0 1,030 - Inerte
89,4 0,700 - No tóxico
- No inflamable
68,7 0,300 - Alto valor de Hfg (a 100ºC es
12,7 0,015 de 2256 kJ/kg)
6,7 0,010 Inconveniente: Es difícil conseguir estas
presiones con compresores normales

- Además, congela a 0ºC (a presión

atmosférica).
 Mezcla zeotrópica

Mezcla zeotrópica

Mezcla azeotrópica

Problema: Aumento de la presión en altas temperaturas

Puede ser el refrigerante del FUTURO, si se vence el reto tecnológico de
las altas presiones de trabajo
 SISTEMAS DE BOMBA DE CALOR

Coeficiente de operación: Q C
m TC (s a  s b ) TC
(Ciclo de Carnot) COPBC   

W 
WTURB (TC  TF )(s a  s b ) TC  TF
COMP 
m m

• Fuentes de energía: aire atmosférico, agua y suelo

• El principal inconveniente de los sistemas que emplean aire es la
formación de escarcha en el evaporador
• Tanto la capacidad como la eficiencia de la bomba de calor
disminuyen de manera significativa a bajas temperaturas
•Las bombas de calor son más competitivas en áreas con veranos
calurosos e inviernos suaves
Problema 7.3
Bomba de calor ideal con R134a, conocidas las presiones y el caudal másico de
refrigerante. Determinar (a) QH, (b) Caudal volumétrico en 1, (c) COP
Ciclos avanzados en sistemas de
refrigeración

Profesor:
Joaquín Zueco Jordán
Área de Máquinas y Motores Térmicos
 SISTEMAS DE REFRIGERACIÓN EN CASCADA
Características
Fuente caliente
• Adecuados para aplicaciones a
QC
temperaturas bajas, que suponen 7 6
niveles de presión inadecuados y un  Condensador

pobre del compresor (tóxico)

WCB
• En cada ciclo puede emplearse Compresor
refrigerantes con las características
más adecuadas
8 5
• Permiten obtener mejoras de hasta un Evaporador
15% en eficiencia energética respecto
Q
a los sistemas clásicos 3 2
QF
T COPMF  Condensador
6 WCA  WCB
WCA
7 2 Disminución Compresor
del Wc
8
3 5 4 1
4 1
Evaporador

QF
Aumento del QF
s Fuente fría
Problema 7.4
Ciclo refrigeración en cascada ideal con R134a, conocidas las presiones y el
caudal másico de refrigerante de mayor presión. Determinar (a) Caudal másico
del otro refrigerante, (b) Trabajo de compresión, (c) COP
SISTEMAS DE REFRIGERACIÓN MULTIETAPA
Características COPMF 
QF
• Por idénticas razones que los Fuente caliente WCA  WCB
sistemas en cascada, adecuados QC
5 4
para aplicaciones a temperaturas Condensador
bajas
• Emplean un único refrigerante, y WCB
Compresor
el intercambiador entre etapas es
sustituido por una cámara de 6
mezcla (mejores características 9
3
de transferencia de calor)
T 4 7
Cámara
de mezcla 2

WCA
5 2 Compresor

6
7 9 3 8 1
Evaporador
8
1 QF
s Fuente fría
 Problema 7.5
Ciclo refrigeración multietapa ideal con R134a, conocidas las presiones.
Determinar (a) Título de vapor x6, (b) qL y trabajo de compresión, (c) COP
 SISTEMAS DE REFRIGERACIÓN DE
PROPÓSITO MÚLTIPLE CON UN COMPRESOR

Características QC
• Sistema económico y práctico 3 2
Condensador
cuando se requiere refrigeración
a más de una Tª
WC
• En el caso de un refrigerador Compresor
doméstico
Q FR  Q FC
COPMF 
WC 4
2 Refrigerador
T 5 6
QFR

3 1
4 5 5ºC 7
Congelador
-15ºC
7 6 1 QFC
s
 SISTEMAS DE REFRIGERACIÓN PARA
LICUEFACCIÓN DE GASES
Características QC
• Muchos procesos de ingeniería a 3
temperaturas criogénicas (<-100ºC)
dependen de gases licuados Intercambiador
de calor WC
• Las temperaturas críticas de helio, 4
Compresor
hidrógeno y nitrógeno son -268,
8
-240 y -147ºC respectivamente 2
1
• El ciclo Linde-Hampson, es uno de
Q Regenerador Gas compuesto
los más empleados,
T 5
3 7 Vapor
recirculado
4 6
5
1
2
8 9 Líquido
9 6 7 extraído

s
 APLICACIONES DE LOS SISTEMAS DE
LICUEFACCIÓN DE GASES

i) N2 líquido: para enfriado del ii) Preparación de propulsores

hormigón, para la conservación líquidos para cohetes
de semen, embriones, sueros,
enzimas, tejidos, células y
productos químicos iii) Producción de O2 (O2
clínico, industria metalúrgica
-sopletes, producción de
acero,…-) y de N2
(refrigerante, producción de
iv) Estudio de propiedades amoniaco), mediante un ciclo
de materiales a Tª bajas: Linde y destilación
La superconductividad
La superconductividad es el fenómeno que se produce en algunas sustancias (ciertos metales y cerámicas) en las que
desaparece su resistividad eléctrica cuando la Tª se reduce por debajo de una Tª crítica, Tc.
 CICLOS DE REFRIGERACIÓN DE GAS

Ciclo Brayton invertido QC

3 2
• El refrigerante permanece todo el Condensador
tiempo en fase gaseosa
• La principal diferencia respecto al Compresor
Wneto
Turbina
ciclo de Carnot invertido es que los
procesos de transferencia de calor
no son isotérmicos. 4
Evaporador
1

• Inconveniente: menor COP y son

QF
más caros que la compresión de
vapor.
• Ventajas:
T 2
– componentes simples más ligeros
(aptos para enfriamiento de
aviones)
– ciclos adecuados para licuefacción 3 TC
de gases al permitir la TF
incorporación de regeneración 1
4
s
 APLICACIONES ADICIONALES DE LA
REFRIGERACIÓN DE GAS
QC
3 2
Condensador
2´
Wneto
Turbina Compresor
Regener
4´
4 1
Evaporador

QF

No están limitados por la Pv del

fluido de trabajo. Por ello, pueden
conseguir una Tª del espacio
refrigerado muy baja. De hecho,
con un regenerador, un gas
refrigerador puede alcanzar
temperaturas criogénicas utilizadas
para licuar muchos gases.
 Problema 7.6
Ciclo regenerativo de refrigeración ideal con He(g), conocidas temperaturas (ver
figura) y la relación de compresión=3. Determinar (a) T4, (b) COP, (c) Trabajo
de compresión
Sistema de aire
acondicionado en un avión.
 CICLOS DE REFRIGERACIÓN POR
ABSORCIÓN (Simple efecto)
Características Ferdinand Carré (1846)
• Sistemas accionados mediante calor
• Competitivos cuando la fuente de
energía térmica es barata (solar,
residual de cogeneración, etc.)
• Ventaja: se comprime un líquido en
lugar de un vapor. El consumo de
energía es muy pequeño.

QF Q
COPMF   F
QG  WB QG
8. Mezcla de Gases y 
Psicrometría

Profesor:
Joaquín Zueco Jordán
Área de Máquinas y Motores Térmicos
Descripción de la composición de la mezcla

mi
Mi =  Masa molecular
ni

m = m1 + m2 + m3 +...+ mj
V

n = n1 + n2 + n3 +...+ nj
T P
ni j
yi =  Fracción molar  yi = 1
n i=1

j
mi  fmi = 1
fmi =  Fracción másica
m i=1
 Descripción de la composición de la mezcla

m1 + m2 +........+ mj j
M =   yi Mi = M
n i=1

V n1 M1 + n2 M2+........+ nj Mj
M = 
n

T P Aire seco
Componente Fracción molar (%)
N2 78,08
O2 20,95
M= 28,97 kg/kmol
Ar 0,93
CO2 0,03
Otros 0,01
 Ejemplos de mezclas de gases
•Aire (mezcla de 21% de oxígeno y 79% de nitrógeno, más otros 
gases en pequeñas proporciones) 
• Cronigón (mezcla de 99% argón y 1% oxígeno), empleado como 
gas de protección en soldadura MAG de aceros inoxidables  
•Trimix (mezcla de oxígeno, helio y nitrógeno) usado en buceo 
técnico y en buceo recreacional avanzado, para evitar los problemas 
derivados de respirar nitrógeno bajo presión 
•Desodorante en aerosol 
•Mezcla de neón, argón y xenón 
•Mezcla de 85% metano, 9% etano, 4% propona y 2% butano. 
•Hexafluoruro de azufre y aire 
•Insecticida en aerosol, para el control de diversas plagas de 
interiores 
•Aire y helio, en globos 
•Nitrox (mezcla de aire, enriquecido en oxígeno) empleado en 
buceo técnico y recreativo
Relaciones P‐V‐T en mezclas de gases reales

Reglas de presión y volumen aditivo

P = P1 + P2 + P3 +...+ Pj V = V1 + V2 + V3 +...+ Vj
T,P
T,V
Zi ni RT Zi ni RT
Pi =  Vi = 
V P

j
Z =  yi Zi
i=1

En general la regla de vol. aditivo es mejor excepto a presiones muy bajas

Relaciones P‐V‐T en mezclas de gases ideales
Son casos particulares de las reglas de presión y volumen aditivos
Modelo de Dalton
ni R T
P = P1 + P2 + P3 +...+ Pj Pi = 
V Ambas son inexactas
Pi ni porque desprecian el
yi =  =  efecto perturbador
P n de un componente en
una mezcla sobre
las variables de
estado de los otros
Modelo de Amagat componentes. Es
ni R T
V = V1 + V2 + V3 +...+ Vj Vi =  decir, se supone que
P no reaccionan unos
gases con otros.
Vi ni
yi =  = 
V n
La ley de Dalton es una aproximación adecuada hasta alrededor de 100
atm, pero falla a partir de entonces. La ley de Amagat, ofrece una
excelente aproximación para presiones incluso superiores a 1000 atm.
Propiedades en mezclas de gases ideales
Ley de Gibbs-Dalton (es el modelo de Dalton reformulado por Gibbs)
Las propiedades de una mezcla de gases ideales se pueden
calcular a partir de las propiedades de los gases constituyentes

j
u =  y i ui
i=1
nu = n1 u1 + n2 u2 +...+ nj uj j
h =  y i hi
nh = n1 h1 + n2 h2 +...+ nj hj i=1
j
ns = n1 s1 + n2 s2 +...+ nj sj s =  y i si
i=1
ncp= n1 cp1 + n2 cp2 +...+ nj cpj
j
cp =  yi cpi
i=1
Procesos con mezclas a composición constante
Estado 1 Estado 2

n1, n2, ...,nj n1, n2, ...,nj

Tini, Pini Tfin, Pfin

j
u =  yi [ ui(Tfin) – ui(Tini)]
i=1

j
h =  yi [ hi(Tfin) – hi(Tini)]
i=1
j
s =  yi [ si(Tfin,pi,fin) - si(Tini,pi,ini)]
i=1
 pi,fin
Tablas de gases ideales si = sio(Tfin) - sio(Tini) – R ln ----
pi,ini

Pi,fin yi pfin pfin pfin

----- = -------- = ---- si = sio(Tfin) - sio(Tini) – R ln ----
Pi,ini yi pini pini pini

Calores específicos constantes

u = cv ( Tfin – Tini) ui = cv,i ( Tfin – Tini)

h = cp ( Tfin – Tini) hi = cp,i ( Tfin – Tini)

Tfin pfin Tfin pfin
s = cp ln ---- - R ln ---- si = cp,i ln ---- - R ln -----
Tini pini Tini pini
 Mezcla de gases ideales
T2 P2
T1 P1 Tf Pf

n = n1+n2

V1 V2 Vf=V1+V2

Estado inicial Estado final

U = Q - W Ui = Uf

Ui = n1 u1(T1) + n2 u2(T2) = Uf = n1 u1(Tf) + n2 u2(Tf)

(proceso iterativo) Tf
 Mezcla de gases ideales
Calores específicos constantes
n1 cv,1 T1 + n2 cv,2 T2
n1 cv,1 (T1 - Tf) + n2 cv,2 (T2 - Tf) = 0 Tf = 
n1 cv,1 + n2 cv,2

n1 R T1 n2 R T2
Vf = V1+ V2 Vf =  
P1 P2

(n1 + n2) Tf
Pf = 
Vf Pf = (n1+ n2) R Tf n1 T1 n2 T2
 
P1 P2
Segundo Principio
0
Qi Si = n1 s1(T1,P1 ) + n2 s2 (T2,P2)
Sf - Si=   SG
T Sf = n1 s1(Tf, y1 Pf ) + n2 s2 (Tf, y2 Pf)

SG = n2 s2(Tf, y2 Pf ) - s2(T2,P2) + n1 s1(Tf, y1 Pf ) – s1(T1,P1)

Tf y2 Pf Tf y1 Pf
SG = n2 cp,2 ln   R ln -------  + n1  cp,1 ln   R ln -------- 
T2 P2 T1 P1

Factores que contribuyen a la generación de entropía

- Diferencia de Tª entre los gases

- Diferencia de presiones entre los gases
- Los gases son distintos
Problema 8.1
A un compresor entran 2 kg/s de una mezcla de gases ideales (de composición
volumétrica 60% de N2 y 40% de CO2). Las condiciones a la entrada son de 1,5 bar y
30ºC, mientras a la salida son de 4,5 bar y 142ºC. Determinar:
a) Potencia consumida
b) Rendimiento isoentrópico

2
4,5 bar
Ecuación de la 142ºC
w
energía aplicada a
compresores y
bombas

q=0
1
1,5 bar
30ºC
 .
Mezcla en volumen: N2 (0,6) y CO2 (0,4), m = 2 kg/s

j
M =  yi Mi = 34,41 kg/kmol
i=1
j
h1-h2 =  yi [ hi(T2) – hi(T1)]/M
i=1
. .
(a) Potencia W = m (h1-h2) = -217,5 kW
j
s2s – s1 =  yi [ si(T2, pi2) - si(T1, pi1)] = 0
i=1
yN2 [ s2so – s1o – R ln P2/P1]N2 + yCO2 [ s2so – s1o – R ln P2/P1]CO2 = 0

(b) Rendimiento isoentrópico

Mediante proceso iterativo: T2s = 400K
h1 – h2s
SC =  
h1(T1) – h2s (T2s) h1 – h2
 Psicrometría
Ciencia que estudia las propiedades termodinámicas del
aire húmedo

Aire Aire seco Vapor de agua

húmedo = R = 0,287 kJ/kgK + R =0,4615 kJ/kgK
Aire húmedo
a v

P = Pa + Pv PaV = ma RaT PvV = mv Rv T

na, ma: aire seco P, T, V

nv, mv: vapor de agua
n, m: mezcla
 Temperatura de rocío TR
Mínima Tª que puede tener el aire húmedo sin que el
vapor de agua condense

T A) Pv, T > TR Aire húmedo

Pg no saturado
C
Pv R) Pv, T = TR Aire húmedo
TA saturado
G A G) P > Pv, T > TR Aire húmedo
TR saturado
R
s
T < TR Aire húmedo
sobresaturado

La cantidad de vapor de agua en el aire húmedo varía desde

cero en el aire seco hasta un máximo, en mezcla saturada
 mv kg v. Pv
Humedad específica  =    = 0,622 
ma kg a. P-Pv

yv Pv Aire saturado 100 

Humedad relativa  = ----- 
yv,sat P,T Pg P,T Aire seco 0 

T
Pg Pv = y v P Pg = yv,sat P
C
Pv
T Grado de humedad
G A

φ = 
s g
¿Cúal es la densidad del vapor de agua en g/m3 en el desierto donde la humedad relativa es del 10% y la Tª ambiente 40ºC?

Pg(T) = 7,3851 kPa

T(ºC)
Pv =  Pg(T) = 0,73851 kPa
Pg
C
Pv v = R T/ Pv = (8,314/18) (40+273)/0,73851 
40 v = 195,76 m3 /kg = 0,19576  m3 /g

Densidad:

Ρ = 1/v = 5,108 g/m3
v
 Aire húmedo en contacto con agua liquida
IDEALIZACIONES:
- Aire seco y el vapor se comportan como gases ideales independientes
- El equilibrio entre ambas fases no se altera por la presencia de aire,
por ello la Pv=Pg(T)

T
Pv = y v P Pg

Aire húmedo
Pv
T
Evaporación
2 1

Agua liquida
s
En realidad la presencia de aire altera la Pv respecto a Pg en una pequeña cantidad
Problema 8.2

Un tanque de 3 m3 contiene aire húmedo a 100ºC, 4,4 bar y Φ1=40%. Si el tanque

es enfriado a 20ºC. Determinar si existe condensación y la cantidad de ésta en
kg/kg aire seco, así como la transferencia de calor.
Entalpía del aire húmedo H = maha + mvhv
H
h =  ha+  hv Origen de referencia (0ºC, 1 atm)
ma
hv  hg(T)
Origen de referencia Tabla A-2
ha = cpaT(ºC)
(0 K, 1 atm)
Ejemplo: Aire húmedo a 40ºC y =0.015 kg v/kg a
ha - tabla aire seco
h = cpa T +  hg =  1.006 ּ 40 + 0.015 ּ 2573.5 =  = 
hv - tabla H2O (g) 78.59 kJ/kg a

Tablas de gases ideales h = ha+  hv =  313.26 + 0.015 ּ (10404.35/18.015)

= 321.9 kJ/kg a

Energía interna del aire húmedo uv  ug(T)

Uso de las mismas tablas o ua = cvaT
Entropía del vapor de agua Sv(T, Pv) = Sg Rv ln
Punto de rocío: Enfriamiento a P=cte
P

T1 y Pv, Pa

T
Pv = y v P Pg
Estado inicial Pv
1
P Tr

3 2
Q
T s

Estado final
 Conservación de la masa y la energía
.
. . Q . .
ma mv1 ma mv2 . .
m a1 = m a2
Aire Aire . . .
m v1 + m w = m v2
húmedo húmedo . .
mw = ma (2 - 1)
1 3 . 2
mw
Agua líquida o vapor

. . . . . .
0 = Q + (ma ha1 + mv1 hv1) + mw hw - (ma ha2 + mv2 hv2)

. .
0 = Q + ma  (ha1 – ha2) + 1 hg1 + (2 - 1) hw - 2 hg2) 

ha1 – ha2 = cpa (T1 - T2) o con tabla A-16

 Técnica de saturación adiabática   (para obtener )
T T1  Tsa ´
Aire no Aire
C Pg(T2)
saturado saturado
T Pv
Tsa 2 1 1 2
3 Tsa, hf
TR
R Agua líquida
s
En la práctica resulta difícil lograr un estado saturado sin utilizar un canal muy largo, resultando
impráctico. Se suele emplear una Tª equivalente a la Tsa, denominada Tª de bulbo húmedo

Hent = Hsal cpa (Tsa - T1) + ´ (hg,sa – hf,sa)

 = 
h1 + ( ´ - ) hf,sa = hg,sa (hg,1– hf,sa)
h1 = cpa T1 +  hg,1 Pg(Tsa)
´ = 
hsa = cpa Tsa+ ´ hg,sa P- Pg(Tsa)
 Psicrómetro  TBS Tª de bulbo seco
TBS TBH  TBH Tª de bulbo húmedo

TBS = TBH aire saturado

TBS - TBH
TBS - TBH aire no saturado
Aire Mirando en tablas 
gasa
humedecida TBS >>> TBH (TBS - TBH)
 Disminuye (mucha evaporación)
TBS > TBH (TBS - TBH)
Psicrómetro normal  Aumenta (poca evaporación)
 Psicrómetro de cabestrillo
 Se obtiene a partir de la medida de TBS y TBH 

Se produce evaporación de agua en torno al bulbo del termómetro húmedo con el consiguiente
enfriamiento. La Tª del agua que humedece el tejido desciende, se produce una ∆T que genera
una TC entre el aire y el agua del termómetro .

Al cabo de cierto tiempo la pérdida de calor del agua por evaporación se iguala a la ganancia
de calor del aire, se alcanza un estado estacionario, alcanzándose el punto de lectura de la TBH
  Humedad relativa
90 70 60 50 40 30 20
Carta psicrométrica
30

 Humedad absoluta kg/kg aire seco

0.025

25 0.020

0.015
20

15 0.010

10

5 0.005
0
-5
-10
0.000
-10 -5 0 5 10 Tª bulbo seco ºC 35 40 45 50 55 60
Procesos de acondicionamiento de aire

Humidificación

Enfriamiento Calentamiento
Deshumidificación
Calentamiento y enfriamiento sensible
. .
Q = ma (h2 - h1)



h2

h1

1 2
 1=  2

TBS T1 T2
Enfriamiento con deshumidificación

1-2 Deshumidificación
. . .
Qe = ma (h1 – h2) - ma ( 1 –2) hf2


2-3 Calentamiento h1 
. .
Qc = ma (h3 – h2) 1
h3 1
h2
3
2 2,3

TBS
 Problema 8.3

Aire húmedo a 25ºC, 1 atm y con una humedad relativa del

50% es enfriado hasta alcanzar la temperatura de 10ºC.
Determinar si existe condensación y la cantidad de ésta en
kg/kg aire seco.
 Aire húmedo a 25ºC, 1 atm, =50% es enfriado, p=cte y 10ºC
- Determinar si ocurre condensación y su valor
T  
Pv = y v P
Pg
Pv 50%

1
Tr a
1 1= a
2 a
3
2 10%

TBS(ºC) 10 14 25
. . .
mw = mv1 - mv2
. .
mw/ma = 1 - 2 = 0,004 – 0,0076 = 0,0024 kg (v)/kg (a)
 Problema 8.4

Un deshumidificador entra aire húmedo a temperatura T1 y

humedad relativa Φ1. A la salida se alcanza el estado de
saturación a la temperatura T2. Sabiendo que el agua
condensada sale a la temperatura de T3.
Determinar:
‐ Cantidad de agua que condensa en kg/kg aire seco
‐ Calor intercambiado
 Deshumidificador, T1, 1, T2, 2(saturación), P1=P2=P3, T3
- Cantidad de agua que condensa y calor intercambiado
.
Q

1 2

3
Agua
. .
mw/ma = 1 - 2

. . . . . .
0 = Q + (ma ha1 + mv1 hv1) - mw hw - (ma ha2 + mv2 hv2)

. .
Q/ ma = (ha +  hg)2 – (ha +  hg)1 + (1 – 2) hf3
Enfriamiento con deshumidificación

.
QE =
 Humidificación Adición de vapor

h2

2
h1 2

1 1

TBS T1 T2

Humidificación con calentamiento adicional, T1 < T2

 Problema 8.5

Un sistema de aire acondicionado trabaja con los datos de la figura.

Determinar:
‐ Calor aportado en la sección de calentamiento
‐ Cantidad de vapor aportado
Humidificación
Inyección de agua líquida
h2 (pulverizada)

h1 2
2

1
1

TBS T2 T1

Humidificación con enfriamiento, T1 > T2

El calor latente de vaporización es suministrado por la

propia corriente de aire húmedo, por lo que ésta se enfría
1 agua 2 Enfriamiento evaporativo
(climas cálidos y algo secos)

Tela mojada h 1=h 2
2´ 2
h1 + (2 – 1) hf = h2 2
h1>> (2 – 1) hf 1
1
h1 h2
TBS T2 T1
Mezcla adiabática de dos corrientes 
.
m1 h1 1
3
.
. m3 h3
m2 h2 2

• Balance materia 
.
ma1 3- 2 

. =  h2
ma2  1 -  3
h3
2
h1 2
3
3
1
1

TBS T1 T3 T2
 • Balance energía
.
ma1 h 3 - h 2 T 3 - T 2

. =   
ma2 h 1 - h 3 T1- T3
Problema 8.6

Una mezcla adiabática de dos corrientes de la figura.

Determinar el estado de salida
Torres de refrigeración “Similar al enfriamiento evaporativo”
COOLING TOWERS La masa de agua evaporada extrae el calor 
latente de vaporización del propio líquido. Este 
calor latente es cedido al aire, por lo que se 
enfría el agua y se calienta el aire.

Balance de materia
. .
ma ( 4 – 3) = m5

h3 = cpa T3+ 3 hg3

h4 = cpa T4+ 4 hg4

Balance de energía
. . . .
h4 + m2/ma hf2 = ( 4 – 3) hf5 + m2/ma hf1 + h3
Torres de refrigeración
 Problema 8.7
Una torre de refrigeración posee los datos de la figura.
Determinar el caudal volumétrico de aire y el caudal másico de
agua de reposición.

Debido a la evaporación del agua, debe aportarse

agua para mantener constante ese caudal, y que
el sistema trabaje estacionariamente.
Balance de materia . .
m5 = ma ( 4 – 3)

Balance de energía
. . . .
h4 + m2/ma hf2 = ( 4 – 3) hf5 + m2/ma hf1 + h3

. .
m5 = ma ( 4 – 3) = 96.9 (0.0273 – 0.0087) = 1.80 kg/s
9. Combustión

Profesor:
Joaquín Zueco Jordán
Área de Máquinas y Motores Térmicos
Combustible + Comburente Productos + energía

masa (PRODUCTOS) = masa (REACTIVOS)

nº de moles (PRODUCTOS)  nº de moles (REACTIVOS)

Energía de activación

Comburente Combustible

Reacción en cadena
• Teórica o neutra Se quema todo el combustible con el aire mínimo

Combustible + ( O2+ 3.76 N2) CO2 + H2O + N2

• Completa Se quema todo el combustible

Combustible + (O2+ 3.76 N2) CO2 + H2O + N2 + O2

• Incompleta No se quema todo el combustible

Combustible + ( O2+ 3.76 N2) CO2 + H2O + N2 +CO

Composición Composición
teórica técnica
21 % O2 M=28.96 kg/ kmol
78 % N2 21 % O2
1 % Ar 79 % N2
Comburente = aire
O2+ 3,76 N2

Ley de Dalton
 C, H, S, H2O,cenizas, etc.
Combustible
Clasificación Composición
Sólido Líquido Gaseoso % peso = % masa
madera gasolina gas natural C. gravimétrica
turba gasoil propano
Cx Hy Hidrocarburo
carbón
Gasolina: C8H18
Gasoil: C12H26
 Características de familias de hidrocarburos
Familia Fórmula Estructura Saturado Terminación y ejemplos
Parafínica CnH2n+2 En cadena Sí “ano” propano, octano
Olefínica CnH2n En cadena No “ileno” o “eno” propileno
Diolefínica CnH2n-2 En cadena No “dieno” butadieno
Nafténica CnH2n En anillo Sí prefijo “ciclo” ciclopentano
Aromática
Benceno CnH2n-6 En anillo No
Naftaleno CnH2n-12 En anillo No

Isómero: Mismo nº de H y C con distinta estructura

H H
Alcoholes: Sustitución de H por OH, CH3OH (metanol) C
H H
H H H H H H C C
H– C–CC–C–H H C H
H– C–C–C–C–H
H H H H H H H H H H
En cadena, saturada En cadena, no saturada
En anillo, saturada
 Separación por destilación fraccionada

Combustible Intervalo
destilación
GLP (butano y 5-20ºC
propano)
gasolinas 20-200ºC

queroseno 180-315ºC

gasóleos 220-350ºC

Fuel-oil 400ºC

Poder comburivoro: Cantidad mínima de aire seco para una

combustión neutra, unidades Nm3/kg comb. o Nm3/Nm3 comb.
Poder fumigeno: Cantidad de productos de combustión para una
combustión neutra, unidades Nm3/kg comb.
• La inflamación del combustible depende:
- Límites de inflamabilidad

- Temperatura de autoinflamación: Tª a la que el vapor se enciende sin una

fuente de ignición
- Mínima energía para su activación

Límites de inflamabilidad: Porcentaje de combustible en la mezcla

para asegurar la combustión:
L i < L < Ls

FACTORES DE LOS QUE DEPENDEN:

T aumenta el campo de inflamabilidad

P aumenta el campo de inflamabilidad

Los gases inertes y el vapor de agua lo disminuyen

 Sustancia Tª inflamación (ºC) Tª de ignición (ºC) Li Ls

Alcohol etílico 392 558 3.28 18.95

Butano 415 482 1.86 8.41
Metano 535 650/750 5.0 15.0
Etano 472 520/630 3.1 12.45
Etileno 490 542/547 2.7 28.6
Benceno 580 740 1.35 6.75
Hidrógeno 4.1 75
Acetileno 2.3 82

Mínima Tª con la que la llama originada

Tª ignición es persistente y duradera, aún retirando la
(flash point)  fuente de ignición.
Tª fuente de ignición > Tª ignición
Máxima Tª a la que puede calentarse un
Tª inflamación  combustible sin riesgo de incendio.
Tª de ignición: 493ºC
 Ecuaciones químicas de la combustión
Reacción estequiométrica o kg aire
ACe= 
teórica kg comb.
Combustible + ( O2+ 3.76 N2) CO2 + H2O + N2
kg aire AC
Reacción real AC=  = 
kg comb. ACe
Con exceso de aire >1
Combustible + ( O2+ 3.76 N2) CO2 + H2O + N2 + O2

Con defecto de aire <1

Combustible + ( O2+ 3.76 N2) CO2 + H2O + N2 +CO

Combustible + ( O2+ 3.76 N2) Productos + SO2

% azufre Termotecnia.- Joaquín Zueco Jordán
 Ecuaciones químicas de la combustión

Reacción estequiométrica o teórica

Combustible + ( O2+ 3.76 N2) CO2 + H2O + N2

kg aire Maire
ACe=  = ACe 
kg comb. Mcomb.

AC AC 1
=    
ACe ACe FR

FR>1--- (<1): Mezcla rica, defecto de aire

FR<1--- (>1): Mezcla pobre, exceso de aire

 Ecuaciones para un hidrocarburo

Reacción real  ( x+y/4)  4.76 28.97

AC= 
(12x +y)
Con exceso de aire >1

CXHY +( x+y/4)( O2+ 3.76 N2)

xCO2 + y/2 H2O + 3.76 ( x+y/4)N2 +( -1)( x+y/4)O2

Con defecto de aire <1

CXHY +( x+y/4)( O2+ 3.76 N2)

aCO2 +bCO +y/2 H2O + 3.76 ( x+y/4)N2
 Análisis de los productos de combustión
Analizador de Orsat
Análisis volumétrico Composición en
base seca
% moles
%H2O %CO2 %CO %N2 %O2

Analizador de gases
electrónico
Equipo que realiza análisis de
gases de combustión. CO2, O2,
CO, Eficiencia , temperatura y
también análisis de SO2 y NOx
Termotecnia.- Joaquín Zueco Jordán
 Problema 9.1 H2 (2 kg/h) se quema con aire (120 m3/h, 30ºC y 1 atm)

H2 - Determinar el coeficiente de exceso de aire

2 kg/h
1
Productos
Cámara ACe = 2.38 kmol a./kmol c.
Aire
30ºC
2 combustión 3
1 atm
AC=120 m3/h H2 + 0.5 ( O2+ 3.76 N2) H2O + 1.88 N2
AC2 AC2 P2
ma =   kg/h
va R/Ma T2
139.8 kg a./h
AC =  kg a./kg c.
2 kg c./h
2.018 kg c./kmol c.
AC = 69.92 kg a./kg c.  kmol a./ kmol c.
28.97 kg a./kmol a.

F 1 AC 4.87
FR =    =   
Fe  ACe 2.38
Problema 9.2 Análisis de los gases del CH4 a 22ºC y 1atm
Combustión neutra

CH4 + 2 ( O2+ 3,76 N2 ) CO2+ 2 H2O + 7,52 N2

Ley de Dalton, Pi = Yi P
nvapor
Psat(22ºC)= 0.026424 bar =  P(1 atm)
(1+nvapor+7.52)
nvapor = 0.23 kmol/kmol c.

Gases de escape: CO2+ 0.23 H2O + 7.52 N2

Agua liquida = 2 – 0.23 = 1.77 kmoles H2O(l)/kmol c.

 Análisis másico del combustible: 0.85 C; 0.13 H2, 0.02 O2
Problema 9.3 - Determinar la cantidad de vapor de agua que condensa
mi Mi ni
Combustible Enfriamiento
Cámara a 25ºC y 1 atm C 85 12 7.083
Aire combustión H2 13 2 6.500

Condensación
O2 2 32 0.063

(7.083 C + 6.5 H2 +0.063 O2) + 10.27 ( O2+ 3.76 N2)

7.083 CO2 + 6.5 H2O + 38.615 N2
nv Pv = Pg(25ºC) = 0.03169 bar
Pv =  P
(nv+ na) na = 7.083 + 38.615 = 45.698 kmol
Cantidad de vapor de agua que condensa:
Despejando: nv = 1.4509 kmol
(6.5 –1.4509)
 kmol /kmol c.

 Influencia de la humedad del aire en la combustión

Combustión completa humedad del aire

nv
CXHY + ( x+y/4)( O2+ 3.76 N2 + 4.76  H2O )
na
nv
x CO2 + y/2 H2O + 4.76( x+y/4)  H2O + 3.76 ( x+y/4)N2
na
Humedad absoluta
nv mv Ma 28.97 moles v mv kg v
 =       =  
na ma Mv 18 mol a ma kg a
Productos de combustión
y/2 H2O  debido al H del combustible
nv
4.76( x+y/4)  H2O  debido a la humedad del aire
na
 Problema 9.4 CH4 se quema con =1.5, HR%=75%
- Determinar la reacción de combustión y Tª rocío prod.

CH4
Productos Mezcla de aire húmedo
Cámara
Aire
25ºC
combustión 0.945 bar Pv =  Pg(25ºC) = 0.0237bar
0.945 bar na = 3  4.76 = 14.28 kmol a.
=75%

Reacción con aire seco

CH4 + 3 ( O2+ 3.76 N2) CO2 + 2 H2O + 11.28 N2 + O2

nv
Pv =  P Despejando: nv = 0.367 kmol
(nv+ na)
 Reacción con aire húmedo
CH4 + 3( O2+ 3.76 N2)+0.367 H2O CO2+ 2.367 H2O+11.28 N2+ O2

Pv en los
productos

14.3 kPa

(hay que interpolar)

Solución mediante THERMOCombustion
 Entalpía de formación
Cambio de energía relacionado con la formación de un _ _
compuesto, a partir de sus elementos constituyentes, q = hf o
en las condiciones de referencia estándar

Tref = 25ºC
Pref = 1 atm
_
H2 a 25ºC q
_ MJ
f , agua =-241.82 
h o
O2 a 25ºC Combustión
H2O a 25ºC kmol

h (T, p)  h f
o
 
 h (T, p)  h (T , p )  h f
ref ref
o
 Δh
 Entalpía de formación
- o
Sustancia 25ºC Formula Estado hf ( kJ/kmol)
Acetileno C2H2 gas +226730
Dióxido de carbono CO2 gas -393520
Vapor de agua H2O gas -241820
Agua (liquida) H2O líquido -285830
Metano CH4 gas -74870
Etano C2H6 gas -84670
Propano C3H8 gas -103840
Butano C4H10 gas -126140
Heptano C7H16 líquido -224390
Octano C8H18 líquido -249950
Oxígeno O2 gas 0
Nitrógeno N2 gas 0

Transferencia de calor desde el reactor:

Reacción exotérmica (entalpía de formación negativa)
Aporte de calor al reactor:
Reacción endotérmica (entalpía de formación positiva)
 Combustión en un flujo estacionario

Reactivos 1 Q
_
1 kmol comb. Cámara 2 q = h2 – h1
combustión Productos _
Aire q = HP – HR

_ _ _
HR = hcomb + nO2 hO2 + nN2 hN2 kJ
_ _ _ 
HP = nCO2 hCO2 + nH2O hH2O + nN2 hN2 kmol
 EVALUACIÓN DE LOS TÉRMINOS DE ENTALPÍA:
Productos de la combustión
Gas ideal (considerado cada componente)  h (T)
h o (tabla A-26 del Cengel) h  h f o  h(T )  h(T )
f  P ref 

Se puede evaluar de dos formas:

- Con las tablas A-17 a A-25
- Integración del C p en la tabla A-2
Aire de la combustión 0
h  h f o  h(T )  h(T )
 A ref 

Combustible
h  h f o  h(T )  h(T )
 C ref 
Entalpía de los gases de combustión en función de la Tª
Tª Oxígeno Nitrógeno Dióxido Vapor de
de carbono agua
K kJ/kmol kJ/kmol kJ/kmol kJ/kmol
298 0 0 D 9364 9904
400 3028 2972 i 13372 14734
s
600 9249 8895 t 22280 20402
800 15838 15045 i 32179 27896
n
1000 22701 21459 t 42769 35882
o
1200 29758 28110 53848 44380
1400 36956 34941 o 65271 53351
r
1600 44269 41913 i 76944 62748
1800 51679 48992 g 88806 72513
e
2000 59189 56156 n 100804 82593
2200 66792 63380 H
112939 92940
2400 74484 70661 125152 83553
Termotecnia.- Joaquín Zueco Jordán
 Cámara de combustión
Aire primario, Aire secundario Aire terciario
25% 35% 45%

Productos de
combustión
Aire
Mezcla
Combustión Dilución

CombustibleTubo de llama 45% Reducir la Tª hasta la

25% Se introduce alrededor 35% Se introduce a través requerida a la entrada de la
del chorro de combustible de orificios practicados en el turbina, por las condiciones
creando una mezcla de aire- tubo de llama para completar impuestas por el punto
combustible relativamente y refrigerar la combustión metalúrgico a la entrada de los
rica con el objeto de obtener (eliminación de NOx). álabes de la turbina.
una Tª elevada, necesaria
para una combustión rápida.
 Entalpía de combustión
Cantidad de energía térmica liberada durante un _ _
proceso de combustión completa (neutra) q = hC
_
R q
Comb.
Cámara P
combustión Productos
Aire
_
h c   ns h s   ne h e hC= HP – HR
P R
_ _ _
HR = hcomb(Tcomb.) + nO2 hO2 (Taire) + nN2 hN2(Taire)
kJ
_ _ _ 
HP = nCO2 hCO2(TP) + nH2O hH2O (TP) + nN2 hN2 (TP) kmol
 Entalpía de combustión
Caso particular (estado estándar): _
R q
1 kmol
Comb. Cámara P
combustión Productos 25ºC
Aire 25ºC
_ _ _
HR = hcomb + nO2 hO2 + nN2 hN2 kJ
_ _ _ 
kmol
HP = nCO2 ( hfo) CO2 + nH2O (hfo) H2O + nN2 hN2

_ _ _ _ _
hC=q = HP – HR = nCO2 ( hfo) CO2 + nH2O (hfo) H2O – hfocomb
 calor liberado _
 Poder calorífico =   hC
kg comb.

P.C.I. Poder calorífico inferior P.C.S. Poder calorífico superior

No se tiene en cuenta el calor Calor liberado cuando los

liberado para vaporizar el agua productos de la combustión son
formada por la combustión del H enfriados hasta la Tª referencia

P.C.I. = P.C.S. – hfg mH2O

hfg = 2441.7 kJ/kg a 25ºC

 - Obtener la entalpía de combustión del etano, C2H6(g), a 25ºC
Problema 9.5 y 1atm.

C2H6
Productos
25ºC Cámara
Aire
25ºC
combustión 25ºC,
1 atm
1 atm

C2H6 + 3.5 ( O2+ 3.76 N2) 2 CO2 + 3 H2O(l) + 13.16 N2

PCS = - hc= 1559.85 kJ/mol C2H6

 - Obtener el PCI a 25ºC y 1atm:

C2H6 + 3.5 ( O2+ 3.76 N2) 2 CO2 + 3 H2O(g) + 13.16 N2

PCI = - hc= 1427.82 kJ/mol C2H6

PM(etano) = 30.07 kg/kmol

PCS = 1559850 kJ/kmol C2H6 /30.07 kg/kmol = 51873.9 kJ/kg etano

PCI = 1427820 kJ/kmol C2H6 /30.07 kg/kmol= 47483.2 kJ/kg etano
 Temperatura adiabática de la llama

R Q=0

Comb. Cámara Productos

0= HP – HR
combustión TP
Aire HP = HR
TP  Tª adiabática de llama

Energía de los productos

Energía de la combustión:
Cesión de calor al entorno
- Combustible
- Tª y P de los reactivos
CONDICIONES
- Combustión completa
- Cantidad estequiométrica de aire
 h h
P R

 n s (h f o  Δ h ) s   n e (h f o  Δ h )
P R e
O
 n s ( Δh )s   n e (Δh)   n e h feo   n s h fso
P R e R P
_ _ _
nCO2 hCO2(TP) + nH2O hH2O (TP) + nN2 hN2 (TP) =
_ _ _
= nCO2 hCO2(298K) + nH2O hH2O (298K) + nN2 hN2 (298K)
_ _ _
+ hcomb(Tcomb.) + nO2 hO2 (Taire) + nN2 hN2(Taire)
- El aporte de más aire diluye los productos, bajando la Tª
- Combustión incompleta tiende a disminuir la Tª
- Existen fugas de calor, no se cumple la condición adiabática
- Tª crece, hay disociaciones de algunos productos, reacciones endotérmicas
Termotecnia.- Joaquín Zueco Jordán
Problema 9.6

Una mezcla en volumen de un 10% de etanol (C2H6O) y 90%

de octano (C8H18) con exceso de aire del 10% produce una
combustión incompleta en una cámara de combustión.
La relación CO2/CO molar es 9. Productos y reactivos se
encuentran a 25ºC y 1 atm. Determinar:
- Temperatura de rocío en ºC
- Calor intercambiado en kJ/kg comb.
- Considerando la condensación del H2O (v)
nv Pv = Pg(25ºC) = 3.1698 kPa
Pv =  P
(nv+ na) na = 6.66 + 0.74 + 1.525 + 47.77 = 56.695 kmoles
Despejando: nv = 1.8559 kmol H2O (v)
nliq = 8.4 - 1.8559 = 6.5441 kmol H2O (liq)
Evolución de Tp con el exceso de aire Observe la
diferencia de
pendiente para
<1 y >1.
Problema 9.7

Hidrógeno (H2) es quemado en una cámara de combustión

adiabática con un exceso de aire del 50%. Determinar la
temperatura adiabática de la llama
Problema 9.8

Una mezcla de O2 (0.12 kg) y CH4 (0.6 kg) contenido en un

depósito es quemada a V=cte hasta alcanzar la Tª de 1200K
de los productos de combustión. Obtener:
- Presión en el estado final en kPa
- Calor intercambiado en kJ/kmol combustible

Sistema
cerrado
7.5 CH4 + 18.75 O2 7.5 CO2 + 15 H2O + 3.75 O2
u=h–pv=h–RT Sistema CERRADO
 Contaminación ambiental

- Afectan al equilibrio terrestre:

Efecto invernadero
Desaparición de la capa de ozono
Alteración de la meteorología, etc.
- Afectan a la salud humana
- Afectan a la vida animal y vegetal
- Ensucian y deterioran los materiales: edificios, vehículos
 Principales causas de contaminación
• CO2 es el principal causante del efecto
invernadero
Carbono • CO es un gas muy tóxico
CO2, CO • Producir combustiones completas, para no
producir CO
•Consumo mínimo de combustible, para
producir menos CO2
• Gas natural, menor producción CO2
• Lluvia ácida, H2SO4
Azufre • Condensación del H2SO4 en chimeneas
SOx • Combustibles con bajo contenido de S
• Destrucción de la capa de ozono
Nitrógeno • Incrementa la lluvia ácida
NOx • El NO2 es un gas venenoso
• A Tªs inferiores a 1300ºC apenas existen
• El exceso
Termotecnia.- de aire favorece su formación
Joaquín Zueco Jordán
COCs – compuestos volátiles orgánicos
Smog – niebla contaminante
 Es un gas que se encuentra
(similar al CO) comúnmente durante la
perforación y producción de
petróleo crudo y gas natural,
además de en instalaciones de
tratamiento de aguas residuales.
Se nota por su olor a huevo
podrido.