Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 249 vistas

    Asignación 1

    Cargado por

    Emerson Sivincha
    Este documento presenta la resolución de dos problemas de modelado matemático mediante métodos numéricos. El primer problema resuelve un sistema de ecuaciones diferenciales ordinarias que modela la concentración de un reactante en tres reactores en serie, usando el método de Runge-Kutta. El segundo problema resuelve una ecuación en derivadas parciales que modela la transferencia de calor en una dimensión, usando los métodos explícito e implícito de diferencias finitas.

    Copyright:

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

    Asignación 1

    Cargado por

    Emerson Sivincha
    249 vistas12 páginas
    Este documento presenta la resolución de dos problemas de modelado matemático mediante métodos numéricos. El primer problema resuelve un sistema de ecuaciones diferenciales ordinarias que modela la concentración de un reactante en tres reactores en serie, usando el método de Runge-Kutta. El segundo problema resuelve una ecuación en derivadas parciales que modela la transferencia de calor en una dimensión, usando los métodos explícito e implícito de diferencias finitas.

    Título original

    ASIGNACIÓN 1

    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 dos problemas de modelado matemático mediante métodos numéricos. El primer problema resuelve un sistema de ecuaciones diferenciales ordinarias que modela la concentración de un reactante en tres reactores en serie, usando el método de Runge-Kutta. El segundo problema resuelve una ecuación en derivadas parciales que modela la transferencia de calor en una dimensión, usando los métodos explícito e implícito de diferencias finitas.

    Copyright:

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

    Asignación 1

    Cargado por

    Emerson Sivincha
    Este documento presenta la resolución de dos problemas de modelado matemático mediante métodos numéricos. El primer problema resuelve un sistema de ecuaciones diferenciales ordinarias que modela la concentración de un reactante en tres reactores en serie, usando el método de Runge-Kutta. El segundo problema resuelve una ecuación en derivadas parciales que modela la transferencia de calor en una dimensión, usando los métodos explícito e implícito de diferencias finitas.

    Copyright:

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

    UNIVERSIDAD NACIONAL DE SAN AGUSTIN

    FACULTAD DE INGENIERÍA PROCESOS


    Escuela Profesional de Ingeniería Metalúrgica

    CURSO: MODELAMIENTO Y SIMULACION

    DOCENTE: Ing. FEDERICO CALISAYA MAMANI

    INTEGRANTE:

    SIVINCHA PAUCAR OSCAR EMERSON

    AREQUIPA
    2021
    ASIGNACIÓN N°1

    RESOLVER POR EL METODO DE NEWTON-RAPHSON EL SIGUIENTE SISTEMA


    DE ECUACIONES NO LINEALES
    f ( x , y , z )=3 x−cos ( yz ) −0.5=0
    f ( x , y , z )=x 2−625 y 2=0
    10 π −3
    f ( x , y , z )=e− xy +20 z+ =0
    3
    : Obtener la matriz de derivadas parciales

    ∂f 2 ∂f1 ∂f 1
    =3 =zsen ( yz ) = ysen ( yz )
    ∂x ∂y ∂z

    ∂f 2 ∂f 2 ∂f 2
    =2 x =−1250 y =0
    ∂x ∂y ∂z
    ∂f 3 ∂f 3 ∂f 3
    =− y e−xy =−x e−xy =20
    ∂x ∂y ∂z

    Obtener la matriz jacobiana formado por h1, h2 y h3 que son incognitas de esta matriz.
    B

    3 zsen ( yz ) ysen ( yz ) h 1 −f 1
    2x −1250 y 0 h 2=−f 2
    A X
    − y e−xy −x e−xy 20 h 3 −f 3

    Expresar la matriz en forma [x]=[A^(-1)][B] para encontrar los valores de [hi]


    siendo
    Este último lo correspondiente a pequeños aumentos para valores de
    xi+1=xi+h

    Primero, tomando x0=[25,1,1] se obtiene la matriz jacobiana

    Matriz inversa de [A]

    0.32963496 0.0002219 0.01386891


    0.0131854 -0.00079112 0.00055476
    4.578E-13 -1.358E-14 0.05

    Los valores de hi son:


    Siendo el valor nuevo
    Para X
    x1+1=x0+h

    x1+1=25 -24.2336895=0.76631051

    Realizando el mismo procedimiento descrito anteriormente se obtiene el siguiente


    cuadro.
    FINALMENTE LA SOLUCION DE LA MATRIZ AMPLIADA DA COMO RESULTADO:

    Resolver el siguiente sistema de ecuaciones por el método de punto fijo


    multivariable método interactivo de Jacobo y gaus seidal

    F1(x,y)= x2-10+y2+8=0
    F2(x,y)= xy2+x+10y+8=0

    Solución:
    Xk+1 = 0.1(xk)2+0.1(yk)2+0.8
    Yk+1 = 0.1xk+1(yk)2 +0.1 xk+1 +0.8
    Aplicando los criterios de convergencia se requiere derivar parcialmente:

    Evaluando para x0 = 0, y0 = 0
    1er iteración:
    x1= 0.8
    y1= 0.1(0.8)+0.8=0.88
    2da iteración:
    x2 = 0.1(0.8)2 +0.1(0.88)2 + 0.8 = 0.9414
    y2 = 0.1(0.9414)(0.88)2+0.1(0.9414)+0.8 = 0.9670
    3er iteración:
    x3= 0.1(0.9414)2 +0.1(0.9670)2 + 0.8 =0.9821
    y3= 0.1(0.9821)(0.9670)2+0.1(0.9821)+0.8 =0.9900
    Los resultados del proceso iterativo son:

    K Xk yk
    0 0.0000 0.0000
    1 0.8000 0.8800
    2 0.9414 0.9670
    3 0.9821 0.9900
    4 0.9944 0.9969
    * * *
    11 1.0000 1.0000
    1. RESOLVER

    x , =2 x + 4 y
    y ´ ´ =−x +6 y

    POR EL METODO DE RUNGE-KUTTA APROXIME A X(0.6), Y(0.6), COMPARE LOS


    RESULTADOS OBTENIDOS CON h1 Y h2.
    1
    zi+1=zi+ ( k 1+2 k 2+2 k 3+k 4 )∗h
    6
    k 1=f ( t ; z )

    1 1
    (
    k 2=f t+ h ; z+ k 1h
    2 2 )
    1 1
    k 3=f ( t+ h ; z + k 2 h)
    2 2

    k 4=f ( t+ h; z +k 3 h )
    Para la solucion de este problema se utilizó propiedades de matrices
    2 x +4 y
    f ( t , z )= … … … … … … … … … … ….(1)
    −x+ 6 y
    PRIMERA INTERACCIÓN (h=0.1) = VALOR INICIAL z=(1;1)
    PARA H=0.1

    2 ( 1 )+ 4 ( 1 )
    =6
    k 1=
    ( )()
    −( 1 ) +6 ( 1 ) 5

    1 1
    k 2=f 0+ ∗0.1; 1 + ∗0.1∗ 6 =f 0.05 ; 1.3 = 7.6
    ( 2 1 2 () 5 ( )) ( ( )) ( )
    1.25 6.2

    ( 12 1 12 6.2 ) ( 1.31 ) 6.48


    k 3=f 0+ ∗0.1; ( 1 )+ ∗( 7.6 )∗0.1 =f 0.05 ; ( 1.38 ) =( 8 )

    k 4=f 0+0.1 ; ( 1 ) +( 8 )∗0.1 =f 0.1 ; ( 1.8 ) =( 10.192 )


    ( 1 6.48 ) ( 1.648 ) 8.088
    1
    zi+1=( 1) + ( 6 ) +2 ( 7.6 ) +2 ( 8 ) + ( 10.192) 0.1=( 1.789)
    1 6( 5 6.2 6.48 8.088 ) 1.641
    Se sigue el mismo procedimiento hasta x(0.6) y y(0.6) es decir t=0.6 cuadro de
    resumen de la solución de las ecuaciones

    PARA H=0.2
    El error de truncamiento originado por la naturaleza de las técnicas empleadas para
    aproximar los valores de y es menor a un tamaño de paso de H=0.1 que a H=0.2. Con
    valores de H suficientemente pequeña, los errores normalmente disminuyen, en tanto
    f ´ ´ (x , y )h4
    aumenta el orden. Et =
    2!

    RESUELVA EL SIGUIENTE PROBLEMA DE VALORES EN LA FRONTERA POR EL


    METODO EXPLICITOS E IMPLICITOS DE DIFERENCIAS DIFINITAS

    ∂T ∂2 T
    =α
    ∂t ∂ x2
    T ( x , o )=20 senx
    T ( 0 , t )=100 ° C
    T ( 1 , t ) =50° C
    pi e2
    α =1.0 L=1.0 pie tmax=1 h
    h

    METODO EXPLICITO

    En las esquinas de la malla, debido a que está gobernado por dos condiciones de
    frontera, en este caso se tomó el promedio
    La segunda derivada se aproxima mediante una diferencia dividida finita centrada y
    una diferencia dividida finita hacia adelante sirve para aproximar a la derivada con
    respecto al tiempo.
    0.52 100 99.598 99.431 99.598 100
    0.54 100 99.673 99.538 99.673 100
    0.56 100 99.734 99.624 99.734 100
    0.58 100 99.784 99.695 99.784 100
    0.6 100 99.825 99.752 99.825 100
    0.62 100 99.858 99.799 99.858 100
    0.64 100 99.884 99.836 99.884 100
    0.66 100 99.906 99.867 99.906 100
    0.68 100 99.924 99.892 99.924 100
    0.7 100 99.938 99.912 99.938 100
    0.72 100 99.950 99.929 99.950 100
    0.74 100 99.959 99.942 99.959 100
    0.76 100 99.967 99.953 99.967 100
    0.78 100 99.973 99.962 99.973 100
    0.8 100 99.978 99.969 99.978 100
    0.82 100 99.982 99.975 99.982 100
    0.84 100 99.985 99.979 99.985 100
    0.86 100 99.988 99.983 99.988 100
    0.88 100 99.990 99.986 99.990 100
    0.9 100 99.992 99.989 99.992 100
    0.92 100 99.994 99.991 99.994 100
    0.94 100 99.995 99.993 99.995 100
    0.96 100 99.996 99.994 99.996 100
    0.98 100 99.997 99.995 99.997 100
    1 100 99.997 99.996 99.997 100
    METODO IMPLICITO
    En este método la derivada espacial se aproxima en un nivel de tiempo anterior l-1 y la
    segunda derivada mediante diferencia finita centrado, con ello se consigue obtener un
    conjunto de ecuaciones algebraicas, por lo que el método se reduce a la solucion de
    un sistema de ecuaciones en cada punto en el tiempo.
    El resultado se muestra en la tabla.

    2. Un sistema de 3 reactores en serie es discreta por las siguientes


    EDO

    DCa 1 1
    = x ( Ca0−Ca1 )−hxCa1
    Dt ῖ
    DCa 2 1
    = x ( Ca0−Ca2 ) −hxCa2
    Dt ῖ
    DCa 3 1
    = x ( Ca2−Ca3 ) −hxCa3
    Dt ῖ

    Donde
    Ca= Concentración del reactante A
    T= Tiempo
    t= V/a= 2 minutos tiempo de residencia promedio
    h=Constante de la cinética de reacción 0.5 min-1
    Ca0= Concentración de entrada se mantiene constante en 1.8 kg.mol de
    A/m3

    Resolver para Ca.(T),Ca2(T), Ca3(T) por el método de Runge kutta Con


    AT=0.1 minutos hasta Ca1(2.9). Ca2(2.9),Ca3 (2.9 Con Ca1(0)=0.4 Kg.mol/m3,
    Ca2(0) = 0.2 km.mol/m3 Ca3(0)= 0.1 kg.mol/m3

    dCa
    =0.5 ( 1.8−Ca )−0.5Ca
    dt
    0 Ca1+ Ca2 =0.9
    Ca1 ( 0+ 1 )=0.9

    d Ca 2
    =0.5 ( Ca 1−Ca 2 )−0.5Ca 2
    dt
    dCa2−0.5Ca1 +Ca2=0
    d Ca 3
    =0.5 ( Ca 2−Ca 3 ) −0.5Ca 3
    dt
    dCa3−0.5Ca2 +Ca3=0
    ¿

    d Ca 1
    + Ca1 =0.9
    dt
    t d Ca 1 t t
    e + e Ca1=e 0.9
    dt
    1
    ∫ d Ca1= e t ∫ e t 0.9 dt
    0.9
    Ca1= t (e t +a)
    e

    Ca1=0.9+ k e−t concentracion 1


    ¿

    ( D+1 )3 Ca2=−0.45(D+1)
    d 2 Ca 2 2 d Ca 2
    + +Ca 2=−0.45
    dt2 dt
    r 2 +2 r +r=0
    Ca2= A 0 C−t + A 1 Tc−t
    C−t T e−t
    W=
    −e(
    −t −t
    e −T e −t =e
    −2 t
    )
    −T e−2 t +T e−2 t =e−2 t

    0 T e−t
    U ( t )= ( −t
    0.45 e −T e −t =∫ 0.45 e)
    −2t

    ∫ 0.45 Te−t dt =0.45 ¿ ¿


    ∫ 0.45 Te−t dt =0.45 Tet−0.45 e t
    e−t 0.45 t

    U ( t )=
    ( −e−t 0.45 )
    e−2 t
    U ( t )=∫ −0.45 e t dt=−0.45 e t
    Ca2=0.45 T + 0.45−0.45 t
    Ca2=−0.45+ A 0 et + A 1 T e−t Concentracion 2

    D+1 0 0.9
    (
    ( D+1 )3 Ca3 = −0.5 D+1 0
    0 −0.5 1 )
    ( D+1 )3 Ca3 =0.225( D+1)
    d 2 Ca 3 2 d Ca 3
    + +Ca 3 =0.225
    dt2 dt
    r 2 +2 r +r=0
    Ca2= A1 e−t + A 2 T e−t
    Ca2=0.225+ A 1 e−t + A 2 T e−t Concentracion 3
    Ca1=0.9+ k e−t
    Ca2=−0.45+ A 0 et + A 1 T e−t Concentracion 2
    Ca2 ( 0 )=0.2
    0.2=−0.45+ A0
    A0 =−0.65
    Ca3 ( 0 )=0.1
    0.1=0.225− A1
    A1=0.325
    Ca0=0.4
    0.4=0.9+ k
    k =0.4−0.9=−0.5
    Ca1=0.9−0.5 e−t

    Ca1 ( 0.1 ) =0.9−0.5 e−0.1

    0.44 kg
    Ca1 ( 0.1 ) =
    m3
    8.187 kg
    Ca1 ( 2.9 ) =
    m3
    7.75 kg
    Ca1 ( 2.9 ) −Ca1 ( 0.1 )=
    m3

    También podría gustarte