Captulo 4

Ecuaci on de Laplace
4.1. Condiciones de frontera
4.2. Metodo de la funci on de Green
Si los problemas electrosaticos solo envolviesen cargas (o densidades de
carga) sin condiciones de frontera entonces =

k(

k)
|xx

|
dv

sera sucien-
temente.
Pero muchos problemas electrostaticos envuelven regiones nitas del es-
pacio, con o sin carga adentro, y con condiciones de frontera.
Esas condiciones de frontera pueden ser simuladas por una apropiada dis-
tribucion de cargas fuera de la region de interes, pero solo es util en casos
muy simples( Ej: metodo de imagenes).
Para lidiar con este problema. Suponga dos funciones escalares: tal
que

A = . Donde

A es una funcion bien comportada, denida dentro
de un volumen V y una supercie S. Ademas

A sin discontinuidades e
integrable.
Usando el teorema de divergencia:

Adv =

Ads

A = () =
2
+

A n = n =

n
Primera identidad de Green . George Green (1824):

V
(
2
+ )dV =

n
dS (4.1)
Intercambiando

V
(
2
+ )dV =

n
dS (4.2)
31
32 CAP

ITULO 4. ECUACI

ON DE LAPLACE
Ahora restando (3.1)-(3.2), se obtiene la segunda identidad de Green:

V
(
2

2
)dV =

n
)dS (4.3)
Usando la ecuacion de Poisson:

2
=

0
Ademas se escoge =
1
|rr

|
=
1
R
1
y del hecho que
2
(
1
R
) = 4(r r

). En
(3.1):

V
(
2
1
R
)dV

