Termo Aplicada

Descargar como docx, pdf o txt
Descargar como docx, pdf o txt
Está en la página 1de 78

Fundamentos de Termodinámica

FUNDAMENTOS DE TERMODINÁMICA

CARLOS ARTURO BELLO BLANCO

Ingeniero Químico, Magíster en Ingeniería Química

UNIVERSIDAD DEL ATLANTICO


FACULTAD DE INGENIERIA
PROGRAMA DE INGENIERÍA QUÍMICA
Fundamentos de Termodinámica

INTRODUCCIÓN

Al ingeniero químico le corresponde en su perfil ocupacional, el análisis, el estudio, la


investigación y la aplicación de los principios físicos, químicos y biológicos que brinda la
naturaleza, con el fin de diseñar, simular, controlar y optimizar procesos fisicoquímicos y
bioquímicos con el fin de obtener productos valiosos a partir de materias primas en
armonía con la conservación de los recursos naturales y el medio ambiente.

Los siguientes principios y leyes fundamentales son herramientas para el ingeniero


químico:

- Principio de la conservación de la materia (balance de materia)


- Principio de la conservación de la energía (balance de energía) o primera ley de
la Termodinámica
- Principio del incremento de la entropía (balance de entropía) o segunda ley de la
Termodinámica
- Análisis exergético (balance de exergía)
- Leyes del equilibrio de fases y equilibrio químico (balance de potenciales químicos)
- Principio de la Conservación de la cantidad de movimiento (segunda ley de
Newton)
- Leyes de la cinética química
- Correlaciones empíricas y modelos matemáticos
- Principios de diseño, control y simulación de procesos.

Las primeras herramientas arriba mencionadas están relacionadas con la termodinámica.


La Termodinámica es la ciencia que estudia la energía, sus transformaciones y las
propiedades involucradas en dichas transformaciones; sin embargo, es frecuente
analizar los procesos desde el punto de vista energético y se descuida el enfoque basado
en la segunda ley. La energía tiene cantidad y calidad y este último concepto es más
importante que el primero. La termodinámica es una herramienta fundamental en la
interpretación, análisis, diseño, control y simulación de procesos químicos (reactores
químicos) y procesos de separación (operaciones unitarias) tales como absorción,
destilación, secado, evaporación, extracción líquida, lixiviación, etc.

El objetivo de la Termodinámica es la energía, pero surge la pregunta: ¿qué es la


energía? Inmediatamente se relaciona con la capacidad para realizar un trabajo. ¿Y si se
usa para calentar o enfriar un fluido? Una definición más general sería que la energía es
la capacidad para causar un cambio de las propiedades de un sistema o cambio de
estado o proceso. Esta capacidad se mide en cantidad, pero esta cantidad de energía
puede tener menor o mayor capacidad de producir trabajo dependiendo del valor de sus
propiedades con respecto a las del entorno. Surgen otras preguntas: ¿qué es propiedad
de un sistema? ¿Qué es sistema? ¿Qué es estado? ¿Qué es proceso? ¿Qué diferencia
hay entre un proceso reversible y un proceso irreversible? ¿Qué dirección tiene un
Fundamentos de Termodinámica

proceso? ¿Qué restricciones debe cumplir un proceso? además de la primera ley, se


debe satisfacer la segunda ley.

Considérese los siguientes procesos:


1) En una turbina adiabática entra un fluido a una presión y temperatura altas donde
se expande hasta una presión y temperatura bajas. ¿Cómo se determina la
máxima potencia que puede generar y la temperatura mínima de salida?

2) En un intercambiador de calor, por ejemplo, vapor de agua saturado a 150° C


intercambia energía calorífica con un fluido frio que entra a 80° C y se calienta
hasta 140° C. El vapor se condensa y sale como líquido saturado a 150° C o como
líquido a una temperatura menor a 150° C pero mayor que 80° C. El proceso
mencionado es posible, pero el proceso contrario, es imposible. ¿Cómo se explica
esto? Estos conceptos se tratarán más adelante.

3) En una válvula entra líquido saturado de un refrigerante a alta presión donde se


expande adiabáticamente hasta una presión baja. ¿Por qué no es posible que el
refrigerante a la salida no alcanza el estado de la entrada si el contenido de energía
de la salida es igual al de la entrada?

La primera ley no da criterios para responder esta pregunta; sin embargo, la segunda ley,
sí. Esto indica que no basta realizar un análisis energético porque no arroja elementos
para determinar qué tan eficiente es el proceso.
Fundamentos de Termodinámica

1. CONCEPTOS BÁSICOS Y PRINCIPIOS FUNDAMENTALES

Un sistema es una porción del universo sometido a análisis termodinámico que tiene sus
fronteras que lo separa de sus alrededores o resto del universo. Las fronteras pueden
ser reales o imaginarias, rígidas o móviles y a través de ellas se puede transferir energía
en forma de trabajo o calor. Un sistema cerrado es aquel que no intercambia materia
con sus alrededores, aunque puede intercambiar calor y/o trabajo con sus alrededores.
Este tipo de sistema se denomina masa de control. Un sistema abierto es aquel que
intercambia materia con sus alrededores además de intercambiar energía de la materia,
calor y trabajo. Este sistema se denomina volumen de control. Un sistema adiabático
es aquel en donde no hay transferencia de calor con los alrededores y puede ser cerrado
a abierto.

Vapor de agua
a 3500 kPa y
450° C

Producción de
Calor trabajo
Q
Vapor a 10
kPa y 46° C

Sistema cerrado: cilindro-émbolo Sistema abierto: turbina de vapor

El sistema cerrado mostrado arriba es un dispositivo cilindro-émbolo que contiene un


fluido con una determinada masa constante al que se transfiere una cantidad de calor Q
a través de la pared del cilindro y como consecuencia hay un aumento de la energía del
fluido y una producción de trabajo de expansión debido al aumento del volumen, dado
2
por 𝑊 = ∫1 𝑃𝑑𝑉; mientras que el otro sistema es una turbina que tiene una corriente de
entrada de un fluido con un alto nivel de energía donde se produce trabajo de eje como
consecuencia de la expansión del fluido el cual sale a un nivel bajo de energía. La turbina
puede ser adiabática.

El estado termodinámico es la condición misma de un sistema y está especificado por


el valor de sus propiedades que lo diferencian de otro estado.

Las propiedades de un sistema son las características particulares de un estado


termodinámico, medidas directa o indirectamente por instrumentos a través de señales
analógicas o digitales.

La regla de las fases de Gibbs establece los grados de libertad de un sistema de NC


componentes y NF fases homogéneas. Un componente es una sustancia pura o especie
con una composición química invariable que hace parte de una mezcla o solución. Una
fase homogénea es el estado de una sustancia pura o una mezcla de componentes que
puede ser sólida, líquida o gaseosa. Los grados de libertad representan las variables o
Fundamentos de Termodinámica

propiedades intensivas independientes de un sistema para estar completamente


especificado.

Según la regla de las fases de Gibbs los grados de libertad están dados por:

π = NC − NF + 2

Una sustancia pura, NC = 1, en una fase homogénea (sólida, líquida o gaseosa), NF = 1,


tiene dos grados de libertad, es decir que hay que especificar dos variables o propiedades
intensivas independientes. Este resultado está en concordancia con el postulado de
estado.

El agua es una sustancia pura tiene un estado termodinámico definido por su presión y
su temperatura, las cuales determinan su fase. A las condiciones de entrada de la turbina
de vapor mostrada anteriormente (3500 kPa y 450° C), el agua se encuentra en fase
gaseosa o vapor y su estado y el valor de las otras propiedades estarán establecidos por
estas dos propiedades especificadas. La variación de la presión a 10 kPa y la temperatura
a 46° C, establecerá otro estado y otras propiedades.

El estado termodinámico de un sistema se expresa matemáticamente como:

z = f(P, T, y)

La representación gráfica de estas variables en el espacio son los denominados


diagramas de fases, donde z es una propiedad termodinámica.

El estado termodinámico de una mezcla binaria (dos componentes) en una fase


homogénea se establece completamente especificando tres propiedades intensivas
independientes, por ejemplo, la presión, la temperatura y la fracción molar o másica de
uno de los componentes.

Las propiedades extensivas son aquellas que dependen de la cantidad de materia del
sistema, mientras que las propiedades intensivas son independientes de la materia.
Las propiedades pueden ser intrínsecas o extrínsecas. Las propiedades intrínsecas son
aquellas que están ligadas a la naturaleza química y física de la materia.

Las propiedades termodinámicas intensivas son matemáticamente exactas, es decir son


independientes de la trayectoria y son funciones de estado o funciones punto, es decir
que su valor está dado por la condición del estado.

Así, el cambio de una propiedad termodinámica está dado por:


2

∆z = z2 − z1 = ∫ dz
1
Fundamentos de Termodinámica

Existen unas cantidades termodinámicas que no son propiedades del sistema y son
funciones de trayectoria porque no son matemáticamente exactas. El calor y el trabajo
son formas de energía en transición y surgen durante un proceso; su cambio está dado
por:
2

Q12 = ∫ δQ
1
2

W12 = ∫ δW
1

Las propiedades termodinámicas extensivas e intensivas más importantes se denotan


con letras mayúsculas y minúsculas, respectivamente y son las siguientes:

Propiedades extensivas Propiedades intensivas


Volumen, V m3, pies3 Volumen específico, v m3/kg, pies3/lbm
------------------ ------------------ Presión, P kPa, psia
------------------ ------------------ Temperatura, T K (° C), R (° F)
Energía interna, U kJ, BTU Energía interna específica, u kJ/kg, BTU/lbm
Entalpía, H kJ, BTU Entalpía específica, h kJ/kg, BTU/lb
Entropía, S kJ/K, BTU/R Entropía específica, s kJ/(kg-K), BTU/(lbm-R)
------------------ ------------------ Capacidad térmica, Cp, Cv kJ/(kg-K), BTU/(lbm-R)
Energía de Gibbs, G kJ, BTU Energía de Gibbs específica, g kJ/kg, BTU/lbm

Un proceso termodinámico es la serie de etapas que sufre un sistema desde un estado


inicial hasta un estado final. El camino que recorre el sistema durante el proceso es la
trayectoria. Algunos procesos se caracterizan porque el valor de una propiedad
permanece constante.

a) Proceso isovolumétrico o isocórico: volumen constante; dV = 0


b) Proceso isotérmico: temperatura constante; dT = 0
c) Proceso isobárico: presión constante; dP = 0
d) Proceso isentálpico: entalpía constante; dH = 0
e) Proceso isentrópico: entropía constante; dS = 0
f) Proceso adiabático: la transferencia de calor es igual a cero; dQ = 0

Energía del sistema. Un proceso termodinámico se caracteriza porque la energía del


sistema se transforma de una forma a otra y surgen formas dinámicas como el trabajo
(W) y el calor (Q) como resultado de las interacciones con el resto del universo o
alrededores. Esto significa que la energía de un sistema es la capacidad que tiene para
causar un proceso. La energía total de un sistema puede aumentar o disminuir como
consecuencia de las interacciones energéticas con los alrededores en forma de calor y/o
trabajo.

La energía de un sistema (E) está constituida de energía térmica, energía mecánica,


energía cinética, energía potencial, energía magnética, energía eléctrica, energía química
Fundamentos de Termodinámica

y energía nuclear. Estas formas de energía se clasifican en energías extrínsecas y


energía intrínsecas.

Las formas de energía extrínsecas son aquellas independientes de la naturaleza química


de la materia y es la suma de la energía cinética y energía potencial, principalmente y de
energía magnética, en menor grado.

Las formas de energía intrínsecas son las energías térmica, mecánica, química y nuclear,
y están relacionadas con la estructura molecular de la materia y están representadas por
la energía interna (U).

Despreciando los efectos magnéticos, eléctricos y de tensión superficial, la energía de un


sistema está representado por:

u̅2
E = U + Ec + Ep = mu + m + mgZ
2

La energía interna de un sistema es el resultado de las interacciones moleculares,


vibracionales, rotacionales y cinéticas en función de la temperatura y la presión. Una
forma de la energía interna es la energía sensible que está asociada con la variación de
la temperatura y presión sin cambio de fase, mientras que la energía latente está ligada
al cambio de fases.

Un fluido sometido a calentamiento sufre un aumento de su energía sensible debido al


aumento de la temperatura. Una sustancia pura en fase líquida que se convierte en vapor
sufre un cambio de fase debido a la energía latente de vaporización suministrada.

1) Energía latente de vaporización (de líquido a vapor)


2) Energía latente de condensación (de vapor a líquido)
3) Energía latente de fusión (de sólido a líquido)
4) Energía latente de solidificación (de líquido a sólido)
5) Energía latente de sublimación (de sólido a vapor)
6) Energía latente de sublimación inversa (de vapor a sólido)

Trabajo. Es la energía que surge al aplicar una fuerza F con un desplazamiento dX o


cuando se eleva un peso venciendo el efecto de la fuerza gravitacional:
2

W = ∫ FdX
1

Hay diferentes tipos de trabajo:


- Trabajo mecánico que puede ser de expansión o compresión.
Fundamentos de Termodinámica

V2

Wexp = ∫ PdV ; V2 > V1


V1

- Trabajo de eje o de torsión (torque τ) de N revoluciones


2πN

W𝜏 = ∫ τdθ
0

- Trabajo eléctrico debido a un diferencial de voltaje y una intensidad de corriente:


t t

Welec = ∫ VIdt = ∫ I2 Rdt


0 0
Welec en J, V en voltios, I en amperios, R en ohm, t en s

Calor. El calor es una forma de energía que surge debido a una diferencia finita de
temperatura entre dos cuerpos o entre el sistema y sus alrededores. Los mecanismos de
transferencia de calor son: conducción, convección y radiación. Estos mecanismos se
tratan con más detalle en los cursos de transferencia de calor.

1.1. Principio de la Conservación de la materia.

El principio de la conservación de la materia establece que la materia del universo no se


crea ni se destruye sino que se conserva aunque puede transformarse cuando hay
reacciones químicas. Sin embargo, la variación de la materia de un sistema es igual a la
variación de la materia del resto del universo (alrededores). Es decir:

δm δm
( ) +( ) =0
dt sistema dt alrededores

Para fluidos en procesos de flujo este principio se expresa mediante la ecuación de


continuidad. La ecuación de continuidad de un elemento de volumen de control está dada
por:
Fundamentos de Termodinámica

Velocidad de acumulación velocidad de flujo velocidad de flujo


{ de materia en el } + { de materia que sale } − { de materia que entra } = 0
volumen de control del volumen de control al volumen de control

Considerando un volumen de control al cual entran NCE corrientes y salen NCS


corrientes de flujo másico, la ecuación de continuidad queda expresada por:

NCS NCE
δm
[ ] + ∑ ṁj − ∑ ṁi = 0 (1 − 1)
∂t vc
j=1 i=1

Donde ṁj , ṁi = flujos másicos de las corrientes de salida (j) y entrada (i),
respectivamente.

En general,
Ak u̅k
ṁk =
vk

Ak es el área de flujo, u̅k es la velocidad media y vk es el volumen específico de la


corriente k.

Si el proceso ocurre a estado estacionario, el término de acumulación de materia en el


volumen de control es cero y la ecuación de continuidad se convierte en:

NCS NCE

∑ ṁj = ∑ ṁi
j=1 i=1

Si el proceso es de carga, esto es, no hay corrientes de salida, la ecuación (1-1) se


transforma en:
NCE
δm
[ ] − ∑ ṁi = 0
∂t vc
i=1

Si el proceso es de descarga, esto es, no hay corrientes de entrada, la ecuación (1-1) se


transforma en:
NCS
δm
[ ] + ∑ ṁj = 0
∂t vc
j=1

La ecuación (1-1) para un sistema cerrado se reduce a:

δm
[ ] =0
∂t mc
Fundamentos de Termodinámica

1.2. Principio de la conservación de la cantidad de movimiento.

El principio de la conservación de la cantidad de movimiento representa también la


segunda ley de Newton y establece que:

velocidad de acumulación de suma de fuerzas


{ }={ }
cantidad de movimiento en el volumen de control externas

Con la aplicación de este principio a un elemento de volumen diferencial se obtiene la


siguiente ecuación de variación de la cantidad de movimiento en forma vectorial:

Dv
ρ = −∇p + ∇. [τ] + ρ. g
Dt

El término de la izquierda es la derivada sustancial del vector de la velocidad v con


respecto al tiempo y representa la acumulación de cantidad de movimiento debido a la
aceleración y las fuerzas de inercia debido al movimiento. El primer término de la derecha
representa la fuerza de presión. El segundo término es la matriz del esfuerzo cortante y
representa las fuerzas viscosas. El tercer término representa la fuerza gravitacional
debido al peso del elemento diferencial.

1.3. Primera ley de la termodinámica.

El principio de la conservación de la energía establece que la energía del universo se


conserva. La forma general de este principio se expresa de la siguiente manera:
Velocidad de velocidad de flujo velocidad de flujo velocidad neta de
acumulación de energía que sale de energía que entra flujo de energía que
de energía en el + con la masa − con la masa − entra desde los =0
volumen de desde el volumen hacia el volumen alrededores a través
{ control } { de control } { de control } { de las fronteras }

Los términos de transporte de energía son la energía interna, la energía cinética, la


energía potencial que acompañan a la masa que fluye hacia o desde sistema. La
velocidad neta de energía que entra desde los alrededores está representada por el flujo
de calor, la potencia de flujo y la potencia debido al trabajo de eje a través de la frontera
física del sistema.

Para el mismo volumen de control, la ecuación del principio de la conservación de la


energía está dada por:

NCS NCE
dU
[ ] + ∑ [ṁj (h + ec + ep ) ] − ∑ [ṁi (h + ec + ep ) ] − (Q̇ neto − Ẇneto ) = 0 (1 − 2)
dt vc j i
j=1 i=1

Donde


Ecuación de variación de cantidad de movimiento en Fenómenos de Transporte
Fundamentos de Termodinámica

u = energía interna, en kJ/kg


h = entalpía, en kJ/kg
̅2
𝑢
ec = 2 = energía cinética, en kJ/kg
ep = gz = energía potencial, en kJ/kg
Q̇neto = δQ⁄dt = flujo de calor neto suministrado, kW
Ẇneto = δW⁄dt = potencia de eje neta producida, en kW

Para un proceso estacionario, el término de acumulación de energía se reduce a cero:

NCS NCE

∑ [ṁj (h + ec + ep ) ] − ∑ [ṁi (h + ec + ep ) ] − (Q̇neto − Ẇneto ) = 0 (1 − 3)


j i
j=1 i=1

Si el proceso es de carga, esto es, no hay corrientes de salida, la ecuación (1-2) se


transforma en:
NCE
dU
[ ] − ∑ [ṁi (h + ec + ep ) ] − (Q̇neto − Ẇneto ) = 0
dt vc i
i=1

Si el proceso es de descarga, esto es, no hay corrientes de entrada, la ecuación (1-2) se


transforma en:
NCS
dU
[ ] + ∑ [ṁj (h + ec + ep ) ] − (Q̇neto − Ẇneto ) = 0
dt vc j
j=1

La ecuación (1-2) para un sistema cerrado se reduce a:

dU δQ δW
[ ] = − ó ∆U = Q − W
dt vc dt dt

1.4. Ecuación de la energía mecánica.

De la combinación del principio de la conservación de la cantidad de movimiento y de la


energía se deduce la ecuación de la energía mecánica, la cual está dada por:

P2
u̅2 gz
∫ vdP + ∆ ( ) + ∆ ( ) + weje + lw = 0 (1 − 4)
2g c gc
P1

Donde lw es el trabajo perdido debido a las irreversibilidades.

La forma más simple de la ecuación de energía mecánica es la ecuación de Bernoulli


aplicada a fluidos incompresibles:
Fundamentos de Termodinámica

∆P ∆u̅2 g∆z
( )+( )+( )=0
ρ 2g c gc

La ecuación de energía mecánica se usa frecuentemente en la solución de sistemas de


flujo en tuberías, conductos y accesorios como válvulas, codos, etc.

1.5. Segunda ley de la termodinámica.

Con la primera ley se estima la cantidad de la energía que se transforma en un proceso


determinado, pero no limita la transformación de calor en trabajo útil.

La importancia de la segunda ley consiste principalmente en:


a) restringe el uso de la primera ley porque establece criterios para determinar la
calidad de la energía y la energía disponible, es decir, la energía que puede
transformarse en trabajo útil
b) determina la dirección de los procesos
c) establece cuando un proceso es posible, espontáneo o está en equilibrio
termodinámico
d) proporciona los criterios para determinar la idealidad de un proceso
e) permite evaluar la eficiencia de un proceso.

Como introducción a la segunda ley se definen los conceptos de depósito térmico y


máquinas cíclicas.

Por definición, un depósito térmico es aquel espacio de dimensiones relativamente


grandes con respecto al sistema que intercambia calor con éste. La característica más
importante de un depósito térmico es que sufre procesos isotérmicos e internamente
reversibles. Un depósito térmico es una fuente térmica cuando suministra calor al sistema
y es un sumidero térmico cuando recibe calor del sistema.

Por otra parte, una máquina térmica es aquel proceso cíclico cuya finalidad es producir
trabajo a partir de una transferencia de calor desde al menos un depósito térmico a alta
temperatura hacia al menos otro depósito térmico a baja temperatura. La eficiencia de
una máquina térmica se define como el cociente entre el trabajo producido y el calor
suministrado desde los depósitos térmicos a alta temperatura o fuentes, es decir:

Wneto ∑ Qfuentes − ∑ Qsumideros ∑ Qsumideros


ηt = = =1−
Qneto ∑ Qfuentes ∑ Qfuentes

El calor que se transfiere hacia el o los depósitos térmicos a baja temperatura o sumideros
limita la conversión de la totalidad del calor suministrado desde las fuentes en trabajo
neto producido. Esto quiere decir que no es posible operar una máquina térmica con una
eficiencia de 100%.
Fundamentos de Termodinámica

Un refrigerador es aquel proceso cíclico cuya finalidad es mantener un espacio a


temperatura por debajo de la del entorno (refrigeración) mediante la transferencia de calor
desde ese espacio refrigerado hacia otro depósito a temperatura alta (entorno) y el
suministro de trabajo desde el entorno. El coeficiente de rendimiento de un refrigerador
(COP) se define como el cociente entre el calor extraído (Qrefrig) desde el depósito a
temperatura baja y el trabajo consumido, es decir:

|Qrefrig | |Qrefrig | 1
(COP)REF = = =
|Wsum | |Qentorno | − |Qrefrig | |Qentorno ⁄Qrefrig | − 1

Téngase en cuenta que el entorno está a una temperatura por encima de la del espacio
a refrigerar.

Una bomba de calor es aquel proceso cíclico cuya finalidad es mantener un depósito a
temperatura por encima de la del entorno (calefacción) mediante la transferencia de calor
desde el entorno hacia aquel depósito a temperatura alta y el suministro de trabajo desde
el entorno. El coeficiente de rendimiento de una bomba de calor (COP) se define como
el cociente entre el calor suministrado al depósito (Qcalef) a temperatura alta y el trabajo
consumido, es decir:

|Qcalef | |Qcalef | 1
(COP)BC = = =
|Wsum | |Qcalef | − |Qentorno | 1 − |Qentorno ⁄Qcalef |

Téngase en cuenta que el entorno está a una temperatura por debajo de la del espacio a
calentar.

La primera ley aplicada a una máquina térmica o un refrigerador reversible o irreversible


que intercambia calor con una fuente térmica y un sumidero térmico, está expresada por:

∮ δQ − ∮ δW = ∮ dEsist = 0

Pero ∮ δQ = Qsum − Qelim ; ∮ δW = Wneto

Reemplazando en la ecuación de la primera ley queda:

Qsum − Qelim = Wneto

El resultado anterior indica que en una máquina térmica no es posible convertir el calor
suministrado en trabajo neto.

Los esquemas de una máquina térmica y un refrigerador que intercambian calor con una
fuente térmica y un sumidero térmico, son los siguientes:
Fundamentos de Termodinámica

