Resumen Edp

Descargar como docx, pdf o txt
Descargar como docx, pdf o txt
Está en la página 1de 24

Captulo 7

Ecuacin de Laplace
En el captulo anterior estudiamos los problemas valor inicial y valor de lmite inicial para las
ecuaciones de la onda y de calor. Vimos que estos problemas son bien representados por la
ecuacin de la onda y de calor en el sentido de que los problemas tienen soluciones nicas, las
cuales dependen continuamente de los datos iniciales. En este ltimo captulo, tenemos que
tener cuidado con ciertos problemas que son bien representados para las ecuaciones de Laplace
en dos variables independientes.

2 2
+ =0
2 2

Esta ecuacin sirve como el prototipo para las ecuaciones diferenciales elpticas de la misma
forma que las ecuaciones de la onda y calor lo hacer para las ecuaciones diferenciales
hiperblicas y parablicas, respectivamente.

1. Problemas de Valor Limite

En orden de ver como las ecuaciones de Laplace pueden surgir en un problema fsico,
consideremos la distribucin de temperatura en una delgada placa uniforme ocupando un
dominio D en el plano (Fig. 7.1). Asumimos que D tiene un lmite continuo C. Haciendo
suposiciones similares a esos que se hicieron derivando ecuaciones de flujo de calor
unidimensional, estos pueden ser mostrados como la distribucin de temperatura en la placa
cumpliendo la ecuacin homognea en D.

2 2
(1.1) ( 2 + 2) = 0

Fig. 7-1 Dominio de plano y curva lmite

Esta ecuacin esa llamada la ecuacin de calor bidimensional. Suponiendo que el flujo de calor
en D es continuo: entonces la distribucin de temperatura es independiente del tiempo (la
cual es, / = 0) y por lo tanto la ecuacin (1.1) se convierte en:

2 2
(1.2) + =0
2 2
La cual es la ecuacin de Laplace bidimensional.

Si hay una fuente de calor en D, la funcin de la temperatura va a satisfacer la ecuacin


de calor no homognea:

2 2
(1.3) ( 2 + 2 ) = (, , )

Consecuentemente, cuando el flujo de calor es continuo entonces no depende de , (1.3) se


reduce a la ecuacin de Laplace no homognea.

2 2 (, )
(1.4) + = , (, ) =
2 2
La ecuacin (1.4) es conocida como la ecuacin de Poisson. Ambas ecuaciones, (1.2) y (1.4), son
ejemplos prominentes de ecuaciones diferenciales elpticas que son de gran importancia en
fsica matemtica. A menudo escribimos estas ecuacin en la forma compacta = 0 y =
, donde denota el operador de Laplace:

2 2
= +
2 2

De la discusin precedente, parece que en un estado estacionario, la distribucin de la


temperatura en un cuerpo conductor de calor uniforme satisface la ecuacin de Laplace.

Correspondientemente, cuando hay una fuente de calor conocida, la funcin de temperatura


satisface la ecuacin de Poisson. A partir de esta consideracin fsica, parece natural considerar
problemas de valores lmite para el operador de Laplace en lugar de problemas de valores
iniciales o valores de lmite inicial. De hecho, los problemas de los valores lmite son tpicos para
las ecuaciones diferenciales elpticas porque los problemas de valor inicial o de valor lmite inicial
no estn bien presentados para este tipo de ecuacin. A este respecto, examinemos un ejemplo
ideado por Hadamard para la ecuacin de Laplace. Considere el siguiente problema de valor
inicial:
+ = 0 ( < < , > 0)
(1.5) (, 0) = 0 ( < < )
(, 0) = sin ( < < )

Por el mtodo de separacin de variables, obtenemos la solucin explicita de este problema.

(1.6) (, ) = sin sinh


Ahora observamos que el parmetro se convierte en infinito, el dato inicial del problema
tiende a cero, mientras que para 0, la solucin (1.6) se vuelve arbitrariamente larga. Esto
muestra que la solucin del problema no depende continuamente del dato inicial, y por lo tanto
el problema no est bien fundamentado.

Un valor lmite tpico para una ecuacin diferencial elptica consiste en encontrar una
solucin de la ecuacin diferencial en un dominio D que satisfaga ciertas condiciones de
contorno en C de D. En este captulo, estaremos conteniendo problemas de valor lmite para la
ecuacin de Laplace con condiciones de contorno que son tpicas para los operadores
diferenciales elpticos. Nuestra discusin ser limitada a el plano bidimensional, si bien
muchos de estos resultados sern obtenidos y los procesos usados pueden ser inmediatamente
extendidos al caso de tres dimensiones o tridimensional.

El problema de valor lmite de primer tipo para la ecuacin de Laplace consiste en encontrar
una solucin de la ecuacin en un dominio D, el cual asume valores en los lmites C de D. Esto
es:
= 0
(1.7)
=
Donde es una funcin conocida en C. Este tipo de problema es llamado un problema de
Dirichlet, y la condicin lmite es siempre llamada condicin Dirichlet o condicin lmite del
primer tipo. Fsicamente, el problema puede ser interpretado como encontrado la distribucin
de temperatura de equilibrio en D cuando una distribucin de temperatura arreglada est en el
lmite de D.

