function f=ELVnPropisoPropanolf(w)

global z
global P
global P1vap
global P2vap

En este caso deben satisfacerse las ecuaciones del balance de masa y del equilibrio líquido-vapor:

  1. $$ x_1*P^o_1+(1-x_1)P^o_2 -P =0$$
  2. $$ x_1*P^o_1-y_1P=0 $$
  3. $$ f_{vap}y_1+(1-f_{vap})x_1-z_1=0 $$

Se tienen tres ecuaciones con tres incógnitas: x_1, y_1, fvap w: vector de incógnitas

  1. w(1)=x_1 fracción molar del más volátil en el líquido
  2. w(2)=y_1 fracción molar del más volátil en el vapor
  3. w(3)= fracción de vaporización de la alimentación
f(1)=w(1)*P1vap+(1-w(1))*P2vap-P;
f(2)=w(1)*P1vap-w(2)*P;
f(3)=w(3)*w(2)+(1-w(3))*w(1)-z(1);

Error using ELVnPropisoPropanolf (line 21)
Not enough input arguments.


Published with MATLAB® R2012b