DEPOSITO TERMICO DEPOSITO TERMICO A


A TEMPERATURA TEMPERATURA
ALTA ALTA
O FUENTE O SUMIDERO

MÁQUINA TÉRMICA Wneto REFRIGERADOR Wneto

DEPÓSITO TERMICO DEPOSITO TERMICO A


A TEMPERATURA BAJA TEMPERATURA BAJA O
O SUMIDERO FUENTE

Los siguientes enunciados expresan las limitaciones de las máquinas térmicas y los
refrigeradores.

1- El enunciado de Kelvin y Planck enmarca la limitación de la primera ley:

Es imposible construir una máquina térmica cuyo único efecto sea


el intercambio de calor con una sola fuente inicialmente en equilibrio
y la producción de trabajo neto.

De acuerdo con este enunciado es imposible alcanzar una eficiencia térmica del 100%
en la transformación de calor suministrado en trabajo neto.

2- El enunciado de Clausius establece que:

Es imposible operar una máquina cíclica de tal manera que el único


efecto externo a la máquina sea la transferencia de calor desde un
depósito a otro depósito a una temperatura mayor.

Así, es imposible transferir calor desde una temperatura menor a una mayor a menos que
se suministre trabajo al sistema.

3- Teorema de Carnot dice:

La eficiencia térmica de una máquina térmica irreversible que funcione entre


dos temperaturas, es siempre menor que la de una máquina de Carnot
(reversible) que funcione entre las mismas temperaturas

4- La desigualdad de Clausius establece que en un proceso cíclico:


Fundamentos de Termodinámica

La integral cíclica de la cantidad Q/T es menor o igual a cero

δQ
∮ ≤0
T

La igualdad se cumple para procesos cíclicos reversibles y la desigualdad para procesos


cíclicos irreversibles. Un proceso cíclico que no cumple con la desigualdad de Clausius
es un proceso imposible.

1.5.1. Procesos reversibles. Un proceso reversible es aquel que al retornar a su estado


inicial no deja consecuencias en los alrededores.

La máquina de Carnot es un ejemplo de máquina térmica reversible, por lo tanto,


cumple la igualdad de Clausius.

Reemplazando las dos etapas no adiabáticas de la máquina térmica de Carnot, en la


ecuación de Clausius, ésta se convierte en:
2 1
δQ δQH δQL
∮ =∫ ( ) +∫ ( ) =0
T 1 TH fuente 2 TL sumidero

La primera integral del término derecho representa el proceso de intercambio de calor QH


con el depósito térmico a temperatura alta TH o fuente, y la segunda, el proceso de
intercambio de calor QL con el depósito térmico a temperatura baja TL o sumidero.

Debido a que los procesos en los depósitos térmicos son reversibles e isotérmicos, las
integrales de la ecuación anterior se resuelven y se obtiene:

QH QL
( ) +( ) =0
TH fuente TL sumidero

Para la máquina térmica, QL = −|QL |

|QH | |QL |
( ) =( )
TH fuente TL sumidero

El resultado anterior es una característica de la máquina de Carnot. Esta característica


de los procesos reversibles es necesaria para la definición de la propiedad termodinámica
denominada entropía.

Considere un ciclo que consta de dos etapas internamente reversibles A y B. En este


caso el proceso cumple la igualdad de la integral cíclica de Q/T:
Fundamentos de Termodinámica

2 1
δQ δ𝑄 δQ 2
∮ =∫ ( ) +∫ ( ) =0
T 1 T A 2 T B A

2 2
δ𝑄 δQ
∫ ( ) =∫ ( ) B
1 T A 1 T B 1

Este resultado permite ver que el valor de la integral es independiente de la trayectoria A


o B o cualquiera otra y solo depende de los estados inicial y final; por lo tanto, esta integral
representa el cambio de una propiedad termodinámica o variable de estado. Esta
propiedad es la entropía designada por S y se define como:
2
δQ δQ
dS = ( ) ∆S = S2 − S1 = ∫ ( ) (A)
T reversible 1 T reversible

1.5.2. Procesos irreversibles. Un proceso irreversible, contrario a uno reversible, deja


consecuencias en los alrededores para que pueda retornar a su punto de partida.

Considérese ahora un ciclo irreversible donde la trayectoria A es irreversible mientras que


la B es reversible. La integral cíclica de Q/T es menor que cero de acuerdo a la ecuación
de Clausius, de tal manera que:
2 1 2
δQ δQ δQ δQ
∮ =∫ ( ) +∫ ( ) =∫ ( ) + (S1 − S2 ) < 0
T 1 T A,irrev 2 T B,rev 1 T A,irrev

2
δQ
∆S = (S2 − S1 ) > ∫ ( ) (B)
1 T irrev

Las ecuaciones (A) y (B) representan la segunda ley de la Termodinámica. La


combinación de estas ecuaciones da el cambio de entropía de cualquier proceso es:
2
δQ δQ
∆S = (S2 − S1 ) ≥ ∫ ( ) ó dS ≥ ( ) (C)
1 T T

Una máquina térmica puede operar si satisface la siguiente desigualdad:

QH QL
( ) +( ) ≤ 0
TH fuente TL sumidero

En general:
k=DT
Qk
∑ ≤0 (D)
Tk
k=1
Fundamentos de Termodinámica

La ecuación (D) es una expresión algebraica, de manera que debe tenerse en cuenta el
signo del calor de cada depósito térmico.
Ejemplo 1.1. Una máquina térmica intercambia calor con un depósito (fuente) a 1000 K y un depósito
(sumidero) a 300 K. Desde la fuente se suministra 2000 kJ de calor (Q H). Calcule la eficiencia térmica, el
trabajo producido y el calor eliminado al sumidero, en kJ, si:
a) la máquina opera como una máquina de Carnot,
b) la máquina es irreversible y su eficiencia térmica es el 75% de la eficiencia de la máquina de Carnot.
Demuestre si cada máquina satisface la desigualdad de Clausius.

Solución:
QH = 2000 kJ; TH = 1000 K; TL = 300 K
Balance de energía W neto = QH – QL;
La eficiencia térmica está definida por t = 1 – QL/QH = W neto/QH
a) máquina de Carnot: De la ecuación (D), (tCarnot) = 1 – QL/QH = 1 – TL/TH

t = 1 – QL/QH = 1 – TL/TH = W neto/QH = 1 – 300/1000 = 0.7;

W neto = t *QH = 0.7 *2000 = 1400 kJ; QL = 600 kJ.

Desigualdad de Clausius: (QH/TH) + (QL/TL) = (2000/1000) + (-600/300) = 0 (Ok)

b) máquina irreversible: t = 0.75*(tCarnot) = 0.75*0.7 = 0.525 = 1 – (QL/QH);

(QL/QH) = 0.475; QL = 950 kJ; W neto = 1050 kJ;

Desigualdad de Clausius: (QH/TH) + (QL/TL) = (2000/1000) + (-950/300) = -1.167 < 0 (Ok)

Los resultados demuestran que ambas máquinas cumplen con la desigualdad de Clausius, por lo tanto, el
proceso es posible.

1.5.3. Principio del Incremento de la entropía del universo. Por definición, un sistema
es una porción del universo que se somete a un análisis termodinámico. Los alrededores
junto con el sistema forman un sistema compuesto adiabático que es el universo. El
cambio de entropía del universo, según la ecuación (C), es mayor o igual cero. Esto
significa que la entropía del universo no se conserva, sino que se incrementa debido a
las irreversibilidades de los procesos reales, es decir:

(∆S)universo = Sgen = (∆S)sistema + (∆S)alrededores ≥ 0 (1 − 5)

Los procesos ideales son aquellos que no incrementan entropía del universo y suelen
compararse con los procesos reales bajo las mismas condiciones para medir la eficiencia
de éstos. Así, la eficiencia de un proceso, proceso se define como:

Energía útil esperada en el proceso real


ηproceso =
Energía útil esperada en el proceso si fuera reversible

Considerando el mismo volumen de control de la primera ley, la segunda ley establece:


Fundamentos de Termodinámica

Velocidad de velocidad de flujo velocidad de flujo velocidad de


Velocidad de
acumulación de entropía que sale de entropía que entra flujo de entropía
generación de
{ } = de entropía en el + con la masa − con la masa + debido a la transferencia ≥ 0
entropía en
volumen de desde el volumen hacia el volumen de calor con los
el universo { } { } { } { }
control de control de control alrededores

La expresión de la segunda ley, conforme a la ecuación anterior, es la siguiente:

dS Q̇k
Ṡgen = ( ) + ∑[ṁj sj ] − ∑[ṁi si ]ent + ∑ ≥0 (1 − 6)
dt vc sal Tk
j=1 i=1 k=1

Los tres primeros términos de la derecha de la ecuación (1-6) representan el cambio de


entropía en el sistema, siendo positivo cuando el proceso es de calentamiento o
evaporación y negativo para enfriamiento o condensación; el último término es la suma
del cambio de entropía de los depósitos térmicos, incluyendo los alrededores, con los que
el sistema intercambia calor. Obsérvese que la producción o el consumo de trabajo no
están asociados a la generación de entropía.

La suma de los cambios de entropía de todos los depósitos térmicos incluyendo los
alrededores está dada por:

Q̇k Q̇0 Q̇n


∑ = +∑
Tk T0 Tn
k=1 n=1

En forma general, para procesos reales la ecuación (1-6) indica que la entropía del
universo siempre se incrementa y no es conservativa. Esto representa el principio del
incremento de la entropía del universo.

1.6. Análisis exergético (de disponibilidad) de procesos.

La combinación de la primera y segunda leyes de la Termodinámica permite encontrar la


potencia o el trabajo útil en un proceso en función del cambio de la energía, el cambio de
la entropía y la entropía generada. Para un proceso reversible, esta potencia o trabajo
útil sería máximo si se produce o mínimo si se suministra.

Cuando este proceso se lleva a cabo reversiblemente desde el estado inicial del sistema
hasta el estado de los alrededores a P0 y T0, la potencia útil se denomina exergía. La
exergía se define como la máxima potencia o el trabajo máximo (o el mínimo) que se
produciría o se consumiría si se llevara a cabo dicho proceso reversiblemente; así, la
exergía es la disponibilidad o energía disponible que conforma la totalidad de la
energía de un sistema:

{Energía total} = {Energía disponible} + {Energía no disponible}


Fundamentos de Termodinámica

La energía no disponible es aquella parte de la energía total que no puede convertirse en


trabajo útil y representa el trabajo perdido o potencial destruido durante el proceso. Este
trabajo perdido se conoce también como la irreversibilidad del proceso, I, de modo que:

Ẇperdida = Ẇreversible − Ẇreal útil = İ

El estado de los alrededores a P0 y T0 se define como el estado muerto en donde el


sistema no tiene exergía o capacidad de realizar un trabajo. El estado de los alrededores
con energía cinética y potencial despreciables, se puede considerar como el estado
muerto.
Fundamentos de Termodinámica

2. PROPIEDADES TERMODINÁMICAS DE UNA SUSTANCIA PURA

Las propiedades termodinámicas específicas intrínsecas que dependen de la naturaleza


química de la materia son la presión (P), la temperatura (T), el volumen (v), la energía
interna (u) y la entropía (s).

La entalpía (h), la energía de Gibbs (g) y la energía de Helmhotz (a) son propiedades
termodinámicas por definición en función de las propiedades intrínsecas.
Asi:
h = u + pv; g = h − Ts; a = u − Ts

Otras propiedades derivadas de las anteriores son:


El calor específico a presión constante, el calor específico a volumen constante, el
coeficiente de Joule-Thompson, el coeficiente de expansión volumétrica y el coeficiente
de compresibilidad isotérmica

∂h ∂u ∂T 1 ∂v 1 ∂v
Cp = ( ) ; Cv = ( ) ; μJT = ( ) ; β= ( ) ; k=− ( )
∂T P ∂T v ∂P h v ∂T P v ∂P T

Se puede observar que todas las propiedades son funciones de dos propiedades
intensivas independientes, de acuerdo con la regla de las fases de Gibbs. Si estas dos
son la presión y la temperatura, la determinación del volumen, la entalpía y la entropía,
permite calcular las demás.

v = v(P, T); h = h(P, T); s = s(P, T)

La representación del comportamiento de estas propiedades con respecto a la presión y


temperatura da origen a los diagramas de fases Presión-volumen (P-v), Presión-entalpía
(P-h) y Presión-entropía (P-s) a temperatura constante, a los diagramas Temperatura-
volumen (T-v), Temperatura-entalpía (T-h) y Temperatura-entropía (T-s) a presión
constante y al diagrama Presión-Temperatura (P-T) a volumen constante.
Fundamentos de Termodinámica

Diagrama P-v a temperatura constante Diagrama T-v a presión constante

Diagrama Presión-entalpía a temperatura constante Diagrama Temperatura-entropía a presión constante

En los diagramas de arriba se puede observar dos curvas cuyo punto máximo es el punto
crítico. Este punto crítico es la condición en la que por encima de él no es posible
distinguir la fase líquida de la fase vapor. La curva que está a la izquierda representa los
estados de líquido saturado y las propiedades se denotan por el subíndice f; la curva
derecha representa los estados de vapor saturado y las propiedades se denotan con el
subíndice g. La región que está a la izquierda de la curva de líquido saturado es la región
líquida comprimida o sub-enfriada. La región que está a la derecha de la curva de vapor
saturado es la región vapor sobrecalentado.

La presión que corresponde al líquido saturado y al vapor saturado en equilibrio a una


determinada temperatura, se denomina presión de saturación o presión de vapor. A
su vez, a una determinada presión le corresponde su temperatura de saturación.

P sat = P sat ( T) ó T sat = T sat ( P)

Propiedades del vapor saturado: a = ag (P sat , T) ó a = ag (P, T sat )


Propiedades del líquido saturado: a = af (P sat , T) ó a = af (P, T sat )
Fundamentos de Termodinámica

La región comprendida entre estas dos curvas representa estados de mezclas en


equilibrio del líquido saturado y el vapor saturado.

Considerando una mezcla de una masa mf de líquido saturado y una masa mg de vapor
saturado, la relación entre la masa del vapor y la masa de la mezcla, se denomina calidad
del vapor (x),
mg
x=( ) (2 − 1)
mf + mg

Las propiedades de la mezcla líquido-vapor se determinan a partir de las del líquido


saturado y del vapor saturado, mediante la siguiente ecuación:

mf mg
a=( ) af + ( ) a = (1 − x)af + xag = af + x(ag − af ) (2 − 2)
mf + mg mf + mg g

Donde la propiedad a puede ser volumen, energía interna, entalpía, entropía o cualquiera
propiedad termodinámica intensiva.

Los criterios para establecer los estados termodinámicos de una sustancia pura son los
siguientes:

1. Líquido comprimido o sub-enfriado: P > P sat ; T < T sat ; a < af


2. Líquido saturado: P = P sat ; T = T sat ; a = af ; x = 0
3. Mezcla líquido-vapor: P = P sat ; T = T sat ; af < a < ag ; 0 < x < 1
4. Vapor saturado: P = P sat ; T = T sat ; a = ag ; x = 1
5. Vapor sobrecalentado: P < P sat ; T > T sat ; a > ag

El comportamiento PvT, PhT y PsT se representa a través de los diagramas o de tablas


denominadas tablas de vapor de saturación y tablas de líquido comprimido o vapor
sobrecalentado. Las tablas de saturación a su vez se dividen en tablas de temperatura y
tablas de presión.

Tabla de temperatura del vapor de agua (presión de saturación)


T Psat vf vg uf ug hf hg sf sg
(° C) (kPa) (m3/kg) (m3/kg) (kJ/kg) (kJ/kg) (kJ/kg) (kJ/kg) (kJ/kg-K) (kJ/kg)
120 198,67 0,001060 0,8912 503,601 2528,86 503,812 2705,93 1,5279 7,1291
125 232,24 0,001065 0,7700 524,827 2534,27 525,074 2713,10 1,5816 7,0770

Tabla de presión del vapor de agua (temperatura de saturación)


P Tsat vf vg uf ug hf hg sf sg
(kPa) (° C) (m3/kg) (m3/kg) (kJ/kg) (kJ/kg) (kJ/kg) (kJ/kg) (kJ/kg-K) (kJ/kg-K)
190 118,6 0,001059 0,9292 497,649 2527,32 497,85 2703,88 1,5127 7,1440
200 120,21 0,001060 0,8857 504,492 2529,09 504,70 2706,23 1,5302 7,1269
Fundamentos de Termodinámica

Tabla de líquido comprimido o sub-enfriado


P = 200 kPa (Tsat = 120,91° C) P = 300 kPa (Tsat = 133,52° C)
T V U h s v u h s
(° C) (m3/kg) (kJ/kg) (kJ/kg) (kJ/kg-K) (m3/kg) (kJ/kg) (kJ/kg) (kJ/kg-K)
50 0,001012 209,301 209,504 0,7037 0,001012 209,286 209,59 0,7037
100 0,001043 419,032 419,24 1,3071 0,001043 419,002 419,32 1,3071

Tabla de vapor sobrecalentado


P = 200 kPa (Tsat = 120,91° C) P = 300 kPa (Tsat = 133,52° C)
T V U h s v u h s
(° C) (m3/kg) (kJ/kg) (kJ/kg) (kJ/kg-K) (m3/kg) (kJ/kg) (kJ/kg) (kJ/kg-K)
135 0,9228 2553,27 2737,84 7,2058 0,6083 2545,72 2728,22 6,9998
140 0,9352 2561,28 2748,33 7,2313 0,6170 2554,27 2739,36 7,0269
Fuente: termodinamicaparaiq.blogspot.com/p/los-programas-que-se-muestran.html

En la web se puede descargar software gratuito para la determinación de las propiedades


termodinámicas del vapor.

Esta aplicación tiene tres pestañas. Una pestaña de estados de saturación donde se especifica la temperatura o la
presión y la calidad entre cero y uno. La segunda pestaña es de estados de vapor sobrecalentado o líquido
comprimido/sub-enfriado en donde se puede especificar 1) presión y temperatura; 2) temperatura o presión y cualquiera
propiedad (v, u, h, s). En la tercera pestaña aparecen las constantes del agua como las constantes críticas, masa
molecular, temperatura y presión del punto triple.
Fundamentos de Termodinámica

Pantallazos de un software gratuito de tablas de vapor del agua

Ejemplo.
Un sistema cilindro-émbolo conectado a un resorte lineal contiene 2,5 kg de agua a 150 kPa y calidad de
50%. El resorte apenas toca el émbolo de tal manera que no se realiza ningún trabajo. Posteriormente se
aplica calor al cilindro-émbolo hasta que el volumen y la presión se duplican.
Determine:
a) El trabajo sobre el resorte, en kJ
b) El calor suministrado, en kJ

Solución:
La fuerza del resorte lineal es directamente proporcional al desplazamiento: Fresorte = k resorte ∆X
La presión total que ejerce el fluido es la suma de la presión sobre el émbolo más la presión sobre el
resorte.
Fresorte k resorte ∆X k resorte (X − X1 )
P = P1 + Presorte = P1 + = P1 + = P1 +
A A A

En donde A es el área del émbolo, X es el desplazamiento.


El volumen está relacionado con el desplazamiento: V = V1 + A ∗ (X − X1 )
k(V−V1 ) k (P −P1 )
La presión queda expresada en términos del volumen: P = P1 + 2 , donde = (V2
A A2 2 −V1 )
El trabajo de frontera (expansión sobre los alrededores y sobre el resorte) está dado por:

V2 V2 V2 (V2 −V1 )
k(X − X1 ) k(V − V1 )
W = ∫ PdV = ∫ [P1 + ] dV = ∫ P1 dV + ∫ d(V − V1 )
V1 V1 A V1 0 A2

k(V2 − V1 )2 (P2 − P1 )(V2 − V1 )


W = P1 (V2 − V1 ) + 2
= P1 (V2 − V1 ) +
2A 2
(P2 − P1 )(V2 − V1 )
Wresorte =
2

Balance de energía en un sistema cerrado: Q − W = ∆U


Estado 1:
P1 = 150 kPa, Calidad, x = 0,50; m = 2,5 kg
Fundamentos de Termodinámica

v1 = vf + x*(vg – vf) = 0,58017 m3/kg; V1 = m*v1 = 1,450425 m3


u1 = uf + x*(ug – uf) = 1493,09 kJ/kg; U1 = m*u1 = 3732,73 kJ

Estado 2:
P2 = 2*P1 = 300 kPa, m = 2,5 kg; V2 = 2*V1 = 2,90085 m3
v2 = V2/m = 2*v1 = 1,16034 m 3/kg;
A 300 kPa, vg = 0,6058 m3/kg; por lo tanto, el estado 2 es vapor sobrecalentado
Por interpolación, u2 = 3102,13 kJ/kg; U2 = m*u2 = 7755,33 kJ

Trabajo total, W = 150*(2,90085-1,450425)+(300-150)*( 2,90085-1,450425)/2 = 326,35 kJ

Trabajo sobre el resorte: Wresorte = (300-150)*(2,90085-1,450425)/2 = 108,78 kJ

Cambio de la energía interna: ∆U = 7755,33 − 3732,73 = 4022,6 kJ

Calor necesario: 𝐐 = 𝐖 + ∆𝐔 = 𝟒𝟑𝟒𝟖, 𝟗𝟓 𝐤𝐉

Conclusiones:
Del total del trabajo realizado por el sistema, el trabajo sobre el resorte es el trabajo útil y representa apenas
el 2,5% del calor suministrado (W resorte/Q), mientras que el 92,5% del calor es entregado al fluido para
aumentar su energía interna (U/Q). El resto es entregado al émbolo.

Distribución del calor suministrado:


Aumento de la energía interna ∆U = 4022,60 kJ (92,5%)
Trabajo total: Wtotal = 326,35 kJ (7,5%)
Trabajo sobre el resorte (útil) Wresorte = 108,78 kJ (2,5%)

Ejemplo.
Establezca las condiciones de los estados termodinámicos del agua de acuerdo a las especificaciones
dadas (en negrilla) en la tabla
P T Psat Tsat Calidad v u h s
Condición
(kPa) (° C) (kPa) (° C) (%) (m3/kg) kJ/kg) (kJ/kg) (kJ/kg-K)
1000 75 38,59 179,88 ----- Liq comp 0,00103 313,78 314,81 1,0152
800 250 3976,2 170,41 ----- Vap sob 0,29320 2715,87 2950,44 7,0401
400 143,61 400 143,61 75 Mezcla 0,34706 2065,88 2204,70 5,6157
550 175 892,6 155,46 ----- Vap sob 0,36177 2599,04 2798,01 6,8928
198,67 120 198,67 120 100 Vap sat 0,89121 2528,86 2705,93 7,1291
2500 223,95 2500 223,95 0 Liq sat 0,00120 958,91 961,91 2,5543
3000 233,85 3000 233,85 91,49 Mezcla 0,06110 2467,16 2650,45 5,8844
12,35 50 12,35 50 25,41 Mezcla 3,05633 776,75 814,50 2,5765
1500 131,05 1500 131,05 ----- Liq comp 0,00107 550,07 551,68 1,6445
4500 600 ------ 257,44 ----- Vap sob 0,08766 3276,38 3670,86 7,3127

Ejemplo.
Vapor de agua a 220 kPa y 250° C fluye a través de una tubería de 2 pulg (5,08 cm) de diámetro interno
con un flujo volumétrico de 2,5 m 3/s saliendo en dos corrientes, cada una a 210 kPa y 150° C. Una de las
ramas tiene un diámetro interno de 1 pulg (2,54 cm) y la otra tiene un diámetro interno de ½ pulg (1,27 cm).
Si la velocidad de la rama de ½ pulg es igual de la velocidad de la rama de 1 pulg, determine: a) el flujo
másico de cada corriente; b) el cambio de entalpía del agua; c) el cambio de entropía del agua.

2
V2 = ¿? m3/s
P2 = 210 kPa
T2 = 150° C
2 = 2,54 cm
1
V1 = 1,5 m3/s
P1 = 220 kPa 3
T1 = 250° C V3 = ¿? m3/s
1 = 5,08 cm P3 = 210 kPa
T3 = 150° C
3 = 1,27 cm
Fundamentos de Termodinámica