El problema de segundo tipo para la ecuacin (1.2) consiste en encontrar una solucin de la
ecuacin en el dominio D satisfaciendo la condicin lmite:

(1.8) =

En el lmite C de D. Aqu / la derivada normal de en C. Este problema es llamado el
problema de Neumann, y la condicin (1.8) es llamada condicin o condicin lmite de segundo
tipo de Neumann. Por ejemplo, cuando determinamos la temperatura de equilibrio en un
cuerpo uniforme D para el cual una suma dada de calor es suplida al lmite C de D.

El problema de valor lmite del tercer tipo tiene la siguiente condicin:



(1.9) + =

En el lmite C donde y son funciones conocidas. Como un problema llamado el problema de
Robin o problema de los valores limites mesclados. La condicin lmite (1.9) sugiere cuando, para
algunos casos, aceptamos por radicacin de calor desde el lmite del cuerpo dentro del medio
circundante.

Esto es de inters al observar que las condiciones de contorno (1.7), (1.8) y (1.9)
corresponden a la condicin lmite del tercer tipo, encontramos en conexin con nuestra
discusin de la cuerda vibrante y la conduccin de calor en la vara delgada.

Problemas anlogos puedes, por supuesto, ser considerados por las ecuaciones de Poisson
(1.4). Siempre, como sea, ser posible reducir un problema usando las ecuaciones de Poisson a
uno que involucre la ecuacin de Laplace. El proceso consiste en primero encontrar una solucin
particular de la ecuacin de Poisson = y = . Luego, satisface la ecuacin de
Laplace y su condicin limite correspondiente que envolver . La funcin devera ser tal que
esta es definida y diferenciada en el lmite del dominio en consideracin. En el caso general, el
problema es la mejor aproximacin las funciones certeras Kernel, de las cuales se habla en la
Seccin 9 de este captulo.

Cada uno de los diversos problemas antes mencionados se denomina problema de valor
lmite interior y su solucin es requerida para satisfacer la ecuacin diferencial dentro del
dominio D. Si, en cambio, buscamos una solucin de la ecuacin diferencial fuera del dominio D
con condicin de frontera en C, entonces el problema se denomina problema de valor lmite
exterior. En este trabajo nos ocuparemos principalmente de los problemas interiores, aunque
en los ejercicios se presentarn algunos problemas exteriores.

2. Teorema de Green y singularidad de las soluciones

En esta seccin se tratar la cuestin de la singularidad de las soluciones del problema de


Dirichlet, el problema de Neumann y el problema de Robin para la ecuacin de Laplace, con la
ayuda de la identidad de Green. Consideremos primero el problema de Dirichlet (1.7).
Supongamos que 1 y 2 son dos soluciones de este problema y establecemos = 1 2 .
Est claro que:

(2.1) = 1 2 = 0
Y
(2.2) = =0

Vamos a considerar la primera identidad de Green sobre el dominio D:


(2.3) = ( + + )

Reemplazando = en (2.3), tenemos:



(2.4) = ( + 2 + 2 )

Si sustituimos la funcin = 1 2 por en (2.4) y usando (2.1) y (2.2) obtenemos:

(2.5) (2 + 2 ) = 0

Dado que el integrando en (2.5) es continuo y no negativo, concluimos que 2 + 2 = 0


idnticamente en D. Esto implica que debe ser constante en D + C. Pero = 0 en C. Por lo
tanto, resulta que = 0 identicamente en D + C: esto es 1 = 2. As, hemos propuesto nuestro
primer teorema en esta seccin.

Teorema 2.1. Una solucin del problema de Dirichlet (1.7) est determinada de manera nica.

En el caso o el problema de Neumann



(2.6) = 0 , =

Si 1 y 2 son dos soluciones, entonces su diferencia = 1 2 ser una solucin del


problema homogneo:

(2.7) = 0 , =0

Aplicando la identidad (2.4) a la funcin y usando (2.7), se nos lleva de nuevo a la
conclusin de que debe ser constante en D + C: es decir (, ) = . En este caso, sin
embargo, no puede establecerse que = 0, puesto que no desaparece necesariamente en
C. Por lo tanto, solo podemos tener:

1 = 2 +
Que dice que cualquier dos soluciones del problema (2.6) difieren por una constante.
Declaramos este resultado como nuestro segundo teorema.

Teorema 2.2. Una solucin del problema de Neumann (2.6) es nica hasta una constante aditiva.

Finalmente, para el problema de Robin



(2.8) = 0 , + =

Se puede establecer la unicidad de la solucin si asumimos que 0 y no es idnticamente
cero en C. De hecho, supongamos 1 y 2 son dos soluciones del problema (2.8). Entonces =
1 2 satisface la condicin

+ = 0

O

(2.9) =

Por lo tanto, por la identidad (2.4), tenemos:

(2.10) 2 = (2 + 2 )

