Universidad nacional “José Faustino Sánchez Carrión”

EAP de Ingeniería Química

PROBLEMAS

1.11 ¿cuál es el cambio en la entropía cuando 0.7m3 de CO2 0.3m3 de N2 cada uno
a 1 bar y 25 oC se combinan para formar una mezcla de gas?
Solución:
Datos:
T= 25+273.15 = 298.15K
P=1 bar
Para un gas ideal fracción molar=fracción volumen
VCO2 =0.7m3 Y1 = 0.7
VN2 =0.3m3 Y2 = 0.3
Reemplazando en la ecuacion:

PV t 1 ( 1 ¯m 3 )
nt = =
RT 83.14∗10−6 m3 ¯¿
∗298.15 k =40.342 mol ¿
mol k

La fórmula para mezcla de gases ideales el cambio de entropía es:

S gi −∑ Yi S gi =−nR ∑ Yi ln(Yi )
i i

S gi −∑ Yi S gi
i

J
−8.314 ∗40.342 mol( 0.7 ln 0.7+0.3 ln 0.3)
mol k

J
∆ S=204.8859
K

1.12 Un recipiente dividido en dos partes por un separador contiene 4 moles de

gas N2 a 75OC Y 30bar en un lado y 2.5moles de gas argón a 130 OC y 20 bar en el
otro. Si el separador se elimina y los gases se mezclan en forma adiabática y
completa ¿Cuál es el cambio en la entropía? suponga que el N2 y Ar son gases
ideales con Cv=5/2R ; Cv=3/2R respectivamente.

Solución:
Datos:
N2 Ar
n=4mol n=2.5mol
 Universidad nacional “José Faustino Sánchez Carrión”
EAP de Ingeniería Química

Hallando la temperatura:

T=75+¿273.15=348.15k T=130+273.15=403.15k

P=30bar P=20 bar

Cv=5/2 R Cv=3/2 R
Cv=3/2 R Cv=3/2 R

Si: nt =4 mol+2.5 mol=6.5 mol

4 2.5
Y N 2= =0.615 ; Y Ar = =0.385
6.5 6.5

Hallando con la ecuacion del cambio de entropía:

∆ S=∆ S N 2 +∆ S Ar + ∆ S mezcla

Calculo de la Teq haciendo un balance de energía en una mezcla:

∆ E N 2 =−∆ E Ar

n N 2∗Cv N 2∗( T eq −T N 2 ) =n Ar∗Cv Ar∗( T Ar −T eq )

mol∗5 mol∗3
4 R∗( T eq −348.15 ) =2.5 R∗(403.15−T eq )
2 2

27.5 T eq=9986.625 K

T eq=363.15 K

Calculo de la presión según: PV=nRT

V T =V N 2+ V Ar

nt∗R∗T n N 2∗R∗T N 2 n Ar∗R∗T Ar

= +
P PN2 P Ar

6.5∗363.15 4∗348.15 2.5∗403.15

= +
P 30 20

P=24.38 ¯¿

Calculo de la entropía para gases ideales (Cp = cte)

 Universidad nacional “José Faustino Sánchez Carrión”
EAP de Ingeniería Química

T2 P2
∆ S=nCpln ( )
T1
−nRln ( )
P1

mol∗7 J 363.15 J 24.38

∆ S N 2=4
2
∗8.314
mol k
∗ln (
348.15 )
−4 mol∗8.314
mol k
∗ln
30( )
J
∆ S N 2=11.806
mol k

mol∗5 J 363.15 J 24.38

∆ S Ar =2.5
2
∗8.314
mol k
∗ln (
403.15 )
−2.5 mol∗8.314
mol k
∗ln
20( )
J
∆ S Ar =−9.547
mol k

∆ S mezcla=nT (−R ∑ Yi ln (Yi))

i

J
∆ S mezcla=−6.5 mol∗8.314 ¿
mol k

J
∆ S mezcla=35.99
mol k

