Laplace

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 31

Capı́tulo 7

Ecuación de Laplace

Problema 7.1 Sea u(x, y) una función armónica para x, y ∈ (0, π) que es nula
en el borde de dicho cuadrado salvo u(x, 0) = 3 sin x, x ∈ (0, π).

Solución:

Sabemos que la ecuación se puede resolver usando series de Fourier,



X sinh n(π − y)
u(x, y) = fn sin nx,
n=1
sinh nπ

si la función u(x, 0) admite desarrollo de Fourier en el intervalo (0, π),


∞ π
2
X Z
u(x, 0) = fn sin nx, fn = f (x) sin nx dx.
n=1
π 0

El desarrollo de f (x) = 3 sin x en serie de senos no precisa integrales, ya que


podemos leer directamente que todos los coeficientes fn son nulos salvo f1 = 3.
Por tanto, la solución del problema es

3 sinh(π − y)
u(x, y) = sin x. ✷
sinh π
Problema 7.2 Hallar en forma de serie la solución del problema para la ecua-
ción de Laplace ∆u(x, y) = uxx + uyy = 0, u(0, y) = 0, u(L, y) = 0, u(x, 0) =
f (x), uy (x, 0) = g(x) para x ∈ (0, L), y > 0. Aplicarlo al caso L = π, u(x, 0) =
0, uy (x, 0) = N −1 sin N x, donde N es un número natural. Tomar el lı́mite de la
solución cuando N tiende a infinito. Interpretar el resultado. Aplicarlo también
al caso L = π, u(x, 0) = 0, uy (x, 0) = x.

Solución:

Se trata de un problema mixto, ya que las dos primeras condiciones son de


contorno y las dos últimas, iniciales, pensando en y como en una coordenada
temporal.
Intentaremos resolver la ecuación de Laplace por separación de variables,
buscando soluciones de la forma u(x, y) = X(x)Y (y), que sustituidas en la

163
ecuación,

X ′′ (x) Y ′′ (y)
0 = uxx + uyy = X ′′ (x)Y (y) + X(x)Y ′′ (y) ⇒ − = = λ,
X(x) Y (y)

nos permiten separar las variables de la ecuación y reducirla a dos ecuaciones


ordinarias,
X ′′ (x) + λX(x) = 0, Y ′′ (y) − λY (y) = 0.
Las condiciones de contorno se traducen en condiciones de contorno para la
primera ecuación ordinaria,

X(0) = 0, X(L) = 0.

Resolvemos primero el problema para X(x),

X ′′ (x) + λX(0) = 0, X(0) = 0, X(L) = 0,

que corresponde a condiciones de contorno de extremos fijos y ya lo hemos


resuelto con anterioridad y tiene por autovalores y autofunciones
 nπ 2 nπx
λn = , n = 1, 2, . . . , Xn (x) = sin , n = 1, 2, . . .
L L
y sus múltiplos.
Por tanto, por superposición lineal de soluciones u(x, y) = X(x)Y (y), ob-
tendremos soluciones más generales

X nπx
u(x, y) = Yn (y) sin .
n=1
L

Abordamos la otra ecuación, lineal ordinaria y homogénea, con valores de λ


restringidos ya a los autovalores del problema de contorno,
 nπ 2
Yn′′ − Yn = 0,
L
cuya solución general, dado que tiene autovalores ±nπ/L es
nπy nπy
Yn (y) = an enπy/L + bn e−nπy/L = An cosh + Bn sinh ,
L L
para n = 1, 2, . . .. Usaremos la forma hiperbólica de la solución general. Por
tanto,
∞ 
X nπy nπy  nπx
u(x, y) = An cosh + Bn sinh sin ,
n=1
L L L

resuelve tanto la ecuación como las condiciones de contorno.


Falta utilizar las condiciones iniciales, u(x, 0) = f (x), uy (x, 0) = g(x). Co-
menzamos con la primera,

X nπx
f (x) = u(x, 0) = An sin .
n=1
L

164
Por tanto, si f (x) tiene desarrollo en serie de senos en el intervalo [0, L],
∞ L
nπx 2 nπx
X Z
f (x) = fn sin , fn = f (x) sin dx,
n=1
L L 0 L

comparando ambas expresiones, como el desarrollo de Fourier es único, obtene-


mos los coeficientes An = fn , n = 1, 2, . . .
Los coeficientes Bn los obtenemos con la otra condición, para lo cual nece-
sitamos la derivada,

X nπ  nπy nπy  nπx
uy (x, y) = An sinh + Bn cosh sin ,
n=1
L L L L


X nπ nπx
g(x) = uy (x, 0) = Bn sin .
n=1
L L

Ası́ pues, si g(x) tiene desarrollo en serie de senos en el intervalo [0, L],
∞ L
nπx 2 nπx
X Z
g(x) = gn sin , gn = g(x) sin dx,
n=1
L L 0 L

comparando ambas expresiones, como el desarrollo de Fourier es único, obtene-


mos los coeficientes Bn = Lgn /nπ, n = 1, 2, . . .
Ya tenemos identificados todos los coeficientes de nuestra solución, ası́ que
podemos afirmar:
“El problema de la ecuación de Laplace uxx + uyy = 0 para x ∈ (0, L), y > 0
con condiciones mixtas u(x, 0) = f (x), uy (x, 0) = g(x), u(0, y) = 0 = u(L, y)
para x ∈ (0, L), y > 0, tiene solución
∞  
X nπy Lgn nπy nπx
u(x, y) = fn cosh + sinh sin ,
n=1
L nπ L L

si las funciones f (x), g(x) admiten desarrollo convergente de Fourier en el in-


tervalo (0, L)”
∞ L
nπx 2 nπx
X Z
f (x) = fn sin , fn = f (x) sin dx,
n=1
L L 0 L

∞ L
nπx 2 nπx
X Z
g(x) = gn sin , gn = g(x) sin dx. ✷
n=1
L L 0 L

En el caso, L = π, las expresiones se simplifican


∞ 
X gn 
u(x, y) = fn cosh ny + sinh ny sin nx,
n=1
n


X ∞
X
f (x) = fn sin nx, g(x) = gn sin nx.
n=1 n=1

165
Si f (x) = 0, es obvio que todos los coeficientes fn son nulos. Si g(x) =
N −1 sin N x, la propia función es su desarrollo en serie, ası́ que gN = 1/N y el
resto de coeficientes son nulos. Por tanto, la solución en este caso es
sinh N y sin N x
u(x, y) = .✷
N2
Observamos que, si tomamos el lı́mite N tendiendo a infinito, el dato inicial
es trivial u(x, 0) = 0, uy (x, 0) = 0, ya que el seno está acotado por la unidad,
ası́ que g(x) decrece como 1/N . La solución en ese caso es la trivial u(x, y) ≡
0, como es fácil comprobar. Sin embargo, para y 6= 0, la solución u(x, y) =
N −2 sinh N y sin N x tiende a infinito como una exponencial eN y , en lugar de
tender a cero.
Por tanto, las condiciones iniciales conducen a malos problemas para la ecua-
ción de Laplace. Esto es razonable, ya que la ecuación de Laplace aparece en
problemas estacionarios y, por ello, no tiene mucho sentido considerar la coor-
denada y como un tiempo. ✷
En el caso f (x) = 0 todos los coeficientes fn son nulos. Y g(x) = x tiene
coeficientes

2 π 2(−1)n+1

2 sin nx x cos nx
Z
gn = x sin nx dx = − = ,
π 0 π n2 n 0 n
con lo cual la serie de la solución es

X 2(−1)n+1
u(x, y) = sinh ny sin nx. ✷
n=1
n2

La serie de Fourier en senos en [0, π] de g(x) = x, al ser de clase C ∞ , converge


a x para x ∈ (0, π) y converge a cero en los extremos. Por tanto, la serie sólo
converge a un valor distinto del de la función en x = π. ✷

Problema 7.3 Hallar en forma de serie la solución del problema de Dirichlet


para la ecuación de Laplace ∆u(x, y) = 0 para x ∈ (0, L), y ∈ (0, M ) con
condiciones de contorno u(x, M ) = f (x), u(x, 0) = 0 para x ∈ (0, L), u(0, y) =
0 = u(L, y) para y ∈ (0, M ). Aplicar el resultado al caso en el que L = π = M ,
f (x) = 5 sin x cos x. Comprobar explı́citamente que la función obtenida satisface
la ecuación y las condiciones de contorno. Aplicar el resultado al caso en el que
L = π = M , f (x) = 1.

Solución:

Consideremos funciones armónicas u(x, y) es un rectángulo, x ∈ (0, L), y ∈


(0, M ). Intentaremos resolver la ecuación de Laplace por separación de variables,
buscando soluciones de la forma u(x, y) = X(x)Y (y), que sustituidas en la
ecuación,
X ′′ (x) Y ′′ (y)
0 = uxx + uyy = X ′′ (x)Y (y) + X(x)Y ′′ (y) ⇒ − = = λ,
X(x) Y (y)
nos permiten separar las variables de la ecuación y reducirla a dos ecuaciones
ordinarias,
X ′′ (x) + λX(x) = 0, Y ′′ (y) − λY (y) = 0.