Desde que 0, el termino en la izquierda de (2.10) no es positivo y el termino a la derecha no


es negativo. Adems, ambos trminos deben ser igual a 0, entonces:

(2 + 2 ) = 0

Lo que implica que = . en D+C. De (2.9) se sabe que = 0 en D+C; esto da 1 = 2 .


Esto prueba el prximo teorema.

Teorema 2.3. Una solucin del problema de valor de lmite mixto (2.8) con 0 y no es
idnticamente a cero el cual se determina de forma exclusiva.

Es importante sealar que la identidad de Green (2.3) tambin nos permite establecer una
condicin necesaria para la existencia de una solucin del problema de Neumann (2.6). Si
intercambiamos el papel de y en (2.3), y ponemos = 1, obtenemos.


(2.11) =

As, si el problema (2.6) tiene una solucin , entonces la aplicacin de (2.11) produce.

(2.12) = = 0

Esto demuestra que si el problema (2.6) tiene una solucin, entonces la integral de sobre el
lmite C debe desaparecer. En trminos de nuestra interpretacin fsica del problema, la
condicin (2.12) simplemente dice que en un flujo de calor en estado estacionario el suministro
neto de calor en el lmite de D debe ser cero.

En el caso en el que el problema (2.6) involucra la ecuacin de Poisson, = en D, la


aplicacin de (2.11) da la condicin necesaria

(2.13) (, ) =

para la existencia de una solucin de tal problema. En consecuencia, si la condicin (2.12) o


(2.13) no se cumple, entonces el problema de Neumann para la ecuacin de Laplace o la
ecuacin de Poisson no puede tener una solucin.

3. Principio mximo para funciones armnicas


Una funcin se denomina funcin armnica en un dominio D si tiene derivadas continuas
de segundo orden y satisface la ecuacin de Laplace en cada punto de D. En esta seccin
probaremos un principio mximo para las funciones armnicas, que es muy similar al principio
mximo que probamos para las funciones que satisfacen la ecuacin del calor.

Teorema 3.1. (Principio mximo) Sea una funcin armnica en un dominio delimitado D con
lmite C. Si es continua en D + C, entonces el valor mximo de sobre D + C se alcanza en el
lmite C de D.

Prueba: Demostramos primero que cualquier funcin continua en D + C y que satisface la


desigualdad diferencial
(3.1) + > 0

debe alcanzar su valor mximo en el lmite C. Sea una funcin continua en D + C que satisface
(3.1) y supongamos que u Asumi su valor mximo en un punto interior (0 , 0 ) de D. Entonces,
en tal punto, tendramos.

(0 , 0 ) 0 (0 , 0 ) 0

Entonces

(0 , 0 ) + (0 , 0 ) 0

Esto contradice el hecho de que satisface (3.1) en D. Por lo tanto, una funcin continua en D
+ C y satisfaciendo (3.1) alcanza su valor mximo en C.

Supongamos ahora que es una funcin armnica en D y contina en D + C. Sea M su valor


mximo en C y considere la funcin:

(, ) = (, ) + ( 2 + 2 )
Donde > 0 es un nmero arbitrario. Es claro que es continuo en D + C y

= + 4 = 4 > 0
Por lo tanto, por el resultado anterior, el mximo valor de es alcanzado en C. As, en cualquier
punto en D, tenemos:

(, ) max (, ) +

Donde = max( 2 + 2 ) para todo (, ) en D + C. Desde (, ) = (, ) ( 2 + 2 )


(, ), esto conlleva a:

(, ) (, ) +
Para todo (, ) en D + C, y para cualquier numero > 0. Hacemos que 0, entonces
encontramos (, ) , el cual estable el teorema.

Aplicando el principio mximo a la funcin armnica = , obtenemos un resultado


anlogo para el valor mnimo de . As, una funcin armnica en un dominio limitado que es
continua a lo largo del cierre del dominio alcanza sus valores mximos y mnimos en el lmite de
ese dominio. Por supuesto, si es una constante. Su mximo (que coincide con su mnimo) se
alcanza en cada punto de D. De hecho, una constante es la nica funcin armnica que puede
alcanzar un mximo en un punto interior de un dominio. Limitada o sin lmites. Esta es la versin
fuerte del principio mximo. Que estableceremos en la Seccin 6.

Usando el principio mximo para las funciones armnicas, es posible dar una prueba
alternativa del Teorema 2.1 para la unicidad de la solucin del problema de Dirichlet para la
ecuacin de Laplace. De hecho, si 1 y 2 son dos soluciones del problema (1.7), entonces su
diferencia = 1 2 es una funcin armnica que desaparece en el lmite C. Entonces, por el
principio mximo y mnimo, no puede ser mayor o menor que cero en D. Por lo tanto, debe
desaparecer idnticamente en D + C; Es decir, 1 = 2 .

El principio mximo tambin nos permite establecer una dependencia continua de la solucin
del problema de Dirichlet para la ecuacin de Laplace en los datos de frontera. Declaramos este
resultado como un teorema.

