Método Numérico
Método Numérico
Método Numérico
1
1. Introducción
El método de diferencias finitas es un clásica aproximación para encontrar
la solución numérica de las ecuaciones diferenciales parciales. El método
consiste básicamente en que las derivadas son reemplazadas por
aproximaciones en diferencias finitas, convirtiendo el problema de
ecuaciones diferenciales en un problema algebraico Ax = b que puede
resolverse por métodos matriciales.
1.1 Aproximaciones polinómicas
O en su forma cuadrática:
𝑓′′(𝑥0 )
𝑓 𝑥 ≈ 𝑓 𝑥0 + ℎ 𝑓’ 𝑥0 + ℎ2
2!
Mientras menor sea h más aproximada será la solución.
1.2 Error de truncamiento
𝑓 𝑥 + ℎ − 𝑓(𝑥)
𝑓′ 𝑥 = − 𝑂𝑝 (ℎ)
ℎ
𝑓 𝑥 −𝑓(𝑥−ℎ)
𝑓 ′ (𝑥) = − 𝑂𝑟 (h)
ℎ
ℎ2
Siendo 𝑂𝑟 ℎvez
una = mas 𝑓′′(∈) el error de truncamiento.
2!
2.2 Aproximación de la primera derivada
ℎ2 ′′ ℎ3
𝑓 𝑥 − ℎ = 𝑓 𝑥 − ℎ𝑓 ′ 𝑥
+ 𝑓 ∈ − 𝑓′′′(∈𝑝 )
2! 3!
2.2 Aproximación de la primera derivada
𝑓 𝑥+ℎ −𝑓 𝑥−ℎ
𝑓 ′ (𝑥) = − 𝑂𝑐 (ℎ2 )
2ℎ
2.2 Aproximación de la primera derivada
Donde:
Derivada de Orden Uno en Una Dimensión
Robin.
1
Para 𝜽 = se obtiene un esquema de diferencias finitas centradas en el tiempo,
2
conocido como esquema Crank-Nicolson, este esquema también es incondicionalmente
estable y es el más usado por tal razón.
Entonces en el esquema de Crank-Nicolson se toman una diferencia
progresiva para el tiempo y se promedian las diferencias progresivas
y regresivas en el tiempo para las derivadas espaciales.
𝒋
𝒖𝒋+𝟏 −𝒖𝒊
𝒖𝒕 ≈
𝒌
Para el tiempo, donde k=∆𝒕 𝒀
𝑗 𝑗
además de 𝑔𝑢𝑖 𝑦 ℎ𝑖 para el espacio.
𝒋
𝒖𝒋+𝟏 −𝟐𝒖𝒊 +𝒖𝒋−𝟏
𝒖𝒕𝒕 ≈
𝒌𝟐
𝒋
𝒖𝒋+𝟏 −𝟐𝒖𝒊 +𝒖𝒋−𝟏
𝒖𝒕𝒕 ≈
𝒌𝟐
𝒖𝒕 + 𝒂𝒖𝒙 = 𝟎
Error global: Se define como U-u, sea U las soluciones del método y u las
soluciones analíticas.
Norma infinito: 𝐸 ∞ = 𝑚á𝑥𝑖 𝑒𝑖
Norma 1: 𝐸 1 = σ𝑖 ℎ𝑖 𝑒𝑖
2 1Τ2
Norma 2: 𝐸 2 = ( 𝑖 ℎ𝑖 𝑒𝑖 )
σ
La norma infinito es la más usada para calcular los errores en la
discretización.
El error local de truncamiento sirve para medir que tanto se acerca la
discretización de diferencias finitas a la ecuación diferencial y se define
𝑢 𝑥−ℎ −2𝑢 𝑥 +𝑢(𝑥+ℎ)
como: 𝑇 𝑥 = 𝑃𝑢 − 𝑃ℎ 𝑢 = 𝑓 𝑥 − 2
ℎ
Convergencia, consistencia y estabilidad
Para determinar la condición de estabilidad de los métodos con paso de tiempo, se utiliza el
análisis de estabilidad de von Neumann, que puede describirse en los siguientes pasos:
Esquema discreto ⇒ transformada discreta de Fourier ⇒ crecimiento del factor g (ξ) ⇒
estabilidad si |g (ξ) | ≤ 1
10. Estabilidad de métodos con paso de
tiempo.
10.1 Análisis de von Neumann Simplificado para métodos con Paso
de Tiempo 𝑈 𝑘+1 = 𝑓(𝑈 𝑘 , 𝑈 𝑘+1 )
Asúmase que se tiene un método de paso temporal
Se tiene el siguiente teorema con el cual se puede conocer la
estabilidad del esquema de diferencias finitas
𝑔, 𝜃, ∆𝑡, ℎ ≤ 1 + 𝐾∆𝑡
10. Estabilidad de métodos con paso de
tiempo.
10.1 Análisis de von Neumann Simplificado para métodos con Paso de Tiempo
para todo Θ, 0 < h ≤ h0 y 0 < ∆t ≤ ∆t0. Si g (Θ, ∆t, h) es independiente de h y ∆t, la
condición para la estabilidad es
𝑔(𝜃) ≤ 1
3. Se resuelve 𝑔 ξ y se determina si |𝑔 ξ | ≤ 1
Conclusión