166
Tenemos un problema de Dirichlet simplificado, en el que sólo hay condicio-
nes de contorno no triviales en uno de los lados del rectángulo,
u(x, M ) = f (x), u(x, 0) = 0, u(0, y) = 0 = u(L, y),
para x ∈ (0, L), y ∈ (0, M ).
Las condiciones de contorno se traducen en condiciones de contorno para las
ecuaciones ordinarias,
X(0) = 0 = X(L), Y (0) = 0.
Resolvemos primero el problema para X(x),
X ′′ (x) + λX(x) = 0, X(0) = 0 = X(L).
Este problema ya lo hemos resuelto con anterioridad y tiene por autovalores
y autofunciones,
 nπ 2 nπx
λn = , Xn (x) = sin , n∈N
L L
Abordamos la otra ecuación,
 nπ 2 nπy nπy
Yn′′ − Y = 0 ⇒ Yn (y) = An sinh + Bn cosh ,
L L L
para n ∈ N.
Imponemos la condición de contorno Yn (0) = 0,
0 = Yn (0) = Bn , n ∈ N,
y, por tanto, a falta de la condición inhomogénea, la solución se puede expresar
en forma de serie,

X nπy nπx
u(x, y) = An sinh sin .
n=0
L L

Falta por imponer la condición inhomogénea u(x, M ) = f (x),



X nπM nπx
f (x) = u(x, M ) = An sinh sin ,
n=1
L L

que es un desarrollo en serie de senos.


Por tanto, si la función f (x) admite desarrollo de Fourier en serie de senos,

nπx 2 L nπx nπM
X Z
f (x) = fn sin , fn = f (x) sin dx = An sinh ,
n=0
L L 0 L L

tenemos la solución del problema de contorno:


“El problema de la ecuación de Laplace ∆u(x, y) = 0 para x ∈ (0, L), y ∈
(0, M ) con condiciones de contorno u(x, M ) = f (x), u(x, 0) = 0 para x ∈ (0, L),
u(0, y) = 0 = u(L, y) para y ∈ (0, M ) tiene solución
nπy
∞ sinh
L sin nπx ,
X
u(x, y) = fn
nπM L
n=1 sinh
L

167
si la función f (x) admite desarrollo de Fourier en el intervalo (0, L) en serie de
senos”

nπx 2 L nπx
X Z
f (x) = fn sin , fn = f (x) sin dx, n ∈ N.
n=1
L L 0 L

En el caso L = π = M se simplifica el resultado anterior,


∞ ∞
X sinh ny X
u(x, y) = fn sin nx, f (x) = fn sin nx,
n=1
sinh nπ n=1

π
2
Z
fn = f (x) sin nx dx, n ∈ N.
π 0

La función f (x) = 5 sin x cos x tiene un desarrollo de Fourier con sólo un


término,
5 5
5 sin x cos x = sin 2x ⇒ f2 = , fn = 0, n 6= 2,
2 2
5 sinh 2y
u(x, y) = sin 2x,
2 sinh 2π
como se comprueba explı́citamente,
5
u(x, 0) = 0 = u(0, y) = u(π, y), u(x, π) = sin 2x,
2
sinh 2y
uyy (x, y) = 10 sin 2x = −uxx(x, y). ✷
sinh 2π
Para la función f (x) = 1, en cambio, el desarrollo de Fourier tiene infinitos
términos,

2 π

2 h cos nx iπ 2
Z
0 n par
fn = sin nx dx = − = (1 − (−1)n ) =
π 0 π n 0 nπ 4/nπ n impar

4 X sinh ny
u(x, y) = sin nx ✷.
π n impar
n sinh nπ

Problema 7.4 Hallar en forma de serie la solución del problema de Neumann


para la ecuación de Laplace ∆u(x, y) = 0 para x ∈ (0, L), y ∈ (0, M ) con
condiciones de contorno uy (x, 0) = −g(x), uy (x, M ) = 0 para x ∈ (0, L),
ux (0, y) = 0 = ux (L, y) para y ∈ (0, M ).

Solución:

Consideremos funciones armónicas u(x, y) es un rectángulo, x ∈ (0, L), y ∈


(0, M ). Intentaremos resolver la ecuación de Laplace por separación de variables,
buscando soluciones de la forma u(x, y) = X(x)Y (y), que sustituidas en la
ecuación,

X ′′ (x) Y ′′ (y)
0 = uxx + uyy = X ′′ (x)Y (y) + X(x)Y ′′ (y) ⇒ − = = λ,
X(x) Y (y)

168
uy
M
y
-uX uX

x L
-uy

Figura 7.1: Problema de Neumann en el rectángulo [0, L] × [0, M ]

nos permiten separar las variables de la ecuación y reducirla a dos ecuaciones


ordinarias,
X ′′ (x) + λX(x) = 0, Y ′′ (y) − λY (y) = 0.

Tenemos un problema de Neumann simplificado, en el que sólo hay condi-


ciones de contorno no triviales en uno de los lados del rectángulo,

uy (x, 0) = −g(x), uy (x, M ) = 0, ux (0, y) = 0 = ux (L, y),

para x ∈ (0, L), y ∈ (0, M ).


Las condiciones de contorno se traducen en condiciones de contorno para las
ecuaciones ordinarias,

X ′ (0) = 0 = X ′ (L), Y ′ (M ) = 0.

Resolvemos primero el problema para X(x),

X ′′ (x) + λX(x) = 0, X ′ (0) = 0 = X ′ (L).

Este problema ya lo hemos resuelto con anterioridad y tiene por autovalores


y autofunciones,
 nπ 2 nπx
λn = , Xn (x) = cos , n = 0, 1, 2, . . .
L L
Abordamos la otra ecuación,
 nπ 2
Yn′′ − Y = 0 ⇒ Y0 (y) = A0 + B0 y, Yn (y) = An enπy/L + Bn e−nπy/L ,
L
para n ∈ N.
Imponemos la condición de contorno Yn′ (M ) = 0,

0 = Y0′ (M ) = B0 ⇒ Y0 (y) = A0 ,
nπ nπM/L nπ
0 = Yn′ (M ) = An e − Bn e−nπM/L ⇒ Bn = An e2nπM/L , n ∈ N.
L L

169
 
Yn (y) = An enπy/L + e−nπy/L e2nπM/L
  nπ(M − y)
= An enπM/L enπ(M−y)/L + e−nπ(M−y)/L = an cosh
L
denotando a0 = A0 , an = 2An enπM/L , n ∈ N, y, por tanto, a falta de la
condición inhomogénea, la solución se puede expresar en forma de serie,

X nπ(M − y) nπx
u(x, y) = an cosh cos .
n=0
L L

Falta por imponer la condición inhomogénea uy (x, 0) = −g(x),



X nπ nπM nπx
g(x) = −uy (x, 0) = an sinh cos ,
n=1
L L L

que es un desarrollo en serie de cosenos al que le falta el coeficiente g0 .


Por tanto, si la función g(x) admite desarrollo de Fourier en serie de cosenos,

g0 X nπx 2 L nπx nπan nπM
Z
g(x) = + gn cos , gn = g(x) cos dx = sinh ,
2 n=1
L L 0 L L L

2 L
Z
g0 = g(x) dx = 0,
L 0
tenemos la solución del problema de contorno:
“El problema de la ecuación de Laplace ∆u(x, y) = 0 para x ∈ (0, L), y ∈
(0, M ) con condiciones de contorno uy (x, 0) = −g(x), uy (x, M ) = 0 para x ∈
(0, L), ux (0, y) = 0 = ux (L, y) para y ∈ (0, M ) tiene solución
∞ nπ(M−y)
L X gn cosh L nπx
u(x, y) = K + nπM
cos ,
π n=1 n sinh L L

si la función g(x) admite desarrollo de Fourier en el intervalo (0, L) en serie de


cosenos,

nπx 2 L nπx
X Z
g(x) = gn cos , gn = g(x) cos dx,
n=1
L L 0 L

con coeficiente g0 nulo,”


L
2
Z
g0 = g(x) dx = 0.
L 0

La constante K = a0 aparece ya que el problema de Neumann tiene solución


única salvo constante. Es la constante a0 que el desarrollo de Fourier no puede
identificar, ya que desaparece al calcular la derivada.
La condición de tener el coeficiente g0 era esperable, ya que el dato del
problema de Neumann tiene que tener integral nula en el borde del recinto.
En nuestro caso, como el dato se anula en tres de los lados del rectángulo, la
condición necesaria se reduce a la anulación de la integral de la función g(x) en
el cuarto lado. ✷

170
Problema 7.5 Hallar en forma de serie la solución del problema de Neu-
mann para la ecuación de Laplace ∆u(x, y) = 0 para x ∈ (0, L), y ∈ (0, M )
con condiciones de contorno uy (x, M ) = g(x), uy (x, 0) = 0 para x ∈ (0, L),
ux (0, y) = 0 = ux (L, y) para y ∈ (0, M ). Aplicar el resultado al caso en el que
L = π = M , g(x) = 2 sin2 3x. Lo mismo para g(x) = 2 sin2 3x − 1.