Teorema 3.2. Si y son dos funciones armnicas en un dominio limitado D con lmite C tal que
| | en C, luego:

| |
La prueba es simple y se deja como un ejercicio.

En las siguientes secciones, se determinarn soluciones del problema (1.7) en dominios


rectangulares y circulares. Por lo tanto, para estos tipos de dominios, los Teoremas 3.1 y 3.2
implican que el problema de Dirichlet est bien planteado para la ecuacin de Laplace.

4. Problema de Dirichlet en un Rectngulo

A partir de esta seccin consideraremos las soluciones de los problemas de Dirichlet y


Neumann para la ecuacin de Laplace en dominios simples tales como un rectngulo y un disco.
Se observar que para este tipo de dominio, el mtodo de separacin de variables puede
utilizarse para obtener una solucin del problema de Dirichlet o Neumann. Consideremos
primero en esta seccin el problema de Dirichlet en un dominio rectangular R (Fig. 7.2) 0 < <
, 0 < < . Buscamos una solucin de la ecuacin de Laplace.

(4.1) + = 0

La cual satisface las condiciones de contorno

(0, ) = 0 (, ) = 0 (0 )
(4.2)
(, 0) = 0 (, ) = () (0 )

Este problema puede ser interpretado fsicamente como encontrar la temperatura de equilibrio
en todo un cuerpo rectangular R delgado cuando sus lados = 0, = , = 0 se mantienen a
temperatura cero y el lado = se mantiene a una distribucin dada de temperatura (), 0
.
Fig. 7-2 El problema de Dirichlet en un Rectngulo

De acuerdo con el mtodo de separacin de variables, asumimos soluciones no triviales de (4.1)


de la forma:

(4.3) (, ) = ()()
La cual satisface las tres condiciones de contorno homogneas en (4.2). Sustituyendo (4.3) en
(4.1) y separando variable, obtenemos:

= =

Esto conduce a las dos ecuaciones diferenciales ordinarias.

(4.4) + = 0 = 0
Que deben ser satisfechas por las funciones X e Y, respectivamente. Para que (4.3) satisfaga las
tres condiciones de frontera homogneas para (4.2), debemos tener:
(4.5) (0) = 0 () = 0 (0) = 0
Por lo tanto, la funcin X e Y debe determinarse a partir del problema de valor propio:
(4.6) + = 0 (0) = 0 () = 0
Y el problema de valor inicial:
(4.7) = 0 (0) = 0
Cuando encontramos el problema (4.6) antes, encontramos los valores propios =
/2 con las correspondientes funciones propias:
2 2


() = sin ( = 1,2, )

Para cada = 2 2 /2, una solucin de (4.7) es dada por:
1 (/)
() = [ (/) ]
2

() = sinh ( = 1,2, )

Esto sigue la funcin:

(, ) = sin sinh ( = 1,2, )

Son todas soluciones particulares de la ecuacin de Laplace que se desvanecen en = 0, = ,
e = 0. Estas funciones se llaman armnicos rectangulares. Para obtener una solucin del
problema (4.1), (4.2) ahora formamos la serie


(4.8) (, ) = sin sinh

=1

donde los coeficientes son constantes an por determinar. Definiendo formalmente = y


usando la condicin de contorno en (4.2), obtenemos:

() = (, )


(4.9) = ( sinh ) sin

=1


= ( ) sin

=1

Donde:

(4.10) = sinh ( = 1,2, )

La ltima expresin en (4.9) es la expansin en serie de Fourier de sobre (0, ), y por lo tanto:
2
(4.11) = () sin ( = 1,2, )
0

Sustituyendo de (4.10) en (4.8), obtenemos:



sinh
(4.12) (, ) = sin

=1 sinh

Para la solucin del problema (4.1), (4.2). Los coeficientes estn dados por (4.11).

Si asumimos que es continua, por partes lisa en 0 , y se desvanece en = 0 y =


, podemos verificar que (4.12) da la solucin de nuestro problema (4.1), (4.2) por el mismo
mtodo utilizado en la Seccin 4 del Captulo 6. Observe que la desigualdad:

sinh / /
= /
/
sinh

1 2/ ()/
= ()/
1 2/ 1 2/
Puede utilizarse para probar la convergencia uniforme de la serie (4.12) y las derivadas parciales
de segundo orden de esa serie con respecto a e para 0 , 0 0 , con cualquier
0 < .

El problema general de Dirichlet

+ = 0

(, 0) = 1 () (, ) = 2 () (0 < < )
(4.13)
(0, ) = 3 () (, ) = 4 () (0 < < )

En el dominio rectangular R puede ser solucionado por una superposicin de las cuatro
funciones 1 , 2 , 3 , y 4 , donde cada (1 4) denota la solucin de (4.13) donde toda
dato lmite, excepto son cero.

Ejemplo 4.1 Encontrar la solucin del problema (4.1), (4.2) cuando : 0 < < , 0 < < , y
() = sin3 .

Solucin: Recordamos de la trigonometra que


