Metodo de Diferencias Finitas para EDP de Segundo Orden.

Universidad Nacional Aut
onoma de
Honduras
UNAH -VS
Departamento de Fsica
Laboratorio de electricidad y magnetismo 2
LF-415
Apuntes Clase #2
Diferencias finitas para E.D.P.

En lo sucesivo consideramos las notaciones para E.D.P. siguientes: Supongamos una E.D.P.
de orden 2 para la funcion V (x, y) dada por:
AVxx + Buxy + CVyy = f (x, y, V, Vx , Vy )
Elptica: B 2 4AC < 0.
Ejemplo : Ecuacion de Poisson uxx + uyy = f (x, y).
Parabolica : B 2 4AC = 0.
Ejemplo : Ecuacion del calor ut = uxx + f (x, t).
Hiperbolica : B 2 4AC > 0.
Ejmplo : Ecuacion de ondas utt uxx = f (x, y).
Mallado
Denotamos por la region del plano sobre la que se define V (x, y), R es la frontera de .
Considerada una red o malla de puntos
{(xi , yj ) : i = 0, 1, . . . , N + 1, j = 0, 1, . . . , M + 1}.
Para la region , llamaremos:
Vi,j = V (xi , yj ) sera el valor exacto de la solucion de la E.D.P.
Ui,j al valor numerico proporcionado por un metodo de resolucion aproximada de la E.D.P.
Supongamos que es como se muestra y que hacemos un mallado rectangular
Notemos que algunos puntos de nuestro mallado estaran fuera de , en cambio si es un

rectangulo esto no sucede.
El metodo de diferencias finitas es mucho mejor cuando las regiones son regulares. Trabajaremos considerando que = [a, b][c, d], es decir un rectangulo como el de la figura anterior.
Podemos definir
xi = a + (i 1)k
con
yj = c + (j 1)h
con
ba
N 1
dc
h=
M 1
k=
k es el espaciado del mallado en el eje x y h es el espaciado en y.
F
ormulas de aproximaci
on de derivadas parciales
Podemos encontrar formulas de aproximacion de funciones de varias variables al considerar
una serie de Taylor para funciones en varias variables o tambien extender las que ya tenemos, recordando que al derivar parcialmente con respecto a una variable las restantes se
toman como constantes.
De nuestra clase anterior recordemos que
f 0 (x) =
f (x + h) f (x)
+ O(h),
h
por lo que
V (x + h, y) V (x, y)
V (x, y)
=
+ O(h)
x
h
V (x, y)
V (x, y + h) V (x, y)
=
+ O(h).
y
h
Ademas recordemos que
f 00 (x) =
f (x + h) 2f (x) + f (x h)
+ O(h2 ),
h2
entonces
2 V (x, y)
V (x + h, y) 2V (x, y) + V (x h, y)
=
+ O(h2 )
x2
h2
2 V (x, y)
V (x, y + h) 2V (x, y) + V (x, y h)
=
+ O(h2 ).
2
2
y
h
Siguiendo la misma idea podemos extender todas las formulas de aproximacion con las que
trabajamos en la clase anterior.
Problema de Laplace Dirichlet

Resolver usando diferencias finitas.
2 V (x, y) 2 V (x, y)
+
= 0,
x2
y 2
V (x, y) = g(x, y),
(x, y) = [a, b] [c, d]

(x, y) .
Mallado bidimensional
(xi , yj ) = (a + (i 1)k, b + (j 1)h).
Tomando h = k.
(xi , yj ) = (a + (i 1)h, b + (j 1)h).
Aproximaciones de segundo orden
Ui+1,j 2Ui,j + Ui1,j
2 Vi,j
2
x
h2
Ui,j+1 2Ui,j + Ui,j1
2 Vi,j
.
y 2
h2
2 Vi,j 2 Vi,j
+
=0
x2
y 2
Ui+1,j 2Ui,j + Ui1,j Ui,j+1 2Ui,j + Ui,j1
+
=0
h2
h2
Agrupando obtenemos
4Ui,j + Ui+1,j + Ui1,j + Ui,j+1 + Ui,j1 = 0,
i, j = 1, 2, 3, . . . , M .
La molecula computacional del metodo es
Por la condicion de contorno sabemos el valor de V en ciertos nodos, estos nodos estan
marcados en rojo en la figura siguiente
De lo anterior
4Ui,j + Ui+1,j + Ui1,j + Ui,j+1 + Ui,j1 = 0,
3
i, j = 2, 3, . . . , M 1.
Forma matricial del problema