Solución:

Se puede abordar con un cambio de variable usando el resultado anterior,


pero lo resolveremos directamente.
Consideremos funciones armónicas u(x, y) es un rectángulo, x ∈ (0, L), y ∈
(0, M ). Intentaremos resolver la ecuación de Laplace por separación de variables,
buscando soluciones de la forma u(x, y) = X(x)Y (y), que sustituidas en la
ecuación,

X ′′ (x) Y ′′ (y)
0 = uxx + uyy = X ′′ (x)Y (y) + X(x)Y ′′ (y) ⇒ − = = λ,
X(x) Y (y)

nos permiten separar las variables de la ecuación y reducirla a dos ecuaciones


ordinarias,
X ′′ (x) + λX(x) = 0, Y ′′ (y) − λY (y) = 0.
Tenemos un problema de Neumann simplificado, en el que sólo hay condi-
ciones de contorno no triviales en uno de los lados del rectángulo,

uy (x, M ) = g(x), uy (x, 0) = 0, ux (0, y) = 0 = ux (L, y),

para x ∈ (0, L), y ∈ (0, M ).


Las condiciones de contorno se traducen en condiciones de contorno para las
ecuaciones ordinarias,

X ′ (0) = 0 = X ′ (L), Y ′ (0) = 0.

Resolvemos primero el problema para X(x),

X ′′ (x) + λX(x) = 0, X ′ (0) = 0 = X ′ (L).

Este problema ya lo hemos resuelto con anterioridad y tiene por autovalores


y autofunciones,
 nπ 2 nπx
kn = , Xn (x) = cos , n = 0, 1, 2, . . .
L L
Abordamos la otra ecuación,
 nπ 2 nπy nπy
Yn′′ − Y = 0 ⇒ Y0 (y) = A0 +B0 y, Yn (y) = An cosh +Bn sinh ,
L L L
para n ∈ N.
Imponemos la condición de contorno Yn′ (0) = 0,

0 = Y0′ (0) = B0 ⇒ Y0 (y) = A0 ,



0 = Yn′ (0) = Bn
L

171
y observamos que todos los coeficientes Bn tienen que ser nulos.
Por tanto, a falta de la condición inhomogénea, la solución se puede expresar
en forma de serie,

X nπy nπx
u(x, y) = An cosh cos .
n=0
L L

Falta por imponer la condición inhomogénea uy (x, M ) = g(x),



X nπ nπM nπx
g(x) = uy (x, M ) = An sinh cos ,
n=1
L L L

que es un desarrollo en serie de cosenos al que le falta el coeficiente g0 .


Por tanto, si la función g(x) admite desarrollo de Fourier en serie de cosenos,

nπx 2 L nπx nπAn nπM
X Z
g(x) = gn cos , gn = g(x) cos dx = sinh ,
n=0
L L 0 L L L

L
2
Z
g0 = g(x) dx = 0,
L 0
tenemos la solución del problema de contorno:
“El problema de la ecuación de Laplace ∆u(x, y) = 0 para x ∈ (0, L), y ∈
(0, M ) con condiciones de contorno uy (x, M ) = g(x), uy (x, 0) = 0 para x ∈
(0, L), ux (0, y) = 0 = ux (L, y) para y ∈ (0, M ) tiene solución

L X gn cosh nπ
L y nπx
u(x, y) = K + nπM
cos ,
π n=1 n sinh L L

si la función g(x) admite desarrollo de Fourier en el intervalo (0, L) en serie de


cosenos,

nπx 2 L nπx
X Z
g(x) = gn cos , gn = g(x) cos dx,
n=1
L L 0 L

con coeficiente g0 nulo,”


L
2
Z
g0 = g(x) dx = 0.
L 0

La constante K = a0 aparece ya que el problema de Neumann tiene solución


única salvo constante. Es la constante a0 que el desarrollo de Fourier no puede
identificar, ya que desaparece al calcular la derivada.
La condición de tener el coeficiente g0 era esperable, ya que el dato del
problema de Neumann tiene que tener integral nula en el borde del recinto.
En nuestro caso, como el dato se anula en tres de los lados del rectángulo, la
condición necesaria se reduce a la anulación de la integral de la función g(x) en
el cuarto lado. ✷
Usando trigonometrı́a elemental,

g0 X
g(x) = 2 sin2 3x = 1 − cos 6x = + gn cos nx,
2 n=1

172
la función g(x) tiene desarrollo de Fourier con coeficientes nulos salvo g0 = 2,
g6 = −1. Por tanto, este dato de contorno no cumple la condición necesaria
g0 = 0, por lo que el problema de Neumann en este caso no tiene solución. ✷
En cambio, en el otro caso,

g0 X
g(x) = 2 sin2 3x − 1 = − cos 6x = + gn cos nx,
2 n=1

la función g(x) tiene desarrollo de Fourier con coeficientes nulos salvo g6 = −1.
Por tanto, la solución de problema de Neumann en este caso es
cosh 6y cos 6x
u(x, y) = K − .✷
6 sinh 6π
Problema 7.6 Sea u(x, y) una función armónica para x, y ∈ (0, π) cuya de-
rivada normal es nula en el borde de dicho cuadrado salvo uy (x, 0) = 5 cos 3x,
x ∈ (0, π).

Solución:

En primer lugar, comprobamos la condición necesaria para que el problema


de Neumann tenga solución,
Z π Z π Z π
du
Z
ds = − uy (x, 0) dx + ux (π, y) dy − uy (x, π)
∂Ω dν
Z π0 0
π
0
5
Z
π
+ ux (0, y) dy = −5 cos 3x dx = − [sin 3x]0 = 0.
0 0 3
Luego el problema de Neumann tiene solución, que podemos obtener usando
series de Fourier,

X gn cosh n(π − y)
u(x, y) = K + cos nx,
n=1
n sinh nπ

si la función g(x) = −uy (x, 0) admite desarrollo de Fourier en el intervalo (0, π)


en serie de cosenos,

X nπ
g(x) = gn cos x,
n=1
L
con coeficiente g0 nulo.
El desarrollo de g(x) = −5 cos 3x en serie de cosenos no precisa integrales,
ya que podemos leer directamente que todos los coeficientes gn son nulos salvo
g3 = −5.
Por tanto, la solución del problema es
5 cosh 3(π − y)
u(x, y) = K − cos 3x. ✷
3 sinh 3π
Problema 7.7 Hallar en forma de serie la solución del problema para la ecua-
ción de Laplace ∆u(x, y) = 0 para x ∈ (0, L), y ∈ (0, M ) con condiciones de con-
torno u(x, 0) = 0, u(x, M ) = f (x) para x ∈ (0, L), ux (0, y) = 0 = ux (L, y) para
y ∈ (0, M ). Aplicar el resultado al caso en el que L = π = M , f (x) = cos2 2x.
Lo mismo para f (x) = cos(x/2).

173
Solución:

Consideremos funciones armónicas u(x, y) es un rectángulo, x ∈ (0, L), y ∈


(0, M ). Intentaremos resolver la ecuación de Laplace por separación de variables,
buscando soluciones de la forma u(x, y) = X(x)Y (y), que sustituidas en la
ecuación,

X ′′ (x) Y ′′ (y)
0 = uxx + uyy = X ′′ (x)Y (y) + X(x)Y ′′ (y) ⇒ − = = λ,
X(x) Y (y)

nos permiten separar las variables de la ecuación y reducirla a dos ecuaciones


ordinarias,
X ′′ (x) + λX(x) = 0, Y ′′ (y) − λY (y) = 0.
Tenemos un problema de contorno simplificado, en el que sólo hay condicio-
nes de contorno no triviales en uno de los lados del rectángulo,

u(x, M ) = f (x), u(x, 0) = 0, ux (0, y) = 0 = ux (L, y),

para x ∈ (0, L), y ∈ (0, M ).


Las condiciones de contorno se traducen en condiciones de contorno para las
ecuaciones ordinarias,

X ′ (0) = 0 = X ′ (L), Y (0) = 0.

Resolvemos primero el problema para X(x),

X ′′ (x) + λX(x) = 0, X ′ (0) = 0 = X ′ (L).

Este problema ya lo hemos resuelto con anterioridad y tiene por autovalores


y autofunciones,
 nπ 2 nπx
kn = , Xn (x) = cos , n = 0, 1, 2, . . .
L L
Abordamos la otra ecuación,
 nπ 2 nπy nπy
Yn′′ − Yn = 0 ⇒ Y0 (y) = A0 +B0 y, Yn (y) = An cosh +Bn sinh ,
L L L
para n ∈ N.
Imponemos la condición de contorno Yn (0) = 0,

0 = Y0 (0) = A0 ⇒ Y0 (y) = B0 y,

0 = Yn (0) = An
y observamos que todos los coeficientes An tienen que ser nulos.
Por tanto, a falta de la condición inhomogénea, la solución se puede expresar
en forma de serie,