3 1
sin3 = sin sin 3
4 4
De modo que desde (4.9) vemos que 1 = 3/4, 3 = 1/4 y = 0 para cada otro valor de
. Por lo tanto, por (4.12), la solucin es:
3 sinh 1 sinh 3
(, ) = sin sin 3
4 sinh 4 sinh 3
El problema general de Neumann para la ecuacin de Laplace en el rectngulo R,

+ = 0

(, 0) = () (, ) = () (0 < < )
(4.14)
(0, ) = () (, ) = () (0 < < )

Puede tambin ser solucin de la misma manera, siempre que los datos lmites satisfagan la
condicin necesaria (2.12), lo que significa que en el presente caso

(4.15) [() ()] + [() ()] = 0
0 0

Ilustramos la situacin en este caso por el significado del siguiente problema.

Ejemplo 4.2 Encontrar una solucin al problema de Neumann


= 0 (0 < < , 0 < < )
(, 0) = /2 (0 )
(, ) = 0 (0 )
(0, ) = (, ) = 0 (0 )

Solucin: Por lo tanto notamos que la condicin (4.15) cumple con el dato de lmite () =
/2. El mtodo de separacin de las variables (, ) = ()() conduce al problema de
Sturm-Loiuville

+ = 0 (0) = 0 () = 0
Para la funcin X. Esto tiene los valores propios = 2 y la correspondiente funcin propia
() = cos , = 0, 1, 2, . Para cada valor propio = 2, el problema de valor inicial
para la funcin Y:

+ 2 = 0 () = 0
Dada () = cosh ( ). Adems, para solucionar el problema, asumimos una solucin en
serio en la forma:

(, ) = 0 + cosh ( ) cos

Consiste de la forma , 1. Diferenciando esta frmula con respecto ha y aplicando la


condicin limite = 0, encontramos:


= sinh cos
2
=1

= cos
=1

Donde = sinh , = 1, 2, . En vista del factor que 0 ( /2) = 0, vimos
que el termino constante 0 en la expansin en serie de coseno de Fourier de la funcin () =
/2 es cero. Por lo tanto, la serie dada anteriormente es precisamente la serie de coseno
de Fourier de la funcin, y as:
2
= ( ) cos
0 2
2 (1) 1
= ( = 1,2, )
2
Por lo tanto,
2 (1) 1
=
3 sinh
O
4
21 = ( = 1,2, )
(2 1)3 sinh(2 1)
Entonces una solucin del problema dado de Neumann es:

4 cosh(2 1)( )
(, ) = 0 + cos
(2 1)3 sinh(2 1)
=1

Y 0 es una constante arbitrario.

5. Problema de Dirichlet en un Disco


Sea D y C denotado respectivamente el interior en el contorno del disco 2 + 2 (Fig.
7.3) con radio y centro en el origen. Buscamos una funcin armnica en D, que asume el
valor prescrito = en el contorno C; esto es:

(5.1) + = 0 ,
=

Fig. 7-3 Problema de Dirichlet en una regin circular

En orden de ser capaz de aplicar el mtodo de separacin de variables a este problema,


introducimos las coordenadas polares

= cos , = sin (0 < < 2, 0 < < )


Para transformar la ecuacin de Laplace en la forma:
1 1
(5.2) + + 2 = 0

La condicin de contorno en (5.1) se vuelve:

(5.3) (, ) = () (0 2)

As, el problema (5.1) es reducido al problema equivalente (5.2), (5.3), donde ahora el mtodo
de separacin de variables puede ser aplicado. Asumimos la solucin particular de (5.2) en la
forma:

(, ) = () ()

Sustituyendo esto en (5.2) y separando variables, obtenemos


() () ()
2 + = =
() () ()
Donde es nuestra separacin constante. Esto gobierna a las dos ecuaciones diferenciales
ordinarias

(5.4) 2 () + () () = 0
Y

(5.5) () + () = 0
Para las funciones y , respectivamente.

En orden para la solucin del problema (5.2), (5.3) es necesario requerir que la solucin sea
peridica en , peridica en 2; esto es,

(, + 2) = (, )

Este requerimiento es claramente conocido si es peridico, peridica en 2. Por lo tanto,


impusimos la condicin peridica de contorno:

(5.6) () () = 0, () () = 0
Para la funcin. La ecuacin (5.5) junto con la condicin (5.6) constituye un problema de valor
propio [Problema 15, Ejercicio 4.4] para cada valor propio son = 2 , = 1, 2, , y la
correspondiente funcin propia son:

() = C cos + D sin
Cuando = 0, ecuacin (5.4) tiene la solucin general

R 0 () = A0 + B0 ln

Y cuando = 2 , = 1, 2, esto tiene la solucin general

R () = A + B
Adems, todas las funciones
(5.7) u (, ) = A0 + B0 ln + (A + B )(C cos + D sin )

= 0, 1, 2, , > 0 que son peridicas del perodo 2, satisfacen la ecuacin de Laplace (5.2)
en coordenadas polares. Estas son llamadas las funciones armnicas circulares.