Solución:
El flujo másico está relacionado con el flujo volumétrico y el volumen específico
V̇ Au̅
ṁ = =
v v
Estado 1:
A 220 kPa y 250° C, v1 = 1,08914 m3/kg, h1 = 2970,55 kJ/kg; s1 = 7,665 kJ/kg-K
V̇1 A1 u̅1 1,5
ṁ1 = = = = 1,3772 kg/s
v1 v1 1,08914
Estado 2:
A 210 kPa y 150° C, v2 = 0,91332 m3/kg, h2 = 2768,33 kJ/kg; s2 = 7,257 kJ/kg-K;
Área de flujo: A2 = (/4)*(2)2 = (/4)*(23)2 = 4A3 = 5,067*10-4 m2
Estado 3:
A 210 kPa y 150° C, v3 = 0,91332 m3/kg, h3 = 2768,33 kJ/kg; s3 = 7,257 kJ/kg-K;
Área de flujo: A3 = (/4)*(3)2 = 1,267*10-4 m2

Balance de materia:
A2 u̅2 A3 u̅3 4A3∗ u̅3 A3 u̅3 5A3 u̅3
ṁ1 = ṁ2 + ṁ3 = + = + = = 5ṁ3
v2 v3 v3 v3 v3
3
1,3772 kg m kg m3
ṁ3 = = 0,2754 ; V̇3 = 0,2516 ; ṁ2 = 1,1018 ; V̇2 = 1,0063
5 s s s s

Cambio de entalpía:
∆Ḣ = (Ḣ2 + Ḣ3 ) − Ḣ1 = (ṁ2 h2 + ṁ3 h3 ) − ṁ1 h1 = (ṁ2 + ṁ3 )h2 − ṁ1 h1 = ṁ1 (h2 − h1 )
∆Ḣ = 1,3772 ∗ (2768,33 − 2970,55) = −278,50 kW

Cambio de entropía:
∆Ṡ = (Ṡ2 + Ṡ3 ) − Ṡ1 = (ṁ2 s2 + ṁ3 s3 ) − ṁ1 s1 = (ṁ2 + ṁ3 )s2 − ṁ1 s1 = ṁ1 (s2 − s1 )
kW
∆Ṡ = 1,3772 ∗ (7,257 − 7,665) = −0,5619
K

Ejemplo.
Un cilindro aislado está dividido en dos compartimientos de 10 pies3 cada uno mediante un pistón
inicialmente bloqueado, tal como se muestra en la figura de abajo. El compartimiento A tiene aire a 2 atm
y 600 R y el compartimiento B contiene aire a 10 atm y 2000 R. Ahora el pistón se desbloquea de tal manera
que queda libre para moverse y transferir calor de modo que la temperatura final es uniforme para los dos
compartimientos, es decir T A = TB. Determine: a) la masa de cada compartimiento; b) la temperatura y
presión final; c) el cambio de entropía del aire
Fuente bibliográfica: Tablas del aire en Cengel (A-17E)

Solución:
Constante específica del aire: R = 0,3704 psia-pies3/(lbm-R) = 0,06855 BTU/(lbm-R)

En la siguiente tabla aparecen las propiedades calculadas y obtenidas de las tablas A-17E del Cengel
Las propiedades en negrillas son las especificaciones, las demás fueron calculadas.
Fundamentos de Termodinámica

Condiciones iniciales Condiciones finales Final total


A B A B
T (R) 600 2000 1475,18 1475,18 1475,18
P (psia) 29,4 147 90,355 90,355 90,355
u (BTU/lbm) 102,34 367,61 261,50 261,50
s°T (BTU/lb-R) 0,62607 0,93205 0,84975 0,84975
V (pies3) 10 10 8 12 20
m (lb) 1,323 1,984 1,323 1,984 3,307

La masa de cada compartimiento se encuentra con la ecuación de los gases ideales: m = PV⁄(RT)
29,4 ∗ 10
mA = P1A V1A ⁄(RT1A ) = = 1,323 lb
0,3704 ∗ 600
147 ∗ 10
mB = P1B V1B ⁄(RT1B ) = = 1,984 lb
0,3704 ∗ 2000

Masa total: m = 1,323 + 1,984 = 3,307 lb. Volumen total: V = 20 pies 3

Balance de energía:
mA (u2 − u1A ) + mB (u2 − u1B ) = 0;
Despejando u2:
u2 = (mA u1A + mB u1B )⁄(mA + mB ) = (1,323 ∗ 102,34 + 1,984 ∗ 367,61)⁄(1,323 + 1,984) = 261,50 BTU/lb

Usando las tablas A-17E se interpola con u2 para hallar 𝐓𝟐 = 𝟏𝟒𝟕𝟓, 𝟏𝟖 𝐑 y s°T2 = 0,84975 BTU/lb-R

De la ecuación de los gases ideales:

𝐦𝐑𝐓𝟐 𝟑, 𝟑𝟎𝟕 ∗ 𝟎, 𝟑𝟕𝟎𝟒 ∗ 𝟏𝟒𝟕𝟓, 𝟏𝟖


𝐏𝟐 = = = 𝟗𝟎, 𝟑𝟓𝟓 𝐩𝐬𝐢𝐚
𝐕 𝟐𝟎

El volumen final de cada compartimiento se halla con:

mA RT2 1,323 ∗ 0,3704 ∗ 1475,18


V2A = = = 8 pies 3 ; V2B = 12 pies 3
P2 90,355

Balance de entropía:
Sgen = mA (s2 − s1A ) + mB (s2 − s1B )

El cambio de entropía de cada compartimiento se expresa en términos de la temperatura y la presión:


P2A BTU
s2 − s1A = sT° 2 − sT° 1A − Rln ( ) = 0,84975 − 0,62607 − 0,06855 ∗ ln(90,355⁄29,4) = 0,14672
P1A lb − R
P2B BTU
s2 − s1B = sT° 2 − sT° 1B − Rln ( ) = 0,84975 − 0,93205 − 0,06855 ∗ ln(90,355⁄147) = −0,04894
P1B lb − R

Finalmente,
𝐒𝐠𝐞𝐧 = 𝟏, 𝟑𝟐𝟑 ∗ 𝟎, 𝟏𝟒𝟔𝟕𝟐 + 𝟏, 𝟗𝟖𝟒 ∗ (−𝟎, 𝟎𝟒𝟖𝟗𝟒) = 𝟎, 𝟎𝟗𝟕𝟎 𝐁𝐓𝐔/𝐑

Obsérvese que T1A < T2 < T1B ; P1A < P2 < P1B

Igualmente ∆TA > 0; ∆TB < 0; ∆PA > 0; ∆PB < 0; ∆VA < 0; ∆VB > 0
Fundamentos de Termodinámica

3. RELACIONES TERMODINÁMICAS GENERALIZADAS

A partir de la combinación de la primera y segunda leyes en un sistema cerrado, de la


definición de la entalpía, la energía de Gibbs (g = h – Ts) y la energía de Helmhotz (a = u
– Ts), se obtienen las siguientes relaciones termodinámicas generalizadas:

du = Tds − Pdv (3 − 1)

dh = Tds + vdP (3 − 2)

dg = −sdT + vdP (3 − 3)

da = −sdT − Pdv (3 − 4)

Las relaciones Maxwell se encuentran usando el criterio de la exactitud en cada una de


las relaciones anteriores:

∂T ∂P
( ) = −( ) (3 − 5)
∂v s ∂s v

∂T ∂v
( ) =( ) (3 − 6)
∂P s ∂s P

∂s ∂v
−( ) =( ) (3 − 7)
∂P T ∂T P

∂s ∂P
( ) =( ) (3 − 8)
∂v T ∂T v

Algunas relaciones derivadas de las anteriores son:

1. Efecto de la presión a temperatura constante sobre la entalpía,

∂h ∂s ∂v
( ) = T ( ) + v = −T ( ) + v (3 − 9)
∂P T ∂P T ∂T P

Para gases ideales, la entalpía es independiente de la presión, es decir:

∂h
( ) =0
∂P T

2. Efecto del volumen a temperatura constante sobre la energía interna,


Fundamentos de Termodinámica

∂u ∂s ∂P
( ) = T( ) − P = T( ) − P (3 − 10)
∂v T ∂v T ∂T v

Para gases ideales, la energía interna es independiente del volumen, es decir:

∂u
( ) =0
∂v T

3. Variación de la entropía en función de la temperatura y la presión,

Cp dT v
ds = − dP (3 − 11)
T T

4. Variación de la entropía en función de la temperatura y el volumen,

Cv dT P
ds = + dv (3 − 12)
T T

5. Coeficiente de Joule-Thompson, variación de la temperatura con respecto a la


presión a entalpía constante (expansión isentálpica),

∂h ∂v
∂T ( ) T( ) − v
∂P T ∂T P
μJT =( ) =− = (3 − 13)
∂P h ∂h CP
( )
∂T P

Por ejemplo, para procesos de refrigeración, el coeficiente de Joule-Thompson del


refrigerante es positivo, es decir, su expansión isentálpica causa una disminución de la
temperatura; sin embargo, la expansión isentálpica de un gas ideal es también isotérmica,
porque su coeficiente de Joule-Thompson es igual a cero.

6. Coeficiente de expansión isentrópica, variación de la temperatura con respecto a


la presión a entropía constante,

∂s ∂v
∂T ( ) T( )
∂P T ∂T P
μS = ( ) = − = (3 − 14)
∂P s ∂s CP
( )
∂T P

Por ejemplo, en turbina o compresores, el coeficiente de expansión isentrópica del fluido,


siempre es positivo, es decir, su expansión isentrópica causa una disminución de la
temperatura.
Fundamentos de Termodinámica

4. PROPIEDADES TERMODINÁMICAS DE GASES IDEALES PUROS

Las propiedades termodinámicas de los gases ideales puros se basan en la ecuación de


estado correspondiente:

PV = nR u T
Donde:
n = moles del gas, en kmol
Ru = constante universal de los gases ideales = 8,314 kPa.m3/kmol-K
V = volumen, en m3
P = presión absoluta, en kPa
T = temperatura absoluta, en K

1. Volumen específico. La ecuación para el volumen específico a partir de la ecuación


de estado de los gases ideales es la siguiente:

RT
v= (4 − 1)
P

Donde:
v = volumen específico, en m3/kg
R = constante específica del gas = Ru/PM, en kPa.m3/kg-K
PM = masa molecular del gas, en kg/kmol

La capacidad térmica a presión constante, CP, se define como:

∂h
CP = ( ) = F(T) = a + bT + cT 2 + dT 3
∂T P

Representa el calor que necesita una unidad de masa de una sustancia para aumentar
un grado de temperatura a presión constante.

2. Entalpía específica. La ecuación (3-9) aplicada a un gas ideal da como resultado que
la entalpía de un gas ideal es independiente de la presión y función únicamente de la
temperatura,

∂h ∂v RT
( ) = −T ( ) + v = − +v=0
∂P T ∂T P P
Entonces,
∂h dh
CP = ( ) =
∂T P dT
Fundamentos de Termodinámica

La integración de la ecuación diferencial de la entalpía permite determinar la entalpía en


función de la temperatura:
T
b c d
h − h0 = ∫ CP dT = a(T − T0 ) + (T 2 − T02 ) + (T 3 − T03 ) + (T 4 − T04 ) (4 − 2)
T0 2 3 4

Donde CP está en kJ/kg-K y T en K; a, b, c y d son constantes propias de cada sustancia


como gas ideal; la temperatura T0 y la entalpía h0 tienen valores arbitrarios y se cancelan
en un proceso determinado desde un estado 1 hasta un estado 2.
T2
b c d
h2 − h1 = ∫ CP dT = a(T2 − T1 ) + (T22 − T12 ) + (T23 − T13 ) + (T24 − T14 ) (4 − 2A)
T1 2 3 4

La capacidad térmica a volumen constante, Cv, se define como:

∂u
Cv = ( )
∂T v

Representa el calor que necesita una unidad de masa de una sustancia para aumentar
un grado de temperatura a volumen constante.

3. Energía Interna específica. La ecuación (3-10) aplicada a un gas ideal da como


resultado que la energía interna de un gas ideal es independiente del volumen y función
únicamente de la temperatura,

∂u ∂P RT
( ) = T( ) − P = −P=0
∂v T ∂T v v
Entonces,
∂u du
Cv = ( ) =
∂T v dT

La relación entre Cp y Cv de un gas ideal, se halla a partir de la definición de entalpía:

dh = du + d(Pv) = CP dT = Cv dT + RdT

CP = Cv + R

b c d
u − u0 = a(T − T0 ) + (T 2 − T02 ) + (T 3 − T03 ) + (T 4 − T04 ) − R(T − T0 ) (4 − 3)
2 3 4

El cambio de energía interna en un proceso está dado por:

b c d
u2 − u1 = a(T2 − T1 ) + (T22 − T12 ) + (T23 − T13 ) + (T24 − T14 ) − R(T − T0 ) (4 − 3A)
2 3 4
Fundamentos de Termodinámica

4. Entropía específica. La ecuación (3-11) aplicada a un gas ideal puro permite


determinar la entropía específica, en función de la temperatura y la presión:
T P T P T
CP dT v CP dT R CP dT P
s − s0 = ∫ − ∫ dP = ∫ − ∫ dP = ∫ − Rln ( )
T0 T P0 T T0 T P0 P T0 T P0

T c d P
s − s0 = aln ( ) + b(T − T0 ) + (T 2 − T02 ) + (T 3 − T03 ) − Rln ( ) (4 − 4)
T0 2 3 P0

La temperatura T0, la presión P0 y la entropía s0 tienen valores arbitrarios y se cancelan


en un proceso determinado desde un estado 1 hasta un estado 2.

T2 c d P2
s2 − s1 = aln ( ) + b(T2 − T1 ) + (T22 − T12 ) + (T23 − T13 ) − Rln ( ) (4 − 4A)
T1 2 3 P1

Ejemplo. 20 m3/s de metano entra a un intercambiador de calor a 47° C y 110 kPa y luego sale a 127° C y
105 kPa. Usando la capacidad térmica del metano en términos de la temperatura, determine el flujo
volumétrico en la salida, el cambio de entalpía, en kW y el cambio de entropía, en kW/K.
Datos:
Capacidad térmica, CP = 19,89 + 5,024*10-2T+1,269*10-5T2 -1,101*10-8T3, CP en kJ/kmol-K, T en K
Masa molecular, PM = 18 kg/kmol
Constante específica, R = Ru/PM = 8,314/18 = 0,46189 kPa.m3/kg-K
Solución:
Estado 1 = entrada; Estado 2 = salida
T1 = 47° C = 320 K; P1 = 110 kPa; T2 = 127° C = 400 K; P2 = 105 kPa
El flujo volumétrico está relacionado con el volumen específico y el flujo másico:
V̇1 P1 V̇1 ṁ2 RT2
ṁ 1 = = = ṁ2 = ṁ; V̇2 =
v1 RT1 P2

V̇1 P1 V̇1 110 ∗ 20 kg ṁ2 RT2 14,885 ∗ 0,46189 ∗ 400


ṁ 1 = = = = 14,885 ; V̇2 = = = 26,19 m3 /s
v1 RT1 0,46189 ∗ 320 s P2 105

Cambio de entalpía:
5,024 ∗ 10−2 1,269 ∗ 10−5
h2 − h1 = 19,89(400 − 320) + (4002 − 3202 ) + (4003 − 3203 )
2 3
1,101 ∗ 10−8 kJ kmol
− (4004 − 3204 ) = 3,129 ∗ 103 ∗ = 1,738 ∗ 102 kJ/kg
4 kmol 18 kg

Ḣ2 − Ḣ1 = ṁ(h2 − h1 ) = 2587,11 kW

Cambio de entropía:
400 1,269 ∗ 10−5
s2 − s1 = 19,89 ∗ ln ( ) + 5,024 ∗ 10−2 (400 − 320) + (4002 − 3202 )
320 2
1,101 ∗ 10−8 1,05 kJ kmol kJ
− (4003 − 3203 ) − 8,314 ∗ ln ( ) = 9,095 ∗ = 0,5053
3 1,10 kmol − K 18 kg kg − K

Ṡ2 − Ṡ1 = ṁ(s2 − s1 ) = 7,521 kW/K

Ejemplo. 0,35 kmol/s de CO2 se expande desde 1200 K y 400 kPa (4 bar) hasta 800 K y 100 kPa (1 bar).
Determine el cambio de entalpía, en kW y el cambio de entropía, en kW/K
Fundamentos de Termodinámica

Datos:
Capacidad térmica, CP = 22,26 + 5,981*10-2T -3,501*10-5T2 +7,469*10-9T3, CP en kJ/kmol-K, T en K
Masa molecular, PM = 44 kg/kmol
Solución:
T1 = 1200 K; P1 = 4,0 bar; T2 = 800 K; P2 = 1,0 bar; flujo másico = m 1 = m2 = 0,35*44 = 15,40 kg/s
Cambio de entalpía:
5,981 ∗ 10−2 3,501 ∗ 10−5
h2 − h1 = 22,26(800 − 1200) + (8002 − 12002 ) − (8003 − 12003 )
2 3
7,469 ∗ 10−9 kJ kmol
+ (8004 − 12004 ) = −2,174 ∗ 104 ∗ = −494,19 kJ/kg
4 kmol 44 kg

Ḣ2 − Ḣ1 = ṁ(h2 − h1 ) = −7610,53 kW

Cambio de entropía:
800 3,501 ∗ 10−5
s2 − s1 = 22,26 ∗ ln ( ) + 5,981 ∗ 10−2 (800 − 1200) − (8002 − 12002 ) +
1200 2
7,469 ∗ 10−9 1,0 kJ kmol kJ
(8003 − 12003 ) − 8,314 ∗ ln ( ) = 1,078 ∗ = 0,0245
3 4,0 kmol − K 44 kg kg − K

Ṡ2 − Ṡ1 = ṁ(s2 − s1 ) = 0,377 kW/K

Ejemplo. Un tanque rígido de 1 m 3 contiene aire a 2500 kPa y 300 K conectado mediante una válvula a un
cilindro-émbolo como se muestra en la figura. El pistón de área de 0,1 m 2 requiere 250 kPa para moverse.
Ahora la válvula se abre hasta que el pistón se eleva 2 m lentamente y luego la válvula se cierra. Durante
el proceso la temperatura del aire permanece constante en 300 K. Determine la presión final en el tanque,
en kPa y el calor suministrado, en kJ

Tanque A

Cilindro B

Pistón

Solución:
Sistema: tanque A y cilindro-émbolo B
Estado inicial: 1A y 1B; estado final: 2A y 2B
Proceso en el émbolo: isobárico, P1B = P2B = PB = 250 kPa
Proceso en el tanque: isovolumétrico, V1A = V2A = 1 m3
Estado 1A: V1A = 1 m3; P1A = 2500 kPa; T1A = 300 K
P1A V1A 2500 ∗ 1,0
m1A = = = 29,036 kg
RT1A 0,287 ∗ 300
Estado 2A: V2A = 1 m3; T2A = 300 K; m2A = m1A – m2B
m2A RT2A
P2A =
V2A
Estado 1B: V1B = 0; m1B = 0
Estado 2B: V2B = 0,1m2*2m = 0,2 m3; P2B = 250 kPa; T2B = 300 K
P2B V2B 250 ∗ 0,2
m2B = = = 0,581 kg; m2A = m1A − m2B = 28,455 kg
RT2B 0,287 ∗ 300
Fundamentos de Termodinámica

Trabajo de expansión en el cilindro: 𝐖𝐁 = 𝐏𝐁 (𝐕𝟐𝐁 − 𝐕𝟏𝐁 ) = 𝟐𝟓𝟎 ∗ 𝟎, 𝟐 = 𝟓𝟎 𝐤𝐉


𝟐𝟖, 𝟒𝟓𝟓 ∗ 𝟎, 𝟐𝟖𝟕 ∗ 𝟑𝟎𝟎
𝐏𝟐𝐀 = = 𝟐𝟒𝟓𝟎 𝐤𝐏𝐚
𝟏, 𝟎

Cambio de energía interna total:


∆U = (U2A − U1A ) + (U2B − U1B ) = (m2A u2A − m1A u1A ) + (m2B u2B − m1B u1B )
Pero m1B = 0; u1A = u2A = u2B = u2 = F(T)

∆U = (m2A u2 − m1A u1A ) + m2B u2 = (m2A + m2B )u2 − m1A u1A = (m2A + m2B − m1A )u2 = 0

𝐐 = 𝐖𝐁 + ∆𝐔 = 𝟓𝟎 𝐤𝐉

Ejemplo. Un tanque rígido está dividido en dos compartimientos separados por un pistón buen conductor
del calor y que se puede mover libremente. Inicialmente el compartimiento A contiene 1 kg de N2 a 800 kPa
y 400 K. El compartimiento B tiene un volumen que es doble del compartimiento A y contiene O2 a 800 kPa
y 600 K. Ahora el sistema intercambia calor con los alrededores hasta que el N2 y el CO alcanzan una
temperatura de 300 K.
Usando las tablas A-9 del libro Fundamentals of Thermodynamics de Sonntag-van Wylen, 7ª edición,
determine:
a) el calor total transferido, en kJ
b) la presión final, en kPa

Compartimiento A Compartimiento B
mA = 1 kg de N2 O2
T1A = 400 K T1B = 600 K
P1A = 800 kPa P1B = 800 kPa
V1A = 2V1B

pistón

Solución:
Constantes específicas: RN2 = RA = 0,2968 kJ/kg-K; RO2 = RB = 0,2598 kJ/kg-K
La masa en cada compartimiento es constante.

Balance de materia:

mA R A T1A 1,0 ∗ 0,2968 ∗ 400


V1A = = = 0,1484 m3 ; V1B = 0,2968 m3
P1A 800

P1B V1B 800 ∗ 2 ∗ 0,1484


mB = = = 1,523 kg
R B T1B 0,2598 ∗ 600

Trabajo neto a los alrededores es igual a cero, W = 0

Ecuación de energía:

Q = W + ∆U = ∆UN2 + ∆UCO = mN2 (u2 − u1 )N2 + mCO (u2 − u1 )CO

A partir de las Tablas A-9 se busca la energía interna del nitrógeno y del oxígeno

Para el N2, a 400 K, uN2 = 297,09 kJ/kg; a 300 K, uN2 = 222,63 kJ/kg
Para el O2, a 600 K, uO2 = 404,46 kJ/kg; a 300 K, uO2 = 195,20 kJ/kg
∆UN2 = 1.0 ∗ (222,63 − 297,09) = −74,46 kJ
∆UN2 = 1.523 ∗ (195,20 − 404,46) = −318,75 kJ
Fundamentos de Termodinámica

𝐐 = −𝟕𝟒, 𝟒𝟔 − 𝟑𝟏𝟖, 𝟕𝟓 = −𝟑𝟗𝟑, 𝟐𝟏 𝐤𝐉

Presión final P2A = P2B = P2; temperatura final T2A = T2B = T2 = 300 K;

Volumen total V2 = V1A + V1B = 0,1484*3 = 0,4452 m3

(𝐦𝐀 𝐑 𝐀 + 𝐦𝐁 𝐑 𝐁 )𝐓𝟐 (𝟏, 𝟎 ∗ 𝟎, 𝟐𝟗𝟔𝟖 + 𝟏, 𝟓𝟐𝟑 ∗ 𝟎, 𝟐𝟓𝟗𝟖) ∗ 𝟑𝟎𝟎


𝐏𝟐 = = = 𝟒𝟔𝟔, 𝟔𝟕 𝐤𝐏𝐚
𝐕𝟐 𝟎, 𝟒𝟒𝟓𝟐

mA R A T2 1,0 ∗ 0,2968 ∗ 300


V2A = = = 0,1908 m3 ; V2B = V2 − V2A = 0,2544 m3
P2 466,67
Fundamentos de Termodinámica

