Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 74 vistas

    Reactores No Isotérmicos. Ej 7

    Cargado por

    Samanta Silva
    Este documento presenta la resolución de un reactor tipo batch no isotérmico mediante el uso de MATLAB. Se desarrolla un modelo matemático que incluye ecuaciones para la velocidad de reacción, balance de energía y temperatura en función de la conversión. Luego, se presenta un algoritmo y su programación en MATLAB para calcular el tiempo requerido para alcanzar un 90% de conversión, el cual resulta ser de 146 minutos. Finalmente, se grafica la conversión en función del tiempo.

    Copyright:

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

    Reactores No Isotérmicos. Ej 7

    Cargado por

    Samanta Silva
    74 vistas13 páginas
    Este documento presenta la resolución de un reactor tipo batch no isotérmico mediante el uso de MATLAB. Se desarrolla un modelo matemático que incluye ecuaciones para la velocidad de reacción, balance de energía y temperatura en función de la conversión. Luego, se presenta un algoritmo y su programación en MATLAB para calcular el tiempo requerido para alcanzar un 90% de conversión, el cual resulta ser de 146 minutos. Finalmente, se grafica la conversión en función del tiempo.

    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 la resolución de un reactor tipo batch no isotérmico mediante el uso de MATLAB. Se desarrolla un modelo matemático que incluye ecuaciones para la velocidad de reacción, balance de energía y temperatura en función de la conversión. Luego, se presenta un algoritmo y su programación en MATLAB para calcular el tiempo requerido para alcanzar un 90% de conversión, el cual resulta ser de 146 minutos. Finalmente, se grafica la conversión en función del tiempo.

    Copyright:

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

    Reactores No Isotérmicos. Ej 7

    Cargado por

    Samanta Silva
    Este documento presenta la resolución de un reactor tipo batch no isotérmico mediante el uso de MATLAB. Se desarrolla un modelo matemático que incluye ecuaciones para la velocidad de reacción, balance de energía y temperatura en función de la conversión. Luego, se presenta un algoritmo y su programación en MATLAB para calcular el tiempo requerido para alcanzar un 90% de conversión, el cual resulta ser de 146 minutos. Finalmente, se grafica la conversión en función del tiempo.

    Copyright:

    © All Rights Reserved
    Descargue como DOCX, PDF, TXT o lea en línea desde Scribd
    Descargar como docx, pdf o txt
    Está en la página 1de 13

    ESCUELA POLITCNICA NACIONAL

    FACULTAD DE INGENIERA QUMICA


    Y AGROINDUSTRA
    CARRERA DE INGENIERA QUMICA
    INGENIERA DE LA REACCIN I

    Resolucin de un reactor no
    isotrmico en MATLAB
    INTEGRANTES:
    Barreto Paola
    Delgado Cindy
    Proao Daro
    Silva Samanta
    Fecha de Entrega: 18 Febrero 2015

    1. OBJETIVO
    Resolver el ejercicio planteado acerca de un reactor tipo batch no isotrmico,
    mediante el uso del software de modelado matemtico MATLAB, desarrollando un
    algoritmo y su respectiva programacin.
    Objetivos secundarios.

    Calcular el tiempo que toma la reaccin en un reactor tipo batch no


    isotrmico para las condiciones dadas.
    Determinar grficamente la funcin entre la conversin y el tiempo para el
    reactor planteado.
    Aplicar los conocimientos adquiridos durante todo el curso de Ingeniera de
    la Reaccin I.

    2. PLANTEAMIENTO DEL PROBLEMA


    Calcular el tiempo en el cual se obtiene una conversin del 90% para la reaccin
    2 A B +C , que se lleva a cabo en fase gaseosa en un reactor tipo batch, en
    el cual la temperatura no permanece constante y se trata de un sistema
    adiabtico.
    Evaluar de manera grfica el cambio de la conversin en dicho reactor, con
    respecto al tiempo, hasta que la conversin sea total. Conociendo que:
    =250 C

    P=2 atm
    x=0,9

    k =31 0 e

    8000
    RT

    L
    molmin

    H rx a 25 C =17000

    R=1,987

    Kcal
    Kmol

    cal
    molK

    Tabla 1. Datos de capacidad calorfica a presin constante

    Compuestos
    A
    B
    C

    Cp (cal/mol C)
    4
    3
    2

    3. ESQUEMA

    Figura 1. Esquema del reactor BATCH adiabtico

    4. DESARROLLO DEL MODELO MATEMTICO


    Ecuacin del tiempo de reaccin en un reactor batch en funcin de la conversin.
    x

    t=N Ao
    0

    dx
    r AV

    [1]

    Concentracin inicial del gas A, debemos tener en cuenta que el volumen


    permanece constante y que la fraccin molar de A es 1.
    C Ao =

    N Ao
    V

    [2]

    C Ao =

    C Ao =

    P
    RT

    2 atm
    Latm
    0,082
    523 K
    Kmol

    C Ao =0,047

    mol
    L

    Velocidad de reaccin del gas A y concentracin de A en funcin de la conversin.


    r A =kC A 2
    C Ao( A + Ax )T o
    P
    T
    CA=
    Po

    A=

    C Ao
    =1
    C Ao

    A=

    2
    =1
    2

    CA=

    C Ao( 1x )T o
    T
    C ( 1x )T o
    r A =k Ao
    T

    Reemplazar las ecuaciones 3 y 2 en 1

    [3]

    t=C Ao
    0

    dx

    C ( 1x )T o
    k Ao
    T

    [4]

    Balance de energa para la reaccin.


    n

    QWsF Ao i C pidT H rx a T F Aox=0


    i=1

    i C pidT = H rx a T x
    i =1

    [5]

    Sabiendo que el calor (Q) y el trabajo (Ws) son iguales a 0 debido a que el sistema
    es adiabtico, se tiene la resolucin de la entalpia de reaccin en funcin de la
    temperatura final de reaccin.
    T

    H rx a T = H rxa 25 C + CpdT
    TR

    iC p sale iC p entra
    i=1
    n

    Cp=
    i=1

    1
    1
    cal
    Cp= 3+ 214=1,5
    2
    2
    molK
    T

    cal
    cal
    H rx a T =17000
    + 1,5
    dT
    mol 298
    mol

    H rx a T =1,5T 16553

    cal
    mol

    [6]

    Al reemplazar la ecuacin 6 en la ecuacin 5 se obtiene la ecuacin 7, en donde


    se muestra la temperatura en funcin de la conversin.

    4
    523

    cal
    cal
    dT = 1,5T 16553
    x
    molK
    mol

    T=

    0,67(16553x +2092)
    x2,67

    [7]

    Finalmente se reemplaza la ecuacin 7 en la ecuacin 4.


    x

    t=

    C AoT o2 0

    0,67(16553x+2092)
    dx
    x2,67
    2

    31 0 e

    0,9

    t=

    0,047523 2 0

    8000
    R0,67(16553 x+2092)
    x2,67

    ( 1x )

    0,67(16553x +2092)
    dx
    x2,67

    31 02e

    8000
    1,9870,67(16553 x+2092)
    x2,67

    ( 1x )

    Las ecuaciones detalladas sern resueltas en el software.

    5. ALGORITMO
    El algoritmo que se utiliza para la programacin en MATLAB se presenta a
    continuacin.

    INICIO

    ,,,,,,,

    FIN

    6. PROGRAMACIN
    La programacin que se realiz en el software MATLAB se presenta a
    continuacin.

    % Reactor "batch" no isotermico y adiabatico


    clear;
    clc;
    syms CAo CBo CCo x xA t To T rA k P R R1 CpA CpB CpC thA thB thC a b c ga
    gb gc CA Q W FAo y DTr DT DCp
    % Constantes
    P=2; % atm
    x=0.9;
    To=523; % K
    R=0.082; % (L*atm)/(mol*K)
    R1=1.987; % cal/(mol*K)
    DTr=-17000; % cal/mol
    CpA=4; % cal/(mol*K)
    CpB=3; % cal/(mol*K)
    CpC=2; % cal/(mol*K)
    a=2;
    b=1;
    c=1;
    Q=0;
    W=0;
    CBo=0;
    CCo=0;
    % Resolucin
    thA=CAo/CAo
    thB=CBo/CAo
    thC=CCo/CAo
    gA=-a/a
    gB=b/a
    gC=c/a
    DCp=gB*CpB+gC*CpC+gA*CpA % cal/(mol*K)
    DT=DTr+int(DCp,298,T) %cal/mol
    y=thA*int(CpA,To,T)+thB*int(CpB,To,T)+thC*int(CpC,To,T) % cal/(mol*K)
    T=solve(Q-W-FAo*sum(y)-DT*xA*FAo==0,T)
    CAo=P/(R*To) % mol/L
    CA=(CAo*(thA+gA*xA)*To)/T
    k=3*10^2*exp(-8000/(R1*T)) % L/(mol*min)
    rA=k*CA^2
    t=CAo*int(1/rA,0,x);
    t=subs(t);
    t=vpa(t,'3') % min
    % Grfico
    syms tx t1 t2 t3 t4 t5 t6 t7 t8 t9 t10 xt x1 x2 x3 x4 x5 x6 x7 x8 x9 x10
    xB
    xt=0.01:0.1:1
    x1=0.01;
    x2=0.1;
    x3=0.2;
    x4=0.3;
    x5=0.4;
    x6=0.5;
    x7=0.6;
    x8=0.7;

    x9=0.8;
    x10=0.9;
    t1=CAo*int(1/((300*exp((8000000*(3*xB - 8))/(1987*(33106*xB +
    4184)))*(3*xB - 8)^2*((1000*xB)/41 - 1000/41)^2)/(33106*xB +
    4184)^2),0,x1);
    t2=CAo*int(1/((300*exp((8000000*(3*xB - 8))/(1987*(33106*xB +
    4184)))*(3*xB - 8)^2*((1000*xB)/41 - 1000/41)^2)/(33106*xB +
    4184)^2),0,x2);
    t3=CAo*int(1/((300*exp((8000000*(3*xB - 8))/(1987*(33106*xB +
    4184)))*(3*xB - 8)^2*((1000*xB)/41 - 1000/41)^2)/(33106*xB +
    4184)^2),0,x3);
    t4=CAo*int(1/((300*exp((8000000*(3*xB - 8))/(1987*(33106*xB +
    4184)))*(3*xB - 8)^2*((1000*xB)/41 - 1000/41)^2)/(33106*xB +
    4184)^2),0,x4);
    t5=CAo*int(1/((300*exp((8000000*(3*xB - 8))/(1987*(33106*xB +
    4184)))*(3*xB - 8)^2*((1000*xB)/41 - 1000/41)^2)/(33106*xB +
    4184)^2),0,x5);
    t6=CAo*int(1/((300*exp((8000000*(3*xB - 8))/(1987*(33106*xB +
    4184)))*(3*xB - 8)^2*((1000*xB)/41 - 1000/41)^2)/(33106*xB +
    4184)^2),0,x6);
    t7=CAo*int(1/((300*exp((8000000*(3*xB - 8))/(1987*(33106*xB +
    4184)))*(3*xB - 8)^2*((1000*xB)/41 - 1000/41)^2)/(33106*xB +
    4184)^2),0,x7);
    t8=CAo*int(1/((300*exp((8000000*(3*xB - 8))/(1987*(33106*xB +
    4184)))*(3*xB - 8)^2*((1000*xB)/41 - 1000/41)^2)/(33106*xB +
    4184)^2),0,x8);
    t9=CAo*int(1/((300*exp((8000000*(3*xB - 8))/(1987*(33106*xB +
    4184)))*(3*xB - 8)^2*((1000*xB)/41 - 1000/41)^2)/(33106*xB +
    4184)^2),0,x9);
    t10=CAo*int(1/((300*exp((8000000*(3*xB - 8))/(1987*(33106*xB +
    4184)))*(3*xB - 8)^2*((1000*xB)/41 - 1000/41)^2)/(33106*xB +
    4184)^2),0,x10);
    t1=subs(t1);
    t1=vpa(t1,'3'); % min
    t2=subs(t2);
    t2=vpa(t2,'3'); % min
    t3=subs(t3);
    t3=vpa(t3,'3'); % min
    t4=subs(t4);
    t4=vpa(t4,'3'); % min
    t5=subs(t5);
    t5=vpa(t5,'3'); % min
    t6=subs(t6);
    t6=vpa(t6,'3'); % min
    t7=subs(t7);
    t7=vpa(t7,'3'); % min
    t8=subs(t8);
    t8=vpa(t8,'3'); % min
    t9=subs(t9);
    t9=vpa(t9,'3'); % min
    t10=subs(t10);
    t10=vpa(t10,'3'); % min
    tx=[t1, t2, t3, t4, t5, t6, t7, t8, t9, t10]
    plot(tx, xt)
    legend('xA')
    grid on

    7. RESULTADOS Y DISCUSIN
    Los resultados obtenidos en el software MATLAB se presentan a continuacin.
    thA =
    1
    thB =
    0
    thC =
    0

    Nota: thA=A, thB=B y thC=C

    gA =
    -1
    gB =
    0.5000
    gC =
    0.5000

    Nota: gA=A, gB=B y gC=C

    DCp =
    -1.5000

    Nota: DCp=Cp

    DT =
    - (3*T)/2 16553

    H rx a T

    Nota: DT=

    y=
    n

    4*T 2092

    T=

    Nota: y=

    i=1

    i C p idT

    -(33106*xA + 4184)/(3*xA - 8)

    Nota: Temperatura en funcin de la conversin


    T=f(x)

    CAo =
    0.0466
    CA =
    ((3*xA - 8)*((1000*xA)/41 - 1000/41))/(33106*xA + 4184)
    k=
    300*exp((8000000*(3*xA - 8))/(1987*(33106*xA + 4184)))
    rA =
    (300*exp((8000000*(3*xA - 8))/(1987*(33106*xA + 4184)))*(3*xA 8)^2*((1000*xA)/41 - 1000/41)^2)/(33106*xA + 4184)^2
    t=
    146.0

    Nota: tiempo (min) que se demora en alcanzar la conversin de 90%

    xt =
    0.0100 0.1100
    0.8100 0.9100

    0.2100

    0.3100

    0.4100

    0.5100

    0.6100

    0.7100

    tx =
    [ 1.2834032854298129677772521972656,
    4.9751985737821087241172790527344,
    6.5017320675542578101158142089844,
    7.9485056351404637098312377929688,
    9.8487003440968692302703857421875,
    12.785931311314925551414489746094,
    17.918300170917063951492309570312,
    28.170443774666637182235717773438, 53.18107000552117824554443359375,
    146.1510930545628070831298828125]

    Tabla 2. Resultados obtenidos en el software MATLAB


    Tiempo para una conversin del 90%
    146. 0 min

    ( t
    Grfica obtenida por el software MATLAB

    Figura 2. Grafico conversin (xA) versus tiempo (min).


    Discusin de resultados
    De los resultados obtenidos con el uso del software utilizado, como se observa en
    la tabla 2 se tiene que una conversin del 90% para la reaccin dada se alcanza
    en un tiempo de aproximadamente 2 horas y 25 minutos, lo que evidencia que es
    una reaccin, tericamente, un tanto lenta.
    En la figura 2, por otra parte se muestra la dependencia de la conversin con el
    tiempo, en donde, en los primeros 50 minutos se tiene un gran aumento de la
    conversin, en aproximadamente un 80%. A partir de este tiempo, hasta
    aproximadamente 150 minutos, la conversin aumenta cerca de un 10%, y
    superado este tiempo, al acercase la conversin al 100%, el tiempo de reaccin
    tiende al infinito.
    Es importante tener en cuenta que las ecuaciones resultantes del desarrollo
    matemtico son bastante complejas y pueden presentar ciertos errores en su
    resolucin debido a los mtodos de aproximacin que usan los softwares de
    simulacin matemtica.

    CONCLUSIONES Y RECOMENDACIONES

    El tiempo en el reactor no isotrmico para una conversin del 90 % es de 146


    min.
    El tiempo aumenta con la conversin hasta que tiende al infinito.

    Se recomienda utilizar otros mtodos numricos ms avanzados debido a que


    en el software utilizado las ecuaciones de integral indefinida que involucran
    expresiones exponenciales son susceptibles de error y no expresa el resultado
    exacto.

    También podría gustarte