Formaremos una combinacin infinita adecuada de las funciones (5.7) para nuestra solucin
del problema (5.2), (5.3). Puesto que la solucin debe ser continua en D, el cual contiene el
origen = 0, debemos excluir de nuestra consideracin el trmino logartmico y los trminos
que implican las potencias negativas de . Por lo tanto, elegimos B = 0, = 0, 1, 2, en (5.7)
y consideremos la serie infinita:

0
(5.8) u(r, ) = + ( cos + sin )
2
=1

Donde hemos cambiado de smbolo las constantes. Para determinar las constantes 0 , y ,
1, ponemos = y aplicamos la condicin de contorno (5.3) para obtener:

0
(5.9) () = + ( cos + sin )
2
=1
Esta es la expansin de la serie de Fourier de en el intervalo [, ], para la cual los
coeficientes estn dados por:
1
= = () cos ( = 0, 1, 2, )

(5.10)
1

= = () sin ( = 1, 2, )

Sustituyendo por y en (5.8), obtenemos



0
(5.11) u(r, ) = + ( ) ( cos + sin )
2
=1

Donde y son los coeficientes de Fourier en [, ], dados por las integrales en (5.10).

La frmula (5.11) junto con (5.10) da la solucin de nuestro problema (5.2), (5.3). Para
verificar esto, asumimos ahora que es continua, por partes lisa, y peridica del perodo 2 en
[, ]. Sea
1
= |()|

Entonces | | y | | , y define:

(r, ) = ( ) ( cos + sin ) ( = 1, 2, )

Con 0 = 0 /2. Luego, por 0 con cualquier 0 < ,

(5.12) | (r, )| ( ) [| | + | |] 2 ( )


La serie que consiste en los trminos 2 ( ) converge uniformemente para 0. Por el

teorema 6.1 del Captulo 5, la serie (5.11) tambin converge uniformemente para = . Por lo
tanto, la serie (5.11) converge uniformemente a (, ) para , y por lo tanto (, ) es
continua para 0 , 0 2. Dado que:

0
(, ) = + ( cos + sin ) = ()
2
=1

se cumple la condicin de contorno (5.3).

Para demostrar que u satisface la ecuacin de Laplace (5.2), observamos que la serie (5.11) y
sus derivadas parciales primera y segunda estn dominadas por algn mltiplo constante de la
serie 22 (/)2 . Esta serie converge uniformemente para 0 < , por lo que tiene
derivadas continuas de primer y segundo orden para < , que puede obtenerse diferenciando
a lo largo de la serie (5.11). Resulta que:

2 1 1 2 1
2
+ + 2 2 = 2 (, )[( 1) + 2 ] = 0

=1

As se verifica que satisface (5,2).


Ejemplo 5.1. Solucionar el problema de Dirichlet
1 1
+ + 2 = 0 ( < )

(, ) = cos 2 (0 2)

Solucin: Notamos que


1 1
cos2 = ( + cos 2)
2 2
Adems, por (5.11), tenemos

1 2
u(r, ) = ( + cos 2)
2

Ejemplo 5.2 Encontrar la solucin del problema de Dirichlet

+ = 0 2 + 2 < 1

= 2 2 + 2 = 1
Solucin: Considerando el problema en coordenadas polares, notamos que las condiciones de
1
contorno se vuelve (1, ) = sin2 = (1 cos 2). As, por (5.1), tenemos
2

1
u(r, ) = (1 r 2 cos 2)
2
Volviendo a las coordenadas rectangulares, encontramos
1
(, ) = (1 2 + 2 )
2
La cual es la solucin del problema.

6. Frmula Integral de Poisson

La solucin (5.11) del problema de Dirichlet en un disco puede expresarse en forma de una
integral a partir de la cual podemos deducir varias propiedades importantes y tiles de las
funciones armnicas. Esto se realiza sustituyendo las frmulas (5.10) por los coeficientes de
Fourier y en (5.11). Obtenemos:

1 1
(, ) = () + ( ) ()(cos cos + sin sin )
2
=1

1 1
(6.1) = () + ( ) () cos ( )
2 =1

1 1
= () [ + ( ) cos ( )]
2
=1

Donde hemos intercambiado el orden de sumatoria e integracin. Este intercambio est


justificado, ya que por (5.12) la serie

1
(6.2) + ( ) cos ( )
2
=1

converge uniformemente para < y .

Calculemos la suma de las series (6.2). Hacemos = = (cos + sin ), donde || =


< 1, y consideremos la serie

1
+
2
=1

Ya que || < 1, tenemos:



1 1 1+
+ = + =
2 2 1 2(1 )
=1

1 + cos + sin
=
2(1 cos sin )
(1 + cos + sin )(1 cos sin )
(6.3) =
2((1 cos )2 + ( sin )2 )
1 2 + 2 sin
=
2(1 2 cos + 2 )
Por el otro lado, ya que

= = (cos + sin )
Tenemos

1 1
(6.4) + = + (cos + sin )
2 2
=1 =1

As, al igualar las partes reales de los lados derechos de (6.3) y (6.4), encontramos