5. PRIMERA Y SEGUNDA LEYES DE LA TERMODINÁMICA

El principio de Conservación de la Energía y el principio del Incremento de la Entropía del


Universo están representados mediante las siguientes ecuaciones para un sistema
abierto:

NCS NCE
∂U
[ ] + ∑ [ṁj (h + ec + ep ) ] − ∑ [ṁi (h + ec + ep ) ] − (Q̇neto − Ẇneto ) = 0
∂t vc j i
j=1 i=1

dS Q̇k
Ṡgen = [ ] + ∑[ṁj sj ] − ∑[ṁi si ]ent + ∑ ≥0
dt vc sal Tk
j=1 i=1 k=1

Donde, Q̇ neto = δQ⁄∂t y Ẇneto = δW⁄∂t

Para un sistema cerrado, las ecuaciones anteriores se transforman en:

∂U δQ δW
[ ] −( − ) = 0; ∆U = Q − W (5 − 1)
∂t vc ∂t ∂t

δSgen dS 1 δQk Qk
= [ ] +∑ ≥ 0; Sgen = ∆S + ∑ ≥0 (5 − 2)
dt dt vc Tk dt Tk
k=1 k=1

Esquema de un sistema abierto


Fundamentos de Termodinámica

APLICACIONES DE LA PRIMERA LEY EN SISTEMAS ABIERTOS

PROCESOS DE CARGA:

En los procesos de carga se considera que no hay corrientes másicas de salida y que
hay una corriente másica de entrada; sin embargo, puede haber intercambio de energía
en forma de calor y/o trabajo.
Suposiciones:
1. Flujo uniforme de calor y trabajo
2. Entalpía y entropía de la corriente de entrada constantes

Balance de materia:

m2 − m1 = mi

Balance de energía:

dU δmi δQ δW
[ ] − hi = −
∂t vc ∂t ∂t ∂t

Balance de entropía:

δSgen dS δmi 1 δQk


= [ ] − si +∑ ≥0
dt ∂t vc ∂t Tk ∂t
k=1

Eliminando el diferencial de tiempo en los balances de energía y entropía e integrando,


2
∫ dU − hi ∫ δmi = ∫ δQ − ∫ δW
1

2
1
∫ δSgen = ∫ dS − si ∫ δmi + ∑ ∫ δQk ≥ 0
1 Tk
k=1

(m2 u2 − m1 u1 ) − (m2 − m1 )hi = Q − W (5 − 4)

Qk
Sgen = (m2 s2 − m1 s1 ) − (m2 − m1 )si + ∑ ≥0 (5 − 5)
Tk
k=1

Si el proceso de carga es adiabático, no hay trabajo realizado y el volumen de control


está inicialmente evacuado (m1 = 0), las ecuaciones (5-4) y (5-5) quedan reducidas a:

u2 = hi

Sgen = mi (s2 − si )
Fundamentos de Termodinámica

Ejemplo (Problema 6.111 Sonntag-van Wylen). Un tanque de 25 L inicialmente evacuado está conectado
mediante una válvula a una línea de suministro de aire a 20° C y 800 kPa. La válvula se abre y el aire fluye
al interior del tanque hasta que la presión alcanza 600 kPa. Determine la temperatura, la masa final en el
tanque y la entropía generada. Suponga que el proceso es adiabático. Encuentre una expresión que
relacione la temperatura de la línea de suministro de aire con la temperatura final en el tanque, usando
calores específicos constantes.

Solución:
T TC
p dT
h = ∫ Cp dT; s= ∫ − Rln(P) = sT0 − Rln(P)
T
Raire = 0,287 kPa.m3/kg-K
V1 = V2 = 0,025 m3; P2 = 600 kPa
Ti = 20° C (293 K); Pi = 800 kPa
Consideraciones: m1 = 0; m2 = mi; Q = 0; W = 0
Ecuación de energía: u2 = hi
Ecuación de entropía: Sgen = mi (s2 − si )

De la tabla del aire (Tabla A7-1 del libro Fundamentals of Thermodynamics de Sonntag-van Wylen), por
interpolación a 20° C (293 K), hi = 293,5 kJ/kg, siT° = 6,84565 kJ/kg-K

si = siT° - Raire*ln(Pi) = 6,84565 – 0,287*ln(8) = 6,24885 kJ/kg-K

u2 = hi = 293,5 kJ/kg

Debe haber una temperatura a la cual, u2 sea igual a 293,5 kJ/kg; usando la misma tabla A7-1, por
interpolación,

T2 = 405,51 K (132,51° C)

A T2 = 405,51 K, por interpolación usando la tabla A7-1, s2T° = 7,17290 kJ/kg-K

s2 = s2T° - Raire*ln(P2) = 7,17290 – 0,287*ln(6) = 6,65866 kJ/kg-K

De la ecuación de los gases ideales,

𝐏𝟐 𝐕𝟐 𝟔𝟎𝟎 ∗ 𝟎, 𝟎𝟐𝟓
𝐦𝟐 = 𝐦𝐢 = = = 𝟎, 𝟏𝟐𝟖𝟗 𝐤𝐠
𝐑 𝐚𝐢𝐫𝐞 𝐓𝟐 𝟎, 𝟐𝟖𝟕 ∗ 𝟒𝟎𝟓, 𝟓𝟏

Cambio de entropía: 𝐒𝐠𝐞𝐧 = 𝐦𝐢 (𝐬𝟐 − 𝐬𝐢 ) = 𝟎, 𝟏𝟐𝟖𝟗 ∗ (𝟔, 𝟔𝟓𝟖𝟔𝟔 − 𝟔, 𝟐𝟒𝟖𝟖𝟓) = 𝟎, 𝟎𝟓𝟐𝟖 𝐤𝐉/𝐊

A partir de la ecuación de la energía para este proceso específico,


T2 Ti T2 T2
u2 = hi = ∫ Cv dT = ∫ Cp dT = ∫ Cp dT − RT2 ; ∫ Cp dT = RT2
0 0 0 Ti

𝐂𝐩 𝐓𝐢
Si Cp es constante, 𝐂𝐩 (𝐓𝟐 − 𝐓𝐢 ) = 𝐑𝐓𝟐 ; 𝐓𝟐 =
𝐂𝐩 −𝐑

Ejemplo (Problema 6.108 Sonntag-van Wylen). Un tanque rígido de acero (Cacero = 0,46 kJ/kg-K) de 40
kg y 1 m3 contiene aire a 500 kPa. Tanto el tanque como el aire están a 20° C. El tanque está conectado a
una línea de flujo de aire a 2MPa y 20° C. La válvula se abre permitiendo que el aire fluya hacia el tanque
hasta que la presión alcanza 1,5 MPa y luego se cierra. Suponga que el tanque y el aire están siempre a
la misma temperatura. Si la temperatura final es de 35° C, determine la masa final en el tanque y el calor
transferido.
Tome Cpaire = 1,004 kJ/kg-K (constante)
Fundamentos de Termodinámica

Solución.
Raire = R = 0,287 kJ/kg-K
Vtanque = V1 = V2 = 1 m3
Sistema: aire y tanque
Proceso de carga.
La masa y el volumen del tanque son constantes. mtanque = 40 kg
Trabajo de frontera realizado es igual a cero
Consideraciones: Flujos uniformes. Las propiedades termodinámicas de entrada son constantes
kJ
Condición de entrada: Ti = 20° C (293 K), hi = Cpi Ti = 1,004 ∗ 293 = 294,17
kg−K

P1 V1 500∗1,0
Condición inicial (1): T1 = 20° C (293 K), P1 = 500 kPa, m1 = = = 5,946 kg
RT1 0,287∗293

𝐏𝟐 𝐕𝟐 𝟏𝟓𝟎𝟎∗𝟏,𝟎
Condición final (2): T2 = 35° C (308 K), P2 = 1500 kPa, 𝐦𝟐 = = = 𝟏𝟔, 𝟗𝟕 𝐤𝐠
𝐑𝐓𝟐 𝟎,𝟐𝟖𝟕∗𝟑𝟎𝟖
Balance de materia del aire:
m2 − m1 = mi = 11,023 kg
Balance de energía del sistema:
(m2 u2 − m1 u1 ) + mtanque Ctanque (T2 − T1 ) − mi hi = Q

Pero u = CvT, donde Cv = Cp – R = 1,004 – 0,287 = 0,717 kJ/kg-K

(m2 Cv T2 − m1 Cv T1 ) + mtanque Ctanque (T2 − T1 ) − mi Cpi Ti = Q

𝐐 = 𝟎, 𝟕𝟏𝟕 ∗ (𝟏𝟔, 𝟗𝟕 ∗ 𝟑𝟎𝟖 − 𝟓, 𝟗𝟒𝟔 ∗ 𝟐𝟗𝟑) + 𝟒𝟎 ∗ 𝟎, 𝟒𝟔 ∗ (𝟑𝟎𝟖 − 𝟐𝟗𝟑) − 𝟏𝟏, 𝟎𝟐𝟑 ∗ 𝟐𝟗𝟒, 𝟏𝟕 = −𝟒𝟔𝟖, 𝟏𝟗 𝐤𝐉

Ejemplo. Si el proceso del ejemplo anterior fuese adiabático, con las mismas condiciones de entrada e
iniciales, determine la presión y la masa final en el tanque.

Solución.
Raire = R = 0,287 kJ/kg-K
Vtanque = V1 = V2 = 1 m3
Sistema: aire y tanque
Proceso de carga.
La masa y el volumen del tanque son constantes. m tanque = 40 kg
Trabajo de frontera realizado es igual a cero, W = 0
Proceso adiabático, Q = 0
Consideraciones: Flujos uniformes. Las propiedades termodinámicas de entrada son constantes
kJ
Condición de entrada: Ti = 20° C (293 K), hi = Cpi Ti = 1,004 ∗ 293 = 294,17
kg−K

P1 V1 500∗1,0
Condición inicial (1): T1 = 20° C (293 K), P1 = 500 kPa, m1 = = = 5,946 kg
RT1 0,287∗293

m2 RT2
Condición final (2): T2 = 35° C (308 K), P2 = ¿?, P2 =
V2
Balance de materia del aire:
m2 − m1 = mi = 11,023 kg
Balance de energía del sistema:
(m2 u2 − m1 u1 ) + mtanque Ctanque (T2 − T1 ) − mi hi = 0

Pero u = CvT, donde Cv = Cp – R = 1,004 – 0,287 = 0,717 kJ/kg-K

(m2 Cv T2 − m1 Cv T1 ) + mtanque Ctanque (T2 − T1 ) − (m2 − m1 )Cpi Ti = 0

mtanque Ctanque (T2 − T1 ) + m1 Cpi Ti − m1 Cv T1 = m2 Cpi Ti − m2 Cv T2


Fundamentos de Termodinámica

mtanque Ctanque (T2 − T1 ) + m1 (Cpi Ti − Cv T1 ) = m2 (Cpi Ti − Cv T2 )

[𝐦𝐭𝐚𝐧𝐪𝐮𝐞 𝐂𝐭𝐚𝐧𝐪𝐮𝐞 (𝐓𝟐 − 𝐓𝟏 ) + 𝐦𝟏 (𝐂𝐩𝐢 𝐓𝐢 − 𝐂𝐯 𝐓𝟏 )]


𝐦𝟐 = = 𝟏𝟎, 𝟓𝟖 𝐤𝐠
(𝐂𝐩𝐢 𝐓𝐢 − 𝐂𝐯 𝐓𝟐 )
𝐦𝟐 𝐑𝐓𝟐 𝟏𝟎, 𝟓𝟖 ∗ 𝟎, 𝟐𝟖𝟕 ∗ 𝟑𝟎𝟖
𝐏𝟐 = = = 𝟗𝟑𝟓, 𝟐𝟑 𝐤𝐏𝐚
𝐕𝟐 𝟏, 𝟎

PROCESOS DE DESCARGA:

En los procesos de descarga se considera que no hay corrientes másicas de entrada y


que hay una corriente másica de salida; sin embargo, puede haber intercambio de
energía en forma de calor y/o trabajo.
Suposiciones:
1. Flujo uniforme de calor y trabajo
2. Entalpía y entropía de la corriente de salida a las condiciones promedio entre los
estados inicial y final
(h1 + h2 ) (s1 + s2 )
hj promedio = ; sj promedio =
2 2

Balance de materia:
m1 = m2 + mj

Balance de energía:
d(mu) δmj δQ δW
[ ] + hj promedio = −
∂t vc ∂t ∂t ∂t

Balance de entropía:

δSgen d(ms) δmj 1 δQk


=[ ] + sj promedio + ∑ ≥0
∂t ∂t vc ∂t Tk ∂t
k=1

Eliminando el diferencial de tiempo en los balances de energía y entropía e integrando,

(m2 u2 − m1 u1 ) + (m1 − m2 )hj promedio = Q − W (5 − 6)

Qk
Sgen = (m2 s2 − m1 s1 ) + (m1 − m2 )sj promedio + ∑ ≥0 (5 − 7)
Tk
k=1

Ejemplo. Una planta de vapor está basada en una turbina y una caldera. El tanque de la caldera tiene un
volumen de 100 L e inicialmente contiene líquido saturado con una cantidad despreciable de vapor saturado
a 100 kPa. Ahora de suministra calor al quemador de la caldera desde una fuente térmica que está a 177°
C. El regulador de presión mantiene la presión constante y no se abre antes de que se alcance una presión
de 700 kPa en la caldera. El vapor saturado a 700 kPa entra a la turbina que lo descarga como vapor
saturado a 100 kPa. El quemador se apaga cuando no haya líquido en el tanque de la caldera. Determine
Fundamentos de Termodinámica

el trabajo producido por la turbina, el calor transferido al quemador de la caldera y la entropía generada en
todo el proceso.

Solución:
Proceso de descarga en la caldera no estacionario
Producción de trabajo por una turbina, W T en estado estacionario
Se requiere un calor en el quemador de la caldera, Q C desde una fuente que está a TF = 177° C (450 K)

vapor saturado a 700 kPa


,

Regulador de
presión

Liquido saturado

Turbina adiabática vapor saturado a 100 kPa


Calor

Estado 1:
Líquido saturado a 100 kPa; Volumen constante, V = 100 L = 0,1 m 3
De las tablas de vapor, a 100 kPa, Tsat = 99,606° C, vf = v1 = 1,043*10-3 m3/kg, uf = u1 = 417,4 kJ/kg
sf = s1 = 1,30276 kJ/kg-K
m1 = V/v1 = 0,1/1,043*10-3 = 95,8773 kg

Estado 2:
Vapor saturado a 700 kPa
De las tablas de vapor, a 700 kPa, T sat = 164,946° C, vg = v2 = 0,27278 m3/kg, ug = u2 = 2571,81 kJ/kg
sg = s2 = 6,70706 kJ/kg-K
m2 = V/v2 = 0,1/0,27278 = 0,3666 kg

Salida de la caldera:
Vapor saturado a 700 kPa, hg = hj = 2762,75 kJ/kg, sg = sj = 6,70706 kJ/kg-K

Salida de la turbina:
Vapor saturado a 100 kPa, hg = 2674,95 kJ/kg, sg = 7,35885 kJ/kg-K

Balance en la caldera:
mj = m1 – m2 = 95,5107 kg

El regulador de presión controla que la entalpía y la entropía de salida sea la del vapor saturado a 700 kPa

Q = (m2 u2 − m1 u1 ) + mj hj

𝐐 = 𝟎, 𝟑𝟔𝟔𝟔 ∗ 𝟐𝟓𝟕𝟏, 𝟖𝟏 − 𝟗𝟓, 𝟖𝟕𝟕𝟑 ∗ 𝟒𝟏𝟕, 𝟒 + 𝟗𝟓, 𝟓𝟏𝟎𝟕 ∗ 𝟐𝟕𝟔𝟐, 𝟕𝟓 = 𝟐𝟐𝟒𝟕𝟗𝟓, 𝟕𝟕 𝐤𝐉 = −𝐐𝐅

QF
(Sgen ) = (m2 s2 − m1 s1 ) + mj sj +
caldera TF
224795,77 kJ
(Sgen ) = (0,366 ∗ 6,70706 − 95,8773 ∗ 1,30276) + 95,5107 ∗ 6,70706 − = 18,60
caldera 450 K

Balance en la la válvula y la turbina:


mi = mj = 95,5107 kg
Fundamentos de Termodinámica

Wturbina = mj (hi − hj )

𝐖𝐭𝐮𝐫𝐛𝐢𝐧𝐚 = 𝟗𝟓, 𝟓𝟏𝟎𝟕 ∗ (𝟐𝟕𝟔𝟐, 𝟕𝟓 − 𝟐𝟔𝟕𝟒, 𝟗𝟓) = 𝟖𝟑𝟖𝟓, 𝟖𝟒 𝐤𝐉

kJ
(Sgen ) = mj (si − sj ) = 95,5107 ∗ (7,35885 − 6,70706) = 62,253
turbina+válvula K

Balance de entropía en el proceso (caldera, válvula y turbina):

QF
(Sgen ) = (m2 s2 − m1 s1 ) + mj sj + = (Sgen ) + (Sgen )
proceso TF caldera t+v

(𝐒𝐠𝐞𝐧 ) = (𝟎, 𝟑𝟔𝟔 ∗ 𝟔, 𝟕𝟎𝟕𝟎𝟔 − 𝟗𝟓, 𝟖𝟕𝟕𝟑 ∗ 𝟏, 𝟑𝟎𝟐𝟕𝟔) + 𝟗𝟓, 𝟓𝟏𝟎𝟕 ∗ 𝟕, 𝟑𝟓𝟖𝟖𝟓
𝐩𝐫𝐨𝐜𝐞𝐬𝐨
𝟐𝟐𝟒𝟕𝟗𝟓, 𝟕𝟕 𝐤𝐉
− = 𝟖𝟎, 𝟖𝟓𝟑
𝟒𝟓𝟎 𝐊

APLICACIONES DE LA PRIMERA Y SEGUNDA LEYES EN SISTEMAS ABIERTOS EN


ESTADO ESTACIONARIO

Las ecuaciones generales de la primera y segunda leyes para sistemas abiertos en


estado estacionario son las siguientes:

NCS NCE

∑ [ṁj (h + ec + ep ) ] − ∑ [ṁi (h + ec + ep ) ] − (Q̇neto − Ẇneto ) = 0 (5 − 8)


j i
j=1 i=1

Q̇k
Ṡgen = ∑[ṁj sj ] − ∑[ṁi si ]ent + ∑ ≥0 (5 − 9)
sal Tk
j=1 i=1 k=1

Los equipos o dispositivos donde se pueden aplicar estas ecuaciones son: toberas,
difusores, intercambiadores de calor, mezcladores de corrientes o equipos de mezclado,
separador de corrientes o equipos de separación (columnas de destilación, columnas de
absorción, etc), turbinas, compresores y bombas, unidades de procesos, ciclos de
potencia, ciclos de refrigeración, etc.

Cada equipo o dispositivo tiene un objetivo y contribuye a alcanzar el objetivo de un


proceso o unidad de proceso. Para alcanzar su objetivo necesita un suministro de energía
o una interconversión de su energía.

Toberas
El objetivo de una tobera es convertir la energía térmica y mecánica de un fluido en
energía cinética. No hay potencia y el flujo de calor es despreciable. Hay una entrada (1)
y una salida (2). Se desprecia el cambio de la energía potencial.

Balance de materia:
Fundamentos de Termodinámica

A1 u̅1 A2 u̅2 A1 v1 u̅2


ṁ1 = = ṁ2 = ; = ( )( )
v1 v2 A2 v2 u̅1
Balance de energía:

2
u̅2 = √[2(h1 − h2 ) + u̅12 ]

Balance de entropía:
Ṡgen = ṁ1 (s2 − s1 ) ≥ 0

En el caso que la velocidad de entrada sea despreciable, la ecuación se reduce a:

2
u̅2 = √2(h1 − h2 )

Ejemplo. Vapor sobrecalentado de amoniaco entra a una tobera aislada a 20° C y 500 kPa a velocidad
despreciable y un flujo estable de 0,01 kg/s. El amoniaco sale a 200 kPa y una velocidad de 450 m/s.
Determine: a) la temperatura (o calidad si está saturado); b) la entropía generada

0,01 kg/s de amoniaco


20° C 0,01 kg/s de amoniaco
500 kPa 200 kPa
1 450 m/s
2

Solución:
h en kJ/kg; s en kJ/kg-K
Se usan las tablas del amoniaco (Tabla A.7.7 de Sonntag-van Wylen)
A 20° C y 500 kPa, h1 = 1488,3; s1 = 5,4244
Balance de materia: ṁ1 = ṁ2 = 0,01 kg/s
Balance de energía:
u̅22 4502
h2 = h1 − = 1488,3 − = 1387,05 kJ/kg
2 2000

A 200 kPa, hg = 1419,6 > h2. El estado de la salida es mezcla líquido-vapor.

T2 = temperatura de saturación a 200 kPa = -18,86° C)

hf = 93,86 kJ/kg; hfg = 1325,74; sf = 0,3844; sfg = 5,2135


La calidad del vapor a la salida: x2 = (h2 – hf)/hfg = (1387,05 – 93,86)/1325,74 = 0,9754

Calidad del vapor, x2 = 97,54%

s2 = sf + x2*sfg = 0,3844 + 0,9754*5,2135 = 5,4696

Balance de entropía:
𝐒̇𝐠𝐞𝐧 = 𝐦̇𝟏 (𝐬𝟐 − 𝐬𝟏 ) = 𝟎, 𝟎𝟏 ∗ (𝟓, 𝟒𝟔𝟗𝟕 − 𝟓, 𝟒𝟐𝟒𝟒) = 𝟒, 𝟓𝟐 ∗ 𝟏𝟎−𝟒 𝐤𝐖/𝐊 ≥ 𝟎

Ejemplo. Vapor saturado de refrigerante R134a entra a una tobera aislada a 800 kPa a velocidad de 10
m/s y un flujo estable de 0,5 kg/s. El refrigerante sale a 200 kPa. Si la relación entre el área de entrada y
de salida es de 5, determine: a) la velocidad de salida; b) la temperatura de salida (o calidad si está
saturado); c) la entropía generada
Fundamentos de Termodinámica

0,5 kg/s de R134a


Vapor saturado 0,5 kg/s de R134a
800 kPa 200 kPa
10 m/s A1/A2 = 5
2
1

Solución:
v en m3/kg; h en kJ/kg; s en kJ/kg-K
Se usan las tablas del R134a (Tablas A.11-13 de Cengel)
A 800 kPa, vapor saturado, v1 = 0,02562; h1 = 267,29; s1 = 0,91835
kg A u
̅ A u
̅
Balance de materia: ṁ1 = ṁ2 = 0,01 = 1 1 = 2 2
s v1 v2
A 1 v2 v2
u̅2 = u̅1 ( ) ( ) = 10 ∗ 5 ∗ ( ) = 1951,6 ∗ v2
A2 v1 0,02562
Balance de energía:
u̅12 u̅22 102 u̅22 u̅22 (1951,6)2 ∗ v22
h2 = h1 + − = 267,29 + − = 267,34 − = 267,34 −
2 2 2000 2000 2000 2000

h2 = 267,34 − 1904,37v22

A 200 kPa, suponer calidad o temperatura, estimar volumen específico y entalpía. Verificar la ecuación
anterior.

La temperatura de salida es aproximadamente -6,35° C

Las propiedades son: v2 = 0,1017; h2 = 247,64; s2 = 0,9496

La velocidad de salida es: 𝐮


̅ 𝟐 = 𝟏𝟗𝟓𝟏, 𝟔 ∗ 𝟎, 𝟏𝟎𝟏𝟕 = 𝟏𝟗𝟖, 𝟒𝟕 𝐦/𝐬

Balance de entropía:
𝐒̇𝐠𝐞𝐧 = 𝐦̇𝟏 (𝐬𝟐 − 𝐬𝟏 ) = 𝟎, 𝟓 ∗ (𝟎, 𝟗𝟒𝟗𝟔 − 𝟎, 𝟗𝟏𝟖𝟑𝟓) = 𝟎, 𝟎𝟏𝟓𝟔𝟐𝟓 𝐤𝐖/𝐊 ≥ 𝟎