Reemplazando en la ecuación (1) las entropías halladas tenemos:

J J J
∆ S=11.806 −9.547 +35.99
mol k mol k mol k

J
∆ S=38.27
mol k

1. Una corriente de N2 que circula con una relación de 2kg/s y otra de H2 que
circula con una relación de 0.5 kg/s se mezclan de manera adiabática en un
proceso de flujo uniforme .Si los gases se suponen ideales, ¿cuál es la rapidez
con que aumenta la entropía como resultado del proceso?
Solución:
Datos:
N2 H2
ṁ=2Kg/s ṁ=0.5Kg/s
PM=28.014Kg/kmol PM=2.016 kg/Kmol
Calculo del flujo molar del N2:
 Universidad nacional “José Faustino Sánchez Carrión”
EAP de Ingeniería Química

2 kg s−1 Kmol
ṅ= =0.0714
28.014 Kg kmol−1
S

Calculo del flujo molar del H2:

0.5 kg s−1 Kmol
ṅ= =0.248
2.016 Kg kmo l−1
s
Flujo molar total
Kmol
ṅt =n ˙N 2 + n ˙H 2=0.0714+ 0.248=0.3194
s
n ˙N 2 0.0714
Y 1= = =0.224
ṅt 0.3194
n ˙H 2 0.248
Y 2= = =0.776
ṅt 0.3194
∆ S=−R∗ṅt ∑ Yi ln (Yi)
i
kJ kmol
∆ S=−8.314 ∗0.3194 ∗(0.224∗ln ( 0.224 ) +0.776∗ln ⁡( 0.776))
kmol k s
kJ
∆ S=1.412
sK

2. ¿cuál es el trabajo ideal requerido para separar una mezcla equimolar de

metano y etano a 175oC y 3 bar en un proceso de flujo uniforme en dos
corrientes de gases puros a 35oC y 1 bar si T σ =300 K ?

SOLUCION:

Para el cálculo de trabajo ideal se aplica la ecuación siguiente:

W ideal =∆ H −T σ ∆ S … .(1)

Tabla de capacidades caloríficas de los componentes:

moles A B C D
metano 0,5 1,702 9,081*10-3 -2,164*10-6 0
etano 0,5 1,31 19,225*10-3 -5,561*10-6 0

Para gases ideales:

T
∆ Cp ∆ B 2( 2 ) ∆ C 3( 3 )
∆ H =R ∫ dT =( ∆ A ) T 0 ( τ−1 ) + T τ −1 + T τ −1 … .(2)
T0 R 2 0 3 0

T
T0=448.15K ; T=308.15K ; τ = =0.6876
T0
 Universidad nacional “José Faustino Sánchez Carrión”
EAP de Ingeniería Química

∆ A=0.5∗1.702+ 0.5∗1.31=1.4165

∆ B=0.5∗9.081∗10−3 +0.5∗19.225∗10−3=14.153∗10−3

∆ C=0.5∗(−2.164∗10−6 ) +0.5∗(−5.561∗10−6 )=−3.8625∗10−6

Reemplazando los datos obtenidos en la ecuación (2):

T

∫ ∆ RCp dT =−198.3129−749.2829+78.2098=−869.386
T0

J J
∆ H =8.314 (−869.386 k )=−7228.075
mol k mol
T
∆S ∆Cp dT P
=∫ −ln + ∑ Yi ln (Yi ) … … .(3)
R T0 R T P0 i
T
⟨ ∆ Cp ⟩ lnT
∫ ∆ RCp dT
T
=
R T0
… … …..( 4)
T0

⟨ ∆Cp ⟩ ∆D ( τ+ 1 )
R ( (
=∆ A+ ∆ B T 0 + ∆ C T 0 +
2
T
2 2 0
τ T0 )( 2 ))( τ−1
lnτ )
… ..(¿)

Reemplazando en (*)

