Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 473 vistas

    Tarea - 2 Ramos Requena F.

    Cargado por

    florencio ramos montano
    La oxidación de sulfuro sódico Na2S hasta tiosulfato Na2S2O3 fue estudiada por investigadores japoneses en un reactor intermitente durante 100 minutos. Se graficaron las concentraciones de cada componente en función del tiempo para analizar la reacción. Se desarrolló un programa en MATLAB para simular numéricamente la reacción usando el método de Runge-Kutta de cuarto orden.

    Copyright:

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

    Tarea - 2 Ramos Requena F.

    Cargado por

    florencio ramos montano
    473 vistas5 páginas
    La oxidación de sulfuro sódico Na2S hasta tiosulfato Na2S2O3 fue estudiada por investigadores japoneses en un reactor intermitente durante 100 minutos. Se graficaron las concentraciones de cada componente en función del tiempo para analizar la reacción. Se desarrolló un programa en MATLAB para simular numéricamente la reacción usando el método de Runge-Kutta de cuarto orden.

    Título original

    TAREA -2 RAMOS REQUENA F.

    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?

    La oxidación de sulfuro sódico Na2S hasta tiosulfato Na2S2O3 fue estudiada por investigadores japoneses en un reactor intermitente durante 100 minutos. Se graficaron las concentraciones de cada componente en función del tiempo para analizar la reacción. Se desarrolló un programa en MATLAB para simular numéricamente la reacción usando el método de Runge-Kutta de cuarto orden.

    Copyright:

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

    Tarea - 2 Ramos Requena F.

    Cargado por

    florencio ramos montano
    La oxidación de sulfuro sódico Na2S hasta tiosulfato Na2S2O3 fue estudiada por investigadores japoneses en un reactor intermitente durante 100 minutos. Se graficaron las concentraciones de cada componente en función del tiempo para analizar la reacción. Se desarrolló un programa en MATLAB para simular numéricamente la reacción usando el método de Runge-Kutta de cuarto orden.

    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 5

    INGENIERIA DE LAS REACCIONES QUIMICAS II

    Ingeniero: Salvador Ore Vidalon


    Alumno: RAMOS REQUENA, Frans Florencio
    Eacuela: Ing. Química

    TAREA No.2

    La oxidación de sulfuro sódico Na 2S (A) hasta tiosulfato Na2S2O3 (R) ha sido


    seguida por investigadores japoneses en un reactor intermitente. como se muestra
    a continuación.

    CA(0) = 185 mol/litro


    Cs(0) = 0
    CT(0) = 0
    CR(0) = 0

    Tiempo = 100 minutos

    Graficar tiempo vs. concentración de cada componente.


    SOLUCION:
    function RANGEKUTTAconforywhileparaREACCIONESCOMPLEJAS
    clc,clear
    n=input('Ingrese valor de n:');
    a=input('Ingrese valor inicial de "t":');
    b=input('Ingrese valor final de "t":');
    h=(b-a)/n;
    m=1;
    xo(m)=a;
    to(m)=a;
    CAo(m)=input('Ingresa valor inicial de "CA":');
    CSo(m)=input('Ingresa valor inicial de "CS":');
    CTo(m)=input('Ingresa valor inicial de "CT":');
    CRo(m)=input('Ingresa valor inicial de "CR":');
    syms t CA CS CT CR
    FA=input('Ingrese la función FA:');
    FS=input('Ingrese la función FS:');
    FT=input('Ingrese la función FT:');
    FR=input('Ingrese la función FR:');
    while m < (n+1)
    t=to(m);
    CA=CAo(m);
    CS=CSo(m);
    CT=CTo(m);
    CR=CRo(m);
    k1A(m)=eval(FA);
    k1S(m)=eval(FS);
    k1T(m)=eval(FT);
    k1R(m)=eval(FR);
    t=to(m)+h/2;
    CA=CAo(m)+h/2*k1A(m);
    CS=CSo(m)+h/2*k1S(m);
    CT=CTo(m)+h/2*k1T(m);
    CR=CRo(m)+h/2*k1R(m);
    k2A(m)=eval(FA);
    k2S(m)=eval(FS);
    k2R(m)=eval(FR);
    k2T(m)=eval(FT);
    t=to(m)+h/2;
    CA=CAo(m)+h/2*k2A(m);
    CS=CSo(m)+h/2*k2S(m);
    CT=CTo(m)+h/2*k2T(m);
    CR=CRo(m)+h/2*k2R(m);
    k3A(m)=eval(FA);
    k3S(m)=eval(FS);
    k3R(m)=eval(FR);
    k3T(m)=eval(FT);
    t=to(m)+h;
    CA=CAo(m)+h*k3A(m);
    CS=CSo(m)+h*k3S(m);
    CT=CTo(m)+h*k3T(m);
    CR=CRo(m)+h*k3R(m);
    k4A(m)=eval(FA);
    k4S(m)=eval(FS);
    k4R(m)=eval(FR);
    k4T(m)=eval(FT);
    to(m+1)=to(m)+h;
    CAo(m+1)=CAo(m)+h/6*(k1A(m)+2*k2A(m)+2*k3A(m)+k4A(m));
    CSo(m+1)=CSo(m)+h/6*(k1S(m)+2*k2S(m)+2*k3S(m)+k4S(m));
    CTo(m+1)=CTo(m)+h/6*(k1T(m)+2*k2T(m)+2*k3T(m)+k4T(m));
    CRo(m+1)=CRo(m)+h/6*(k1R(m)+2*k2R(m)+2*k3R(m)+k4R(m));
    m=m+1;
    end
    disp(' n t CA CS CT CR ')
    for a=1:m
    fprintf(' %3.4f %5.6f %5.6f %5.6f %5.6f
    %5.6f\n',a,to(a),CAo(a),CSo(a),CTo(a),CRo(a));
    figure (1)
    plot(to,CAo,'r-');
    plot(to,CSo,'g-');
    plot(to,CTo,'m-');
    plot(to,CRo,'k-');
    axis([0 100 0 200],'square')
    grid on
    hold on
    title('DIAGRAMA t-C ( "OXIDACIÓN DE SULFATO SÓDICO Na2S(A) HASTA
    TIOSULFATO Na2S2O3(R)"'),xlabel('t(min)'), ylabel('C (mol/L)')
    legend('t','CA','CS','CT','CR')

    plot(to(a),CAo(a),'r-',to(a),CSo(a),'g-',to(a),CTo(a),'m-',to(a),CRo(a),'
    k-');
    hold on
    end

    EJECUTANDO PROGRAMA:
    Ingrese valor de n:30
    Ingrese valor inicial de "t":0
    Ingrese valor final de "t":100
    Ingresa valor inicial de "CA":185
    Ingresa valor inicial de "CS":0
    Ingresa valor inicial de "CT":0
    Ingresa valor inicial de "CR":0
    Ingrese la función FA:-0.0205*CA
    Ingrese la función FS:0.0070*CA-0.0099*CS
    Ingrese la función FT:0.0027*CA-0.0163*CT
    Ingrese la función FR:0.0108*CA+0.0099*CS+0.0163*CT
    n t CA CS CT CR
    1.0000 0.000000 185.000000 0.000000 0.000000 0.000000
    2.0000 3.333333 172.780587 4.103616 1.565960 6.549837
    3.0000 6.666667 161.368277 7.802975 2.945674 12.883074
    4.0000 10.000000 150.709759 11.129103 4.155822 19.005315
    5.0000 13.333333 140.755247 14.110836 5.211753 24.922164
    6.0000 16.666667 131.458239 16.774970 6.127580 30.639212
    7.0000 20.000000 122.775306 19.146398 6.916280 36.162016
    8.0000 23.333333 114.665889 21.248244 7.589772 41.496095
    9.0000 26.666667 107.092106 23.101981 8.159003 46.646910
    10.0000 30.000000 100.018577 24.727543 8.634023 51.619857
    11.0000 33.333333 93.412262 26.143434 9.024047 56.420257
    12.0000 36.666667 87.242300 27.366824 9.337524 61.053352
    13.0000 40.000000 81.479869 28.413641 9.582198 65.524293
    14.0000 43.333333 76.098052 29.298656 9.765155 69.838138
    15.0000 46.666667 71.071708 30.035564 9.892881 73.999847
    16.0000 50.000000 66.377359 30.637058 9.971307 78.014276
    17.0000 53.333333 61.993076 31.114898 10.005849 81.886177
    18.0000 56.666667 57.898379 31.479976 10.001453 85.620193
    19.0000 60.000000 54.074140 31.742375 9.962629 89.220856
    20.0000 63.333333 50.502495 31.911429 9.893488 92.692588
    21.0000 66.666667 47.166761 31.995770 9.797770 96.039700
    22.0000 70.000000 44.051355 32.003380 9.678878 99.266387
    23.0000 73.333333 41.141724 31.941639 9.539903 102.376735
    24.0000 76.666667 38.424277 31.817361 9.383649 105.374714
    25.0000 80.000000 35.886319 31.636839 9.212655 108.264186
    26.0000 83.333333 33.515996 31.405882 9.029221 111.048901
    27.0000 86.666667 31.302235 31.129844 8.835425 113.732497
    28.0000 90.000000 29.234695 30.813661 8.633138 116.318507
    29.0000 93.333333 27.303717 30.461881 8.424047 118.810355
    30.0000 96.666667 25.500283 30.078687 8.209669 121.211361
    31.0000 100.000000 23.815967 29.667930 7.991363 123.524741

    GRAFICO:

    También podría gustarte