AU = B
Donde
U = {U22 , U32 , . . . , UM 1,2 , U23 , U33 , . . . , UM 1,3 . . . , U2,M 1 , U3,M 1 , . . . , UM 1,M 1 }T
D I
I D I
I D I
A=
,
.
.
.
.
.
.
.
.
.
I D I
I D
donde I denota a matriz identidad de dimension M 2, ademas
4 1
1
4 1
1
4 1
D=
.
.
.
.
.
.
.
.
.
1
4 1
1
4
BT =
{g12 + g21 , g31 , . . . , gM 2,1 , gM 1,1 + gM,2 , g13 , . . . , gM,3 , . . . , g1,M 1, + g2,M , g3,M , . . . , gM 1,M, + gM,M 1 },
donde g(xi , yi ) = gi,j .
Caso M = 5
AU = B
4 1
0 1
0
0
0
0
0
U22
U32
1
4 1
0 1
0
0
0
0
0
0
4
0
0 1
0
0
0
U42
1
0
0
4 1
0 1
0
0
U23
0 1
0 1
4 1
0 1
0
U33
0
0 1
0 1
4
0
0 1
U43
0
0
0 1
0
0
4 1
0
U24
0
0
0
0 1
0 1
4 1 U34
0
0
0
0
0 1
0 1
4
U44
g12 + g21
g31
g41 + g52
g13
0
g53
g14 + g25
g35
g45 + g54
1.
Ejercicios propuestos
Resolver en forma clara y ordenada dejando evidencia de sus procedimientos.

1. Obtener la expresion matricial y la molecula computacional para el metodo de diferencias finitas en los siguientes problemas de contorno. Considerar N = M = 4
2 V (x, y)
V (x, y)
, (x, y) = [a, b] [c, d]
=
2
x
y
v(x, y) = g(x, y), (x, y) .
Usando cada una de las formulas de aproximacion,
a)
2 Vi,j
Vi+1,j 2Vi,j + Vi1,j
=
+ O(h2 ),
2
x
h2
Vi,j+1 Vi,j
Vi,j
=
+ O(h).
y
h
b)
2 Vi,j
Vi+1,j 2Vi,j + Vi1,j
=
+ O(h2 ),
2
x
h2
Vi,j
Vi,j+1 Vi,j1
=
+ O(h2 ).
y
2h
c)
2 Vi,j
Vi+2,j 2Vi+1,j + Vi,j
=
+ O(h),
2
x
h2
Vi,j
Vi,j Vi,j1
=
+ O(h).
y
h
2. Aplicando diferencias finitas para cada caso y con la ayuda de una calculadora programable encontrar la solucion en cada nodo; usar las aproximaciones:
Ui+1,j 2Ui,j + Ui1,j
2 Vi,j
2
x
k2
Ui,j+1 2Ui,j + Ui,j1
2 Vi,j
.
2
y
h2
a) 2 V = 0
b) 2 V = x(y 1),
= [0, 1] [0, 1].

Metodo de Diferencias Finitas para EDP de Segundo Orden.

Cargado por

Copyright:

Formatos disponibles

Metodo de Diferencias Finitas para EDP de Segundo Orden.

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Metodo de Diferencias Finitas para EDP de Segundo Orden.

Cargado por

Copyright:

Formatos disponibles

Universidad Nacional Aut

Diferencias finitas para E.D.P.

Notemos que algunos puntos de nuestro mallado estaran fuera de , en cambio si es un

k es el espaciado del mallado en el eje x y h es el espaciado en y.

Problema de Laplace Dirichlet

(x, y) = [a, b] [c, d]

La molecula computacional del metodo es

Forma matricial del problema

Resolver en forma clara y ordenada dejando evidencia de sus procedimientos.

= [0, 1] [0, 1].

También podría gustarte