(4(r r

) +
(r

0
R
)
= 4(r

) +

(r

dV

0
|r r

(
1
R
)
1
R

dS

Si r esta dentro de V:
(r) =
1
4
0

V
(r

)dV

R
+
1
4

1
R

n
(
1
R
)

dS

1
4
H
S
h
(r)
|rr

|

(r)(rr

)
|rr

|
3
i
ndS

2
Comentarios:
Que son los otros terminos que aparecen?. Son terminos introducidos por
la libertad que hay en la ecuacion de Poisson. Se deben jar usando con-
diciones de frontera. Aqu no asumimos que sea la carga del universo.
de todas maneras esto tiene demasiada libertad. Notar qu estamos cons-
truyendo a partir de , de su valor (de ) en la frontera, y de sus
derivadas en la frontera. No debemos jar estas dos cantidades (fronte-
ra) y (frontera), arbitrariamente y esperar una solucion consistente.
que podemos hacer es jar una de ellas y pedir
2
= 0
Lo mejor que podems hacer es introducir una funcion G(r,r

)que tenga la
propiedad:

2
G(r,r

) = 4(r r

)
Ej:
G(r,r

) =
1
|r

r

|
+F(r r

)
1
No es bien comportada pero podemos usar lmites o excluir
vectr

deV
2
Condiciones de frontera: Dirichlet = 0sobresuperficieS; Newman

n
= 0 sobre
supercie S
4.2. M

ETODO DE LA FUNCI

ON DE GREEN 33
Tal que F
3
:
2
F(r,r

) = 0 La funcion F se encarga de manejar las condiciones

de frontera: usemos G(r,r

) por en: ( )
_
V
(
2

2
)dV =
_
S
(

n
)dS
entonces:
_
V
(G(r,r

)
2

2
G(r,r

))dV =
_
S
(G(r,r

) G(r,r

)) ndS
_
_
G(r,r

)
2
+ 4(r

)(

r r

)
_
dV

=
(r) =
_
V
G(r,r

)
(r

)
4
0
dV

+
1
4
_ _
G(r,r

)
(r

)
n

(r

)
G(r,r

)
n

_
dS

Ahora podemos tratar de especicar las condiciones de frontera, apro-

vechando que aun no hemos decidido que usar como funcion de Green.
Podemos por ejemplo, pedir que:
G
Dirichlet
(r,r

)|
s
= 0
En ese caso:
(r) =
_
V
dV

G
D
(r,r

)(r

)
1
4
_
S
dS

(r

)
G(r,r

)
n

Ahora podemos intentar hacer lo contrario, es decir tomar

G(r,r

)
n

|
S
= 0
en la supercie. Pero
2
G(r,r

) = 4(r r

). Entonces:
_

2
G(r,r

)dV =
_
G(r,r

)dV = 4
_
S
dSG(r,r

) n =
_
S
dS

G(r,r

)
n

= 4 = 0
Por lo tanto, no podemos tomar
G(r,r

)
n

|
s
= 0, lo mas que podemos hacer
es tomar:
G
V on Newman
(r,r

)
n
=
4
S
Se tendra entonces:
(r) =
_
dV

G
N
(r,r

)
(r

)
4
0
+
1
4
_
S
dS

G
N
(r,r

)
(r

)
n

+
1
S
_
S
dS

(r

)
=
_
dV

G
N
(r,r

)
(r

)
4
0
+
1
4
_
S
dS

G
N
(r,r

)
(r

)
n

+ < (r) >

S
. .
V alorpromediodesobrelasuperficieS
Nota: La solucion de ecuaciones homogeneas (o para condiciones de frontera
homogenea) pueden multiplicarse por un factor y aun ser solucion. Las solucio-
nes de ecuaciones inhomogeneas (o para condiciones de frontera inhomogeneas)
3
F =
1
|r

|
es s olo una clase de funcion que dependen de r y r

y que satisfacen

2
F(r,r

)(r

)
34 CAP

ITULO 4. ECUACI

ON DE LAPLACE
no pueden ser modicadas. La funcion de Green es una solucion para el caso
que es homogenea en todas partes excepto en un punto . Cuando este punto
esta en la supercie (frontera), la funcionde Green puede usarse para satisfacer
las condiciones de frontera inhomogeneas. Cuando ese punto esta afuera en el
espacio, puede ser usada para satisfacer ecuaciones inhomogeneas.
Ecuacion por un campo en presencia de fuentes es una ecuacion diferencial
inhomogenea (ej: ec. de Poisson
2
=

0
).
Ecuacion sin fuentes presentes es una ecuacion difertencial homogenea (ej:
ec LaPlace
2
= 0)
Analogamente:
Condiciones de frontera requiriendo campo zero en la supercie : Condi-
ciones de frontera homogeneas(condicion de Dirichlet).
Condiciones de frontera requiriendo gradiente normal zero en la supercie:
condiciones de frontera homogenea (condicion de Von Newman).
Condiciones de frontera requiriendo a + b

n
= 0: condicion de frontera
mixta.
Principio de Reciprocidad:
El intercambio de fuente (r

) y observado (r), no cambia G(r,r

) = G(r

,r).
Esta funcion va a satisfacer ciertas condiciones denidas en S y S y tendra una
discontinuidad en r =r

Para recapitular
Una alternativa a la tecnica de autofunciones (ej: armonicos esfericos, problema
de Legendre; polinomios de Bessel, que veremos mas adelante) para resolver
problemas con condiciones de frontera. Resolver la ecuacion inhomogenea para
las funciones homogeneas para una fuente puntual en alg un punto de r

dentro
de la supercie. La funcion de Green resultante tiene una descontinuidad en
r =r

4.3. Metodos de las imagenes

Pasemos a otro problema:
V (r) =
1
k

q
|r d

k|

q
|r + d

k|

4.3. M

ETODOS DE LAS IM

AGENES 35
= k

x
2
+ y
2
+ (z d)
2

x
2
+ y
2
+ (z + d)
2

=
0
V
n
=
0
V
z
|
z=0
|
V
z
= k

q(z d)
( )
3
2
+
q(z + d)
( )
3
2

|
z=0
(x, y) =
qd
2(x
2
+ y
2
+ d
2
)
3
2
Q =

da = q
4.3.1. Teorema de la unicidad(uniqueness)
Dos soluciones de la ecuacion de LaPlace que satisfacen las misma condicio-
nes de frontera dieren cuando mucho en una constante aditiva.
=
kq
|r z

k|
+
kq

|r z

k|
=
kq
|r n z

k|
+
kq

|r n z

k|
36 CAP

ITULO 4. ECUACI

ON DE LAPLACE
Si hacemos :r = R y = 0, entonces obtendremos una relacion entre q, q, R,
z:
(r = R) =
kq
R| n
z
R

k|
+
kq

| n
R
z

k|
= 0
Entonces:
q
R
=
q

z
R
=
R
z

con lo que se obtendra:

q

=
R
z
q
z

=
R
2
z
Regresando a la ecuacion preliminar
=
kq
|r n z

k|

k
R
z
q
|r n
R
2
z

k|
= kq
1
r

1
| n
z
r

k|

R
z
| n
R
2
zr

k|

entonces
=
0

r
|
r=R
=
q
4R
2

R
z

1
R
2
z
2
(1 +
R
2
z
2
2
R
z
cos )
3
2
Solucion a la ecuacion de LaPlace en dos dimensiones (coordenadas
cartesianas)

2
=

2

x
2
+

2

y
2
4.4. M

ETODO DE SEPARACI

ON DE VARIABLES 37
Por separacion de variables: = X(x)Y (y):
=

2
(XY )
x
2
+

2
(XY )
y
2
= 0
= Y
d
2
X
dx
2
+ X
d
2
Y
dy
2
= 0
=
1
X
d
2
X
dx
2

a
2
+
1
Y
d
2
Y
dy
2

a
2
= 0
Con esto las soluciones tendr an la forma de :
X = A
1
e
ax
+ A
2
e
ax
o X = A
3
cosh ax + A
4
senh ax
Y = B
1
e
iay
+ B
2
e
iay
o Y = B
3
cos ay + B
4
sen ay
Entonces la solucion para , sera:

a
(x, y) = (A
1
e
ax
+ A
2
e
ax
)(B
1
e
iay
+ B
2
e
iay
)
o de otra forma:

a
(x, y) = (A
3
cosh ax + A
4
senh ax)(B
3
cos ay + B
4
sen ay)
Analogamente para 3 dimensiones: El exito de este metodo de separacion de
variables esta en dos de sus propiedades: Completitud y ortogonalidad.
Un conjunto de funciones es completo si cualquier otra funcion f(x) puede
ser expresada como f(x) =

n=1
c
n
f
n
(x)

f
n
(x)f
m
(x)dx = 0 con n = m
4.4. Metodo de Separaci on de Variables
4.4.1. Coordenadas esfericas y polinomios de Legendre
simetra azimutal
1
r
2

r
2

+
1
r
2
sen

sen

= 0
38 CAP

ITULO 4. ECUACI

ON DE LAPLACE
Tecnica de separacion de variables
= z(r) P()
1
r
2
P()
d
dr

r
2
dz
dr

+
z(r)
r
2
sin
d
d

sen
dP
d

= 0
1
z
d
dr

r
2
dz
dr

=
1
P sin
d
d

sen
dP
d

= k (cte)

1
sin
d
d

sen
dP
d

+ kP = 0
sabemos dx = d cos = sin d

d
dx

(1 x
2
)
d
dx

+ n(n + 1)

P = 0
d
dr

r
2
dz
dr

= n(n + 1)z
entonces:
z
n
=

r
n
r
(n+1)
z(r) = Ar
n
+
B
r
n+1

n
(r, ) =

A
n
r
n
+
B
n
r
n+1

P
n
(cos )
=

n
Ejemplo:
=?
2
= 0
(r, ) =

A
0
r
0
+
B
0
r

P
0
+

A
1
r +
B
1
r
2

cos +

A
2
r
2
+
B
2
r
3

1
2
(3 cos
2
1)+
4.4. M

ETODO DE SEPARACI

ON DE VARIABLES 39

E(r, )

r
=

E
0
= E
0

k
(r, )|
r
= E
0
z + cte = E
0
r cos + cte
=
_

E d

l
(r, )|
r
= (A
0
) + (A
1
r) cos +
_
A
2
r
2
_
1
2
(3 cos
2
1) +
Entonces:
A
1
= E
0
A
0
= cte
A
n
= 0, n 2
(r, ) =
_
A
0
+
B
0
r
_
+
_
E
0
r +
B
1
r
2
_
cos +
_
B
2
r
3
_
1
2
(3 cos
2
1) +
(a, ) =
a
A
0
=
a
, B
1
= E
0
a
3
yB
n
= 0n 2
(r, ) =
a
E
0
r cos +
E
0
a
3
cos
r
2
, r a
(r, ) =
a
, r a
E
r
=

r
= E
0
(1 + 2
a
3
r
3
) cos
E

r
= E
0
(1
a
3
r
3
) sin
() =
0
E
r
|
r=a
= 3
0
E
0
cos
4.4.2. Funciones Asociadas de Legendre y Arm onicos Esferi-
cos

2
(r, , ) =
1
r

2
r
2
(r) +
1
r
2
sen

_
sen

_
+
1
r
2
sen
2

2
= R(r)P()Q()
1
rR
d
2
dr
2
(rR) +
1
Pr
2
sen
d
d
_
sen
d
d
_
+
1
Qr
2
sen
2

d
2
Q
d
2
= 0
r
2
sen
2

_
1
r
2
R
d
dr
_
r
2
dR
dr
_
+
1
Pr
2
sen
d
d
_
sen
dP
d
__
+
1
Q
d
2
Q
d
2
. .
Q=e
im
= 0
1
R
d
dr
_
r
2
dR
dr
_
. .
+
+
1
P sen
d
d
_
sen
dP
d
_

m
2
sen
2

. .

= 0
De (*)
d
dr
_
r
2
dR
dr
_
= n(n + 1)R
40 CAP

ITULO 4. ECUACI

ON DE LAPLACE
R =
_
r
n
r
(n+1)
De (**)
1
sen
d
d
_
sen
dP
d
_

m
2
sen
2

= 0
Para m
2
= 0 P
n
(cos ) Polinomios de Legendre. PolinomiosAsociadosdeLegendre
P
m
n
(x) = (1)
m
(1 x
2
)
m
2
d
m
dx
m
P
n
(x)
Se sabe que:
P
n
(x) =
1
2
n
n!
_
d
dx
_
n
(x
2
1)
n
P
m
n
(x) =
(1)
m
2
n
n!
(1 x
2
)
m
2
d
n+m
dx
n+m
(x
2
1)
n
P
m
n
= (1)
m
(n m)!
(n +m)!
P
n
(x)
_
1
1
P
n
(x)P
n
(x)dx =
2
2n + 1

nn

_
1
1
P
m
n
(xP
m
n
)dx =
2
2n + 1
(n +m)!
(n m)!

nn

(r, , ) = R
n
(r)Y
nm
(, )

2
=
_
p
2
r
+

L
2
r
2
_
donde:

L =
1
2
r
Y
nm
(, ) =

2n + 1
4
(n m)!
(n +m)!
P
m
n
(cos )e
im
Y
nm
(, ) = (1)
m
Y

nm
(, )
Ortogonalidad:
_
2
0
d
_

0
sen d
. .
R
d
Y

m
(, )Y
nm
(, ) =
nn
mm

Completitud:

l=0
l

m=l
Y

nm
(

)Y
nm
(, ) = (

)(cos ) cos

L
2
Y
nm
(, ) = n(n + 1)Y
nm
(, )

LY
nm
(, ) = mY
nm
(, )
4.4. M

ETODO DE SEPARACI

ON DE VARIABLES 41
Teorema de Adicion:
P
n
(cos ) =
4
2n + 1
n

m=n
Y

nm
(

)Y
nm
(, )
donde cos = cos cos

+ sen sen

cos(

1
|r r

|
= 4

n=0
n

m=n
1
2n + 1
r
n
r
n+1
Y

nm
(

, )Y
nm
(, )

n=0

r
n
r
n+1
P
n
(cos )
Solucion General
1
|r r

|
=
1
r

n=0
_
r

r
_
n
P
n
(cos )(r r

n=0
r
n
r
n+1
P
n
(cos )
(r, , ) =

n=0
n

m=n
_
A
nm
r
n
+ B
nm
r
(n+1)
_
Y
nm
(, )
si conocemos en algun valor de r(r )
_
d(a, , ) = Y

nm
(, ) = A
nm
a
n
+
B
nm
a
n+1
Y
lm
(, )
l = 0 Y
00
=
1

4
=
_
1
4
P
0
(cos )
l = 1
_
_
_
Y
11
=
_
3
4
sen e
i
Y
10
=
_
1
4
cos =
_
3
4
P
1
(cos )
l = 2
_

_
Y
22
=
1
4
_
15
2
sen
2
e
2i
Y
21
=
_
15
8
sen cos e
i
Y
20
=
_
5
4
_
3
2
cos
2

1
2
_
=
_
5
4
P
2
(cos )
l = 3
_

_
Y
33
=
1
4
_
35
4
sen
3
e
3i
Y
32
=
1
4
_
105
2
sen
2
cos e
2i
Y
31
=
1
4
_
21
4
sen (5 cos
2
1)e
i
Y
30
=
_
7
4
_
5
2
cos
3

3
2
cos
_
=
_
7
4
P
3
(cos )
1
|r r

|
=
1
r

n=0
_
r

r
_
n
P
n
(cos )
Ejemplo:
42 CAP

ITULO 4. ECUACI

ON DE LAPLACE
(R, ) = () = V
0
(10 cos
3
6 cos ) = V
0
4P
3
(cos )
=

A
n
r
n
+
Bn
r
n+1

P
n

int
=

n=1
A
n
r
n
P
n
+A
0

e
xt =

n+1
B
n
r
n+1
P
n
+A
0

int
(R, ) =
ext
(R, )

n+1
A
n
R
n
P
n
=

n+1
B
n
R
n+1
P
n
B
n
= A
n
R
2n+1
Ademas:

ref
= () = 0 A
0
= 0

int
(R, ) =

n+1
A
n
R
n
P
n
= V
0
4P
3
(cos )
A
n
= 0n = 3
A
3
R
3
P
3
= V
0
4
A
3
=
4V
0
R
3
B
3
= A
3
R
7
= 4V
0
R
4

int
(r, ) =
4V
0
R
3
r
3
P
3
(cos )

ext
(r, ) =
4V
0
R
4
r
4
P
3
(cos )
4.4. M

ETODO DE SEPARACI

ON DE VARIABLES 43
Ejemplo: Un dipolo electrico pse encuentra en el interior de un cascaron esferico
conducor (radio R) conectado a tierra. Hallar el potencial en todo el espacio y
la densidad de carga en la supercie.

int
(r, ) =

A
n
r
n+1
+ B
n
r
n

P
n
cerca del origen:

int
(r 0, ) = k
p r
r
3
= k
p cos
r
2
=
A
1
= pk, A
n
= 0n = 1
En la supercie, hay conexion a tierra.
ext
= 0

int
(R, ) =
kp
R
2
P
1
(cos ) +

n=0
B
n
R
n
P
n
= 0
B
n
= 0 n = 1
kp
R
2
P
1
+ B
1
RP
1
= 0 B
1
=
kp
R
3

int
(r, ) = kp

1
r
2

1
R
3

cos
=
0

r
|
r
= R = kp
0

2
R
3

1
R
3

= 3
kp
0
R
3
cos
44 CAP

ITULO 4. ECUACI

ON DE LAPLACE
4.5. Coordenadas cilndricas y funciones de Bes-
sel
4.5.1. Soluci on de la ecuacion de LaPlace en coordenadas
cilndricas (Arm onicos Cilndricos)
En dos dimensiones
1
r

+
1
r
2

2
= 0
= R(r)P()
r
R
d
dr

r
dR
dr

=
1
P
d
2
P
d
2
= k
Entonces con P:
d
2
P
d
2
= kP
P =

cos

k
sen

k
Con R:
d
dr

r
dR
dr

=
kR
r
Si n 0:
R =

r
n
r
n
Si n=0:
R =

ln r
cte
Armonicos Cilndricos

ln r
r
n
cos nr
n
sen n
r
n
sen nr
n
cos n
En tres dimensiones

2
=
1
r

+
1
r
2

2
+

2

z
2
= 0
= RZ
z
r
d
dr

r
dR
dr

+
RZ
r
2
d
2

d
2
+ R
d
2
Z
dz
2
= 0
1
R
d
2
R
dr
2
+
1
Rr
dR
dr
+
1
r
2

d
2

d
2
=
1
Z
d
2
Z
dz
2
= b
2

z=Acosh bz+Bsenh bz
r
2
R
d
2
R
dr
2
+
r
R
dR
dr
+ b
2
r
2
=
1

d
2

d
2
= a
2

=C cos a+Dsen a
4.5. COORDENADAS CIL

INDRICAS Y FUNCIONES DE BESSEL 45

Con lo que se obtiene la ecuacion diferencial de Bessel
d
2
R
dr
2
+
1
r
dR
dr
+b
2

a
2
r
2
R = 0
R = EJ
a
(br) +FN
a
(br)
Donde:
J
a
(br) =

m=0
(1)
m
(
br
2
)
a+2m
m!(a +m+ 1)
N
a
(br) =
(cos a)J
a
(br) J
a
(br)
sen a
Si a=n y br
J
n
(x) =

2
x
cos(x

n
2
)
N
n
(x) =

2
x
sen(x

n
2
)
Ejemplo:
(r, ) = k

r
da = k

a sen 5
r
Rddz =
kaRl
r

sen 5d =
2kR

al

r
=
2kR
r
(r, ) = A
0
+A
1
ln r + (B
n
r
n
+C
n
r
n
) cos n + (D
n
r
n
+E
n
r
n
) sen n
Sir R,A
1
= C
n
= E
n
= 0

int
(r, ) = A
0
+B
n
r
n
cos n +D
n
r
n
sen n

int
(r, ) = A
0
+B
1
r cos +B
2
r
2
cos 2 + +D
1
r sen +D
2
r
2
sen 2 +
Sir R,A

1
= B

n
= D

n
= 0

ext
(r, ) = A

0
+C

n
r
n
cos n +E

n
r
n
sen n
46 CAP

ITULO 4. ECUACI

ON DE LAPLACE

ext
(r, ) = A
0
+
C
1
r
cos +
C
2
r
2
cos 2 + +
E
1
r
sen +
E
2
r
2
sen 2 +

ext
r


int
r

|
r=R
=
()

n=1

nC
n
R
n+1
cos n
nE
n
R
n+1
sen n

n=1

nB
n
R
n1
cos n +nD
n
R
n1
sen n

=
1

0
a sen 5

5E

5
R
6
sen 5 5D
5
R
4
sen 5 =
1

0
a sen 5
E

5
+D
5
R
1
0 =
1
5
0
aR
6

int
= A
0
+D
5
r
5
sen 5

ext
= A
0
+E
5
r
5
sen 5 = A
0
+

aR
6

0
D
5
R
10

r
5
sen 5

int
(r = R) =
ext
(r = R)
D
5
R
5
=
aR
5
0
D
5
R
5
2D
5
R
5
=
aR
5
0
D
5
=
a
10
0
r
4
E
5
=
aR
6
5
0

aR
6
10
0
=
aR
6
10
0
4.6. Soluciones especiales en dos dimensiones