Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 57 vistas

    Primer Examen

    Cargado por

    jvalencia_196986
    Este documento presenta la resolución de cuatro ecuaciones no lineales utilizando el método de la secante en MATLAB. También resuelve un sistema de ecuaciones químicas aplicando el método de Gauss-Seidel y calcula un factor de compresibilidad usando Newton-Raphson.

    Copyright:

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

    Primer Examen

    Cargado por

    jvalencia_196986
    57 vistas9 páginas
    Este documento presenta la resolución de cuatro ecuaciones no lineales utilizando el método de la secante en MATLAB. También resuelve un sistema de ecuaciones químicas aplicando el método de Gauss-Seidel y calcula un factor de compresibilidad usando Newton-Raphson.

    Título original

    primer examen

    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 cuatro ecuaciones no lineales utilizando el método de la secante en MATLAB. También resuelve un sistema de ecuaciones químicas aplicando el método de Gauss-Seidel y calcula un factor de compresibilidad usando Newton-Raphson.

    Copyright:

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

    Primer Examen

    Cargado por

    jvalencia_196986
    Este documento presenta la resolución de cuatro ecuaciones no lineales utilizando el método de la secante en MATLAB. También resuelve un sistema de ecuaciones químicas aplicando el método de Gauss-Seidel y calcula un factor de compresibilidad usando Newton-Raphson.

    Copyright:

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

    UNIVERSIDAD NACIONAL DE TRUJILLO

    PROFESOR:
    Dr. EVANGELISTA BENITES, GUILLERMO
    DAVID
    ALUMNO:
    VALENCIA CHAVARRIA, JHON ALONSO
    CURSO:
    MTODOS NUMRICOS
    CICLO:
    VI

    2015
    Ejercicio 2.16 del libro de Nieves Domnguez
    Resuelva por el mtodo de la secante, posicin falsa o biseccin, las siguientes
    ecuaciones:
    a) ex + 2-x + 2cosx 6 = 0

    b) xlogx 10 = 0
    c) ex + x3 + 2x2 + 10x 20 = 0
    d) senx cscx + 1 = 0
    Eh escogido el mtodo de la secante para resolver estas ecuaciones:
    PROGRAMACIN EN MATLAB:
    Para la ecuacin (a):
    clc, clear all
    fprintf('\tPOBLEMA 2.16 PARTE a MTODO DE LA SECANTE\n')
    fprintf('\t=========================================\n')
    x0=0; x1=1; i=0;
    fprintf('
    ----------------------------\n')
    fprintf('
    i
    Xi
    |Xi+1 - Xi|\n')
    fprintf('
    ----------------------------\n')
    fprintf('
    %d
    %8.5f\n',i,x0)
    i=i+1;
    fprintf('
    %d
    %8.5f
    %8.5f \n',i,x1,x1-x0)
    for i=2:14
    f0=exp(x0)+2^(-x0)+2*cos(x0)-6;
    f1=exp(x1)+2^(-x1)+2*cos(x1)-6;
    x2=x1-(x1-x0)*f1/(f1-f0);
    dist=abs(x2-x1);
    fprintf('
    %d
    %8.5f
    %8.5f \n',i,x2,dist)
    x0=x1;
    x1=x2;
    end
    fprintf('
    ----------------------------\n')
    fprintf('\n
    El valor de X es: %8.5f\n\n',x2)

    El valor de X es: 1.82938

    Para la ecuacin (b):


    clc, clear all
    fprintf('\tPOBLEMA 2.16 PARTE b MTODO DE LA SECANTE\n')
    fprintf('\t=========================================\n')
    x0=1; x1=2; i=0;
    fprintf('
    ----------------------------\n')
    fprintf('
    i
    Xi
    |Xi+1 - Xi|\n')

    fprintf('
    ----------------------------\n')
    fprintf('
    %d
    %8.5f\n',i,x0)
    i=i+1;
    fprintf('
    %d
    %8.5f %8.5f \n',i,x1,x1-x0)
    for i=2:7
    f0=x0*log(x0)-10;
    f1=x1*log(x1)-10;
    x2=x1-(x1-x0)*f1/(f1-f0);
    dist=abs(x2-x1);
    fprintf('
    %d
    %8.5f %8.5f \n',i,x2,dist)
    x0=x1;
    x1=x2;
    end
    fprintf('
    ----------------------------\n')
    fprintf('\n
    El valor de X es: %8.5f\n\n',x2)

    El valor de X es: 5.72893


    Para la ecuacin (c):
    clc, clear all
    fprintf('\tPOBLEMA 2.16 PARTE c MTODO DE LA SECANTE\n')
    fprintf('\t=========================================\n')
    x0=0; x1=1; i=0;
    fprintf('
    ----------------------------\n')
    fprintf('
    i
    Xi
    |Xi+1 - Xi|\n')
    fprintf('
    ----------------------------\n')
    fprintf('
    %d
    %8.5f\n',i,x0)
    i=i+1;
    fprintf('
    %d
    %8.5f
    %8.5f \n',i,x1,x1-x0)
    for i=2:6
    f0=exp(x0)+x0^3+2*x0^2+10*x0-20;
    f1=exp(x1)+x1^3+2*x1^2+10*x1-20;
    x2=x1-(x1-x0)*f1/(f1-f0);
    dist=abs(x2-x1);
    fprintf('
    %d
    %8.5f
    %8.5f \n',i,x2,dist)
    x0=x1;
    x1=x2;
    end
    fprintf('
    ----------------------------\n')
    fprintf('\n
    El valor de X es: %8.5f\n\n',x2)

    El valor de X es: 1.20320


    Para la ecuacin (d):
    clc, clear all
    fprintf('\tPOBLEMA 2.16 PARTE d MTODO DE LA SECANTE\n')
    fprintf('\t=========================================\n')
    x0=2; x1=3; i=0;
    fprintf('
    ----------------------------\n')
    fprintf('
    i
    Xi
    |Xi+1 - Xi|\n')

    fprintf('
    ----------------------------\n')
    fprintf('
    %d
    %8.5f\n',i,x0)
    i=i+1;
    fprintf('
    %d
    %8.5f
    %8.5f \n',i,x1,x1-x0)
    for i=2:9
    f0=sin(x0)-csc(x0)+1;
    f1=sin(x1)-csc(x1)+1;
    x2=x1-(x1-x0)*f1/(f1-f0);
    dist=abs(x2-x1);
    fprintf('
    %d
    %8.5f
    %8.5f \n',i,x2,dist)
    x0=x1;
    x1=x2;
    end
    fprintf('
    ----------------------------\n')
    fprintf('\n
    El valor de X es: %8.5f\n\n',x2)

    El valor de X es: 2.47535

    Ejercicio 2.6.6 del libro de Ledanois


    Una reaccin qumica compleja produce tres componentes B, C y D a partir de A.
    Todas las reacciones son de primer orden (Ki en s -1). Si entran (FA0) 10 Kmol/s de
    A, establecer las concentraciones de cada uno de los componentes en la salida

    del reactor en condiciones estacionarias (V=20L; V0=20L/s; K 1=0,05; K2=0,04;


    K3=0,05; K4=0,02; K5=K6=0,01).
    V[(K1 + K4)CA K6CD] = FA0 V0CA
    20[(0.05 + 0.02)CA 0.01CD] = 10 20CA
    20(0.07CA 0.01CD) = 10 20CA
    1.4CA 0.2CD = 10 20CA
    21.4CA 0.2CD = 10

    (I)

    V[(K2 + K4)CB K1CA] = -V0CB


    20[(0.04 + 0.02)CB 0.05CA] = -20CB
    20(0.06CB 0.05CA) = -20CB
    1.2CB CA = -20CB
    21.2CB CA = 0
    -CA + 21.2CB = 0

    (II)

    V[K3CC K2CB K5CA] = -V0CC


    20[0.05CC 0.04CB 0.01CA] = -20CC
    CC 0.8CB 0.2CA = -20CC
    21CC 0.8CB 0.2CA = 0
    -0.2CA 0.8CB + 21CC = 0 (III)
    V[K6CD K4CB K3CC] = -V0CD
    20[0.01CD 0.02CB 0.05CC] = -20CD
    0.2CD 0.4CB CC = -20CD
    -0.4CB CC + 20.2CD = 0

    (IV)

    21.4
    0 0 0.2
    1 21.2 0
    0
    0.2 0.8 21
    0
    0
    0.4 1 20.2

    ][ ] [ ]
    CA
    10
    CB = 0
    CC
    0
    CD
    0

    PROGRAMACIN EN MATLAB:
    function [x,k]=problema_2_6_6_gausseidel(a,b,x,e,m)
    n=length(x);
    for k=1:m
    t=x;
    for i=1:n
    s=a(i,1:i-1)*x(1:i-1)+a(i,i+1:n)*x(i+1:n);
    x(i)=(b(i)-s)/a(i,i);
    end
    if norm((x-t),inf)<e
    return
    end
    end
    x=[];
    k=m;

    Las concentraciones halladas por el mtodo de Gauss-Seidel


    CA = 0.4673
    CB = 0.0220
    CC = 0.0053
    CD = 0.0007
    k = 3 iteraciones

    Ejercicio 3.6.8 del libro de Ledanois


    La ecuacin de estado de Peng-Robinson (PR) es muy utilizada para la prediccin
    del equilibrio lquido-vapor. Relaciona el volumen especfico v (m 3/Kg) con la
    temperatura T (K) y la presin P (KPa), a travs de la relacin:
    P=

    RT
    a

    V b V ( V +b )+ b(V b)

    Donde R es la constante de los gases (KJ/Kg.K), a=0.45724 (R 2Tc2/Pc)alfa,


    b=0.07780(RTc/Pc), alfa=(1+K(1-Tr0.5))2 y K=0.37464 + 1.54226w 0.26992w 2. Esta
    ecuacin puede ser reescrita en trmino del factor de compresibilidad Z(=PV/RT)
    como sigue:
    Z3 (1-B)Z2 + (A 3B2 2B)Z (AB B2 B3) = 0
    Donde A=Ap/R2T2, B=bP/RT. Resuelva esta ecuacin cbica utilizando el mtodo
    de Newton-Raphson. Se harn los clculos para el metano cuyos datos son:
    R=0.51835 KJ/Kg.K, Tc=191.04K, Pc=4.64MPa w=0.008(factor acntrico). Calcular
    el factor de compresibilidad a 37.73C y 6.804 atm. Se recuerda que la
    temperatura reducida Tr es la relacin entre la temperatura absoluta y la
    temperatura crtica.
    P=

    RT
    a

    V b V ( V +b )+ b(V b)

    Multiplicando por (V b) a ambos lados nos da:

    ( V b ) P=

    RT (V b)
    a(V b)

    (V b)
    V (V + b ) +b(V b)

    dividimos por P a ambos lados y la ecuacin queda de la siguiente manera:

    ( V b )=

    a(V b)
    RT

    P P [V ( V +b )+ b (V b )]

    por ultimo expresamos la ecuacin en funcin del volumen:


    V=

    a ( V b )
    RT

    +b
    P P [ V ( V +b ) + b (V b ) ]

    Este volumen lo usaremos para dar un valor inicial a Z


    Z=PV/RT
    PROGRAMACIN EN MATLAB:
    clc, clear all
    disp('PROBLEMA 3.6.8')
    disp('Obtencin del volumen especfico')
    disp(' ')
    disp('Solucin por el Mtodo de Punto Fijo')
    R=0.51835;
    %constante de gas ideal para el metano en, KJ/Kg.K
    T=310.73;
    %temperatura absoluta en, K
    P=689.42;
    %presin absoluta en, KPa
    Tc=191.04;
    %temperatura crtica en, K
    Pc=4640;
    %presin crtica en, KPa
    Tr=T/Tc;
    %temperatura reducida
    w=0.008;
    %factor acntrico
    K=0.37464+1.54226*w-0.26992*w^2;
    alfa=(1+K*(1-Tr^0.5))^2;
    a=0.45724*(R^2*Tc^2/Pc)*alfa;
    b=0.07780*(R*Tc/Pc);
    A=a*P/(R^2*T^2);
    B=b*P/(R*T);
    V0=1;
    V=R*T/P-a*(V0-b)/(P*(V0*(V0+b)+b*(V0-b)))+b;
    error=0.00001;
    fprintf('\nTolerancia = %7.5f\n\n', error)
    iter=1;
    disp(' Iter
    V ')
    disp('--------------------')
    fprintf(' %2d
    %10.5f\n', iter, V)
    while abs(V-V0) > error
    iter=iter+1;
    V0=V;
    V=R*T/P-a*(V0-b)/(P*(V0*(V0+b)+b*(V0-b)))+b;
    fprintf(' %2d
    %10.5f\n', iter, V)
    end
    disp('--------------------')
    fprintf('\nEl Volumen especfico es: %8.5f\n\n',V)
    fprintf('\tCALCULO DEL FACTOR DE COMPRESIBILIDAD APLICANDO EL METODO DE
    NEWTON-RAPHSON\n')
    fprintf('\t==============================================================
    =============\n')
    i=0;
    Z=1;
    dg=abs(1-((3*Z^2-2*(1-B)*Z+(A-3*B^2-2*B))^2-(Z^3-(1-B)*Z^2+(A-3*B^22*B)*Z-(A*B-B^2-B^3))*(6*Z-2*(1-B)))/...
    (3*Z^2-2*(1-B)*Z+(A-3*B^2-2*B))^2);
    Z0=P*V/(R*T);
    fprintf('-----------------------------------------------\n')
    fprintf(' i
    Zi
    |Zi+1 - Zi|
    |dg(Zi)|\n')
    fprintf('-----------------------------------------------\n')
    fprintf(' %d
    %8.5f
    %8.2e\n',i,Z,dg)
    for i=1:2

    f=Z0^3-(1-B)*Z0^2+(A-3*B^2-2*B)*Z0-(A*B-B^2-B^3);
    df=3*Z0^2-2*(1-B)*Z0+(A-3*B^2-2*B);
    Z=Z0-f/df;
    dist=abs(Z-Z0);
    dg=abs(1-((3*Z^2-2*(1-B)*Z+(A-3*B^2-2*B))^2-(Z^3-(1-B)*Z^2+(A-3*B^22*B)*Z-(A*B-B^2-B^3))*(6*Z-2*(1-B)))/...
    (3*Z^2-2*(1-B)*Z+(A-3*B^2-2*B))^2);
    fprintf(' %d
    %8.5f
    %8.5f
    %8.2e\n',i,Z,dist,dg)
    Z0=Z;
    end
    disp('---------------------------------------------------')
    fprintf('\nEl factor de compresibilidad (Z) es: %8.5f\n\n',Z)

    RESPUESTA: El factor de compresibilidad (Z) es: 0.98677

    También podría gustarte