Difusores
El objetivo de un difusor es aumentar la energía mecánica o aumento de la presión de un
fluido a expensas de su energía cinética. Igual que en la tobera, no hay potencia, ni flujo
de calor, el cambio de energía potencial es despreciable y hay una entrada y una salida.

El balance de materia es igual al de la tobera.

Balance de energía:
u̅12 u̅22
h2 = h1 + [ − ]
2 2
Balance de entropía:
Ṡgen = ṁ1 (s2 − s1 ) ≥ 0

Si la velocidad de salida es despreciable,

u̅12
h2 = h1 +
2
Fundamentos de Termodinámica

Ejemplo. Un motor de propulsión a chorro consta de un difusor, seguido de un compresor, una cámara de
combustión, una turbina y una tobera.
El difusor recibe aire a 900 km/h, -5° C y 50 kPa que sale a 80 m/s antes de entrar al compresor. Si el área
de flujo de la salida se aumenta en un 150% del área de entrada, determine: a) la temperatura y presión a
la entrada del compresor; b) la entropía generada por kg/s de aire.
Considerar Cp del aire constante e igual a 1,004 kJ/kg-K

Solución:
Base de cálculo: m1 = 1,0 kg/s
Cp = 1,004 kJ/kg-K; R = 0,287 kJ/kg-K
T1 = -5° C = 268 K; u1 = 900 km/h = 250 m/s; u2 = 80 m/s
v1 = RT1/P1 = 0,287*268/50 = 1,5383 m 3/kg; A2/A1 = 2,5
A2 u̅2 80 m3
v2 = v1 ( ) ( ) = 1,5383 ∗ 2,5 ∗ = 1,2306
A1 u̅1 250 kg
Balance de energía:
u̅12 u̅22 u̅12 u̅22
h2 = h1 + [ − ] = Cp T2 = Cp T1 + [ − ]
2 2 2 2
𝟐 𝟐
𝟏 𝐮 ̅𝟏 𝐮 ̅𝟐 𝟏 𝟐𝟓𝟎𝟐 𝟖𝟎𝟐
𝐓𝟐 = 𝐓𝟏 + [ − ] = 𝟐𝟔𝟖 + [ − ] = 𝟐𝟗𝟓, 𝟗𝟒 𝐊
𝐂𝐩 𝟐 𝟐 𝟏, 𝟎𝟎𝟒 𝟐𝟎𝟎𝟎 𝟐𝟎𝟎𝟎
𝐑𝐓𝟐 𝟎, 𝟐𝟖𝟕 ∗ 𝟐𝟗𝟓, 𝟗𝟒
𝐏𝟐 = = = 𝟔𝟗, 𝟎𝟏𝟗 𝐤𝐏𝐚
𝐯𝟐 𝟏, 𝟐𝟑𝟎𝟔
Balance de entropía:
Ṡgen = ṁ1 (s2 − s1 ) ≥ 0
𝐓𝟐 𝐏𝟐 𝟐𝟗𝟓, 𝟗𝟒 𝟔𝟗, 𝟎𝟏𝟗
𝐒̇𝐠𝐞𝐧 = 𝟏, 𝟎 ∗ [𝐂𝐩 𝐥𝐧 ( ) − 𝐑𝐥𝐧 ( )] = 𝟏, 𝟎𝟎𝟒 ∗ 𝐥𝐧 ( ) − 𝟎, 𝟐𝟖𝟕 ∗ 𝐥𝐧 ( ) = 𝟎, 𝟎𝟎𝟕𝟎𝟓 𝐤𝐖/𝐊
𝐓𝟏 𝐏𝟏 𝟐𝟔𝟖 𝟓𝟎

Intercambiadores de calor
Un intercambiador de calor permite que dos fluidos, uno caliente y el otro frio,
intercambien su energía térmica; el fluido caliente cede energía mientras que el fluido frio,
gana energía. Si el intercambiador es adiabático, la totalidad de la energía cedida por el
fluido caliente, la gana el fluido frio. En caso contrario, una fracción de la energía cedida
se elimina hacia los alrededores si los niveles de temperatura están por encima de la
temperatura de los alrededores.

Un intercambiador de calor puede clasificarse según algún tipo de criterio.

Desde el punto de vista del fluido:


- Calentador, si el fluido importante se calienta. El fluido caliente puede ser vapor
saturado, gases de combustión o gases de escape
- Enfriador, si el fluido importante se enfría. El fluido frio puede ser agua de
enfriamiento o aire
- Condensador, si el fluido importante se condensa. El fluido frio puede ser agua de
enfriamiento
- Vaporizador o evaporador, si el fluido importante se vaporiza o evapora. El fluido
caliente puede ser gases de combustión

Desde el punto de vista de la disposición de los fluidos:


- Flujo en paralelo. Las corrientes de los fluidos van en el mismo sentido
Fundamentos de Termodinámica

- Flujo en contracorriente. Las corrientes de los fluidos van en sentido contrario


- Flujo cruzado
- Tubo y coraza
- Doble tubo

El análisis termodinámico que se realiza tiene en cuenta las condiciones de entrada y


salida de las corrientes y el resultado final estará de acuerdo al tipo de intercambiador.

Considere un intercambiador de calor en donde se ponen en contacto dos fluidos. El


fluido caliente entra en el estado 1 y sale en el estado 2, mientras que el fluido frio entra
en el estado 3 y sale en el estado 4. Los cambios de energía cinética y potencial son
despreciables.

Balance de materia:

ṁ2 + ṁ4 = ṁ1 + ṁ3 ; pero, ṁ2 = ṁ1 y ṁ4 = ṁ3

Balance de energía:
ṁ1 (h2 − h1 ) + ṁ3 (h4 − h3 ) = Q̇

Balance de entropía:
Q̇ 0
Ṡgen = ṁ1 (s2 − s1 ) + ṁ3 (s4 − s3 ) + ≥0
T0

Si el intercambiador de calor es adiabático,

ṁ1 (h1 − h2 ) = ṁ3 (h4 − h3 )

Si el intercambiador no es adiabático y se elimina calor hacia los alrededores, Q 0,

ṁ1 (h1 − h2 ) = ṁ3 (h4 − h3 ) + Q̇0

Para un condensador, donde el fluido caliente cambia de fase a temperatura constante,


las ecuaciones de energía y entropía son las siguientes:

ṁ1 hfg = ṁ3 (h4 − h3 )


Fundamentos de Termodinámica

Ṡgen = ṁ1 (−sfg ) + ṁ3 (s4 − s3 ) ≥ 0

Donde hfg y sfg son la entalpía y entropía latente de vaporización o de condensación.

Para un evaporador o generador de vapor saturado, donde el fluido frio se vaporiza,

ṁ1 (h1 − h2 ) = ṁ3 hfg

Ṡgen = ṁ1 (s2 − s1 ) + ṁ3 sfg ≥ 0

Téngase en cuenta que Q0 es positivo con respecto a los alrededores cuando los niveles
de temperatura son mayores que T0. En los evaporadores de un ciclo de refrigeración,
Q0 es negativo.

Ejemplo. Un intercambiador de calor recibe 50 kg/s de CO2 a 250 kPa y 620 K donde intercambia energía
térmica saliendo a 230 kPa y 460 K. El fluido frio es aire que entra a 37° C y 110 kPa y sale a 287° C y 105
kPa. Usando las tablas de CO2 y del aire, determine: a) el flujo másico de aire, en kg/s; b) la entropía
generada en el intercambiador. Determine a) y b) si el 5% de la energía suministrada por el fluido caliente
se pierde hacia los alrededores que están a 35° C.
Solución:
Fuente: Tablas del CO2 y aire de Cengel
Para el CO2, h en kJ/kmol, s° en kJ/kmol-K; para el aire, h en kJ/kg, s° en kJ/kg-K
Estado 1: entrada del CO2; Estado 2: salida del CO2; Estado 3: entrada del aire; Estado 4: salida del aire
Fluido caliente: CO2, m1 = 50 kg/s; T1 = 620 K; P1 = 250 kPa; T2 = 460 K, P2 = 230 kPa
Fluido frio: aire, m 2 = ¿?; T3 = 310 K; P3 = 110 kPa; T4 = 560 K, P4 = 105 kPa
Alrededores: T0 = 308 K
De las tablas A-20 del CO2:
A 620 K, h1 = 23231 kJ/kmol; s°1 = 244,758 kJ/kmol-K; s1 = 244,758 – 8,314*ln(2,5) = 237,4 kJ/kmol-K
A 460 K, h2 = 15916 kJ/kmol; s°2 = 231,144 kJ/kmol-K; s2 = 231,144 – 8,314*ln(2,3) = 224,22 kJ/kmol-K
De las tablas A-17 del aire:
A 310 K, h3 = 310,24 kJ/kg; s°3 = 1,73498 kJ/kg-K; s3 = 1,73498 – 0,287¨ln(1,10) = 1,7076 kJ/kg-K
A 560 K, h4 = 565,17 kJ/kg; s°4 = 2,33685 kJ/kg-K; s4 = 2,33685 – 0,287*ln(1,05) = 2,3228 kJ/kg-K
50
Energía suministrada: Ėsum = ṁ1 (h1 − h2 ) = ∗ (23231 − 15916) = 8312,5 kW
44
Balance de energía (intercambiador adiabático)
(𝐡𝟏 − 𝐡𝟐 ) 𝟓𝟎 (𝟐𝟑𝟐𝟑𝟏 − 𝟏𝟓𝟗𝟏𝟔) 𝐤𝐠
𝐦̇𝟏 (𝐡𝟏 − 𝐡𝟐 ) = 𝐦̇𝟑 (𝐡𝟒 − 𝐡𝟑 ); 𝐦̇𝟑 = 𝐦̇𝟏 = ∗ = 𝟑𝟐, 𝟔𝟎𝟕
(𝐡𝟒 − 𝐡𝟑 ) 𝟒𝟒 (𝟓𝟔𝟓, 𝟏𝟕 − 𝟑𝟏𝟎, 𝟐𝟒) 𝐬
Balance de entropía:
50
Ṡgen = ṁ1 (s2 − s1 ) + ṁ3 (s4 − s3 ) = ( ) ∗ (224,22 − 237,44) + 32,607 ∗ (2,3228 − 1,7076)
44
𝐤𝐖
𝐒̇𝐠𝐞𝐧 = 𝟓, 𝟎𝟑𝟕
𝐊
Balance de energía (intercambiador no adiabático)
Q̇ 0 = 0,05 ∗ Ėsum = 0,05 ∗ 8312,5 = 415,625 kW
ṁ1 (h1 − h2 ) = ṁ3 (h4 − h3 ) + Q̇ 0
(𝐡𝟏 − 𝐡𝟐 ) 𝐤𝐠
𝐦̇𝟑 = 𝟎, 𝟗𝟓 ∗ 𝐦̇𝟏 = 𝟑𝟎, 𝟗𝟕𝟕
(𝐡𝟒 − 𝐡𝟑 ) 𝐬
Balance de entropía:
Q̇ 0
Ṡgen = ṁ1 (s2 − s1 ) + ṁ3 (s4 − s3 ) +
T0
𝟓𝟎 𝟒𝟏𝟓, 𝟔𝟐𝟓 𝐤𝐖
𝐒̇𝐠𝐞𝐧 = ( ) ∗ (𝟐𝟐𝟒, 𝟐𝟐 − 𝟐𝟑𝟕, 𝟒𝟒) + 𝟑𝟎, 𝟗𝟕𝟕 ∗ (𝟐, 𝟑𝟐𝟐𝟖 − 𝟏, 𝟕𝟎𝟕𝟔) + = 𝟓, 𝟑𝟖𝟒 ≥𝟎
𝟒𝟒 𝟑𝟎𝟖 𝐊
Fundamentos de Termodinámica

Turbina de vapor o de gas


Una turbina de vapor expande vapor saturado o sobrecalentado desde una presión y
temperatura altas hasta una presión y temperatura bajas, produciendo una potencia de
eje que se convierte posteriormente en potencia eléctrica. La condición de descarga
puede ser vapor húmedo. Una turbina de gas maneja gases a temperaturas más altas
pero presiones más bajas que las turbinas de vapor, pero tienen la misma finalidad:
producir potencia de eje.

Esquemas de una turbina de vapor y de una turbogas (compresor-cámara de combustión y turbina de gas)

Las turbinas pueden clasificarse en turbinas de extracción y de condensación. Las


turbinas de extracción son aquellas que tienen dos o más corrientes de salida a presiones
intermedias y una corriente de salida de baja presión. Las turbinas de condensación son
aquellas en la que la corriente de salida se alimenta a un condensador para el uso
posterior del líquido condensado.

En la figura de arriba la turbina de vapor opera entre los estados 1 y 2. Las ecuaciones
de energía y entropía se convierten en:

Ẇturbina = ṁ1 [(h1 + ec1 ) − (h2 + ec2 )] − Q̇0

Q̇ 0
Ṡgen = ṁ1 (s2 − s1 ) + ≥0
T0
Fundamentos de Termodinámica

Una turbina que maneja altos niveles de presión (2 a 7 MPa) puede tener una etapa de
alta presión, una etapa de presión intermedia y una etapa de baja presión. Si el vapor a
presión intermedia se extrae y se recalienta en una caldera y retorna a alta temperatura,
la potencia total de la turbina se incrementa.

La figura de arriba corresponde a un ciclo de potencia con recalentamiento en donde hay


dos turbinas: una de alta presión (TAP) y una de baja presión (TBP) unidas por el mismo
eje. En este caso, la turbina de alta presión opera entre los estados 2 y 3.

Las ecuaciones de energía y entropía son las siguientes:

Ẇturbina = ṁ1 {[(h2 + ec2 ) − (h3 + ec3 )] + [(h4 + ec4 ) − (h5 + ec5 )]} − Q̇ 0

Ṡgen = ṁ1 [(s3 − s2 ) + (s5 − s4 )] + Q̇0 ⁄T0 ≥ 0

Despreciando los cambios en la energía cinética y potencial,

Ẇturbina = ṁ1 [(h2 − h3 ) + (h4 − h5 )] − Q̇0

Si una fracción (x) del flujo másico del vapor a presión intermedia se extrae y se envía a
otro proceso y el resto se retorna a la turbina de baja presión, la potencia total de la
turbina es la suma de las potencias producidas en cada etapa. Las ecuaciones de energía
y entropía de este tipo de turbina son las siguientes:
Vapor a alta presión y alta
temperatura con flujo m1
1 3

TAP TBP
2
(1-x)m1
xm1
Fundamentos de Termodinámica

Ẇturbina = ṁ1 [(h1 − h2 ) + (1 − x)(h2 − h3 )] − Q̇ 0

Ṡgen = ṁ1 [(s2 − s1 ) + (1 − x)(s3 − s2 )] + Q̇0 ⁄T0 ≥ 0

Ejemplo. Una turbina de vapor recibe vapor de dos calderas. Una de las corrientes tiene un flujo de 5 kg/s
a 3 MPa y 700° C y la otra, un flujo de 15 kg/s a 800 kPa y 500° C. La condición de salida es 10 kPa y
calidad de 96%. Determine la potencia de la turbina adiabática y la entropía generada.

5 kg/s 15 kg/s
1 2
3 MPa, 700° C 800 kPa, 500° C

10 kPa, x3 = 96%
3

Solución:
Las propiedades se estiman usando la aplicación ChemicaLogic ™ (h en kJ/kg, s en kJ/kg-K)
Turbina adiabática, Q = 0
Entrada 1: P1 = 3 MPa, T1 = 700° C; h1 = 3912,24; s1 = 7,75896; m1 = 5 kg/s
Entrada 2: P2 = 800 kPa, T2 = 500° C; h2 = 3481,25; s2 = 7,86922; m2 = 15 kg/s
Salida 3: P3 = 10 kPa, x3 = 96%; h3 = 2488,18; s3 = 7,84884
Balance de materia:
ṁ1 + ṁ2 = ṁ3 = 15 kg/s
Balance de energía:
Ẇturbina = ṁ1 (h1 − h3 ) + ṁ2 (h2 − h3 )
𝐖̇𝐭𝐮𝐫𝐛𝐢𝐧𝐚 = 𝟓, 𝟎 ∗ (𝟑𝟗𝟏𝟐, 𝟐𝟒 − 𝟐𝟒𝟖𝟖, 𝟏𝟖) + 𝟏𝟓 ∗ (𝟑𝟒𝟖𝟏, 𝟐𝟓 − 𝟐𝟒𝟖𝟖, 𝟏𝟖) = 𝟐𝟐𝟎𝟏𝟔, 𝟑𝟓 𝐤𝐖
Balance de entropía:
Ṡgen = ṁ1 (s3 − s1 ) + ṁ2 (s3 − s2 )
𝐒̇𝐠𝐞𝐧 = 𝟓, 𝟎 ∗ (𝟕, 𝟖𝟒𝟖𝟖𝟒 − 𝟕, 𝟕𝟓𝟖𝟗𝟔) + 𝟏𝟓 ∗ (𝟕, 𝟖𝟒𝟖𝟖𝟒 − 𝟕, 𝟖𝟔𝟗𝟐𝟐) = 𝟎, 𝟏𝟒𝟑𝟕 𝐤𝐖/𝐊

Compresores y bombas.
El objetivo de un compresor o de una bomba es incrementar la presión del fluido desde
una presión baja a expensas de la potencia suministrada. Los compresores pueden ser
compresores de vapor o compresores de gas, mientras que las bombas operan con
líquidos.
Si el compresor hace parte de una turbogas, el fluido de entrada es aire atmosférico. Si
hace parte en un ciclo de refrigeración, el fluido es vapor saturado o ligeramente
sobrecalentado de un refrigerante.

En los compresores los cambios de energía cinética y potencial son despreciables. En


cambio en las bombas el cambio de energía potencial si es importante.

Ecuaciones de energía y entropía de un compresor no adiabático:

|Ẇcompresor | = ṁ1 (h2 − h1 ) + Q̇0


Fundamentos de Termodinámica

Ṡgen = ṁ1 (s2 − s1 ) + Q̇0 ⁄T0 ≥ 0

Considerando líquido incompresible, la relación entre la entalpía y la presión en las


bombas, además, es la siguiente:

h2 = h1 + (P2 − P1 )v1

Ecuación de energía de una bomba que opera desde una altura Z 1 hasta Z2:

|ẆB | = ṁ1 (h2 + gZ2 − h1 − gZ1 ) = ṁ1 [(P2 − P1 )v1 + g(Z2 − Z1 )]

Donde v1 = volumen específico a la entrada de la bomba

Si el cambio de energía potencial es despreciable, la ecuación para el trabajo de la bomba


específico es la siguiente:

|𝑤B | = (h2 − h1 ) = (P2 − P1 )v1

Ejemplo. Un compresor de vapor toma refrigerante R-134a desde 150 kPa y -10° C hasta 1200 kPa y 50°
C. Si las pérdidas de calor a los alrededores se estiman en 40 kW y la potencia de entrada es de 150 kW,
determine: a) el flujo másico de refrigerante; b) la entropía generada si los alrededores están a 35° C.
Solución:
Fluido: Refrigerante R-134a
Potencia de entrada: Ẇc = 150 kW
Flujo de calor a los alrededores: Q̇ 0 = 40 kW
Entrada: P1 = 150 kPa; T1 = -10° C; h1 = 393,84 kJ/kg; s1 = 1,7606 kJ/kg-K
Salida: P2 = 1200 kPa; T2 = 50° C; h2 = 426,84 kJ/kg; s2 = 1,7237 kJ/kg-K
Temperatura de los alrededores, T 0 = 308 K
Balance de energía:

|Ẇcompresor | = ṁ1 (h2 − h1 ) + Q̇ 0 ;


(|𝐖̇𝐜𝐨𝐦𝐩𝐫𝐞𝐬𝐨𝐫 | − 𝐐̇𝟎 ) (𝟏𝟓𝟎 − 𝟒𝟎) 𝟏𝟎
𝐦̇𝟏 = = = 𝐤𝐠/𝐬
(𝐡𝟐 − 𝐡𝟏 ) (𝟒𝟐𝟔, 𝟖𝟒 − 𝟑𝟗𝟑, 𝟖𝟒) 𝟑

Balance de entropía:
𝟏𝟎 𝟒𝟎 𝐤𝐖
𝐒̇𝐠𝐞𝐧 = 𝐦̇𝟏 (𝐬𝟐 − 𝐬𝟏 ) + 𝐐̇𝟎 ⁄𝐓𝟎 = ∗ (𝟏, 𝟕𝟐𝟑𝟕 − 𝟏, 𝟕𝟔𝟎𝟔) + = 𝟔, 𝟖𝟕 ∗ 𝟏𝟎−𝟑 ≥𝟎
𝟑 𝟑𝟎𝟖 𝐊

Mezcladores o calentadores abiertos.


Los mezcladores son aquellos equipos que reciben dos o más corrientes másicas de una
sustancia pura o mezcla multicomponente con sus respectivas condiciones de presión,
temperatura y flujo, donde se ponen en contacto directamente e intercambian su energía
térmica, saliendo una sola corriente con una condición termodinámica que depende del
intercambio de calor con los alrededores.

Las ecuaciones de energía y entropía para un mezclador de NK corrientes de entrada


son las siguientes:
Fundamentos de Termodinámica

NK NK
Q̇k
∑ ṁki (hj − hki ) = Q̇ Ṡgen = ∑ ṁki (sj − ski ) + >0
Tk
k=1 k=1

Donde hj y sj son la entalpía y entropía de la corriente de salida, hki y ski son la entalpía y
entropía de la corriente de entrada k.

Si el mezclador es adiabático, la temperatura de salida estará entre las más baja y la más
alta de las corrientes de entrada. Esta temperatura permite determinar cuáles corrientes
ceden energía y cuáles ganan energía.
NK NK

∑ ṁki (hj − hki ) = 0 Ṡgen = ∑ ṁki (sj − ski ) > 0


k=1 k=1

Si (hj − hki ) < 0, la corriente k es una corriente caliente, por lo tanto cede energía; en
caso contrario, gana energía.

Es posible que el mezclador elimine calor a los alrededores; en este caso,


NK NK
Q̇ 0
∑ ṁki (hj − hki ) + Q̇0 = 0 Ṡgen = ∑ ṁki (sj − ski ) + >0
T0
k=1 k=1

Si una corriente externa suministra (retira) calor al (del) mezclador, las ecuaciones de
energía y entropía son:
NK NK

∑ ṁki (hj − hki ) + ṁce (hcej − hcei ) = 0; Ṡgen = ∑ ṁki (sj − ski ) + ṁce (scej − scei ) > 0
k=1 k=1

Donde hce y sce son la entalpía y entropía de la corriente externa (j = salida, i = entrada)

Ejemplo. Un mezclador o calentador abierto de una planta generadora de potencia calienta 4 kg/s de agua
a 45° C y 100 kPa mezclándolo con vapor proveniente de una turbina a 100 kPa y 250° C. Asuma que el
flujo de salida es líquido saturado a la presión dada. Si el mezclador es adiabático, determine: a) el flujo
másico del vapor proveniente de la turbina; b) la entropía generada en el proceso.
Solución:
Fluido: Agua
Proceso: Mezclado isobárico y adiabático de dos corrientes (1 y 2) y una salida (3)

Balance de materia: ṁ3 = ṁ1 + ṁ2

Balance de energía: ṁ1 (h3 − h1 ) + ṁ2 (h3 − h2 ) = 0

Balance de entropía: Ṡgen = ṁ1 (s3 − s1 ) + ṁ2 (s3 − s2 ) > 0


Fundamentos de Termodinámica

Los valores de las propiedades de cada corriente fueron obtenidos de las tablas de vapor del agua
(Sonntag-van Wylen) y aparecen en la siguiente tabla.

Tabla de propiedades:
Corriente P (kPa) T (° C) Condición h (kJ/kg) s (kJ/kg-K)
1 100 45 líq sub 188,42 0,6386
2 100 250 vap sob 2974,33 8,0332
3 100 99,62 líq sat 417,44 1,3025

