Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 75 vistas

    Programas de Equilibrio Liquido Vapor ISOBÁRICOS

    Cargado por

    Ronaldo Orihuela
    Este documento presenta el análisis termodinámico del equilibrio líquido-vapor de una mezcla ideal de acetato de etilo y etanol a presión constante utilizando los métodos de Van Laar, Margules y Wilson. Se proporcionan las ecuaciones y parámetros para cada modelo y se grafican los resultados frente a datos experimentales para evaluar su precisión.

    Copyright:

    © All Rights Reserved
    Descargue como DOCX, PDF, TXT o lea en línea desde Scribd

    Programas de Equilibrio Liquido Vapor ISOBÁRICOS

    Cargado por

    Ronaldo Orihuela
    75 vistas7 páginas
    Este documento presenta el análisis termodinámico del equilibrio líquido-vapor de una mezcla ideal de acetato de etilo y etanol a presión constante utilizando los métodos de Van Laar, Margules y Wilson. Se proporcionan las ecuaciones y parámetros para cada modelo y se grafican los resultados frente a datos experimentales para evaluar su precisión.

    Descripción original:

    deeeee

    Título original

    Programas-de-Equilibrio-liquido-vapor-ISOBÁRICOS

    Derechos de autor

    © © All Rights Reserved

    Formatos disponibles

    DOCX, PDF, TXT o lea en línea desde Scribd

    ¿Le pareció útil este documento?

    ¿Este contenido es inapropiado?

    Este documento presenta el análisis termodinámico del equilibrio líquido-vapor de una mezcla ideal de acetato de etilo y etanol a presión constante utilizando los métodos de Van Laar, Margules y Wilson. Se proporcionan las ecuaciones y parámetros para cada modelo y se grafican los resultados frente a datos experimentales para evaluar su precisión.

    Copyright:

    © All Rights Reserved
    Descargue como DOCX, PDF, TXT o lea en línea desde Scribd
    Descargar como docx, pdf o txt
    75 vistas7 páginas

    Programas de Equilibrio Liquido Vapor ISOBÁRICOS

    Cargado por

    Ronaldo Orihuela
    Este documento presenta el análisis termodinámico del equilibrio líquido-vapor de una mezcla ideal de acetato de etilo y etanol a presión constante utilizando los métodos de Van Laar, Margules y Wilson. Se proporcionan las ecuaciones y parámetros para cada modelo y se grafican los resultados frente a datos experimentales para evaluar su precisión.

    Copyright:

    © All Rights Reserved
    Descargue como DOCX, PDF, TXT o lea en línea desde Scribd
    Descargar como docx, pdf o txt
    de 7

    “AÑO DE LA DIVERSIFICACION PRODUCTIVA Y

    DEL FORTALECIMIENTO DE LA EDUCACION”

    UNIVERSIDAD NACIONAL DEL CENTRO DEL


    PERÚ
    FACULTAD DE INGENIERIA QUIMICA
    E.A.P. DE INGENIERIA QUIMICA

    TERMODINÁMICA DE LOS PROCESOS QUÍMICOS II

    Equilibrio Líquido-Vapor en un proceso


    isobárico

    CATEDRÁTICO:

    Ms. ORÉ VIDALÓN, Salvador Teódulo

    ALUMNA:

    CRISTOBAL ESPINOZA, Andrea Cristina

    SEMESTRE: VI SECCION: “A”

    Huancayo – Perú

    2015
    MÉTODO DE VAN LAAR
    PROCESO ISOBÁRICO: (ACETATO DE ETILO (1) – ETANOL (2))

    function VAN_LAAR_ISOBARICO
    clc, clear all
    % Ideal Isobárico
    % Equilibrio Liquido Vapor: Acetato de Etilo (1) – Etanol (2)
    syms T x1 x2
    disp('Sistema Isobarico:VAN LAAR Acetato de etilo(1) - Etanol(2)')
    X1=[0.0001,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1];
    X2=[1,0.9,0.8,0.7,0.6,0.5,0.4,0.3,0.2,0.1,0];
    % Constantes de antoine para el acetato de etilo y etanol:
    a1=7.09808; a2=8.04494 ;
    b1=1238.71; b2=1554.3 ;
    c1=217.0; c2=222.65 ;
    % Parámetros de Wilson para el acetato de etilo(1) - etanol(2):
    a12=0.3311;
    a21=0.3972;
    % Presión constante en mmHg:
    disp('a presion de: ')
    P=760
    f1=10^(a1-b1/(c1+T));
    f2=10^(a2-b2/(c2+T));
    f3=10^(a12/(1+x1*a12/(x2*a21)))^2;
    f4=10^(a21/(1+x2*a21/(x1*a12)))^2;
    s=12;r=1;p=1;
    while s>r
    x1=X1(p);
    x2=X2(p);
    q=1;v=0;
    tmin=40;
    while v<0.999
    T=tmin(q);
    pvp1=eval(f1);
    pvp2=eval(f2);
    gam1=eval(f3);
    gam2=eval(f4);
    Y1(q+1)=x1*pvp1*gam1/P;
    Y2(q+1)=x2*pvp2*gam2/P;
    v(q+1)= Y1(q+1)+Y2(q+1);
    tmin(q+1)=tmin(q)+0.05;
    q=q+1;
    end
    t(p)=tmin(q);
    y1(p)=Y1(q);
    y2(p)=Y2(q);
    p1(p)=pvp1;
    p2(p)=pvp2;
    u(p)=v(q);
    g1(p)=gam1;
    g2(p)=gam2;
    r(p+1)=r(p)+1;
    p=p+1;
    end
    disp(' x T pvp1 pvp2 P gamma1 gamma2
    y1 y2 y1+y2')
    for a=1:1:s-1;
    fprintf('%0.1f %5.2f %5.3f %5.3f %5.1f %5.3f %5.3f
    %5.3f %5.3f
    %5.2f\n',X1(a),t(a),p1(a),p2(a),P,g1(a),g2(a),y1(a),y2(a),u(a))
    end
    plot(X1,y1)
    subplot(2,1,1),plot(X1,y1,'c',X1,y1,'r*',X1,X1,'k')
    title('ISOBARICO: Equilibrio Liquido Vapor: Acetato de etilo – Etanol')
    xlabel('X1')
    ylabel('Y1')
    legend('VAN LAAR','DATOS EXPERIMENTALES','DIAGONAL')
    subplot(2,1,2),plot(X1,t,'c',X1,t,'r*',y1,t,'k',y1,t,'b+')
    title('ISOBARICO: Equilibrio Liquido Vapor: Acetato de etilo – Etanol')
    xlabel('X1,Y1')
    ylabel('T')
    legend('VAN LAAR','DATOS EXPERIMENTALES X1','VAN LAAR','DATOS
    EXPERIMENTALES Y1')
    MÉTODO DE MARGULES
    PROCESO ISOBÁRICO: (ACETATO DE ETILO (1) – ETANOL (2))
    function MARGULES_ISOBARICO
    clc, clear all
    % Ideal Isobárico
    % Equilibrio Liquido Vapor: Acetato de Etilo (1) – Etanol (2)
    syms T x1 x2
    disp('Sistema Isobarico: MARGULES Acetato de etilo(1) - Etanol(2)')
    X1=[0.0001,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1];
    X2=[1,0.9,0.8,0.7,0.6,0.5,0.4,0.3,0.2,0.1,0];
    % Constantes de antoine para el acetato de etilo y etanol:
    a1=7.09808; a2=8.04494 ;
    b1=1238.71; b2=1554.3 ;
    c1=217.0; c2=222.65 ;
    % Parámetros de Margules para el acetato de etilo(1) - etanol(2):
    a12=0.3313;
    a21=0.3925;
    % Presión constante en mmHg:
    disp('a presion de: ')
    P=760
    f1=10^(a1-b1/(c1+T));
    f2=10^(a2-b2/(c2+T));
    f3=10^(a12/(1+x1*a12/(x2*a21)))^2;
    f4=10^(a21/(1+x2*a21/(x1*a12)))^2;
    s=12;r=1;p=1;
    while s>r
    x1=X1(p);
    x2=X2(p);
    q=1;v=0;
    tmin=40;
    while v<0.999
    T=tmin(q);
    pvp1=eval(f1);
    pvp2=eval(f2);
    gam1=eval(f3);
    gam2=eval(f4);
    Y1(q+1)=x1*pvp1*gam1/P;
    Y2(q+1)=x2*pvp2*gam2/P;
    v(q+1)= Y1(q+1)+Y2(q+1);
    tmin(q+1)=tmin(q)+0.05;
    q=q+1;
    end
    t(p)=tmin(q);
    y1(p)=Y1(q);
    y2(p)=Y2(q);
    p1(p)=pvp1;
    p2(p)=pvp2;
    u(p)=v(q);
    g1(p)=gam1;
    g2(p)=gam2;
    r(p+1)=r(p)+1;
    p=p+1;
    end
    disp(' x T pvp1 pvp2 P gamma1 gamma2
    y1 y2 y1+y2')
    for a=1:1:s-1;
    fprintf('%0.1f %5.2f %5.3f %5.3f %5.1f %5.3f %5.3f
    %5.3f %5.3f
    %5.2f\n',X1(a),t(a),p1(a),p2(a),P,g1(a),g2(a),y1(a),y2(a),u(a))
    end
    plot(X1,y1)
    subplot(2,1,1),plot(X1,y1,'c',X1,y1,'r*',X1,X1,'k')
    title('ISOBARICO: Equilibrio Liquido Vapor: Acetato de etilo – Etanol')
    xlabel('X1')
    ylabel('Y1')
    legend('MARGULES','DATOS EXPERIMENTALES','DIAGONAL')
    subplot(2,1,2),plot(X1,t,'c',X1,t,'r*',y1,t,'k',y1,t,'b+')
    title('ISOBARICO: Equilibrio Liquido Vapor: Acetato de etilo – Etanol')
    xlabel('X1,Y1')
    ylabel('T')
    legend('MARGULES','DATOS EXPERIMENTALES X1','MARGULES','DATOS
    EXPERIMENTALES Y1')
    MÉTODO DE WILSON
    PROCESO ISOBÁRICO: (ACETATO DE ETILO (1) – ETANOL (2))

    function WILSON_ISOBARICO
    clc, clear all
    % Ideal Isobárico
    % Equilibrio Liquido Vapor: Acetato de Etilo (1) – Etanol (2)
    syms T x1 x2
    disp('Sistema Isobarico: WILSON Acetato de etilo(1) - Etanol(2)')
    X1=[0.0001,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1];
    X2=[1,0.9,0.8,0.7,0.6,0.5,0.4,0.3,0.2,0.1,0];
    % Constantes de antoine para el acetato de etilo y etanol:
    a1=7.09808; a2=8.04494 ;
    b1=1238.71; b2=1554.3 ;
    c1=217.0; c2=222.65 ;
    % Parámetros de Wilson para el acetato de etilo(1) - etanol(2):
    v1=95.3;
    v2=60.356;
    d1=-178.81;
    d2=844.69;
    % Presión constante en mmHg:
    disp('a presion de: ')
    P=760
    a12=(v2/v1)*exp(-(d1)/(1.987*T));
    a21=(v1/v2)*exp(-(d2)/(1.987*T));
    f1=10^(a1-b1/(c1+T));
    f2=10^(a2-b2/(c2+T));
    s=12;r=1;p=1;
    while s>r
    x1=X1(p);
    x2=X2(p);
    q=1;v=0;
    tmin=40;
    while v<0.999
    T=tmin(q);
    a12=(v2/v1)*exp(-(d1)/(1.987*T));
    a21=(v1/v2)*exp(-(d2)/(1.987*T));
    gam1=exp(-log(x1+a12*x2)+x2*(a12/(x1+a12*x2)-a21/(a21*x1+x2)));
    gam2=exp(-log(x2+a21*x1)-x1*(a12/(x1+a12*x2)-a21/(a21*x1+x2)));
    pvp1=eval(f1);
    pvp2=eval(f2);
    Y1(q+1)=x1*pvp1*gam1/P;
    Y2(q+1)=x2*pvp2*gam2/P;
    v(q+1)= Y1(q+1)+Y2(q+1);
    tmin(q+1)=tmin(q)+0.05;
    q=q+1;
    end
    t(p)=tmin(q);
    y1(p)=Y1(q);
    y2(p)=Y2(q);
    p1(p)=pvp1;
    p2(p)=pvp2;
    u(p)=v(q);
    g1(p)=gam1;
    g2(p)=gam2;
    r(p+1)=r(p)+1;
    p=p+1;
    end
    disp(' x T pvp1 pvp2 P gamma1 gamma2
    y1 y2 y1+y2')
    for a=1:1:s-1;
    fprintf('%0.1f %5.2f %5.3f %5.3f %5.1f %5.3f %5.3f
    %5.3f %5.3f
    %5.2f\n',X1(a),t(a),p1(a),p2(a),P,g1(a),g2(a),y1(a),y2(a),u(a))
    end
    plot(X1,y1)
    subplot(2,1,1),plot(X1,y1,'c',X1,y1,'r*',X1,X1,'k')
    title('ISOBARICO: Equilibrio Liquido Vapor: Acetato de etilo – Etanol')
    xlabel('X1')
    ylabel('Y1')
    legend('WILSON','DATOS EXPERIMENTALES','DIAGONAL')
    subplot(2,1,2),plot(X1,t,'c',X1,t,'r*',y1,t,'k',y1,t,'b+')
    title('ISOBARICO: Equilibrio Liquido Vapor: Acetato de etilo – Etanol')
    xlabel('X1,Y1')
    ylabel('T')
    legend('WILSON','DATOS EXPERIMENTALES X1','WILSON','DATOS EXPERIMENTALES
    Y1')

    También podría gustarte