⟨ ∆Cp ⟩
=1.4165+ ( 6.3427−0.6545 )∗0.8341=6.161
R

∆S 308.15 1
R
=6.161∗ln
448.15( 3 )
−ln +2∗0.5∗ln 0.5=−1.902

J J
∆ S=8.314 (−1.902 )=−15.813
mol k mol k

Reemplazando en la ecuación (1):

W ideal =−7228.075−300∗(−15.813)

J
W ideal =−2484.1
mol

3. ¿Cuál es el trabajo requerido para separar aire (21% mol de oxigeno y 79%
mol de nitrógeno) a 25oC y 1 bar, si la eficiencia termodinámica del proceso es
5% y si T σ =300 K ?
 Universidad nacional “José Faustino Sánchez Carrión”
EAP de Ingeniería Química

Solución:
Base de cálculo: 1mol de aire
Y1=0.21 ; Y2=0.21; n= 0.05

W ideal
W= ….( 1) ; W ideal =∆ H −T σ ∆ S …(2)
ni
Para gases ideales: H gi =∑ Yi H gi
H gi−∑ Yi H gi =0 ; ∆ H =0
Para el cálculo de la entropía aplicamos la siguiente expresión:
S gi =∑ Yi S gi −R ∑ Yi lnYi
J J
∆ S=R ∑ Yi lnYi=−8.314 ( 0.21 ln 0.21+0.79 ln 0.79 )=−4.273
mol k molk

Reemplazando en (2) y (1):

J J
(
W ideal =−300 K −4.273
mol K )
=1281.9
mol

1281.9 J J
W= =25638
0.05 mol mol

4. A partir de la información siguiente para el factor de compresibilidad de CO 2 a

150oC elabore graficas de la fugacidad y del coeficiente de fugacidad del CO 2 en
función de P para presiones hasta 500 bar. Compare sus resultados con los que
se obtengan a partir de la correlación generalizada representada por la
ecuación (11.68).
P
[
∅=exp r ( Bo + ω B1 )
Tr ]
P (bar) Z
10 0,985
20 0,970
40 0,942
60 0,913
80 0,885
100 0,869
200 0,765
300 0,762
400 0,824
500 0,910

Solución:
 Universidad nacional “José Faustino Sánchez Carrión”
EAP de Ingeniería Química

Para una P=0 sabemos que:

ln ∅1=0 ∅ 1=1

Utilizamos la regla del trapecio para integrar la ecuación:

P
dP
ln ∅1=∫ (¿ Z i−1) ¿
0 P

Z i−1
F i= Para i=2, 3,4…
Pi

F i−F i−1
Ai = . ( Pi−P i−1 )
2

ln ∅i=ln ∅ i−1 + Ai

Pi Zi Fi Fi-1 Ai ln∅i-1 ln∅i ∅i fi (bar)

(bar)
0 0.000 0.00000 0.00000 0.00000 0.00000 0.00000 1.000 0.000
10 0.985 -0.00150 0.00000 -0.00750 0.00000 -0.00750 0.993 9.925
20 0.970 -0.00150 -0.00150 -0.01500 -0.00750 -0.02250 0.978 19.555
40 0.942 -0.00145 -0.00150 -0.02950 -0.02250 -0.05200 0.949 37.973
60 0.913 -0.00145 -0.00145 -0.02900 -0.05200 -0.08100 0.922 55.332
80 0.885 -0.00144 -0.00145 -0.02888 -0.08100 -0.109875 0.896 71.676
100 0.869 -0.00131 -0.00144 -0.02748 -0.10988 -0.137350 0.872 87.167
200 0.765 -0.00118 -0.00131 -0.12425 -0.13735 -0.261600 0.770 153.96
4
300 0.762 -0.00079 -0.00118 -0.09842 -0.26160 -0.36001667 0.698 209.29
9
400 0.824 -0.00044 -0.00079 -0.06167 -0.36002 -0.42168333 0.656 262.37
7
500 0.910 -0.00018 -0.00044 -0.03100 -0.42168 -0.45268333 0.636 317.96
0

