% ob_03_tr2. C´alculo de destilaci´on binaria por el m´etodo de Sorel con SIMOBIQ
  
% Opción de cálculo de propiedades físicas

propfisop=0

% Opción de cálculo de equilibrio

flashop=0
%flashop=1 % Para usar Tb o Tr modificada

% Opción de cálculo para bloques de balance

modbalop=0

% Opción de cálculo para cálculo de destilación

destilacionop=0

% Opción de cálculo para dimensionamiento de equipos

dimeneqop=0

% Llamada a simobiq

  simobiq;
 
% Componentes 
comp=[C6H6,C7H8];
comp_Henry=[];

% Datos:
%  106 kmol/h benceno
%   35 kmol/h tolueno

f=[106, 35]
F=sum(f)
z=f/F

% La presión de la alimentaci´on es la atmosférica
PF=101.3e3 % Pa

% Temperatura  de la alimentación (líquido saturado o parcialmente vaporizado)
TF=tb(PF,z) % K
%TF=390 %K

CF={F,z,PF,TF}

% Las purezas deseadas son:
  xD = 0.995
  xB = 0.05
    
% Condensador
 FVD = 0 % Total
 %FVD =1 % Parcial
 
 % Presi´on de la columna
 P=101.3e3 % Atmosf´erica
 
 % Relación de reflujo (externa)
 R=3
 
tic
% C´alculo de la destilaci´on
[N,nf,DL,DV,B,QC,QH,PER] = sorel(CF,xB,xD,P,R,FVD)

tf=toc
   


