TD 5 GC

TD 5 : EQUILIBRE FLASH ISOTHERME LIQUIDE VAPEUR
LIEN THERMODYNAMIQUE NIST (Pour Vérification)

https://webbook.nist.gov/chemistry/fluid/
EXERCICE 1 (Traité au cours )

1/ Reprendre l’exercice traité au cours en utilisant Excel
A B C D E F G H
1 v = 0.425837
2 Phase
3 Equil. Terml = Term2 = Ratio =
4 zi K-Value (Ki−l)zi (Ki−l)v+l Terml/Term2 xi
yi
5 Propane 0.1 6.8 0.58 3.4699 0.1672 0.0288 0.1960
6 n-Butane 0.3 2.2 0.36 1.5110 0.2383 0.1985 0.4368
7 n-Pentane 0.4 0.8 −0.08 0.9148 −0.0874 0.4372 0.3498
8 n-Octane 0.2 0.052 −0.1896 0.5963 −0.3180 0.3354 0.0174
1 f (v) = 8.2719E-07 0.9999996 1.0
2/ Donner le polynôme du troisième degré solution au problème et le résoudre par la méthode de Newton.
3/ Ecrire le programme Matlab suivant. Analyser le programme et indiquer ce qu’il réalise. Vérifier les résultats précédents Quel est l’estimé initial qui produit le meilleur résultat.
function Y=flash_equ_1(K,z,x0) set(gca,'fontsize',10,'fontweight','bold')