Cálculo del flujo másico del vapor:


(𝐡𝟑 − 𝐡𝟏 ) (𝟒𝟏𝟕, 𝟒𝟒 − 𝟏𝟖𝟖, 𝟒𝟐)
𝐦̇𝟐 = 𝐦̇𝟏 = 𝟒, 𝟎 ∗ = 𝟎, 𝟑𝟓𝟖𝟑 𝐤𝐠/𝐬
(𝐡𝟐 − 𝐡𝟑 ) (𝟐𝟗𝟕𝟒, 𝟑𝟑 − 𝟒𝟏𝟕, 𝟒𝟒)

Cálculo de la entropía generada:


𝐒̇𝐠𝐞𝐧 = 𝟒, 𝟎 ∗ (𝟏, 𝟑𝟎𝟐𝟓 − 𝟎, 𝟔𝟑𝟖𝟔) + 𝟎, 𝟑𝟓𝟖𝟑 ∗ (𝟏, 𝟑𝟎𝟐𝟓 − 𝟖, 𝟎𝟑𝟑𝟐) = 𝟎, 𝟐𝟒𝟒𝟎 𝐤𝐖/𝐊

ANÁLISIS DE LOS GRADOS DE LIBERTAD DE UN PROCESO

A continuación, se analiza un proceso que está integrado por un número determinado de


equipos o dispositivos unidos por un número determinado de corrientes. El análisis de los
grados de libertad permite hallar las variables independientes, las ecuaciones de balance
y las especificaciones para que el proceso pueda resolverse y determinar los
requerimientos energéticos cumpliendo los principios o leyes de la termodinámica.

Considere un proceso conformado por NE equipos y NCO corrientes másicas.

Los equipos o dispositivos pueden ser válvulas, toberas, difusores, turbinas,


compresores, bombas, mezcladores, separadores de fases, separadores de corrientes,
intercambiadores de calor, etc. Cada equipo puede tener dos corrientes energéticas
(Flujo de calor y potencia)

Los balances de materia y energía en cada equipo representan el número de ecuaciones


independientes. Se puede hacer los balances de materia y energía del proceso global y
balances de materia y energía parciales de (NE-1) equipos.

Las corrientes másicas pueden ser de entrada y de salida al proceso o a un equipo. Cada
corriente tiene tres variables (presión, temperatura y flujo másico o volumétrico).

Las variables totales del proceso son: NV = 3 ∗ NCO + 2 ∗ NE


Las ecuaciones independientes son: NEC = 2 ∗ NE
Los datos o especificaciones independientes son: ND = NV − NEC = 3 ∗ NCO

El siguiente esquema representa un proceso donde hay tres (3) corrientes de entrada,
dos (2) corrientes de salida y tres (3) corrientes intermedias que unen uno o varios
equipos. El primer equipo es un mezclador, seguido de un separador de donde salen dos
corrientes. La corriente 5 entra al compresor y sale en el punto 6. La corriente 7 entra a
un intercambiador de calor y sale en el punto 8.
Fundamentos de Termodinámica

mezclador separador
1
compresor
5
2
4
3

6
7
8
intercambiador

Las variables totales del proceso son: NV = 3 ∗ NCO + 2 ∗ NE = 3 ∗ 8 + 2 ∗ 4 = 32


La presión de cada corriente = 8 (P1, P2, P3, P4, P5, P6, P7, P8)
La temperatura de cada corriente = 8 (T1, T2, T3, T4, T5, T6, T7, T8)
El flujo másico de cada corriente = 8 (m1, m2, m3, m4, m5, m6, m7, m8)
El flujo de calor de cada equipo = 4 (QM, QS, QC, QI)
La potencia de cada equipo = 4 (W M, W S, W C, WI)

Las ecuaciones independientes son: NEC = 2 ∗ NE = 2 ∗ 4 = 8


Balance de materia de cada equipo = 4
Balance de energía de cada equipo = 4

Los datos a especificar son: ND = NV − NEC = 3 ∗ NCO = 3 ∗ 8 = 24


Los datos a especificar pueden tener muchas alternativas, por ejemplo:
Flujos de las corrientes de entrada = 3 (m1, m2, m3)
Flujo de calor de tres equipos = 3 (QM, QS, QC)
Las potencias de tres equipos = 3 (W M = 0, WS = 0, W I = 0)
La presión de cada corriente = 8 (P1, P2, P3, P4, P5, P6, P7, P8)
La temperatura de siete corrientes = 7 (T1, T2, T3, T5, T6, T7, T8)

Ejemplo. Un proceso usa energía geotérmica proveniente de agua caliente extraída desde las
profundidades de la tierra para generar potencia en una pequeña turbina. El agua a alta presión a 1,5 MPa
y 180° C se estrangula en una cámara de vaporización flash, la cual forma vapor saturado y líquido saturado
a 40 kPa. El líquido se descarga en el fondo mientras que el vapor sale en el tope y se alimenta a una
turbina adiabática, de donde sale a 10 kPa y calidad del 90%. Si la turbina produce 1 MW, determine: a) el
flujo másico de agua geotérmica; b) la entropía generada en el proceso.

Agua caliente
1,5 MPa
1
180° C
vapor saturado
40 kPa

2
Cámara flash
40 kPa

4 vapor
10 kPa, Calidad 90%
Líquido saturado
40 kPa
Fundamentos de Termodinámica

Solución:
Las variables totales del proceso son: NV = 3 ∗ NCO + 2 ∗ NE = 3 ∗ 4 + 2 ∗ 2 = 16
La presión de cada corriente = 4 (P1, P2, P3, P4)
La temperatura de cada corriente = 4 (T1, T2, T3, T4)
El flujo másico de cada corriente = 4 (m1, m2, m3, m4)
El flujo de calor de cada equipo = 2 (QC, QT)
La potencia de cada equipo = 2 (W C, W T)
Total = 16

Las ecuaciones independientes son: NEC = 2 ∗ NE = 2 ∗ 2 = 4


Balance de materia de cada equipo = 2
Balance de energía de cada equipo = 2
Total = 4

Los datos a especificar son: ND = NV − NEC = 3 ∗ NCO = 3 ∗ 4 = 12


Los datos especificados son:
Flujo de calor de los equipos = 2 (QC = 0, QT = 0)
Las potencias de los equipos = 2 (W C = 0, W T = 1000 kW)
La presión de cada corriente = 4 (P1 = 1500 kPa, P2 = P4 = 40 kPa, P3 = 10 kPa)
La temperatura de la corriente 1 = 1 (T1 = 180° C)
Condición de las corrientes 2, 3, 4) = 3 (x2 = 1, x3 = 0,9 y x4 = 0)
Total = 12

Las siguientes propiedades fueron extraídas de la aplicación ChemicaLogic SteamTab.


A P1 = 1500 kPa y T1 = 180° C, h1 = 763,304 kJ/kg; s1 = 2,13856 kJ/kg-K
A P2 = 40 kPa, vapor saturado, T2 = Tsat =,75,857° C, h2 = hg = 2636,05 kJ/kg; s2 = sg = 7,669 kJ/kg-K
A P4 = 40 kPa, líquido saturado, T4 = Tsat =,75,857° C, h4 = hf = 317,624 kJ/kg; s4 = sf = 1,02607 kJ/kg-K
A P3 = 10 kPa, calidad x3 = 0,90, T3 = Tsat =45,806° C, h3 = 2344,65 kJ/kg; s3 = 7,39886 kJ/kg-K

Balances en la turbina
Balance de materia: m 2 = m3
Balance de energía: W T = m2*(h2 – h3) = 1000 kW; m2 = m3 = 1000/(2636,05 – 2344,65) = 3,4317 kg/s

Balances en la cámara:
Balance de materia: m 2 + m4 = m1
Balance de energía: m 2*h2 + m4*h4 = m1*h1 = (m2 + m4)*h1;
m4 = m2*(h2 – h1)/(h1 – h4) = 3,4317*(2636,05 – 763,304)/(763,304 – 317,624) = 14,378 kg/s
m1 = 3,4317 + 14,378 = 17,81 kg/s
m2 = m3 = 3,4317 kg/s
m4 = 14,378 kg/s

Balance de entropía en el proceso: Sgen = m3*s3 + m4*s4 – m1*s1


Sgen = 3,4317*7,39886 + 14,378*1,02607 – 17,81*2,13856 = 2,056 kW/K

CONSECUENCIAS DE LA SEGUNDA LEY

Los diferentes equipos y procesos vistos anteriormente cumplen los principios de


conservación de materia y energía y el principio del incremento de la entropía del
universo.

Sin embargo, ¿surge el interrogante: cómo determinar si el proceso es eficiente y cuál es


el grado de eficiencia?
Fundamentos de Termodinámica

Considere el siguiente proceso que se lleva a cabo en una turbina de vapor adiabática:

A una turbina adiabática entran 50 kg/s de vapor a 5 MPa y 500° C donde se expande
hasta 120 kPa y a las condiciones que se muestran en la siguiente tabla. Para cada
condición de salida, determine: a) la potencia generada, en kW; b) la entropía generada,
en kW/K.
Se usa la aplicación ChemicaLogic Steam Tab™
A 120 kPa, las propiedades de saturación son:
hg = 2683,05 kJ/kg, hf = 439,359 kJ/kg; sg = 7,29767 kJ/kg-K, sf = 1,36093 kJ/kg-K
Las condiciones en negrillas son las especificaciones en cada caso
Entrada Salida isent Salida 1 Salida 2 Salida 3 Salida 4
P (kPa) 5000 120 120 120 120 120
T (° C) 500 104,784 110
Calidad vap sob 0,9462 0,95 0,98 1,00 vap sob
h (kJ/kg) 3434,66 2562,27 2570,87 2638,18 2683,05 2693,95
s (kJ/kg-K) 6,97808 6,97808 7,00083 7,17894 7,29767 7,3263
W T = m*(h1 - h2) = 43619,7 43189,7 39824,2 37580,5 37035,5
Sgen = m*(s2 - s1) = 0,00 1,14 10,04 15,98 17,41

Los resultados anteriores indican que para cada condición de salida, se satisfacen las
dos leyes de la termodinámica y que la potencia generada por la turbina disminuye a
medida que la entropía generada aumenta. Esto quiere decir que la eficiencia de la
turbina también disminuye. También se puede observar que la máxima potencia
generada ocurre cuando la entropía generada es igual a cero.

Procesos reversibles. Un proceso reversible es aquel en que un sistema puede retornar


a su condición inicial sin dejar huellas o consecuencias en el resto del universo. Considere
1 kg de agua como líquido saturado a 150 kPa (Tsat = 111,35° C) que se vaporiza
isobáricamente hasta vapor saturado con gases de combustión (Cpgc = 1,10 kJ/kg-K) que
entran a 250° C y salen a 120° C. Posteriormente el vapor saturado se condensa
igualmente a presión constante de 150 kPa hasta líquido saturado necesita un fluido de
enfriamiento, por ejemplo agua de enfriamiento (Cpagua = 4,184 kJ/kg-K) desde 35° C
hasta 80° C. Los resultados obtenidos después de aplicar las leyes de la termodinámica
se presentan en la siguiente tabla.
Agua gases de combustión Agua agua de enfriamiento
P (kPa) 150 150 150 150
T (° C) 111,35 111,35 250 120 111,35 111,35 35 80
Condición liq sat vap sat ----- ------ vap sat liq sat ----- ------
Flujo (kg/s) 100 100 1556,63 1556,63 100 100 1182,27 1182,27
Cp (kJ/kg-K) 1,10 1,10 4,184 4,184
h (kJ/kg) 467,13 2693,11 h (kJ/kg) = -143,00 2693,11 467,13 h (kJ/kg) = 188,28
s (kJ/kg-K) 1,4337 7,223 s (kJ/kg-K) = -0,31435 7,223 1,4337 s (kJ/kg-K) = 0,57056
H (kW) 222598,30 -222598,3 -222598,3 222598,3
S (kW/K) 578,93 -489,33 -578,93 674,56
Sgen (kW/K) 89,60 95,63
Fundamentos de Termodinámica

Se puede observar que los cambios de entalpía y entropía del agua son iguales a cero;
es decir que el agua sufre un proceso internamente reversible; sin embargo, los
alrededores (gases de combustión y agua de enfriamiento) no podrían regresar a sus
respectivos estados iniciales con el único contacto entre ellos. Esto quiere decir que el
agua de enfriamiento a 80° C no puede calentar a los gases de 120 a 250° C. El efecto
total es que la entropía generada en este proceso es igual a 185,24 kW/K (89,60 + 95,63)
representando las consecuencias en los alrededores.

Un aumento de la entropía generada implica un aumento de la irreversibilidad y una


disminución de la eficiencia del proceso. Un proceso reversible se considera un proceso
ideal y es el punto de referencia de un proceso irreversible. En el caso de la turbina
analizada anteriormente, la potencia máxima se obtendría si la turbina fuese reversible.
De esta manera, una forma de determinar la eficiencia de un proceso es comparar el
objetivo real con el objetivo que se obtendría si operara reversiblemente.

Eficiencia isentrópica. Los equipos como turbinas, bombas, compresores, toberas y


difusores tienen su respectiva eficiencia isentrópica, como a continuación se relaciona,
teniendo en cuenta que debe cumplirse la condición de la conservación de la entropía del
universo, es decir:
Ṡgen = 0

Turbina. Considere una turbina que tiene una entrada 3 y una salida 4 (según el diagrama
Temperatura-entropía (Ts) mostrado abajo); además hay una pérdida de calor a los
alrededores, Q0. Los balances de energía y entropía son:

ẆT = ṁ(h3 − h4 ) − Q̇ 0
Q̇0
Ṡgen = ṁ(s4 − s3 ) + ≥0
T0

Si la turbina fuese reversible la condición de salida estaría dada por:


Q̇ 0 ⁄ṁ
s4 = s3 −
T0
Obsérvese lo siguiente: a) la entropía de salida puede ser menor que la entropía de
entrada y por tanto, una disminución de la entropía del fluido; b) un alto valor de las
pérdidas de calor a los alrededores causa una disminución de la potencia de la turbina.

Ahora, si la turbina fuese adiabática y además reversible, sería isentrópica y se obtendría


la máxima potencia.
s4 = s3 = constante

ẆTS = ẆTmáx = ṁ(h3 − h4s )

Donde h4s es la entalpía de salida se evalúa a la presión de salida y a entropía constante.


Fundamentos de Termodinámica

Una salida real debe estaría a la derecha del punto 4s y la entalpía sería mayor que h 4s.
En este caso la potencia sería:
ẆT = ṁ(h3 − h4 )

La eficiencia isentrópica o adiabática de la turbina en particular está dada por:

ẆT (h3 − h4 )
ηT = =
ẆTmáx (h3 − h4s )

Si la turbina tiene más de una o dos entradas o salidas, el análisis es el mismo, teniendo
en cuenta el balance de materia.

En el ejemplo de arriba de la turbina con diferentes condiciones de salida, las eficiencias


isentrópicas son:

Salida isent Salida 1 Salida 1 Salida 1 Salida 1


P (kPa) 120 120 120 120 120
T (° C) 104,78 104,78 104,78 104,78 110
Calidad 0,9462 0,95 0,98 1,00 vap sob
W T (kW) 43619,7 43189,7 39824,2 37580,5 37035,5
T = 99,0% 91,3% 86,2% 84,9%

Compresores. Considere un compresor no adiabático de una entrada 1 y salida 2. Los


balances de energía y entropía son:

|ẆC | = ṁ(h2 − h1 ) + Q̇ 0
Q̇0
Ṡgen = ṁ(s2 − s1 ) + ≥0
T0

Obsérvese lo siguiente: a) la entropía de salida puede ser menor que la entropía de


entrada y por tanto, una disminución de la entropía del fluido; b) un alto valor de las
pérdidas de calor a los alrededores causa un aumento de la potencia consumida por el
compresor.
Fundamentos de Termodinámica

Ahora, si el compresor fuese adiabático y además reversible, sería isentrópico y


consumiría la mínima potencia.

s2 = s1 = constante

|ẆCS | = ẆCmín = ṁ(h2s − h1 )

Donde h2s es la entalpía de salida se evalúa a la presión de salida y a entropía constante.

Una salida real debe estaría a la derecha del punto 2s y la entalpía sería mayor que h 2s.

En este caso la potencia sería:


ẆC = ṁ(h2 − h1 )

La eficiencia isentrópica o adiabática del compresor, en particular está dada por:

ẆCmín (h2s − h1 )
ηC = =
ẆC (h2 − h1 )

Tobera. En una tobera reversible y adiabática, el aumento de la energía cinética sería el


máximo y la velocidad de salida igualmente, sería la máxima para las condiciones de
entrada dadas. Los balances de energía y entropía para la tobera isentrópica son:

2
u̅2máx = √2(h1 − h2s ) + u̅12

s2 = s1 = constante

La expresión para la eficiencia isentrópica es la siguiente:

u̇ 22 2(h1 − h2 ) + u̅12
ηtobera = =
u̇ 22máx 2(h1 − h2s ) + u̅12

PROCESOS ISENTRÓPICOS DE GASES IDEALES


El análisis de un proceso isentrópico aplicado a un gas ideal, permite encontrar la relación
entre la temperatura y la presión, a saber:
T2 C
p P2
(s2 − s1 ) = ∫ dT − Rln ( ) = 0
T1 T P1
Si se define, sT° = sT° (T)
TC
p
sT° = ∫ dT
T
y presión relativa, Pr = P⁄P0
Fundamentos de Termodinámica

P2 Pr2 P2 Pr2 sT° 2s − sT° 1


sT° 2s − sT° 1 = Rln ( ) = Rln ( ) ó ( ) = ( ) = exp ( )
P1 Pr1 P1 Pr1 R

De esta manera, la presión relativa en un proceso isentrópico está relacionado con la


temperatura, Pr = Pr(T) o viceversa, T = T(Pr)

Si el Cp es constante y además definiendo k = Cp/Cv, la relación entre la temperatura y la


presión estará dada por:

k−1
T2s P2 T2s P2 k
Cp ln ( ) = Rln ( ) ; ( )= ( )
T1 P1 T1 P1

Por otra parte, la relación entre la temperatura y el volumen se puede hallar a partir del
cambio de la entropía en términos de estas dos variables:
T2
Cv v2
(s2 − s1 ) = ∫ dT + Rln ( ) = 0
T1 T v1
Si Cv es constante,
R
T2 v1 T2 v1 Cv v1 k−1
Cv ln ( ) = Rln ( ) ; ( )=( ) =( )
T1 v2 T1 v2 v2

De la ecuación de los gases ideales se obtiene la relación entre la presión y el volumen:

P2 v1 k
( )=( ) ; Pv k = constante
P1 v2

Ejemplo. A un compresor adiabático se alimentan continuamente 1,4 m 3/s de hidrógeno a 1,2 bar y 320 K.
El hidrógeno sale a 7,2 bar y 545 K. Determine: a) la potencia requerida en el compresor; b) la eficiencia
isentrópica del compresor; c) la entropía generada. Use las tablas del hidrógeno (Tabla A-22 del Cengel)
Solución:
Flujo molar, n1 = P1V1/(RuT1) = 120*1,4/(8,314*320) = 0,06315 kmol/s
Relación de presión, rP = P2/P1 = 7,2/1,2 = 6
Relación isentrópica entre la presión y la temperatura:
P2
sT° 2s − sT° 1 = R u ln ( )
P1
De la tabla A-22, a 320 K, h1 = 9100 kJ/kmol; s°T1 = 132,621 kJ/kmol-K
P2 P2 kJ
sT° 2s − sT° 1 = R u ln ( ) ; sT° 2s = sT° 1 + R u ∗ ln ( ) = 132,621 + 8,314 ∗ ln(6) = 147,518
P1 P1 kmol − K

Interpolando con s°T2s para hallar T2s y h2s, los resultados son: T2s = 533,70 K y h2s = 15336,29 kJ/kmol

De la tabla A-22, a 545 K, h2 = 15667,5 kJ/kmol; s°T2 = 148,13 kJ/kmol-K

Potencia real del compresor: WC = n*(h2 – h1) = 0,06315*(15667,5 – 9100) = 414,74 kW

Potencia isentrópica del compresor: W CS = n*(h2s – h1) = 0,06315*(15336,29 – 9100) = 393,82 kW


Fundamentos de Termodinámica

Eficiencia isentrópica del compresor, C = WCS/WC = 393,82/414,74 = 0,9496 ( 95%)

Cálculo de la entropía generada,

Sgen = n[s°T2 – s°T1 – Ru*ln(P2/P1)] = 0,06315*[148,13 – 132,621 – 8,314*ln(6)] = 0,39kW/K

Procesos politrópicos. Un proceso politrópico es aquel proceso reversible que obedece


las siguientes relaciones:

Pv n = P1 v1n = P2 v2n = constante


n−1
P2 v1 n T2 P2 n T2 v1 n−1
( )=( ) ( )= ( ) ( )=( )
P1 v2 T1 P1 T1 v2

Donde n indica la trayectoria:


0, proceso isobárico
1, proceso isotérmico
n=
k = Cp ⁄Cv , proceso isentrópico
{ ∞, proceso isovolumétrico }

En el diagrama Pv se pueden observar las curvas para varios valores de n. El área bajo
la curva representa el trabajo reversible

Sistemas cerrados: Wrev = ∫ PdV = m ∫ Pdv

Sistemas abiertos: Ẇeje,rev = ∫ V̇dP = ṁ ∫ vdP

Ejemplo. 0,22 kg de CO2 se comprimen en un dispositivo cilindro-émbolo desde una presión de 120 kPa y
volumen de 0,1 m3 hasta 480 kPa y 363,72 K (90,72° C). Si el proceso es politrópico, determine: a) el orden
de la trayectoria; b) el trabajo suministrado, en kJ; c) el calor transferido, en kJ; d) la entropía generada si
los alrededores están a 15° C. Use el Cp del CO2 a 300 K de la Tabla A-2 del libro Termodinamica, Cengel,
6ª edición
P

n=1

n = k = Cp/Cv
n>k

v
1 2
Fundamentos de Termodinámica

Solución:
Constante específica del CO2, R = 0,18895 kPa.m3/kg-K = 0,18895 kJ/kg-K
Para un proceso politrópico se cumple: Pv n = P1 v1n = constante, donde n es el orden de la trayectoria.
El CO2 cumple la ecuación de los gases ideales: PV = mRT
El trabajo de compresión se halla con la integral de PdV y se necesita la relación entre P y V a través de la
trayectoria n.
V2
P2 V2 − P1 V1 mR(T2 − T1 )
W = P1 V1n ∫ V −n dV = = ; n ≠1
V1 1−n 1−n

En el estado inicial, T1 = P1V1/(mR) = 120*0,1/(0,22*0,18895) = 288,68 K


En el estado final, V2 = mRT2/P2 = 0,22*0,18895*363,72/480 = 0,0315 m 3
El orden de la trayectoria se halla con la siguiente expresión a partir de los datos de presión y volumen en
los estados iniciales y finales:
P 120
P1 V2 n ln ( 1 ) ln ( )
P2 480 = 1,2
=( ) ; n= =
P2 V2 V 0,0315
ln ( 2 ) ln ( )
V2 0,1
El trabajo suministrado es:

𝐦𝐑(𝐓𝟐 − 𝐓𝟏 ) 𝟎, 𝟐𝟐 ∗ 𝟎, 𝟏𝟖𝟖𝟗𝟓 ∗ (𝟑𝟔𝟑, 𝟕𝟐 − 𝟐𝟖𝟖, 𝟔𝟖)


𝐖= = = −𝟏𝟓, 𝟔 𝐤𝐉
𝟏−𝐧 𝟏 − 𝟏, 𝟐

Para un pequeño rango de temperatura se puede asumir que Cp es constante.


De la tabla A-2, Cp = 0,846 kJ/kg-K y Cv = 0,657 kJ/kg-K
El cambio de energía interna es:

∆U = mCv (T2 − T1 ) = 0,22 ∗ 0,657 ∗ (363,72 − 288,68) = 10,846 kJ

