Universidad de Talca Escuela de Bioinformtica

Facultad de Ingeniera Claudio Tenreiro

Tpico Especial: Ecuacin de Poisson y Laplace

En el aula presencial se mostr un resultado importante, cual es la derivacin
de una ecuacin para el potencial electrosttico en funcin de la densidad de la
distribucin de la carga espacial. La derivacin resulta fundamentalmente de la
aplicacin de la Ley de Gauss para el flujo del campo elctrico, teniendo en
cuenta el teorema de la divergencia del clculo vectorial. Veamos esto
nuevamente.

El flujo del campo vectorial elctrico,
a travs de una superficie S cuya
forma es arbitraria, como la
mostrada en el dibujo y que
encierra una carga total q es:
0

q
da
S S
= = =

n E ds E

Evidentemente si se trata de una
distribucin de carga espacial y no
de una carga puntual q, tendremos
que integrar la carga total sobre el
volumen que ella ocupa, es decir:

=
V S
dV da

0
1
n E

En este caso S es la superficie que encierra al volumen V, el cual contiene la
carga. Esta expresin es la conocida forma integral de la ley de Gauss, la cual
podemos re escribir, en su forma diferencial, si hacemos uso del teorema de la
divergencia aplicndolo al campo vectorial electrosttico. El teorema de la
divergencia permite relacionar integrales de superficie y de volumen. Dicho lo
anterior tendremos entonces que:

= =
V V S
dV dV da

0
1
E n E y dado que el volumen considerado es
totalmente arbitrario o general, es posible deducir que hay una nueva relacin,
entre los argumentos de las integrales, involucrando esta vez slo un operador
diferencial, es decir

0
1
= E
Esta expresin es la que se denomina forma diferencial de la Ley de Gauss.
Ahora bien, si consideramos que para mover una cierta cantidad de carga en
presencia de un campo elctrico se requiere un trabajo, lo que es evidente ya
que el campo ejerce una fuerza sobre la carga que se desea mover, entonces
tendremos que dicho trabajo representa una energa, en este caso, una
energa potencial elctrica proveniente de la interaccin electrosttica entre las
cargas. Podemos denotar dicha energa como U, como una variacin de
energa potencial, la cual es posible evaluar a partir del clculo del trabajo
suponiendo que slo la fuerza elctrica es la nica que se ejerce sobre el
objeto cargado durante el desplazamiento, a saber:
q
S
E
n

da
Universidad de Talca Escuela de Bioinformtica
Facultad de Ingeniera Claudio Tenreiro


= = =
b
a
b
a
b
a
q U q W dr E dr E dr F y el potencial electrosttico se define
como la energa potencial elctrica por unidad de carga
q
U
= donde hemos
supuesto, arbitrariamente, que este potencial para distribuciones separadas
una distancia muy grande, de hecho una distancia infinita, es cero.
En este caso se obtiene una relacin bastante simple entre el potencial y el
campo elctrico, dada por una integral de lnea

= dr E
El valor de esta integral no depende de la trayectoria usada al mover la carga
en presencia del campo
1
, por ser este un campo de fuerzas conservativo
(teora de potencial) al igual que el caso de la interaccin gravitacional. Un
anlisis simple, en una dimensin de esta expresin, por ejemplo en x, nos
muestra algo notable, veamos;

= =
x x
E
dx
d
dx E


Es decir que si se conoce los valores del potencial electrosttico (campo
escalar) entonces es posible encontrar el campo elctrico (campo vectorial) por
aplicacin simple de una derivada. Este resultado se generaliza a la situacin
3-D mediante el uso del operador diferencial vectorial gradiente, el campo
apunta entonces en la direccin del gradiente y las lneas equipotenciales son
perpendiculares a l. En 3-D se tiene:
= E Donde en coordenadas cartesianas z
z
y
y
x
x

=
Luego, si el campo se expresa como el gradiente del potencial, = E ,
entonces, al substituir esta expresin para el campo en la forma diferencial de
la Ley de Gauss se obtiene una nueva ecuacin,
( )

0
2
0 0
1 1 1
= = = E
Esta ltima expresin se llama ecuacin de Poisson y el operador diferencial es
llamado Laplaciano.
En coordenadas cartesianas este operador es bastante simple,
2
2
2
2
2
2
2
z y x

=

Entonces tenemos ahora una manera de encontrar el
potencial en un problema electrosttico a partir de una ecuacin diferencial y se
distinguen dos situaciones, cuando la distribucin de carga es no nula, que
implica resolver la ecuacin de Poisson, y el caso en que se busca la solucin
en una regin del espacio libre de carga, en cuyo caso la ecuacin se llama de
Laplace.
) ( 0
) (
1
2
0
2
Laplace
Poisson
=
=

1
Como ejemplo suponga que desea mover un objeto cargado en el campo generado por una carga puntual
y verifique que esa integral se anula si el punto inicial y final de la trayectoria del recorrido es el mismo.
Universidad de Talca Escuela de Bioinformtica
Facultad de Ingeniera Claudio Tenreiro
z
r

Esto significa que disponemos de una alternativa adicional para calcular el
potencial electrosttico, no quiere decir que sea ms fcil, ya que estas son
ecuaciones diferenciales de segundo orden y requieren que sean especificadas
de manera correcta las condiciones iniciales o de contorno del problema.
Observacin: En general el problema ser bastante ms fcil si se identifican
simetras, por ejemplo que el potencial no deba depender de ciertas
coordenadas, eso reduce drsticamente la forma y complejidad de la ecuacin
diferencial resultante. Por ejemplo, el Laplaciano en coordenadas esfricas es:
2
2
2 2 2
2
2
2
1 1 1

=
sen r
sen
sen r r
r
r r


Con lo que si se trata de un problema que presente
simetra esfrica, slo el primer trmino es el que
interesa ya que el potencial resultante no puede
depender de los ngulos.




Finalmente, tambin es interesante el caso de coordenadas cilndricas (en este
caso r es la distancia perpendicular al eje del cilindro, ver figura):
2
2
2
2
2
2
1 1
z r r
r
r r

Ejemplos de este tipo de clculos:

Flujo de Campo:
Dada una expresin del campo, en forma vectorial, como por ejemplo
( ) ( ) k z j y i x A E

1

2
+ + = encontrar el flujo de ese campo por un disco de radio b,
en el plano x-y y en z=0 y orientado en la direccin k (podra ser la tapa
inferior de un cilindro de ese radio). Por definicin de flujo de un vector

= dS E
E
, entonces )