X nπ nπx
u(x, y) = B0 y + Bn sinh y cos .
n=1
L L

174
Falta por imponer la condición inhomogénea u(x, M ) = f (x),

X nπM nπx
f (x) = u(x, M ) = B0 M + Bn sinh cos ,
n=1
L L

que es un desarrollo en serie de cosenos.


Por tanto, si la función f (x) admite desarrollo de Fourier en serie de cosenos,
∞ L
f0 X nπx 2 nπx nπM
Z
f (x) = + fn cos , f (x) cos dx = fn = Bn sinh ,
2 n=1
L L 0 L L

L
2
Z
f (x) dx = f0 = 2B0 M,
L 0

tenemos la solución del problema de contorno:


“El problema de la ecuación de Laplace ∆u(x, y) = 0 para x ∈ (0, L), y ∈
(0, M ) con condiciones de contorno u(x, M ) = f (x), u(x, 0) = 0 para x ∈ (0, L),
ux (0, y) = 0 = ux (L, y) para y ∈ (0, M ) tiene solución

f0 X sinh nπy
L nπx
u(x, y) = y+ fn nπM
cos ,
2M n=1
sinh L L

si la función f (x) admite desarrollo de Fourier en el intervalo (0, L) en serie de


cosenos,
∞ L
f0 X nπx 2 nπx
Z
f (x) = + fn cos , fn = f (x) cos dx.
2 n=1
L L 0 L

Este problema no es ni de Dirichlet, ni de Neumann, ya que mezcla ambas


condiciones, según los lados del rectángulo. La solución en forma de series es
única y existe si f (x) admite desarrollo en serie de cosenos. ✷
En el caso L = π = M ,

f0 X sinh ny
u(x, y) = y+ fn cos nx,
2π n=1
sinh nπ

∞ π
f0 X 2
Z
f (x) = + fn cos nx, fn = f (x) cos nx dx.
2 n=1
π 0

Usando trigonometrı́a elemental,



1 + cos 4x f0 X
f (x) = cos2 2x = = + fn cos nx,
2 2 n=1

la función f (x) tiene desarrollo de Fourier con coeficientes nulos salvo f0 = 1,


f4 = 1/2. Por tanto, la solución de problema en este caso es

y sinh 4y cos 4x
u(x, y) = + .✷
2π 2 sinh 4π

175
En el otro caso, f (x) = cos(x/2), sı́ que nos vemos abocados a calcular el
desarrollo de Fourier,

2 π 1 π
Z      
1 1
Z
fn = cos(x/2) cos nx dx = cos n + x + cos n − x dx
π 0 π 0 2 2
x=π
(−1)n (−1)n
  
1 sin (n + 1/2) x sin (n − 1/2) x 1
= + = −
π n + 1/2 n − 1/2 x=0 π n + 1/2 n − 1/2
1 (−1)n+1
= , n = 0, 1, . . .
π n2 − 1/4

con lo cual la solución del problema es



2 X 1 (−1)n+1 sinh ny
u(x, y) = y + cos nx. ✷
π2 n=1
π n2 − 1/4 sinh nπ

Problema 7.8 Una función u(r, φ) satisface la ecuación de Laplace en un


cı́rculo de radio R centrado en el origen, y verifica u(R, φ) = cos2 φ. Hallar
u(r, φ) en el interior del cı́rculo.

Solución:

La función u(r, φ) es periódica en el ángulo φ, u(r, 0) = u(r, 2π), por lo


que este problema se puede abordar por series de Fourier, que verifican esta
propiedad. Sabemos que las funciones armónicas en el disco de radio a se pueden
expresar como

a0 X  r  n
u(r, φ) = + (an cos nφ + bn sin nφ) ,
2 n=1
R

siendo an , bn los coeficientes del desarrollo de Fourier del valor de la función en


el borde,

a0 X 1 cos 2φ
u(R, φ) = + (an cos nφ + bn sin nφ) = cos2 φ = + ,
2 n=1
2 2

que en nuestro caso no es preciso realizar, ya que es patente que todos los
coeficientes bn son nulos, lo mismo que los an , salvo a0 = 1, a2 = 1/2. Por
tanto, la solución del problema es

1 r2
u(r, φ) = + cos 2φ. ✷
2 2R2
Problema 7.9 La función u(r, φ) es armónica en el interior del cı́rculo de
radio a centrado en el origen, siendo u(a, φ) = 5 + sin φ cos φ. Hallar u(r, ϕ) en
el interior del cı́rculo.

Solución:

La función u(r, φ) es periódica en el ángulo φ, u(r, 0) = u(r, 2π), por lo


que este problema se puede abordar por series de Fourier, que verifican esta

176
propiedad. Sabemos que las funciones armónicas en el disco de radio a se pueden
expresar como

a0 X  r  n
u(r, φ) = + (an cos nφ + bn sin nφ) ,
2 n=1
a

siendo an , bn los coeficientes del desarrollo de Fourier del valor de la función en


el borde,

a0 X sin 2φ
u(a, φ) = + (an cos nφ + bn sin nφ) = 5 + ,
2 n=1
2

que en nuestro caso no es preciso realizar, ya que es patente que todos los coe-
ficientes son nulos, salvo a0 = 10, b2 = 1/2. Por tanto, la solución del problema
es
r2
u(r, φ) = 5 + 2 sin 2φ. ✷
2a
Problema 7.10 Hallar en forma de serie la solución del problema de Dirichlet
para la ecuación de Laplace ∆u(r, φ) = 0 en el semicı́rculo de radio R, y ≥ 0,
con condiciones de contorno u(R, φ) = f (φ), u(r, 0) = 0 = u(r, π), r ∈ (0, R),
φ ∈ (0, π). Aplicarlo al caso f (φ) = 6 sin 4x cos 4x.

Solución:

Consideremos funciones armónicas u(r, φ) en un cı́rculo centrado en el origen


y de radio R. Intentaremos resolver la ecuación de Laplace en coordenadas
polares, r ∈ (0, R), φ ∈ (0, π),
ur uφφ
urr + + 2 = 0,
r r
por separación de variables, buscando soluciones de la forma u(r, φ) = R(r)Φ(φ),
que sustituidas en la ecuación,
ur uφφ R′ (r)Φ(φ) R(r)Φ′′ (φ)
0 = urr + + 2 = R′′ (r)Φ(φ) + + ,
r r r r2
r2 R′′ (r) + rR′ (r) Φ′′ (φ)
=− = λ,
R(r) Φ(φ)
nos permiten separar las variables de la ecuación y reducirla a dos ecuaciones
ordinarias,

r2 R′′ (r) + rR′ (r) − λR(r) = 0, Φ′′ (φ) + λΦ(φ) = 0.

Las condiciones de contorno nos dan el valor de u en r = R y sobre el eje Y ,


donde toma el valor nulo. Las condiciones homogéneas u(r, 0) = 0, u(r, π) = 0
conducen a condiciones de contorno

Φ(0) = 0, Φ(π) = 0,

lo cual lleva a un problema de contorno para Φ ya resuelto con anterioridad,


que tiene por autovalores λn = n2 y autofunciones Φn (φ) = sin nφ, n ∈ N.

177
Queda por resolver la ecuación en la coordenada radial,

r2 Rn′′ (r) + rRn′ (r) − n2 Rn (r) = 0,

que es una ecuación de Euler, que se reduce, como hemos visto, a una ecuación
lineal con coeficientes constantes con el cambio de variable independiente r =
es ⇔ s = ln r,

R̈n − n2 Rn = 0 ⇒ Rn (s) = An ens + Bn e−ns = An rn + Bn r−n , n = 1, 2, . . .

denotando por un punto la derivada con respecto a s.


Llegado este punto, observamos que aparecen funciones singulares, las po-
tencias inversas r−n . Imponiendo que las soluciones sean funciones acotadas y
regulares en el cı́rculo, es patente que todos los coeficientes Bn tienen que ser
nulos. Por tanto, las funciones radiales son Rn (r) = An rn , n = 1, . . .
Y podemos escribir la solución de la ecuación de Laplace en el cı́rculo como

X
u(r, φ) = An rn sin nφ
n=1

Usando la condición de contorno inhomogénea, u(R, φ) = f (φ),



X
u(R, φ) = An Rn sin nφ = f (φ),
n=1

y si la función f (φ) admite desarrollo de Fourier en serie de senos,


∞ π
2
X Z
f (φ) = fn sin nφ, fn = f (φ) sin nφ dφ,
n=1
π 0

observamos, por la unicidad del desarrollo en serie de Fourier, que An Rn = fn ,


n ∈ N y, por tanto, podemos identificar los coeficientes de la función armónica
con los coeficientes del desarrollo de Fourier:
“El problema de la ecuación de Laplace ∆u(r, φ) = 0 para r ∈ (0, R), φ ∈
(0, π) con condiciones de contorno u(R, φ) = f (φ), u(r, 0) = 0 = u(r, π) tiene
solución
X∞  r n
u(r, φ) = fn sin nφ (7.1)
n=1
R

si la función f (φ) admite desarrollo de Fourier en serie de senos,”