El calor transferido es: 𝐐 = 𝐖 + ∆𝐔 = −𝟏𝟓, 𝟔 + 𝟏𝟎, 𝟖𝟒𝟔 = −𝟒, 𝟕𝟓𝟒 𝐤𝐉


El calor hacia los alrededores es: Q 0 = 4,574 kJ

La entropía generada se halla con el cambio de entropía del CO 2 y de los alrededores.

T2 P2 363,72 480 kJ
(∆S)CO2 = m [Cp ln ( ) − Rln ( )] = 0,22 ∗ [0,846 ∗ ln ( ) − 0,18895 ∗ ln ( )] = −0,0146
T1 P1 288,68 120 K
4,574 kJ
(∆S)0 = = 0,0159
288 K
𝐤𝐉
(∆𝐒)𝐠𝐞𝐧 = (∆𝐒)𝐂𝐎𝟐 + (∆𝐒)𝟎 = −𝟎, 𝟎𝟏𝟒𝟔 + 𝟎, 𝟎𝟏𝟓𝟗 = 𝟏, 𝟐𝟔 ∗ 𝟏𝟎−𝟑 > 𝟎
𝐊

Ejemplo. Una turbina recibe un flujo estable de 2,5 kg/s de aire a 9 bar y 1500 K que se expande
politrópicamente saliendo a 1,5 bar y un flujo volumétrico de 5,0141 m 3/s.
Considerando Cp del aire constante e igual a 1,004 kJ/kg-K, determine: a) la potencia de la turbina, en kW;
b) el flujo de calor, en kW; c) la entropía generada, en kW/K, si el calor es suministrado desde una fuente
a 1600 K.
Solución:
El aire se considera gas ideal, por lo tanto cumple la ecuación: PV̇ = ṁRT
La turbina y la fuente forman el universo, de tal manera que Q̇ F = −Q̇ rev ; además TF = 1600 K

La relación termodinámica de la entalpía es coherente con la primera ley:


∆h = ∫ Tds + ∫ vdP = qrev − weje,rev
Para un sistema abierto,
∆Ḣ = ∫ ṁCp dT = ∫ TdṠ + ∫ V̇dP = Q̇ rev − Ẇeje,rev
Fundamentos de Termodinámica

La potencia se halla integrando la relación entre la presión y el volumen: Pv n = P1 v1n = C

n−1 n−1
n P2 n n P2 n
Ẇeje,rev = −ṁ ∫ vdP = − P V̇ [( ) − 1] = − ṁRT1 [( ) − 1]
n − 1 1 1 P1 n−1 P1

Estado 1: P1 = 9 bar (900 kPa); T1 = 1500 K; V1 = mRT1/P1 = 2,5*0,287*1500/900 = 1,1958 m3/s


Estado 2: P2 = 1,5 bar (150 kPa); V2 = 5,0141 m3/s; T2 = P2V2/(mR) = 150*5,0141/(2,5*0,287) = 1048,24 K

De la relación entre la temperatura y la presión de un proceso politrópico,


P 150
P2 v1 n ln ( 2 ) ln ( )
P1 900 = 1,25
( )=( ) ; n= v =
P1 v2 ln ( 1 ) ln (1,1958)
v2 5,0141
Cálculo de la potencia reversible:
𝐧−𝟏 𝟏,𝟐𝟓−𝟏
𝐧 𝐏𝟐 𝐧 𝟏, 𝟐𝟓 𝟏 𝟏,𝟐𝟓
𝐖̇𝐞𝐣𝐞,𝐫𝐞𝐯 =− 𝐦̇𝐑𝐓𝟏 [( ) − 𝟏] = − ∗ 𝟐, 𝟓 ∗ 𝟎, 𝟐𝟖𝟕 ∗ 𝟏𝟓𝟎𝟎 [( ) − 𝟏] = 𝟏𝟔𝟐𝟎, 𝟔𝟒 𝐤𝐖
𝐧−𝟏 𝐏𝟏 𝟏, 𝟐𝟓 − 𝟏 𝟔
Cálculo del cambio de entalpía:
∆Ḣ = ∫ ṁCp dT = ṁCp (T2 − T1 ) = 2,5 ∗ 1,004 ∗ (1048,24 − 1500) = −1133,92 kW
Cálculo del flujo de calor:

𝐐̇𝐫𝐞𝐯 = ∆𝐇̇ + 𝐖̇𝐞𝐣𝐞,𝐫𝐞𝐯 = −𝟏𝟏𝟑𝟑, 𝟗𝟐 + 𝟏𝟔𝟐𝟎, 𝟔𝟒 = 𝟒𝟖𝟔, 𝟕𝟐 𝐤𝐖

Cálculo de la entropía generada:


𝐓𝟐 𝐏𝟐 𝐐̇𝐅 𝟏𝟎𝟒𝟖, 𝟐𝟒 𝟏 𝟒𝟖𝟔, 𝟕𝟐
𝐒̇𝐠𝐞𝐧 = 𝐦̇ [𝐂𝐩 𝐥𝐧 ( ) − 𝐑𝐥𝐧 ( )] + = 𝟐, 𝟓 ∗ [𝟏, 𝟎𝟎𝟒 ∗ 𝐥𝐧 ( ) − 𝟎, 𝟐𝟖𝟕 ∗ 𝐥𝐧 ( )] −
𝐓𝟏 𝐏𝟏 𝐓𝐅 𝟏𝟓𝟎𝟎 𝟔 𝟏𝟔𝟎𝟎
𝐤𝐖
𝐒̇𝐠𝐞𝐧 = 𝟎, 𝟎𝟖𝟐
𝐊

PROCESOS DE COMPRESIÓN MULTIETAPAS CON INTERENFRIAMIENTO


Es frecuente realizar trabajo de compresión con relaciones de presiones altas; si el
proceso es adiabático, la temperatura de descarga es elevada. Una manera de alcanzar
altas relaciones de presión y bajas temperaturas es la compresión en varias etapas con
enfriamiento entre etapas. El objetivo del enfriamiento es disminuir la temperatura antes
de la siguiente etapa. La relación de presión de cada etapa se determina a partir de la
presión mínima y la presión máxima.

Considere una compresión de N etapas desde una presión mínima de entrada P0 y


temperatura de entrada mínima T0, hasta una presión máxima de salida PN, la presión de
descarga de la etapa j (j = 1, N) se obtiene de tal manera que la relación de presión es
igual para cada etapa con el fin de consumir el mínimo trabajo de compresión.

N
j N−j
Pj = √(PN ∗ P0 ) j = 1, 2, 3, N

La relación de presión de cualquiera etapa está dada por:


Fundamentos de Termodinámica

N PN
rP = √
P0

Si el fluido se enfría hasta la temperatura de entrada de primera etapa, es decir T0, la


temperatura de salida también es la misma y la potencia de compresión de las etapas
son iguales. Como consecuencia la potencia total y el flujo de calor total eliminado serán,

- Con Cp constante:

|Ẇtotal | = N|Ẇetapaj | = N ∗ ṁCp T0 [(rP )(k−1)/k − 1]

Q̇total = (N − 1)Q̇etapaj = (N − 1) ∗ ṁCp T0 [(rP )(k−1)/k − 1]

- Con Cp variable:

|Ẇtotal | = N|Ẇetapaj | = N ∗ ṁ[hj − h0 ]

Q̇total = (N − 1)Q̇etapaj = (N − 1) ∗ ṁ[hj − h0 ]

Donde hj es la entalpía de salida de la etapa j que se halla a partir de s° Tj

sT° j − sT° 0 = R u ln(rP )

Ejemplo. 75 moles/s de hidrógeno se comprime desde 1 bar y 300 K hasta 16 bar en cuatro etapas
isentrópicas con enfriamiento entre etapas de compresión. Cada etapa de enfriamiento baja la temperatura
del hidrógeno hasta 300 K antes de entrar a la siguiente etapa de compresión. Usando las tablas del
hidrógeno (Tabla A-22 del Cengel), determine: a) la potencia total del compresor, en kW; b) la temperatura
de salida de cada etapa de compresión, en K; c) el flujo de calor total eliminado, en kW
Solución:
Cálculo de la relación de presión de cada etapa:
N PN 4 16,0
rP = √ =√ =2
P0 1,0
T P4 P3 P2
P1 P0

Tj

T0

s
Fundamentos de Termodinámica

De la Tabla A-22 del Cengel, a 300 K, h0 = 8522 kJ/kmol; s°T0 = 130,754 kJ/kmol-K
Cálculo de s°Tj:
1. sT° j = sT° 0 + R u ln(rP ) = 130,754 + 8,314 ∗ ln(2) = 136,517 kJ/kmol − K

La entalpía y la temperatura de salida cada etapa se halla interpolando con s°Tj. Los resultados son:
hj = 10438,55 kJ/kmol, Tj = 366,08 K

Cálculo de la potencia de cada etapa:


kmol kJ
|Ẇetapaj | = ṅ [hj − h0 ] = 0,075 ∗ (10438,55 − 8522,00) = 143,74 kW
s kmol

Cálculo de la potencia total:


|𝐖̇𝐭𝐨𝐭𝐚𝐥 | = 𝐍|𝐖̇𝐞𝐭𝐚𝐩𝐚𝐣 | = 𝟒 ∗ 𝟏𝟒𝟑, 𝟕𝟒 = 𝟓𝟕𝟒, 𝟗𝟔 𝐤𝐖

Cálculo del flujo de calor total eliminado:


𝐐̇𝐭𝐨𝐭𝐚𝐥 = (𝐍 − 𝟏)𝐐̇𝐞𝐭𝐚𝐩𝐚𝐣 = 𝟑 ∗ 𝟏𝟒𝟑, 𝟕𝟒 = 𝟒𝟑𝟏, 𝟐𝟐 𝐤𝐖

Ejemplo. Repita el ejemplo anterior pero la compresión isentrópica se lleva a cabo en una etapa.
Solución:
16,0
Cálculo de la relación de presión: rP = = 16
1,0
Cálculo de s°Tj:
sT° j = sT° 0 + R u ln(rP ) = 130,754 + 8,314 ∗ ln(16) = 153,805 kJ/kmol − K

La entalpía y la temperatura de salida se hallan interpolando con s°Tj. Los resultados son:
hj = 19075,19 kJ/kmol, Tj = 661,14 K

Cálculo de la potencia total (una etapa):


𝐤𝐦𝐨𝐥 𝐤𝐉
|𝐖̇𝐭𝐨𝐭𝐚𝐥 | = 𝐧̇ [𝐡𝐣 − 𝐡𝟎 ] = 𝟎, 𝟎𝟕𝟓 ∗ (𝟏𝟏𝟗𝟎𝟕𝟓, 𝟏𝟗 − 𝟖𝟓𝟐𝟐, 𝟎𝟎) = 𝟕𝟗𝟏, 𝟒𝟗 𝐤𝐖
𝐬 𝐤𝐦𝐨𝐥

De los resultados anteriores se puede concluir que la compresión multietapa consume menos potencia que
la compresión en una sola etapa. Si el número de etapas de compresión tiende a infinito, la compresión
tiende a ser isotérmica y la potencia mínima se halla con la siguiente expresión:

𝐏𝐍 𝟏𝟔
|𝐖̇𝐢𝐬𝐨𝐭 | = 𝐧̇ (𝐑 𝐮 𝐓𝟎 )𝐥𝐧 ( ) = 𝟎, 𝟎𝟕𝟓 ∗ 𝟖, 𝟑𝟏𝟒 ∗ 𝟑𝟎𝟎 ∗ 𝐥𝐧 ( ) = 𝟓𝟏𝟖, 𝟔𝟓 𝐤𝐖
𝐏𝟎 𝟏
Fundamentos de Termodinámica

EFECTOS TÉRMICOS EN SISTEMAS REACCIONANTES

En los procesos que presentan reacciones químicas, las leyes de la termodinámica son
las mismas pero las propiedades como entalpía, energía interna y entropía deben tener
valores absolutos por cuanto las moles de las especies químicas no se conservan.

Considere la siguiente reacción química que ocurre en un reactor:

aA + bB ⇔ rR + sS

Donde A y B son especies químicas reactivas, R y S son especies químicas productos;


a, b, r y s son los coeficientes estequiométricos (positivos para productos y negativos
para los reactivos).

Una especie química puede ser un inerte, el cual no interviene en la reacción y su


coeficiente estequiométrico en dicha reacción es igual a cero.

El alimento al reactor está constituido de nAi moles de A, nBi moles de B, nRi moles de R,
nSi moles de S (téngase presente que el subíndice i es para la entrada). Las moles
alimentadas de R y S pueden ser iguales a cero.

Las moles totales a la entrada es la suma de las moles de todas las especies químicas
a la entrada del reactor expresada como ni:
NEQ

ni = ∑ nki = nAi + nBi + nRi + nSi


k=1

NEQ = número de especies químicas

El grado de conversión , se define como el alcance de la reacción y es una medida


del grado en que se convierte un reactivo. Por ejemplo, si A es el reactivo límite, el
máximo valor de  es niA.

El balance de materia de cada especie química k, en una reacción simple como la de


arriba queda expresado como:

nkj = nki + νk α

Donde:
nkj , nki = moles a la salida y a la entrada de la especie k
νk = coeficiente estequiométrico de la especie química k en la reacción química
α = grado de conversión de la reacción expresado en moles

Balance de A: nAj = nAi − aα


Balance de B: nBj = nBi − bα
Fundamentos de Termodinámica

Balance de R: nRj = nRi + rα


Balance de S: nSj = nSi + sα

Las moles totales a la salida del reactor nj es la suma de todas las especies química que
salen del reactor:
NEQ NEQ

nj = ∑ nkj = nAj + nBj + nRj + nSj = ni + ∑ νk α = ni + Δnα


k=1 k=1

Δn es el cambio de moles de la reacción y puede ser menor que cero, cero o mayor que
cero.

Si las reacciones ocurren en fase gaseosa homogénea las fracciones molares de cada
especie química a la entrada y a la salida del reactor, se pueden determinar,
respectivamente con:
nki nkj
yki = ykj =
ni nj

La conversión fraccional de un reactivo A, por ejemplo, es la relación entre sus moles


convertidas y sus moles alimentadas y expresa como,

(nAi − nAj ) aα
XA = = 0 ≤ XA ≤ 1
nAi nAi

El grado de conversión de la reacción está relacionado con la conversión fraccional de la


especie química A, a partir de la ecuación anterior:

nAi XA
α=
a

Si el reactivo A es un reactivo límite, su conversión fraccional puede ser del 100%

Ejemplo. La siguiente reacción ocurre en fase gaseosa en un reactor continuo:

CO2(g) + 3H2(g) ⇔ CH3 OH(g) + H2 O(g)

Al reactor se alimentan 1 kmol/h de CO 2 y 4 kmol/h de H2. La conversión fraccional del CO2 es del 75%.
Calcule la composición molar de los productos y la producción de metanol, en kmol/h.
Solución:
Sea A: CO2(g) B : H2(g) R : CH3OH(g) S : H2O(g)
Los flujos en la entrada, en kmol/h son: ṅ Ai = 1,0 ṅ Bi = 4,0 ṅ Rj = 0 ṅ Sj = 0
El grado de conversión se halla a partir de la conversión fraccional del CO2:
ṅ Ai XA
α= = XA = 0,75 kmol/h
a
El balance de materia por especie química por especie en kmol/h, es el siguiente:
nAj = 1 − α = 0,25 nBj = 4 − 3α = 1,75 nRj = 0 + α = 0,75 nSj = 0 + α = 0,75
Las moles totales a la salida son: nj = ni + Δnα = 5 − 2 ∗ 0,75 = 3,5 kmol/h
Fundamentos de Termodinámica

La composición molar de los productos aparece en la siguiente tabla:


Especie química
CO2(g) H2(g) CH3OH(g) H2O(g) Sumas
nki (kmol/h) 1,0 4,0 0,0 0,0 5,00
k -1 -3 +1 +1 -2
nkj (kmol/h) 0,25 1,75 0,75 0,75 3,50
ykj 0,0714 0,5000 0,2143 0,2143 1,000

La producción de metanol es de 0,75 kmol/h

Los balances de energía y entropía de sistemas reaccionantes se basan en las


propiedades termodinámicas absolutas.

Debido a que los reactivos y los productos son mezclas de especies químicas, la entalpía
y entropía absoluta de la mezcla están en función de la temperatura, la presión y la
composición molar.
NEQ NEQ

Ḣ = ∑ ṅ k h̅k Ṡ = ∑ ṅ k s̅ k
k=1 k=1

Donde h̅k y h̅k son entalpía y entropía absoluta de la especie química k expresadas como:
T T
Cp dT
h̅k = ∆h°f,298 + ∫ Cp dT s̅ k = s°298 + ∫ − R u ln(pk )
T
298 298

∆h°f,298 = Entalpía de formación de la especie química a 298 K y 1 bar


s°298 = Entropía absoluta de la especie química a 298 K y 1 bar
pk = yk P = Presión parcial de la especie química, en bar

La entalpía de formación de una especie química k, es la energía necesaria para


formarlas a partir de sus elementos estables. Los elementos estables como el C(s), S(s),
H2(g), O2(g), N2(g), Cl2(g) están libres en la naturaleza y su entalpía de formación es cero.
Por ejemplo, para formar un mol de CO2 a 298 K y 1 bar, la reacción es:

C(s) + O2(g) ⇒ CO2)g) ∆h°kf,298 = −393520 J/mol

La energía que se libera cuando se quema un mol de carbono sólido con un mol de
oxígeno para formar un mol de dióxido de carbono a 298 K y 1 bar, es de 393520 J.

Reemplazando las entalpías y entropías absolutas de las especies químicas en la


correspondiente propiedad de una mezcla en un sistema reaccionante, las expresiones
son las siguientes:

- Mezcla de NEQ especies químicas que entra al reactor (en este caso, un reactor
continuo) a temperatura Ti, presión Pi y flujos nki,
Fundamentos de Termodinámica

NEQ Ti NEQ

Ḣi = ∑(ṅ ki ∆h°kf,298 ) + ∫ [∑(ṅ ki Cpk )] dT


k=1 298 k=1

NEQ Ti NEQ NEQ


dT
Ṡi = ∑(ṅ ki s°298 ) + ∫ [∑(ṅ ki Cpk )] − R u ∑(ṅ ki ln(yki Pi ))
T
k=1 298 k=1 k=1

- Mezcla de NEQ especies químicas que sale al reactor (en este caso, un reactor
continuo) a temperatura Tj, presión Pj y flujos nkj,

NEQ Tj NEQ

Ḣj = ∑(ṅ kj ∆h°kf,298 ) + ∫ [∑(ṅ kj Cpk )] dT


k=1 298 k=1

NEQ Tj NEQ NEQ


dT
Ṡj = ∑(ṅ kj s°298 ) + ∫ [∑(ṅ kj Cpk )] − R u ∑ (ṅ kj ln(ykj Pj ))
T
k=1 298 k=1 k=1

Primera Ley aplicada a sistemas reaccionantes. El balance de energía de un sistema


reaccionante queda expresado así:

NEQ Tj NEQ Ti NEQ

Q̇ = Ḣj − Ḣi = ∑(Δṅ k ∆h°kf,298 ) + ∫ [∑(ṅ kj Cpk )] dT − ∫ [∑(ṅ ki Cpk )] dT


k=1 298 k=1 298 k=1

Donde:
Δṅ k = ṅ kj − ṅ ki = νk α = cambio de moles de la especie química k

> 0, si k es producto
Δṅ k = { = 0, si k es un inerte }
< 0, si k es un reactivo
NEQ NEQ

∑(Δṅ k ∆h°kf,298 ) = α ∑(νk ∆h°kf,298 ) = ∆Ḣ298


°

k=1 k=1

 Si Cpk es función de la temperatura, Cpk = ak + bk T + ck T 2 + dk T 3 (Cengel, van


Wylen)
Fundamentos de Termodinámica

NEQ NEQ NEQ NEQ NEQ

∑(ṅ k Cpk ) = ∑(ṅ k ak ) + ∑(ṅ k bk ) T + ∑(ṅ k ck ) T 2 + ∑(ṅ k dk ) T 3


k=1 k=1 k=1 k=1 k=1

Cpk
 Si Cpk es función de la temperatura, = ak + bk T + ck T 2 + dk /T 2 (Smith van
R
Ness)

NEQ NEQ NEQ NEQ NEQ

∑(ṅ k Cpk ) = R [∑(ṅ k ak ) + ∑(ṅ k bk ) T + ∑(ṅ k ck ) T 2 + ∑(ṅ k dk ) /T 2 ]


k=1 k=1 k=1 k=1 k=1

 Si se dispone de los valores de la entalpía en función de la temperatura absoluta,


Tj NEQ NEQ Ti NEQ NEQ

∫ [∑(ṅ kj Cpk )] dT = ∑ ṅ kj (hkj − hk298 ) ; ∫ [∑(ṅ ki Cpk )] dT = ∑ ṅ ki (hki − hk298 )


298 k=1 k=1 298 k=1 k=1

Donde hkj , hki y hk298 son las entalpías de la especie k evaluadas a T j (temperatura de
los productos), Ti (temperatura de los reactivos o alimento) y 298 K, respectivamente.

Se denomina calor o entalpia estándar de reacción a 298 K, al cambio de entalpía


cuando se convierte completamente un kmol de un reactivo de coeficiente
estequiométrico igual a menos uno, a 298 K y 1 bar; en este caso T i = Tj = 298 K. La
ecuación general de energía para 1 kmol de un reactivo se reduce a:

NEQ

∆h°298 = ∑(νk ∆h°kf,298 )


k=1

Ejemplo. Determine la entalpía estándar de la reacción del ejemplo anterior si el CO 2 se convierte


completamente. Repita el problema si el metanol y el agua salen en fase líquida. Compare los resultados.
Solución.
Los datos de entalpia de formación de las especies químicas se tabulan en la siguiente tabla.

Especie química
CO2(g) H2(g) CH3OH(g) H2O(g)
k -1 -3 +1 +1
h°kf,298 (kJ/kmol) -393520 0 -200670 -241820
Fuente: Termodinámica de Cengel

Cálculo de la entalpía estándar de la reacción de conversión de 1 kmol de CO 2 para formar 1 kmol de


metanol gaseoso:
NEQ
kJ
∆h°298 = ∑(νk ∆h°kf,298 ) = (−1)(−393520) + (−3)(0) + (1)(−200670) + (1)(−241820) = −48970
kmol
k=1

Si el metanol y el agua salen en fase líquida, la reacción es:


Fundamentos de Termodinámica

CO2(g) + 3H2(g) ⇔ CH3 OH(l) + H2 O(l)

Ahora las entalpías de formación del metanol y el agua serán diferentes debido al cambio de fase.

∆h°kf,298(l) = ∆h°kf,298(g) − hfg,25° C

Los nuevos valores de las entalpías de formación para el metanol y el agua aparecen en la tabla de abajo.

Especie química
CO2(g) H2(g) CH3OH(l) H2O(l)
k -1 -3 +1 +1
h°kf,298 (kJ/kmol) -393520 0 -238660 -285830
Fuente: Termodinámica de Cengel

Cálculo de la entalpía estándar de la reacción de conversión de 1 kmol de CO 2 para formar 1 kmol de


metanol líquido:
NEQ
kJ
∆h°298 = ∑(νk ∆h°kf,298 ) = (−1)(−393520) + (−3)(0) + (1)(−238660) + (1)(−285830) = −130970
kmol
k=1

La diferencia en los resultados se debe a que hay que retirar el calor de vaporización al metanol y al agua
para condensarlos.

Ecuación de energía de casos especiales:


 Temperatura de los reactivos o alimento y de los productos iguales a 298 K
(Ti = Tj = 298 K). Caso de un reactor isotérmico a 298 K

NEQ NEQ

Q̇ = ∑(Δṅ k ∆h°kf,298 ) = α ∑(νk ∆h°kf,298 ) = ∆Ḣ298


°

k=1 k=1

 Temperatura de los reactivos o alimento y de los productos iguales


