Universidad Nacional Autónoma de

Honduras
UNAH -VS
logo.jpg

Departamento de Fı́sica
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 función V (x, y) dada por:
AVxx + Buxy + CVyy = f (x, y, V, Vx , Vy )
Elı́ptica: B 2 − 4AC < 0.
Ejemplo : Ecuación de Poisson uxx + uyy = f (x, y).
Parabólica : B 2 − 4AC = 0.
Ejemplo : Ecuación del calor ut = uxx + f (x, t).
Hiperbólica : B 2 − 4AC > 0.
Ejmplo : Ecuación de ondas utt − uxx = f (x, y).

Mallado
Denotamos por Ω la región 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 región Ω, llamaremos:

Vi,j = V (xi , yj ) será el valor exacto de la solución de la E.D.P.

Ui,j al valor numérico proporcionado por un método de resolución aproximada de la E.D.P.

Supongamos que Ω es como se muestra y que hacemos un mallado rectangular

malladoEDP1.jpg

Notemos que algunos puntos de nuestro mallado estarán fuera de Ω, en cambio si Ω es un

rectángulo esto no sucede.

malladoEDP2.jpg

El método de diferencias finitas es mucho mejor cuando las regiones son regulares. Trabaja-
remos considerando que Ω = [a, b]×[c, d], es decir un rectángulo como el de la figura anterior.

Podemos definir
b−a
xi = a + (i − 1)k con k=
N −1
d−c
yj = c + (j − 1)h con h =
M −1
k es el espaciado del mallado en el eje x y h es el espaciado en y.

1
Fórmulas de aproximación de derivadas parciales
Podemos encontrar fórmulas de aproximación de funciones de varias variables al considerar
una serie de Taylor para funciones en varias variables o también extender las que ya tene-
mos, recordando que al derivar parcialmente con respecto a una variable las restantes se
toman como constantes.

De nuestra clase anterior recordemos que

f (x + h) − f (x)
f 0 (x) = + O(h),
h
por lo que

∂V (x, y) V (x + h, y) − V (x, y)
= + O(h)
∂x h
∂V (x, y) V (x, y + h) − V (x, y)
= + O(h).
∂y h

Además recordemos que

f (x + h) − 2f (x) + f (x − h)
f 00 (x) = 2
+ O(h2 ),
h
entonces
∂ 2 V (x, y) V (x + h, y) − 2V (x, y) + V (x − h, y)
2
= + O(h2 )
∂x h2
∂ 2 V (x, y) V (x, y + h) − 2V (x, y) + V (x, y − h)
2
= + O(h2 ).
∂y h2

Siguiendo la misma idea podemos extender todas las fórmulas de aproximación con las que
trabajamos en la clase anterior.

2
Problema de Laplace – Dirichlet
Resolver usando diferencias finitas.
∂ 2 V (x, y) ∂ 2 V (x, y)
+ = 0, (x, y) ∈ Ω = [a, b] × [c, d]
∂x2 ∂y 2
V (x, y) = g(x, y), (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

∂ 2 Vi,j Ui+1,j − 2Ui,j + Ui−1,j

2
≈
∂x h2
∂ 2 Vi,j Ui,j+1 − 2Ui,j + Ui,j−1
2
≈ .
∂y h2

∂ 2 Vi,j ∂ 2 Vi,j
+ =0
∂x2 ∂y 2

Ui+1,j − 2Ui,j + Ui−1,j Ui,j+1 − 2Ui,j + Ui,j−1

+ =0
h2 h2
Agrupando obtenemos

−4Ui,j + Ui+1,j + Ui−1,j + Ui,j+1 + Ui,j−1 = 0, i, j = 1, 2, 3, . . . , M .

La molécula computacional del método es

molecula.jpg

Por la condición de contorno sabemos el valor de V en ciertos nodos, estos nodos están
marcados en rojo en la figura siguiente

malladoEDP3.jpg

De lo anterior

−4Ui,j + Ui+1,j + Ui−1,j + Ui,j+1 + Ui,j−1 = 0, i, j = 2, 3, . . . , M − 1.

3
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, además

 
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 g12 + g21

 −1 4 −1 0 −1 0 0 0 0 

 U32 


 g31 


 0 0 4 0 0 −1 0 0 0 
  U42
 


 g41 + g52 


 −1 0 0 4 −1 0 −1 0 0 
  U23
 


 g13 


 0 −1 0 −1 4 −1 0 −1 0 

 U33
 = 
  0 


 0 0 −1 0 −1 4 0 0 −1  

 U43



 g53 


 0 0 0 −1 0 0 4 −1 0 
  U24
 


 g14 + g25 

 0 0 0 0 −1 0 −1 4 −1   U34   g35 
0 0 0 0 0 −1 0 −1 4 U44 g45 + g54

4
1. Ejercicios propuestos
Resolver en forma clara y ordenada dejando evidencia de sus procedimientos.

1. Obtener la expresión matricial y la molécula computacional para el método de dife-

rencias finitas en los siguientes problemas de contorno. Considerar N = M = 4

∂ 2 V (x, y) ∂V (x, y)
2
= , (x, y) ∈ Ω = [a, b] × [c, d]
∂x ∂y
v(x, y) = g(x, y), (x, y) ∈ ∂Ω.

Usando cada una de las fórmulas de aproximación,

∂ 2 Vi,j Vi+1,j − 2Vi,j + Vi−1,j ∂Vi,j Vi,j+1 − Vi,j

a) 2
= 2
+ O(h2 ), = + O(h).
∂x h ∂y h

∂ 2 Vi,j Vi+1,j − 2Vi,j + Vi−1,j ∂Vi,j Vi,j+1 − Vi,j−1

b) = + O(h2 ), = + O(h2 ).
∂x2 h2 ∂y 2h

∂ 2 Vi,j Vi+2,j − 2Vi+1,j + Vi,j ∂Vi,j Vi,j − Vi,j−1

c) 2
= + O(h), = + O(h).
∂x h2 ∂y h

2. Aplicando diferencias finitas para cada caso y con la ayuda de una calculadora pro-
gramable encontrar la solución en cada nodo; usar las aproximaciones:

∂ 2 Vi,j Ui+1,j − 2Ui,j + Ui−1,j

≈
∂x2 k2
∂ 2 Vi,j Ui,j+1 − 2Ui,j + Ui,j−1
2
≈ .
∂y h2

a) ∇2 V = 0

proble1.jpg

b) ∇2 V = x(y − 1), Ω = [0, 1] × [0, 1].

proble2.jpg