∞ π
2
X Z
f (φ) = fn sin nφ fn = f (φ) sin nφ dφ. ✷
n=1
π 0

En nuestro caso particular f (φ) = 3 sin 8φ, con lo cual el desarrollo de Fourier
en senos sólo tiene no nulo el término f8 = 3 y, por tanto, la solución del
problema de Dirichlet queda

3r8
u(r, φ) = sin 8φ. ✷
R8

178
Problema 7.11 Sea u una función armónica en el disco de radio a centrado
en el origen y que se anula en el diámetro y = 0. Demostrar que u(x, −y) =
−u(x, y).

Solución:

Sabemos que las funciones armónicas en el disco de radio a se pueden ex-


presar como

a0 X  r  n
u(r, φ) = + (an cos nφ + bn sin nφ) .
2 n=1
a

El diámetro y = 0 es el eje X corresponde a las semirrectas φ = 0, φ = π. Por


tanto, la función verifica u(r, 0) = 0 = u(r, π), lo que trasladado al desarrollo
supone

a0 X  r  n
0 = u(r, 0) = + an ,
2 n=1
a

serie nula de términos positivos, de donde se deduce que los coeficientes an son
todos nulos. Con la condición u(r, π) se llega a la misma conclusión.
Por tanto, la condición de anulación impone que el desarrollo sea sólo de
senos,

a0 X  r  n
u(r, φ) = + bn sin nφ.
2 n=1
a

Comprobamos la condición de imparidad u(x, −y) = −u(x, y). Como en


coordenadas polares cambiar y por −y es una reflexión respecto al eje X, supone
cambiar φ por −φ,
∞ ∞
a0 X  r n a0 X  r  n
u(r, −φ) = + bn sin n(−φ) = − + bn sin nφ = −u(r, φ),
2 n=1 a 2 n=1 a

ya que el seno es una función impar. ✷

Problema 7.12 Sea u(x, y) una función armónica para x ∈ (0, π), y > 0. Se
tiene que u(x, 0) = 1 para x ∈ (0, π), u(0, y) = 0 = u(π, y) para y > 0. Hallar
u(x, y).

Solución:

Dado que u(0, y) = 0 = u(π, y), hacemos un desarrollo de Fourier en serie


de senos que cumpla dichas condiciones de contorno,

X
u(x, y) = Yn (y) sin nx,
n=1

que sustituimos en la ecuación de Laplace,



X
Y ′′ (y) − n2 Y (y) sin nx,

0 = uxx + uyy =
n=1

179
y por la unicidad del desarrollo de Fourier, obtenemos un sistema infinito de
ecuaciones ordinarias lineales homogéneas,
Y ′′ − n2 Y = 0 ⇒ Yn (y) = an cosh ny + bn sinh ny, n ∈ N.
Por tanto, la función más general que es solución de la ecuación de Laplace
con dichas condiciones de contorno es

X
u(x, y) = (an cosh ny + bn sinh ny) sin nx.
n=1

Sabemos además que u(x, 0) = 1. Por tanto,



X
1 = u(x, 0) = an sin nx,
n=1

y los coeficientes an no son más que los del desarrollo de la función unidad en
serie de senos,
2 π 2 1 − (−1)n

2 h cos nx iπ
Z
0 n par
an = sin nx dx = − = = .
π 0 π n 0 π n 4/nπ n impar
Pero esto fija sólo los coeficientes an , no los bn . Luego el problema está mal
planteado, ya que es indeterminado, como ya sabı́amos, ya que es un proble-
ma mixto de valores iniciales para la ecuación de Laplace, que sólo permite
problemas de contorno.
El problema tendrı́a solución única si impusiéramos la condición adicional
de que u sea acotada en infinito. Dado que
an + bn ny an − bn −ny
an cosh ny + bn sinh ny = e + e ,
2 2
eliminar los términos divergentes en infinito supone eliminar las exponenciales
crecientes eny , lo que se consigue si bn = −an .
Por tanto, la función buscada tiene desarrollo en serie

X 4 X e−ny
u(x, y) = an e−ny sin nx = sin nx
n=1
π n impar
n

4 X e−(2n−1)y
= sin(2n − 1)x. ✷
π n=1 2n − 1

Problema 7.13 Hallar en forma de serie la solución del problema de Neumann


para la ecuación de Laplace ∆u(r, φ) = 0 en el cı́rculo de radio R con condiciones
de contorno ur (R, φ) = g(φ), r ∈ (0, R), φ ∈ (0, 2π). Aplicar el resultado al caso
en el que g(φ) = 1 y al caso en el que g(φ) = sin φ cos φ.
Solución:

Consideremos funciones armónicas u(r, φ) en un cı́rculo centrado en el origen


y de radio R. Intentaremos resolver la ecuación de Laplace en coordenadas
polares, r ∈ (0, R), φ ∈ (0, 2π),
∂ 2 u 1 ∂u 1 ∂2u
+ + = 0,
∂r2 r ∂r r2 ∂φ2

180
por separación de variables, buscando soluciones de la forma u(r, φ) = R(r)Φ(φ),
que sustituidas en la ecuación,

ur uφφ R′ (r)Φ(φ) R(r)Φ′′ (φ)


0 = urr + + 2 = R′′ (r)Φ(φ) + + ,
r r r r2

r2 R′′ (r) + rR′ (r) Φ′′ (φ)


=− = λ,
R(r) Φ(φ)
nos permiten separar las variables de la ecuación y reducirla a dos ecuaciones
ordinarias,

r2 R′′ (r) + rR′ (r) − λR(r) = 0, Φ′′ (φ) + λΦ(φ) = 0.

Independientemente del problema de contorno que tengamos que resolver,


las propias coordenadas polares exigen condiciones de periodicidad,

Φ(0) = Φ(2π), Φ′ (0) = Φ′ (2π),

lo cual conduce a un problema de contorno para Φ ya resuelto con anterioridad,


que tiene por autovalores λn = n2 y autofunciones Φ0 (φ) = 1, Φn (φ) = cos nφ,
Φ̃n (φ) = sin nφ, n ∈ N.
Queda por resolver la ecuación en la coordenada radial,

r2 Rn′′ (r) + rRn′ (r) − n2 Rn (r) = 0,

que es una ecuación de Euler, que se reduce, como hemos visto, a una ecuación
lineal con coeficientes constantes con el cambio de variable independiente r =
es ⇔ s = ln r,

R̈n − n2 Rn = 0 ⇒ Rn (s) = An ens + Bn e−ns = An rn + Bn r−n , n 6= 0,

R̈0 = 0 ⇒ R0 (s) = A0 + B0 s = A0 + B0 ln r,
denotando por un punto la derivada con respecto a s.
Llegado este punto, observamos que aparecen funciones singulares, las poten-
cias inversas r−n y ln r. Imponiendo que las soluciones sean funciones acotadas
y regulares en el cı́rculo, es patente que todos los coeficientes Bn tienen que ser
nulos. Por tanto, las funciones radiales son Rn (r) = rn , n = 0, 1, . . .
Y podemos escribir la solución de la ecuación de Laplace en el cı́rculo como

X  
u(r, φ) = A0 + rn An cos nφ + Ãn sin nφ .
n=1

El problema de Neumann supone conocer la derivada de u normal al borde


del recinto. En este caso el recinto es un cı́rculo, ası́ que la dirección normal
es la radial, por lo que la condición de contorno de Neumann es de la forma
ur (R, φ) = g(φ). Sustituyendo,

X  
ur (R, φ) = nRn−1 An cos nφ + Ãn sin nφ = g(φ),
n=1

181
y si la función g(φ) admite desarrollo de Fourier en serie de senos y cosenos,

a0 X
g(φ) = + (an cos nφ + bn sin nφ) ,
2 n=1

2π 2π
1 1
Z Z
an = g(φ) cos nφ dφ, bn = g(φ) sin nφ dφ,
π 0 π 0

observamos, por la unicidad del desarrollo en serie de Fourier, que A0 = K está


indeterminado, ya que tenemos sólo condiciones sobre las derivadas de u, algo
que sucede siempre con los problemas de Neumann.
El coeficiente a0 del desarrollo del dato de Neumann tiene que ser nulo, lo
cual era de esperar, ya que los problemas de Neumann tienen como condición
necesaria para tener solución que la integral de la derivada normal de u a lo
largo del borde del recinto se anule y, en nuestro caso,
Z Z 2π
0= ur (R, φ) ds = R g(φ) dφ = Rπa0 .
Γ 0

Por su parte, nRn−1 An = an , nRn−1 Bn = bn , n ∈ N y, por tanto, pode-


mos identificar los coeficientes de la función armónica con los coeficientes del
desarrollo de Fourier:
“El problema de la ecuación de Laplace ∆u(r, φ) = 0 para r ∈ (0, R), φ ∈
(0, 2π) con condición de contorno ur (R, φ) = g(φ) tiene solución

X rn
u(r, φ) = K + (an cos nφ + bn sin nφ) ,
n=1
nRn−1

unica salvo una constante K, si la función g(φ) admite desarrollo de Fourier en