(Ti = Tj = T  298 K). Caso de un reactor isotérmico a una temperatura T

NEQ T NEQ

Q̇ = α {∑(νk ∆h°kf,298 ) + ∫ [∑(νk Cpk )] dT}


k=1 298 k=1

 Temperatura de los reactivos o alimento igual a 298 K


(Ti = 298 K, Tj  298 K). Caso de un reactor no isotérmico

NEQ Tj NEQ

Q̇ = α ∑(νk ∆h°kf,298 ) + ∫ [∑(ṅ kj Cpk )] dT


k=1 298 k=1

 Temperatura de los reactivos o alimento y de los productos diferentes a 298 K


(Ti  Tj  298 K). Caso de un reactor no isotérmico
Fundamentos de Termodinámica

NEQ Tj NEQ Ti NEQ

Q̇ = α ∑(νk ∆h°kf,298 ) + ∫ [∑(ṅ kj Cpk )] dT − ∫ [∑(ṅ ki Cpk )] dT


k=1 298 k=1 298 k=1

Ejemplo. La siguiente reacción ocurre en un reactor continuo en fase gaseosa ideal:


CO2(g) + 3H2(g) ⇔ CH3 OH(g) + H2 O(g)
Al reactor entran 1 kmol/h de CO2 y 4 kmol/h de H2 a 298 K. Los productos salen a 400K. Si la conversión
fraccional del CO2 es de 95%, determine el flujo de calor transferido, en kJ/h. Los datos de entalpía del
CO2, H2 y H2O se pueden obtener de las tablas A-17 a A-22 del Cengel. Para el metanol, el Cp está dado
por: Cp = 19,0 + 9,152 ∗ 10−2 T − 1,22 ∗ 10−5 T 2 − 8,039 ∗ 10−9 T 3
Solución:
 (kmol/h) = 0,95 Tj (K) = 400 T (K)
Especie química
CO2(g) H2(g) CH3OH(g) H2O(g)
Sumas
(A) (B) (R) (S)
nki (kmol/h) 1 4 0 0 5,00 ni (kmol/h)

k -1 -3 1 1 -2 Δn
nkj (kmol/h) 0,05 1,15 0,95 0,95 3,10 nj (kmol/h)
ykj 0,0161 0,3710 0,3065 0,3065 1,00 ∑(ykj )
Δh°kf,298 (kJ/kmol) -393520 0 -200670 -241820 -48970,00 Δh°298
h°k,298 (kJ/kmol) 9364 8468 9602 9904
h°k,400 (kJ/kmol) 13372 11426 14610 13356
Δh°k (kJ/kmol) 4008 2958 5008 3452 11639,10 ∑(nkj Δh°k )
Q (kJ/h) -34882,40

El balance de materia por especie química por especie en kmol/h, es el siguiente:


nAj = 1 − α = 0,05 nBj = 4 − 3α = 1,15 nRj = 0 + α = 0,95 nSj = 0 + α = 0,95
El cambio de entalpía del metanol es: hRj − hR298 = 5008 kJ/kmol
NEQ

∑(νk ∆h°kf,298 ) = ∆h°298 = −48970,0 kJ/kmol


k=1
NEQ

Q̇ = α∆h°298 + ∑ ṅ kj (hkj − hk298 )


k=1
Q̇ = 0,95 ∗ (−48970,0) + 0,05 ∗ 4008 + 1,15 ∗ 2958 + 0,95 ∗ 5008 + 0,95 ∗ 3452 = −34882,4 kJ/h

Ejemplo. Si el reactor del ejemplo anterior operara isotérmicamente a 400 K (Ti = Tj = 400 K), determine el
flujo de calor.
Solución:
En este caso la ecuación de energía será:
NEQ

Q̇ = α [∆h°298 + ∑ νk (hkj − hk298 )]


k=1

kJ
Q̇ = 0,95 ∗ [(−48970,0) + (−1) ∗ 4008 + (−3) ∗ 2958 + (1) ∗ 5008 + (1) ∗ 3452] = −50722,4
h
Fundamentos de Termodinámica

Segunda Ley aplicada a sistemas reaccionantes. El balance de entropía de un


sistema reaccionante queda expresado así:

NEQ Tj NEQ Ti NEQ


dT dT
(ΔṠ)reacción = ∑ Δṅ k s °298 + ∫ [∑ ṅ kj Cpk ] − ∫ [∑ ṅ ki Cpk ]
T T
k=1 298 k=1 298 k=1
NEQ NEQ

− R u [ṅ j ln(Pj ) − ṅ i ln(Pi )] − R u [∑ ṅ kj ln(ykj ) − ∑ ṅ ki ln(yki )]


k=1 k=1

El cambio de entropía del depósito térmico (DT) con el cual el sistema intercambia calor,
es (Q̇DT = −Q̇):
Q̇ DT
(ΔṠ)DT =
TDT

La entropía generada es:


Q̇DT
Ṡgen = (ΔṠ)reacción + ≥0
TDT
Ejemplo. A un reactor continuo isotérmico que opera a 550° C y 10 atm (10,1325 bar), entran 1 kmol/h de
CO y 3 kmol/h de vapor de agua en donde reaccionan para formar H2 y CO2 en fase gaseosa ideal:

CO + H2O <=====> CO2 + H2

Calcular el flujo de calor y la entropía generada, si la conversión fraccional del CO es de 70% y el calor es
eliminado a los alrededores que están a TDT = 300 K

Solución: (ver la Hoja H2 - 550°C del archivo Excel PROBLEMAS DE EFECTOS TERMICOS 20172.xls)
Fuente bibliográfica: Introducción a la Termodinámica en Ingeniería Química. Smith-van Ness

La capacidad térmica de las especies químicas están dadas por: Cpk = R[a k + bk T + ck T 2 + dk /T 2 ]

Datos de entalpía de formación (kJ/kmol), constantes de Cpk/R y entropía absoluta (kJ/kmol-K)


h°f298 ak bk ck dk s°298
CO(g) -110525 3,376 5,57*10-4 0 -3,10*103 197,65
H2O(g) -241820 3,470 1,45*10-3 0 1,21*104 188,83
-4 3
H2(g) 0 3,249 4,22*10 0 8,30*10 130,68
CO2(g) -393520 5,457 1,05*10-3 0 -1,16*105 213,80

Cálculo del grado de conversión,


ṅ Ai XA 1,0 ∗ 0,70 kmol
α= = = 0,70
a 1 h

Balance de materia para cada especie química: nkj = nki + νk α

En la siguiente tabla se presentan los resultados del balance de materia:


Fundamentos de Termodinámica

k nki (kmol/h) nkj (kmol/h) nk (kmol/h) yki = nki/ni ykj = nkj/nj
CO(g) -1 1,00 0,30 -0,70 0,250 0,075
H2O(g) -1 3,00 2,30 -0,70 0,750 0,575
H2(g) 1 0,00 0,70 0,70 0,000 0,175
CO2(g) 1 0,00 0,70 0,70 0,000 0,175
Sumas 0 ni = 4,00 nj = 4,00 1,000 1,000

Balance de energía:

NEQ Tj NEQ Ti NEQ

Q̇ = α ∑(νk ∆h°kf,298 ) + ∫ [∑(ṅ kj Cpk )] dT − ∫ [∑(ṅ ki Cpk )] dT


k=1 298 k=1 298 k=1

NEQ

α ∑(νk ∆h°kf,298 ) = (∆nk )T (∆h°f298 ) = −28822,50 kJ/h


k=1

Para los reactivos:

NEQ NEQ NEQ NEQ

∑(ṅ ki Cpk ) = R [∑(ṅ ki a k ) + ∑(ṅ ki bk ) T + ∑(ṅ ki dk ) /T 2 ] = R[13,786 + 4,907 ∗ 10−3 T + 3,32 ∗ 104 /T 2 ]
k=1 k=1 k=1 k=1

823 NEQ

∫ [∑(ṅ ki Cpk )] dT = 72769,66 kJ/h


298 k=1

Para los productos:

NEQ NEQ NEQ NEQ

∑(ṅ kj Cpk ) = R [∑(ṅ kj a k ) + ∑(ṅ kj bk ) T + ∑(ṅ kj dk ) /T 2 ] = R[15,088 + 4,529 ∗ 10−3 T − 4,83 ∗ 104 /T 2 ]
k=1 k=1 k=1 k=1

823 NEQ

∫ [∑(ṅ kj Cpk )] dT = 76077,80 kJ/h


298 k=1

𝐤𝐉 𝐤𝐉
𝐐̇ = −𝟐𝟖𝟖𝟐𝟐, 𝟓𝟎 + 𝟕𝟔𝟎𝟕𝟕, 𝟖𝟎 − 𝟕𝟐𝟕𝟔𝟗, 𝟔𝟔 = −𝟐𝟓𝟓𝟏𝟒, 𝟑𝟔 ; 𝐐̇𝐃𝐓 = −𝐐̇ = 𝟐𝟓𝟓𝟏𝟒, 𝟑𝟔
𝐡 𝐡

Balance de entropía:

Cambio de entropía de la reacción:


NEQ Tj NEQ Ti NEQ
dT dT
(ΔṠ)reacción = ∑ Δṅ k s °298 + ∫ [∑ ṅ kj Cpk ] − ∫ [∑ ṅ ki Cpk ] − R u [ṅ j ln(Pj ) − ṅ i ln(Pi )]
T T
k=1 298 k=1 298 k=1
NEQ NEQ

− R u [∑ ṅ kj ln(ykj ) − ∑ ṅ ki ln(yki )]
k=1 k=1
Fundamentos de Termodinámica

NEQ Tj NEQ Ti NEQ


dT dT kJ
∑ Δṅ k s °298 = (∆nk )T (s298
° )
= −29,40; ∫ [∑ ṅ kj Cpk ] = 145,24; ∫ [∑ ṅ ki Cpk ] = 139,20
T T h−K
k=1 298 k=1 298 k=1

R u [ṅ j ln(Pj ) − ṅ i ln(Pi )] = 77,01 − 77,01 = 0

NEQ NEQ
kJ
R u [∑ ṅ kj ln(ykj ) − ∑ ṅ ki ln(yki )] = (−37,33) − (−18,70) = −18,63
h−K
k=1 k=1

kJ
(ΔṠ)reacción = −29,40 + (145,24 − 139,20) − (−18,63) = −4,74
h−K

Cambio de entropía de los alrededores:

Q̇ DT 25514,36 kJ
(ΔṠ)DT = = = 85,05
TDT 300 h−K

Cálculo de la entropía generada:

𝐤𝐉
𝐒̇𝐠𝐞𝐧 = −𝟒, 𝟕𝟒 + 𝟖𝟓, 𝟎𝟓 = 𝟖𝟎, 𝟑𝟏 >𝟎
𝐡−𝐊

Ejemplo. El gas amoniaco entra al reactor de una planta de ácido nítrico mezclado con 25% más de aire
seco del que se requiere para la conversión completa de amoniaco en óxido nítrico (NO) y vapor de agua.
Los reactivos entran a 75° C. Si la conversión del amoniaco es de 80% y la reacción es adiabática,
determine la temperatura de los gases que salen del reactor. Suponga gases ideales y C p constante.
La reacción que ocurre en fase gaseosa ideal y es la siguiente:

NH3(g) + (5/4)O2(g) --------> NO(g) + (3/2)H2O(g)

Datos de entalpía de formación a 298 K y Cp/R


NH3(g) O2(g) N2(g) NO(g) H2O(g)
h°f298 (kJ/kmol) -46110 0 0 90250 -241818
Cpk/R 4,269 3,535 3,502 3,590 4,038

Solución. (ver la Hoja NO – NH3 del archivo Excel PROBLEMAS DE EFECTOS TERMICOS 20172.xls)

Ejemplo. Un gas que consiste de solo CO y N2 se obtiene al pasar una mezcla de gases de combustión y
aire a través de un lecho de carbón incandescente (suponga que el carbón es puro). Las dos reacciones
que ocurren en forma completa son:

CO2(g) + C(s) ⇒ 2CO(g) y 2C(s) + O2(g) ⇒ 2CO(g)

El gas de combustión tiene la siguiente composición molar: 12,80% de CO 2, 3,70% de CO, 5,40% de O2 y
78,10% de N2. La mezcla de gas de combustión y aire se proporciona de tal manera que los calores de las
reacciones se cancelan y en consecuencia la temperatura del lecho de carbón es constante. Si esta
temperatura es 867° C (1140 K), si la corriente de alimentación se precalienta hasta 867° C y si el proceso
es adiabático, determine: a) la relación entre las moles de gas de combustión y aire que se requiere; b) la
composición molar del gas producto.
Solución:
Favor consultar la Hoja CO-N2 del archivo Excel PROBLEMAS DE EFECTOS TERMICOS 20172
Fundamentos de Termodinámica

Base de cálculo = 100 kmol/h de gases de combustión


Defínase, a = kmol/h de O2 y 3,76a kmol/de N2 para 3,76a kmol/h de aire seco.
CpC(s) = 8,496 kJ/kmol-K
Ecuación de energía:
NEQ NEQ

Q̇ = ∑ Δnk ∆h°kf,298 + ∑ Δnk (hT − h298 ) = 0


k=1 k=1
Balance de materia de las especies químicas:

nCO2 j = 12,8 − α1 = 0; nCOj = 3,7 + 2α1 + 2α2 ; nO2j = 5,4 + a − α2 = 0; nN2j = 78,1 + 3,76a

α1 = 12,8; α2 = 5,4 + a

∆nCO2 = 0; ∆nCO = 2α1 + 2α2 = 36,4 + 2a; ∆nO2 = 0; ∆nN2 = 0; ∆nC = −α1 − 2α2 = −(23,6 + 2a)

En la siguiente tabla aparecen los datos, los cálculos y resultados


nki k1 k2 nkj nk h°f298 h298 hTj hk nkj
CO2(g) 12,8 -1 0 0 -12,8 -393520 9364 50484 41120,00 0
CO(g) 3,7 2 2 40,1 + 2a 36,4 + 2a -110530 8669 35054 26385,00 50,85
O2(g) 5,4 + a 0 -1 0 -(5,4 + a) 0 8682 36314 27632,00 0
N2(g) 78,1 + 3,76a 0 0 78,1 + 3,76a 0 0 8669 34760 26091,00 98,31
C(s) --- -1 -2 ---- -(23,6 + 2a) 0 0 7153,63 7153,63 0
149,16

Reemplazando en la ecuación de energía,


(−12,8)(−393520 + 41120) + (36,4 + 2a)(−110530 + 26385) − (5,4 + 2a)(0 + 27632)
− (23,6 + 2a)(0 + 7153,63) = 0

Resolviendo, a = 5,374 kmol/h Flujo molar de aire seco: ṅ aire seco = 4,76 ∗ 5,374 = 25,58 kmol/h

Los resultados son:

1) Relación gas/aire: 𝐑 = 𝟏𝟎𝟎⁄𝟐𝟓, 𝟓𝟖 = 𝟑, 𝟗𝟏

2) Composición molar del gas producto: 𝐲𝐂𝐎 = 𝟎, 𝟑𝟒𝟏; 𝐲𝐍𝟐 = 𝟎, 𝟔𝟓𝟗

Ejemplo. Una caldera genera 36000 kg/h de vapor a 5 MPa y 500° C a partir de agua a 80° C y 5 MPa.
Para generar el vapor se usa un gas combustible que consiste de 75% de CH 4; 10% de H2; 5% de CO2 y
10% de CO que se quema con aire atmosférico. El gas combustible y el aire entran a la caldera a 25° C y
1,2 bar. El aire entra con una humedad relativa de 80% y un exceso de 10% con respecto a la combustión
completa. Los gases de combustión salen a 460 K y contienen 2,5% (mol) de O2 en base seca. Determine:
a) el flujo molar de mezcla combustible; b) la relación molar aire/combustible; c) la composición molar de
los gases de combustión secos.
Base de cálculo: 1,0 kmol/h
yki ykj
Esp Qca nki nkj nk h°f298 h298 h460
(%) (%)
combustibl

CH4 (1) 75,0% 0,7500 0,0000 -0,750 -74850,0 0,0 0,0 0,000%
mezcla

H2 (2) 10,0% 0,1000 0,0000 -0,100 0,0 8468,0 13179,0 0,000%


e

CO2 (3) 5,0% 0,0500 0,8345 0,7845 -393520,0 9364,0 15916,0 10,823%
CO (4) 10,0% 0,1000 0,0655 -0,0345 -110530,0 8669,0 13412,0 0,850%
N2 (5) 79,0% 6,6176 6,6176 0 0,0 8669,0 13399,0 85,827%
Aire

O2 (6) 21,0% 1,7600 0,1928 -1,5672 0,0 8682,0 13525,0 2,500%


H2O (7) 0,0% 0,1809 1,7809 1,6000 -241820,0 9904,0 15428,0 0,000%
Fundamentos de Termodinámica

NCC = número de carbonos en el combustible; NHC = número de hidrógenos en el combustible;


NOC = número de oxígenos en el combustible; NNC = número de nitrógenos en el combustible;

Sea GS = kmol/h de gases secos

A 25° C, psat = 3,17 kPa


Humedad molar, wm = HR*psat/(P – HR*psat) = 0,80*3,17/(120 – 0,80*3,17) = 0,0216 kmol de vapor/kmol as

Combustión completa:

0,75CH4 +0,10H2 + 0,05CO2 + 0,10CO + atO2 --------> BCO2 + dH2O

Balance de Carbono: NCC = B; NCC = 0,75 ∗ 1 + 0,1 ∗ 0 + 0,05 ∗ 1 + 0,1 ∗ 1 = 0,90; 𝐁 = 𝟎, 𝟗𝟎

Balance de Hidrogeno: NHC = 2d; NHC = 0,75 ∗ 4 + 0,1 ∗ 2 + 0,05 ∗ 0 + 0,1 ∗ 0 = 3,20 = 2d; 𝐝 = 𝟏, 𝟔𝟎

Balance de Oxígeno: NOC + 2at = 2B + d; NOC = 0,75 ∗ 0 + 0,1 ∗ 0 + 0,05 ∗ 2 + 0,1 ∗ 1 = 0,20

𝟐𝐁 + 𝐝 − 𝐍𝐎𝐂 𝐤𝐦𝐨𝐥
𝐚𝐭 = = 𝟏, 𝟔𝟎; 𝐚 = 𝐚𝐭 (𝟏 + ) = 𝟏, 𝟔 ∗ (𝟏 + 𝟎, 𝟏) = 𝟏, 𝟕𝟔 𝐝𝐞 𝐎𝟐 𝐞𝐧 𝐞𝐥 𝐚𝐢𝐫𝐞 𝐬𝐞𝐜𝐨
𝟐 𝐡

Combustión incompleta:

0,75CH4 + 0,10H2 + 0,05CO2 + 0,10CO + aO2 + 3,76aN2 --------> bCO2 + cCO + dH2O + eO2 + fN2

Balance de Carbono: NCC = b + c = 0,90; b + c = B = (yCO2 + yCO)GS; (yCO2 + yCO) = B/GS

Balance de Oxígeno: NOC + 2a = 2b + c + d + 2e; e = yO 2*G; yO2 = e/GS

Balance de Nitrógeno: NNC + 2*(3,76a) = 2f; NNC = 0

𝐟 = 𝟑, 𝟕𝟔𝐚 = 𝟔, 𝟔𝟏𝟕𝟔 𝐤𝐦𝐨𝐥/𝐡 𝐝𝐞 𝐍𝟐 𝐞𝐧 𝐞𝐥 𝐚𝐢𝐫𝐞; 𝐲𝐍𝟐 = 𝐟/𝐆𝐒

yCO2 + yCO + yN2 + yO2 = 1; (1 − yO2 ) = (B + f)⁄G S

(𝐁 + 𝐟) (𝟎, 𝟗𝟎 + 𝟔, 𝟔𝟏𝟕𝟔) 𝐤𝐦𝐨𝐥


𝐆𝐒 = = = 𝟕, 𝟕𝟏𝟎𝟒 𝐝𝐞 𝐠𝐚𝐬𝐞𝐬 𝐝𝐞 𝐜𝐨𝐦𝐛𝐮𝐬𝐭𝐢ó𝐧 𝐬𝐞𝐜𝐨𝐬
(𝟏 − 𝐲𝐎𝟐 ) (𝟏 − 𝟎, 𝟎𝟐𝟓) 𝐡

Con el valor de GS:

yN2 = 6,6176/7,7104 = 0,85827;

e = 0,025*7,7104 = 0,1928 kmol/h de O2 en los gases secos;

Del balance de carbono: b + c = 0,90


Del balance de oxígeno: 2b + c = NOC + 2a – d – 2e = 0,20 + 2*1,76 – 1,60 – 2*0,1928 = 1,7345
2b + c = 1,7345;
Resolviendo simultáneamente las ecuaciones anteriores: b = 1,7345 – 0,900 = 0,8345; c = 0,0655 kmol/h

yCO2 = b/G = 0,8345/7,7104 = 0,10823; yCO = 0,0085

Flujo molar de aire seco = 4,76a = 8,378 kmol/h

Flujo de vapor con el aire = 4,76a*wm = 8,378*0,0216 = 0,1809 kmol/h


Fundamentos de Termodinámica

Flujo de vapor en los productos: 0,1809 + d = 0,1809 + 1,60 = 1,7809 kmol/h

Los flujos de cada especie química en los reactivos y productos y la composición molar de los gases secos,
resultado del balance de materia, aparecen en la tabla de arriba. Estos resultados están basados en una
base de cálculo de 1 kmol/h de combustible.

Relación molar aire seco/gas combustible = 8,378

Balance de energía: Q̇ = Ḣj − Ḣi

Las entalpías de cada especie química a 298 y 460 K aparecen en la tabla de arriba.
kJ
Ḣj = ∑ nkj [∆h°f298 + (h460 − h298 )] = − 7,184 ∗ 105
h
k

kJ
Ḣi = ∑ nki [∆h°f298 + (hTR − h298 )] = − 1,306 ∗ 105
h
k

𝐤𝐉 𝐤𝐦𝐨𝐥
𝐐̇ = 𝐇̇𝐏 − 𝐇̇𝐑 = −𝟓, 𝟖𝟕𝟖 ∗ 𝟏𝟎𝟓 𝐩𝐨𝐫 𝟏, 𝟎 𝐝𝐞 𝐜𝐨𝐦𝐛𝐮𝐬𝐭𝐢𝐛𝐥𝐞
𝐡 𝐡

Para calcular el flujo molar de gas combustible se debe hacer un balance de energía en la caldera.

A 80° y 5,0 MPa, hagua = 338,955 kJ/kg


A 500° y 5,0 MPa, hvapor = 3434,66 kJ/kg

𝐤𝐉
𝐐̇ = 𝐦̇𝐯𝐚𝐩 (𝐡𝐯𝐚𝐩 − 𝐡𝐚𝐠𝐮𝐚 ) = 𝟑𝟔𝟎𝟎 ∗ (𝟑𝟒𝟑𝟒, 𝟔𝟔 − 𝟑𝟑𝟖, 𝟗𝟓𝟓) = 𝟏, 𝟏𝟏𝟒 ∗ 𝟏𝟎𝟕
𝐡

Para generar 36000 kg/h de vapor se requiere un flujo de gas combustible de:
𝐤𝐦𝐨𝐥
𝐧̇ 𝐜𝐨𝐦𝐛 = 𝟏, 𝟏𝟏𝟒 ∗ 𝟏𝟎𝟕 ⁄𝟓, 𝟖𝟕𝟖 ∗ 𝟏𝟎𝟓 = 𝟏𝟖, 𝟗𝟔
𝐡

Flujo molar de gases de combustión secos:


𝐤𝐦𝐨𝐥
𝐧̇ 𝐠𝐚𝐬𝐞𝐬 𝐬𝐞𝐜𝐨𝐬 = 𝟏𝟖, 𝟗𝟔 ∗ 𝟕, 𝟕𝟏𝟎𝟒 = 𝟏𝟒𝟔, 𝟏𝟖𝟏
𝐡

También podría gustarte