Scribd
Upload
34 views1 page

%parametros de Interaccion: For For End End For For en

This document contains constants and parameters for calculating properties of a two-component mixture using an activity coefficient model. It defines critical temperatures and pressures for the components, as well as Antoine constants for vapor pressure calculations. Parameters for an interaction model including Aij, Bij, and Cij values are provided to calculate activity coefficients. Mole fractions and a temperature are given as inputs, and the document outlines loops to calculate binary interaction coefficients and activity coefficients for each component in the mixture.

Uploaded by

Sergio Luna
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
34 views1 page

%parametros de Interaccion: For For End End For For en

This document contains constants and parameters for calculating properties of a two-component mixture using an activity coefficient model. It defines critical temperatures and pressures for the components, as well as Antoine constants for vapor pressure calculations. Parameters for an interaction model including Aij, Bij, and Cij values are provided to calculate activity coefficients. Mole fractions and a temperature are given as inputs, and the document outlines loops to calculate binary interaction coefficients and activity coefficients for each component in the mixture.

Uploaded by

Sergio Luna
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 1

Constancia de la reacción   

exotérmica Tcr=[521.64 562.05] %K 


Pcr=[8097 4895]  % kPa 
R= 8.314 % kJ/kmolK 
  
P=1000 %kPa 
  
X1=0 

% constantes Antoine 
Ai=[16.5938 13.8594]  %PARAMETROS DE INTERACCION 
Bi=[3644.3 2773.78]    
Ci=[239.76 220.07] 
Aij=[0 -0.3547 ;-2.8744 0] 
Bij=[0 389.079 ;1603.93 0] 
Cij=[0 0.4;0.4 0] 
  
Xi=[0.1 0.9] 
T= 60 %°C 
  
for i =1:NC 
   for j=1:NC 
     Tij(i,j)=Aij(i,j)+Bij(i,j)/(T+273.15) 
     Gij(i,j)=exp(-Cij(i,j)*Tij(i,j)) 
   end    
end 
        
for i =1:NC       
            gi(i)=exp((sum(TjiGjiXj
    for j=1:NC  )/sum(GkiXk))+sum(sumatoria)) 
       TjiGjiXj(j)=Tij(j,i)*Gi      
j(j,i)*Xi(j)       
         en
       for k=1 : NC 
          TkjGkjXk(k)=Tij(k,j)
*Gij(k,j)*Xi(k) 
          GkjXk(k)=Gij(k,j)*Xi
(k) 
          GkiXk(k)=Gij(k,i)*Xi
(k) 
       end 
       sumatoria(j)=((Xi(j)*Gi
j(i,j))/sum(GkjXk))*(Tij(i,j)-
sum(TkjGkjXk)/sum(GkjXk)) 
        
       end 

You might also like