PROBLEMA 1: Separacin FlashUna cmara flash opera a 50 C y 200 kPa, separando 1000 kmol/h de una alimentacin que tiene 30 % mol de propano, 10 % n-butano, 15 % n-pentano, y 45 % n-hexano. Cules son las composiciones de los productos y los flujos. A travs de la Carta Depriester se determina que los valores k son: k1 (propano) = 7.0; k2 (n-butano) = 2.4; k3 (n-pentano) = 0.80; k4 (n-hexano) = 0.302.

Solucin

disp(' Calculo de la fraccion de vaporizacion')disp(' con el metodo de Newton Raphson ')disp(' ')n=input(' Ingrese el numero de componentes: ');z=input(' Ingrese la composicion de la alimentacion, [ ]: ');k=input(' Ingrese las constantes de equilibrio, [ ]: ');a=input(' Ingrese el valor inicial: ');suma1=0; iter=0;for i=1:n suma1=z(i)*(k(i)-1)/(1+a*(k(i)-1))+suma1;endfprintf(' N x fx\n');while abs(suma1)>1e-4dsuma1=0; for i=1:n dsuma1=z(i)*(k(i)-1)^2/((1+a*(k(i)-1)))^2+dsuma1; end ak=a+suma1/dsuma1; a=ak; suma1=0; for i=1:n suma1=z(i)*(k(i)-1)/(1+a*(k(i)-1))+suma1; end iter=iter+1; fprintf('%2d \t %11.8f \t %7.4e \n',iter,a,suma1);end disp(' ') fprintf(' La raiz es: %5.8f\n',a) fprintf(' El numero de iteraciones es: %3.0f\n',iter)