serie de senos y cosenos con término independiente a0 nulo”

X
g(φ) = (an cos nφ + bn sin nφ) ,
n=1

2π 2π
1 1
Z Z
an = g(φ) cos nφ dφ, bn = g(φ) sin nφ dφ.
π 0 π 0

Esta solución en serie de Fourier se puede sumar, al igual que se hizo para
el problema de Dirichlet,

X r0n
u(r0 , φ0 ) − K = (an cos nφ0 + bn sin nφ0 )
n=1
nRn−1

1 2π r0n
Z X
= g(φ) (cos nφ0 cos φ + sin nφ0 sin φ) dφ
π 0 n=1
nRn−1

1 2π r0n
Z X
= g(φ) cos n(φ − φ0 ) dφ,
π 0 n=1
nRn−1

ya que las series de Fourier se pueden integrar bajo condiciones muy generales.

182
El integrando es la integral de una serie geométrica, como se comprueba
desarrollando el coseno en exponenciales imaginarias,
∞ ∞ ∞
X r0n X r0n inγ
X rn
0
2 cos nγ = e + e−inγ
n=1
nRn−1 n=1
nR n−1
n=1
nR n−1

iγ n−1
Z r0 ∞  ∞  n−1 !
reiγ


X re −iγ
X
= e +e dr
0 n=1
R n=1
R
Z r0 
Reiγ Re−iγ

= + dr
0 R − reiγ R − re−iγ
r0
−R ln R − reiγ + ln R − re−iγ 0
 
=
2R ln R − R ln R2 + r02 − 2Rr0 cos γ ,

=

del cual podemos omitir el término constante, ya que lo podemos incorporar a


la constante K,

R
Z
g(φ) ln R2 + r02 − 2Rr0 cos(φ − φ0 ) dφ. ✷

u(r0 , φ0 ) = K −
2π 0

Procedemos con los casos particulares: el caso g(φ) = 1 no tiene solución, ya


que no se cumple la condición necesaria,
Z 2π
g(φ) dφ = 2π 6= 0. ✷
0

El caso g(φ) = sin φ cos φ es sencillo, ya que su desarrollo de Fourier,

sin 2φ
g(φ) = ,
2
tiene todos los coeficientes nulos, salvo b2 = 1/2. Por tanto, la solución del
problema es
r2
u(r, φ) = K + sin 2φ. ✷
4R
Problema 7.14 Hallar en forma de serie la solución del problema de Neu-
mann para la ecuación de Laplace ∆u(x, y) = 0 para x ∈ (0, L), y ∈ (0, M )
con condiciones de contorno uy (x, M ) = g(x), uy (x, 0) = 0 para x ∈ (0, L),
ux (0, y) = 0 = ux (L, y) para y ∈ (0, M ). Aplicar el resultado al caso en el que
L = π = M , g(x) = 2 sin2 3x. Lo mismo para g(x) = 2 sin2 3x − 1.

Solución:

Consideremos funciones armónicas u(x, y) es un rectángulo, x ∈ (0, L), y ∈


(0, M ). Intentaremos resolver la ecuación de Laplace por separación de variables,
buscando soluciones de la forma u(x, y) = X(x)Y (y), que sustituidas en la
ecuación,

X ′′ (x) Y ′′ (y)
0 = uxx + uyy = X ′′ (x)Y (y) + X(x)Y ′′ (y) ⇒ − = = λ,
X(x) Y (y)

183
nos permiten separar las variables de la ecuación y reducirla a dos ecuaciones
ordinarias,
X ′′ (x) + λX(x) = 0, Y ′′ (y) − λY (y) = 0.
Tenemos un problema de Neumann simplificado, en el que sólo hay condi-
ciones de contorno no triviales en uno de los lados del rectángulo,

uy (x, M ) = g(x), uy (x, 0) = 0, ux (0, y) = 0 = ux (L, y),

para x ∈ (0, L), y ∈ (0, M ).


Las condiciones de contorno se traducen en condiciones de contorno para las
ecuaciones ordinarias,

X ′ (0) = 0 = X ′ (L), Y ′ (0) = 0.

Resolvemos primero el problema para X(x),

X ′′ (x) + λX(x) = 0, X ′ (0) = 0 = X ′ (L).

Este problema ya lo hemos resuelto con anterioridad y tiene por autovalores


y autofunciones,
 nπ 2 nπx
kn = , Xn (x) = cos , n = 0, 1, 2, . . .
L L
Abordamos la otra ecuación,
 nπ 2 nπy nπy
Yn′′ − Yn = 0 ⇒ Y0 (y) = A0 +B0 y, Yn (y) = An cosh +Bn sinh ,
L L L
para n ∈ N.
Imponemos la condición de contorno Yn′ (0) = 0,

0 = Y0′ (0) = B0 ⇒ Y0 (y) = A0 ,



0 = Yn′ (0) = Bn
L
y observamos que todos los coeficientes Bn tienen que ser nulos.
Por tanto, a falta de la condición inhomogénea, la solución se puede expresar
en forma de serie,

X nπy nπx
u(x, y) = An cosh cos .
n=0
L L

Falta por imponer la condición inhomogénea uy (x, M ) = g(x),



X nπ nπM nπx
g(x) = uy (x, M ) = An sinh cos ,
n=1
L L L

que es un desarrollo en serie de cosenos al que le falta el coeficiente g0 .


Por tanto, si la función g(x) admite desarrollo de Fourier en serie de cosenos,
∞ L
g0 X nπx 2 nπx nπAn nπM
Z
g(x) = + gn cos , gn = g(x) cos dx = sinh ,
2 n=1
L L 0 L L L

184
L
2
Z
g0 = g(x) dx = 0,
L 0
tenemos la solución del problema de contorno:
“El problema de la ecuación de Laplace ∆u(x, y) = 0 para x ∈ (0, L), y ∈
(0, M ) con condiciones de contorno uy (x, M ) = g(x), uy (x, 0) = 0 para x ∈
(0, L), ux (0, y) = 0 = ux (L, y) para y ∈ (0, M ) tiene solución
∞ nπy
L X gn cosh L nπx
u(x, y) = K + nπM
cos ,
π n=1 n sinh L L

si la función g(x) admite desarrollo de Fourier en el intervalo (0, L) en serie de


cosenos,

nπx 2 L nπx
X Z
g(x) = gn cos , gn = g(x) cos dx,
n=1
L L 0 L

con coeficiente g0 nulo,”


L
2
Z
g0 = g(x) dx = 0.
L 0

La constante K = a0 aparece ya que el problema de Neumann tiene solución


única salvo constante. Es la constante a0 que el desarrollo de Fourier no puede
identificar, ya que desaparece al calcular la derivada.
La condición de tener el coeficiente g0 era esperable, ya que el dato del
problema de Neumann tiene que tener integral nula en el borde del recinto.
En nuestro caso, como el dato se anula en tres de los lados del rectángulo, la
condición necesaria se reduce a la anulación de la integral de la función g(x) en
el cuarto lado. ✷
Usando trigonometrı́a elemental,

g0 X
g(x) = 2 sin2 3x = 1 − cos 6x = + gn cos nx,
2 n=1

la función g(x) tiene desarrollo de Fourier con coeficientes nulos salvo g0 = 2,


g6 = −1. Por tanto, este dato de contorno no cumple la condición necesaria
g0 = 0, por lo que el problema de Neumann en este caso no tiene solución. ✷
En cambio, en el otro caso,

g0 X
g(x) = 2 sin2 3x − 1 = − cos 6x = + gn cos nx,
2 n=1

la función g(x) tiene desarrollo de Fourier con coeficientes nulos salvo g6 = −1.
Por tanto, la solución de problema de Neumann en este caso es
cosh 6y cos 6x
u(x, y) = K − .✷
6 sinh 6π
Lo comprobamos,
cosh 6y sin 6x 6 cosh 6y cos 6x
ux (x, y) = , uxx (x, y) = ,
sinh 6π sinh 6π

185
sinh 6y cos 6x 6 cosh 6y cos 6x
uy (x, y) = − , uyy (x, y) = − .
sinh 6π sinh 6π
Obviamente, ux (0, y) = 0 = ux (π, y), porque el seno se anula en cero y π.
Igualmente uy (x, 0) = 0, porque se anula el seno hiperbólico. A su vez,

uy (x, π) = − cos 6x = g(x).

Y la ecuación de Laplace se verifica, ya que


6 cosh 6y cos 6x
uxx (x, y) = = −uyy (x, y). ✷
sinh 6π
Problema 7.15 Hallar en forma de serie la solución del problema para la
ecuación de Laplace ∆u(z, φ) = uzz + uφφ = 0, u(0, φ) = 0, u(L, φ) = f (φ),
u(z, −π) = u(z, π), uφ (z, −π) = uφ (z, π) para z ∈ (0, L), φ ∈ (−π, π). Aplicar
el resultado al caso en el que L = π, f (φ) = 4 sin 2φ cos 2φ + cos2 2φ.

Solución:

Intentaremos resolver la ecuación de Laplace por separación de variables,


