Calculo Flash
Calculo Flash
Calculo Flash
Temperatura°
R 600
Presion(psia) 1000
2
Coeficientes de interacción binaria C1-C7 = 0.03812
Inicio Calculo de equilibrio de fases (bifásico)
Calculo de
Whitsol & Brulé, 2000
kEquilibrio Phase Behavior, Whitson, Curtis, H. &
Brulé, Michael R.
2. Encontrar los
valores kEquilibrio Monograph volumen 20 SPE Henry L.
Max y MIn Doherty Series
3. Resolver Rachford-Rice con
Fvmin y Fvmax
9. Actualiza kEquilibrio con
relación de fugacidades y/o
4. Calcular las composiciones xi GDEM
y yi (fase líquida y vapor)
8. Evaluar
condición de
5. Calcular los factores de
igualdad de
compresibilidad Zl y Zv de EOS
fugacidades. 𝑁
2
Converge? 10. 𝑙𝑛𝑘𝐸𝑞𝑢𝑖𝑙𝑖𝑏𝑟𝑖𝑜𝑖 < 1𝑥10−4
𝑖=1
3
CALCULO DE EQUILIBRIO DE FASE
PARA LA ECUACION DE ESTADO DE PENG-ROBINSON:
Z 3 (1 B)Z 2 ( A 2B 3B 2 )Z ( AB B 2 B 3 ) 0
para cada componente dentro de la mezcla:
fi Bi A Bi 2
Z B 1 2
l n l ni Z 1 l n Z B
son independientes
y j Ai j l n
cuando ysep estableceB el equilibrio n+1 variables
i
especificadas anterior al calculo.
2 2 B B A j
Z y Bdeben
1 2ser
𝑁 𝑁
2 1 m 1 Tr1/2
1
donde:
5
https://es.wikipedia.org/wiki/Ecuaci%C3%B3n_de_tercer_grado
Paso 1. Valor inicial para Ki y 2.
Identificar valor máx. y min.
Componente Pc, psia Tc, °R Pr Tr Ki Paso 1 identificar máx. y min.
C1 667.78 343.1 1.497 1.749 6.85878 Ki
C1 6.85878
C7+ 371.43 1082 2.692 0.555 0.00084
C7 0.00084
Ki
Min 0.00084
𝑒𝑥𝑝 5.37 1 + 𝜔𝑖 1 − 𝑇𝑟𝑖−1 Max 6.85878
𝐾𝑖 =
𝑝𝑟𝑖
Paso 2 calculo de Ki
Ki
C1 4.2793
C7 0.0045
Ki
Min 0.0045
Max 4.2793
6
3. Ecuación de Rachford-Rice, modificada por Muskat y
McDowell
Frac
C Ki 1
Mol 𝐹𝑉𝑚𝑖𝑛 = Paso 3.1 Ciclo cerrado para condición de equilibrio
1 − 𝐾𝑚𝑎𝑥
C1 0.4 Encontrar el valor Fvmin y Fvmax Valor medio
4.2793
1
C7+ 0.6 𝐹𝑉𝑚𝑎𝑥 =
0.0045 1 − 𝐾𝑚𝑖𝑛 Fvmin -0.30495 0.34981
Fvmax 1.00456
7
Fracción de vapor en la mezcla a P y T, Fv = 0.21886
4. Calculo de las fracciones molares
𝑧𝑖 𝑧𝑖 𝐾𝑖
𝑥𝑖 = 𝑦𝑖 = = 𝐾𝑖 𝑥𝑖
𝐹𝑉 𝐾𝑖 − 1 + 1 𝐹𝑉 𝐾𝑖 − 1 + 1
Component
Frac. Mol
e
Ki Xi Yi
C1 0.4 4.2793 0.2329 0.9965
C7+ 0.6 0.0045 0.7671 0.0035
8
Paso 5.0
• Calcular valores que no se modificaran a lo largo del
algoritmo
Frac.
Componente PM Pc, psia Tc, °R Tri Pri mi alfai Ai Bi
Mol
C1 0.4 16.043 667.78 343.1 0.013 1.749 1.497 0.395 0.76169034 0.17052088 0.067
C7+ 0.6 128 371.43 1082 0.413 0.555 2.692 0.966 1.55377693 6.21565304 0.378
2 1 m 1 Tr1/2
1
Términos constantes
𝑎 𝑇 𝑝 𝑝𝑟𝑖
𝐴𝑖 = 2 2
= 0.45724 2 𝛼 (AiAj)^0.5*(1-kij)
𝑅 𝑇 𝑇𝑟𝑖
i=1,j=1 0.17052
𝑏𝑝 𝑝𝑟𝑖 i=1,J=2 0.99027
𝐵𝑖 = = 0.0778
𝑅𝑇 𝑇𝑟𝑖 i=2,j=2 6.21565
𝑁 𝑁
i=1,J=2 0.99027
i=2,j=2 6.21565
𝑁
𝐵 = 𝑦𝑖 𝐵𝑖
𝑖=1
10
5. Calcular ZL y Zg
• P-R Z 3 (1 B)Z 2 ( A 2B 3B 2 )Z ( AB B 2 B 3 ) 0
Liq Gas Variable Líquido Gas
A 4.02088 0.1763 a1 -0.6948 -0.9323
B 0.30517 0.0677 a2 3.1311 0.0271
a3 -1.1055 -0.0070
Q= 0.9901 -0.0875
R= 0.2026 0.0293
D= 1.01155 0.000189
D >0 Una solución real, 2 complejas
S1= 1.0651 0.3505
S2= -0.9295 0.2497
xreal 0.3672 0.9110
Z 0.3672 0.9110
Considerando el
discriminante
se tiene
1.una de las raíces es real y
dos de ellas son complejas si
D>0
11
fi
6. Cálculo de las fugacidades
Bi A Bi 2 Z B 1 2
l n l ni Z 1 l n Z B y j Ai j l n
yi p B 2 2 B B A j
Z B 1 2 +
Liq Gas
xi yi fi Liq fi gas gL gV
916.806
0.2329 0.9965
C1 9 906.4146 1.58839333 6.78577995
C7 0.7671 0.0035 0.801961 0.87135 -0.16930119 -0.00047964
1.41909214 6.78530031
gmezcla 2.59351888
Liq Gas
xi yi fi Liq fi gas (fLi/fvi-1)^2
916.806
0.2329 0.9965
C1 9 906.4146 0.000131
0.80196
C7
0.7671 0.0035 1 0.87135 0.006341
0.006473
14
9. De no alcanzarse convergencia
actualizar Ki
𝒏
𝒏+1 𝒏 𝒇𝑳𝒊
𝑲𝒊 = 𝑲𝒊 𝒏
,
𝒇𝑽𝒊
15
10. Evaluar la convergencia para la solución trivial
(𝐾𝑖 → 1) con la condición:
𝑵
2
𝒍𝒏 𝑲𝒊 < 1𝒙10−4 .
𝒊=1
ki(n+1) (ln(Ki))^2
C1 4.3283 2.1467
C7 0.0042 30.0060
32.1527
16
En las siguientes iteraciones….
Valor para convergencia: 5.3256E-14
A 4.04537919 0.17582578
B 0.30607749 0.06761541
Z 0.36799325 0.91132737
Ma(lb-mol) 102.254472 16.4021621
Densidad(lb/pie3) 43.160999 2.79560441
VolumenMolar(pie3/lb
mol) 2.36914053 5.86712558
17