Universidad de las Fuerzas Armadas ESPE

Extensi on Latacunga
Matem atica Superior

Ecuaciones de Laplace

Kevin Ivan Barrera Llanga y Christyan Mario Cruz Ulloa

4o B Ingenier a Mecatr onica

15 de diciembre de 2013

 Indice
I Objetivos 2
2 2

1. Objetivo General 2. Objetivos Espec cos

II III IV V VI VII

Introducci on Marco teorico Ejercicio Resuelto Ejercicio de Aplicaci on en Matlab Ejercicio Propuesto Concluciones Bibliograf a

2 3 5 8 11 11 11

VIII

Parte I

Objetivos
1. Objetivo General
Determinar los conceptos y las ecuaciones de Laplace en su forma Cil` ndrica, Esf` erica y Rectangular.

2.

Objetivos Espec cos

c) Comprobar el comportamiento te orico del potencial electroest atico en el interior de un sistema de placas a partir de gr acas simuladas en matlab.

Parte II

Introducci on
La ecuaci on de Laplace describe siempre procesos estacionarios en los que el tiempo no es una de las variables independientes. Esto signica, que el tipo de problemas que se considera para esta ecuaci on no va a incluir condiciones iniciales, lo que provocar a que la estructura de las soluciones que encontraremos sea distinta a la que presentan las ecuaciones de ondas y difusi on. Las ecuaciones el pticas se obtienen, en general, cuando se estudian procesos estacionarios. La ecuaci on de este tipo que aparece m as frecuentemente es la ecuaci on de Laplace, cuya expresi on general es: u = 0 donde representa el operador Laplaciano.Otra clase de problemas en los que la ecuaci on de Laplace juega un papel importante es el relacionado con los campos vectoriales que derivan de un potencial (electro-magn etico, gravitatorio, etc.). En estos casos, si en la regi on donde se quiere obtener el potencial no hay fuentes de campo, este verica la ecuaci on de Laplace. Si por el contrario, existe una cierta densidad de estas, el potencial verica la correspondiente ecuaci on de Poisson.

Parte III

Marco teorico
La ecuaci on de Laplace es una ecuaci on en derivadas parciales de segundo orden de tipo el ptico, que recibe ese nombre en honor al f sico y matem atico Pierre-Simon Laplace. Introducida por las necesidades de la mec anica newtoniana, la ecuaci on de Laplace aparece en muchas otras ramas de la f sica te orica como la astronom a, la electrost atica, la mec anica de uidos o la mec anica cu antica. Para obtener la ecuaci on de Laplace partimos del operador Nabla. d d d + y + z dx dy dz

= x

Laplaciano es igual al producto escalar de nabla:

. . .

= = (i

.(

V)

d d d d d d +j + k ).(i +j +k ) dx dy dz dx dy dz d2 V d2 V d2 V = + + dx2 dy 2 dz 2 3

Igualando a cero el operador nabla obtenemos la ecuaci on de laplace que en coordenadas cartesianas la podemos escribir como:
2

d2 V d2 V d2 V + + =0 2 2 dx dy dz 2

E. Laplace en coordenadas Rectangulares

d2 u d2 u d2 u + 2 + 2 =0 dx2 dy dz Podemos ampliar esta ecuaci on en coordenadas cil ndricas y coordenadas esf ericas, para esto es necesario los siguientes cambios de variables: E. Laplace en coordenadas Cil ndricas

x = r.cos y = r.sen z=z

1 d2 u 1 d du d2 u . 2 + . (r ) + 2 = 0 2 r d r dr dr dz E. Laplace en coordenadas Esf ericas

x = p.sen.cos y = p.sen.cos z = p.cos 1 d 1 d du 1 d2 u 2 du . . ( p ) + ( sen. ) + . =0 p2 dp dp p2 sen d d p2 sen2 d2

Parte IV

Ejercicio Resuelto
Fuente: Ecuaciones Diferenciales-Deniss Zill, Cap. 11, Vol. 1, P ag. 545. C alculo de la temperatura u en una placa rectangular

