Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 15 vistas

    Actividad 2

    Cargado por

    jorge.ramirez9331
    Numerical methods report

    Copyright:

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

    Actividad 2

    Cargado por

    jorge.ramirez9331
    15 vistas8 páginas
    Numerical methods report

    Título original

    Actividad_2

    Derechos de autor

    © © All Rights Reserved

    Formatos disponibles

    PDF, TXT o lea en línea desde Scribd

    ¿Le pareció útil este documento?

    ¿Este contenido es inapropiado?

    Numerical methods report

    Copyright:

    © All Rights Reserved
    Descargue como PDF, TXT o lea en línea desde Scribd
    Descargar como pdf o txt
    15 vistas8 páginas

    Actividad 2

    Cargado por

    jorge.ramirez9331
    Numerical methods report

    Copyright:

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

    UNIVERSIDAD DE GUADALAJARA

    CENTRO UNIVERSITARIO DE CIENCIAS EXACTAS E INGENIERÍAS

    Cómputo Cientı́fico II

    Generalidades

    Alumnos: Jorge Ramı́rez López, (222793314)


    Profesor: Dr. Diego Armando Pantoja González

    6 de febrero de 2024
    Índice
    1. Introducción II

    2. Metodologı́a III

    3. Desarrollo III

    4. Resultados III

    5. Conclusiones III

    6. Referencias IV

    Apéndices V

    A. Código fuente V

    B. Imágenes VI

    I
    Resumen
    En la práctica de esta semana, se empleó el Método de Euler para resolver ecuaciones diferen-
    ciales ordinarias, para modelar la transferencia de calor en un cuerpo de agua. Posteriormente, se
    generó un conjunto de gráficas que mostraban tanto la solución exacta como la solución numérica
    del modelo, junto con la radiación solar que afectaba al cuerpo de agua.

    1. Introducción Con el punto Po se puede evaluar la primera de-


    rivada de F (x), la curva solución del PVI, en ese
    Se deduce del análisis de un sistema la si- punto; a saber
    guiente ecuación diferencial ordinaria
    dy
    F ′ (x) = . (4)
    dT dx Po =f (xo ,yo )
    + q(T − Ti ) = aJs sin ωt, (1)
    dt
    Con esta información se procede a trazar una
    Q
    en donde las expresiones de q y a son q = V , recta, aquella que pasa por Po y de pendiente
    a = VAρc , y V es el volumen, ρ la densidad, c el f (xo , yo ). Esta recta aproxima F(x) en una ve-
    calor especı́fico, T la temperatura, Q es el flu- cindad de xo . Tómese la recta como remplazo de
    jo, A es el área, y Js es la radiación solar. La F (x) y localı́cese en ella (la recta) el valor de y
    cual modela como la temperatura T del sistema correspondiente a x1 . Entonces,
    cambia con el tiempo debido al flujo del calor
    y1 − yo
    q(T − Ti ) y la radiación solar aJs sin ωt. La ra- = f (xo , yo ). (5)
    diación solar afecta la temperatura del sistema, x1 − xo
    y la transferencia de calor también contribuye a Se resuelve para y1
    esta variación.
    y1 = yo + (x1 − xo )f (xo , yo ) = yo + hf (xo , yo ).
    La solución de esta ecuación permite enton- (6)
    ces predecir como evoluciona la temperatura del Aplicando el método de Euler a nuestra ecua-
    sistema en función del tiempo y de los paráme- ción de transferencia de calor (Ec. 1), obtenemos
    tros dados. la solución numérica

    Es posible encontrarla solución analı́tica de la Tnum (k+1) = Tnum (k)+dt[aJs sin ωt(k)−q(Tnum (k)−Ti )]
    ecuación (Ec. 1), la cual es (7)
    donde t va en intervalos de dt dı́as.
    aJs
    T (t) = Ti +(To −Ti )e−qt + 2 (ωe−qt +ω cos ωt+q sin ωt).
    ω + q2
    (2) Es posible componer un script en software de
    El método de Euler es el más simple de los programación para resolver las operaciones nece-
    métodos numéricos para resolver un problema de sarias y graficar la función, tanto analı́tica, como
    valor inicial (Hurtado y Sánchez, 2012, pp. 539- numérica.
    541) como el que se planteó anteriormente. Este
    consiste en dividir el intervalo que va de xo a xf
    en n subintervalos de ancho h; o sea

    xf − xo
    h= (3)
    n
    de manera que se obtiene un conjunto discreto
    de (n + 1) puntos del intervalo de interés [xo , xf ].

    II
    2. Metodologı́a 3. Desarrollo
    Se siguieron los siguientes pasos para cons- Para la simulación presentada, se utilizaron
    truir el script en Matlab: los siguientes parámetros:

    1. Se definen las constantes: Variable Valor


    A 2.5 × 106 m2
    Se escribieron los parámetros como ρ 1000 kg/m3
    constantes. V 5.0 × 108 m3
    C 4182 J/kg◦ C
    Se corrigieron las unidades.
    Q 80 m3 /s
    Se definió el vector de tiempo. Ti 10◦ C
    Js 290 J/s m2
    2. Generar los vectores de las solucio- To 15◦ C
    nes: 2π
    ω 2×365 rad/d
    dt 60 dı́as
    Para la solución exacta, se definió un
    vector de igual dimensión que el vec- donde Q y Js tienen que multiplicarse por
    tor de tiempo, a partir de la Ec. 2. 86400 para corregir sus unidades (conversión de
    segundos a dı́as).
    Para la solución numérica, se cons-
    Cabe destacar que el intervalo de tiempo es dt =
    truyó un bucle for tal que se completó
    60 dı́as, un valor relativamente grande.
    un vector de igual dimensión que el
    Se graficaron en un plano las soluciones analı́tica
    vector del tiempo, a partir de la Ec.
    y numérica, y en otro, la radiación solar. Todo
    7.
    en función del tiempo y exportado en una misma
    figura.
    3. Graficación:

    Para cada gráfica se definió un eje 4. Resultados


    (“Create Cartesian Axes - MATLAB
    Axes - MathWorks América Latina”, La gráfica obtenida, en donde se pueden apre-
    s.f.) de tal manera que estuvieran or- ciar las curvas de la solución exacta y numérica
    ganizados en dos columnas. se puede consultar en el Apéndice B.
    El error relativo porcentual encontrado fue de
    Se exportó la figura como archivo png
    2.7573 %.
    para la presentación de los resultados.

    4. Cálculo del error: 5. Conclusiones


    Se calcula el error relativo porcentual.Se concluye que el método de Euler resul-
    ta muy sencillo de utilizar y es una muy buena
    De esta manera se construyó el script dispo- aproximación, pues da un error pequeño, incluso
    nible en el Apéndice A. cuando el intervalo de tiempo fue muy grande.

    III
    6. Referencias

    Referencias
    Create Cartesian Axes - MATLAB Axes - MathWorks América Latina. (s.f.). https://la.mathworks.
    com/help/matlab/ref/axes.html
    Hurtado, A. N., & Sánchez, F. C. D. (2012, 1 de enero). Métodos numéricos aplicados a la ingenierı́a
    (4.a ed.).

    IV
    Apéndices
    A. Código fuente

    Código MATLAB.

    clc
    clear a l l

    % Constantes
    A = 2 . 5 E6 ;
    V = 5E8 ;
    Q = 80 ∗ 8 6 4 0 0 ;
    J s = 290 ∗ 8 6 4 0 0 ;
    omega = ( 2 ∗ pi ) / ( 2 ∗ 3 6 5 ) ;
    rho = 1 0 0 0 ;
    C = 4182;
    Ti = 1 0 ;
    To = 1 5 ;

    a = A/ (V∗ rho ∗C ) ;
    q = Q/V;

    dt = 6 0 ;
    t = 0 : dt : 3 6 5 ;

    % S o l u c i n Exacta
    T e x a c t a = Ti + (To − Ti ) ∗ exp(−q . ∗ t ) +
    ( ( a ∗ J s ) / ( omega ˆ2 + q ˆ 2 ) ) ∗ ( omega ∗ exp(−q . ∗ t )
    − omega ∗ cos ( omega . ∗ t ) + q ∗ sin ( omega . ∗ t ) ) ;

    % S o l u c i n Num rica
    T numerica ( 1 ) = To ;
    for i = 1 : 6
    T numerica ( i + 1 ) = T numerica ( i ) +
    dt . ∗ ( a . ∗ J s . ∗ sin ( omega . ∗ t ( i ) ) − q . ∗ ( T numerica ( i ) − Ti ) ) ;
    end

    V
    Código MATLAB

    % Figura y p l o t t i n g
    figure ( 1 ) , c l f (1)

    ax1 = axes ( ’ P o s i t i o n ’ , [ 0 . 1 0 . 5 5 0 . 8 0 . 3 5 ] ) ;

    hold on

    plot ( ax1 , t , T exacta , Marker = ” diamond ” , C o l o r = ’#0d98ba ’ )


    plot ( ax1 , t , T numerica , Marker = ” o ” , C o l o r = ’#cb4154 ’ )
    legend ( ’ S o l u c i n Exacta ’ , ’ S o l u c i n N u m r i c a ’ )
    xlabel ( ’ Tiempo [ d a s ] ’ )
    ylabel ( ’ Temperatura [ C ] ’ )
    t i t l e ( ’ F l u j o s de c a l o r ’ )

    hold o f f

    grid on , box on , axis t i g h t

    ax2 = axes ( ’ P o s i t i o n ’ , [ 0 . 1 0 . 1 0 . 8 0 . 3 5 ] ) ;
    plot ( ax2 , t , sin ( omega . ∗ t ) )
    xlabel ( ’ Tiempo [ d a s ] ’ )
    ylabel ( ’ R a d i a c i n S o l a r [W/mˆ 2 ] ’ )

    grid on , box on , axis t i g h t

    s a v e a s ( 1 , ’ f i g u r a . png ’ )

    ERP = abs (norm( T numerica ) − norm( T e x a c t a ) ) /norm( T e x a c t a ) ∗ 1 0 0 ;

    B. Imágenes

    VI
    VII

    También podría gustarte