Scribd
Upload
43 views

Para A Fase Biodiesel

This document contains code for calculating the flash calculations for a three-component system (glycerol, biodiesel, water) using the NRTL model. It provides the component compositions (xg and xb) in each phase, as well as performing a flash calculation to determine the vapor (z) and liquid (xl) phase compositions and the distribution coefficient (K).

Uploaded by

Rafael Cavalcante
Copyright
© © All Rights Reserved
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
43 views

Para A Fase Biodiesel

This document contains code for calculating the flash calculations for a three-component system (glycerol, biodiesel, water) using the NRTL model. It provides the component compositions (xg and xb) in each phase, as well as performing a flash calculation to determine the vapor (z) and liquid (xl) phase compositions and the distribution coefficient (K).

Uploaded by

Rafael Cavalcante
Copyright
© © All Rights Reserved
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 4

alpha = [0,0.2,0.2;0.2,0,0.2;0.2,0.

2,0];
A = [0,-624.48,352.24;1585.1,0,1313.3;2115.4,-737.53,0];
T = 298.15;
xg = [0.000499966
0.240702466 0.758797568; 0.000782409
0.349511685
0.649705906; 0.001301873 0.447998527 0.5506996; 0.002322981 0.537063359
0.460613659; 0.003281255 0.608428964 0.388289781; 0.005116127
0.670549514 0.324334359; 0.007173415 0.731493332 0.261333252];
xg = xg.';
xb = [0.772545845, 0.177495642,
0.049958513; 0.692971221,
0.265522817,
0.041505963; 0.623659809,
0.340525776,
0.035814415;
0.546742391,
0.422633039,
0.030624569; 0.489489813,
0.482824941,
0.027685246; 0.433515126,
0.541089406,
0.025395468; 0.381274393,
0.594541766,
0.024183841];
xb = xb.';
G = exp((-(alpha*A))/T);
lambda = (A*G)/T;
for i = 1:7;
E = lambda*diag((1./(G.'*xb(:,i))));
L = G*diag((1./(G.'*xb(:,i))));
lNphi = ((E + E.') - (L*diag(xb(:,i))*E.'))*xb(:,i);
PHIxb(:,i) = exp(lNphi);
end
PHIxb

Para a fase biodiesel:


PHIxb =
Columns 1 through 4
0.961129472247431
60929.682622926
0.00947827887009364

1.75951390689812
8537.11119321678
0.00695559497590399

2.81841316228288
2912.15090098583

4.40326150210996
1368.49707650677

0.00969515877150202

0.0201155080616873

Columns 5 through 7
5.79097378884541

7.22587257052674

8.49947812689779

955.083314871389

752.216245536153

648.468049184562

0.0400074454908165

0.0833189172600538

0.170504068274346

0.0234187340032305

0.061093731281491

Para a fase glicerol:


PHIxg =

Columns 1 through 4
0.00770797932904329
547.778827482095

487.119749168873

22.903332544756

23.991355245119

0.143433985237334

470.8181452788

470.252937407505

24.4774853358013

24.4193841827135

Columns 5 through 7
0.284009087104227

0.522826451545729

0.962996439434078

473.768905149596

479.216388880707

484.246293798811

24.0614837939054

23.4008927079974

22.5297502278963

Flash LL

z = [0.068902937; 0.712605254; 0.218491809];


alpha = [0,0.2,0.2;0.2,0,0.2;0.2,0.2,0];
A = [0,-624.48,352.24;1585.1,0,1313.3;2115.4,-737.53,0];
T = 298.15;
G = exp((-(alpha*A))/T);
lambda = (A*G)/T;
Ez = lambda*diag((1./(G.'*z)));
Lz = G*diag((1./(G.'*z)));
lNphiz = ((Ez + Ez.') - (Lz*diag(z)*Ez.'))*z;
PHIz = exp(lNphiz);
xl = (1./(PHIz.'*z))*(diag(PHIz)*z);
Exl = lambda*diag((1./(G.'*xl)));
Lxl = G*diag((1./(G.'*xl)));
lNphixl = ((Exl + Exl.') - (Lxl*diag(xl)*Exl.'))*xl;

PHIxl = exp(lNphixl);
kl = diag(1./PHIxl)*PHIz;
klmax = max(kl(:));
klmin = min(kl(:));
N1 = diag((klmax*(ones(3,1)))-kl)*z;
N2 = diag(kl-(klmin*(ones(3,1))))*z;
x1 = 1./((ones(1,3))*N1)*N1;
x2 = 1./((ones(1,3))*N2)*N2;
Ex1 = lambda*diag((1./(G.'*x1)));
Lx1 = G*diag((1./(G.'*x1)));
lNphix1 = ((Ex1 + Ex1.') - (Lx1*diag(x1)*Ex1.'))*x1;
Ex2 = lambda*diag((1./(G.'*x2)));
Lx2 = G*diag((1./(G.'*x2)));
lNphix2 = ((Ex2 + Ex2.') - (Lx2*diag(x2)*Ex2.'))*x2;
BETA = (1./(klmax - klmin))*((ones(1,3))*N2);
K = exp((lNphix1)-(lNphix2));
BETA1 = 0;
dBETA = BETA - BETA1;
BETA1 = BETA;
K1 = [0;0;0];
dK = K - K1;
K1 = K;
while dBETA > 0,01;
a = diag(1./((ones(3,1))+BETA1*(K1-(ones(3,1)))))*(K1-(ones(3,1)));
BETA = BETA1 + (1./(z.'*diag(a)*a))*(z.'*a);
dBETA = BETA - BETA1;
BETA1 = BETA;
end
BETA1
while dK > 0,01;
XI = diag(1./((ones(3,1))+BETA1*(K1-(ones(3,1)))))*z;
XII = diag(K1)*XI;
EXI = lambda*diag((1./(G.'*XI)));
LXI = G*diag((1./(G.'*XI)));
lNphiXI = ((EXI + EXI.') - (LXI*diag(XI)*EXI.'))*XI;
EXII = lambda*diag((1./(G.'*XII)));
LXII = G*diag((1./(G.'*XII)));
lNphiXII = ((EXII + EXII.') - (LXII*diag(XII)*EXII.'))*XII;
K = exp((lNphiXI)-(lNphiXII));
dK = K - K1;
K1 = K;
end
K1
XI
XII

BETA1 =
0.8564

K1 =
0.0373
54.5585
0.2492

XI =
0.4247
0.0139
0.5616

XII =
0.0092
0.8298
0.1610

You might also like