Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 42 vistas

    Problemas Con Condiciones de Frontera

    Cargado por

    luishernandez142857
    Este documento describe diferentes problemas relacionados con condiciones de frontera en ecuaciones diferenciales parciales resueltas numéricamente mediante diferencias finitas. Explica conceptos como convergencia, consistencia y estabilidad de los métodos numéricos, y cómo tratar diferentes tipos de condiciones de frontera como Dirichlet, Neumann y Robin. También presenta métodos para resolver problemas de autovalores como el método polinomial y el método de potencias. Finalmente, propone algunos ejercicios prácticos para aplicar estos conceptos.

    Copyright:

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

    Problemas Con Condiciones de Frontera

    Cargado por

    luishernandez142857
    42 vistas27 páginas
    Este documento describe diferentes problemas relacionados con condiciones de frontera en ecuaciones diferenciales parciales resueltas numéricamente mediante diferencias finitas. Explica conceptos como convergencia, consistencia y estabilidad de los métodos numéricos, y cómo tratar diferentes tipos de condiciones de frontera como Dirichlet, Neumann y Robin. También presenta métodos para resolver problemas de autovalores como el método polinomial y el método de potencias. Finalmente, propone algunos ejercicios prácticos para aplicar estos conceptos.

    Descripción original:

    Título original

    Problemas con condiciones de frontera

    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?

    Este documento describe diferentes problemas relacionados con condiciones de frontera en ecuaciones diferenciales parciales resueltas numéricamente mediante diferencias finitas. Explica conceptos como convergencia, consistencia y estabilidad de los métodos numéricos, y cómo tratar diferentes tipos de condiciones de frontera como Dirichlet, Neumann y Robin. También presenta métodos para resolver problemas de autovalores como el método polinomial y el método de potencias. Finalmente, propone algunos ejercicios prácticos para aplicar estos conceptos.

    Copyright:

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

    Problemas Con Condiciones de Frontera

    Cargado por

    luishernandez142857
    Este documento describe diferentes problemas relacionados con condiciones de frontera en ecuaciones diferenciales parciales resueltas numéricamente mediante diferencias finitas. Explica conceptos como convergencia, consistencia y estabilidad de los métodos numéricos, y cómo tratar diferentes tipos de condiciones de frontera como Dirichlet, Neumann y Robin. También presenta métodos para resolver problemas de autovalores como el método polinomial y el método de potencias. Finalmente, propone algunos ejercicios prácticos para aplicar estos conceptos.

    Copyright:

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

    Problemas con

    condiciones de
    frontera
    PhD. Alejandro Paredes
    Problemas con condiciones de frontera
    Boundary Value Problems (BVP)

    Consideremos la ecuación del calor de una barra con extremos


    en los puntos 𝑥 = 𝑎 y 𝑥 = 𝑏.

    • 𝑢(𝑥, 𝑡): Temperatura de la barra en un punto e instante


    determinado.
    • 𝜅(𝑥): Coeficiente de conducción dependiente de la
    posición.
    • 𝜓 𝑥, 𝑡 : Fuente de calor (𝜓 > 0) o sumidero (𝜓 < 0).
    Problemas con condiciones de frontera

    • Consideremmos el estado estacionario, 𝜅 𝑥 = 𝜅,


    temperatura fija en los bordes (condiciones de
    frontera de Dirichlet) y extremos 𝑎 = 0, 𝑏 = 1. 𝑢 0 = 𝑢𝑎 , 𝑢 1 = 𝑢𝑏
    • Consideremos la discretización (malla)

    𝑢𝑎 𝑢𝑏

    𝑥0 = 0 𝑥𝑖−1 𝑥𝑖 𝑥𝑖+1 𝑥𝑛 = 1
    𝑒𝑟𝑟𝑜𝑟: Error local o de truncamiento.
    Problemas con condiciones de frontera
    Definimos 𝑈𝑖 como la solución en diferencias finitas. 𝑈𝑖 es una aproximación de
    𝑢(𝑥) en el punto 𝑥𝑖 . Si 𝑈𝑖 existe, entonces verifica el sistema lineal

    𝐴 𝑈 𝐹

    𝐴 : Es la representación matricial del operador 𝑑2 Τ𝑑𝑥 2 .


    Convergencia, consistencia y estabilidad
    Error global

    • Si 𝑼 = 𝑈1 , … , 𝑈𝑛−1 𝑇 es la solución en DF y 𝒖 = 𝑢 𝑥1 , … , 𝑢 𝑥𝑛−1 𝑇 es la


    solución analítica evaluada en los puntos de la malla con paso ℎ, entonces
    definimos el error global
    𝑬=𝑼−𝒖
    • Definición: el método de diferencias finitas es convergente si se tiene una norma
    sobre la malla y lim 𝑬 = 0.
    ℎ→0
    Norma infinita: p-norma : • Definición: el método DF tiene un
    1 orden de precisión 𝑝, si
    𝑛−1 𝑞 𝑞

    𝑬 = max {𝐸𝑖 } 𝑬 = ℎ ෍ 𝐸𝑖 𝑬 ≤ 𝐶 ℎ𝑝 , 𝑝>0


    ∞ 𝑝
    1≤𝑖≤𝑛−1
    𝑖=1 dodne 𝐶 no depende de ℎ.
    Convergencia, consistencia y estabilidad
    Error local

    𝑇
    • En general la solución analítica 𝒖 = 𝑢 𝑥1 , … , 𝑢 𝑥𝑛−1 no satisface exactamente
    la ecuación A𝐮 ≠ 𝑭, mas bien

    𝝉 = 𝜏1 , … , 𝜏𝑛−1 𝑇
    𝝉 = 𝐴𝒖 − 𝑭

    • El vector 𝝉 recibe el nombre de error local. Si definimos los operadores :


    𝑑2𝑢 𝑥 𝑢 𝑥 − ℎ − 2𝑢 𝑥 + 𝑢(𝑥 + ℎ) 𝜏𝑗 = 𝑃ℎ 𝑢 𝑥𝑗 − 𝑃𝑢(𝑥𝑗 )
    𝑃𝑢 𝑥 = 2
    , 𝑃ℎ 𝑢(𝑥) =
    𝑑𝑥 ℎ2
    • Definición: el método de diferencias finitas es consistente si: lim 𝜏𝑗 = 0 ∀ 𝑗 = 1, … , 𝑛 − 1
    ℎ→0
    • Consistencia significa que la ecuación discretizada tiende a la ecuación diferencial en
    cada punto cuando ℎ → 0. En nuestro caso 𝝉 = 𝑂 ℎ2 .
    Convergencia, consistencia y estabilidad
    En general tenemos: Para una malla determinada con paso ℎ

    𝑬=𝑼−𝒖
    𝐴ℎ 𝑬ℎ = −𝝉ℎ , 𝑑𝑖𝑚 𝐴 = 𝑛 − 1 = 1Τℎ − 1
    𝐴𝑼 = 𝑭 𝐴𝑬 = −𝝉
    Donde 𝐴ℎ es la matriz asociada al operador
    𝝉 = 𝐴𝒖 − 𝑭 diferencial 𝑑2 Τ𝑑𝑥 2 y su dimensión depende de ℎ.
    −1
    Sea 𝐴ℎ la inversa de 𝐴ℎ , Definición: un método DF genera una secuencia de
    entonces ecuaciones matriciales de la forma 𝐴ℎ 𝑼ℎ = 𝑭ℎ donde ℎ es el
    −1
    paso de la malla. El método DF es estable si 𝐴ℎ existe
    para todo ℎ suficientemente pequeño (ℎ < ℎ0 ) y si existe una
    constante C, independiente de ℎ, tal que
    −1
    𝐴ℎ ≤ 𝐶 ∀ ℎ < ℎ0 .
    Convergencia, consistencia y estabilidad

    Teorema: Un método DF consistente y estable es convergente.

    Consistencia + Estabilidad ⇒ Convergencia

    cuando

    Error local por truncamiento 𝑂(ℎ2 ) + Estabilidad ⇒ Error global 𝑂(ℎ2 )


    Condiciones de frontera y puntos
    fantasmas
    Condiciones de Dirichlet Condiciones de Robin

    Condiciones de Neumann Condiciones de Cauchy

    𝑢 𝑎 = 𝑢𝑎 , 𝑢 𝑏 = 𝑢𝑏
    𝑢′ 𝑎 = 𝛼, 𝑢′ 𝑏 = 𝛽
    Condiciones de frontera y puntos
    fantasmas
    Condiciones de Neumann Ecuación para el punto 𝑥0

    (Orden 1)

    Para los puntos interiores

    𝑢(𝑎)
    ' 𝑢(𝑏)
    𝑥0 = 𝑎 𝑥𝑖−1 𝑥𝑖 𝑥𝑖+1 𝑥𝑛 = 𝑏
    Condiciones de frontera y puntos
    fantasmas
    Para mantener la precisión al segundo orden necesitamos extender la solución al
    intervalo [𝑎 − ℎ, 𝑎]

    Punto fantasma :

    Condición de frontera al segundo orden:


    Condiciones de frontera y puntos
    fantasmas

    Pendiente ≈ 1

    Sol. analítica

    Pendiente ≈ 2
    Problemas autoadjuntos con condiciones
    de frontera
    Ecuación conducción del calor con conductivad variable (problema de Sturm-Liouville)

    𝑑 𝑑𝑢
    𝐿 : operador auto adjunto −𝐿 𝑢 𝑥 =𝑓 𝑥 , 𝐿 𝑢 = − (𝜅(𝑥) )
    𝑑𝑥 𝑑𝑥

    Dado un producto interno se cumple que: 𝐿𝑢, 𝑣 = (𝑢, 𝐿𝑣)

    Si además escogemos una base, la matriz asociada a 𝐿 (𝐿𝑚𝑛 ) es una matriz simétrica.
    Problemas autoadjuntos con condiciones
    de frontera

    Discretizamos con diferencias centradas de orden 2

    Matriz no simétrica !!!!


    Problemas autoadjuntos con condiciones
    de frontera

    Matriz simétrica !!!!


    Ejercicio 1

    Condiciones de Dirichlet

    • Encontrar la solución analítica. • Graficar el error global por refinamiento.


    • Resolver numéricamente considerando:

    𝐿 = 10𝑚 , ℎ’ = 0.01𝑚−2 ∆𝑥 = 1𝑚 𝑇 0 = 40, 𝑇 10 = 200°𝐶, 𝑇𝑎 = 20°𝐶


    Ejercicio 2

    • Encontrar la solución analítica.


    • Resolver numéricamente considerando:

    𝐿 = 10𝑚 , ℎ’ = 0.01𝑚−2 ∆𝑥 = 1𝑚
    Condiciones de Neumann y Dirichlet 𝑇∞ = 40°𝐶, 𝑇 10 = 200°𝐶, 𝑇′𝑎 = 10°𝐶/𝑚

    y las aproximaciones para la condición


    de Neumann.
    • Graficar el error global por refinamiento.
    Problema de autovalores

    • La variable 𝜔 es un valor propio necesario


    para hallar 𝐴1 y 𝐴2.
    • Se determina 𝜔 y luego se resuelve un
    sistema lineal.
    Problema de autovalores

    Para cuatro puntos dentro del intervalo


    Método polinomial

    A partir del determinante se obtine la ecuación característica que se resuelve


    Numericamente con e método de la bisección, newtonrhapson, etc.

    Los valores de 𝑦1, 𝑦2, 𝑦3, 𝑦4 se obtiene resolviendo el


    sistema lineal
    Método de potencias

    • Se utiliza cuando la matriz A (dim(A)=n) tiene n vectore propios LI {V1, ...Vn}.


    • Los autovalores deben verifivar

    • El método calcula el autovalor dominante λ1.


    • Un vector caulquiera X0 cumple
    Método de potencias

    Multiplicando sucesivamente por la matriz A, tenemos

    Y dividiendo entre λ1m, tenemos Para m grande, tenemos


    Método de potencias

    Para cualquier vector Y que no sea perpendicular a V1


    Método de potencias

    Ejemplo: Estimar el autovalor dominante de

    Comenzamos con cualquiera y calculamos las potencias:


    Método de potencias

    Escogemos y m=5

    Se puede escoger

    Respuesta real:
    Método de potencias
    Consideremos con

    Autovalor

    Autovector
    Ejercicio 4

    Usando el método de potencias encontrar en cada caso el


    valor propio más elevado y su correspondiente vector propio

    También podría gustarte