close all h=xlabel('v');
F=@(x) sum(((K-1).*z)./(1+(K-1)*x)); set(gca,'fontsize',10,'fontweight','bold')
[v,feval]=fsolve(F,x0); h=ylabel('f(v)')
x=z./(1+(K-1)*v); set(gca,'fontsize',10,'fontweight','bold')
Y.v=v; hold on
Y.x=x; V=axis;
Y.y=K.*x; h=plot(v,feval,'og',[0 v],[feval feval],':r',[v v],[ feval V(3)],':r');
Y.FX=feval; set(h,'MarkerSize',10,'MarkerFaceColor','g','linewidth',2)
xx=linspace(0,1,100); shg
for ii=1:100 h=text(.5,-1,['\nu = ',num2str(v,10)]);
yy(ii)=sum(((K-1).*z)./(1+(K-1)*xx(ii))); set(h,'fontsize',12,'fontweight','bold','color','k','FontAngle','Itali
end c')
h=plot(xx,yy,'b'); h=text(.5,-1.5,['f(\nu) = ',num2str(feval,10)]);
set(h,'linewidth',2); set(h,'fontsize',12,'fontweight','bold','color','k','FontAngle','Itali
grid on c')
EXERCICE 2
Un mélange est à la pression de 100 psia et une température de 178◦F. Calculer la fraction vapeur en résolvant l’équation de Rachford–Rice
utilisant
a) Excel
b) MATLAB.
Composant Débit (lb mol/h) Coefficient de partage (K)
Propane 20 3.7
n-Butane 30 1.4
n-Pentane 50 0.6
Solution : v: 0.70789, x: [0.06869 0.233798 0.6975] y: [0.2541 0.327317 0.41850] fX=1.45 e-6
EXERCICE 3
Considérant un mélange de 78% N2, 21% O2, et 1% Ar à une pression de 2 atmosphères et une température de 86 °K. Barron (1985) donne
pour K-les valeurs respectives de 1.23, 0.37, et 0.47, Déterminer les compositions liquide- vapeur à l’équilibre sous ces conditions.
Solution : v: 0.29027, x: [0.7311 0.2569 0.01181] y: [0.8993 0.09508 0.00554] fX=-1.5339 e-8
v: 2.902742037903515e-001
x: [7.311839191127075e-001 2.569979150563136e-001 1.181817028377174e-002]
y: [8.993562205086303e-001 9.508922857083604e-002 5.554540033372715e-003]
FX: -1.533995388353082e-008
EXERCICE 4
Tracé de la courbe V=f(P) pour un mélange, à 60 °C : Les coefficients K sont donnés par la relation :
𝑃𝑃 𝑦𝑦
𝐾𝐾𝑖𝑖 = 𝑃𝑃𝑖𝑖 = 𝑥𝑥𝑖𝑖 Avec Pi pression de vapeur du composant i, et P la pression du mélange. Pi est donné par la relation :
𝑖𝑖
Log10(Pi)=Ai-Bi/(Ci+T)
Où Ai, Bi, Ci sont les coefficients d’Antoine ( tabulés)
A partir du tableur Excel et de son solveur, déterminer la fraction vapeur v=V/F
T 60
ETHYLENE ETHANE PROPANE N-BUTANE
A 3,8669 3,93264 3,977221 3,84431
B 584,146 659,739 819,296 909,65
C -18,307 -16,719 -24,417 -36,146
z 0,1 0,3 0,4 0,2
logPi 2,0115435 1,847702051 1,3234846 0,781556636
P k1 k2 k3 k4 x1 x2 x3 x4 y1 y2 y3 y4 v f(v)
18 5,7052019 3,912276549 1,170070675 0,335957338 0,0298 0,122 0,369 0,299 0,17 0,4779 0,431 0,1006 0,5 -0,35993421
19 5,4049281 3,706367257 1,108488008 0,318275373 0,0312 0,127 0,379 0,303 0,17 0,4725 0,421 0,0966 0,5 -0,31688327
20 5,1346817 3,521048894 1,053063607 0,302361604 0,0326 0,133 0,39 0,307 0,17 0,4673 0,41 0,0929 0,5 -0,27577994
21 4,890173 3,353379899 1,002917721 0,287963432 0,034 0,138 0,399 0,311 0,17 0,4622 0,401 0,0894 0,5 -0,23647213
22 4,6678924 3,20095354 0,957330552 0,274874186 0,0353 0,143 0,409 0,314 0,16 0,4572 0,391 0,0862 0,5 -0,19882484
INDICATIONS !
1/ Fonctions à introduire
Cellule B7 (logP1) : =B3-B4/(B5+$B$1+273,15), Cellule B9 K1=10^B$7/$A9, Cellule E9 x1 =B$6/(1+(B9-1)*$N9) ; Cellule J9 y1=B9*F9
2/ Vérifier que les solutions v pour les pressions de 10 à 21 atmosphères sont données par :
P k1 k2 k3 k4 x1 x2 x3 x4 y1 y2 y3 y4 v f(v)
18 5,70520187 3,91227655 1,17007067 0,33595734 0,0180023 0,0785501 0,3434548 0,5599486 0,1027066 0,3073096 0,4018664 0,1881188 0,96804674 -4,5711E-05
19 5,40492808 3,70636726 1,10848801 0,31827537 0,0199221 0,0864654 0,3639683 0,5292159 0,1076775 0,3204726 0,4034545 0,1684364 0,91251275 -0,00046928
20 5,13468168 3,52104889 1,05306361 0,3023616 0,0219113 0,0945478 0,3825051 0,5016599 0,1125076 0,3329076 0,4028022 0,1516827 0,86194159 0,00072413
21 4,89017303 3,3533799 1,00291772 0,28796343 0,024076 0,1031732 0,3990562 0,4730384 0,117736 0,3459791 0,4002205 0,1362178 0,81063439 -0,00080948
22 4,66789244 3,20095354 0,95733055 0,27487419 0,0263134 0,1119247 0,4134696 0,4480445 0,1228281 0,3582657 0,3958271 0,1231559 0,7634755 -0,00032453
3/ Vérifier les contraintes sur les fractions molaires des phase liquide et vapeur
4/ Tracer la courbe v=f(P)
5/ Quelle loi de dépendance suggérez vous ?
6/ Donner les coefficients de régression
I. SEPARATION FLASH ISOTHERME —EQUATIONS DEVELOPPEES
Considérons le schéma de traitement présenté à la figure 3.1.
Figure 3.1 : Séparateur de Phase Flash
Supposons que la température, la pression, et la composition générale du flux d'entrée sont connues. Les fractions molaires des produits chimiques dans
l'entrée sont notées par {zi}. Dans le séparateur de phases, cependant, le liquide et la vapeur sont séparés. Les fractions molaires des produits chimiques en
phase vapeur sont désignées par {yi} et celles du liquide phase par {xi}. Lorsque la vapeur et le liquide sont en équilibre, on peut relier les fractions molaires
de chaque produit chimique dans la vapeur et le liquide par l'équation :
yi = K i xi (3.1)
Les ensembles de {Ki} sont appelés valeurs K et peuvent être prédits à partir de la thermodynamique comme indiqué dans la section suivante. Pour l'instant,
supposons que leurs valeurs sont connues.
Pour calculer l'équation régissant le phénomène, on doit d'abord ajouter les fractions molaires de vapeur et liquide de tous les composants:
NC OMP NC OMP
� yi = 1, � xi = 1
i=1 i=1
(3.2)
Ensuite, soustrayant ces deux équations:
NC OMP NC OMP
� yi − � xi = 0
i=1 i=1
(3.3)
Ensuite, après avoir substitué l'équation expression d'équilibre et réarrangé, il vient :
NC OMP NC OMP NC OMP
� K i xi − � xi = � (K i − 1)xi = 0
i=1 i=1 i=1
(3.4)
Ensuite, effectuons un bilan de masse pour chaque composant sur le séparateur de phases. F est le débit molaire total, V est le débit molaire de vapeur et L est
le débit molaire de liquide. La balance molaire est alors :
zi 𝑭𝑭 = yi 𝑽𝑽 + xi L
Et :
𝐹𝐹 = 𝑉𝑉 + 𝐿𝐿
(3.5)
𝑉𝑉
Divisant par F et notant par v la fraction du débit molaire vapeur (𝜈𝜈 = )
𝐹𝐹
𝐳𝐳𝐢𝐢 = 𝐲𝐲𝐢𝐢 𝛎𝛎 + 𝐱𝐱𝐢𝐢 (𝟏𝟏 − 𝛎𝛎)

