Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 19 vistas

    Metodos - Aplicaciones

    Cargado por

    Alexandra Huertas Rodriguez
    Este documento presenta tres ejercicios de aplicación numérica relacionados con temas de cálculo. El primer ejercicio resuelve numéricamente la distribución de temperatura en una barra delgada utilizando diferencias finitas. El segundo calcula la cantidad de trabajo realizado al mover una partícula entre dos puntos usando métodos de integración numérica. El tercer ejercicio aproxima valores de una función tabulada usando interpolación polinómica de Lagrange.

    Copyright:

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

    Metodos - Aplicaciones

    Cargado por

    Alexandra Huertas Rodriguez
    19 vistas6 páginas
    Este documento presenta tres ejercicios de aplicación numérica relacionados con temas de cálculo. El primer ejercicio resuelve numéricamente la distribución de temperatura en una barra delgada utilizando diferencias finitas. El segundo calcula la cantidad de trabajo realizado al mover una partícula entre dos puntos usando métodos de integración numérica. El tercer ejercicio aproxima valores de una función tabulada usando interpolación polinómica de Lagrange.

    Descripción original:

    Solución Numérica de una Ecuación Diferencial de Calor utilizando el Método Numérico de las Diferencias Finitas

    Título original

    METODOS - APLICACIONES

    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 tres ejercicios de aplicación numérica relacionados con temas de cálculo. El primer ejercicio resuelve numéricamente la distribución de temperatura en una barra delgada utilizando diferencias finitas. El segundo calcula la cantidad de trabajo realizado al mover una partícula entre dos puntos usando métodos de integración numérica. El tercer ejercicio aproxima valores de una función tabulada usando interpolación polinómica de Lagrange.

    Copyright:

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

    Metodos - Aplicaciones

    Cargado por

    Alexandra Huertas Rodriguez
    Este documento presenta tres ejercicios de aplicación numérica relacionados con temas de cálculo. El primer ejercicio resuelve numéricamente la distribución de temperatura en una barra delgada utilizando diferencias finitas. El segundo calcula la cantidad de trabajo realizado al mover una partícula entre dos puntos usando métodos de integración numérica. El tercer ejercicio aproxima valores de una función tabulada usando interpolación polinómica de Lagrange.

    Copyright:

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

    I.

    PROBLEMA APLICADO

    Determine numéricamente la distribución de temperatura en un barra delgada de 2.5cm


    de longitud, la cual se encuentra inicialmente a una temperatura de 0ºC, los extremos de
    la barra se someten a temperaturas de T(0,t)=100ºC en el lado derecho y T(2,t)=50ºC en
    el lado izquierdo constantes con el tiempo.
    Asuma un ∆ x=¿0.5cm y ∆ t=0.5 s

    ∂T ∂2 T
     Ecuación de Calor =k donde T =T ( x , t)
    ∂t ∂ x2

     Con Diferencias Finitas

    ∂T ' T i+1−T i
    =f ( x i )=
    ∂t h
    2
    ∂T '' T i+1 −2T i +T i−1
    2
    =f ( xi ) =
    ∂x h2
    En ( x i ,t j )
    T i , j+1 −T i , j T i+ 1, j−2T i , j+T i−1 , j
    =k
    ∆t ∆ x2
    T i=T i + λ (T i+1−2 T i +T i−1)

    Δt
    λ=K
    ( ∆ x )2
    0.5
    λ=0.075 2
    =0.15
    (0.5)

     Ecuación

    T i=T i + 0.15(T i +1−2 T i +T i−1)

    Condiciones de frontera
    T(0,t)=100ºC
    T(2,t)=50ºC
    Tenemos:
    1 :T 0 =100; x 0=0
    2 :T 1=0 ; x 1=0.5
    3 :T 2 =0 ; x2 =1
    4 :T 3 =0 ; x3 =1.5
    5 :T 4=0 ; x 4 =2
    6 :T 5=50 ; x 4=2.5
     i=1
    T 1+ 0.015 ( T 2−2 T 1+T 0 )=T 1
    reemplazando
    T 1+ 0.15T 2−0.3 T 1+ 0.15 ( 100 )=T 1

     i=2
    T 2+ 0.015 ( T 3−2 T 2+T 1 )=T 2
    reemplazando
    T 2+ 0.15T 3−0.3 T 2+ 0.15T 1=T 2

     i=3
    T 3+ 0.015 ( T 4 −2 T 3 +T 2 )=T 2
    reemplazando
    T 3+ 0.15T 4 −0.3 T 3 +0.15 T 2=T 2

     i=4
    T 4 +0.015 ( T 5−2 T 4 + T 3 ) =T 3
    reemplazando
    T 4 +0.15(50)−0.3T 4 + 0.15T 3=T 3

     Sistema de ecuaciones

    0.7 T 1+ 0.15T 2 +15=T 1


    0.7 T 2+ 0.15T 3 +0.15 T 1 =T 2
    0.7 T 3+ 0.15 T 2 +0.15 T 4=T 3
    0.7 T 4 +0.15 T 3+7.5=T 4

    Con las condiciones iniciales:


    T 1=T 2 =T 3 =T 4 =0
    APLCANDO EN MATLAB

    II. EJERCICIO APLICATIVO DE INTEGRACIÓN


    Teniendo en cuanta que cuando una partícula está a x pies, actúa sobre ella una fuerza
    igual a ex/2 +2x-1 libras. ¿Cuánto trabajo se efectúa al moverla de x=1 a x=3?

    Solución:
    Sabemos
    Formula de trabajo: W=F*d
    b

    Adaptado a integración W =∫ F x . dx
    a

    Gráfico y= ex/2 +2x-1

    Integración numérica:
    Método Trapecio
    Asumiendo n=4

    Ecuación:
    b
    b−a
    W =∫ f x . dx= f ( x ) +2∗f ( x 1 ) +2∗f ( x 2) + 2∗f ( x 3 ) + …+f ( x n ) ]
    a 2∗n [ 0

    Reemplazando:
    b−a
    ∆ x= ; x0 =1; x 1=x 0 +∆ x ; x2 =x1 + ∆ x ; x3 =x 2+ ∆ x ; x 4= x3 + ∆ x
    n
    3−1
    ∆ x= =0.5 ; x 0 =1; x 1=1.5 ; x 2=2 ; x 3=2.5; x 4 =3
    4
    3 x
    2 3−1
    ∫e +2 x−1. dx=
    2∗4
    [ f ( 1 )+ 2∗f ( 1.5 )+2∗f ( 2 ) +2∗f ( 2.5 ) + f ( 3 ) ]
    1

    1
    ¿ [ 2.64872+ 2 ( 4.117 ) +2 ( 5.71828 ) +2 ( 7.49034 ) +9.48169 ]
    4
    1
    ¿ [ 38.54765 ] =9.63691
    4
    Método Simpson 1/3:
    Asumiendo n=4
    b

    ∫ f x . dx= ∆3x [ f ( x 0 ) + 4∗f ( x 1 ) +2∗f ( x 2 ) +4∗f ( x 3 )+ 2∗f ( x 4 ) +…+ f ( x n ) ]


    a

    3−1
    ∆ x= =0.5 ; x 0 =1; x 1=1.5 ; x 2=2 ; x 3=2.5; x 4 =3
    4
    3 x
    0.5
    ∫ e 2 +2 x−1. dx= 3
    [ f (1 ) + 4∗f ( 1.5 ) +2∗f ( 2 ) +4∗f ( 2.5 ) +f ( 3 ) ]
    1

    1
    ¿ [ 2.64872+4 ( 4.117 )+ 2 ( 5.71828 ) + 4 ( 7.49034 ) + 9.48169 ]
    6
    1
    ¿ [ 69.99633 ] =11.66606
    6
    Método Simpson 3/8:
    Asumiendo n=3
    b

    ∫ f x . dx= 3∗∆ x
    f ( x ) +3∗f ( x 1 ) +3∗f ( x 2 ) +f ( x 3 ) ]
    8 [ 0
    a

    3−1 2 5 7
    ∆ x= = ; x0 =1; x 1= ; x 2= ; x 3=3
    3 3 3 3
    2
    3 x
    3 5 7
    2
    [ ()
    ∫ e +2 x−1. dx= 3 f ( 1 ) +3∗f 3 +3∗f 3 + f ( 3 )
    1
    () ]
    3 x
    2
    ∫ e 2 +2 x−1. dx= 9
    [ 2.64872+3(4.63431)+3 (6.87794)+9.48169 ]
    1

    2
    ¿ [ 46.66716 ] =10.37048
    9
    Valor por integración: 11.66593559927587, el método de Simpson 1/3 se aproxima a la
    respuesta con 11.66606 lb/ft
    Respuesta: Mover la partícula desde x=1 hasta x=3 requiere un trabajo de 11.6
    libras/pie.

    III. EJERCICIO DE APLICACIÓN DE INTERPOLACION


    Supongamos que tenemos la función tabular sin una función conocida
    i 0 1 2 3
    F(Xi) -3 0 5 7
    Xi 0 1 3 6

    a) Determinar la aproximación polinomial de LaGrange usando todos los puntos


    (función f (x)).
    b) Determinar el valor aproximado de f (x) para x = 1.8
    Solución:
    Como se sabe la tabla contiene 4 puntos, por ende se deduce a que existe un polinomio
    de tercer orden.
    P3 ( x )=L0 ( x )f (x 0 )+L1 ( x ) f ( x1 )+ L2 ( x ) f ( x2 )+L3 ( x ) f ( x3 )
    P3 ( x )=L0 (x )(−3 )+L1 ( x )(0 )+L2 ( x )(5 )+L3 ( x )(7 )
    ( x−x 1 )(x −x 2 )( x−x 3 ) (x −1)( x−3 )(x −6 ) ( x−1 )( x−3)( x−6)
    L0 ( x )= = =
    ( x 0 −x 1 )( x 0 −x 2 )( x 0 −x 3 ) (0−1 )(0−3 )(0−6) −18
    ( x−x 0 )(x −x 2 )( x−x 3 ) ( x−0 )(x −3)( x−6 ) x( x−3 )( x−6 )
    L1 ( x )= = =
    ( x 1−x 0 )( x 1−x 2 )( x 1−x 3 ) (1−0 )(1−3 )(1−6 ) 10
    ( x−x 0 )( x−x 1 )( x−x 3 ) ( x−0 )( x−3 )( x−6) x (x −1)( x−6 )
    L2 ( x )= = =
    ( x 2−x 0 )( x 2 −x1 )( x 2−x 3 ) (3−0)(3−1)(3−6) −18
    ( x−x 0 )( x−x 1 )( x−x 2 ) x (x −1)( x−3 ) x ( x−1)( x−3 )
    L3 ( x )= = =
    ( x 2−x 0 )( x 3 −x1 )( x 3 −x 2 ) 6 (6−1)(6−3 ) 90
    Operando tenemos la ecuación de la recta:

    x 3 x 2 46
    P3 =− − + x −3
    30 30 15
    Reemplazando hallamos el valor aproximado de la función cuando x = 1.8

    (1. 8)3 (1. 8 )2 46


    P3 (1 . 8)=− − + (1 . 8)−3=2 . 2176
    30 30 15

    También podría gustarte