buscando soluciones de la forma u(z, φ) = Z(z)Φ(φ), que sustituidas en la ecua-
ción,

Z ′′ (z) Φ′′ (φ)


0 = uzz + uφφ = Z ′′ (z)Φ(φ) + Z(z)Φ′′ (φ) ⇒ =− = λ,
Z(z) Φ(φ)

nos permiten separar las variables de la ecuación y reducirla a dos ecuaciones


ordinarias,
Z ′′ (z) − λZ(z) = 0, Φ′′ (φ) + λΦ(φ) = 0.
Tenemos un problema de Dirichlet,

u(0, φ) = 0, u(L, φ) = f (φ), u(z, −π) = u(z, π), uφ (z, −π) = uφ (z, π),

para z ∈ (0, L), φ ∈ (−π, π).


Las condiciones de contorno se traducen en condiciones de contorno para las
ecuaciones ordinarias,

Z(0) = 0, Φ(−π) = Φ(π), Φ′ (−π) = Φ′ (π).

Resolvemos primero el problema para Φ(φ),

Φ′′ (φ) + λΦ(φ) = 0, Φ(−π) = Φ(π), Φ′ (−π) = Φ′ (π),

que corresponde a condiciones de contorno periódicas y ya lo hemos resuelto


con anterioridad y tiene por autovalores y autofunciones,

kn = n2 , n = 0, 1, . . . ,

Φ0 (φ) = 1, Φn (φ) = cos nφ, Φ̃n (φ) = sin nφ, n = 1, 2, . . .


Abordamos la otra ecuación,

Zn′′ − n2 Zn = 0 ⇒ Z0 (z) = A0 + B0 z, Zn (z) = An cosh nz + Bn sinh nz,

186
para n ∈ N.
Imponemos la condición de contorno Zn (0) = 0,

0 = Z0 (0) = A0 ⇒ Z0 (z) = B0 z,

0 = Zn (0) = An ⇒ Zn (z) = Bn sinh nz,


y observamos que todos los coeficientes An tienen que ser nulos.
Por tanto, a falta de la condición inhomogénea, la solución se puede expresar
en forma de serie,
∞ 
X 
u(z, φ) = B0 z + Bn cos nφ + B̃n sin nφ sinh nz.
n=1

Falta por imponer la condición inhomogénea u(L, φ) = f (φ),


∞ 
X 
f (φ) = u(L, φ) = B0 L + Bn cos nφ + B̃n sin nφ sinh nL,
n=1

que es un desarrollo en serie de senos y cosenos para f (φ).


Por tanto, si la función f (φ) admite desarrollo de Fourier en serie de senos
y cosenos,
∞ ∞
f0 X
f˜n sin nφ,
X
f (φ) = + fn cos nφ +
2 n=1 n=1

1 π 1 π
Z Z
fn = f (φ) cos nφ dφ, f˜n = f (φ) sin nφ dφ,
π −π π −π
comparando ambos desarrollos, concluimos

f0 fn f˜n
B0 = , Bn = , B̃n = , n = 1, 2, . . .
2L sinh nL sinh nL
“El problema de la ecuación de Laplace uzz + uφφ = 0 para z ∈ (0, L),
φ ∈ (−π, π) con condiciones de contorno u(0, φ) = 0, u(L, φ) = f (φ), u(z, −π) =
u(z, π), uφ (z, −π) = uφ (z, π) tiene solución
∞ 
f0  sinh nz
fn cos nφ + f˜n sin nφ
X
u(z, φ) = z+ ,
2L n=1
sinh nL

si la función f (φ) admite desarrollo de Fourier en el intervalo (−π, π) en serie


de senos y cosenos”
∞ ∞
f0 X
f˜n sin nφ,
X
f (φ) = + fn cos nφ +
2 n=1 n=1
π π
1 1
Z Z
fn = f (φ) cos nφ dφ, f˜n = f (φ) sin nφ dφ. ✷
π −π π −π
Teniendo en cuenta que un cilindro circular de altura L y radio R, parame-
trizado en coordenadas cilı́ndricas por g(z, φ) = (R cos φ, R sin φ, z), z ∈ (0, L),
φ ∈ (−π, π) tiene por laplaciano
uφφ
∆u = uzz + ,
R2

187
observamos que este problema es un problema de Dirichlet para un cilindro de
radio unidad y altura L con datos de contorno sobre las circunferencias situadas
en los planos z = 0 y z = L.
Usando trigonometrı́a elemental,

1 1
f (φ) = 4 sin 2φ cos 2φ + cos2 2φ = 2 sin 4φ + + cos 4φ,
2 2
observamos que la función f (φ) tiene desarrollo de Fourier con coeficientes nulos
salvo f0 = 1, f4 = 1/2, f˜4 = 2. Por tanto, la solución de problema de Dirichlet
en este caso es
 
1 1 1
u(z, φ) = z+ 2 sin 4φ + cos 4φ sinh 4z. ✷
2L sinh 4L 2

Lo comprobamos. Obviamente, u(0, φ) = 0, porque se anula el seno hiperbóli-


co, y u es una función periódica de periodo 2π por serlo el seno y el coseno.
Además,
1 1
u(L, φ) = + 2 sin 4φ + cos 4φ = f (φ),
2 2
y también se verifica la ecuación de Laplace,
 
16 1
uzz + uφφ = 2 sin 4φ + cos 4φ sinh 4z
sinh 4L 2
 
16 1
− 2 sin 4φ + cos 4φ sinh 4z = 0. ✷
sinh 4L 2

Problema 7.16 Hallar en forma de serie la solución del problema para la ecua-
ción de Laplace R2 uzz + uφφ = 0, uz (0, φ) = 0, uz (h, φ) = g(φ), u(z, −π) =
u(z, π), uφ (z, −π) = uφ (z, π) para z ∈ (0, h), φ ∈ (−π, π), siendo constantes
h, R > 0. Aplicarlo al caso h = 1 = R, g(φ) = φ. Aplicarlo al caso h = 1 = R,
g(φ) = φ2 .

Solución:

Se trata de un problema simplificado de Neumann para la ecuación de La-


place en un rectángulo con condiciones de contorno periódicas, es decir, en un
cilindro de radio R y altura h. Los datos de contorno se dan en las dos circun-
ferencias del borde, situadas en z = 0, z = h.
Intentaremos resolver la ecuación de Laplace por separación de variables,
buscando soluciones de la forma u(z, φ) = Z(z)Φ(φ), que sustituidas en la ecua-
ción,

R2 Z ′′ (z) Φ′′ (φ)


0 = R2 uzz + uφφ = R2 Z ′′ (z)Φ(φ) + Z(z)Φ′′ (φ) ⇒ =− = λ,
Z(z) Φ(φ)

nos permiten separar las variables de la ecuación y reducirla a dos ecuaciones


ordinarias,
λ
Z ′′ (z) − 2 Z(z) = 0, Φ′′ (φ) + λΦ(φ) = 0.
R

188
Las condiciones de contorno se traducen en condiciones de contorno para las
ecuaciones ordinarias,

Z ′ (0) = 0, Φ(−π) = Φ(π), Φ′ (−π) = Φ′ (π).

Resolvemos primero el problema para Φ(φ),

Φ′′ (φ) + λΦ(φ) = 0, Φ(−π) = Φ(π), Φ′ (−π) = Φ′ (π),

que corresponde a condiciones de contorno periódicas y ya lo hemos resuelto


con anterioridad y tiene por autovalores y autofunciones,

kn = n2 , n = 0, 1, . . . ,

Φ0 (φ) = 1, Φn (φ) = cos nφ, Φ̃n (φ) = sin nφ, n = 1, 2, . . .


Abordamos la otra ecuación,

n2 nz nz
Zn′′ − Zn = 0 ⇒ Z0 (z) = A0 + B0 z, Zn (z) = An cosh + Bn sinh ,
R2 R R
para n = 1, 2, . . .
Necesitamos las derivadas,
nAn nz nBn nz
Z0′ (z) = B0 , Zn′ (z) = sinh + cosh ,
R R R R
para imponer la condición de contorno Zn′ (0) = 0,

0 = Z0′ (0) = B0 ⇒ Z0 (z) = A0 ,

nBn nz
0 = Zn′ (0) = ⇒ Zn (z) = An cosh ,
R R
y observamos que todos los coeficientes Bn tienen que ser nulos.
Por tanto, a falta de la condición inhomogénea, haciendo uso del principio
de superposición lineal, la solución se puede expresar en forma de serie,
∞ 
X  nz
u(z, φ) = A0 + An cos nφ + Ãn sin nφ cosh ,
n=1
R

y su derivada, que necesitaremos,



X n  nz
uz (z, φ) = An cos nφ + Ãn sin nφ sinh .
n=1
R R

Falta por imponer la condición inhomogénea uz (h, φ) = g(φ),



X n  nh
g(φ) = uz (h, φ) = An cos nφ + Ãn sin nφ sinh ,
n=1
R R

que es un desarrollo en serie de senos y cosenos para g(φ), sin término constante,
lo que impone una condición sobre el dato de contorno, como sucede en los
problemas de Neumann.