De la correlación generalizada representada por la ecuación:

Pr o
∅=exp
[ Tr
( B + ω B1 )
]
Para CO2 a 150ºC; T=423.15 K

Datos: Tc= 304.2 K Pc= 73.83bar ω=0.224

 Universidad nacional “José Faustino Sánchez Carrión”
EAP de Ingeniería Química

T 423.15
T r= = =1.3910
T c 304.2
P
Pr =
Pc

Expresaremos la ecuación anterior en función de la P:

P
P
[
∅ ( P)=exp c (B o +ω B1)
Tr ]
0.422
Bo=0.083−
T 1.6
r
0.422
Bo=0.083− =−0.16588
1.39101.6

0.172
B1=0.139−
T 4.2
r
0.172
B1=0.139− =0.09599
1.3910 4.2

Dato Nº Pi (bar) ∅i fi
1 0 1.000 0.000
2 10 0.986 9.860
3 20 0.972 19.445
4 40 0.945 37.813
5 60 0.919 55.146
6 80 0.894 71.490
7 100 0.869 86.885
8 200 0.755 150.980
9 300 0.656 196.768
10 400 0.570 227.949
11 500 0.495 247.566

Graficas:
 Universidad nacional “José Faustino Sánchez Carrión”
EAP de Ingeniería Química

Coeficiente de fugacidad vs P
1

Coeficiente de fugacidad
0.8 vs P
Logarithmic (Coeficiente
de fugacidad vs P)
Power (Coeficiente de
fugacidad vs P)
Correlacion

0.6

0.4
0 100 200 300 400 500 600

Fugacidad vs P
350.000

300.000

250.000
Fugacidad vs P
200.000 Correlacion
150.000

100.000

50.000

0.000
0 100 200 300 400 500 600
 Universidad nacional “José Faustino Sánchez Carrión”
EAP de Ingeniería Química

5. Para el SO2 a 600k y 300 bar establezca aproximaciones convenientes de la

fugacidad y de G R /RT .

Solución:

De la ecuación:

∅=(∅ ¿¿ 0) ¿¿

ω= 0.245 Tc=430.8 K Pc= 78.84 bar

Tr\Pr 3,0000 3,8052 5,0000

1,30 0,6383   0,5383
1,39   ∅0  
1,40 0,7145   0,6237

∅ 0=¿ 0,67659842

Tr\Pr 3,0000 3,8052 5,0000

1,30 1,2853   1,3868
1,39   ∅1  
1,40 1,2942   1,4888

∅1 1,3826429
= 1

∅=(∅ ¿¿ 0) ¿¿ 0,73249546

f =∅ P=( 0,73249546 ) ¿219,748637 bar.

GR
=ln ∅
RT

GR
=ln ( 0,73249546 )=¿-0,31129814
RT

6. calcule la fugacidad de isobutileno como un gas:

a) A 280oC y 20 bar ; b) A 280oC y 100 bar

Solución:

a) T=553.15 K P=20 bar

Datos: Tc= 417.9 K Pc= 40 bar ω=0.194
 Universidad nacional “José Faustino Sánchez Carrión”
EAP de Ingeniería Química

T 553.15
T r= = =1.3236
T c 417.9
P 20
Pr = = =0.5
Pc 40

Pr o
∅=exp
[ Tr
( B + ω B1 )
]
0.422
Bo=0.083−
T 1.6
r

o 0.422
B =0.083− =−0.18640
1.32361.6

0.172
B1=0.139−
T 4.2
r

1 0.172
B =0.139− =0.0860
1.3236 4.2

0.5
∅=exp [ 1.3236
(−0.18640+ 0.194 ×0.0860) ]
∅=0.9379

f =∅ P=( 0,9379 ) ¿18.76 bar.

b) T=553.15 K P=100 bar

