10.
9 Para el sistema acetona (1) /acetonitrilo (2), la presión de vapor de los
componentes puros está dada por las siguientes ecuaciones de Antoine:
2940.46
Acetona=ln P 1sat =14.5463− (1)
t +237.22
2945.47
Acetonitrilo=ln P2sat =14.2724− (2)
t +224
Donde T está en °C y la presión de vapor en kPa. Considerando que la ley de
Raoult describe los estados de equilibrio vapor/líquido para este sistema,
determine:
a. x1 y y1 para las fases en equilibrio a 54°C y 65 kPa.
Diagrama x1 vs y1 para el sistema acetona(1)/ acetonitrilo (2) a 65 KPa y 54°C
1
0.9
0.8
0.7
0.6
0.5
y1
0.4
0.3
0.2
0.1
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
x1
Solución:
Datos:
2940.46
ln P1sat =14.5463−
54+237.22
Ecuación 1:
+B
( )
P1sat =exp T +C
2940.46
(14.5463− )
sat 54 +237.22
P1 =exp
P1sat =85.5636
Ecuación 2:
2945.47
ln P2sat =14.2724−
54+ 224
Ecuación:
2945.47
(4.2724− )
P2sat =exp 54 +224
P2sat =39.5349
Igualamos y resolvemos:
y 1 P+ y 2 P=x 1 P 1sat + x 2 P2sat
P( y ¿ ¿ 1+ y 2 )=x 1 P1sat +(1−x) P2sat ¿
P=x 1 P 1sat + P 2sat −x 1 P2sat
P=P2 sat + x 1 ( P¿¿ 1¿¿ sat−P 2sat ) ¿ ¿
Obtenemos X1
P−P 2sat
x 1= sat
P 1 −P2sat
65−39.53
x 1=
85.5636−39.5349
x 1=0.5533
Despejamos y obtenemos X2
x 2=1−x 1
x 2=1−0.5533
x 2=0.4467
Igualamos las ecuaciones para obtener Y1
y 1 P ¿ x 1 P1sat
x 1 P1sat
y 1=
P
0.5533∗85.5636
y 1=
65
y 1=0.7283
Despejamos y obtenemos Y2
y 2=1− y 1
y 2=1−0.7283
y 2=0.2717
b. T y y1 para P = 65 kPa y x1 = 0.4
Diagrama T vs x1, y1 para el sistema acetona(1)/ acetonitrilo (2) a 65 KPa.
70
65
60
55
T(°C)
50
45
40
35
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
x1, y1
Solución:
sat 2940.46
De la ecuación de la acetona: ln P1 =14.5463−
t+237.22
Despejamos T1sat
2940.46
sat
=14.5463−ln P 1sat
T 1 +237.22
2940.46
sat
=T 1sat + 237.22
14.5463−ln P1
2940.46
T 1sat = −237.22
14.5463−ln P1sat
2940.46
T 1sat = −237.22
14.5463−ln 65
T 1sat =46.2821 ° C
Calculamos T2 saturada.
B
T 2sat = −C
ln P− A
−2945.47
T 2sat = −224
ln (65)−14.2724
T 2sat =67.6981
Calculamos el valor de To:
T 0=x 1 T 1sat + x 2 T 2sat
T 0=0.4∗46.2821+ 0.6∗67.6981
T 0=59.1317 ° C
Calculamos el valor de la volatilidad:
B1 B2
ln α =0.2739+ −
T 0+ C1 T 0 +C 2
2940.46 2945.47
ln α =0.2739− +
59.1317+ 237.22 59.1317+224
α =2.1273
Calculamos P2 saturación:
P
P2sat =
x 1∗α + x 2
65
P2sat =
0.4∗2.1273+0.6
P2sat =44.7991
Calculamos T2 Saturación:
sat B2
T2 = −C 2
A 2−ln P2sat
2945.47
T 2sat = −224
14.2724−ln( 44.7991)
T 2sat =57.3190 ° C
Con la siguiente ecuación:
ln α =ln P1sat −ln P2sat
P 1sat
ln α =ln
( ) P 2sat
P1 sat
α=
( )
P2 sat
Despejamos para tener P1 saturación:
P1sat =α∗P 2sat
P1sat =2.1273∗44.7991
P1sat =95.3011
Con la siguiente ecuación:
y 1 P ¿ x 1 P1sat
Calculamos y1:
x 1 P1sat
y 1=
P
0.4∗95.3011
y 1=
65
y 1=0.5864
Calculamos y2:
y 1 + y 2=1
y 2=1− y 1
y 2=1−0.5864
y 2=0.4136
c. P y y1 para T = 54°C y x1 = 0.4
Solución:
Diagrama P vs x1, y1. Sistema acetona (1)/ Acetronitrilo (2), a 54 °C
90
80
70
60
Presión (KPa)
50
40
30
20
10
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
x1, y1
Con la gráfica obtenemos el valor de x1=0.4 y x2= 0.6
P1sat =85.5636
P2sat =39.5349
Calculamos P:
P=x 1 P 1sat + x2 P2sat
P=(0.4∗85.5636)+(0.6∗39.5349)
P=57.9463
Utilizamos la siguiente ecuación:
y 1 P=x 1 P1sat
Despejamos y1:
x 1 P1sat
y 1=
P
0.4∗85.5636
y 1=
57.9463
y 1=0.5906
De la siguiente ecuación:
y 1 + y 2=1
Despejamos y2
y 2=1− y 1
y 2=1−0.5906
y 2=0.4094
d. T y x1 para P = 65 kPa y y1 = 0.4
Diagrama T vs x1, y1 para el sistema acetona(1)/ acetonitrilo (2) a 65 KPa.
70
65
60
55
T(°C)
50
45
40
35
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
x1, y1
Solución:
sat 2940.46
De la ecuación de la acetona: ln P1 =14.5463−
t+237.22
Despejamos T1sat
2940.46
sat
=14.5463−ln P 1sat
T 1 +237.22
2940.46
sat
=T 1sat + 237.22
14.5463−ln P1
2940.46
T 1sat = −237.22
14.5463−ln P1sat
2940.46
T 1sat = −237.22
14.5463−ln(65)
T 1sat =46.2821 ° C
Calculamos T2 saturada.
B
T 2sat = −C
ln P− A
−2945.47
T 2sat = −224
ln (65)−14.2724
T 2sat =67.6981
Como:
y 1=0.4
y 2=0.6
Calculamos el valor de To:
T 0=x 1 P1sat + x 2 P2 sat
T 0=0.4∗46.2821+ 0.6∗67.6981
T 0=59.1317 ° C
Calculamos el valor de la volatilidad:
B1 B2
ln α =0.2739+ −
T 0+ C1 T 0 +C 2
2940.46 2945.47
ln α =0.2739− +
59.1317+ 237.22 59.1317+224
α =2.1273
Calculamos P1 saturación:
P1sat =P( y 1+ y 2∗α )
P1sat =65∗(0.4+ 0.6∗2.1273)
P1sat =108.9647
Calculamos la temperatura con los datos anteriores:
B1
T= −C 1
A 2−ln P1sat
2940.46
T= −237.22
14.2724−ln(108.9647)
T =69.6732 ° C
Con la siguiente ecuación:
y 1 P=x 1 P1sat
Calculamos x1:
y1 P
x 1=
y 1 P1sat
0.4∗65
x 1=
108.9647
x 1=0.2386
Calculamos x2:
x 1+ x2=1
x 2=1−x 1
x 2=1−0.2386
x 2=0.7614
e. P y x1 para T = 54°C y y1 = 0.4
Solución:
Diagrama P vs x1, y1. Sistema acetona (1)/ Acetronitrilo (2), a 54 °C
90
80
70
60
Presión (KPa)
50
40
30
20
10
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
x1, y1
Sabemos que y1=0.4 entonces y2=0.6
P1sat =85.5636
P2sat =39.5349
Tenemos las siguientes ecuaciones:
y1 P
x 1= sat
P1
y2 P
x 2= sat
P2
Como:
x 1+ x2=1
Entonces:
y1 P y2 P
+ =1
P 1sat P2sat
y1 y2
P∗
( P1 sat
+ sat
P2 )
=1
Despejamos P:
1
P=
y1 y
(
P1 sat
+ 2sat
P2 )
1
P=
0.4 0.6
( 85.5636 +
39.5349 )
P=50.3744
Utilizando la siguiente ecuación obtenemos x1:
y1 P
x 1= sat
P1
0.4∗50.3744
x 1=
85.5636
x 1=0.2354
f. La fracción del sistema que es líquida, x1 y y1 a 54°C y 65kPa, cuando
la composición promedio del sistema es de un 70% molar de acetona.
Solución:
z1=0.70 por lo tanto la z2= 0.30
P1sat =85.5636
P2sat =39.5349
Ubicación el punto de burbuja:
z 1=x 1
Por tanto:
z 1=x 1=0.70
z 2=x 2=0.30
P burb=x1 P1sat + x 2 P2sat
P burb=(0.70∗85.5636)+( 0.30∗39.5349)
P burb=71.7549
Ubicación en el punto de rocío:
z 1= y 1
Por tanto:
z 1= y 1=0.70
z 2= y 2=0.30
1
P de rocío=
y1 y2
P1(sat
+ sat
P2 )
1
P de rocío=
0.7 0.3
( +
85.5636 39.5349 )
P de rocío=63.4144
Por tanto, deducimos que:
P de rocío< 65 KPa< p burbuja
Calculamos K1:
P1sat
k 1=
P
85.5636
k 1=
65
k 1=1.3163
Calculamos K2
P2 sat
k 2=
P
39.5349
k 2=
65
k 2=0.6082
Igualando:
z 1∗k 1 z 2∗k 2
(1+v ( k 1−1)
+
)(
1+v ( k 2−1)
=1
)
( 1+v0.7∗1.3163
(1.3163−1) ) +
( 0.3∗0.6082
1+ v (0.6082−1) )
=1
v=0.8382
L=0.1618
Calculamos y1:
z1∗k 1
y 1=
(
1+ v (k 1 −1) )
0.7∗1.3163
y 1=
(
1+ 0.8382∗(1.3163−1) )
y 1=0.7283
Calculamos x1:
y1
x 1=
k1
0.7283
x 1=
1.3163
x 1=0.5532
Por lo que:
x 1+ x2=1
x 2=1−x 1
x 2=1−0.5532
x 2=0.4468
Se utilizó el método numérico de newton rhapson, programado en Matlab
Código principal
clc
clear all
%variables
z1=0.7;
z2=0.3;
k1=1.3163;
k2=0.6082;
%variable supuesta
vsup=0.6;
itermax=10;
erreq=1e-5;
for iter=1:itermax
[j,f,vcal]=solucion_m(z1,k1,z2,k2,vsup);
[ERR1]=error_1(vsup,vcal);
if norm(ERR1) < erreq break
else
ERR1=vsup-vcal;
vsup=vcal;
end
end
iter
ERR1
vcal
código para las ecuaciones
function[j,f,vcal]=solucion_m(z1,k1,z2,k2,vsup)
j= (z1*k1/(1+vsup*(k1-1)))+(z2*k2/(1+vsup*(k2-1)))-1; %funcion objetivo
f=(-z1*k1*(k1-1)/(1+vsup*(k1-1))^2)-(z2*k2*(k2-1)/(1+vsup*(k2-1))^2); %derivada
de la funcion obejetivo
vcal=vsup-(j/f); %newton rhapson
captura del valor obtenido
g. La fracción del sistema que es líquida, x1 y y1 a 54°C y 65kPa, cuando
la composición promedio del sistema es de un 60% molar de acetona.
Solución:
z1=0.60 por lo tanto la z2= 0.40
P1sat =85.5636
P2sat =39.5349
Ubicación el punto de burbuja:
z 1=x 1
Por tanto:
z 1=x 1=0.60
z 2=x 2=0.40
P burb=x1 P1sat + x 2 P2sat
P burb=(0.60∗85.5636)+( 0.40∗39.5349)
P burb=67.1521
Ubicación en el punto de rocío:
z 1= y 1
Por tanto:
z 1= y 1=0.60
z 2= y 2=0.40
1
P de rocío=
y1 y2
P1(sat
P2 )
+ sat
1
P de rocío=
0.6 0.4
( +
85.5636 39.5349 )
P de rocío=58.3772
Por tanto, deducimos que:
P de rocío< 65 KPa< p burbuja
Calculamos K1:
P1sat
k 1=
P
85.5636
k 1=
65
k 1=1.3163
Calculamos K2
P2 sat
k 2=
P
39.5349
k 2=
65
k 2=0.6082
Igualando:
z 1∗k 1 z 2∗k 2
(1+v ( k 1−1)
+
)(
1+v ( k 2−1) )
=1
( 1+v0.6∗1.3163
(1.3163−1) ) +
( 0.4∗0.6082
1+ v (0.6082−1) )
=1
v=0.2668
L=0.7332
Calculamos y1:
z1∗k 1
y 1=
(
1+ v (k 1 −1) )
0.6∗1.3163
y 1=
( 1+ 0.2668∗(1.3163−1) )
y 1=0.7283
Calculamos x1:
y
x 1= 1
k1
0.7283
x 1=
1.3163
x 1=0.5532
Por lo que:
x 1+ x2=1
x 2=1−0.5532
x 2=0.4468
Se utilizó el método numérico de newton rhapson, programado en Matlab
Código principal
clc
clear all
%variables
z1=0.6;
z2=0.4;
k1=1.3163;
k2=0.6082;
%variable supuesta
vsup=0.6;
itermax=10;
erreq=1e-5;
for iter=1:itermax
[j,f,vcal]=solucion_m(z1,k1,z2,k2,vsup);
[ERR1]=error_1(vsup,vcal);
if norm(ERR1) < erreq break
else
ERR1=vsup-vcal;
vsup=vcal;
end
end
iter
ERR1
vcal
código para las ecuaciones
function[j,f,vcal]=solucion_m(z1,k1,z2,k2,vsup)
j= (z1*k1/(1+vsup*(k1-1)))+(z2*k2/(1+vsup*(k2-1)))-1; %funcion objetivo
f=(-z1*k1*(k1-1)/(1+vsup*(k1-1))^2)-(z2*k2*(k2-1)/(1+vsup*(k2-1))^2); %derivada
de la funcion obejetivo
vcal=vsup-(j/f); %newton rhapson
10.10 Para el sistema n-pentano (1) / n-heptano (2), la presión vapor de los
componentes puros está dada por.
2477.07
LnP 1sat =13.8183−
t +233.21
2911.32
LnP 2sat =13.8587−
t +216.64
Donde t está en °C y la presión vapor en kpa. Considerando que la ley de Raoult
describe los estados de equilibrio vapor/ liquido de este sistema, determine:
a) X1, y y1 para las fases en equilibrio a 63°c y 95 kpa
Presión constante 95kpa
Solución:
De la ecuación de n-pentano
Datos:
Escriba aquí laecuación .Ecuación 1:
+B
( )
P1sat =exp T +C
2477.07
(13.8183− )
P1sat =exp 63+ 233.21
P1sat =234.10 kpa
Ecuación 2:
2911.32
LnP 2sat =13.8587−
t +216.64
Ecuación:
2911.32
(13.8587− )
sat 63 +216.64
P2 =exp
P2sat =31.42 kpa
Igualamos y resolvemos:
y 1 P+ y 2 P=x 1 P 1sat + x 2 P2sat
P( y ¿ ¿ 1+ y 2 )=x 1 P1sat +(1−x) P2sat ¿
P=x 1 P 1sat + P 2sat −x 1 P2sat
P=P2 sat + x 1 ( P¿¿ 1¿¿ sat−P 2sat ) ¿ ¿
Obtenemos X1
P−P 2sat
x 1= sat sat
P 1 −P2
95−31.42
x 1=
234.10−31.34
x 1=0.31
Despejamos y obtenemos X2
x 2=1−x 1
x 2=1−0.31
x 2=0.69
Igualamos las ecuaciones para obtener Y1
y 1 P ¿ x 1 P1sat
x 1 P1sat
y 1=
P
0.31∗234.10
y 1=
95
y 1=0.7639
Despejamos y obtenemos Y2
y 2=1− y 1
y 2=1−0.7639
y 2=0.2361
Diagrama de equilibrio Y1 vs X1 para el sistema n-pentano
(1)/ n-heptano (2) a 95kpa
1
0.9
0.8
0.7
0.6
0.5
Y1
0.4
0.3
0.2
0.1
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
X1
B) t y Y1 para P=95 kpa y x1=0.34
Solución de la ecuación de n-pentano
2477.07
LnP 1sat =13.8183−
t +233.21
2477.07
=13.8183−lnP 1 sat
T 1 sat +233.21
2477.07
T 1 sat= −233.21
13.8183−lnP 1 sat
T 1 sat=34.16° C
Calculamos T2 saturada.
B
T 2sat = −C
ln P− A
−2911.32
T 2sat = −216.64
ln (65)−13.8587
T 2sat =83.98
Calculamos el valor de To:
T 0=x 1 T 1sat + x 2 T 2sat
T 0=0.34∗34.16+ 0.66∗83.94
T 0=67.01 °C
Calculamos el valor de la volatilidad:
B1 B2
ln α =0.2739+ −
T 0+ C1 T 0 +C 2
2477.07 2911.32
ln α =0.2739− +
67.01+233.26 67.01+216.64
α =9.8571
Calculamos P2 saturación:
P
P2sat =
x 1∗α + x 2
95
P2sat =
0.34∗9.8571+0.66
P2sat =23.68 kpa
Calculamos T2 Saturación:
B2
T 2sat = −C 2
A 2−ln P2sat
2911.32
T 2sat = −216.64
13.8586−ln(23.66)
T 2sat =55.57 ° C
Con la siguiente ecuación:
ln α =ln P1sat −ln P2sat
P 1sat
( )
ln α =ln
P 2sat
P1 sat
α=
( )
P2 sat
Despejamos para tener P1 saturación:
P1sat =α∗P 2sat
P1sat =9.8571∗23.68
P1sat =233.41 kpa
Con la siguiente ecuación:
y 1 P ¿ x 1 P1sat
Calculamos y1:
x 1 P1sat
y 1=
P
0.34∗233.41
y 1=
95
y 1=0.83
Calculamos y2:
y 1 + y 2=1
y 2=1− y 1
y 2=1−0.83
y 2=0.17
Diagrama T vs x1, y1 para el sistema n-pentano( 1) / n-
heptano (2) a 95 kpa
100
90
80
70
60
50
40
30
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Diagrama P vs x1, y1 para el sistema n-pentano( 1) / n-heptano (2) a 60°c
300
280
260
240
220
200
180
160
P
140
120
100
80
60
40
20
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
x1,y1
C) P y Y1 para t=60°c y x1=0.44
Con la gráfica obtenemos el valor de x1=0.44 y x2= 0.56
P1sat =214.90 kpa
P2sat =28.0739 kpa
Calculamos P:
P=x 1 P 1sat + x2 P2sat
P=(0.44∗214.90)+(0.56∗28.0739)
P=110.27
Utilizamos la siguiente ecuación:
y 1 P=x 1 P1sat
Despejamos y1:
x 1 P1sat
y 1=
P
0.44∗214.90
y 1=
110.27
y 1=0.85 74
De la siguiente ecuación:
y 1 + y 2=1
Despejamos y2
y 2=1− y 1
y 2=1−0.8574
y 2=0.1426
d) t y x1 para P 85 kpa y y1=0.86
Diagrama T vs x1, y1 para el sistema n-pentano( 1) / n-heptano (2) a 85 kpa
100
90
80
70
60
50
40
30
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Solución:
sat 2477.07
De la ecuación de la n-pentano: LnP 1 =13.8183−
t +233.21
Despejamos T1sat
2477.07
sat
=13.8183−ln P 1sat
T 1 +233.21
2477.07
sat
=T 1sat + 233.21
13.8183−ln P1
2477.07
T 1sat = −233.21
13.8183−ln P1sat
2477.07
T 1sat = −233.21
13.8183−ln ( 85 )
T 1sat =30.99 ° C
Calculamos T2 saturada.
B
T 2sat = −C
ln P− A
−2911.32
T 2sat = −216.64
ln (85)−13.8587
T 2sat =92.5470 ° C
Como:
y 1=0.86
y 2=0.14
Calculamos el valor de To:
T 0= yT 1sat + y 2 T 2sat
T 0=0.86∗30.99+0.14∗92.5470
T 0=39.6079
Calculamos el valor de la volatilidad:
B1 B2
ln α =0.2739+ −
T 0+ C1 T 0 +C 2
2477.07 2911.32
ln α =0.2739− +
39.6070+ 216.64 39.6070+216.64
α =7.1603
Calculamos P1 saturación:
P1sat =P( y 1+ y 2∗α )
P1sat =85∗(0.86+0.14∗7.1603)
P1sat =158.3075 kpa
Calculamos la temperatura con los datos anteriores:
B1
T= −C 1
A 2−ln P1sat
2477.07
T= −233.21
13.8587−ln(158.3075)
T =48.46 ° C
Con la siguiente ecuación:
y 1 P=x 1 P1sat
Calculamos x1:
y1 P
x 1=
P 1sat
0.86∗85
x 1=
158.3075
x 1=0.4617
Calculamos x2:
x 1+ x2=1
x 2=1−x 1
x 2=1−0.4617
x 2=0.5383
e) P y x1 para t=70°C y y1=0.08
Diagrama P vs x1, y1 para el sistema n-pentano( 1) / n-heptano (2) a 70°c
300
280
260
240
220
200
180
160
P
140
120
100
80
60
40
20
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
x1,y1
Sabemos que y1=0.08 entonces y2=0.92
P1sat =283.9543
P2sat =40.5277
Tenemos las siguientes ecuaciones:
y P
x 1= 1 sat
P1
y P
x 2= 2 sat
P2
Como:
x 1+ x2=1
Entonces:
y1 P y2 P
+ =1
P 1sat P2sat
y1 y2
P∗
( P1 sat
+ sat
P2 )=1
Despejamos P:
1
P=
y1 y
(
P1 sat
+ 2sat
P2 )
1
P=
0.08 0.92
( 283.9543 +
40.5277 )
P=43.5118
Utilizando la siguiente ecuación obtenemos x1:
y1 P
x 1= sat
P1
0.08∗43.5118
x 1=
283.9543
x 1=0.01
f) La fracción del sistema es líquida x1 y y1 a 60 °c y 115 kpa, cuando la
composición promedio del sistema es equimolar
z1=0.50 por lo tanto la z2= 0.50
P1sat =214.9039
P2sat =28.0739
Ubicación el punto de burbuja:
z 1=x 1
Por tanto:
z 1=x 1=0.50
z 2=x 2=0.50
P burb=x1 P1sat + x 2 P2sat
P burb=(0.50∗214.9083)+(0.50∗28.0739)
P burb=121.4911
Ubicación en el punto de rocío:
z 1= y 1
Por tanto:
z 1= y 1=0.50
z 2= y 2=0.50
1
P de rocío=
y1 y2
( P1 sat
+ sat
P2 )
1
P de rocío=
0.5 0.5
( 214.9039 +
28.0739 )
P de rocío=49.6604
Por tanto, deducimos que:
P de rocío< 115 KPa< p burbuja
Calculamos K1:
P1sat
k 1=
P
214.9039
k 1=
115
k 1=1.8687
Calculamos K2
P2 sat
k 2=
P
28.0739
k 2=
115
k 2=0.2441
x 2=1−x 1
x 2=1−0.5532
x 2=0.4468
g) La fracción del sistema que es líquida, x1 y y1 a 60°C y 115 cuando la
composición del sistema es de un 60% molar de n-pentano
z1=0.60 por lo tanto la z2= 0.40
P1sat =214.9039
P2sat =28.0739
Ubicación el punto de burbuja:
z 1=x 1
Por tanto:
z 1=x 1=0.50
z 2=x 2=0.50
P burb=x1 P1sat + x 2 P2sat
P burb=(0.60∗214.9083)+(0.40∗28.0739)
P burb=140.1745
Ubicación en el punto de rocío:
z 1= y 1
Por tanto:
z 1= y 1=0.60
z 2= y 2=0.40
1
P de rocío=
y1 y2
P1(sat
+ sat
P2 )
1
P de rocío=
0.6 0.4
( +
214.9039 28.0739 )
P de rocío=58.6852
Por tanto, deducimos que:
P de rocío< 115 KPa< p burbuja
Calculamos K1:
P1sat
k 1=
P
214.9039
k 1=
115
k 1=1.8687
Calculamos K2
P2 sat
k 2=
P
28.0739
k 2=
115
k 2=0.2441
Igualando:
z 1∗k 1 z 2∗k 2
(1+v ( k 1−1)
+
)(
1+v ( k 2−1)
=1
)
0.5∗1.8687 0.5∗0.2441
(1+v (1.8687−1)
+
)(
1+v (0.2441−1)
=1
)
v=0.3333
L=0.6667
Calculamos y1:
z1∗k 1
y 1=
( 1+ v (k 1 −1) )
0.6∗1.8687
y 1=
(
1+ 0.3333∗(1.8687−1) )
y 1=0.9680
Calculamos x1:
y1
x 1=
k1
0.9680
x 1=
1.8687
x 1=0.5180
Por lo que:
x 1+ x2=1
x 2=1−x 1
x 2=1−0.5532
x 2=0.4820
Se utilizó un método numérico newton rhapson, programado en Matlab
Código principal
%variables
z1=0.6;
z2=0.4;
k1=1.8687;
k2=0.2441;
%variable supuesta
vsup=0.6;
itermax=5000;
erreq=1e-5;
for iter=1:itermax
[j,f,vcal]=solucion_m(z1,k1,z2,k2,vsup);
[ERR1]=error_1(vsup,vcal);
if norm(ERR1) < erreq break
else
ERR1=vsup-vcal;
vsup=vcal;
end
end
iter
ERR1
vcal
Código de ecuaciones
function[j,f,vcal]=solucion_m(z1,k1,z2,k2,vsup)
j= (z1*k1/(1+vsup*(k1-1)))+(z2*k2/(1+vsup*(k2-1)))-1; %funcion objetivo
f=(-z1*k1*(k1-1)/(1+vsup*(k1-1))^2)-(z2*k2*(k2-1)/(1+vsup*(k2-1))^2); %derivada
de la funcion obejetivo
vcal=vsup-(j/f); %newton rhapson
captura de resultado