189
Por tanto, si la función g(φ) admite desarrollo de Fourier en serie de senos
y cosenos,
∞ ∞
g0 X X
g(φ) = + gn cos nφ + g̃n sin nφ,
2 n=1 n=1
1 π 1 π
Z Z
gn = f (φ) cos nφ dφ, g̃n = g(φ) sin nφ dφ,
π −π π −π
comparando ambos desarrollos, concluimos
Rgn Rg̃n
g0 = 0, An = , Ãn = , n = 1, 2, . . .
n sinh nh
R n sinh nh
R

El coeficiente A0 no queda determinado, como corresponde a un problema


de Neumann, que tiene solución única salvo una constante K = A0 .
“El problema de la ecuación de Laplace R2 uzz + uφφ = 0 para z ∈ (0, h), φ ∈
(−π, π) con condiciones de contorno uz (0, φ) = 0, uz (h, φ) = g(φ), u(z, −π) =
u(z, π), uφ (z, −π) = uφ (z, π) tiene solución
∞ 
cosh nz

X gn g̃n R
u(z, φ) = K + R cos nφ + sin nφ nh
,
n=1
n n sinh R

si la función g(φ) admite desarrollo de Fourier en el intervalo (−π, π) en serie


de senos y cosenos
∞ ∞
g0 X X
g(φ) = + gn cos nφ + g̃n sin nφ,
2 n=1 n=1
π π
1 1
Z Z
gn = g(φ) cos nφ dφ, g̃n = g(φ) sin nφ dφ,
π −π π −π
con coeficiente g0 nulo.” ✷
En el caso h = 1 = R las expresiones se simplifican,
∞  
X gn g̃n cosh nz
u(z, φ) = K + cos nφ + sin nφ .
n=1
n n sinh n

El caso g(φ) = φ2 no tiene solución, ya que


1 π 2 2π 3
Z
g0 = φ dφ = 6= 0. ✷
π −π 3
El caso g(φ) = φ, al tratarse esta de una función impar, su desarrollo no
tiene términos en coseno. Por tanto, todos los coeficientes gn son nulos, incluido
g0 . Luego habrá solución.
Calculamos el resto de coeficientes,

1 π 2(−1)n+1

sin nφ φ cos nφ
Z
g̃n = φ sin nφ dφ = 2
− = ,
π −π n n −π n
y, por tanto, la solución de este problema es

X 2(−1)n+1 cosh nz
u(z, φ) = K + 2
sin nφ .✷
n=1
n sinh n

190
Problema 7.17 Sea u(x, y) una función tal que uxx + uyy = −2 sin x sin y para
x, y ∈ (0, π) que satisface u = 0 en el borde de dicho cuadrado. Calcular u(x, y).
Solución:

La función satisface las condiciones de contorno u(0, y) = 0 = u(π, y),


u(x, 0) = 0 = u(x, π). Podemos realizar un desarrollo en serie de senos de
la función,

X
u(x, y) = Yn (y) sin nx,
n=1
que ya satisface las condiciones u(0, y) = 0 = u(π, y).
El término inhomogéneo F (x, y) = −2 sin x sin y, ya está desarrollado y tiene
sólo el coeficiente n = 1.
Sustituyendo el desarrollo en la ecuación de Poisson,

X
Yn′′ (y) − n2 Yn (y) = sin nx,

−2 sin y sin x = uxx + uyy =
n=1

obtenemos, usando la unicidad del desarrollo de Fourier, el siguiente sistema de


problemas de contorno,

′′ 2 −2 sin y n = 1
Yn (y) − n Yn (y) = , Y (0) = 0 = Y (π),
0 n 6= 1
haciendo uso de las condiciones restantes, u(x, 0) = 0 = u(x, π).
La solución general de las ecuaciones homogéneas Yn′′ − n2 Yn = 0 es
Ynh (y) = An cosh ny + Bn sinh ny.
Para n > 1 no hay término inhomogéneo, por lo que las condiciones de
contorno Y (0) = 0 = Y (π) imponen
0 = An , 0 = An cosh nπ + Bn sinh nπ ⇒ An = 0 = Bn , n > 1,
con lo que estas ecuaciones sólo tienen solución trivial Yn (y) = 0.
Para n = 1 necesitamos una solución particular para la fuerza f (y) =
−2 sin y. Sabemos que la solución tiene que ser de la forma Y1p (y) = a cos y +
b sin y, que sustituida en la ecuación,
′′
0 = 2 sin y + Y1p − Y1p = (−2a) cos y + (2 − 2b) sin y ⇒ a = 0, 2 − 2b = 0,
nos proporciona los coeficientes indeterminados a = 0, b = 1, con lo cual la
solución general de la ecuación es
Y1 (y) = A1 cosh y + B1 sinh y + sin y,
y le podemos imponer las condiciones de contorno,
0 = Y1 (0) = A1 , 0 = Y1 (π) = A1 cosh π + B1 sinh π ⇒ A1 = 0 = B1 ,
con lo cual llegamos a la conclusión de que los coeficientes del desarrollo de
Fourier son Y1 (y) = sin y, Yn (y) = 0, n > 1, que sustituidos en el desarrollo nos
proporcionan la solución del problema,

X
u(x, y) = Yn (y) sin nx = sin x sin y. ✷
n=1

191
Problema 7.18 Sea u una función armónica y no negativa en el cı́rculo de
radio R centrado en el origen, y continua en su borde. Demostrar las desigual-
dades
R − ||x|| R + ||x||
u(0) ≤ u(x) ≤ u(0).
R + ||x|| R − ||x||
Utilizar la fórmula integral de Poisson.
Solución:

Sabemos que una función armónica en el cı́rculo de radio R se puede construir


como
R2 − r02 2π f (φ)
Z
u(r0 , φ0 ) = 2 + r2 − 2Rr cos(φ − φ )
dφ,
2π 0 R 0 0 0

siendo f (φ) = u(R, φ) su valor en la circunferencia del borde.


El denominador de esta integral se puede acotar y sacar fuera de la integral,
dado que el coseno toma sus valores en el intervalo [−1, 1],
(R − r0 )2 ≤ R2 + r02 − 2Rr0 cos(φ − φ0 ) ≤ (R + r0 )2 ,
1 R − r0 2π 1 R + r0 2π
Z Z
f (φ) dφ ≤ u(r0 , φ0 ) ≤ f (φ) dφ.
2π R + r0 0 2π R − r0 0
La demostración se completa con el cálculo del valor de u en el centro del
cı́rculo, Z 2π
1
u(0) = u(0, φ0 ) = f (φ) dφ,
2π 0
valor que se puede obtener también con la ley de la media.
Juntando ambos resultados, concluimos
R − r0 R + r0
u(0) ≤ u(r0 , φ0 ) ≤ u(0),
R + r0 R − r0
resultado conocido como desigualdad de Harnack. ✷
Problema 7.19 Obtener la función de Green para el problema de Dirichlet en
el semiespacio z > 0 en R3 .
Solución:

Sabemos que la solución fundamental en el espacio con polo en x0 es


1 1
ψx0 (x, y, z) = − =− p .
4πkx − x0 k 4π (x − x0 ) + (y − y0 )2 + (z − z0 )2
2

El borde del semiespacio z > 0 es el plano XY (z = 0). Sobre este plano,


1
ψx0 (x, y, 0) = − p .
4π (x − x0 ) + (y − y0 )2 + z02
2

Para que se anule Gx0 (x, y, 0) habrá que añadir a ψx0 una función armónica
Ux0 cuya contribución sobre el eje X cancele la de la solución fundamental. Por
tanto, es preciso que
1
Ux0 (x, y, 0) = p .
4π (x − x0 ) + (y − y0 )2 + z02
2

192
La manera de conseguirlo es situar una carga negativa en el punto simétrico
de x0 respecto al plano XY , que es x̃0 = (x0 , y0 , −z0 ). Ası́, como los puntos
del plano XY son equidistantes de x0 y de x̃0 , las contribuciones de ambos se
cancelarán, al tener signo opuesto,
1 1
Ux0 (x, y, z) = = p ,
4πkx − x̃0 k 4π (x − x0 ) + (y − y0 )2 + (z + z0 )2
2

que es una función armónica, ya que el polo (x0 , y0 , −z0 ) no pertenece al semi-
espacio z > 0.
Por tanto, la función de Green del semiespacio z > 0 es
1 1
Gx0 (x, y, z) = − .✷
4πkx − x̃0 k 4πkx − x0 k

La normal exterior al borde del semiespacio z > 0 es ν = −uz . Por tanto,

dGx0 (x, y, 0) dGx0 (x, y, 0) 1 z0


=− = ,
dν dz 2π ((x − x0 )2 + (y − y0 )2 + z 2 )3/2

y la solución del problema de Dirichlet, ∆u = 0, u(x, y, 0) = f (x, y), será

z0 ∞
Z ∞
f (x, y)
Z
u(x0 , y0 ) = dx dy 3/2
.
2π −∞ −∞ ((x − x0 ) + (y − y0 )2 + z 2 )
2

193

También podría gustarte