Datos: Tc= 417.9 K Pc= 40 bar ω=0.194

T 553.15
T r= = =1.3236
T c 417.9
P 100
Pr = = =2.5
Pc 40

Pr o
∅=exp
[ Tr
( B + ω B1 )
]
0.422
Bo=0.083−
T 1.6
r
 Universidad nacional “José Faustino Sánchez Carrión”
EAP de Ingeniería Química

0.422
Bo=0.083− =−0.18640
1.32361.6

0.172
B1=0.139−
T 4.2
r
0.172
B1=0.139− =0.0860
1.3236 4.2

2.5
∅=exp [ 1.3236
(−0.18640+ 0.194 ×0.0860) ]
∅=0.725659

f =∅ P=( 0.725659 ) ¿72.5659 bar.

7. Calcule la fugacidad de uno de los siguientes:

a) ciclo pentano a 110oC y 275 bar. A 110 oC la presión de vapor del ciclo
pentano es 5.267 bar.
b) 1-buteno a 120 oc y 34 bar .A 120oC la presión de vapor del 1-buteno es
25.83 bar.

Solución:

Para el cálculo de la fugacidad a presiones altas empleamos la siguiente expresión:

V sat sat
i ( P−P i )
sat
f =∅ P exp
i
sat
i ( RT ) …(1)

a) Para el ciclopentano tenemos:

ω=0.196 ; T =110+ 273.15=383.15 K
3
T c =511.8 k ; Pc =45.02 ¯¿ ; Z c =0.273 ; V c =258 c m
mol
383.15 5.2675
T r= =0.7486 ; Psatr = =0.1169
511.8 45.02

Calculo del coeficiente de fugacidad de la presión de vapor empleando la ecuación:

P sat
∅ =exp[ r ¿(B 0+ ω B 1)]¿
sat
i
Tr

0.422 0.422
B0=0.083− 1.6
=0.083− =−0.5877
Tr 0.74861.6
 Universidad nacional “José Faustino Sánchez Carrión”
EAP de Ingeniería Química

0.172 0.172
B1=0.139− 4.2
=0.139− =−0.4413
Tr 0.74864.2

0.1169
∅ sat
i =exp [ 0.7486
(−0.5877+ 0.196∗0.4413 ) =0.900 ]
c m3 c m3
2
2 /7
sat (1−T r ) 7
0.2514
V =V C Z C
i =258 ∗0.273 =107.546
mol mol

Reemplazando en (1) los datos hallados tenemos:

f =0.9∗5.267 ¯¿ exp ¿¿

f =0.9∗5.267 ¯¿ exp ( 0.9106 )=11.783 ¯¿

b) Para 1-buteno tenemos:

ω=0.191 ; T =120+273.15=393.15 K
3
T c =420 k ; Pc =40.43 ¯¿ ; Z c =0.277 ; V c =239.3 c m
mol
393.15 25.83
T r= =0.9361 ; Psatr = =0.6389
420 40.43

Calculo del coeficiente de fugacidad con la presión de vapor empleando la ecuación:

sat P sat
r 0 1
∅ 2 =exp[ ¿(B + ω B )]¿
Tr

0.422 0.422
B0=0.083− 1.6
=0.083− =−0.3860
Tr 0.93611.6

0.172 0.172
B1=0.139− 4.2
=0.139− =−0.088
Tr 0.93614.2

0.6389
∅ sat
i =exp [ 0.9361
(−0.3860−0.191∗0.088 ) =0.7596 ]
c m3 c m3
2
2 /7
sat (1−T r ) 0.039 7

V =V C Z C
i =239.3 ∗0.277 =133.31
mol mol

Reemplazando en (1) los datos hallados tenemos:

f =0.7596∗25.83 ¯¿ exp ¿¿
 Universidad nacional “José Faustino Sánchez Carrión”
EAP de Ingeniería Química

f =0.7596∗25.83 ¯¿ exp ( 0.0333 )=20.284 ¯¿