( k da dS = y para Z=0 k A k A E

) 1 ( = = por lo que el
integrando queda ( ) ( ) da A k da k A dS E = =

y con eso, la evaluacin del

flujo se hace como
2
b A da A da A dS E
disco disco
E
= = = =


r


Universidad de Talca Escuela de Bioinformtica
Facultad de Ingeniera Claudio Tenreiro
Trabajo:
Dado un campo de la forma j y E

= encontrar el trabajo al desplazar una carga

q
0
desde el punto A=(0,0,0) al punto B=(1,1,0) por medio de la trayectoria
0 ,
2
= = z x y . Esta es una trayectoria parablica en el plano x-y. En primer
lugar, el trabajo se define como una integral de lnea, es decir

B
A
B A
dr E q W
0

y el desplazamiento ocurre en el mismo plano de la trayectoria, es decir
j dy i dx dr

+ = por lo tanto el integrando queda ( ) ( ) dy y j dy i dx j y dr E = + =


pero la trayectoria vincula x e y, dx x xdx x ydy xdx dy x y
3 2 2
2 2 2 = = = =
Con lo cual la integral es bastante ms simple y queda como
0
1
0
4
0
1
0
3
0
2
1
2
1
2 q x q dx x q W
B A
= = =

Campo a partir de potencial:
Dado un potencial electrosttico de la forma

+
+ =
2 2
0
1 ) (
a x
x
V x se puede
encontrar el campo elctrico asociado gracias a la relacin entre ellos por
medio del gradiente. En este caso es un problema unidimensional, con lo que
la expresin general = E queda simplemente
dx
d
E
x

= en la direccin de
x (o componente x del gradiente en coordenadas cartesianas). Derivando se
obtiene que el mdulo del campo es
( )

+
=
2
3
2 2
2
0
a x
a
V E
x

Ecuacin de Laplace:
Hay un conjunto de casos en que resolver la ecuacin de Laplace es simple,
por ejemplo cuando las simetras indican que el problema slo depende de una
variable, por ejemplo en coordenadas cartesianas y el problema slo depende
de x, se le ocurre un ejemplo?, piense en el caso de placas paralelas.
constantes b a b x a x a
dx
d
dx
d
, , ) ( 0
2
2
+ = = =


Es decir se tiene un potencial que depende linealmente de x (o sea un campo
constante como en el caso de placas paralelas), sin embargo el problema no es
tan simple, como ya fue mencionado, ya que es necesario incorporar las
condiciones de borde del problema, por ejemplo supongamos que
2 1
2 1 1 2
2 1
2 1
2 2
1 1
2 2
1 1
,
) (
) (
x x
x x
b
x x
a
b x a
b x a
x
x

=
+ =
+ =

=
=

Con esto es posible entonces escribir la solucin general para el potencial
como
( ) ( ) { }
1 2 2 1
2 1
1
) ( x x x x
x x
x

=


Universidad de Talca Escuela de Bioinformtica
Facultad de Ingeniera Claudio Tenreiro
x
d

Supongamos que tenemos una
configuracin de dos placas separadas
una distancia d, la de la izquierda en el
origen, es decir en x
1
=0, y supongamos
que el potencial en ese punto es cero y
en x
2
=d es V
0
, entonces
x
d
V
x
V d x
x
0
0 2 2
1 1
) (
,
0 , 0
=
= =
= =

Y el campo elctrico ser, entonces la
derivada de este potencial, es decir
x
d
V
dx
x d
E
x

) (
0
= =




Otro caso simple es en coordenadas cilndricas, donde el problema es resolver
la ecuacin
b r ln a r a
dr
d
r
dr
d
r
dr
d
r
r
r r
+ = = =

) ( ) ( 0 0
1

a y b se obtienen de las condiciones de borde, de manera semejante al caso
anterior. Por ejemplo se podra especificar el potencial en dos puntos dados por
radios diferentes,

+ =
+ =

= > =
= =
b r ln a
b r ln a V
r r r
V r r
) ( 0
) (
0 ) (
) (
2
1 0
1 2
0 1

De estas ecuaciones se despeja a y b para encontrar el potencial en
cualquier punto r del espacio.

Esta es la tcnica usual para obtener la solucin de una distribucin de cargas
en el espacio (como parte de un objeto de geometra definida) casi en cualquier
condicin. La fortaleza de la tcnica radica que en especial que hay mtodos
numricos para resolver estas ecuaciones diferenciales con cualquier tipo de
condicin de borde.

Recomendado: Busque y aplique soluciones de estos problemas, disponibles
en la web, usando algn software de su gusto y conocimiento, como MathCad,
Matlab, Mathematica o el que sea que Ud. domine.
Discuta la convergencia de las soluciones y restricciones de las mismas en
casos generales.
Aplique el mtodo a problemas que Ud. haya reconocido como de inters en el
mbito de los polmeros y monmeros u otras molculas complejas pero en
aproximaciones de geometras simples, como lneas o curvas, as como
membranas, las cuales pueden ser representadas por figuras como planos.