Suponga que deseamos encontrar la temperatura de estado estable u(x, y ) en una placa rectangular cuyas orillas verticales x = 0x = a se encuentran aisladas, mientras las orillas superior e inferior y = by = 0 se mantienen a temperaturas f (x) y 0, respectivamente. Consulte la gura cuando no escapa calor desde las supercies laterales de la placa, resolvemos el siguiente problema de valores en la frontera: d2 u d2 u + 2 = 0, 0 < x < a, 0 < y < b, (1) dx2 dy du du |x = 0 = 0, |x = a = 0, 0 < y < b, (2) dx dx u(x, 0) = 0, u(x, b) = f (x), 0 < x < a, (3) Soluci[on del problema de valores en la frontera. Con u(x, y ) = X (x)Y (y ), la separaci on de variables en (1) conduce a X = Y Y = X X + .X = 0, (4) Y + .Y = 0, (5) En (2) y (3), las tres condiciones de frontera homog eneas se traducen en X (0) = 0, X (a) = 0 y Y (0) = 0. El problema de Sturm-Liouville asociado con la ecuaci on (4) es entonces X + .X = 0, X (0) = 0, X (a) = 0 (6) El an alisis de los casos correspondientes a = 0, = 2 < 0y = 2 > 0, donce > 0, Por comodidad, a continuaci on presentamos una versi on sintetizada de dicho an alisis. Para = 0, (6) se convierte en X = 0, X (0) = 0, X (a) = 0 La soluci on de la ecuaci on diferencial ordinaria es X = C 1 + C 2x. La condici on de frontera X (0) = 0 entonces, implica que c2 = 0, por lo que X = c1. Observe que para cualquier c1, esta soluci on constante satisface la segunda condici on de frontera X (a) = 0. Haciendo que c1 = 0, X = c1 es una soluci on no trivial del problema de valores en la frontera (6). Para = 2 < 0, (6) no tiene una soluci on no trivial. Para = 2 > 0, (6) se convierte en X + 2 X = 0, X (0) = 0, X (a) = 0

Al aplicar la condici on de frontera X (0) = 0, la soluci on X = c1.cos(.x) + c2.sen(.x) implica que c2 = 0, por lo que X = c1.cos(.x). La segunda condici on de frontera X (a) = 0 aplicada a esta u ltima expresi on nos da entonces c1.sen.a. Debido a que > 0, la u ltima ecuaci on se satisface cuando .a = n. 0 = n , n = 1, 2, ...Los valores propios de (6) son entonces 0 y a 2 2 2 n = n = n 2 , n = 1, 2, ... Por la correspondiente 0 = 0 con n = 0, las a funciones propias de (6) son X = c1, n = 0 n X = c1.cos x, n = 1, 2... a Ahora debemos resolver la ecuaci on (5) sujeta a la u nica condici on de frontera homog enea Y (0) = 0. Primero, para 0 = 0, la ecuaci on diferencial en (5) es simplemente Y = 0 y, por lo tanto, su soluci on es Y = c3 + c4y . Sin embargo, Y (0) = 0 implica que c3 = 0, en consecuencia, Y = c4y . Segundo, para n2 2 2 on diferencial en (5) es Y 2 Y = 0. Como 0 < y < b n = n2 2 , la ecuaci a a es un intervalo nito, escribimos la soluci on general en t erminos de las funciones hiperb olicas: y y Y = c3.cosh(n. ) + c4senh(n ) a a A partir de esta soluci on podemos observar que Y (0) = 0 de nuevo implica y c3 = 0, en consecuencia Y = c4senh(n ). Las soluciones producto un = a X (x)Y (y ) que satisfacen la ecuaci on de Laplace (1) y las tres condiciones de frontera homog eneas dadas en (2) y (3) son A0 y, n = 0 n An senh a n cos x, n = 1, 2.... a donde hemos escrito nuevamente c1c4 como A0 para n = 0 y como An para n = 1, 2... El principio de superposici on da otro resultado

u(x, y ) = A0 y +
n=1

An senh

n n ycos x a a

(7) Por u ltimo, sustituyendo y = b en (7) observamos que u(x, b) = f (x) = A0 b + n n b.cos b n=1 An senh a a Donde A0 = f (x)dx 2 n a An = n 0 f (x)cos a xdx a.senh b a 7 1 ab
a 0

Parte V

Ejercicio de Aplicaci on en Matlab

en donde las placas paralelas al plano XY son placas semiinnitas con potencial igual a 0 y la placa en el plano XZ se encuentra un potencial b0 . Las condiciones de frontera son: a) En y=0 el potencial el ectrico es: V (y = 0, z ) = V0 b) Cuando la variable y tiende al innito el potencial electrico es: V (y , z ) = 0 c) En z=0 el potencial el ectrico es V (y, z = 0) = 0 d) En z=b el potencial el ectrico es V (y, z = b) = 0 Soluci on general: V (y, z ) = Aey + Bey Csen(z ) + Dcos(z ) Programacion en matlab

Archivo algoritmo.m

Condiciones de la frontera Archivo f1.m

Archivo f2.m

Archivo f3.m

Archivo f4.m

Programa que se ejecuta Archivo corre.m

La respuesta de este nos da un valor numerico en matlab que gracamente se representa como:

Donde nos indica que el valor nal es 1.4:

10

Parte VI

Ejercicio Propuesto
Ecuaciones diferenciales Dennis Zill, ejercicios 11.5, pag 549 Resuelva la ecuaci on de Laplace para una placa rectangular sujeta a condiciones inciales de frontera que se proporciona. 1. u(0, y ) = 0, u(a, y ) = 0 u(x, 0) = 0, u(x, b) = f (x)

Parte VII

Concluciones
1. En la vida Real la Ecuaci on de Laplace se puede aplicar en el an alisis del comportamiento te orico del potencial electrost atico en el interior de un sistema de placas. 2. Se pueden resolver ecuaciones de Laplace en el software Matlab, seg un los comandos indicados.

11

Parte VIII

Bibliograf a
Referencias
[1] http://hyperphysics.phy-astr.gsu.edu/hbasees/electric/ laplace.html [2] http://www.youtube.com/watch?v=uj3NO3ydCtE [3] http://es.scribd.com/doc/33933327/Ecuacion-de-Laplace-en-Matlab [4] http://www.tecnun.es/asignaturas/metmat/Texto/En_web/ Ecuacion_de_Laplace/Ecuacion_de_Laplace.pdf [5] http://pitagora.usach.cl/~vguinez/apuntes/cap12.pdf

12