1 1 2 + 2 sin
(6.5) + cos =
2 2(1 2 cos + 2 )
=1

Se deduce que (6.2) tiene la suma



1 2 2
+ ( ) cos ( ) =
2 2(2 2 cos( ) + 2 )
=1

De modo que (6.1) se convierte

1 2 2
(6.6) (, ) = 2 ()
2 2 cos( ) + 2

Esto se conoce como frmula integral de Poisson para la solucin del problema de Dirichlet en
un disco.

Recordamos que en (5.11), se asumi continuo y por partes lisas adems de ser peridico
de periodo 2. En la frmula integral (6.6), basta con exigir que sea continua y peridica. Para
entonces la funcin definida por (6.6) para < y por (, ) = () para = es continua
en todo el disco . Adems, tiene derivadas continuas de todos los rdenes con respecto
a y para < , que pueden obtenerse por diferenciacin dentro del signo integral. Por un
clculo directo, se ve fcilmente que la funcin

1 2 2
(6.7) (, , ) =
2 2 2 cos( ) + 2
satisface (5.2) para < ,0 2. As, la frmula integral (6.6) con (, ) = ()
representa la solucin nica del problema de Dirichlet (5.2), (5.3). Con frecuencia, esta frmula
es ms prctica de usar que la serie (5.11).

La funcin (6.7) a menudo se llama ncleo de Poisson. Esta funcin es positiva para 0
, . De hecho, puesto que 2 2 > 0, vemos que

2 2 cos( ) + 2 = ( )2 + 2[1 cos( )]


( )
= ( )2 + 4 sin2
2
>0
Adems, el ncleo de Poisson satisface la propiedad interesante:

1 2 2
(6.8) 1= 2
2 2 cos( ) + 2

Este resultado se deduce de (6.6) cuando () = 1 y el hecho de que en tal caso (, ) = 1


por el principio mximo.

Si ponemos = 0 en (6.6), entonces obtenemos el resultado muy importante


1
(0,0) = ()
2

O

1
(6.9) (0,0) = (, )
2

Ya que (, ) = () y = . Esta frmula expresa la bien conocida propiedad de valor


medio de las funciones armnicas. Digamos esta propiedad como un teorema en un contexto
ms general.

Teorema 6.1. Sea una funcin armnica en un dominio D. Entonces el valor de en el centro
de cualquier disco que est en D es igual a la media o media de los valores de en el lmite del
disco.

Cabe sealar que la propiedad de valor medio de una funcin armnica puede establecerse
utilizando simplemente la armona de la funcin sin recurrir a la frmula integral de Poisson.
Para ver esto, sea D, un disco de radio y centro (, ) que yace en D. Entonces cada punto
(, ) en D, puede expresarse en coordenadas polares

= + cos (0 )

= + sin (0 2)
Con el polo en (, ). Aplicando la identidad (2.11) al dominio , obtenemos.
2

(6.10) ( + cos , + sin ) =
0

Puesto que es armnico, = 0 y as (6.10) se convierte


2

(6.11) ( + cos , + sin ) = 0
0

Configuremos
2
(6.12) (, ; ) = ( + cos , + sin )
0

Puesto que y sus derivadas son continuas en , as es junto con sus primeras derivadas. Se
deduce que podemos calcular / diferenciando (6.12) bajo el signo integral. En vista de
(6.11), encontramos
2

= ( + cos , + sin ) = 0
0

Esto implica que es independiente de ; en particular, tenemos


(6.13) (, ; 0) = (, ; ) (0 )

Ahora, desde (6.12), vemos que


2
(, , 0) = lim ( + cos , + sin )
0 0

2
= (, )
0

= 2(, )
Por lo tanto, por (6.12) y (6.13), finalmente obtenemos

1 2
(6.14) (, ) = ( + cos , + sin )
2 0

Lo que demuestra el teorema 6.1.

Esto produce el caso especial (6.9) cuando el punto (, ) pasa a ser el origen.

Es interesante notar que la inversa del teorema 6.1 es tambin verdadera. Es decir, si es
continua y satisface la propiedad de valor medio (6.4) para cualquier en un dominio D,
entonces es armnica en D. De hecho, bajo el requisito de continuidad se puede demostrar
que , como se expres en la frmula (6.14), tiene derivadas continuas de todos los rdenes en
D que tambin satisfacen la propiedad del valor significativo. Esto implica que una funcin
armnica necesariamente tiene derivados de todos los rdenes. Aqu, slo mostraremos que
= 0 en D.
Supongamos que hay un punto (, ) en D en el que 0 y > 0. Puesto que es
continua, hay un disco D, alrededor de (, ) con radio y situado en D tal que > 0 en todo
. Entonces. Aplicando la frmula (2.11) a , encontramos

0 <

2
=
0 0
2
u
=
0 r

2
= ( + cos , + sin )
0

Que, por (6.14), produce:


0 < = [(, )] = 0

Una contradiccin. Por lo tanto, = 0 en todo D.

Debido al teorema 6.1 estamos ahora en posicin de establecer la declaracin hecha en la


