Problemas Con Condiciones de Frontera
Problemas Con Condiciones de Frontera
Problemas Con Condiciones de Frontera
condiciones de
frontera
PhD. Alejandro Paredes
Problemas con condiciones de frontera
Boundary Value Problems (BVP)
𝑢𝑎 𝑢𝑏
𝑥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
𝐴 𝑈 𝐹
𝑇
• En general la solución analítica 𝒖 = 𝑢 𝑥1 , … , 𝑢 𝑥𝑛−1 no satisface exactamente
la ecuación A𝐮 ≠ 𝑭, mas bien
𝝉 = 𝜏1 , … , 𝜏𝑛−1 𝑇
𝝉 = 𝐴𝒖 − 𝑭
𝑬=𝑼−𝒖
𝐴ℎ 𝑬ℎ = −𝝉ℎ , 𝑑𝑖𝑚 𝐴 = 𝑛 − 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
cuando
𝑢 𝑎 = 𝑢𝑎 , 𝑢 𝑏 = 𝑢𝑏
𝑢′ 𝑎 = 𝛼, 𝑢′ 𝑏 = 𝛽
Condiciones de frontera y puntos
fantasmas
Condiciones de Neumann Ecuación para el punto 𝑥0
(Orden 1)
𝑢(𝑎)
' 𝑢(𝑏)
𝑥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 :
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 −𝐿 𝑢 𝑥 =𝑓 𝑥 , 𝐿 𝑢 = − (𝜅(𝑥) )
𝑑𝑥 𝑑𝑥
Si además escogemos una base, la matriz asociada a 𝐿 (𝐿𝑚𝑛 ) es una matriz simétrica.
Problemas autoadjuntos con condiciones
de frontera
Condiciones de Dirichlet
𝐿 = 10𝑚 , ℎ’ = 0.01𝑚−2 ∆𝑥 = 1𝑚
Condiciones de Neumann y Dirichlet 𝑇∞ = 40°𝐶, 𝑇 10 = 200°𝐶, 𝑇′𝑎 = 10°𝐶/𝑚
Escogemos y m=5
Se puede escoger
Respuesta real:
Método de potencias
Consideremos con
Autovalor
Autovector
Ejercicio 4