Potencial en la ranura por medio de la ecuación de

Laplace y su gráfica en Python

D.M Jiménez Sarmiento, M zapata

Estudiante del curso de Teorı́a electromagnética 21 , Programa de Fı́sica, Universidad del Atlántico.
Orientador del curso de Teorı́a electromagnética; Francisco Racedo 22 , Universidad del Atlántico.

Mayo 21 del 2024

Resumen
En este proyecto se hallo por medio de la ecuación de laplace el potencial en la ranura de dos
placas, cuya solución se hizo por separación de varibales el cual es el métodos adecuado para la
solución de ecuaciones diferenciales parciales,haciendo todo el desarrollo matemático se hallo el
potencial en la ranura y con ayuda del lenguaje de programación python, se hicieron códigos los
cuales graficaron el comportamiento de el potencial en el contorno n=1,y el comportamiento del
potencial en la ranura con una gráfica 3D para n=1 y para n=47. estos códigos se corrieron y se
demostró que si funcionaba correctamente en la presentación.

1. introducción rrollo, luego las ecuaciones diferenciales parciales

se reducen a ecuaciones difrenciales ordinarias y
La ecuación de laplace. se le da ese nombre en se aplican las condiciones del problema llegando a
honor al astrónomo, fı́sico y matemático Pierre- la solución separables de la ecuación de laplace. se
Simon Laplace, ya que este hizo un gran aporte utiliza este método para para el problema que se
para el posterior descubrimiento de la transfor- tiene como proyecto, ya que El método de sepa-
mada de Laplace, y descubrió la ecuación de La- ración de variables es para encontrar una solución
place. esta ecuación es de segundo orden de tipo completa particular para ciertos problemas que in-
elı́ptico representa las ecuaciones diferenciales par- volucran ecuaciones en derivadas parciales. Este
ciales. la ecuaciones diferenciales ocupan un lugar es uno de los métodos más productivos de la fı́si-
en la solución de varios problemas en deferentes ca matemática para buscar soluciones a problemas
campos como electrostática, fluidos y otros mas, fı́sicos descritos mediante ecuaciones diferenciales
la la mayorı́a de problemas se pueden solucionar de derivadas parciales.
por medio de ecuaciones diferenciales, la solución
muchas veces no son exactas y son restringidas
a a través de condiciones de contorno y y pro-
2. Objetivo
blemas de valor inicial.Las soluciones de la ecua- Hallar el potencial en la ranura y hacer un codigo
ción de Laplace en un abierto de Rn . a partir de en phython, el cual muestre la gráfica del potencial
la ecuaciones diferenciales ordinarias y ecuaciones en esa ranura.
diferenciales parciales, surge la necesidad de los
métodos numéricos, una herramientas necesarias
para estas ya que es capaz de resolver estas ecua- 3. Planteamiento del problema
ciones. una forma para la solución de la ecuación
de laplace, es el método de separación de variables, Dos placas de metal conectadas a tierra permane-
para esto se tiene que tener en cuenta la geometrı́a cen paralelas en el plano XZ; en y=0 y y=a como
para elegir las coordenadas correcta para el desa- se muestra en la figura 1. El extremo en x se cie-

1
rra con una tira infinita aislada de las dos placas se puede escribir de esta manera:
y mantiene al potencial especificoV0 (y), encuentre
el potencial en la ranura.
d2 X(x) d2 Y(y)
Y(y) + X(x) (4)
dx2 dy 2

dividimos por la ecuación (3)

Y(y) d2 X(x) X(x) d2 Y(y)

+ (5)
X( X)Y( y) dx2 X( X)Y( y) dy 2

Figura 1: ilustración de las dos placas enunciadas

en el problema cancelando términos tenemos lo siguiente

4. Desarrollo 1 d2 X(x) 1 d2 Y(y)

+ (6)
X(x) dx2 Y(y) dy 2
Partiendo de que la configuración es independiente
de z es un problema entonces de dos dimensiones.
hacemos uso de la fórmula de Laplace: tenemos una solución de la forma f (x) + g(y) = 0
donde f (X) = −g(y) e igualamos a constantes
∇2 V = 0 (1) f(X) y g(y)

∂2 ∂2 ∂2
( + + )V = 0
∂x2 ∂y 2 ∂z 2
1 d2 X(x)
como es independiente de z entonces queda: = c1 (7)
X(x) dx2
∂2 ∂2
( + )V = 0
∂x2 ∂y 2
1 d2 Y(y)
= c2 (8)
∂ 2 V( x, y) ∂ 2 V( x, y) Y(y) dy 2
( + )=0 (2)
∂x2 ∂y 2

1) V=0 cuando y=0 donde la suma de las constante es cero porque

f (x) + g(Y ) = 0 entonces c1 + c2 = 0 en nues-
tro problema necesitamos una constante positiva
2)V=0 cuando y=a y una negativa, por ende: c1 = k 2 Y c2 = −k 2
entonces ahora tenemos:
3) V=V0 (y) cuando x=0

d2 X(x)
4) V→ 0 cuando x → ∞ = k 2 X(x) (9)
dx2

se supone una solución de la forma:

d2 Y(y)
V( X, Y ) = X( x)Y( y) (3) = −k 2 Y(y) (10)
dy 2
reemplazamos en la ecuación (2)
tenemos dos ecuaciones diferenciales ordinarias
∂ 2 [X(x) Y(y) ] ∂ 2 [X(x) Y(y) ] para que sea mas fácil de resolver, cuyas soluciones
+ =0 son:
∂x2 ∂y 2

2
4.1. Desarrollo de las ecuaciones dife- entonces
renciales X(X) = Aekx + Be−kx (15)
tenemos: y
Y(y) = Csin(ky) + Dcos(ky) (16)
d2 X(x)
− k 2 X(x) = 0 (11) son las solución de la separable apropiada de la
dx2 ecuación de Laplace.
donde:
se plantea
y = erk V( X, Y ) = X( x)Y( y) = (Aekx +Be−kx )(Csin(ky)+Dcos(ky))
(17)
y
d2 X(x) utilizando las condiciones planteadas anteriormen-
− k 2 X(x) = r2 − k 2
dx2 te, se tiene :
tenemos entonces una diferencia de cuadrados: utilizando la condición (4) se requiere que A = 0,
ya que V→ 0 cuando x → ∞
r2 − k 2 = (r − k)(r + k)
y se tiene que:
entonces r − k = 0

V(x,y) = e−kx (csin(ky) + Dcos(ky)) (18)

r+k =0
ahora haciendo uso de la condición (1), la cual dice
que V=0 cuando y=0, tenemos que D=0
tenemos: y1 = Aekx y y2 = Be−kx

V(x,y) = e−kx (csin(ky)) (19)

por ende la solución general es:
haciendo uso de la condición (2), la cual dice que
nπ
X(X) = Ae kx
+ Be −kx
(12) cuando V=0 y =a tenemos que K = a

para donde n toma valores 1,2,3...

d2 Y(y)
+ k 2 Y(y) = 0 (13)
dy 2
La ecuacion de Laplace es lineal por ende se tiene
se plantea V = V1 + V2 + V3 .... = 0, entonces V = C1 V1 +
y = erk C2 V2 + C3 V3 .... = 0

tenemos entonces: ∞
X nπx nπy
V(x,y) = Cn e a sin( ) (20)
2 2 a
r +k =0 n=0

entonces r2 = −k 2 utilizando la ecuación (3), la cual dice V=V0 (y)

cuando x=0, entonces
√
r= −1K, lo cual es ∞
X nπy
V(0,y) = Cn sin( ) = V0 (y) (21)
a
n=0
r = ±ik
utilizamos el truco de fourier, multiplicamos am-
bos lasos por sin( n′πy
a ) e integramos ambos lados.
r = a + ib entonces

y = eay (csin(by) + dsin(by) donde a=0, por ende ∞

X Z a
nπy n′πy
la solución general es: Cn sin( )sin( )dy =
0 a a
n=0
Ra
Y(y) = Csin(ky) + Dcos(−ky) (14) 0 V0 (y)sin( n′πy
a )dy(22)

3
por ortogonalidad sabemos que:
Z a
Z a
nπy n′πy 2V0 (y) n′πy 2V0 (y)
sin( )sin( )dy Cn = sin( )dy = (1−cos(nπ)
a a a 0 a nπ
0 (25)
̸ n
sera 0, si n′ =
entonces es 0 si n es par
a
2 si n′ = n cuando n′ = n queda la siguiente
expresión: 4V0 (y)
cuando n es impar nπ
Z a
a n′πy
Cn = V0 (y)sin( )dy (23)
2 0 a ∞
4V0 (y) X 1 −nπx nπy
V( X, Y ) = e a sin( ) (26)
ahora bien,¿ que pasa n- par o impar? π n a
n=1,2,3..
Z a
2 n′πy
Cn = V0 (y)sin( )dy (24) la forma del potencial en la ranura de las placas.
a 0 a

4
5. gráficos por el programador
para esto se decargaron los paquetes: import numpy as np y import matplotlib, pyplot as plt.

Figura 2

Figura 3: potencial en el contorno n=1

5
Figura 4: potencial en las placas paralelas n=1 en 3d y el 2d se ilustra en el mismo gráfico

Figura 5: potencial en las placas paralelas para n =47

Figura 6: potencial en las placas paralelas para n =1 hasta 101

6
6. conclusiones
se hizo el respectivo desarrollo matemático y analı́tico.la solución que se propone tiene que ser distinto
de cero porque si es igual a cero seria la solución trivial .

en la gráfica de contorno para n=1, la parte mas roja es donde se mide el potencial, en unas condición
se dice que el voltaje es cero cuando x tiende a infinito y eso nos muestra esa imagen como el potencial
va disminuyendo hasta llegar a cero.

en la gráfica 3D entre mayor sea el n el potencial va tendiendo a ser constante ya que el potencial
depende de x,y y n.

Referencias
[1] https://www.amazon.com/-/es/David-J-Griffiths-Introducci%C3%B3n-Electrodin%C3%
A1mica/dp/B081CYQ751
[2] https://link.springer.com/book/10.1007/978-3-319-22569-2
[3] [x]MamaniY.ChacaliazaJ.PelezA.LanchipaO.RamosJ.-Solucinanalticaynumricadelaecuacindelaplaceu
10.15381/rif.v24i2.20730