Seccin 3 en el sentido de que una constante es la nica funcin armnica que puede alcanzar
un mximo dentro de un dominio. Declararemos y demostraremos este hecho como un
teorema.

Teorema 6.2. Sea una funcin armnica en un dominio D. Si alcanza su valor mximo en un
punto dentro de D, entonces es una constante.

Fig. 7-4 El principio mximo fuerte


Demostracin: Sea el valor mximo de en D (Fig. 7.4), que se alcanza en un punto
0 : (0 , 0 ) dentro de . Se demostrar que = idnticamente en . Sea cualquier otro
punto dentro de y unir 0 y mediante una curva suave situada en . A partir del punto 0
y movindose hacia el punto , construiremos un nmero finito de discos situados en y con
centros en tales que el centro de cada disco se encuentra dentro del disco anterior.
Probaremos que = en cada uno de los discos. As, sea 0 un disco con centro en 0 y un
radio 0 tal que 0 est en . Entonces, por la propiedad de valor medio de , tenemos:

1 2
(6.15) = (0 ) = (0 , )
2 0

Ahora, en la circunferencia 0 de 0, tenemos (0 , ) . Si (0 , ) < para un cierto


punto en 0 , entonces por la continuidad tendramos < sobre un cierto arco de 0 que
contiene ese punto; por consiguiente>

1 2
(0 , ) ) <
2 0

Que contradice (6.15). Por lo tanto, = idnticamente en 0 . Pero el argumento precedente


se sostiene tambin para cualquier disco concntrico con radio para 0 0.

Por lo tanto, = idnticamente en 0.

Ahora, sea 1 el punto de interseccin de 0 con la curva de modo que (1 ) = . Sea 1


un disco sobre 1 con radio 1 tal que 1 est en . Repitiendo el argumento precedente,
concluimos que = en todo 1. Procediendo de esta manera, se sigue que despus de un
nmero finito de pasos llegaremos a un disco que contiene el punto y para el cual =
idnticamente en , as () = y dado que el punto es arbitrario en , se deduce que
= en todo .
DEDICATORIA

La dedicatoria de este trabajo grupal est dedicado a nuestros padres,

pilares fundamentales en nuestras vidas. Sin ellos, jams hubisemos

podido conseguir lo que hasta ahora estamos logrando como hijos y

estudiantes universitarios. Su tenacidad y lucha insaciable han hecho de

ellos el gran ejemplo a seguir y destacar, no solo para nosotros como

estudiantes universitarios, sino para nuestros hermanos y familia en

general. Tambin dedicamos este proyecto a nuestro maestro Jos Simen

Quique Broncano, sin su gua y conocimientos sobre el curso, no seramos

capaces de entender el maravilloso mundo que nos abre impartiendo el

curso de Ecuaciones Diferenciales Parciales.


PRLOGO

Los estudiantes de la Escuela Acadmico Profesional de Ingeniera


Mecnica de Fluidos pertenecientes a la Universidad Nacional Mayor de San
Marcos (UNMSM) saluda a los lectores del presente trabajo denominado
La Ecuacin de Laplace, el cual tiene por finalidad contribuir a la
educacin y la formacin integral de los estudiantes en el aspecto
intelectual, moral y fsico. Solo as el hombre estar en la capacidad de
explicarse objetivamente la realidad, asumir una actitud crtica ante los
problemas, plantear soluciones, expresar sensibilidad ante el dolor
humano, profesar y aplicar principios nobles.

El presente texto, La Ecuacin de Laplace, es producto del trabajo


organizado de profesores de matemtica y estudiantes de la Escuela
Acadmico Profesional de Ingeniera Mecnica de Fluidos (Universidad
Nacional Mayor de San Marcos); est dirigido a estudiantes que inician su
formacin universitaria en las especialidades de ingenieras y ciencias de las
diferentes universidades del pas.

Saludamos el aporte de nuestros profesores, el esfuerzo y dedicacin de los


estudiantes de la Escuela Acadmico Profesional de Ingeniera Mecnica de
Fluidos y el apoyo de los padres de familia quienes, a pesar de los problemas
sociales y econmicos, buscan una educacin cientfica y humanista que
beneficie a nuestra sociedad.
UNIVERSIDAD NACIONAL MAYOR DE SAN MARCOS
Universidad del Per Decana de Amrica

FACULTAD DE CIENCIAS FSICAS


E.A.P. INGENIERA MECNICA DE FLUIDOS

LA ECUACIN DE LAPLACE EUTIQUIO C. YOUNG

CURSO : Ecuaciones Diferenciales Parciales

DOCENTE : Mg. QUIQUE BRONCANO, Jos Simen

INTEGRANTES : NGELES DURAND, Kennedy (14130125)


BELITO BARRENECHEA, Jedy John ( )
GUTIRREZ ENCISO, Miguel ngel (14130211)
REYNOSO MARTNEZ, Kevin A. (14130054)
SAMANIEGO ESPINOZA, Eyner A. ( )

LIMA PER

2017

También podría gustarte