Mblopezs,+35221 137844 1 CE

"MOMENTO" Revista Depto. Fisica Univ ersidad Nacional de Colombia,Núm ero 7, Octubre 1992 .
Solución Numérica de la Ecuación de Laplace

en una Región Plana Acotada.
Mariela Mario (*)

Mauricio Garcia Castañeda (**)
Departamento de Fisica
Universidad Nacional de Colombia
RESUMEN: Empleando un algoritmo basado en el método de

Diferencias Finitas, se desarrolló una rutina numérica para
solucionar la Ecuación de Laplace en una región pl ana rectangula r.
As; mismo, se presentan gráficamente las li neas equipotenciales,
calculadas por medio de un esquema de in terpolación de Lagrange.
ABSTRACT: Laplace's equation for a flat rectangular region IS
solved using a numerical routine based on the fini te difference

method. AIso, the graphic representation of the equipotential lines
is presented by means of a Lagrange interpolation scheme.
INTRODUCCION.
La ecuación de Laplace pone de manifiesto que en una región en la cual

no existen cargas eléctricas, el potencial electrostático satisface una ecuación
diferencial parcial de segundo orden.
87
Consideremos la región rectangular de la figura 1. Asi, el problema a
solucionar es la ecuación diferencial parcial:
2 (p (p
\7 V(x'Y) =-2V(x'Y)+-2V(x,y)=O ( 1)
ax ay
con 0< x< aj 0< y < b, sujeta a unas condiciones de frontera que en general
pueden darse por:
v ( x , O ) = f(x) j V ( x , b) = g(x)
V (O , Y ) = p(y) j V ( a , y ) = q(y)
V=g(x)
b
V=p(y) \72 V ( x , y ) =O V=q(y)
V=f(x) a
Figura 1. Definición de la región plana rectangular

para la Ecuación de Laplace.
EL ALGORITMO NUMERICO.
Solucionar numéricamente la ecuación (1) por el método de las

diferencias finitas, implica reemplazar las derivadas presentes en el problema por
sus correspondientes cocientes diferenciables.
Sea la función V (x , y ). Al mantener fija y, la aproximación de la
88
primera derivada de V respecto a x en el intervalo x, x ± ~x , se expresa como
(1,2) :
(2)
Usando la mitad de los incrementos anteriores se puede conseguir la

aproximación para la segunda derivada:
02 V(x,y) :::::::--2
-2 1
(Vx+/::,x,y - 2 Vx y + Vx-~x y) (3)
fu (6~ "
De igual manera se pueden obtener las aproximaciones para los valores

numéricos de las derivadas parciales respecto ay.
Reemplazando (3) en la ecuación de Laplace y generalizando, se obtiene:
v-? V(x,y) (~lX)2 (V x+ óx ,y - 2 Vx,y + Vx_ óx,y)
+~ (Vx,y+ óy - 2 Vx y + Vx y_ óy ) (4)
(~y) "
Si dividimos la región en intervalos Ó x = t:..y = g, la ecuación (4) toma

la forma:
v 2 V (x,y) ::::::: 12 ( Vx+ 6 x,y + Vx,y+ /::, y - 4 Vx,y + V x- /::, x,y + V x,y- /::, y )
g
lo que permite escribir la siguiente expresión:
(5)
89
que proporciona de manera aproximada el potencial electrostático en cualquier
punto x,y dentro de la región, el cual está determinado por el promedio de los
valores en los cuatro puntos adyacentes a lo largo de los respectivos ejes.
U na vez se hayan obtenido los primeros valores para la rejilla, el proceso

se repite tantas veces como sea necesario hasta una última iteración en la cual las
diferencias entre los valores previamente calculados y los hallados, difieren menos
que un valor predeterminado.
CASOS ANALIZADOS
A. Variación Lineal del Potencial Electrostático Sobre las Fronteras.
Para la región R se adoptó la configuración mostrada en la figura 2, la

cual representa un potencial electrostático que varia de manera lineal sobre cada
una de las fronteras.
Al utilizar el programa desarrollado se calcula de manera iterativa los
potenciales electrostáticos punto a punto en el interior de la región aplicando
sucesivamente la ecuación (5). Luego al emplear la rutina de interpolación
mencionada, se encuentran las lineas equipotenciales las cuales se muestran en la
figura 3 para valores entre 1 y 10 voltios en incrementos de 0.5 V.
D. Potenciales Electrostáticos Constantes Sobre las Fronteras.
En este caso, las condiciones de frontera se establecen como aparece en la

figura 4.
Los resultados obtenidos para las equipotenciales se pueden apreciar en la
figura 5 para valores entre -5y 5 voltios con un intervalo de 0.5 V.
90
2 3 4 5 6 7 8 9 10
2 • • • • • • • • 9
3 • • • • • • • • 8
4 • • • • • • • • 7
5 • • • • • • • • 6 Figura 2. Condiciones
6 • • • • • • • 5 Linealmente Variables
T • • • • • • 4
8 • • • • • • • 3
9 • • • • • • • • 2
10 9 6 3
1 10
Jo J. Jo J.
2 9
3 8
\
\,
\
\ ....
---- - - - - - - - -
ó
10
Fig.3
Lineas equipotenciales de 1 a 10 V., para la configuración
91
most rada en la figura 2.
5 5 5 5 5 5 5 5 5 5
5 • • • • • • • • -5
5 • • • • • • • • -5
5 • • • • • • • • -5
5 • • • • • • • • -5
Fig. 4: Condiciones
5 • • • • • • • • -5
Constantes.
5 • • • • • • • • -5
5 • • • • • * • • -5
5 • • • • • • • • -5
.5 -.5 -.5 -5 -5 -5 -.5 -5 -5 -5
5 5 5 5 5 5 5 -5
5 -5 -5 -5 -5 -5 -5 -5 -5
-5
Fig.5
Lineas equipotenciales de -5 a +5 V., para las condiciones
92 de frontera mostradas en la figura 4.
(*) Estudiante del programa de Especialización en Ciencias Fisicas,
departamento de Fisica, Universidad Nacional de Colombia.
(**) Profesor Asociado, departamento de Fisica, Universidad Nacional de
Colombia.
REFERENCIAS.
1. Churchhouse R. , Numerical Methods, J. ''''iley, 1981
2. Merril M., Using Computers in Physics, Houghton Miff1ing, 1976 .
93

Mblopezs,+35221 137844 1 CE

Cargado por

Copyright:

Formatos disponibles

Mblopezs,+35221 137844 1 CE

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

Mblopezs,+35221 137844 1 CE

Cargado por

Copyright:

Formatos disponibles

"MOMENTO" Revista Depto. Fisica Univ ersidad Nacional de Colombia,Núm ero 7, Octubre 1992 .

Solución Numérica de la Ecuación de Laplace

Mariela Mario (*)

RESUMEN: Empleando un algoritmo basado en el método de

ABSTRACT: Laplace's equation for a flat rectangular region IS

solved using a numerical routine based on the fini te difference

La ecuación de Laplace pone de manifiesto que en una región en la cual

V=p(y) \72 V ( x , y ) =O V=q(y)

Figura 1. Definición de la región plana rectangular

Solucionar numéricamente la ecuación (1) por el método de las

Sea la función V (x , y ). Al mantener fija y, la aproximación de la

Usando la mitad de los incrementos anteriores se puede conseguir la

De igual manera se pueden obtener las aproximaciones para los valores

Reemplazando (3) en la ecuación de Laplace y generalizando, se obtiene:

v-? V(x,y) (~lX)2 (V x+ óx ,y - 2 Vx,y + Vx_ óx,y)

Si dividimos la región en intervalos Ó x = t:..y = g, la ecuación (4) toma

lo que permite escribir la siguiente expresión:

U na vez se hayan obtenido los primeros valores para la rejilla, el proceso

A. Variación Lineal del Potencial Electrostático Sobre las Fronteras.

Para la región R se adoptó la configuración mostrada en la figura 2, la

D. Potenciales Electrostáticos Constantes Sobre las Fronteras.

En este caso, las condiciones de frontera se establecen como aparece en la

1. Churchhouse R. , Numerical Methods, J. ''''iley, 1981

2. Merril M., Using Computers in Physics, Houghton Miff1ing, 1976 .

También podría gustarte