(3.6)
Exprimant la valeur de yi en fonction de xi, il vient :
𝐳𝐳𝐢𝐢
𝐱𝐱𝐢𝐢 =
(𝐊𝐊 − 𝟏𝟏)𝛖𝛖 + 𝟏𝟏
(3.7)
Portant cette expression dans (3.4), on aboutit à l’équation finale :
NC OMP
(K i − 1)zi
� =0
(𝐊𝐊 − 𝟏𝟏)𝛖𝛖 + 𝟏𝟏
i=1

TD 5 GC

Transféré par

Droits d'auteur :

Formats disponibles

TD 5 GC

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

TD 5 GC

Transféré par

Droits d'auteur :

Formats disponibles

TD 5 : EQUILIBRE FLASH ISOTHERME LIQUIDE VAPEUR

LIEN THERMODYNAMIQUE NIST (Pour Vérification)

EXERCICE 1 (Traité au cours )

function Y=flash_equ_1(K,z,x0) set(gca,'fontsize',10,'fontweight','bold')

Considérons le schéma de traitement présenté à la figure 3.1.

Figure 3.1 : Séparateur de Phase Flash

𝐳𝐳𝐢𝐢 = 𝐲𝐲𝐢𝐢 𝛎𝛎 + 𝐱𝐱𝐢𝐢 (𝟏𝟏 − 𝛎𝛎)

Portant cette expression dans (3.4), on aboutit à l’équation finale :

Vous aimerez peut-être aussi