Ecuaciones en Derivadas Parciales - Universidad Nacional Del Litoral PDF
Ecuaciones en Derivadas Parciales - Universidad Nacional Del Litoral PDF
Ecuaciones en Derivadas Parciales - Universidad Nacional Del Litoral PDF
net
Presentación
Este apunte contiene las notas del curso Matemática Aplicada, curso de posgrado para
la Facultad de Ingenierı́a Quı́mica y la Facultad de Ingenierı́a y Ciencias Hı́dricas de la
Universidad Nacional del Litoral. Este curso se dicta en paralelo con el curso Ecuaciones
en Derivadas Parciales, optativa de la Licenciatura en Matemática Aplicada, y de pos-
grado para el doctorado y la maestrı́a en Matemática. En la sección Ejercicios de cada
capı́tulo hay problemas marcados con un asterisco. Estos problemas más teóricos están
pensados para los alumnos de estas últimas carreras. Los problemas marcados con dos
asteriscos son también para alumnos de estas carreras, pero son opcionales; si bien están
www.elsolucionario.net
relacionados con los contenidos del capı́tulo correspondiente, no son tan importantes para
el objetivo del curso. Los alumnos del doctorado o la maestrı́a en matemática deberı́an
hacerlos.
Este apunte fue escrito durante el dictado de los cursos en el primer cuatrimestre
de 2011, y será actualizado durante el dictado en el primer cuatrimestre de 2012. Cual-
quier sugerencia u observación sobre errores será agradecida por los autores del apunte y
también por los futuros alumnos del curso.
2
www.elsolucionario.net
Índice general
www.elsolucionario.net
1.7. Ejercicios . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
3
www.elsolucionario.net
www.elsolucionario.net
7.1.2. Un modelo para análisis de poblaciones o inventarios . . . . . . . 84
7.2. EDPs de primer orden con coeficientes variables . . . . . . . . . . . . . . 87
7.2.1. Parametrización preferida . . . . . . . . . . . . . . . . . . . . . . 90
7.2.2. Soluciones en forma paramétrica . . . . . . . . . . . . . . . . . . . 92
7.2.3. Consideraciones globales . . . . . . . . . . . . . . . . . . . . . . . 96
7.3. EDPs de primer orden en más variables . . . . . . . . . . . . . . . . . . . 97
7.4. Ejercicios . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 98
4
www.elsolucionario.net
www.elsolucionario.net
12.Ecuaciones en más variables independientes. 192
12.1. Difusión en dos dimensiones . . . . . . . . . . . . . . . . . . . . . . . . . 192
12.2. Ondas en dos dimensiones . . . . . . . . . . . . . . . . . . . . . . . . . . 194
12.3. Problema de autovalores en un rectángulo . . . . . . . . . . . . . . . . . 196
12.4. Problema de autovalores en el disco. Funciones de Bessel . . . . . . . . . 199
12.4.1. El caso general . . . . . . . . . . . . . . . . . . . . . . . . . . . . 199
12.4.2. El caso con simetrı́a radial . . . . . . . . . . . . . . . . . . . . . . 204
12.5. Ecuación de Laplace en un cilindro circular. Funciones de Bessel modificadas205
12.6. Ejercicios . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 212
5
www.elsolucionario.net
Capı́tulo 1
www.elsolucionario.net
x(t)
es la velocidad media en el intervalo de tiempo [0, t], en metros por segundo.
t
x(b) − x(a)
(con 0 ≤ a < b) es la velocidad media en el intervalo de tiempo [a, b],
b−a
en metros por segundo.
x(a + h) − x(a)
es la velocidad media en el intervalo de tiempo que dura h segundos
h
y comienza en t = a, también en metros por segundo.
x(a + h) − x(a)
lı́m es la velocidad instantánea a tiempo t = a, en metros por
h→0 h
segundo. Si definimos
x(t + h) − x(t)
v(t) = lı́m = x′ (t)
h→0 h
entonces v(t) es la velocidad instantánea a tiempo t.
Por el teorema fundamental del cálculo
Z b Z b
v(t) dt = x′ (t) dt = x(b) − x(a).
a a
6
www.elsolucionario.net
m(R) g
El cociente es la densidad media del cuerpo en la región R, en .
vol(R) cm3
Si B(P, r) denota la bola con centro P y radio r, entonces la densidad puntual f (P )
se define como
m B(P, r)
f (P ) = lı́m ,
r→0 vol B(P, r)
Ahora no resulta tan obvio como en el caso del desplazamiento y la velocidad, pero
en cálculo siempre nos enseñaron (por decreto) que si f (x) es la densidad puntual
www.elsolucionario.net
de un objeto que ocupa una región R, entonces su masa se calcula como
ZZZ
m(R) = f (x) dvol,
R
¿no es cierto? Para afirmar que integrando la f (x) que se obtiene como en el punto
anterior se recupera la masa m(R), harı́a falta un teorema, ¿no les parece?
o sea, Z b
m(b − a) ≤ f (x) dx ≤ M (b − a),
a
y por lo tanto Z b
1
m≤ f (x) dx ≤ M.
b−a a
7
www.elsolucionario.net
Z b
1
Además, si llamamos f = f (x) dx, observemos lo que ocurre al integrar esta
b−a a
constante f sobre el intervalo [a, b]:
Z b Z b Z b
1
f dx = (b − a)f = (b − a) f (x) dx = f (x) dx.
a b−a a a
www.elsolucionario.net
que las dos regiones sombrea-
das tengan la misma área.
a b
En base a estas observaciones:
Z b
1
Llamaremos valor medio de f sobre el intervalo [a, b] a la expresión f (x) dx
b−a a
Si ahora consideramos el valor medio de f en un intervalo [a, a+h] (h > 0), obtenemos
que Z
1 a+h
mh := mı́n f ≤ f (x) dx ≤ máx f =: Mh .
[a,a+h] h a [a,a+h]
Pero como la función es continua, lı́m Mh = f (a) y también lı́m mh = f (a), y por lo
h→0 h→0
tanto: Z a+h
1
lı́m+ f (x) dx = f (a).
h→0 h a
8
www.elsolucionario.net
En palabras, si una sucesión de intervalos está metida en otra sucesión que se encoge
a un punto entonces la sucesión de promedios tiende al valor de la función en ese
punto.
Las demostraciones rigurosas de estos teoremas quedan como ejercicio (ver 1.5– 1.9),
y se sugiere hacerlas utilizando un argumento de tipo ε − δ. El argumento previo con
máximos y mı́nimos fue realizado para ilustración del resultado, pero no es un buen
camino para obtener una demostración rigurosa.
Ejemplos:
Z x0 +h
1
f (x) dx → f (x0 ) cuando h → 0+ .
3h x0 −2h
Z
1 x0 +2h
f (x) dx → f (x0 ) cuando h → 0+ .
h x0 +h
www.elsolucionario.net
Z x0 +h
1
f (x) dx → f (x0 ) cuando h → 0+ .
h − h2 x0 +h2
9
www.elsolucionario.net
Ejemplos:
Si B(x, r) denota la bola con centro en x y radio r en R3 , entonces |B(x, r)| = 34 πr3 .
Si D(x, r) denota el disco con centro en x y radio r en R2 , entonces |D(x, r)| = πr2 .
www.elsolucionario.net
RRR
R
f dvol
mı́n f ≤ ≤ máx f si R es una región de R3
R RR |R| R
f dA
mı́n f ≤ R ≤ máx f si R es una región de R2
R RR |R| R
f dσ
mı́n f ≤ S ≤ máx f si S es una superficie en R3
S R |S| S
f ds
mı́n f ≤ C ≤ máx f si C es una curva en R2 o en R3
C |C| C
y más aún,
RRR ZZZ ZZZ
R
f dvol
Si f := entonces f dvol = f dvol (R región de R3 )
RR |R| ZZ R ZZ R
R
f dA
Si f := entonces f dA = f dA (R región de R2 )
RR |R| Z ZR ZZ R
S
f dσ
Si f := entonces f dσ = f dσ (S superficie en R3 )
R |S| Z S Z S
C
f ds
Si f := entonces f ds = f ds (C curva en R2 o R3 )
|C| C C
En todos estos casos diremos que f es el valor medio de f sobre R, S o C según corres-
ponda.
Un resultado análogo al Teorema 1.1 en más dimensiones es el siguiente (lo enunciamos
para R3 pero lo mismo vale en R2 cambiando B(x, r) por D(x, r) y dvol por dA).
10
www.elsolucionario.net
En otras palabras, si {rn } es una sucesión de radios (rn > 0) que tiende a cero, entonces
ZZZ
1
f dvol −→ f (x0 ) cuando n → ∞.
|B(x0 , rn )| B(x0 ,rn )
entonces necesariamente
m(B(x, r))
f (x) = lı́m+
r→0 vol(B(x, r))
Resultados análogos al Teorema 1.2 son los siguientes:
Teorema 1.5. Sea f una función definida en R3 , continua en x0 . Si {Rn }∞
n=1 es una
3
sucesión de regiones de R tales que Rn ⊂ B(x0 , rn ) y rn → 0, entonces
www.elsolucionario.net
ZZZ
1
f dvol −→ f (x0 ) cuando n → ∞.
|Rn | Rn
Las regiones de integración son cubos Qh de lado h (volumen h3 ) que tienen un vértice en
el punto x0 = (1, 2, 3). El cubo Qh está entonces contenido en la bola B((1, 2, 3), 2h) de
centro (1, 2, 3) y radio 2h. Al hacer tender h a cero, obtenemos el valor de f en el punto
(1, 2, 3),
f (1, 2, 3) = 12 + 22 + 32 = 14.
Teorema 1.7. Sea f una función definida en R3 , continua en x0 . Si {Sn }∞ n=1 es una
sucesión de superficies de R3 tales que Sn ⊂ B(x0 , rn ) y rn → 0, entonces
ZZ
1
f dσ −→ f (x0 ) cuando n → ∞.
|Sn | Sn
donde Qh son los cubos del ejemplo anterior, y ∂Qh denota la superficie que recubre a
Qh . Como el área de ∂Qh es 6h2 tenemos que
ZZ ZZ
1 1
lı́m f (x, y, z) dσ = 6 lı́m 2 f (x, y, z) dσ = 6f (1, 2, 3),
h→0 h2 ∂Qh h→0 6h ∂Qh
puesto que los cubos Qh cumplen las mismas hipótesis que antes.
11
www.elsolucionario.net
1.5. Aplicaciones
En base a las observaciones de la sección anterior, podemos demostrar que:
RRR
Teorema 1.9. Si f es una función continua en una región Ω de R3 y B(x,r)
f dvol = 0
para toda bola B(x, r) contenida en Ω, entonces f ≡ 0.
www.elsolucionario.net
h→0 |Qh | Qh
RRR
Observación 1.11. Es importante notar que si uno sabe que Q
f dvol = 0 sobre una
región Q, esto no implica necesariamente que f RRR = 0 en todos los puntos de Q.
Lo que afirmamos anteriormente es que si Q
f dvol = 0 para todos los cubos Q
contenidos en una región Ω entonces f es cero en todos los puntos de Ω. Por ejemplo,
la función f (x, y, z) = xy tiene integral cero sobre el cubo [−h, h] × [−h, h] × [−h, h] y
sin embargo no es cero en todos los puntos. También la función sen(x) tiene integral cero
sobre el intervalo [−π, π], y sin embargo no es idénticamente cero.
RRR
Observación 1.12. Sin embargo, si uno sabe que Q
f dvol = 0 sobre una región Q, y
también sabe que la función f no cambia de signo, por ejemplo, si se sabe que f (x) ≥ 0
para todo x ∈ Q, entonces sı́ se puede concluir que f (x) = 0 para todo x ∈ Q.
Si definimos (
1
si − h ≤ x ≤ h
χh (x) = 2h
0 en c.o.c.
Z x0 +h Z ∞
1
Entonces f (x) dx = f (x) χh (x0 − x) dx, y luego
2h x0 −h −∞
Z ∞
lı́m f (x) χh (x0 − x) dx = f (x0 ).
h→0 −∞
12
www.elsolucionario.net
5 5
h = 1/8 h = 1/8
h = 1/4 h = 1/4
4 h = 1/2 4 h = 1/2
h=1 h=1
3 3
2 2
1 1
0 0
−1 −1
−2 −1 0 1 2 −2 −1 0 1 2
Figura 1.1: Gráfico de χh (x) (izquierda) y de ρh (x) (derecha) para varios valores de h.
Z ∞
La integral f (x) χh (x0 − x) dx se llama convolución de f y χh , y se denota por
−∞
f ∗ χh (x0 ), es decir
Z ∞
www.elsolucionario.net
f ∗ χh (x0 ) = f (x) χh (x0 − x) dx −→ f (x0 ) (cuando h → 0+ ).
−∞
13
www.elsolucionario.net
Las sucesiones de funciones que se concentran alrededor del origen y tienen esta pro-
piedad se denominan aproximaciones de la identidad, pues al hacer la convolución con
una función f y tomar lı́mite, se recupera el valor de f . También se dice usualmente que
tienden a la delta de Dirac.
1.7. Ejercicios
1.1. Dar ejemplos de funciones continuas que no sean idénticamente nulas en el dominio
indicado, y cuyas integrales den cero.
(a) Intervalo [0, 1] de R.
www.elsolucionario.net
Z 1+h Z 1+h Z 1+h
1
(a) lı́m 3 f (x, y, z) dz dy dx
h→0 h 1 1 1
Z 1+h Z 1+h Z 1+h
1
(b) lı́m 3 f (x, y, z) dz dy dx
h→0 h 1−h 1−h 1−h
Z 1+h Z 1+h Z 1+h
1
(c) lı́m 2 f (x, y, z) dz dy dx
h→0 h 1 1 1
Z 2+2h Z −1+h Z h
1
(d) lı́m 3 f (x, y, z) dz dy dx
h→0 h 2+h −1−h 0
Z −1+h Z 4+h
1
(e) lı́m 2 f (x, y, 0) dy dx
h→0 h −1 4
Z −1+h Z 4+h
1
(f ) lı́m 2 f (x, y, 2) dy dx
h→0 h −1 4
ZZZ
1
(g) lı́m+ 3 f dvol
r→0 r B (−1,0,1),r
Z
1
(h) lı́m+ f ds, donde Cr denota la circunferencia de radio r y centro (1, 1, 0) conte-
r→0 r Cr
nida en el plano xy.
1.3. Demostrar que si f y g son dos funciones continuas en R3 y
ZZZ ZZZ
f dvol = g dvol,
R R
para toda región R del espacio R3 , entonces f (x0 ) = g(x0 ) para todo x0 ∈ R3 .
14
www.elsolucionario.net
* 1.5. Sea f una función definida en R, integrable sobre todo intervalo acodato. Demos-
trar que si f es continua en x0 , entonces
Z x0 +h
1
lı́m+ f (x) dx = f (x0 ).
h→0 2h x0 −h
* 1.6. Sea f una función definida en R, integrable sobre todo intervalo acodato. De-
mostrar que si f es continua en x0 , e {In }∞
n=1 es una sucesión de intervalos para la que
existe una sucesión de números positivos {hn }∞n=1 tales que In ⊂ [x0 − hn , x0 + hn ] para
n = 1, 2, . . . y hn → 0 cuando n → ∞ entonces
www.elsolucionario.net
Z
1
lı́m f (x) dx = f (x0 ).
n→∞ longitud(In ) I
n
* 1.7. Sea f una función definida en R3 , integrable sobre toda bola. Demostrar que si f
es continua en x0 , entonces
ZZZ
1
f dvol −→ f (x0 ) cuando r → 0.
|B(x0 , r)| B(x0 ,r)
* 1.8. Sea f una función definida en R3 , integrable sobre todo conjunto abierto acotado y
continua en x0 . Demostrar que si {Rn }∞
n=1 es una sucesión de conjuntos abiertos no-vacı́os
3
de R tales que Rn ⊂ B(x0 , rn ) para n = 1, 2, . . . y rn → 0 cuando n → ∞, entonces
ZZZ
1
lı́m f dvol = f (x0 ).
n→∞ |Rn | Rn
* 1.9. Sea f una función continua en R3 . Si {Sn }∞ n=1 es una sucesión de superficies
(suficientemente suaves, donde esté definida la integral) de R3 tales que Sn ⊂ B(x0 , rn )
para n = 1, 2, . . . y rn → 0 cuando n → ∞, entonces
ZZ
1
lı́m f dσ = f (x0 ).
n→∞ |Sn | Sn
Interesante: las superficies podrı́an ser bastante raras, y “grandes” de manera que |Sn | →
∞ y aún ası́ el lı́mite ser f (x0 ).
* 1.10. Dada una función f ∈ L2 (Ω), con Ω un subconjunto medible Racotado de Rd
1
(d ∈ N), demostrar que el valor medio de f sobre Ω definido por f = |Ω| Ω
f es el valor
R 2
de t que minimiza la función φ : R → R dada por φ(t) = Ω |f (x) − t| dx.
15
www.elsolucionario.net
* 1.11. Sea Ω unR conjunto abierto no vacı́o de Rd y sea f una función continua en Ω.
Demostrar que si Ω f φ dx = 0 para toda φ : Ω → R de la forma
Demostrar que:
ρ ∈ C ∞ (Rd ),
www.elsolucionario.net
ρ(x) = 0 para todo x ∈ Rd con |x| ≤ 1.
d
** 1.14. Sea Ω unR conjunto abierto no vacı́o de∞R y sea f una función continua en Ω.
Demostrar que si Ω f φ dx = 0 para toda φ ∈ C0 (Ω) entonces f ≡ 0 en Ω.
Aclaraciones y ayuda:
16
www.elsolucionario.net
Capı́tulo 2
El objetivo de este capı́tulo es solamente repasar lo que significa integrar sobre una
lı́nea o sobre una superficie, y cuáles son las fórmulas que se utilizan.
www.elsolucionario.net
Consideremos una curva C parametrizada por una función γ(t), para t ∈ [a, b]:
C
γ(t) γ(b)
γ(a)
a b t
γ ′ (t) = x′1 (t)i + x′2 (t)j + x′3 (t)k
Si por ejemplo f representa una densidad lineal, que es la masa por unidad de longitud
de un alambre que ocupa la curva C, entonces dicha integral es la masa del alambre.
17
www.elsolucionario.net
Cuando se pretende integrar una función vectorial F sobre una curva C, en general
se entiende que se quiere integrar la componente de F tangencial a la curva (que es una
función escalar), es decir Z Z
F · ds = F · τ ds,
C C
donde τ denota un vector tangencial a la curva C con una orientación elegida. Queda
claro que si cambiamos la orientación del vector τ entonces cambia el signo del resultado
de la integral.
Si γ parametriza la curva dada, entonces una elección obvia del vector tangente es
γ′
τ = ′ . Por lo tanto
|γ |
Z Z b Z b
γ ′ (t) ′
F · ds = F (γ(t)) · ′ |γ (t)| dt = F (γ(t)) · γ ′ (t) dt.
C a |γ (t)| a
R
Si F es un campo de fuerzas, entonces C F · ds es el trabajo realizado por el campo
de fuerzas F sobre una partı́cula que recorre la curva C (en la dirección indicada por la
tangente elegida). Si V es un campo de velocidades y C es una curva cerrada, entonces
www.elsolucionario.net
R
C
V · ds es la circulación del campo V a lo largo de la curva C.
Observación 2.1. En muchos casos no hará falta recurrir a la parametrización ni a es-
cribir integrales complicadas, pues todo se simplifica antes de empezar a calcular. Veamos
un par de ejemplos.
Ejemplo 2.2. Consideremos un campo de fuerzas F dado por
F (x1 , x2 , x3 ) = −x2 i + x1 j,
y calculemos el trabajo ejercido por F sobre una partı́cula que recorre la circunferencia C
de radio uno y centro 0 sobre el plano xy, en sentido antihorario cuando lo vemos desde
el lado en que x3 > 0. Observamos que si el punto x1 i + x2 j pertenece a la circunferencia
entonces el vector −x2 i + x1 j es tangente a la misma (y es unitario, pues | − x2 i + x1 j| =
|x1 i + x2 j| = x21 + x22 = 1. Por lo tanto
F (x1 , x2 , x3 ) · τ = (−x2 i + x1 j) · (−x2 i + x1 j) = 1.
Por lo tanto, el trabajo es
Z Z Z
F · ds = F · τ ds = 1 ds = longitud(C) = 2π.
C C C
18
www.elsolucionario.net
v X(u, v)
Xv
Xu
d S
www.elsolucionario.net
c
a b u
Si Xu = xu i + yu j + zu k y Xv = xv i + yv j + zv k (donde los subı́ndices u y v denotan
derivadas con respecto a u y a v respectivamente), la fórmula para calcular el área de
esta superficie es ZZ Z bZ d
área(S) = dσ = |X u × X v | dv du
S a c
La integral de una función escalar f definida sobre S se calcula de la siguiente manera:
ZZ Z bZ d
f dσ = f X(u, v) |X u × X v | dv du
S a c
Z bZ d
= f (x(u, v), y(u, v), z(u, v))|X u × X v | dv du
a c
19
www.elsolucionario.net
Cuando se pretende integrar una función vectorial F sobre una superficie S, en general
se entiende que se quiere integrar la componente de F perpendicular a la superficie, es
decir, se quiere calcular ZZ
F · n dσ,
S
www.elsolucionario.net
través de la superficie
S: x1 = 3, 2 ≤ x2 ≤ 4, 1 ≤ x3 ≤ 6,
Ejemplo 2.6. Flujo del campo de Newton. El campo de Newton está dado por la fórmula
r
V = − 3 , donde r denota el vector posición r = xi + yj + zk y r = |r| es la magnitud
r p
de r, es decir r = x2 + y 2 + z 2 . Calculemos el flujo del campo de Newton hacia afuera
de la superficie esférica S de centro 0 y radio 1. Si r está sobre esta superficie, entonces
r = 1 y además, el vector normal exterior es exactamente el radio r. Luego
ZZ ZZ ZZ ZZ
r r·r
V · n dσ = − 3 · r dσ = − 3 dσ = − 1 dσ = −área(S) = −4π.
S S r S r S
20
www.elsolucionario.net
donde la elección del signo depende de la dirección del flujo que se quiere calcular. Por
lo tanto:
ZZ Z bZ d
F · n dσ = F X(u, v) · n |X u × X v | dv du
S a c
Z bZ d
Xu × Xv
=± F X(u, v) · |X u × X v | dv du
a c |X u × X v |
Z bZ d
=± F X(u, v) · X u × X v dv du
a c
2.3. Ejercicios
2.1. Calcular el trabajo ejercido por la fuerza F del Ejemplo 2.2 sobre una partı́cula que
recorre la circunferencia de centro 0 y radio R que se encuentra sobre el plano xy, en
sentido horario cuando lo vemos desde el lado en que x3 > 0.
2.2. Sea V (x1 , x2 , x3 ) = ω0 (−x2 i + x1 j) con ω0 una constante (velocidad angular). Sea
www.elsolucionario.net
CR = {(R cos t, R sen t, 0) : t ∈ [0, 2π)}, R > 0.
Bibliografı́a complementaria
[Marsden Tromba] Jerrold E. Marsden, Anthony J. Tromba, Cálculo Vectorial, Addison-
Wesley Iberoamericana.
21
www.elsolucionario.net
Capı́tulo 3
www.elsolucionario.net
En esta sección recordamos la definición de los operadores diferenciales más usuales.
Su significado fı́sico quedará claro en las secciones siguientes.
El operador nabla ∇ se define por
∂ ∂ ∂ ∂ ∂ ∂
∇= i+ j+ k o por ∇= i+ j+ k,
∂x ∂y ∂z ∂x1 ∂x2 ∂x3
22
www.elsolucionario.net
2
∂f ∂f ∂f ∂ 2 f ∂ 2f ∂ 2f
∇ f = div grad f = ∇ · (∇f ) = ∇ · i+ j+ k = + + .
∂x1 ∂x2 ∂x3 ∂x21 ∂x22 ∂x23
Si un campo escalar tiene Laplaciano nulo, se dice que es armónico.
www.elsolucionario.net
Si V es un campo vectorial C 2 , entonces la divergencia del rotor es cero, i.e.,
∇ · ∇ × V = 0.
∇ × ∇f = 0.
V = Ψ + Φ, con ∇ · Ψ = 0 y ∇ × Φ = 0.
23
www.elsolucionario.net
C2
(x01 , x02 , x03 )
C4 C5
C3
C1
C6
www.elsolucionario.net
ZZ 6 ZZ
X
V · n dσ = V · n dσ,
S i=1 Ci
ZZ ZZ
1 1 V1 (x01 + h1 , x2 , x3 ) − V1 (x01 , x2 , x3 )
V · n dσ = dx2 dx3 .
vol(R) C1 ∪C2 h2 h3 Q1 h1
24
www.elsolucionario.net
Q1
C2
(x01 , x02 , x03 )
Q2
C5
C3
C4
C1
C6
Q3
www.elsolucionario.net
Si definimos Q2 = [x01 , x01 + h1 ] × [x03 , x03 + h3 ] ⊂ R2 y Q3 = [x01 , x01 + h1 ] × [x02 , x02 + h2 ] ⊂
2
R , haciendo cuentas análogas obtenemos
ZZ ZZ
1 1 V1 (x01 + h1 , x2 , x3 ) − V1 (x01 , x2 , x3 )
V · n dσ = dx2 dx3
vol(R) S h2 h3 Q1 h1
ZZ
1 V2 (x1 , x02 + h2 , x3 ) − V2 (x1 , x02 , x3 )
+ dx1 dx3
h1 h3 Q2 h2
ZZ
1 V3 (x1 , x2 , x03 + h3 ) − V3 (x1 , x2 , x03 )
+ dx1 dx2 .
h1 h2 Q3 h3
25
www.elsolucionario.net
En todos los casos mencionados, el lı́mite se toma cuando (h1 , h2 , h3 ) → (0, 0, 0).
Notemos que entonces valen las siguientes relaciones:
ZZ
0 1
www.elsolucionario.net
∇·V (x ) = lı́m n · V dσ,
vol(R) ∂R
ZZ
0 1
∇ϕ(x ) = lı́m n ϕ dσ,
vol(R) ∂R
ZZ
0 1
∇×V (x ) = lı́m n × V dσ.
vol(R) ∂R
26
www.elsolucionario.net
Pi
Pj
www.elsolucionario.net
(siguiendo ahora con un error menor a ε vol(Ω)). En el miembro izquierdo, todos los
trozos de ∂Ωi que son interiores a Ω aparecen dos veces y con normales opuestas. No
ocurre lo mismo con las caras de ∂Ωi que no son interiores a Ω sino que son parte de
la frontera de Ω. Por consiguiente, en la suma de la izquierda, todas las integrales de
superficie calculadas
RR sobre superficies interiores a Ω se cancelan mutuamente y entonces
esa suma es ∂Ω V · n dσ. En el lı́mite tenemos
ZZ ZZZ
V · n dσ = ∇ · V dvol.
∂Ω Ω
se cumple siempre que V sea diferenciable y sus derivadas sean continuas en Ω = Ω ∪ ∂Ω.
Por ejemplo, si V es el campo de Newton
r
V (r) = − ,
r3
p
con r = (x1 , x2 , x3 ) y r = x21 + x22 + x23 , que es diferenciable en R3 − {0} (tiene una
singularidad en el origen), la igualdad (3.1) vale siempre que Ω sea una región tal que Ω
no toca al origen. Por ejemplo, la igualdad (3.1) vale en la esfera con centro en (1, 1, 1) y
radio 1, pero no vale en la esfera centrada en el origen con radio 1. Tampoco vale en la
esfera centrada en el (1, 0, 0) y radio 1, aunque
sı́ vale en la esfera ahuecada {(x1 , x2 , x3 ) ∈
R3 : 1 < x21 + x22 + x23 < 4} = B (0, 0, 0), 2 − B (0, 0, 0), 1 .
Observación 3.3. Si V es un campo vectorial C 1 (R3 ) y ∇ · V = 0 en R3 , entonces
ZZ ZZZ
V · n dσ = ∇ · V dvol = 0.
∂Ω Ω
27
www.elsolucionario.net
Es decir, para toda superficie cerrada ∂Ω que encierre una región de R3 donde V tiene
divergencia nula, se cumple que el flujo de V total saliente/entrante a través de ∂Ω es nulo.
Esto no nos dice que la velocidad V sea nula, sino que el balance es nulo. Si V representa
la velocidad de un fluido, entonces ∇ · V = 0 indica que el flujo es incompresible.
www.elsolucionario.net
C 1 , entonces
n
ZZ Z
(∇ × V ) · n dσ = V · ds, S
S C
τ C
donde la normal n a S y la dirección de recorrido de C indicada por la elección de τ se
toman de modo que se satisfaga la regla de la mano derecha.
Observación 3.4. La primera observación que se puede hacer de la fórmula de Stokes,
justifica el nombre de rotor o rotacional para el vector ∇ × V , puesto que la anulación
del rotor implica la anulación de la circulación sobre cualquier curva cerrada.
Las consecuencias más profundas del teorema de Stokes se dan en la Teorı́a de Poten-
cial.
Dado un campo vectorial V , se dice que el campo escalar ϕ es un potencial para V si
satisface
∇ϕ = V .
No es cierto que todo campo vectorial admita un potencial, por ejemplo, es fácil ver
que el gradiente de cualquier potencial suficientemente suave ϕ es irrotacional. Es decir
∇ × ∇ϕ = 0.
Esto nos dice que para que un campo vectorial admita un potencial, es necesario que
dicho campo vectorial sea irrotacional:
28
www.elsolucionario.net
www.elsolucionario.net
Z Z 1
F · ds = F (γ(t)) · γ ′ (t) dt
C
Z0 1
= ∇ϕ(γ(t)) · γ ′ (t) dt
0
Z 1
d
= ϕ(γ(t)) dt
0 dt
= ϕ(γ(1)) − ϕ(γ(0))
= ϕ(Q) − ϕ(P ).
Por lo tanto, si F admite un potencial, el trabajo depende solamente de los puntos inicial
y final, y por lo tanto es conservativo.
Donde COP es una curva suave (a trozos) que une O con P . Como F es conservativo,
el resultado da lo mismo a través de cualquier curva COP que se elija. De esta manera
1
El dominio de suavidad de una función es el dominio donde la función está definida y es C 1 . Por
ejemplo, el dominio de suavidad del campo de Newton es R3 − {0}.
29
www.elsolucionario.net
∂ϕ
resulta ∇ϕ = F . Veamos como ejemplo que (x, y, z) = F3 (x, y, z), el resto de las
∂z
igualdades se obtiene de manera análoga (utilizando otra curva C).
(x, y, z)
C3
(0, 0, 0)
C1
(x, 0, 0)
C2 (x, y, 0)
www.elsolucionario.net
Utilicemos como curva C la de la figura, entonces
Z
ϕ(x, y, z) = F · τ ds
C
Z Z Z
= F · i ds + F · j ds + F · k ds
C1 C2 C3
Z Z Z
= F1 ds + F2 ds + F3 ds
C1 C2 C3
Z x Z y Z z
= F1 (s, 0, 0) ds + F2 (x, s, 0) ds + F3 (x, y, s) ds.
0 0 0
Luego,
∂ϕ
(x, y, z) = F3 (x, y, z).
∂z
En conclusión:
Todo campo conservativo es el gradiente de algún campo escalar.
30
www.elsolucionario.net
para toda curva cerrada simple en el dominio. Entonces F es conservativo y por consi-
guiente existe ϕ tal que F = ∇ϕ. El campo escalar ϕ es un potencial para F .
3.5. Ejercicios
3.1. Sea f (r) = 1r el potencial Newtoniano cuyo dominio es R3 \ {0}.
p
(Notación: r = x1 i + x2 j + x3 k y r = x21 + x22 + x23 )
3.2. Demostrar que el gradiente de cualquier función escalar ϕ que tiene derivadas de
segundo orden continuas, es un campo vectorial irrotacional. Más precisamente, demostrar
que
∇ × ∇ϕ = 0,
www.elsolucionario.net
para cualquier función escalar ϕ con derivadas de segundo orden continuas.
r
V (r) = − ,
r3
a través de una superficie esférica centrada en (1, 1, 1) y de radio uno. (Usar Teorema de
Gauss).
¿Cuál es el flujo a través de cualquier superficie cerrada suave que no deje el origen
en el interior?
r
3.5. Sea V (r) = − . Recuerde que según se resolvió en el ejercicio 2.3, el flujo de V
r3
a través de la superficie esférica centrada en 0 = (0, 0, 0) y radio R, hacia el infinito es
−4π.
Calcular el flujo de V a través de cualquier superficie cerrada suave que contenga a 0
en su interior.
3.6. Sea S una superficie cerrada y suave en el espacio, y sea ni la i-ésima componente
del vector normal exterior a S. Demostrar que entonces
ZZ
ni dσ = 0.
S
31
www.elsolucionario.net
3.7. Sea S una superficie cerrada y suave en el espacio, sea n el vector normal exterior a
S, y r el vector posición (o función identidad de R3 ). Demostrar que entonces
ZZ
1
r · n dσ = V,
3 S
donde V es el volumen de la región rodeada por S.
3.8. Sea ϕ(x1 , x2 , x3 ) = ln(x21 + x22 ). Determinar su dominio y el de su gradiente V = ∇ϕ.
¿Es V conservativo? ¿Es incompresible? ¿Es ϕ armónico?
Bibliografı́a complementaria
[Apostol] Tom M. Apostol, Calculus, Vol. 2., Editorial Reverté, 2005.
[Arfken-Weber] Arfken, G.B., Weber, H.J., Mathematical Methods For Physicists, HARCOUT-
Academic Press, 2001.
[Marsden-Tromba] Jerrold E. Marsden, Anthony J. Tromba, Cálculo Vectorial, Addison-
www.elsolucionario.net
Wesley Iberoamericana.
32
www.elsolucionario.net
Capı́tulo 4
Coordenadas Generalizadas en el
Espacio
www.elsolucionario.net
describirse como la intersección de tres superficies: una de cada familia.
33
www.elsolucionario.net
F1 F2 F3
www.elsolucionario.net
P
Q3 = q3
Q1 = q1
34
www.elsolucionario.net
r = 3, 0 ≤ θ < 2π, 0 ≤ ϕ ≤ π,
www.elsolucionario.net
x3 es de π/4 radiantes o 45 grados.
Si consideramos las coordenadas cilı́ndricas, la superficie
www.elsolucionario.net
Un vector tangente a esta curva en el punto (q1 , q2 , q3 ) (t = q1 ) está dado por
d ∂X ∂X1 ∂X2 ∂X3
X(t, q2 , q3 )|t=q1 = (q1 , q2 , q3 ) = (q1 , q2 , q3 )i+ (q1 , q2 , q3 )j + (q1 , q2 , q3 )k.
dt ∂q1 ∂q1 ∂q1 ∂q1
Ası́, los vectores tangentes a las curvas coordenadas que pasan por el punto que tiene
∂X
coordenadas (q1 , q2 , q3 ) son tres, a saber: (q1 , q2 , q3 ), i = 1, 2, 3.
∂qi
Cada par de estos vectores (se pueden armar tres pares diferentes) determina el plano
tangente a una de las tres superficies coordenadas que pasan por dicho punto. Más pre-
cisamente:
∂X ∂X
El plano paralelo a ya que pasa por (q1 , q2 , q3 ) es el plano tangente a la
∂q2 ∂q3
superficie coordenada {q1 = constante} que pasa por dicho punto.
∂X ∂X
El plano paralelo a ya que pasa por (q1 , q2 , q3 ) es el plano tangente a la
∂q1 ∂q3
superficie coordenada {q2 = constante} que pasa por dicho punto.
∂X ∂X
El plano paralelo a ya que pasa por (q1 , q2 , q3 ) es el plano tangente a la
∂q1 ∂q2
superficie coordenada {q3 = constante} que pasa por dicho punto.
∂X
∂q2
∂X
∂q1
Q3 = constante
36
www.elsolucionario.net
En los ejemplos de coordenadas esféricas y cilı́ndricas que consideramos resulta que los
∂X
vectores (q1 , q2 , q3 ), i = 1, 2, 3, son perpendiculares entre sı́ (ver ejercicio 4.2), y por
∂qi
∂X
lo tanto el vector que no es paralelo al plano tangente de la superficie coordenada,
∂qi
es necesariamente perpendicular. Más precisamente
∂X
En cada punto, el vector es perpendicular al plano tangente a la superficie
∂q1
coordenada {q1 = constante} que pasa por dicho punto.
∂X
En cada punto, el vector es perpendicular al plano tangente a la superficie
∂q2
coordenada {q2 = constante} que pasa por dicho punto.
∂X
En cada punto, el vector es perpendicular al plano tangente a la superficie
∂q3
coordenada {q3 = constante} que pasa por dicho punto.
www.elsolucionario.net
En general decimos que un vector es perpendicular a una superficie en un punto si es
perpendicular al plano tangente a dicha superficie en ese punto, es por eso que las tres
afirmaciones anteriores pueden re-escribirse de la siguiente manera (resumida):
∂X
El vector es perpendicular a la superficie {q1 = constante}.
∂q1
∂X
El vector es perpendicular a la superficie {q2 = constante}.
∂q2
∂X
El vector es perpendicular a la superficie {q3 = constante}.
∂q3
∂X
Finalmente observamos que el vector apunta en la dirección en que qi aumenta.
∂qi
37
www.elsolucionario.net
∂X
donde hemos definido hi =
, y como antes ∂X = ∂X1 i + ∂X2 j + ∂X3 k.
∂qi ∂qi ∂qi ∂qi ∂qi
Ejemplo 4.3. Consideremos las coordenadas esféricas y calculemos h2 y a2 . Para realizar
esto debemos utilizar las fórmulas (4.2) y derivar con respecto a q2 = θ:
∂X1
= −r sen ϕ sen θ
∂θ
∂X2 ∂X
= r sen ϕ cos θ =⇒ = r sen ϕ − sen θ i + cos θ j
∂θ ∂θ
∂X3
=0
∂θ
Por lo tanto h2 = r sen ϕ y a2 = − sen θ i + cos θ j.
www.elsolucionario.net
Sea γ(t) una curva en el espacio que describe (por ejemplo) la trayectoria de una
partı́cula. Si representamos la curva en coordenadas cartesianas ortogonales γ(t) = x1 (t)i+
x2 (t)j + x3 (t)k, para calcular la longitud de un arco cualquiera de la curva, necesitamos
dγ
conocer (t) ya que la longitud del arco que une γ(a) con γ(b) está dada por:
dt
Z Z b Z b r
dr dx1 2 dx2 2 dx3 2
ds = dt (t) dt =
+ + dt .
C a a | dt dt
{z dt }
ds
Supongamos ahora que qj (t) = Qj γ(t) son las coordenadas generalizadas de la
dγ 2
curva, y tratemos de hallar una fórmula para (t) en términos de las coordenadas
dt
qj (t). Estamos considerando entonces que
∂X ′ ∂X ′ ∂X ′
γ ′ (t) = q1 (t) + q2 (t) + q (t)
∂q1 ∂q2 ∂q3 3
= h1 a1 q1′ (t) + h2 a2 q2′ (t) + h2 a3 q3′ (t)
p ′ ∂X
Además |γ ′ (t)| = γ (t) · γ ′ (t) y como = hi ai ,
∂qi
γ ′ (t) · γ ′ (t) = h21 (q1′ (t))2 + h22 (q2′ (t))2 + h23 (q3′ (t))2 ,
38
www.elsolucionario.net
www.elsolucionario.net
X i
con hℓ = = t .
∂qℓ i=1
∂q ℓ
ds1 = h1 dq1 .
curva intersección
C
Q3 = constante
Q2 = constante
Para una curva como C en el dibujo, el diferencial de longitud es ds1 = h1 dq1 . Si q1 y
q3 son constantes tenemos ds2 = h2 dq2 , y si q1 y q2 son constantes tenemos ds3 = h3 dq3 .
39
www.elsolucionario.net
Observación 4.5. Notar que para obtener la longitud de cualquier paralelo, se debe
calcular una integral entre 0 y 2π. Sin embargo los resultados dependen de q1 = r y
q3 = ϕ. Aquı́ queda un poco más claro lo que significa h2 : mide cuánto se estira el
segmento [0, 2π] cuando se deforma para describir un paralelo de la esfera.
www.elsolucionario.net
De un modo análogo, si Q2 = constante, dσ13 = ds1 ds3 = h1 h3 dq1 dq3 .
Y si Q1 = constante, dσ23 = ds2 ds3 = h2 h3 dq2 dq3 .
40
www.elsolucionario.net
www.elsolucionario.net
ψ(q1 , q2 , q3 ) = Ψ X1 (q1 , q2 , q3 ), X2 (q1 , q2 , q3 ), X3 (q1 , q2 , q3 ) .
Luego, el gradiente del campo escalar es
∂Ψ ∂Ψ ∂Ψ
∇Ψ(x1 , x2 , x3 ) = (x1 , x2 , x3 ) i + (x1 , x2 , x3 ) j + (x1 , x2 , x3 ) k.
∂x1 ∂x2 ∂x3
Por otro lado, por la regla de la cadena,
∂ψ ∂Ψ ∂X1 ∂Ψ ∂X2 ∂Ψ ∂X3
= + +
∂qi ∂x1 ∂qi ∂x2 ∂qi ∂x3 ∂qi
∂X
= ∇Ψ · = ∇Ψ · (hi ai ),
∂qi
∂X ∂X1 ∂X2 ∂X3 ∂X
donde hemos usado que = i+ j+ k, y hi = (ver Sección 4.3).
∂qi ∂qi ∂qi ∂qi ∂qi
Finalmente,
∂ψ
= hi ∇Ψ · ai ,
∂qi
o, lo que es lo mismo
1 ∂ψ
∇Ψ · ai = .
hi ∂qi
41
www.elsolucionario.net
4.8. Ejercicios
4.1. Considerar las coordenadas esféricas y cilı́ndricas. Para cada caso:
www.elsolucionario.net
Describir F1 , F2 , F3 ;
Calcular h1 , h2 , h3 ;
4.2. Verificar que para las coordenadas cilı́ndricas y esféricas se cumple que
∂X ∂X
· = 0, si i 6= j.
∂qi ∂qj
4.3. Describir todos los tipos de curvas coordenadas en los sistemas cartesianos, esféricos
y cilı́ndricos.
4.4. Utilizando coordenadas cilı́ndricas, calcular la longitud de la hélice θ(t) = t, ρ(t) =
R, z(t) = t, 0 ≤ t ≤ 2π.
4.5. Utilizando coordenadas esféricas. Calcular el área del casquete esférico de apertura
ϕ = π/4.
4.6. Utilizando coordenadas esféricas, calcular el volumen del cono esférico de apertura
ϕ = π/4.
4.7. Sea V (x1 , x2 , x3 ) = ω0 (−x2 i + x1 j) con ω0 una constante (velocidad angular). Sea
42
www.elsolucionario.net
(c) Comparar el resultado del ı́tem anterior con el resultado del ejercicio 2.2 y explicar
la casualidad usando el teorema de Stokes.
Bibliografı́a complementaria
[Arfken-Weber] Arfken, G.B., Weber, H.J., Mathematical Methods For Physicists, HARCOUT-
Academic Press, 2001.
www.elsolucionario.net
43
www.elsolucionario.net
Capı́tulo 5
www.elsolucionario.net
Sea Ω una región del espacio R3 y sea ∂Ω su frontera, que suponemos suave a trozos.
Estamos interesados en describir (cuantitativamente) la evolución espacio-temporal de
una cantidad E en Ω. Por ejemplo, la energı́a térmica, la masa de un compuesto, etc. La
ley de conservación básica establece el hecho casi obvio siguiente:
Supondremos, para fijar ideas, que E es energı́a térmica. Más precisamente, supo-
nemos que si E = E(Ω; t) es la energı́a térmica en la región Ω a tiempo t, entonces
la densidad de energı́a térmica e = e(x; t) es la propiedad intensiva asociada a E y se
relacionan de la siguiente manera:
ZZZ
E(B(x, r); t)
E(Ω; t) = e(x; t) dvol ⇐⇒ e(x; t) = lı́m+ .
Ω r→0 vol(B(x, r))
44
www.elsolucionario.net
φ(x; t) se denomina flujo de energı́a térmica por unidad de área en el punto x a tiempo
t. Por lo tanto, si n es el vector normal a ∂Ω que apunta hacia fuera de Ω,
razón de inmigración ZZ
menos
=− φ(x; t) · n dσ
razón de emigración ∂Ω
a través de ∂Ω
La creación y aniquilación de la cantidad E en Ω (que en el caso de energı́a térmica
puede corresponder a reacciones quı́micas exotérmicas o endotérmicas respectivamente,
o a la circulación de una corriente eléctrica) está dada por una función f (x; t) definida
en Ω. La función f tiene signo arbitrario. En los puntos e instantes donde sea positiva
tendremos fuentes de energı́a, mientras que cuando sea negativa tendremos sumideros
de energı́a. La función f es una tasa de creación/aniquilación de energı́a por unidad de
volumen por unidad de tiempo, de manera que
razón de creación ZZZ
menos
= f (x; t) dvol
razón de desaparición
www.elsolucionario.net
Ω
dentro de Ω
Con estas tres magnitudes básicas, la ley de conservación y el Teorema de Gauss,
podemos establecer una relación cuantitativa precisa entre e, φ y f .
Sea x0 = (x01 , x02 , x03 ) un punto fijo en Ω y sea B = B(x0 , r) una bolita chica centrada
en x0 de manera que B(x0 , r) ⊂ Ω. Consideremos el balance en B.
ZZZ
razón de cambio d d
= E(B; t) = e(x; t) dvol
temporal de E en B dt dt B
ZZ
razón de inmigración/emigración
=− φ(x; t) · n dσ
a través de ∂B ∂B
ZZZ
razón de creación/desaparición
= f (x; t) dvol
dentro de B B
donde ∇ · φ es, claro, la divergencia del flujo térmico. Si e es una función suave, como la
derivación en el miembro izquierdo de (5.2) es con respecto a la variable tiempo, que no
es la de integración, y como estamos considerando una bola B fija que no depende de la
variable temporal, tenemos también que
ZZZ ZZZ
d ∂
e(x; t) dvol = e(x; t) dvol.
dt B B ∂t
45
www.elsolucionario.net
Esta igualdad debe cumplirse para toda bola B tal que B ⊂ Ω. Si e y φ son de clase C 1
y si f es continua, entonces la función entre corchetes es continua, y por lo visto en el
Capı́tulo 1 se tiene que cumplir que
∂
e(x; t) + ∇ · φ(x; t) − f (x; t) = 0
∂t
para todo x ∈ Ω y todo instante de tiempo t.
La ley de conservación (de la energı́a) es entonces la ecuación diferencial en derivadas
parciales de primer orden dada por
∂e
+ ∇ · φ = f. (5.3)
∂t
www.elsolucionario.net
Observación 5.1. Una ley de conservación similar se cumple si C representa una con-
centración, que es también una densidad de masa de una cierta sustancia, por unidad de
volumen:
∂C
+ ∇ · φ = f. (5.4)
∂t
En este caso φ es el flujo de masa por unidad de área, y f es un término fuente de masa
por unidad de volumen por unidad de tiempo.
Entonces
e = cρu.
La ley de conservación toma pues la forma
∂u
cρ + ∇ · φ = f,
∂t
donde c, ρ y f son generalmente datos del problema, y u, φ son las incógnitas.
Es claro que teniendo una sola ecuación, el problema de encontrar u y φ no está sufi-
cientemente determinado. Las ecuaciones de balance, o leyes de conservación tienen que
46
www.elsolucionario.net
www.elsolucionario.net
φ = −K∇u.
Sustituyendo la segunda en la primera, obtenemos la forma general de la ecuación de
difusión o ecuación del calor en el espacio tridimensional:
∂u
cρ − ∇ · (K∇u) = f. (5.5)
∂t
Recordemos que c, ρ y K son funciones del punto.
Cuando la conductividad térmica K es constante, la ecuación toma la forma
∂u 2
cρ − K∇ u = f.
∂t
Cuando no hay fuentes (ni sumideros) y también el calor especı́fico y la densidad son
constantes, tenemos 2
∂u 2 ∂ u ∂ 2u ∂ 2u
= k∇ u = k + +
∂t ∂x21 ∂x22 ∂x23
K
con k = la difusividad del material.
cρ
Observación 5.2. Cuando consideramos una concentración (de masa) C en lugar de
una temperatura, o densidad de energı́a térmica, la relación constitutiva que relaciona C
con el flujo φ se denomina Ley de Fick :
φ = −k∇C,
que es igual a la de Fourier. Incorporando esta relación constitutiva en la ley de conser-
vación obtenemos:
∂C
− ∇ · (k∇C) = f.
∂t
Obtuvimos la misma ecuación para la incógnita C que habı́amos obtenido anteriormente
para la incógnita u. Esta ecuación se denomina ecuación de difusión y modela tanto
difusión de calor como de masa.
47
www.elsolucionario.net
www.elsolucionario.net
∂u
y suele definirse directamente el flujo escalar φ = −K , interpretando que el flujo es
∂x1
hacia la derecha cuando φ > 0 y hacia la izquierda cuando φ < 0.
Llamando x a x1 (ya que es la única variable de interés) obtenemos
∂u ∂ ∂u
cρ − K = f,
∂t ∂x ∂x
que también podemos escribir
cρut − (Kux )x = f.
f
En el caso en que los coeficientes c, ρ y K son constantes, y llamando f˜ a
cρ
∂u ∂2u
− k 2 = f˜,
∂t ∂x
o, lo que es lo mismo
ut − kuxx = f˜.
Ecuación del calor bidimensional: Si consideramos una placa plana delgada, alinea-
da con el plano x1 x2 , aislada en su parte superior e inferior, es válido suponer que todas
las cantidades son constantes sobre rectas paralelas al eje x3 y obtenemos la ecuación:
∂u ∂ ∂u ∂ ∂u
cρ − K + K =f
∂t ∂x1 ∂x1 ∂x2 ∂x2
| {z }
=∇·(K∇u)
48
www.elsolucionario.net
o, lo que es lo mismo
ut − k uxx + uyy = f˜.
www.elsolucionario.net
u(x; 0) = f (x), x ∈ Ω.
Esta condición suele escribirse también de la siguiente manera:
u = f, x ∈ Ω, t = 0.
¿Es esta información suficiente para predecir la temperatura en el futuro? La respuesta
es no. Sabemos la distribución de temperatura inicial y sabemos que la temperatura
cambia acorde a la ecuación diferencial parcial (5.5). Pero aún nos falta saber qué ocurre
en el borde del objeto Ω. Sin esta información, no podemos predecir el futuro.
Los siguientes tipos de condiciones de borde son los más usuales:
49
www.elsolucionario.net
Flujo prescripto. En otras situaciones es posible suponer que se conoce (o está pres-
cripto) el flujo de calor en lugar de la temperatura, es decir
www.elsolucionario.net
Ley de enfriamiento de Newton. Cuando el objeto está en contacto con un fluido
en movimiento (aire, agua, aceite, etc.), entonces las condiciones hasta ahora vistas no
parecen del todo apropiadas. Por ejemplo, imaginemos un objeto que está caliente, y en
contacto con aire en movimiento. El calor saldrá del objeto, calentando el aire. Y el aire
se llevará el calor por el movimiento mismo. Experimentos muestran que el flujo de calor
que sale del objeto es proporcional a la diferencia de temperatura entre el objeto (u(x; t))
y la temperatura del fluido exterior (uB (x, t)). Esta condición de frontera se llama Ley
de enfriamiento de Newton. Si la queremos escribir en términos matemáticos resulta
∂u
− K(x) (x; t) = H[u(x; t) − uB (x, t)], x ∈ ∂Ω, t ≥ 0, (5.9)
∂n
donde la constante de proporcionalidad H se llama coeficiente de transferencia de calor
o coeficiente de convección. Esta condición de borde involucra una combinación lineal de
∂u
uy .
∂n
El coeficiente H en la ley de enfriamiento de Newton se determina experimentalmente.
Depende de las propiedades del objeto como ası́ también de las del fluido (inclusive de
la velocidad con la que el fluido se mueve). Si el coeficiente es muy pequeño, muy poca
energı́a fluye a través del borde. En el lı́mite H → 0, la ley de enfriamiento de Newton,
se convierte en la condición de borde aislado. Podemos pensar que la ley de Newton con
H 6= 0 representa la situación en que el borde no está perfectamente aislado.
Cuando H es grande, mucha energı́a fluye a través del extremo. En el lı́mite H → ∞,
la condición de borde se transforma la condición de temperatura prescripta u(x, t) =
uB (x, t). Esto puede verse fácilmente si dividimos (5.9) por H:
K(x) ∂u
− (x, t) = u(x, t) − uB (x, t),
H ∂n
y tomamos lı́mite cuando H → ∞.
50
www.elsolucionario.net
∂u
La condición de borde en que se relaciona el valor de la derivada normal con la
∂n
diferencia entre u y un dato uB se denomina condición de borde de tipo Robin.
0 u(0, t) = uB (0, t) K0 (0)ux (0, t) = φ(0, t) K0 (0)ux (0, t) = H[u(0, t) − uB (0, t)]
L u(L, t) = uB (L, t) −K0 (L)ux (L, t) = φ(L, t) −K0 (L)ux (L, t) = H[u(L, t) − uB (L, t)]
Es importante observar el cambio de signo: el flujo hacia afuera por unidad de área
es K0 (0)ux (0, t) en x = 0 y es −K0 (L)ux (L, t) en x = L.
www.elsolucionario.net
5.5. La ecuación de Laplace y la ecuación de Poisson
Cuando la conductividad térmica del material y la fuente f son independientes del
tiempo en la ecuación general de conducción del calor
∂u
cρ = ∇ · K∇u + f,
∂t
y nos proponemos encontrar soluciones estacionarias (independientes del tiempo) de esa
ecuación, al hacer ut = 0 tenemos la ecuación de Poisson
∇ · K∇u = −f
2
∇ u = −f˜ si K es constante.
∇ · K∇u = 0
2
∇ u=0 si K es constante.
ut + ∇ · φ = f.
En muchos casos la fuente f es función del punto del espacio, del instante t y también de
u misma: f = f (x, t, u).
51
www.elsolucionario.net
ut + φx = f,
www.elsolucionario.net
Transporte, advección: Relación constitutiva “el flujo es proporcional a la cantidad
que fluye”
φ = cu.
El flujo φ no depende explı́citamente de t ni de x, sólo a través de u. La incógnita u
satisface la ecuación diferencial en derivadas parciales de primer orden
ut + cux = f.
1
V
c x
52
www.elsolucionario.net
u(x, t) = F (x − ct)
para alguna función F de una variable, de clase C 1 (derivable con derivada continua). En
efecto, por la regla de la cadena
ut = F ′ (x − ct)(−c), ux = F ′ (x − ct)1,
y por lo tanto
ut + cux = 0.
www.elsolucionario.net
ut + cux = 0,
u(x, 0) = F (x), x ∈ R,
2
está dada por u(x, t) = F (x − ct). Por ejemplo, is F (x) = e−x , la solución es una
Gaussiana que se traslada a velocidad c:
2
u(x, t) = e−(x−ct) .
Se tiene advección con decaimiento cuando el término fuente f = −λu con λ > 0.
Advección no lineal: Cuando la relación constitutiva φ = ϕ(u) está dada por una
función no lineal ϕ de la variable u. La ecuación diferencial de primer orden toma ahora
la forma
ut + ϕ′ (u)ux = f.
53
www.elsolucionario.net
www.elsolucionario.net
Si V = V1 i + V2 j + V3 k, entonces ∇ V = ∇ V1 i + ∇ V2 j + ∇ V3 k.
Por otra parte, puesto que la divergencia de V es un campo escalar, podemos calcular
su gradiente para obtener ∇(∇ · V ). También está bien definido el rotor del rotor de V :
∇ × (∇ × V ). La fórmula que nos interesa es una relación entre estos tres operadores
vectoriales, que resumimos en el siguiente lema
Lema 5.4. Si V es un campo vectorial C 2 en el espacio, entonces
2
∇ V = ∇(∇ · V ) − ∇ × (∇ × V ).
O sea, el Laplaciano vectorial es el gradiente de la divergencia menos el rotor del rotor.
2 2
Demostración. Supongamos primero que V = V1 i, entonces ∇ V = ∇ V1 i, y por otro
lado
i j k
∂V
1 ∂ ∂ ∂
∇(∇ · V ) − ∇ × (∇ × V ) = ∇ −∇×
∂x1 ∂x1 ∂x2 ∂x3
V1 0 0
i j k
∂V ∂ ∂ ∂
1
=∇ − ∂x1 ∂x2 ∂x 3
∂x1
0 ∂V1 ∂V1
−
∂x3 ∂x2
∂ 2 V1 ∂ 2 V1 ∂ 2 V1
= i + j + k
∂x21 ∂x2 ∂x1 ∂x3 ∂x1
∂ 2V ∂ 2 V1 ∂ 2 V1 ∂ 2 V1
1
+ + i− j− k
∂x22 ∂x23 ∂x2 ∂x1 ∂x3 ∂x1
2 2
= ∇ V1 i = ∇ V .
54
www.elsolucionario.net
Volvamos ahora al sistema de Maxwell con la fórmula del lema. Calculamos el rotor
en ambos miembros de (M.4):
∂B ∂
∇ × (∇ × E) = −∇ × =− (∇ × B),
∂t ∂t
www.elsolucionario.net
∂ ∂ 2E
(∇ × B) = ε0 µ0 2 ,
∂t ∂t
∂ 2E
ε0 µ0 = −∇ × (∇ × E) = −∇ × (∇ × E) + ∇(∇ · E).
∂t2
Luego, por el Lema 5.4
2
−∇ × (∇ × E) + ∇(∇ · E) = ∇ E,
2 ∂ 2E
∇ E = ε0 µ0 .
∂t2
2 ∂ 2B
∇ B = ε0 µ0 .
∂t2
Más precisamente,
2 ∂ 2E
∇ E i = ε0 µ0 2 i ,
∂t i = 1, 2, 3,
2
2 ∂ Bi
∇ B i = ε0 µ0 ,
∂t2
son las ecuaciones de ondas que deben satisfacer las componentes de los campos eléctrico
y magnético.
55
www.elsolucionario.net
La función de onda ψ de una partı́cula bajo la acción de un potencial V tiene que satisfacer
la Ecuación de Schrödinger que, en forma simplificada luce de la siguiente manera
∂ψ 2
i + ∇ ψ = V.
∂t
www.elsolucionario.net
Esta ecuación se parece sólo en su aspecto exterior a una ecuación de difusión, porque el
factor complejo i cambia completamente el comportamiento de las soluciones.
5.9. Ejercicios
* 5.1. Sea Ω un abierto de Rd , y sea B una bola tal que su clausura está contenida en
Ω. Demostrar que si u : Ω × (0, T ) → R es C 1 (todas las derivadas parciales de primer
orden son continuas) entonces
Z Z
d
u(x, t) dx = ut (x, t) dx, ∀t ∈ (0, T ).
dt B B
5.2. Considerar una barra cilı́ndrica de longitud L alineada con el eje x, ocupando el
intervalo [0, L], suponer que su conductividad térmica y temperatura son constantes en
los planos paralelos al plano yz. Llamemos K(x) a la conductividad térmica a distancia
x del extremo izquierdo y u(x; t) a la temperatura a distancia x del extremo izquierdo.
Supongamos que el área de la sección transversal de la barra cilı́ndrica es A.
(a) Escribir sendas fórmulas para el flujo total de calor hacia fuera de la barra en x = 0
y en x = L.
(b) Escribir sendas fórmulas para la cantidad total de calor que fluyó hacia fuera de la
barra en x = 0 y en x = L por el perı́odo de tiempo desde t = 0 hasta t = T .
5.3. Supongamos que se unen dos barras cilı́ndricas de la misma forma, pero diferente
material en el plano x = x0 . Supongamos que la que ocupa el intervalo [0, x0 ] tiene
conductividad térmica K1 y la que ocupa el intervalo [x0 , L] tiene conductividad térmica
K2 . Se dice que estas dos barras están en contacto térmico perfecto si la temperatura es
continua en x = x0 :
u(x− +
0 , t) = u(x0 , t)
56
www.elsolucionario.net
y no se pierde nada de energı́a calórica en x = x0 (es decir, la energı́a calórica que fluye
hacia fuera de una barra fluye hacia adentro de la otra). ¿Qué ecuación matemática
representa la condición
la energı́a calórica que fluye hacia fuera de una barra fluye hacia adentro de la otra
en x = x0 ?
5.4. Considerar un baño conteniendo un volumen V de aceite, con calor especı́fico cf y
densidad de masa ρf , en el que se sumerge un cuerpo de metal con calor especı́fico cm ,
densidad ρm y conductividad térmica Km . Suponer que el baño se agita rápidamente de
modo que la temperatura del baño es aproximadamente uniforme en todo el recipiente
que lo contiene, y ésta es igual a la temperatura del cuerpo de metal en todo su contorno.
Supongamos que el baño está térmicamente aislado excepto donde está en perfecto con-
tacto térmico con el cuerpo de metal, donde puede enfriarse o calentarse dependiendo de
la temperatura del mismo. El objetivo de este problema es determinar una ecuación para
la temperatura del baño, que resultará a su vez una condición de borde para el cuerpo
de metal. Para hacerlo seguiremos los siguientes pasos.
(a) Llamemos T (t) a la temperatura del aceite. Escribir una fórmula que represente la
www.elsolucionario.net
energı́a térmica total del baño de aceite.
(b) Escribir una fórmula que represente la razón de cambio de la energı́a térmica total
del aceite.
(c) Sea Ω la región que ocupa el cuerpo de metal, y sea u(x; t) la temperatura en el
punto x del metal a tiempo t. Escribir una fórmula que represente el flujo de calor
hacia fuera del cuerpo de metal.
(d) Como la cantidad de calor que fluye hacia fuera del cuerpo de metal es la cantidad de
calor que ingresa al fluido, igualando las fórmulas de los ı́tems (b) y (c) obtenemos
una ecuación diferencial (ordinaria) para T (t) que está acoplada con la ecuación del
calor en el metal. Escribirla.
(e) Escribir un sistema de ecuaciones diferenciales para T (t) (la temperatura del aceite)
y u(x; t) (la temperatura del metal).
Bibliografı́a complementaria
[Haberman 1998] Haberman, R., Elementary Applied Partial Differential Equations, Pren-
tice Hall, Upper Saddle River, NJ, 1998.
57
www.elsolucionario.net
Capı́tulo 6
Una ecuación diferencial es una ecuación que involucra una función incógnita y al-
gunas de sus derivadas. Si la incógnita es una función de más de una variable, entonces
la ecuación se llama ecuación en derivadas parciales (EDP) pues las derivadas de la
www.elsolucionario.net
incógnita son derivadas parciales. Ejemplos de EDPs son las ecuaciones que aparecieron
en el capı́tulo anterior cuando vimos las leyes de conservación junto con las ecuaciones
constitutivas, y nuestro objetivo es resolver y comprender propiedades cuantitativas y
cualitativas de las soluciones de estas ecuaciones.
Cuando la función incógnita depende sólo de una variable, la ecuación se llama ecua-
ción diferencial ordinaria (EDO). A continuación veremos un breve repaso de algunos
aspectos básicos de EDOs, que serán necesarios para resolver las EDPs, que son la parte
principal de este curso.
Una EDO de primer orden se dice separable si puede escribirse de la siguiente manera:
dy
f (y) = g(x), (6.1)
dx
donde y es la función incógnita de la variable independiente x.
58
www.elsolucionario.net
dy
=5;
dx
dy 1 dy
= −3y (pues es equivalente a = −3);
dx y dx
dy 1 dy
= xy (pues es equivalente a = x);
dx y dx
dy 1 dy
− y 2 = 0 (pues es equivalente a 2 = 1);
dx y dx
dy
y = −x.
dx
www.elsolucionario.net
Z Z
dy
f (y) dx = f (y) dy = F (y) + C1 ,
dx
′
donde F (y) es una antiderivada, R o primitiva de f (y) (i.e., F (y) = f (y)) y C1 es una
constante arbitraria. El sı́mbolo f (y) dy representa la familia de todas las antiderivadas
de f (y) y como sabemos todas ellas difieren entre sı́ en una constante. En este caso, F (y)
es una cualquiera de ellas, y todas se obtienen sumándole una constante.
Si G(x) denota una antiderivada cualquiera de g(x), integrando el lado derecho de (6.1)
obtenemos
F (y) + C1 = G(x) + C2 .
Llamando C a C2 − C1 obtenemos
F (y) = G(x) + C,
y si pudiéramos despejar y de F (y) obtendrı́amos una fórmula para y(x). En otros casos,
la solución y queda implı́cita.
Veamos cómo resolver los ejemplos mencionados arriba:
dy
= 5. Integrando obtenemos que
dx
Z Z
dy
dx = 5 dx
dx
Z Z
dy = 5 dx
y + C1 = 5x + C2 .
59
www.elsolucionario.net
dy 1 dy
= −3y. Esta ecuación es equivalente a = −3. Integrando obtenemos que
dx y dx
Z Z
1 dy
dx = −3 dx
y dx
Z Z
1
dy = −3 dx
y
ln |y| + C1 = −3x + C2 ,
y(x) = C̃e−3x .
www.elsolucionario.net
• negativa para todo x: C̃ < 0;
• constantemente cero (para todo x): C̃ = 0.
dy 1 dy
= xy. Notemos que esta ecuación es equivalente a = x, y luego
dx y dx
Z Z
1 dy
dx = x dx
y dx
x2
ln |y| + C1 = + C2 ,
2
x2 x2
por lo que ln |y| = + C ó |y| = eC e 2 . La solución general es entonces
2
x2
y(x) = C̃e 2 .
60
www.elsolucionario.net
está desconectada. Por ejemplo, si queremos encontrar la solución que cumple y(0) =
1
4, entonces la constante debe ser C = − y por lo tanto la solución es
4
1 4
y(x) = − 1 = ,
x− 4
1 − 4x
1
que sólo está definida para x 6= . Pero la rama de la solución que cumple la
4 1
condición inicial está definida en − ∞, .
4
dy
y = −x. Integrando a ambos lados obtenemos
dx
Z Z
dy
y dx = −x dx
dx
Z Z
y dy = −x dx
www.elsolucionario.net
y2 x2
+ C1 = − + C2 ,
2 2
es decir que y 2 + x2 = C y la solución está implı́cita. En este caso podemos despejar
y obtener que √
y(x) = ± C − x2
√ √
que sólo existe en el intervalo − C ≤ x ≤ C y para C > 0.
Pero m(x) debe elegirse de manera que el lado izquierdo m(x)y ′ (x) + m(x)p(x)y(x) sea la
derivada de un sólo término, que será el producto de y(x) por otra función. Ahora bien,
′
m(x)y(x) = m(x)y ′ (x) + m′ (x)y(x)
61
www.elsolucionario.net
www.elsolucionario.net
Por lo tanto,
Z
m(x)y(x) = m(x)q(x) dx + C
Z
1
y(x) = m(x)q(x) dx + C .
m(x)
Lo importante no es recordar la fórmula final, sino el procedimiento:
Para resolver
y ′ (x) + p(x)y(x) = q(x)
hay que elegir el factor integrante m(x) de manera que al multiplicar
′ la ecuación por
m(x) resulte que el término de la izquierda sea igual a m(x)y(x) . Para ello, debe
ser
m′ (x) = m(x)p(x).
62
www.elsolucionario.net
Reemplazando lo que se sabe para VA (t) obtenemos la siguiente ecuación para VB (t):
www.elsolucionario.net
VB′ (t) + k VB (t) = k VA0 e−kt .
Multipliquemos por el factor integrante m(t) para ver qué debe satisfacer.
que se simplifica en ′
ekt VB (t) = k VA0 ,
es decir ekt VB (t) = k VA0 t + C, y entonces
63
www.elsolucionario.net
La cantidad de sal que ingresa al tanque B por minuto es la misma cantidad de sal
que drena desde el tanque A:
ℓ g g
20 × 5 = 100 .
min ℓ min
Por otro lado, la cantidad de sal que sale del tanque B por minuto es igual a la cantidad
de salmuera que sale por minuto, por la concentración, es decir
ℓ S(t) g S(t) g
15 × = 15 ,
min V (t) ℓ V (t) min
donde V (t) denota el volumen de salmuera en el tanque B, que fácilmente vemos es
ℓ
V (t) = 100 ℓ + (20 − 15) t min = 100 + 5t ℓ.
min
Luego, el problema a valores iniciales (ecuación diferencial más condición inicial) para
S(t) resulta:
15
S ′ (t) = 100 − S(t), S(0) = 0 .
{z100 + 5t }
| {z }
| condición inicial
ecuación diferencial
La ecuación diferencial es lineal de primer orden, que en forma estándar (y simplificando)
www.elsolucionario.net
se lee
3
S ′ (t) + S(t) = 100.
20 + t
3
El factor integrante debe cumplir m′ (t) = m(t), de manera que
20 + t
Z
3
ln |m(t)| = dt = 3 ln(20 + t) = ln(20 + t)3 ,
20 + t
3
y entonces m(t) = (20+t) . Multiplicando la ecuación diferencial por este factor integrante
obtenemos
(20 + t)3 S ′ (t) + 3 (20 + t)2 S(t) = 100 (20 + t)3
h i′
que es lo mismo que (20 + t)3 S(t) = 100 (20 + t)3 . Luego
64
www.elsolucionario.net
www.elsolucionario.net
Además, dado x0 ∈ R y dos números y0 , y1 , existe una única solución que satisface
y(x0 ) = y0 , y ′ (x0 ) = y1 .
65
www.elsolucionario.net
−b
Caso 2: d = 0. En este caso, la ecuación cuadrática tiene una sola raı́z r = de
2a
multiplicidad doble, que nos provee sólo una solución y1 (x) = erx . Podemos verificar que
otra solución linealmente independiente a esta primera es y2 (x) = xerx . Luego, la solución
www.elsolucionario.net
general de (6.4) es
b
d = b2 − 4ac = 0 =⇒ y(x) = C1 erx + C2 x erx , con r = − .
2a
66
www.elsolucionario.net
www.elsolucionario.net
La solución general es entonces
y(x) = C1 cos(2x) + C2 sen(2x).
67
www.elsolucionario.net
1 3
y(t) = − x′ (t) + x(t), (6.5)
4 4
www.elsolucionario.net
y lo reemplazamos, obteniendo
que resulta
x′′ (t) − 2x′ (t) + x(t) = 0.
x(t) = C1 et + C2 tet .
x(t) = et + C2 tet .
De (6.5)
1 t 3
y(t) = − e + C2 et + C2 tet + et + C2 tet
4 4
1 1 3 1 3
= − − C2 + − C2 t + C2 t et
4 4 4 4 4
1 1 1
= − C2 + C2 t et .
2 4 2
68
www.elsolucionario.net
www.elsolucionario.net
Por ejemplo, si u = u(x, y, z) es una función de tres variables, entonces u ∈ C 2 si u,
ux , uy , uz , uxx , uxy , uxz , uyx , uyy , uyz , uzx , uzy , uzz son continuas. Vemos aquı́ que u ∈ C 2
implica u ∈ C 1 pues u ∈ C 1 cuando u, ux , uy , uz son continuas.
Vemos entonces que en general, u ∈ C k implica que u ∈ C k−1 .
Es también sabido que si u ∈ C 2 entonces uxy = uyx , uyz = uzy , etc. Decimos en
general que si u ∈ C k , las derivadas parciales de orden menor o igual a k conmutan.
Definición 6.10. Dado un entero positivo k, una ecuación en derivadas parciales (EDP)
de orden k es una ecuación que involucra una función incógnita u, tal que k es el mayor
orden de derivada parcial que aparece en la ecuación.
Ejemplo 6.11. Veamos un par de ejemplos:
La ecuación del calor
ut = kuxx u = u(x, t)
es una ecuación de orden dos.
La ecuación de transporte
El dominio donde está planteada la ecuación es D = {(x, t) : 0 < x < L, t > 0}. Este
dominio es un rectángulo “espacio-temporal” semi-infinito.
69
www.elsolucionario.net
Definición 6.12. Una función u es una solución de una EDP de orden k en el dominio
D si es una función C k en la región D y la EDP se satisface en todos los puntos interiores
a D.
www.elsolucionario.net
función f que no depende de la incógnita. Si f ≡ 0 la EDP lineal se dice
homogénea.
u ux + uy = 1
uxx + uy = u2
70
www.elsolucionario.net
u = c 1 u1 + c 2 u2
es solución de
L[u] = c1 f1 + c2 f2 .
En particular, tomando f1 = f2 = 0 obtenemos que si u1 , u2 son soluciones del
problema homogéneo,
L[u1 ] = 0 L[u2 ] = 0
www.elsolucionario.net
entonces u = c1 u1 + c2 u2 es solución de L[u] = 0.
6.5. Ejercicios
6.1. La ley de Torricelly dice que (bajo ciertas circunstancias ideales) un recipiente
cilı́ndrico perderá fluido por un agujero en la base, a una tasa proporcional a la raı́z
cuadrada de la altura de la superficie del fluido. Supongamos que un contenedor cilı́ndri-
co está inicialmente lleno a una altura de 25cm. Si le lleva un minuto perder tres cuartos
del total. ¿Cuánto le llevará perder el total del fluido?
6.2. Resolver las siguientes ecuaciones lineales de primer orden, sujetas a las condiciones
iniciales indicadas
(a) y ′ (x) + 2y(x) = ex , y(0) = 1, (b) x′ (t) + x(t) cos(t) = 0, x(π) = 100.
6.3. Hallar la solución general y(x) de las siguientes EDO homogéneas de segundo orden
(a) y ′′ = 0 (d) y ′′ + 5y ′ + 4y = 0
25
(b) y ′′ + 36y = 0 (e) y ′′ + 4y ′ + 4
y =0
(c) y ′′ − 36y = 0
6.4. Hallar la solución particular que cumple las condiciones dadas:
(a) y ′′ = 0, y(0) = 1, y ′ (0) = 0
71
www.elsolucionario.net
25
(e) y ′′ + 4y ′ + 4
y = 0, y(0) = 0, y ′ (0) = 1
6.5. Hallar la solución particular que cumple las condiciones (homogéneas) dadas:
6.6. Hallar todas las soluciones que cumplen las condiciones dadas:
www.elsolucionario.net
(c) y ′′ − 36y = 0, y(0) = 0, y(π) = 0
ex − e−x ex + e−x
senh(x) = cosh(x) =
2 2
(a) Graficar ambas funciones y determinar el valor de cosh(0), senh(0), cosh′ (0) y senh′ (0).
(b) Determinar todos los valores de x para los que cosh(x) = 0, y también aquellos para
los que senh(x) = 0.
(e) Mostrar que dado un par de números A, B, existe un par de números Ã, B̃ tal que
72
www.elsolucionario.net
(g) Mostrar que si H ∈ R está dado (fijo), y λ > 0 entonces toda solución de y ′′ − λy = 0
se puede escribir
√ √
y(x) = Ø cosh λ (x − H) + B̃
˜
senh λ (x − H)
6.9. Dado λ > 0 fijo, hallar todas las soluciones de y ′′ − λy = 0 que cumplan
Ayuda: usar el ejercicio anterior para hallar la fórmula más sencilla posible.
6.10. Mostrar que para los siguientes operadores L se cumple que
www.elsolucionario.net
(b) L[u] = uxx + uyy + uzz , u = u(x, y, z)
Bibliografı́a complementaria
[Bleecker-Csordas] Bleecker, D., Csordas, G., Basic Partial Differential Equations, Inter-
national Press, Cambridge, Massachusetts, 1996.
[Edwards-Penney] Edwards, C.H. Jr., David E. Penney, Ecuaciones diferenciales elemen-
tales.
73
www.elsolucionario.net
Capı́tulo 7
En este capı́tulo nos dedicaremos a estudiar EDPs de primer orden. Veremos cómo se
resuelven las lineales, y haremos algunas consideraciones sobre la existencia y unicidad
de soluciones.
www.elsolucionario.net
7.1. EDPs de primer orden con coeficientes constan-
tes
Tal vez la EDP más simple es la lineal de primer orden con coeficientes constantes:
a ux + b uy + c u = f (x, y), u = u(x, y), (7.1)
donde a, b, c son constantes dadas tales que a2 + b2 6= 01 y f (x, y) es una función dada.
Como es usual, no se escriben las variables independientes de la función incógnita u. La
expresión u = u(x, y) a la derecha indica que la función u depende de dos variables x e y
y no más.
Como la ecuación es de primer orden, buscamos soluciones u = u(x, y) que sean C 1 ,
y para que esto sea posible, f (x, y) debe ser continua.
74
www.elsolucionario.net
www.elsolucionario.net
Dividimos por 3 y obtenemos la ecuación uy − u = . Multiplicamos por el factor
2
3 3
integrante m(y) = e− 3 y y obtenemos
∂ −2y 2 x + y
e 3 u(x, y) = e− 3 y .
∂y 3
Integrando con respecto a y obtenemos
− 32 y x 2 y 3 2
e u(x, y) = − e− 3 y − + e− 3 y + C(x),
2 2 4
2
y multiplicando por e 3 y obtenemos la solución general
x y 3 2 x+y 3 2
u(x, y) = − − + + C(x) e 3 y = − − + C(x) e 3 y .
2 2 4 2 4
Esto quiere decir que reemplazando C(x) por cualquier función C 1 de la variable x obte-
nemos una solución de (7.2). Ası́, todas las siguientes son soluciones:
u(x, y) C(x)
x + y 3 x2 2 y x2
− − + e3
2 4 8 8
x+y 3 2
− − + cos(πx) e 3 y cos(πx)
2 4
x+y 3 2 2
− − − 2e 3 (x+y) −2e 3 x
2 4
x+y 3 2
− − + e3y 1
2 4
75
www.elsolucionario.net
Está claro que estas no son todas las soluciones, pues hay infinitas. . .
Pasemos ahora al caso más general en que a 6= 0 6= b (esto quiere decir que a 6= 0 y
que b 6= 0, es decir, ninguna de las dos constantes es nula). La clave para resolver este
problema es ver que la parte de la ecuación (7.1) que involucra a las derivadas es:
a ux + b uy = ∇u · (ai + bj)
es decir, el producto escalar entre el gradiente de u y el vector (a, b). Por lo tanto, a ux +b uy
es una derivada direccional en la dirección ai + bj, y la ecuación es entonces una EDO a
lo largo de estas direcciones.
La idea es ahora introducir un cambio de variables de manera que una variable sea
constante a lo largo de las rectas paralelas al vector ai + bj. Estas rectas se llaman rectas
caracterı́sticas y son las de ecuación
b 1 2
y = x+d y =− y=0 y=1 y=
− a a a
− x − bx − a
a bx bx b
www.elsolucionario.net
o también bj
ai +
bx − ay = const.
Usaremos w, z para las nuevas variables. Elegimos w de manera que las caracterı́sticas
sean las rectas w = constante. Elegimos la otra variable z de alguna manera sencilla, pero
que no sea un múltiplo de w. Por ejemplo
x = w + az
(
w = bx − ay
⇐⇒ b
z=y y = z
Definimos entonces una nueva función v = v(w, z) que coincida con u = u(x, y) en
cada punto. Es decir
w + az
v(w, z) = u(x, y) = u( , z) o bien u(x, y) = v(w, z) = v(bx − ay, y)
b
Veamos qué EDP satisface v(w, z):
ux = v w w x + v z z x uy = v w w y + v z z y
ux = v w b + v z 0 y uy = vw (−a) + vz 1
a ux = a b v w b uy = −a b vw + b vz .
a ux + b uy = bvz
76
www.elsolucionario.net
bvz + cv = g(w, z)
www.elsolucionario.net
x = w − 3z
(
w = 2x + 3y
⇐⇒ 2
z=y
y=z
Identificando u(x, y) = v(w, z), veamos qué EDP satisface v(w, z):
ux = v w w x + v z z x uy = v w w y + v z z y
ux = v w 2 + v z 0 y uy = v w 3 + v z 1
3 ux = 6 v w 2uy = 6vw + 2 vz .
w − 3z
Por lo tanto 3ux −2uy = −2vz . Ası́, la ecuación (7.3) se transforma en −2vz +v = ,
2
que es equivalente a
1 w − 3z
vz − v = − .
2 4
z
El factor integrante es m(z) = e− 2 que multiplicándolo nos da
∂ z w − 3z − z
e− 2 v = − e 2.
∂z 4
Integrando obtenemos
z w −z 3 z
e− 2 v = e 2 − (2 + z)e− 2 + C(w),
2 2
z
y finalmente, multiplicando por e 2 , resulta
w 3 z
v= − (2 + z) + C(w)e 2 .
2 2
Reemplazando w, z por sus equivalentes en términos de x, y obtenemos
2x + 3y 3 y
u(x, y) = − (2 + y) + C(2x + 3y)e 2
2 2
77
www.elsolucionario.net
Es importante mencionar aquı́ que C(2x+3y) indica cualquier función C 1 de una variable,
aplicada a la expresión (2x + 3y) como un bloque indivisible. Por ejemplo
C(w) u(x, y)
y
w3 x − 3 + (2x + 3y)3 e 2
y
2 + sen(w) x − 3 + (2 + sen(2x + 3y)) e 2
y
1 e2
x−3+
1 + w2 1 + (2x + 3y)2
y
log(10 + w4 ) x − 3 + log(10 + (2x + 3y)4 )e 2
www.elsolucionario.net
Observación 7.3. Vemos que el problema tiene infinitas soluciones. Si tomamos C ≡ 0
tenemos la solución particular up (x, y) = x − 3, que es una solución de 3ux − 2uy + u = x.
y
La otra parte de la solución general, uh (x, y) = C(2x + 3y)e 2 es la familia de todas las
soluciones del problema homogéneo 3ux − 2uy + u = 0.
Ejemplo 7.4. ¿Qué habrı́a ocurrido en el ejemplo anterior si en lugar de definir z = y
hubiéramos intentado z = x? Veamos. Hacemos el cambio de variables
(
w = 2x + 3y x = z
⇐⇒
z=x y = w − 2z .
3
Identificando u(x, y) = v(w, z), veamos qué EDP satisface v(w, z):
ux = v w w x + v z z x uy = v w w y + v z z y
ux = v w 2 + v z 1 y uy = v w 3 + v z 0
3 ux = 6 vw + 3vz 2uy = 6vw .
Por lo tanto 3ux − 2uy = 3vz . Ası́, la ecuación (7.3) se transforma en 3vz + v = z, que es
equivalente a
1 z
vz + v = .
3 3
z
Multiplicamos por el factor integrante m(z) = e 3 y obtenemos
∂ z z z
e3 v = e3 .
∂z 3
z z
Integrando con respecto a z nos queda e 3 v = (z − 3)e 3 + G(w), o lo que es lo mismo
z
v(z, w) = z − 3 + G(w)e− 3 .
78
www.elsolucionario.net
Conclusiones:
Uno puede tener soluciones generales aparentemente diferentes pero que en realidad
son la misma familia de soluciones.
www.elsolucionario.net
7.1.1. Condición lateral
En muchas ocasiones, uno está interesado en encontrar soluciones que cumplen alguna
condición lateral. Para la ecuación (7.1) una condición apropiada puede ser que u(x, y)
tome ciertos valores especı́ficos en los puntos que están sobre una lı́nea recta o curva.
Ejemplos de tales condiciones son:
79
www.elsolucionario.net
Identificando u(x, y) = v(w, z), veamos qué EDP satisface v(w, z):
ux = v w w x + v z z x uy = v w w y + v z z y
ux = v w 1 + v z 0 y uy = v w 1 + v z 1
ux = v w uy = v w + v z .
www.elsolucionario.net
Una vez hallada la solución general imponemos la condición lateral u(x, 0) = x2 para
todo x ∈ R. Esto implica
1
x2 = u(x, 0) = + C(x), x ∈ R,
2
es decir C(x) = x2 − 21 . Luego, la solución es
1 2 1 2y
u(x, y) = − + (x + y) − e .
2 2
| {z }
C(x+y)
Nuevamente, el primer paso consiste en hallar la solución general. Las caracterı́sticas son
las rectas de pendiente 2, es decir las rectas de la forma y − 2x = const. Viendo que en el
lado derecho de la ecuación aparece ex+y elegimos las variables auxiliares de la siguiente
manera:
( z−w
w = y − 2x x =
⇐⇒ 3
z =x+y
y = 2z +w
3
Identificando u(x, y) = v(w, z), veamos qué EDP satisface v(w, z):
ux = v w w x + v z z x uy = v w w y + v z z y
ux = vw (−2) + vz 1 y uy = v w 1 + v z 1
ux = −2 vw + vz uy = v w + v z .
80
www.elsolucionario.net
www.elsolucionario.net
5x+2
C(2x + 2) = e− 3 . (7.6)
¿Cómo hacemos esto? Definimos r = 2x + 2 (lo que está como argumento de C), y
despejamos x, es decir x = r−2
2
. Luego reemplazamos esta expresión en todos los lugares
en que aparece x en la igualdad (7.6), es decir:
Hasta ahora no nos hemos encontrado con grandes dificultades, veamos un par de
ejemplos más.
Ejemplo 7.7. Hallar la solución de
81
www.elsolucionario.net
es decir
3x−1
C(−1) = e− 3 , x ∈ R.
Aquı́ vemos que estamos pidiendo a la función C(r) que en r = −1 tome diferentes valores
dependiendo de x. Esto es imposible, por lo que el problema plantedo no tiene solución.
La dificultad radica en que la condición lateral está dada sobre una caracterı́stica, la recta
y = 2x − 1 que tiene pendiente 2.
Nos preguntamos por qué tenemos dificultades cuando la condición lateral está dada
sobre una caracterı́stica. La explicación es la siguiente. La EDP resulta una EDO a lo largo
de las caracterı́sticas, y por eso el valor de la solución a lo largo de ellas está determinado
por el valor en un solo punto. Como la condición lateral no cumple la EDO sobre la
caracterı́stica, el problema no tiene solución.
Veamos un último ejemplo:
Ejemplo 7.8. Hallar la solución de
www.elsolucionario.net
Vemos que la ecuación diferencial es la misma, y como ya no somos tan inocentes también
detectamos que la condición lateral está dada sobre una caracterı́stica, la recta y = 2x.
Recordamos que la solución general es
4
u(x, y) = −ex+y + C(y − 2x)e 3 (x+y) ,
es decir
C(0) = 1, x ∈ R.
Aquı́ vemos que estamos pidiendo a la función C(r) que en r = 0 tome el valor 1. Esto
sı́ es posible, por lo que el problema plantedo tiene solución. Lo extraño aquı́ es que
existen infinitas funciones C(r) que valen 1 en r = 0. Por este motivo, este problema
tiene infinitas soluciones. He aquı́ algunos ejemplos:
C(r) u(x, y)
4
(r + 1)3 −ex+y + (y − 2x + 1)3 e 3 (x+y)
4
cos(r) −ex+y + cos(y − 2x) e 3 (x+y)
4 4 4
e− 3 r −ex+y + e− 3 (y−2x) e 3 (x+y) = −ex+y + e4x
4
ln(e + r4 ) −ex+y + ln e + (y − 2x)4 ) e 3 (x+y)
82
www.elsolucionario.net
Lo que ocurre en el ejemplo 7.8 es que la condición lateral está dada sobre una ca-
racterı́stica, pero la misma es compatible con la EDO que se cumple sobre esa recta. Por
lo tanto el problema tiene solución. Ahora uno se pregunta por qué tiene infinitas. La
existencia de infinitas soluciones se debe a que la condición lateral no impone ninguna
exigencia sobre ningún punto de las demás caracterı́sticas.
En base a los ejemplos analizados precedentemente concluimos algo que se puede
demostrar rigurosamente en general, pero cuya demostración queda fuera del alcance de
este curso:
Al resolver EDP lineales de primer orden con coeficientes constantes como (7.1):
Si la condición lateral está dada sobre una caracterı́stica, ocurre una de las
siguientes posiblidades:
www.elsolucionario.net
• El problema tiene infinitas soluciones.
Como mencionamos al principio de esta Sección una condición lateral podrı́a estar
dada sobre una curva, y no necesariamente sobre una lı́nea recta.
Definición 7.9. Si para la solución u = u(x, y) de (7.1) se requiere que tome valores
especı́ficos sobre una curva, tal curva se llama curva lateral.
Se dice que una curva es regular si tiene vector tangente unitario que varı́a de manera
continua con respecto a la longitud de arco. Visualmente es una curva que no tiene esqui-
nas. Una forma de saber si una curva es regular consiste en hallar una parametrización
t → x(t)i + y(t)j con x′ (t), y ′ (t) continuas y tales que (x′ (t))2 + (y ′ (t))2 > 0 para todo t.
Definición 7.10. Se dice que una curva regular intersecta a una lı́nea recta transversal-
mente si en cada punto intersección la recta no es tangente a la curva.
Teorema 7.11. Consideremos la EDP a ux + b uy + c u = f (x, y), u = u(x, y), y supon-
gamos que la condición lateral se impone sobre una curva lateral regular que corta cada
recta caracterı́stica exactamente una vez y lo hace transversalmente. Entonces, si los va-
lores impuestos sobre la curva varı́an de manera suave (C 1 ) el problema tiene solución
única.
Este teorema nos dice que para garantizar que el problema tenga solución única, la
condición lateral debe ser impuesta sobre una curva que satisfaga lo siguiente:
Que corte a todas las caracterı́sticas.
Que corte a cada caracterı́stica una sola vez.
Que en cada punto de la curva la recta tangente no sea caracterı́stica.
En la Figura 7.1 se muestran curvas laterales buenas para la imposición de valores
a la solución (izquierda), y también se muestran curvas laterales que no garantizan la
existencia de solución única (derecha).
83
www.elsolucionario.net
Figura 7.1: Curvas laterales adecuadas e inadecuadas. Las curvas de la izquierda son
adecuadas porque cortan todas las rectas caracterı́sticas, cortan a cada caracterı́stica una
sola vez y lo hacen de manera transversal. Las curvas de la derecha no cumplen estas tres
condiciones a la vez.
www.elsolucionario.net
la edad, a medida que transcurre el tiempo. Los individuos pueden ser tanto organismos
biológicos como objetos manufacturados, por ejemplo lámparas, transistores, productos
alimenticios o medicamentos. Nos referiremos entonces a alguna colección de objetos simi-
lares que se vuelven defectuosos con la edad, o se vencen, acorde a un patrón establecido
estadı́sticamente. Utilizaremos el verbo expirar para decir que un individuo deja de exis-
tir, o de funcionar, o se vence. Y utilizaremos el verbo existir para decir que un individuo
está vivo, o funciona correctamente, o no está vencido.
Supongamos entonces que la variable P (x, t) denota una densidad de población a
tiempo t de individuos de edad x. Es decir,
Z b
P (x, t) dx = cantidad de individuos de edad entre a y b que existen a tiempo t.
a
Si los individuos no expiraran (D ≡ 0), tendrı́amos que P (x, t + ∆t) = P (x − ∆t, t).
Es decir, la densidad de individuos de edad x a tiempo t + ∆t es la misma que la densidad
de individuos de edad x − ∆t a tiempo t.
84
www.elsolucionario.net
www.elsolucionario.net
Z r
P (x, t + r) = P (x − r, t) − D(x − r + s, t + s)P (y − r + s, t + s) ds.
0
85
www.elsolucionario.net
Qz + D(z, z − w)Q = 0,
donde
R hemos reemplazado x por z, y t por z − w. El factor integrante es ahora m(z) =
D(z,z−w) dz
e . Multiplicamos y obtenemos
∂ R D(z,z−w) dz
e Q = 0.
∂z
R
D(z,z−w) dz
Integrando con respecto a z nos queda e Q = C(w). Eligiendo una primitiva
en particular obtenemos
Z z
Q(w, z) = C(w) exp − D(ζ, ζ − w) dζ .
w
www.elsolucionario.net
Z x
P (x, t) = C(x − t) exp − D(ζ, ζ − x + t) dζ .
x−t
Tomando t = 0 obtenemos P (x, 0) = C(x), por lo que C(x) indica la densidad inicial de
población con respecto a la edad, al menos para x ≥ 0. Es decir
Z t
P (x, t) = P (x − t, 0) exp − D(x − t + τ, τ ) dτ .
0
Como P (x, 0) no está definido para x < 0, resulta que P (x, t) no está definido para x < t
(ver Figura 7.2). Es entonces conveniente definir, para x < 0, P (x, 0) como el número
de nacimientos por unidad de tiempo que se producirán a tiempo |x|, definiendo a la vez
D(x, t) = 0 para x < 0.
Consideremos el caso de producción constante C, es decir P (x, 0) = C para todo
x < 0 y tasa de expiración independiente del tiempo, o sea D(x, t) = D(x) (recordando
que D(x) = 0 si x < 0). Analicemos qué ocurre a medida que pasa mucho tiempo. Más
precisamente, calculemos la densidad de población en estado estacionario:
Z t
P∞ (x) = lı́m P (x, t) = lı́m P (x − t, 0) exp − D(x − t + τ ) dτ .
t→+∞ t→+∞ 0
lı́m P (x − t, 0) = C;
t→+∞
86
www.elsolucionario.net
t
D≡0 D>0
x<t
P es constante
sobre estas lı́neas
x>t
x=t
x
además
www.elsolucionario.net
Z t Z x Z 0 Z x Z x
D(x − t + τ ) dτ = D(s) ds = D(s) ds + D(s) ds = D(s) ds.
0 x−t x−t 0 0
| {z }
=0
Más aún, la densidad estacionaria P∞ (x) coincide con la densidad P (x, t) cuando t > x.
Más precisamente
P∞ (x) = P (x, t), si x > t.
Esto es cierto sólo cuando hay producción constante.
87
www.elsolucionario.net
Definición 7.13. Una curva en el plano xy se llama curva caracterı́stica para la EDP (7.7)
si en cada punto (x0 , y0 ) de la curva el vector g(x0 , y0 ) = a(x0 , y0 )i + b(x0 , y0 )j es tangente
a la curva.
Si la curva caracterı́stica es el gráfico de una función y(x), entonces la definición nos
dice que
www.elsolucionario.net
dy b(x, y)
= .
dx
|{z} a(x, y)
| {z }
pendiente de la pendiente del
recta tangente vector g(x, y)
La ecuación
dy b(x, y)
= (7.8)
dx a(x, y)
se llama ecuación caracterı́stica de (7.7) y los gráficos de las soluciones son las curvas
caracterı́sticas de (7.7).
La EDP (7.7) resulta una EDO sobre las curvas caracterı́sticas.
En el caso a, b constante, las caracterı́sticas son las rectas y = ab x + d o también
bx − ay = const. En ese caso definı́amos la variable w como una variable auxiliar que nos
ayudarı́a a resolver el problema. Ahora haremos algo parecido.
Supongamos que ya pudimos resolver (7.8) y que las curvas caracterı́sticas son de la
forma
h(x, y) = const.
Por ejemplo h(x, y) = y − x2 o h(x, y) = x2 + y 2 (ver figura 7.4)
Entonces hacemos el cambio de variables y la identificación
(
w = h(x, y)
v(w, z) = u(x, y)
z=y
La idea es que w será constante en cada caracterı́stica y entonces la EDP se transfor-
mará en una EDO de la variable z. Veamos si esto ocurre
ux = v w w x + v z z x uy = v w w y + v z z y
ux = v w h x + v z 0 y uy = v w h y + v z 1
a ux = ahx vw b uy = bhy vw + b vz .
88
www.elsolucionario.net
www.elsolucionario.net
Por lo tanto
a ux + b uy = (a hx + b hy ) vw + b vz .
Observemos ahora que a hx + b hy = (ai + bj) · ∇h, y ai + bj es tangente a las curvas de
nivel de h, por lo tanto ai + bj es perpendicular a ∇h. Es decir
a hx + b hy = (ai + bj) · ∇h = 0,
−y ux + x uy = 0, u = u(x, y).
El primer paso consiste, como antes, en encontrar las curvas caracterı́sticas. Las mismas
deben satisfacer
dy x
=− ,
dx y
que es una EDO de primer orden separable. Multiplicamos por y e integramos a ambos
lados, obteniendo
dy
y = −x
Z dx Z
y dy = − x dx
y2 x2
=− +C
2 2
2 2
x + y = const.
89
www.elsolucionario.net
Identificando u(x, y) = v(w, z), veamos qué EDP satisface v(w, z):
ux = v w w x + v z z x uy = v w w y + v z z y
ux = vw 2x + vz 0 y uy = vw 2y + vz 1
−yux = −2xy vw x uy = 2xy vw + xvz .
Por lo tanto −yux + xuy = xvz , y la ecuación resulta xvz = 0. Por lo tanto vz debe ser
cero en todos los puntos donde x 6= 0, pero como buscamos soluciones C 1 , vz debe ser
continua y por lo tanto vz ≡ 0. Es decir
v(w, z) = C(w),
www.elsolucionario.net
u(x, y) = C(x2 + y 2 ).
En otras palabras, las soluciones son todas las funciones radiales, o sea, las que son
constantes a lo largo de las circunferencias centradas en el origen. Por ejemplo
C(w) u(x, y)
w3 (x2 + y 2 )3
2 +y 2 )
e−w e−(x
1 1
1+w 1 + x2 + y 2
90
www.elsolucionario.net
Si x(t), y(t) son soluciones del sistema caracterı́stico (7.9) entonces la parametrización
x(t)i + y(t)j se llama parametrización preferida de la caracterı́stica.
www.elsolucionario.net
que nunca.
Resolvamos nuevamente el ejemplo anterior utilizando la parametrización preferida:
Ejemplo 7.15. Utilizando las parametrizaciones preferidas hallar la solución general de
−y ux + x uy = 0, u = u(x, y).
Hallamos ahora las parametrizaciones preferidas de las caracterı́sticas. Para ello plantea-
mos el sistema caracterı́stico correspondiente:
x′ (t) = −y(t), y ′ (t) = x(t),
que es como los estudiados en la Sección 6.3. Para resolverlo derivamos la primera ecuación
y reemplazamos y ′ (t) por x(t) obteniendo
x′′ (t) = −x(t),
con solución general
x(t) = C1 cos t + C2 sen t, y(t) = −x′ (t) = C1 sen t − C2 cos t.
Observamos que x(t)2 + y(t)2 = C12 + C22 , ası́ que las caracterı́sticas son cı́rculos centrados
en el origen, y la misma familia de curvas se obtiene tomando C2 = 0 y C1 constante
arbitraria; pero con una fórmula más simple. Ahora, si u(x, y) es solución de la ecuación
planteada, U (t) = u(x(t), y(t)) debe satisfacer
U ′ (t) = ux x′ + uy y ′ = ux (−y) + uy x = 0,
por lo que U (t) debe ser constante, o lo que es lo mismo, u(x, y) es constante sobre cı́rculos
centrados en el origen. En consecuencia, la solución general es:
u(x, y) = f (x2 + y 2 ), f ∈ C 1.
Observación 7.16. Es interesante notar que el problema anterior no tiene solución
cuando imponemos la condición lateral u(x, 0) = 3x, x ∈ R. Esto ocurre porque el eje
x corta cada caracterı́stica dos veces. Si impusiéramos
p u(x, 0) = 3x, x > 0 entonces el
problema sı́ tendrı́a solución, y única: u(x, y) = 3 x + y 2 .
2
91
www.elsolucionario.net
X(s, t) = s cos t
s > 0, 0 ≤ t < 2π.
Y (s, t) = s sen t
De manera que para cada s0 > 0 fijo la curva X(s0 , t)i + Y (s0 , t)j = s0 cos t i + s0 sen t j,
www.elsolucionario.net
0 ≤ t < 2π describe la circunferencia de radio x0 centrada en el origen. Y para t = 0,
X(s, 0)i + Y (s, 0)j = s cos 0i + s sen 0j = si describe la curva lateral, que en este caso es
el eje x > 0.
Consideremos entonces, en general, la EDP
a(x, y)ux + b(x, y)uy + c(x, y)u = f (x, y), u = u(x, y),
curvas caracterı́sticas
t variable, s fijo
Buscamos entonces una familia de parametrizaciones preferidas X(s, t)i + Y (s, t)j de
las curvas caracterı́sticas, de manera que para t = 0 estemos sobre la curva lateral (ver
92
www.elsolucionario.net
www.elsolucionario.net
Ut + c(X(s, t), Y (s, t))U = f (X(s, t), Y (s, t)).
Si definimos C(s, t) = c(X(s, t), Y (s, t)) y F (s, t) = f (X(s, t), Y (s, t)), U (s, t) debe cum-
plir
Ut (s, t) + C(s, t)U (s, t) = F (s, t), U (s, 0) = g(s), s ∈ I.
Esta última ecuación es entonces una familia de ecuaciones diferenciales ordinarias lineales
de primer orden, una para cada s ∈ I, que sabemos resolver. Veamos un ejemplo concreto:
Ejemplo 7.18. Hallar la solución en forma paramétrica de
Xtt = Xt + Yt = Xt + |{z}
Y −X = Xt + Xt − X − X = 2Xt − 2X,
|{z}
Y −X Xt −X
por lo que
Xtt − 2Xt + 2X = 0.
Proponiendo X = ert obtenemos r2 − 2r + 2 = 0 que tiene soluciones complejas r = 1 ± i.
La solución es entonces X(s, t) = C1 (s)et cos(t) + C2 (s)et sen t; notar que las constantes
pueden depender de la variable s. Las curvas caracterı́sticas están entonces dadas por:
X(s, t) = et C1 (s) cos t + C2 (s) sen t
Y (s, t) = et − C1 (s) sen t + C2 (s) cos t
93
www.elsolucionario.net
Imponiendo las condiciones iniciales X(s, 0) = cos s, Y (s, 0) = sen s tenemos que C1 (s) =
cos s y C2 (s) = sen(s), por lo que
X(s, t) = et cos s cos t + sen s sen t = et cos(s − t)
Y (s, t) = et (− cos s sen t + sen s cos t = et sen(s − t)
Definiendo U (s, t) = u(X(s, t), Y (s, t)), con u(x, y) la solución buscada,
Ut = ux Xt + uy Yt = ux (X + Y ) + uy (Y − X) = u = U,
www.elsolucionario.net
Y (s, t) = et sen(s − t), t ∈ R, 0 ≤ s < 2π. (7.10)
U (s, t) = et ,
A veces, muy pocas veces, puede deducirse la fórmula de u(x, y) en base a la formu-
lación paramétrica. En el ejemplo anterior por ejemplo, es fácil ver que
p
et = U (s, t) = X(s, t)2 + Y (s, t)2
p
y entonces u(x, y) = x2 + y 2 .
Pero en general no es tan fácil hallar la fórmula de u(x, y) aunque la forma paramétri-
ca (7.10) es igualmente útil. Veamos otro ejemplo con la misma ecuación, pero diferente
condición lateral, dada sobre la misma curva lateral.
(y +x)ux +(y −x)uy −u = 0, u = u(x, y), u(cos(s), sen(s)) = cos(2s), 0 ≤ s < 2π.
Las parametrizaciones preferidas son las mismas, y U (t, s) = C(s)et , pero ahora la condi-
ción lateral nos dice que U (0, s) = cos(2s), ası́ que C(s) = cos(2s) y la solución en forma
paramétrica es:
X(s, t) = et cos(s − t),
Y (s, t) = et sen(s − t), t ∈ R, 0 ≤ s < 2π. (7.11)
t
U (s, t) = cos(2s)e ,
Ya no es tan fácil despejar u en términos de x e y, pero sı́ es fácil por ejemplo, graficar
la solución usando un software como MATLAB u OCTAVE; ver Figura 7.6.
94
www.elsolucionario.net
s = [0:0.1:2*pi]; s = [0:0.1:2*pi];
t = [-1:0.05:1]; t = [-1:0.05:1];
www.elsolucionario.net
[s, t] = meshgrid(s, t); [s, t] = meshgrid(s, t);
mesh(x,y,u) mesh(x,y,u)
3 3
2.5 2
2 1
1.5 0
1 -1
0.5 -2
03 -3 3
2 2
3 3
1 1
2 2
0 1 0 1
-1 0 -1 0
-1 -1
-2 -2
-2 -2
-3 -3 -3 -3
95
www.elsolucionario.net
10 10
5 5
0 0
−5 −5
−10 −10
−10 −5 0 5 10 −10 −5 0 5 10
Figura 7.7: Curvas caracterı́sticas y curvas laterales de los problemas (7.12) y (7.13). En
el problema 7.12 la curva lateral determina la solución en el primer y el tercer cuadrante
(izquierda), mientras que en el problema 7.13 la curva lateral determina la solución en el
primer y el segundo cuadrante (derecha).
www.elsolucionario.net
7.2.3. Consideraciones globales
Cuando los valores de u se prescriben sobre una curva regular que intersecta cada
caracterı́stica una sola vez de manera transversal, entonces tenemos una manera de en-
contrar la solución u en todo el plano. Sin embargo, como veremos en el siguiente ejemplo,
puede ocurrir que no existe una curva que corte a todas las caracterı́sticas. En tales ca-
sos, no será posible hallar la solución en todo el plano, sino en algunas regiones, a menos
que impongamos condiciones laterales sobre más de una curva. Esto no es en general
un problema, dado que en muchas ocasiones la región de interés será la cubierta por las
caracterı́sticas que cortan a la curva lateral indicada. Veamos un ejemplo.
Consideremos el problema
xux − yuy + yu = 0, u = u(x, y), u(x, x) = 1. (7.12)
y el problema
xux − yuy + yu = 0, u = u(x, y), u(x, 2 + x2 /10) = ex . (7.13)
dy y
Para ambos problemas, las caracterı́sticas son las curvas que cumplen = −
dx x
1 dy 1
ó = − . Integrando obtenemos ln |y| = − ln |x| + C, o ln |y| + ln |x| = C, que es
y dx x
equivalente a
xy = const.
Las curvas caracterı́sticas son entonces las hipérbolas que se ven en la Figura 7.7. Al
intentar resolver cada uno de los problemas planteados vemos que la condición lateral no
podrá determinar la solución en todo el plano xy sino solamente en un par de cuadrantes;
ver Figura 7.7.
Observación 7.20. Es importante notar que esta particularidad de que ocurra que
una curva no pueda cortar a todas las caracterı́sticas, sólo puede suceder cuando los
96
www.elsolucionario.net
coeficientes a(x, y) y b(x, y) son variables. Cuando los coeficientes son constantes las
caracterı́sticas son todas rectas paralelas entre sı́.
a(x, y, z)ux + b(x, y, z)uy + c(x, y, z)uz + d(x, y, z)u = f (x, y, z), u = u(x, y, z),
para funciones dadas a, b, c, d y f . Las curvas caracterı́sticas x(t)i + y(t)j + z(t)k para-
metrizadas de la manera preferida son las soluciones del sistema
dx dy dz
= a(x(t), y(t), z(t)), = b(x(t), y(t), z(t)), = c(x(t), y(t), z(t)).
dt dt dt
En la práctica es usualmente más conveniente tratar a x como un parámetro en lugar
de t, en cuyo caso el sistema anterior se transforma en un sistema de dos ecuaciones (si
www.elsolucionario.net
suponemos que a(x, y, z) 6= 0)
dy b(x, y(x), z(x)) dz c(x, y(x), z(x))
= , = ,
dx a(x, y(x), z(x)) dx a(x, y(x), z(x))
para las incógnitas y = y(x) y z = z(x).
Veamos un ejemplo con coeficientes constantes.
Ejemplo 7.21. Hallar la solución general de
Identificando u(x, y, z) = ū(x̄, ȳ, z̄), veamos qué EDP satisface ū(x̄, ȳ, z̄):
ux = ūx̄ x̄x + ūȳ ȳx + ūz̄ z̄x uy = ūx̄ x̄y + ūȳ ȳy + ūz̄ z̄y
ux = ūx̄ (−3) + ūȳ (−5) + ūz̄ 0 uy = ūx̄ 2 + ūȳ 0 + ūz̄ 0
2 ux = −6 ūx̄ − 10 ūȳ 3uy = 6 ūx̄ .
97
www.elsolucionario.net
www.elsolucionario.net
2ux + 3uy + 5uz − u = 0, u = u(x, y, z), u(x, y, 0) = x2 sen y, x ∈ R, y ∈ R.
Ya sabemos que la solución general es u(x, y, z) = C(2y − 3x, 2z − 5x)ez/5 . Debemos im-
poner la condición lateral u(x, y, 0) = x2 sen y para determinar la función C. La condición
lateral implica
C(2y − 3x, −5x) = x2 sen y.
Tomemos r = 2y − 3x y s = −5x, entonces x = − 5s , y = 21 (r − 35 s), y entonces
s 2 s 2
1 3 r 3
C(r, s) = sen (r − s) = sen − s .
5 2 5 5 2 10
7.4. Ejercicios
7.1. Encontrar la solución general de cada una de las siguientes ecuaciones lineales de
primer orden con coeficientes constantes:
(a) 2 ux − 3 uy = x (b) ux + uy − u = 0.
(a) Hallar la solución general, y decir cuáles son las curvas caracterı́sticas.
98
www.elsolucionario.net
En cada uno de los siguientes casos hallar (si es posible) una solución de la ecuación
resuelta en (a) que cumpla la condición lateral indicada. Si no hay, indicar por qué. Si
hay más de una, indicar tres soluciones diferentes.
−x + 5
(b) u(−5x, x) = 3ex , (c) u(x, ) = x(1 − x), (d) u(x, x) = x2 .
5
7.3. ¿Qué forma debe tener la función g(x) para que el siguiente problema tenga solución?
(
u x + 3 uy − u = 1 (x, y) ∈ R2
u(x, 3 x) = g(x) x ∈ R.
Si g cumple tal condición. ¿Tiene el problema solución única?
* 7.4. Dar un ejemplo de ecuación lineal de primer orden con coeficientes constantes (en
R2 ), y condición lateral dada sobre una curva lateral que corte todas las caracterı́sticas,
pero que no tenga solución. Explicar por qué no tiene solución.
7.5. Encontrar la solución general, en los dominios indicados, de cada una de las siguientes
ecuaciones lineales de primer orden con coeficientes variables:
www.elsolucionario.net
(a) x ux + 2 y uy = 0, x > 0, y > 0.
(b) y ux − 4 x uy = 2 x y, (x, y) ∈ R2 .
7.6. Encontrar la forma paramétrica de la solución de cada uno de los siguientes proble-
mas:
( (
x ux + 2 y uy = 0, x > 0, y > 0 y ux − 4 x uy = 2 x y, (x, y) ∈ R2
(a) (b)
u(s, e−s ) = sen s, s > 0. u(s, s3 ) = 1, s∈R
Graficar (usando MATLAB) en una parte del dominio para ver cómo lucen las soluciones.
7.7. Resolver el problema
ux + uy + uz = u, u = u(x, y, z), u(x, y, 0) = x2 + y 2 .
7.8. Hallar la solución general de la EDP ux − uy + u = z, u = u(x, y, z).
(a) Hallar la solución que cumple con la condición lateral u(0, y, z) = y 2 ez .
(b) Mostrar que no hay ninguna solución que satisfaga u(x, y, x + y) = 0.
(c) Hallar tres soluciones que satisfagan u(x, y, x + y) = x + y + ey .
(d) Explicar lo ocurrido en (b) y (c) relacionando las rectas caracterı́sticas de la EDP
con el plano z = x + y donde se dan las condiciones laterales.
(e) ¿Qué relación existe entre el plano x = 0 donde se da la condición lateral del item
(a) y las rectas caracterı́sticas?
Bibliografı́a complementaria
[Bleecker-Csordas 1996] Bleecker, D., Csordas, G., Basic Partial Differential Equations,
International Press, Cambridge, Massachusetts, 1996.
99
www.elsolucionario.net
Capı́tulo 8
La ecuación de difusión
unidimensional
Vimos en la Sección 5.3 que la ecuación de difusión unidimensional cuando los coefi-
cientes son constantes resulta
www.elsolucionario.net
ut − kuxx = f, u = u(x, t), (8.1)
Φ(x, t) no está definida para t = 0, pero vemos que Φ(x, t) se aproxima a la Delta
de Dirac cuando t → 0+ .
Z ∞
Φ(x, t) dx = 1, para cualquier t > 0.
−∞
100
www.elsolucionario.net
Si g es continua y acotada, entonces se puede ver que las derivadas con respecto a
x y a t se pueden meter dentro de la integral y
Z ∞
ux (x, t) = g(y)Φx (x − y, t) dy,
Z−∞
∞
uxx (x, t) = g(y)Φxx (x − y, t) dy,
−∞
Z ∞
ut (x, t) = g(y) Φt (x − y, t) dy.
−∞ | {z }
=kΦxx (x,t)
www.elsolucionario.net
Z ∞ Z ∞
u(x, t) = g(y)Φ(x − y, t) dy = g(x − y)Φ(y, t) dy
−∞ −∞
Z ∞
1 y2
=√ g(x − y)e− kt dy
4πkt −∞
Pero el lı́mite de esta última fórmula es g(x) cuando t → 0+ por lo visto en la
Sección 1.6 (ver también ejercicio 1.4).
Por lo tanto, si extendemos la definición de u(x, t) a x ∈ R, 0 ≤ t < ∞ de la
siguiente manera:
Z ∞
g(y)Φ(x − y, t) dy, si t > 0,
u(x, t) = −∞
g(x), si t = 0,
obtenemos una solución del problema
ut = kuxx , x ∈ R, t > 0
u(x, 0) = g(x), x ∈ R (t = 0).
La solución dada por esta fórmula para el dato inicial
(
1, si − 1 ≤ x ≤ 1
g(x) = (8.4)
0, en c.o.c.
Como las soluciones del problema de difusión homogéneo se pueden obtener haciendo
convolución del dato inicial con Φ(x, t) de la fórmula (8.2), esta función Φ(x, t) se
denomina solución fundamental de la ecuación de difusión en 1d.
101
www.elsolucionario.net
2.5
t = 1/4 t=0
t=1 1 t = 1/4
t=4 t=1
2
t = 16 t=4
0.8 t = 16
1.5
0.6
1
0.4
0.5 0.2
0 0
−3 −2 −1 0 1 2 3 −3 −2 −1 0 1 2 3
www.elsolucionario.net
8.2. Solución con CB Dirichlet homogéneas
Consideramos el siguiente problema
∂u ∂2u
EDP: = k 2, 0 < x < L, t > 0,
∂t
( ∂x
u(0, t) = 0, (8.5)
CB: t > 0,
u(L, t) = 0,
CI: u(x, 0) = f (x), 0 < x < L.
Buscamos todas las soluciones posibles que sean producto de una función de x y una
de t y que satisfagan la ecuación diferencial y las condiciones de borde homogéneas.
102
www.elsolucionario.net
∂u dG(t) ∂ 2u d2 φ(x)
(x, t) = φ(x) , (x, t) = G(t).
∂t dt ∂x2 dx2
En consecuencia la ecuación de difusión (8.5) implica que
dG(t) d2 φ(x)
φ(x) =k G(t). (8.7)
dt dx2
www.elsolucionario.net
Ahora nos damos cuenta de que podemos “separar variables” dividiendo ambos miembros
de (8.7) por el producto φ(x)G(t):
1 dG(t) 1 d2 φ(x)
=k .
G(t) dt φ(x) dx2
Se ve que las variables han sido “separadas” en el sentido que el lado izquierdo es sólo
función de la variable t y el lado derecho es sólo función de la variable x. Podemos
continuar trabajando con esta expresión, pero es conveniente dividir también por k, y
luego
1 dG 1 d2 φ
= 2
. (8.8)
kG dt
| {z } φ dx
| {z }
función de t función de x
¿Cómo es posible que una función solamente dependiente de t sea igual a una función que
depende solamente de x? La única posibilidad es que ambas expresiones sean constantes
e independientes tanto de x como de t. Y por la igualdad deben ser iguales a la misma
constante:
1 dG 1 d2 φ
= = −λ, (8.9)
kG dt φ dx2
donde λ es una constante arbitraria conocida como constante de separación. En un
momento explicaremos la aparición del misterioso signo menos, que fue introducido sólo
por conveniencia.
La ecuación (8.9) equivale a dos ecuaciones diferenciales ordinarias, una para G(t) y
otra para φ(x):
d2 φ
= −λφ (8.10)
dx2
dG
= −λkG. (8.11)
dt
103
www.elsolucionario.net
Tenemos ahora tres incógnitas: φ, G, y λ. Para hallar las soluciones producto (8.6)
debemos resolver las dos ecuaciones diferenciales ordinarias que acabamos de deducir, y
ver para qué constantes λ estas soluciones no dan la solucion trivial.
Repetimos que λ es constante, y es la misma en ambas ecuaciones (8.10) y (8.9). Las
soluciones producto también deben satisfacer las dos condiciones de borde. Por ejemplo,
u(0, t) = 0 implica que φ(0)G(t) = 0. Hay aquı́ dos posibilidades. O bien G(t) ≡ 0, o
φ(0) = 0. Si G(t) ≡ 0, entonces u ≡ 0 y obtuvimos la solución trivial. Como estamos
buscando soluciones no-triviales, descartamos este caso y suponemos G(t) 6≡ 0. Si G(t) 6≡
0, es necesario que
φ(0) = 0. (8.12)
φ(L) = 0. (8.13)
www.elsolucionario.net
ecuaciones diferenciales ordinarias (8.10) y (8.11), deben satisfacer las condiciones de
borde (8.12) y (8.13).
d2 φ
= −λφ, φ(0) = 0, φ(L) = 0,
dx2
dG
= −λkG.
dt
Recordemos que λ es la constante de separación, que por el momento pareciera que puede
tomar cualquier valor. Veremos a continuación que λ sólo puede tomar ciertos valores.
104
www.elsolucionario.net
d2 φ
= −λφ
dx2
(8.15)
φ(0) = 0
φ(L) = 0.
www.elsolucionario.net
No hay ninguna teorı́a simple que garantice que la solución existe o es única para este tipo
de problemas. En particular, notamos que φ(x) ≡ 0 satisface la EDO y ambas condiciones
de borde homogéneas, sin importar el valor de λ, aún para λ < 0. Esta solución se conoce
con el nombre de solución trivial del problema a valores de borde. Como u(x, t) =
φ(x)G(t), corresponde a u(x, t) ≡ 0, si las soluciones de (8.15) fueran únicas, φ(x) ≡ 0
serı́a la única solución; y no podrı́amos obtener soluciones no-triviales de una EDP lineal
homogénea por el método de separación de variables. Afortunadamente sı́ existen otras
soluciones de (8.15), aunque no para cualquier valor de λ sino sólo para algunos. Veremos
a continuación que hay ciertos valores de λ, llamados autovalores de (8.15), para los
cuales existen soluciones no-triviales φ(x). Una φ(x) no-trivial, correspondiente a algún
valor de λ se llama autofunción correspondiente al autovalor λ.
105
www.elsolucionario.net
www.elsolucionario.net
Si c2 = 0, obtenemos la solución trivial φ(x) ≡ 0 que no nos interesa. Si c2 6= 0 entonces
tenemos que buscar todos los valores de λ para los que
√
sen λL = 0.
√
En√ otras palabras, λL debe ser un cero de la función seno, es decir un múltiplo entero de
π: λL = nπ. Como hemos supuesto λ > 0, n toma los valores 1, 2, . . . . Los autovalores
λ son entonces nπ 2
λ= , n = 1, 2, 3 . . . . (8.18)
L
2
La autofunción correspondiente al autovalor λ = nπ L
es
√ nπx
φ(x) = c2 sen λx = c2 sen (8.19)
L
donde c2 es una constante arbitraria. A menudo elegimos una autofunción seleccionando
un valor no-nulo para c2 , por ejemplo c2 = 1. Debemos observar que cualquier autofunción
puede ser multiplicada por una constante arbitraria y seguir siendo autofunción del
mismo autovalor pues la ecuación es lineal y homogénea.
106
www.elsolucionario.net
www.elsolucionario.net
similarmente a lo que ocurre con las trigonométricas, pero más fáciles de recordar porque
no cambia el signo al derivar. Otras propiedades importantes son las siguientes:
cosh z 6= 0 para todo z;
senh z = 0 sólo si z = 0.
Como ambas funciones hiperbólicas son sumas de exponenciales y son independientes,
la solución general también se puede expresar como
√ √
φ(x) = c1 cosh −λx + c2 senh −λx, (8.21)
una forma equivalente a (8.20). Aplicamos ahora las condiciones de borde a la solución
√
general (8.21). La condición φ(0) = 0 √implica c1 = 0 y luego φ(x) = c2 senh −λx.
La condición φ(L) = 0 implica c2 senh −λL = 0. Pero senh z = 0 sólo para z = 0 y
entonces, como λ 6= 0 y L > 0, obtenemos c2 = 0 quedándonos nuevamente la solución
trivial. Por lo tanto no hay autovalores negativos para este problema.
107
www.elsolucionario.net
∂u ∂ 2u
= k 2, u(0, t) = 0, u(L, t) = 0,
∂t ∂x
entonces
√ nπ 2
φ(x) = c2 sen λx, G(t) = ce−λkt , y λ= , n = 1, 2, 3, . . . .
L
Escrito de otra forma,
nπx nπ 2
φ(x) = c2 sen , G(t) = ce−( L ) kt
, n = 1, 2, 3, . . . .
L
Por lo tanto, todas las soluciones producto de la ecuación de difusión (8.5) que cumplen
las CB son
www.elsolucionario.net
nπx −( nπ 2
u(x, t) = b sen e L ) kt , n = 1, 2, 3, . . .
L
donde b es una constante arbitraria. Es decir, para cada número real b y cada entero
positivo n, la expresión anterior es una solución de la ecuación de difusión con condiciones
de borde de tipo Dirichlet homogéneas.
Ejemplo 8.1. Las siguientes funciones son soluciones de la ecuación de difusión con
condiciones de borde de u prescripta igual a cero:
3πx − 9π22 kt
u1 (x, t) = sen e L
L
πx − π22 kt
u2 (x, t) = 4,5 sen e L .
L
Y como la ecuación y las condiciones de borde que estamos considerando son lineales y
homogéneas, también
3πx − 9π22 kt πx − π22 kt
u3 (x, t) = sen e L + 4,5 sen e L
L L
es una solución.
¿Qué condición inicial satisfacen estas tres soluciones? Basta reemplazar t por cero
para ver que
3πx
u1 (x, 0) = sen
L
πx
u2 (x, 0) = 4,5 sen
L
3πx πx
u3 (x, 0) = sen + 4,5 sen .
L L
No parece haber mucha libertad para elegir condiciones iniciales generales. Nos ocupare-
mos de este tema a continuación.
108
www.elsolucionario.net
∂u ∂2u
EDP: =k 2
∂t ∂x
CB: u(0, t) = 0
u(L, t) = 0
3πx
CI: u(x, 0) = 4 sen .
L
2
L )
−( nπ kt
La solución producto u(x, t) = b sen nπx L
e satisface la condición inicial u(x, 0) =
b sen nπx
L
. Por lo tanto, eligiendo n = 3 y b = 4, se satisface la condición inicial, y luego
la solución es
3πx −( 3π 2
u(x, t) = 4 sen e L ) kt .
L
Se puede demostrar que este problema fı́sico tiene una única solución, por lo que no
importa el método utilizado para resolverlo. Siempre obtendrı́amos la misma.
www.elsolucionario.net
Principio de superposición. Las soluciones producto son soluciones muy particu-
lares, pues sólo sirven cuando la condición inicial tiene la forma apropiada u(x, 0) =
b sen nπx
L
. Sin embargo, este tipo de soluciones es muy útil también en muchas otras si-
tuaciones. Consideremos la misma EDP con las mismas condiciones de borde pero con la
condición inicial
3πx 8πx
u(x, 0) = 4 sen + 1,7 sen .
L L
La solución de este problema puede obtenerse sumando dos soluciones producto:
3πx −( 3π 2 8πx −( 8π 2
u(x, t) = 4 sen e L ) kt + 1,7 sen e L ) kt .
L L
Inmediatamente vemos que esta función resuelve la EDP, las CB y también la CI.
es también una solución. Como sabemos del método de separación de variables que la
nπ 2
función sen nπx e−( L ) kt es solución de la ecuación de difusión para todo entero n positivo,
L
entonces cualquier combinación lineal de estas soluciones es también solución de la misma
ecuación de difusión con las mismas condiciones de borde homogéneas. Por lo tanto
M
X nπx −( nπ 2
u(x, t) = bn sen e L ) kt (8.22)
n=1
L
109
www.elsolucionario.net
M
X nπx
u(x, 0) = bn sen .
n=1
L
Por lo que podemos resolver la ecuación de difusión para cualquier condición inicial de
la forma
M
X nπx
u(x, 0) = f (x) = bn sen .
n=1
L
¿Qué hacemos si tenemos que resolver la ecuación de difusión con un dato inicial f (x)
que no es de esta forma? La teorı́a de las series de Fourier (Capı́tulo 9) dice que:
1. Cualquier función f (x) (con algunas restricciones razonables) puede ser aproximada
www.elsolucionario.net
(en algún sentido) por combinaciones lineales de sen nπx
L
.
2. La aproximación tal vez no será buena para M pequeño, pero mejora a medida que
M crece.
Afirmamos entonces que cualquier dato inicial f (x) se puede escribir como una com-
binación lineal infinita de sen nπx
L
, conocida como serie de Fourier:
∞
X nπx
f (x) = bn sen . (8.23)
n=1
L
110
www.elsolucionario.net
Haciendo uso de la observación anterior vemos que todos los términos de esta serie son
nulos excepto uno, aquel para el que n = m. Por lo tanto
Z L Z L
mπx mπx mπx L
f (x) sen dx = bm sen sen dx = bm .
0 L 0 L L 2
www.elsolucionario.net
Z L
2 mπx
bm = f (x) sen dx, m = 1, 2, 3, . . . . (8.25)
L 0 L
111
www.elsolucionario.net
1.4
y=f
y=f1 1
1.2
y=f3
1 y=f30 0.8
0.8
0.6
0.6
0.4
0.4
www.elsolucionario.net
0.2
0.2
0 0
0 0.2 0.4 0.6 0.8 1 0 0.1 0.2 0.3 0.4 0.5 0.6
Cuando esto ocurre suele decirse que las CI y las CB son incompatibles. Sin embargo, si
graficamos la solución para diferentes valores de t, a partir de f30 , vemos que se observa
lo que se espera de la solución con condición inical constantemente igual a 1.
En la Figura 8.3 (izquierda) vemos el gráfico de
30
X 2
u30 (x, t) = bn e−(nπ) t sen nπx,
n=1
112
www.elsolucionario.net
como función de x para t = 0,01, 0,03, 0,05, 0,07, 0,09 (también se grafica f (x) ≡ 1 para
referencia, y no se grafica u30 (x, 0)). Se ve que la evolución es muy rápida en el intervalo
de tiempo [0, 0,01] ası́ que graficamos u(x, t) para t = 0,002, 0,004, 0,006, 0,008, 0,01 en el
gráfico de la derecha.
1 1
0.8 0.8
0.6 0.6
0.4 0.4
0.2 0.2
0 0
0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1
www.elsolucionario.net
Figura 8.3: Gráfico de u(x, t) para 0 < x < 1 y t = 0,01, 0,03, 0,05, 0,07, 0,09 (izquierda).
Gráfico de u(x, t) para 0 < x < 1 y t = 0,002, 0,004, 0,006, 0,008, 0,01 (derecha).
En los problemas se pide que usted haga lo mismo para diferentes funciones f y
problemas con otras condiciones de borde.
Supongamos que un alambre delgado se curva para formar una circunferencia, como
en la Figura 8.4. Por cuestiones de conveniencia supongamos que el alambre tiene longitud
2L, es decir, L es la mitad de la longitud del alambre. Como la circunferencia tiene una
2L L
longitud de 2πr, el radio es r = = . Si el alambre es suficientemente delgado y
2π π
113
www.elsolucionario.net
También el flujo de calor deberı́a ser continuo, y como la conductividad térmica es cons-
tante, esto implica que la derivada con respecto a x de la temperatura es continua:
www.elsolucionario.net
∂u ∂u
(−L, t) = (L, t). (8.29)
∂x ∂x
Las dos condiciones de borde para la ecuación en derivadas parciales son entonces (8.28)
y (8.29). La condición inicial es una función dada
El problema matemático consiste luego en la EDP lineal homogénea (8.27), con condicio-
nes de borde lineales y homogéneas (8.28) y (8.29). Es por lo tanto plausible la utilización
del método de separación de variables.
Proponemos una solución producto u(x, t) = φ(x)G(t). Luego obtenemos para φ y G
las ecuaciones diferenciales ordinarias
dG d2 φ
= −λkG = −λφ.
dt dx2
Las condiciones de borde que nos interesan ahora implican que la función φ debe satisfacer
dφ dφ
φ(−L) = φ(L) y (−L) = (L). La constante de separación λ se determina entonces
dx dx
encontrando todos los λ para los que el problema a valores en el borde:
d2 φ dφ dφ
= −λφ, φ(−L) = φ(L) (−L) = (L),
dx2 dx dx
tiene solución no-trivial. Las condiciones de borde con que nos encontramos se denominan
comúnmente condiciones de borde periódicas porque aunque el problema se piensa
definido en −L < x < L, se puede pensar que está definido periódicamente para todo
x ∈ R; la temperatura será periódica (x = x0 es el mismo punto fı́sico que x = x0 + 2L,
y por lo tanto tendrá la misma temperatura). Si λ > 0, la solución general del problema
para φ es √ √
φ(x) = c1 cos λx + c2 sen λx.
114
www.elsolucionario.net
Antes de resolver esta ecuación, analicemos la segunda condición de borde, que involucra
a la derivada,
dφ √ √ √
(x) = λ −c1 sen λx + c2 cos λx .
dx
dφ dφ
Luego (−L) = (L) se cumple si
www.elsolucionario.net
dx dx
√ √ √ √ √ √
λ −c1 sen λ(−L) + c2 cos λ(−L) = λ −c1 sen λL + c2 cos λL
115
www.elsolucionario.net
www.elsolucionario.net
tiplicándola por una constante). En cambio para cada autovalor positivo λ = (nπ/L)2
hay dos autofunciones linealmente independientes, sen nπx
L
y cos nπx
L
. Nuevamente, no hay
ningún autovalor λ < 0 (¿por qué?).
Aquı́ la función f (x) es una combinación lineal de senos y de cosenos. Otra diferencia
crucial es que (8.34) es válida para −L < x < L, mientras que para las series de senos
del ejemplo de la Sección 8.2 la igualdad era válida para 0 ≤ x ≤ L.
116
www.elsolucionario.net
www.elsolucionario.net
Los coeficientes se obtienen de forma similar a como se hizo antes. Si multiplica-
mos (8.34) por cos mπxL
y/o por sen mπx
L
(hacemos las dos cuentas a la vez) y luego inte-
gramos de x = −L a x = L, obtenemos
mπx mπx
Z L
cos
∞ Z L
cos
L X nπx L
f (x) dx = an cos dx
−L
mπx
n=0 −L L mπx
sen sen
L L
mπx
∞ Z L
cos
X nπx L
+ bn sen . dx.
n=1 −L L mπx
sen
L
Ahora utilizamos las propiedades integrales (8.35) y obtenemos
Z L Z L
mπx mπx 2
f (x) cos dx = am cos dx
−L L −L L
Z L Z L
mπx mπx 2
f (x) sen dx = bm sen dx
−L L −L L
Despejando obtenemos
Z L
1
a0 = f (x) dx
2L −L
Z
1 L mπx
am = f (x) cos dx
L −L L
Z
1 L mπx
bm = f (x) sen dx
L −L L
donde las dos últimas fórmulas valen para m = 1, 2, 3, . . . .
117
www.elsolucionario.net
www.elsolucionario.net
A través del principio del máximo.
donde g(x, t), f (x), a(t), b(t) son funciones dadas en [0, L] y en [0, T ], respectivamente.
Si u1 (x, t) y u2 (x, t) son soluciones de (8.36), entonces u1 (x, t) = u2 (x, t) para todo
0 ≤ x ≤ L, t ≥ 0.
Demostración. Sea v(x, t) = u1 (x, t) − u2 (x, t). La demostración consiste en mostrar que
v(x, t) ≡ 0, lo que implicará que u1 (x, t) ≡ u2 (x, t). ¿Qué podemos deducir que cumple
v(x, t)? Claramente
vt = kvxx , 0 < x < L, t > 0
v(0, t) = 0, v(L, t) = 0 t>0 (8.37)
v(x, 0) = 0, 0 < x < L.
Aquı́ hay que tener cuidado con la lógica. Es claro que si tomamos v ≡ 0 obtenemos una
solución de (8.37). Lo que queremos demostrar es que existe una única solución de (8.37),
y es la v ≡ 0.
118
www.elsolucionario.net
Esta función F (t) que depende sólo de t se llama energı́a de la solución v(x, t) a tiempo
t. Veamos qué propiedades cumple F (t).
Z L Z L
2
F (0) = v(x, 0) dx = 0 dx = 0.
0 0
Observemos F ′ (t)
Z
′ d L 2
F (t) = v(x, t) dx
dt 0
Z L
∂ 2
= v(x, t) dx
www.elsolucionario.net
0 ∂t
Z L Z L
= 2v(x, t) vt (x, t) dx = 2 v(x, t) kvxx (x, t) dx por (8.37)
0 0
Z L
= 2k 2v(x, t) vxx (x, t) dx
0
" x=L Z L #
= 2k v(x, t) vx (x, t) − vx (x, t) vx (x, t) dx int. por partes
x=0 0
Z L
2
= 2k v(L, t) vx (L, t) − v(0, t)vx (0, t) −2k vx (x, t) dx .
| {z } |0 {z }
=0 por CB ≥0
| {z }
≤0
Observación 8.4. El método del Teorema 8.3 puede usarse para todas las condiciones
de borde que hemos visto hasta ahora. Basta ver que siempre da menor o igual a cero el
término
v(L, t) vx (L, t) − v(0, t)vx (0, t)
cuando v = u1 − u2 es la resta entre dos soluciones con las mismas condiciones de borde.
Ver ejercicio 8.13
119
www.elsolucionario.net
Observación 8.5. El método del Teorema 8.3 puede usarse para comparar dos soluciones
con diferentes condiciones iniciales. Si u1 , u2 son soluciones de (8.36) con condición inicial
Z L
2
f1 , f2 , respectivamente, entonces definimos v := u1 − u2 y F (t) = v(x, t) dx como
0
antes. Luego F ′ (t) ≤ 0 y F (t) ≤ F (0), es decir
Z L
2
u1 (x, t) − u2 (x, t) dx = F (t)
0
Z L 2
≤ F (0) = u1 (x, 0) − u2 (x, 0) dx
0
Z L
2
= f1 (x) − f2 (x) dx.
0
En palabras, el error cuadrático medio entre u1 (·, t) y u2 (·, t) es decreciente y nunca excede
el error cuadrático medio inicial.
Observación 8.6. La observación anterior nos dice que el error que se comete al apro-
ximar el dato inicial verdadero por uno medido, o por uno dado por una suma finita
www.elsolucionario.net
de senos (por ejemplo) no crece cuando lo medimos como la integral de la diferencia al
cuadrado. Esto no nos dice que el máximo error no crezca, sino que en promedio no crece.
Un resultado más preciso se presenta en la próxima sección.
y
M = máx{f , a, b} = máximo valor de u en Γ.
Aquı́ Γ = {0 ≤ x ≤ L, t = 0} ∪ {x = 0, 0 ≤ t ≤ T } ∪ {x = L, 0 ≤ t ≤ T } es lo que se
denomina frontera parabólica del problema (8.38) (ver Figura 8.5). Entonces
u(x, t) ≤ M, 0 < x < L, 0 < t ≤ T.
O lo que es lo mismo
máx u(x, t) = máx u(x, t), ó máx u = máx u.
0≤x≤L (x,t)∈Γ Ω×[0,T ] Γ
0≤t≤T
120
www.elsolucionario.net
T Ω × (0, T ]
x
L
Ω = [0, L] Ω = (0, L)
www.elsolucionario.net
Observación 8.8. Desde el punto de vista fı́sico, si u denota temperatura, el principio
del máximo nos dice que, en ausencia de fuentes volumétricas de calor, la temperatura
en un punto x a tiempo t es siempre menor o igual a la máxima temperatura alcanzada
entre el tiempo inicial y el borde en el intervalo de tiempo considerado.
Postponemos la demostración del Teorema 8.7 hasta la próxima sección y nos dedi-
camos ahora a ver algunas consecuencias que se desprenden del principio del máximo.
Corolario 8.9 (Principio del mı́nimo). Sea u una solución de (8.38) y definamos
y
m = mı́n{f , a, b} = mı́nimo valor de u en Γ.
Entonces
u(x, t) ≥ m, 0 < x < L, 0 < t ≤ T.
O lo que es lo mismo
Este corolario se demuestra aplicando el principio del máximo a −u, observando que
máx(−u) = − mı́n u, máx(−a) = − mı́n a, etc.
121
www.elsolucionario.net
Ejemplo 8.10. Los principios del máximo y del mı́nimo se pueden utilizar para deducir
cotas ajustadas de la solución aún cuando la misma no pueda calcularse exactamente.
Por ejemplo, si u(x, t) es la solución de
ut = 9uxx , 0 < x < 3, t > 0
u(0, t) = 0, u(3, t) = 0 t>0
x4 (x − 3)4
u(x, 0) = , 0 < x < 3,
100
entonces u(x, t) es siempre un número entre 0 y el máximo del dato inicial. Observemos
4 4
que el dato inicial x (x−3)
100
es simétrico en el intervalo [0, 3] y positivo, y que su máximo se
alcanza en el punto medio del intervalo. Reemplazando x = 1,5 en la fórmula obtenemos
1,58
que el máximo del dato inicial es ≈ 0,25629. Por lo tanto, deducimos que
100
0 ≤ u(x, t) ≤ 0,25629, 0 < x < 3, t > 0.
Una consecuencia importante de los principios del máximo y del mı́nimo es el siguiente
www.elsolucionario.net
teorema.
Teorema 8.11 (Dependencia continua de las CI y las CB). Sean u1 y u2 soluciones de
los problemas
i i
ut − kuxx = 0,
0 < x < L, 0 < t < T
i i i i
u (0, t) = a (t), u (L, t) = b (t) 0<t<T , i = 1, 2.
i i
u (x, 0) = f (x), 0 < x < L,
Si para algún error ǫ > 0 se cumple que
|f1 (x) − f2 (x)| ≤ ǫ, 0<x<L
|a1 (t) − a2 (t)| ≤ ǫ, 0<t<T
|b1 (t) − b2 (t)| ≤ ǫ, 0<t<T
entonces
|u1 (x, t) − u2 (x, t)| ≤ ǫ, 0 < x < L, 0 < t < T.
Demostración. La demostración consiste en aplicar los principios del máximo y mı́nimo
a v(x, t) = u1 (x, t) − u2 (x, t). La función v(x, t) ası́ definida es solución de
122
www.elsolucionario.net
Por lo tanto, aplicando los principios del máximo y del mı́nimo a v(x, t) obtenemos
Observación 8.12. El último teorema es muy importante, pues nos dice que si dispo-
nemos de datos iniciales y de borde aproximados (como los obtenidos por medición),
entonces la solución del problema está igualmente aproximada a la solución con los datos
exactos.
Observación 8.13. Por otro lado, nos permite saber que las soluciones uN , calculadas
con una aproximación fN del dato inicial f (por ejemplo la obtenida con una suma de N
de funciones trigonométricas), están tan aproximadas a la solución exacta u(x, t) como
fN lo está de f .
www.elsolucionario.net
En esta sección procedemos a demostrar el Teorema 8.7, conocido como Principio del
Máximo. Si bien el resultado es obvio si uno piensa en temperaturas y desde el punto de
vista fı́sico, es interesante ver qué propiedades de la fı́sica observables u obvias pueden
deducirse del modelo matemático.
Demostremos primero la siguiente proposición, que nos da el principio del máximo en
un caso particular en el que hay una fuente constantemente negativa.
Proposición 8.14. Sea v(x, t) una solución del problema
Demostración. Sea M = máx v(x, t) y sea M0 = máx v(x, t). Es obvio que M ≤
(x,t)∈Γ (x,t)∈Ω×[0,T ]
M0 porque Γ ⊂ Ω × [0, T ]. Para probar el teorema debemos concluir que M0 ≤ M .
Supongamos entonces que M0 > M y veamos que esto no es posible.
Sea (x0 , t0 ) el punto del dominio espacio temporal Ω × [0, T ] donde se alcanza el
máximo M0 . Como este máximo es mayor estricto que el máximo M que se alcanza sobre
Γ resulta que (x0 , t0 ) ∈
/ Γ, es decir,
123
www.elsolucionario.net
v t v
T Ω × (0, T ]
t0
(x0 , t0 )
x t
0 x0 L 0 t0 T
x
x0 L
Figura 8.6: Si el máximo de v está en (x0 , t0 ) entonces v tiene un máximo en x0 sobre el
segmento horizontal t = t0 y un máximo en t0 sobre el segmento vertical x = x0 . Luego
vx (x0 , t0 ) = 0, vxx (x0 , t0 ) ≤ 0 y vt (x0 , t0 ) ≥ 0.
www.elsolucionario.net
segmento V se da en un t0 que podrı́a ser igual a T , sólo podemos decir que vt (x0 , t0 ) ≥ 0,
pero entonces
vt (x0 , t0 ) − kvxx (x0 , t0 ) ≥ 0,
lo cual contradice la hipótesis de que v es solución de (8.39). Esta contradicción nos dice
que no es posible que M0 > M , y luego M0 ≤ M . Como ya sabı́amos que M ≤ M0
concluimos que M = M0 .
Demostración del Teorema 8.7. Supongamos ahora que u(x, t) satisface las hipótesis del
Teorema 8.7. Para c > 0 (arbitrario) definimos v(x, t) = u(x, t) − c t. Entonces vt =
k vxx − c o, lo que es lo mismo
vt (x, t) − kvxx (x, t) = −c.
Por la Proposición 8.14, para cada (x, t) ∈ Ω × (0, T ]
u(x, t) − c t = v(x, t) ≤ máx v ≤ máx u,
Γ Γ
El Teorema 8.7 se conoce también como Principio del máximo débil. Existe un prin-
cipio del máximo fuerte que enunciamos a continuación, cuya demostración queda fuera
del alcance de este curso (puede encontrarse en [Bleecker-Csordas 1996, Chapter 3].
Teorema 8.15 (Principio del máximo fuerte). Bajo las hipótesis del Teorema 8.7, no
existe ningún (x0 , t0 ) ∈ Ω × (0, T ] tal que u(x0 , t0 ) = M := máxΓ u, a menos que u sea
constante en [0, L] × [0, t0 ], y en ese caso f ≡ a ≡ b ≡ u ≡ M para 0 ≤ x ≤ L y
0 ≤ t ≤ t0 .
124
www.elsolucionario.net
La idea es hallar primero una solución particular que cumpla la Ecuación Diferencial
(ED) y las Condiciones de Borde (CB), ignorando la Condición Inicial (CI). Luego
hallar una que cumpla la ED y las CB homogéneas correspondientes de manera que
sumada a la anterior cumpla todo: la ED, las CB y la CI.
www.elsolucionario.net
con a y b dos constantes y g(x) una función dada. Descomponemos este problema de la
siguiente manera. Primero resolveremos el problema:
up,t − kup,xx = 0, 0 < x < L, t > 0,
up (0, t) = a, up (L, t) = b, t > 0,
up (x, 0) = lo que sea, 0 < x < L.
Una vez resuelto, resolvemos este:
vt − kvxx = 0, 0 < x < L, t > 0,
v(0, t) = 0, v(L, t) = 0, t > 0,
v(x, 0) = g(x) − up (x, 0), 0 < x < L.
De esta manera, la función u(x, t) = up (x, t) + v(x, t) será solución de (8.40). En efecto:
ut − kuxx = (up + v)t − k (up + v)xx
= up,t + vt − kup,xx − kvxx 0 < x < L, t > 0,
= up,t − kup,xx + vt − kvxx = 0,
u(0, t) = up (0, t) + v(0, t) = a + 0 = a, t > 0,
u(L, t) = up (L, t) + v(L, t) = b + 0 = b, t > 0,
u(x, 0) = up (x, 0) + v(x, 0)
= up (x, 0) + (g(x) − up (x, 0)) = g(x), 0 < x < L.
Ahora nos resta resolver los dos (sub)problemas. Para hallar up conviene intentar primero
buscar una solución estacionaria, es decir una que satisfaga up,t = 0. Reemplazando up,t
por 0 en la ED nos da up,xx = 0, o sea que up es lineal en x. Imponiendo las CB obtenemos
b−a
up (x, t) = x + a,
L
125
www.elsolucionario.net
www.elsolucionario.net
L n=1
L
RL
con bn = L2 0 g(x) − b−a L
x + a sen nπx
L
dx, n = 1, 2, . . . . Es importante notar que
los coeficientes bn no se calculan con respecto a la g(x) original, sino con respecto a la
condición inicial que satisface v(x, t), que es el dato inicial original g(x) menos lo que vale
up en t = 0.
Resolvamos un ejemplo más concreto
Ejemplo 8.17. Hallar la solución de
ut − 2uxx = 0, 0 < x < 3, t > 0,
u(0, t) = 2, u(3, t) = −1, t > 0, (8.41)
1 π
u(x, 0) = 2 + sen( x) − x, 0 < x < 3,
10 3
Buscamos una solución estacionaria up = up (x) que cumpla las CB: up = C1 x + C2 , con
up (0) = 2 y up (3) = −1, es decir
up (x) = 2 − x.
Ahora buscamos una solución v(x, t) del problema con CB homogéneas y dato inicial
1 π 1 π
2+ sen( x) − x − 2 − x = sen( x),
10 3 10 3
vt − 2vxx = 0, 0 < x < 3, t > 0,
v(0, t) = 0, v(3, t) = 0, t > 0,
1 π
v(x, 0) = sen( x), 0 < x < 3.
10 3
Por lo visto en la Sección 8.2 la solución de este problema es
1 π π 2
v(x, t) = sen( x)e−( 3 ) 2t ,
10 3
126
www.elsolucionario.net
y finalmente
1 π 2π 2
u(x, t) = up (x, t) + v(x, t) = 2 − x + sen( x)e− 9 t .
10 3
Observamos que
lı́m u(x, t) = 2 − x = up (x, t).
t→∞
Veamos otro ejemplo, con condiciones de borde de diferente tipo en cada extremo:
Ejemplo 8.18. Hallar la solución de
ut − 3uxx = 0, 0 < x < 1, t > 0,
ux (0, t) = 1, u(1, t) = −1, t > 0,
(8.42)
2 3π 5
u(x, 0) = x + cos x − , 0 < x < 1.
4 2
Buscamos, como antes, una solución particular estacionaria up = up (x), que debe ser por
lo tanto lineal: up (x) = C1 x + C2 . Para que satisfaga las CB debe ser C1 = 0 y C2 = −1,
es decir,
www.elsolucionario.net
up (x) = x − 2.
Ahora buscamos una solución v(x, t) del problema con CB homogéneas y dato inicial
2 3π 5 1 2 3π 1 3π
x + cos x − − (x − 2) = − + cos x = cos x ,
4 2 2 4 2 2
donde hemos usado la identidad trigonométrica cos2 (α) = 21 + 12 cos(2α). Es decir, busca-
mos v(x, t) solución de
vt − 3vxx = 0, 0 < x < 1, t>0
vx (0, t) = 0, v(1, t) = 0, t > 0,
1 3π
v(x, 0) = cos x , 0 < x < 3.
2 2
Por lo visto en el Problema 8.9 la solución de este problema es
1 3π 3π 2
v(x, t) = cos( x)e−( 2 ) 3t ,
2 2
y finalmente
1 3π 27π 2
u(x, t) = up (x, t) + v(x, t) = x − 2 + cos( x)e− 4 t .
2 2
Observamos que
lı́m u(x, t) = x − 2 = up (x, t).
t→∞
El truco consiste entonces en hallar la solución de estado estacionario que cumple las
condiciones de borde para up y luego el problema se transforma en uno de los conocidos
con CB homogéneas. Esto funciona bien con cualquier combinación de CB constantes en
el tiempo, excepto cuando las dos condiciones son de tipo Neumann. Veamos un ejemplo
semi-abstracto que ilustra la situación.
127
www.elsolucionario.net
www.elsolucionario.net
siguiente forma
up (x, t) = h(x) + c t.
Veamos si hallamos una solución de esta forma que cumpla la ED y las CB:
up,t (x, t) = c,
up,xx (x, t) = h′′ (x),
up,t (x, t) − kup,xx (x, t) = c − kh′′ (x),
up,x (0, t) = h′ (0),
up,x (L, t) = h′ (L).
La ED y las CB impuestas a up implican una EDO con CB para h(x):
c
h′′ (x) = , 0 < x < L, h′ (0) = a, h′ (L) = b.
k
Intentemos resolverla:
c c 2
h′′ (x) = =⇒ h(x) = x + C1 x + C2
k 2k
c
h′ (x) = x + C1
k
′ c 2
h (0) = a =⇒ C1 = a =⇒ h(x) = x + a x + C2
2k
c
h′ (x) = x + a
k
′ c b−a
h (L) = b =⇒ L + a = b =⇒ c = k
k L
y nos quedó C2 con libertad de elección. Para simplificar elegimos C2 = 0 y nos queda
b−a 2 b−a
up (x, t) = x + |{z}
a x+ kt
| 2L
{z } C1 |L
{z }
c c
2k
128
www.elsolucionario.net
Luego nos falta hallar la solución v(x, t) del problema correspondiente, con condición
inicial g(x) − up (x, 0), que en este caso resulta
vt − kvxx = 0, 0 < x < L, t>0
vx (0, t) = 0, vx (L, t) = 0, t > 0,
b−a 2
v(x, 0) = g(x) − x + ax , 0 < x < L.
2L
La solución de este problema está dada por
∞
X nπx −( nπ 2
v(x, t) = a0 + an cos e L ) kt ,
n=1
L
con
Z L Z L
1 b−a 2 2 b−a 2 nπx
a0 = g(x)− x + ax , an =
g(x) − x + ax cos dx, n = 1, 2, . .
L 0 2L 0 L 2L L
(8.44)
Finalmente, la solución u(x, t) del problema (8.43) es u(x, t) = up (x, t) + v(x, t), o sea
www.elsolucionario.net
∞
b−a 2 b−a X nπx −( nπ 2
u(x, t) = x + ax + k t + a0 + an cos e L ) kt ,
2L L n=1
L
con an , n = 0, 1, 2, . . . dados por (8.44).
129
www.elsolucionario.net
www.elsolucionario.net
Llegamos entonces a la siguiente proposición:
Proposición 8.20. La solución de (8.45) está dada por
b(t) − a(t)
con w(x, t) = x + a(t), u1 (x, t) y u2 (x, t) soluciones de (8.46a) y (8.46b),
L
respectivamente, con
b(0) − a(0)
g(x) = f (x) − w(x, 0) = f (x) − x + a(0) en (8.46a),
L
′
b (t) − a′ (t) ′
h(x, t) = − (wt − kwxx ) = − x + a (t) en (8.46b).
L
El objetivo de esta sección es presentar el Principio de Duhamel que provee una forma de
resolver este tipo de problemas resolviendo una familia de problemas sin término fuente,
pero con dato inicial no-homogéneo.
130
www.elsolucionario.net
www.elsolucionario.net
v1 (0, t) = 0, v1 (L, t) = 0 t > t1 ,
v1 (x, t1 ) = ∆t h(x, t1 ), 0 < x < L, (t = t1 ).
En el tercer intervalo de tiempo (entre t2 y t3 ) pensamos que hay difusión del calor que
ingresó durante los dos primeros intervalos, pero es despreciable la difusión del calor que
ingresa durante el último intervalo, ası́:
Z t3
∼
u(x, t3 ) = v1 (x, t3 ) + v2 (x, t3 ) + h(x, t) dt ∼
= v1 (x, t3 ) + v2 (x, t3 ) + ∆t h(x, t3 ),
t2
131
www.elsolucionario.net
www.elsolucionario.net
Entonces la única solución de
ut − kuxx = h(x, t),
0 < x < L, 0 < t < T,
u(0, t) = 0, u(L, t) = 0, 0 < t < T,
u(x, 0) = 0, 0 < x < L.
está dada por la fórmula Z t
u(x, t) = w(x, t; s) ds.
0
La demostración es sencilla, pero para hacerla debemos ver antes el siguiente lema.
Lema 8.22. Dada una función g(t, s) de dos variables tal que g(t, s) y gt (t, s) son con-
tinuas, se cumple que
Z t Z t
d
g(t, s) ds = g(t, t) + gt (t, s) ds
dt 0 0
Ry
Demostración. Para demostrar este lema definimos G(r, y) = 0 g(r, s) ds y vemos que
Z t
g(t, s) ds = G(t, t).
0
Además Z y
Gr (r, y) = gt (r, s) ds y Gy (r, y) = g(r, y).
0
Ahora bien,
Z t
d
G(t, t) = Gr (t, t) + Gy (t, t) = Gr (r, y) + Gy (r, y) = gt (t, s) ds + g(t, t).
dt r=t r=t 0
y=t y=t
1
Notemos que la condición inicial está dada en t = s y es allı́ donde entra la función h(x, t).
132
www.elsolucionario.net
www.elsolucionario.net
Es decir, ut (x, t) − kuxx (x, t) = h(x, t).
La fórmula puede ser bastante complicada, pero veamos un ejemplo cuando resolvemos
en R usando la solución fundamental.
Entonces
Z tZ ∞
u(x, t) = h(y, s)Φ(x − y, t − s) dy ds, x ∈ R, t > 0.
0 −∞
133
www.elsolucionario.net
es
Z ∞ Z tZ ∞
u(x, t) = g(y)Φ(x − y, t) dy + h(y, s)Φ(x − y, t − s) dy ds, x ∈ R, t > 0.
−∞ 0 −∞
8.8. Ejercicios
8.1. Demostrar que la solución fundamental de la ecuación de difusión en una dimensión
dada por (8.2) es solución de ut = kuxx en R × (0, ∞).
www.elsolucionario.net
* 8.2. Consideremos la solución fundamental de la ecuación de difusión en una dimensión
dada por (8.2). Demostrar que si g : R → R es acotada y de soporte compacto (es decir,
existe M > 0 tal que |g(x)| ≤ M , x ∈ R; y g(x) = 0 si |x| > M ), entonces la función
u(x, t) definida por (8.3) cumple lo siguiente:
(a) u ∈ C 2 (R × (0, ∞)) (Demostrarlo primero para Φ y ver que bajo la hipótesis de g
acotada y de soporte compacto se puede meter la derivada en la integral);
(b) ut = kuxx para (x, t) ∈ R × (0, ∞); (Usar que Φ es solución y usar—nuevamente—
que se puede meter la derivada en la integral);
R 2 √ R 2
R 2
1/2
(i) Demostrar que R
e−x dx = π. Para esto usar que R e−x dx = R2 e−|x| dx y calcular
la integral de R2 usando coordenadas polares.
R
(ii) Demostrar que R Φ(x, t) dx = 1 para todo t > 0.
R
(iii) Demostrar que |x|<δ Φ(x, t) dx → 1 cuando t → 0+ , para cualquier δ > 0 fijo.
Comentario: Se puede demostrar que la ecuación de difusión en R tiene una única solución
acotada para cada dato inicial g. Este problema nos dice que la fórmula para dicha solución
está dada por (8.3).
134
www.elsolucionario.net
www.elsolucionario.net
0
Z L
(n − 21 )πx (m − 21 )πx
(c) sen( ) sen( ) dx, para n, m = 1, 2, 3, . . . .
0 L L
Z L
(n − 12 )πx (m − 21 )πx
(d) cos( ) cos( ) dx, para n, m = 1, 2, 3, . . . .
0 L L
8.4. En base al problema anterior, decir cómo deben calcularse los coeficientes an , bn (en
términos de f ) si se sabe que se cumplen las igualdades siguientes para f definida en
[0, L]:
N N
X nπx X (n − 21 )πx
(a) f (x) = bn sen( ). (c) f (x) = cn sen( ).
n=1
L n=1
L
N N
X nπx X (n − 12 )πx
(b) f (x) = a0 + an cos( ). (d) f (x) = dn cos( ).
n=1
L n=1
L
135
www.elsolucionario.net
www.elsolucionario.net
ut = 5uxx ,
0 < x < 10, t > 0,
ux (0, t) = 0, ux (10, t) = 0, t > 0.
u(x, 0) = f (x), 0 < x < 10,
en cada uno de los siguientes casos
π 10x − x2
(a) f (x) = 10; (b) f (x) = 10 − cos x ; (c) f (x) = 4 + .
10 2
8.9. Hallar todas las soluciones u(x, t) = φ(x)G(t) de variables separadas del problema:
136
www.elsolucionario.net
www.elsolucionario.net
8.14. Usar los principios del máximo/mı́nimo para deducir cotas para la solución de
>> x = [0:.01:5];
>> fe = x .* (5-x);
>> f1 = 200/pi^3*sin(pi/5*x);
>> plot(x, fe, x, f1);
>> legend(’exacta’, ’aproximada’)
137
www.elsolucionario.net
(c) Tomando una grilla de aproximadamente 500 puntos, calcular el máximo error entre
fe y f1 (usar valor absoluto). En MATLAB: >> max(abs(fe-f1))
También se puede hallar el máximo y el mı́nimo de fe (x) − f1 (x) analı́ticamente, pero no vale la
pena en este caso.
(d) ¿Cuál es la solución u1 (x, t) del problema tomando como condición inicial la función
f1 (x)?
(e) Si se toma u1 (x, t) como aproximación de la solución exacta ue (x, t). ¿Cuál es el error
entre ue (x, t) y u1 (x, t) para 0 ≤ x ≤ 5 y t ≥ 0? ¿Sirve u1 (x, t) como aproximación
de ue (x, t)?
(g) Graficar la temperatura en el punto medio del intervalo [0, 5] para t en el intervalo
[0, 7], usando la aproximación u1 (x, t).
8.16. Repetir el ejercicio anterior con f3 (x) = 200
π 3 sen πx
5
+ 200
27π 3 sen 3πx
5
.
www.elsolucionario.net
* 8.17. Sea u una solución C 2 ([0, L] × [0, T ]) del problema
Demostrar que si
M = máx máx a(t), máx f (x)
t∈[0,T ] x∈[0,L]
Demostrar que si
M = máx máx f (x)
x∈[0,L]
138
www.elsolucionario.net
Demostrar que si
M = máx máx a(x, t), máx f (x)
x∈∂Ω,t∈[0,T ] x∈Ω
www.elsolucionario.net
0 < x < 10.
20
8.22. Demostrar la Proposición 8.20.
8.23. Considerar el problema siguiente:
ut = kuxx , 0 ≤ x ≤ L, t ≥ 0,
ux (0, t) = a(t), ux (L, t) = b(t), t ≥ 0, (8.51)
u(x, 0) = f (x), 0 ≤ x ≤ L.
(b) Dar una fórmula para una función w(x, t) que cumpla las condiciones de borde.
(c) ¿Qué problema debe resolver v(x, t) para que u(x, t) = v(x, t) + w(x, t) sea solución
de (8.51)?
(d) Realizar una descomposición del problema descripto en el ı́tem anterior en dos sub-
problemas más sencillos, análogos a (8.46a) y (8.46b). Dar las fórmulas precisas
para los datos no-homogéneos de los problemas resultantes (g(x) y h(x, t) de (8.46a)
y (8.46b) respectivamente).
(e) Demostrar que la suma de w(x, t) y las soluciones u1 (x, t), u2 (x, t) de los subproble-
mas del ı́tem anterior, es solución del problema original.
8.24. Repetir los cinco ı́tems el ejercicio anterior para la siguiente situación:
139
www.elsolucionario.net
Bibliografı́a complementaria
[Bleecker-Csordas 1996] Bleecker, D., Csordas, G., Basic Partial Differential Equations,
International Press, Cambridge, Massachusetts, 1996.
[Haberman 1998] Haberman, R., Elementary Applied Partial Differential Equations, Pren-
tice Hall, Upper Saddle River, NJ, 1998.
www.elsolucionario.net
140
www.elsolucionario.net
Capı́tulo 9
Series de Fourier
En el capı́tulo anterior vimos que las soluciones de ecuaciones de difusión (con coe-
ficientes constantes) sobre intervalos acotados se escriben como combinaciones lineales
de senos y cosenos, multiplicadas por exponenciales. Al tomar t = 0 para verificar la
condición inicial, las exponenciales dan todas igual a uno y desaparecen de la expresión.
www.elsolucionario.net
Afirmamos en ese momento que habı́a una forma de determinar los coeficientes an o bn pa-
ra que la suma de senos o cosenos aproximase al dato inicial f (x) tanto como deseáramos.
Junto con el principio del máximo, vimos que la suma correspondiente una vez vueltas
a incorporar las exponenciales, serı́a una solución igualmente aproximada de la solución
exacta para todo x y para todo t > 0. Nos proponemos ahora estudiar un poco más en
detalle la aproximación de funciones por series de senos y/o de cosenos.
El objetivo de este capı́tulo es estudiar cómo aproximar funciones f (x) por sumas de
alguna de las siguientes formas:
N N N
X nπx X nπx X nπx nπx
a0 + an cos , bn sen , a0 + an cos + bn sen .
n=1
L n=1
L n=1
L L
Estudiaremos primero el caso más general (que engloba de alguna manera a los otros
dos) que consiste en aproximar una función f definida en [−L, L] por una suma de la
forma
N
X nπx nπx
a0 + an cos + bn sen . (9.1)
n=1
L L
141
www.elsolucionario.net
N
X Z L
x̄ · ȳ = xi yi hf, gi = f (x)g(x) dx
i=1 −L
N
X 21 Z L 21
√ p
|x̄| = x̄ · x̄ = x2i kf k = hf, f i = 2
f (x) dx
i=1 −L
Distancia
N
X 21 Z L 21
www.elsolucionario.net
2 2
|x̄ − ȳ| = (xi − yi ) kf − gk = [f (x) − g(x)] dx
i=1 −L
Ortogonalidad o perpendicularidad
f ⊥g ⇐⇒ hf, gi = 0
Z L
x̄ ⊥ ȳ ⇐⇒ x̄ · ȳ = 0
⇐⇒ f (x)g(x) dx = 0
−L
Desigualdad de Cauchy-Schwarz
142
www.elsolucionario.net
Retomemos el estudio de las funciones senos y cosenos del capı́tulo anterior, es decir,
consideremos la sucesión
πx πx 2πx 2πx 3πx 3πx
1 , cos , sen , cos , sen , cos , sen , ...
| {zL} | {zL} | {zL } | {zL } | {zL } | {zL }
|{z}
c0
c1 (x) s1 (x) c2 (x) s2 (x) c3 (x) s3 (x)
Vimos en ese capı́tulo que se cumplen las siguientes relaciones (ver (8.35))
(
0 si n 6= m
hcn , cm i = δn,m L = , n, m = 1, 2, 3, . . . ,
L si n = m
hsn , sm i = δn,m L, n, m = 1, 2, 3, . . . ,
hcn , sm i = 0, n, m = 1, 2, 3, . . . ,
hc0 , sn i = 0, n = 1, 2, 3, . . . ,
hc0 , cn i = 0, n = 1, 2, 3, . . . ,
hc0 , c0 i = 2L.
www.elsolucionario.net
En palabras:
143
www.elsolucionario.net
Si definimos
c0 cn sn
c̃0 = , c̃n = , s̃n = , n = 1, 2, 3, . . . ,
kc0 k kcn k ksn k
entonces la fórmula anterior puede reescribirse como
N
X N
X
f (x) = hf, c̃0 ic̃0 (x) + hf, c̃n ic̃n (x) + hf, s̃n is̃n (x),
n=1 n=1
www.elsolucionario.net
−L L −L L
Figura 9.1: Gráficas de funciones que no son de la forma (9.3).
Sin embargo, casi siempre podemos definir los coeficientes de Fourier, y su correspon-
diente serie.
Definición 9.3 (Serie de Fourier). Sea f (x) definida sobre el intervalo [−L, L] tal que
las siguientes integrales están bien definidas:
Z L
1 1
a0 = f (x) dx = hf, c0 i,
2L −L 2L
Z
1 L nπx 1
an = f (x) cos dx = hf, cn i, n = 1, 2, . . . , N,
L −L L L
Z
1 L nπx 1
bn = f (x) sen dx = hf, sn i, n = 1, 2, . . . , N.
L −L L L
Definimos entonces la serie de Fourier de f por
∞
X nπx nπx
SF f (x) = a0 + an cos + bn sen ,
n=1
L L
144
www.elsolucionario.net
o lo que es lo mismo
SF f (x) = f (x),
o también
www.elsolucionario.net
∞
X nπx nπx
f (x) = a0 + an cos + bn sen .
n=1
L L
Lo que también es interesante es poder medir el error entre SFN f (x) y f (x), que también
estudiaremos en esta sección y en las que siguen.
El primer resultado de aproximación de SFN f a f , usando la norma integral k · k
definida al principio de esta sección está dado en la siguiente proposición.
Proposición 9.6. Sea VN el conjunto de todas las funciones que pueden escribirse como
combinación lineal de c0 , c1 , s1 , c2 , s2 , . . . , cN , sN , entonces
Para m = 1, 2, . . . , N ,
hf, cm i
= hf, cm i − am L = hf, cm i − L = 0,
L
y análogamente,
hf, sm i
= hf, sm i − bm L = hf, sm i − L = 0.
L
www.elsolucionario.net
Resumiendo, hf − SFN f, cn i = 0, hf − SFN f, cn i, y hf − SFN f, sn i, para n = 1, 2, . . . , N .
Si
PNahora v es cualquier función del conjunto VN , entonces se escribe como v = ā0 c0 +
n=1 ān cn + b̄n sn , y por lo tanto
N
X
hf − SF f, vi = ā0 hf − SF f, c0 i + ān hf − SF f, cn i + b̄n hf − SF f, sn i = 0.
n=1
kf − SFN f k2
= hf − SFN f, f − SFN f i definición de norma
= hf − SFN f, f − g + g − SFN f i sumamos y restamos g ∈ VN
= hf − SFN f, f − gi + hf − SFN f, g − SFN f i linealidad del producto escalar
| {z }
=v∈VN
= hf − SFN f, f − gi f − SFN f ⊥ v ∈ VN
≤ kf − SFN f k kf − gk. Cauchy-Schwarz (ver tabla)
kf − SFN f k ≤ kf − gk,
para cualquier g ∈ VN .
146
www.elsolucionario.net
Ejemplo 9.8. Hallar la serie de Fourier de f (x) = x en [−L, L]. Lo que debemos hacer
es calcular los coeficientes de Fourier a0 , a1 , b1 , a2 , b2 , . . .
Z L
1
a0 = x dx = 0,
2L −L
Z
1 L nπx
an = x cos dx = 0, n = 1, 2, . . . (función par por impar es impar)
L −L L
Z
1 L nπx 2L 2L
bn = x sen dx = − cos(nπ) = (−1)n+1 , n = 1, 2, 3, . . . .
L −L L nπ nπ
www.elsolucionario.net
2 2 2
1 1 1
0 0 0
−1 −1 −1
−2 −2 −2
−2 −1 0 1 2 −2 −1 0 1 2 −2 −1 0 1 2
2 2 2
1.5 1.5 1.5
1 1 1
0.5 0.5 0.5
0 0 0
−1 −1 −1
−2 −2 −2
−2 −1 0 1 2 −2 −1 0 1 2 −2 −1 0 1 2
usaron los valores de f (x) para −2 ≤ x ≤ 2 y que sin embargo las funciones SFN f (x) y
SF f (x) están definidas en todo R. En la figura 9.3 se muestra el gráfico de SF25 f (x) y
también el de SF f , comparado con f sobre el intervalo [−6, 6]. El fenómeno de oscilación
de magnitud grande de SFN f (x) en los extremos del intervalo se denomina Fenómeno de
Gibbs, y se debe a que para las condiciones que satisface f (x) = x son de alguna manera
incompatibles con las condiciones periódicas que cumplen los senos y cosenos.
Ejemplo 9.9. Consideremos ahora la función f (x) = 3 + (x − L)2 (x + L)2 , cuyos coefi-
147
www.elsolucionario.net
6 6
4 4
2 2
0 0
−2 −2
−4 −4
−6 −6
−6 −4 −2 0 2 4 6 −6 −4 −2 0 2 4 6
Figura 9.3: Gráficas de SF25 f (x) (izquierda) y SF f (derecha) para f (x) = x en [−2, 2],
pero graficadas sobre el intervalo [−6, 6].
www.elsolucionario.net
1
a0 = (8L4 + 45),
15
48 48
an = − 4 4 cos(nπ) = 4 4 (−1)n+1 , n = 1, 2, . . . ,
nπ nπ
bn = 0, n = 1, 2, . . . .
148
www.elsolucionario.net
4 4
3.8 3.8
3.6 3.6
3.4 3.4
3.2 3.2
3 3
4 4
3.8 3.8
3.6 3.6
www.elsolucionario.net
3.4 3.4
3.2 3.2
3 3
Figura 9.4: Gráficas de SFN f (x) para N = 1, 3, 5 y SF f para f (x) = 3 + (x − L)2 (x + L),
con L = 1.
10 10
9 9
8 8
7 7
6 6
5 5
4 4
3 3
−3 −2 −1 0 1 2 3 −3 −2 −1 0 1 2 3
Figura 9.5: Gráficas de SF10 f (x) (izquierda) y SF f (derecha) para f (x) = 3+(x−L)2 (x+
L) en [−1, 1] (L = 1), pero graficadas sobre el intervalo [−3, 3].
149
www.elsolucionario.net
Ejemplo 9.10. Consideremos ahora la función f (x) = |x|, cuyos coeficientes de Fourier
son
L
a0 =
2
2L 2L
an = 2 2 cos nπ − 1 = 2 2 (−1)n − 1 , n = 1, 2, . . . ,
nπ nπ
bn = 0, n = 1, 2, . . . .
1 1
0.8 0.8
www.elsolucionario.net
0.6 0.6
0.4 0.4
0.2 0.2
0 0
1 1
0.8 0.8
0.6 0.6
0.4 0.4
0.2 0.2
0 0
Figura 9.6: Gráficas de SFN f (x) para N = 1, 3, 5 y SF f para f (x) = |x|, con L = 1.
SF f (x) sólo se usaron los valores de f (x) para −L ≤ x ≤ L y que sin embargo las
funciones SFN f (x) y SF f (x) están definidas en todo R. En la figura 9.7 se muestra el
gráfico de SF5 f (x) y también el de SF f , comparado con f sobre el intervalo [−3, 3].
(
x, si x > 0,
f (x) = (x)+ :=
0, si x ≤ 0,
150
www.elsolucionario.net
3 3
2.5 2.5
2 2
1.5 1.5
1 1
0.5 0.5
0 0
−3 −2 −1 0 1 2 3 −3 −2 −1 0 1 2 3
Figura 9.7: Gráficas de SF10 f (x) (izquierda) y SF f (derecha) para f (x) = |x| en [−1, 1]
(L = 1), pero graficadas sobre el intervalo [−3, 3].
www.elsolucionario.net
L
a0 =
4
L L
an = 2 2 (cos nπ − 1) = 2 2 ((−1)n − 1), n = 1, 2, . . . ,
nπ nπ
L L
bn = − cos nπ = − (−1)n , n = 1, 2, . . . .
nπ nπ
En la figura 9.8 se muestra el gráfico de SFN f (x) para N = 1, 5, 15, y también el de SF f ,
comparado con f sobre el intervalo [−1, 1]. Notar nuevamente que para la definición de
SF f (x) sólo se usaron los valores de f (x) para −L ≤ x ≤ L y que sin embargo las
funciones SFN f (x) y SF f (x) están definidas en todo R. En la figura 9.9 se muestra el
gráfico de SF15 f (x) y también el de SF f , comparado con f sobre el intervalo [−3, 3].
151
www.elsolucionario.net
1 1
0.8 0.8
0.6 0.6
0.4 0.4
0.2 0.2
0 0
−0.2
−1 −0.5 0 0.5 1 −1 −0.5 0 0.5 1
1 1
0.8 0.8
0.6 0.6
www.elsolucionario.net
0.4 0.4
0.2 0.2
0 0
Figura 9.8: Gráficas de SFN f (x) para N = 1, 5, 15 y SF f para f (x) = (x)+ , con L = 1.
3 3
2.5 2.5
2 2
1.5 1.5
1 1
0.5 0.5
0 0
−3 −2 −1 0 1 2 3 −3 −2 −1 0 1 2 3
Figura 9.9: Gráficas de SF15 f (x) (izquierda) y SF f (derecha) para f (x) = (x)+ en [−1, 1]
(L = 1), pero graficadas sobre el intervalo [−3, 3].
152
www.elsolucionario.net
1.2
1
1
www.elsolucionario.net
0.8
0.8
0.6 0.6
0.4 0.4
0.2 0.2
0 0
−0.2
−1 −0.5 0 0.5 1 −1 −0.5 0 0.5 1
1 1
0.8 0.8
0.6 0.6
0.4 0.4
0.2 0.2
0 0
se usaron los valores de f (x) para −L ≤ x ≤ L y que sin embargo las funciones SFN f (x)
y SF f (x) están definidas en todo R. En la figura 9.11 se muestra el gráfico de SF15 f (x)
y también el de SF f , comparado con f sobre el intervalo [−3, 3].
153
www.elsolucionario.net
1 1
0.8 0.8
0.6 0.6
0.4 0.4
0.2 0.2
0 0
−3 −2 −1 0 1 2 3 −3 −2 −1 0 1 2 3
Figura 9.11: Gráficas de SF15 f (x) (izquierda) y SF f (derecha) para f (x) = (x)+ en [−1, 1]
(L = 1), pero graficadas sobre el intervalo [−3, 3].
Observemos primero que dada una función f (x) definida en [−L, L] la suma parcial
www.elsolucionario.net
de N términos de su serie de Fourier es
N
X nπx0 nπx0
SFN f (x0 ) = a0 + an cos + bn sen
n=1
L L
Z L
1
= f (x) dx
−L 2L
N Z L Z L
X 1 nπx nπx0 1 nπx nπx0
+ f (x) cos dx cos + f (x) sen dx sen
n=1 −L
L L L −L L L L
Z " N
#
1 L 1 X nπx nπx0 nπx nπx0
= f (x) + cos cos + sen sen dx
L −L 2 n=1 L L L L
Z " N
#
1 L 1 X nπ(x − x0 )
= f (x) + cos dx.
L −L 2 n=1 L
| {z }
DN (x−x0 )
Es decir, la suma parcial de orden N se puede expresar de manera compacta como una
integral involucrando al llamado núcleo de Dirichlet DN :
Z
1 L
SFN f (x0 ) = f (x)DN (x − x0 ) dx,
L −L
con
N
1 X nπx
DN (x) = + cos .
2 n=1 L
La primera propiedad que podemos ver que satisface el núcleo de Dirichlet es que
DN (x) es periódica de perı́odo 2L, es decir, DN (x + 2L) = DN (x) para todo x ∈ R y por
lo tanto su gráfica se repite en intervalos consecutivos de longitud 2L. Otra propiedad es
154
www.elsolucionario.net
la siguiente:
Z L Z L N Z L
1 X nπx
DN (x) dx = dx + cos dx = L,
−L −L 2 n=1 −L
L
y por lo tanto
Z L
1
DN (x) dy = 1, N = 1, 2, 3 . . . .
L −L
www.elsolucionario.net
valores de N , que puede observarse en la Figura 9.12. Vemos en la gráfica que la función
20 20
15 15
10 10
5 5
0 0
−5 −5
−1 −0.5 0 0.5 1 −1 −0.5 0 0.5 1
20 20
15 15
10 10
5 5
0 0
−5 −5
−1 −0.5 0 0.5 1 −1 −0.5 0 0.5 1
DN (x) es más oscilante para N más grande, pero se concentra alrededor del origen. Uno
RL
tenderı́a entonces a pensar que lı́mN →∞ L1 −L f (x)DN (x) dx = f (0) siempre que f (x).
Debido a las oscilaciones de DN (x), esta afirmación es cierta cuando además de continua,
f es C 1 a trozos, como lo establece la siguiente proposición.
155
www.elsolucionario.net
Proposición 9.13. Sea f (x) una función C 1 a trozos y continua en [−L, L]. Entonces
Z L
1
lı́m f (x)DN (x) dx = f (0).
N →∞ L −L
1
RL
Ahora sı́, dado que L −L
f (x)DN (x − x0 ) dx = SFN f (x0 ), tenemos que
Teorema 9.15. Sea f (x) una función C 1 a trozos y continua en [−L, L], que satisface
www.elsolucionario.net
f (−L) = f (L). Entonces
o también
SF f (x) = f (x), para cada x ∈ [−L, L],
que es lo mismo que decir
∞
X nπx nπx
a0 + an cos + bn sen = f (x), para cada x ∈ [−L, L].
n=1
L L
1 f (0− ) f (0+ )
y las integrales del lado derecho, multiplicadas por convergen a y respec-
L 2 2
tivamente, donde f (0− ) = lı́mx→0− f (x) y f (0+ ) = lı́mx→0+ f (x).
Motivados por esta observación hacemos la siguiente definición.
156
www.elsolucionario.net
Definición 9.16. Sea f (x) una función C 1 a trozos en [−L, L]. Cambiando los valores
de f (x) en algunos puntos obtenemos la función arreglada f¯(x) que se define como sigue:
f (x+ ) + f (x− )
, si − L < x < L,
f¯(x) = 2 (9.4)
f (−L+ ) + f (L− )
, si x = ±L,
2
donde f (x+
0 ) = lı́mx→x+
0
f (x) y f (x−
0 ) = lı́mx→x−
0
f (x).
Observación 9.17. Si f (x) es continua en [−L, L] entonces f¯(x) = f (x) para todo
f (−L+ ) + f (L− )
x ∈ (−L, L) (−L < x < L) y f¯(±L) = .
2
Si f (x) es continua y además f (−L) = f (L), entonces f¯(x) = f (x) para todo x ∈
[−L, L] (−L ≤ x ≤ L).
En la Figura 9.13 vemos algunos ejemplos gráficos de funciones f y sus correspon-
dientes funciones arregladas f¯.
www.elsolucionario.net
−L L −L L −L L
−L L −L L −L L
Figura 9.13: Ejemplos gráficos de funciones f (arriba) y sus correspondientes funciones
arregladas f¯ (abajo).
es decir,
N
X nπx nπx
lı́m a0 + an cos + bn sen = f¯(x), para cada x ∈ [−L, L].
N →∞
n=1
L L
157
www.elsolucionario.net
o,
N
X nπx nπx
lı́m a0 + an cos + bn sen = f˜¯(x), para cada x ∈ R.
N →∞
n=1
L L
En la Figura 9.14 se pueden ver las extensiones periódicas de las funciones arregladas
de la Figura 9.13.
−3L −L L 3L
www.elsolucionario.net
−3L −L L 3L
−3L −L L 3L
Figura 9.14: Gráficos de las extensiones periódicas de las funciones arregladas de la Figu-
ra 9.13, sobre el intervalo [−3L, 3L].
Observación 9.19. El Teorema 9.18 nos dice que para cada x ∈ [−L, L] (fijo) se cumple
que lı́mN →∞ SFN f (x) = f¯(x). Notemos que si f es discontinua en algún punto, o si no se
cumple que f (−L) = f (L), entonces no es cierto que máx−L≤x≤L | SFN f (x) − f¯(x)| → 0.
Esta convergencia es la que se conoce como convergencia uniforme, y se sabe que si una
sucesión de funciones continuas converge uniformemente a otra, entonces la función lı́mite
es también continua. Por lo tanto como las funciones SFN f (x) son continuas y valen lo
mismo en −L y L, si se cumpliera que máx−L≤x≤L | SFN f (x) − f¯(x)| → 0, resultarı́a
necesariamente que f es continua y f (−L) = f (L).
Observación 9.20 (Fenómeno de Gibbs). Cuando la función f tiene una discontinuidad
esencial (no evitable) o cuando f (L− ) 6= f (−L+ ), al aproximarla con su serie de Fourier
ocurre el llamado Fenómeno de Gibbs, por el que las sumas parciales de las series de
158
www.elsolucionario.net
Teorema 9.21. Sea f (x) una función continua y C 1 a trozos en [−L, L] tal que f (−L) =
f (L) entonces
máx |f (x) − SFN f (x)| → 0, cuando N → ∞.
−L≤x≤L
Observación 9.22. Este último teorema nos dice que si f es continua en [−L, L] y C 1
www.elsolucionario.net
a trozos, y además cumple la condición de periodicidad f (−L) = f (L) entonces su serie
de Fourier converge a f uniformemente. El máximo error entre SFN f (x) y f (x) tiende
a cero cuando N → ∞. Vimos que esto no ocurre cuando f es discontinua o cuando
f (−L) 6= f (L) debido al fenómeno de Gibbs.
A continuación vemos un teorema que bajo hipótesis más fuertes sobre f estable-
ce la convergencia uniforme de la serie de Fourier pero además presenta un orden de
convergencia.
Teorema 9.23. Sea f (x) una función de clase C 2 en el intervalo [−L, L] tal que
4L2 M 1
|f (x) − SFN f (x)| ≤ , para todo x ∈ [−L, L].
π2 N
Observación 9.24. La afirmación precedente nos dice que si f tiene derivadas primeras
y segundas continuas en [−L, L] y además cumple las condiciones de compatibilidad
(periodicidad) de igualdad de f y también de f ′ en los extremos, entonces el error entre
f (x) y la suma de N términos de su serie de Fourier es menor que una constante por
1 1
. Esto nos dice que dicho error tiende a cero uniformemente, tan rápido como , y
N N
también nos da una manera de elegir N para aproximar con una tolerancia dada.
Demostración del Teorema 9.23. Por el Teorema 9.21 se cumple que SF f (x) = f (x) para
159
www.elsolucionario.net
www.elsolucionario.net
Veamos ahora si podemos acotar el tamaño de los2 coeficientes an , bn . Recordemos que
cn (x) = cos nπx
L
y que entonces cn (x) = − nπ L
c′′n (x). Luego, integrando por partes
obtenemos Z Z L
1 L L
an = f (x)cn (x) dx = − 2 2 f (x)c′′n (x) dx
L −L n π −L
" x=L Z L #
L
= − 2 2 f (x)c′n (x) − f ′ (x)c′n (x) dx .
nπ x=−L −L
′
x=L
Como c′n (x) = − nπ
L
s n (x), resulta que c ′
n (−L) = c ′
n (L) = 0, entonces f (x)c n (x)
x=−L
=0
y luego "Z #
L
L ′ ′
an = 2 2 f (x)cn (x) dx
nπ −L
" x=L Z L #
L ′
′′
= 2 2 f (x)cn (x) − f (x)cn (x) dx .
nπ x=−L −L
x=L
Como f ′ (−L) = f ′ (L) resulta f ′ (x)cn (x)x=−L = 0 y
Z L
L
an = − 2 2 f ′′ (x)cn (x) dx.
nπ −L
Análogamente,
Z L
L
bn = − 2 2 f ′′ (x)sn (x) dx,
nπ −L
nπx ′′
con sn (x) = sen L
. Como |f (x)| ≤ M tenemos que
160
www.elsolucionario.net
por lo tanto Z L
L 2L2 M
|an | ≤ 2 2 M dx = ,
n π −L n2 π 2
Z L
L 2L2 M
|bn | ≤ 2 2 M dx = ,
n π −L n2 π 2
por lo que finalmente
∞ ∞
X 4L2 M X 1 4L2 M 1
|f (x) − SFN f (x)| ≤ |an | + |bn | ≤ ≤ .
n=N +1
π 2 n=N +1 n2 π2 N
∞ ∞
X 1 1 X 1
Donde hemos usado que 2
≤ , lo que puede verse comparando la suma
n=N +1
n N n=N +1
n2
Z ∞
1
con la integral 2
dx.
N x
www.elsolucionario.net
sirve, además, para saber cuántos términos de la serie de Fourier se deben sumar para
garantizar una aproximación con una tolerancia dada. Veamos un ejemplo.
Ejemplo 9.25. Dada f (x) = x3 −x y L = 1. Determinar N para asegurar que | SFN f (x)−
f (x)| ≤ 0,01. Al ser f un polinomio es de clase C ∞ y por lo tanto también C 2 . Además
f (−1) = f (1) = 0 y f ′ (x) = 3x2 − 1 por lo que f ′ (−1) = f ′ (1) = 2. Luego se cumplen
las hipótesis de f para que valga la afirmación del teorema anterior. Necesitamos saber
M = máx[−1,1] |f ′′ |. Veamos, f ′′ (x) = 6x, y luego M = 6. Queremos saber para qué N ∈ N
se cumple que
4L2 M 1 4×6 1 24
≤ 0,01 ⇐⇒ ≤ 0,01 ⇐⇒ = 243,17 . . . < N.
π2 N π2 N 0,01π 2
Luego, para N ≥ 244, se cumple que |f (x) − SFN f (x)| ≤ 0,01 para todo x ∈ [−1, 1].
Queremos utilizar lo aprendido en la sección anterior para ver si es cierto lo que afirma-
mos en el Capı́tulo 8, que una gran familia de funciones definidas sobre el intervalo [0, L]
se puede aproximar por series de cosenos (para resolver problemas con CB Neumann) y
también por series de senos (para resolver problemas con CB Dirichlet).
Recordemos las siguiente definición.
161
www.elsolucionario.net
Definición 9.26. Una función f (x) se llama par si f (−x) = f (x) para todo x, y una
función f (x) se llama impar si f (−x) = −f (x) para todo x.
La demostración de la siguiente proposición es una consecuencia inmediata del ejer-
cicio 9.4.
Proposición 9.27. Sea f (x) definida en el intervalo [−L, L] con coeficientes de Fourier
Z L Z
1 1 L nπx
a0 = f (x) dx, an = f (x) cos dx, n = 1, 2, . . . ,
2L −L L −L L
Z
1 L nπx
bn = f (x) sen dx, n = 1, 2, . . . .
L −L L
www.elsolucionario.net
Si f (x) es una función impar, entonces
a0 = 0, an = 0, n = 1, 2, . . . ,
Z L
2 nπx
bn = f (x) sen dx, n = 1, 2, . . . .
L 0 L
162
www.elsolucionario.net
Proposición 9.30. Sea f (x) definida en el intervalo [0, L] tal que existen las integrales
de (9.5) y (9.6). Entonces la serie de Fourier de senos de f coincide con la serie de
Fourier de su extensión impar fi y la serie de Fourier de cosenos de f coincide con la
serie de Fourier de su extensión par fp . O sea,
SFs f = SF fi , SFc f = SF fp .
www.elsolucionario.net
Demostración. Veamos que la primera afirmación es cierta, la segunda es análoga. La
serie de Fourier de senos de f es
∞ Z
s
X nπx 2 L nπx
SF f (x) = bn sen , con bn = f (x) sen dx.
n=1
L L 0 L
Por otro lado, si consideramos la extensión impar fi de f , resulta que por la Proposi-
ción 9.27, sus coeficientes de Fourier (de la serie de senos y cosenos) son
a0 = 0, an = 0, n = 1, 2, . . . ,
Z L Z L
2 nπx 2 nπx
bn = fi (x) sen dx = f (x) sen dx, n = 1, 2, . . . .
L 0 L L 0 L
Y por lo tanto, la serie de Fourier de fi es
∞ Z L
Xnπx 2 nπx
SF fi = bn sen , con bn = f (x) sen dx,
n=1
L L 0 L
es decir SFs f = SF fi .
Con lo visto en la última proposición podemos adaptar los resultados de convergencia
de series de Fourier (de senos y cosenos) para obtener resultados de convergencia de las
series de Fourier de senos y de las de cosenos.
Teorema 9.31 (Convergencia de la serie de Fourier de senos). Sea f una función C 1 a
trozos en [0, L]. Entonces las sumas parciales SFsN f de la serie de Fourier de senos de f
convergen puntualmente (en cada x) a
f (x+ ) + f (x− )
, si 0 < x < L,
SFs f (x) = 2
0, si x = 0 ó x = L.
163
www.elsolucionario.net
Demostración. Para demostrar las dos afirmaciones de este teorema basta ver que si f (x)
cumple las hipótesis del mismo, entonces su extensión impar fi cumple las hipótesis de
los Teoremas 9.18 y 9.21, respectivamente.
www.elsolucionario.net
Si f (x) es además continua en [0, L], entonces la suma parcial SFcN f converge uniforme-
mente (no hay fenómeno de Gibbs) a f (x) en [0, L], es decir,
Demostración. Para demostrar las dos afirmaciones de este teorema basta ver que si f (x)
cumple las hipótesis del mismo, entonces su extensión impar fp cumple las hipótesis de
los Teoremas 9.18 y 9.21, respectivamente.
También se puede aplicar el Teorema 9.23 a las extensiones par e impar de una función
f definida en [0, L] para obtener los siguientes teoremas.
Teorema 9.33 (Convergencia con orden de la serie de Fourier de senos). Sea f (x) una
función de clase C 2 en el intervalo [0, L], tal que f (0) = 0 y f (L) = 0, y sea M = máx |f ′′ |.
[0,L]
Entonces, para cada N ≥ 1 se tiene que
4L2 M 1
|f (x) − SFsN f (x)| ≤ , para todo x ∈ [0, L].
π2 N
Demostración. Si f es C 2 en [0, L] y f (0) = f (L) = 0 entonces su extensión impar fi es
C 2 en [−L, L], fi (−L) = fi (L) y fi′ (−L) = fi′ (L) (pues la derivada de una función impar
es una función par). Aplicamos entonces el Teorema 9.23 a fi y obtenemos que
4L2 M 1
|fi (x) − SFN fi (x)| ≤ , para todo x ∈ [−L, L].
π2 N
Por un lado fi (x) = f (x) en [0, L], y por la Proposición 9.30 SFN fi (x) = SFsN f (x) para
4L2 M 1
todo x ∈ [0, L] y todo N ∈ N. Por lo tanto |f (x) − SFsN f (x)| ≤ para todo
π2 N
x ∈ [0, L].
164
www.elsolucionario.net
Teorema 9.34 (Convergencia con orden de la serie de Fourier de cosenos). Sea f (x)
una función de clase C 2 en el intervalo [0, L], tal que f ′ (0) = 0 y f ′ (L) = 0, y sea
M = máx |f ′′ |. Entonces, para cada N ≥ 1 se tiene que
[0,L]
4L2 M 1
|f (x) − SFcN f (x)| ≤ , para todo x ∈ [0, L].
π2 N
Demostración. Si f es C 2 en [0, L] y f ′ (0) = f ′ (L) = 0 entonces su extensión par fp es
C 2 en [−L, L], fp (−L) = fp (L) y fp′ (−L) = fp′ (L). Aplicamos entonces el Teorema 9.23 a
fp y obtenemos que
4L2 M 1
|fp (x) − SFN fp (x)| ≤ , para todo x ∈ [−L, L].
π2 N
Por un lado fp (x) = f (x) en [0, L], y por la Proposición 9.30 SFN fp (x) = SFcN f (x) para
4L2 M 1
todo x ∈ [0, L] y todo N ∈ N. Por lo tanto |f (x) − SFcN f (x)| ≤ para todo
π2 N
x ∈ [0, L].
www.elsolucionario.net
Los últimos teoremas son de aplicación directa para el cálculo aproximado de solucio-
nes de problemas de difusión. Veamos un ejemplo.
Ejemplo 9.35. Considerar el problema
ut = uxx ,
0 < x < 1, t > 0,
u(0, t) = 0, u(1, t) = 0, t > 0.
u(x, 0) = x3 − x, 0 < x < 1,
Hallar una solución aproximada uN (x, t) que cumple |uN (x, t) − u(x, t)| ≤ 0,01 para todo
x ∈ [0, L] y todo t ≥ 0. Las soluciones exactas que conocemos de este problema son de la
forma
N
X 2 2
uN (x, t) = bn sen(nπx)e−n π t .
n=1
Por el principio del máximo, dado que u(x, t) = uN (x, t) = 0 para x = 0 y x = L (t > 0),
sabemos que
máx |u(x, t) − uN (x, t)| ≤ máx |u(x, 0) − uN (x, 0)|.
0≤x≤1 0≤x≤1
t≥0
165
www.elsolucionario.net
Z 1
12 12
con bn = 2 (x3 − x) sen(nπx) dx = 3 3
cos(nπ) = 3 3 (−1)n . De esta manera, la
0 nπ nπ
solución buscada resulta
244
X (−1)n 2 π2 t
uN (x, t) = 12 sen(nπx)e−n .
n=1
n3 π 3
9.5. Ejercicios
(
1 si x > 0
9.1. Para L = 1 y la función f (x) = graficar con la computadora SFN f (x)
0 si x ≤ 0
en [−1, 1], para N = 1, 3, 5, 7, 9, . . . , 101, y verificar que siempre máx SFN f (x) ≥ 1,089
x∈[0,1]
(mientras que máxx∈R f (x) = 1). Para hacerlo, tomar x=[-1:.001:1] y evaluar SFN f
en esa grilla, luego tomar el máximo valor. Este fenómeno, observado experimentalmente
en la computadora, puede demostrarse para todo N ∈ N. Nos dice que la discontinuidad
de f en x = 0 causa un overshooting debido al fenómeno de Gibbs mayor al 8 % del salto
www.elsolucionario.net
en la discontinuidad.
9.2. Para cada una de las siguientes funciones, definir y graficar (a mano) la función
arreglada en [−L, L], con L = 1 y graficar su extensión periódica en [−3L, 3L].
* 9.3. Sea f (x) una función C 3 en el intervalo [−L, L] tal que f (−L) = f (L), f ′ (−L) =
f ′ (L), y f ′′ (−L) = f ′′ (L). Demostrar que si M = máx |f ′′′ (x)|, entonces
x∈[−L,L]
4L3 M 1
(a) máx |f (x) − SFN f (x)| ≤ ,
x∈[−L,L] π3 N 2
4L2 M 1
(c) máx |f ′ (x) − SFN f ′ (x)| ≤ .
x∈[−L,L] π2 N
166
www.elsolucionario.net
9.5. Mostrar que cualquier función f (x) definida en [−L, L] se puede escribir como suma
f (x) + f (−x) f (x) − f (−x)
de una función par y una impar. Ayuda: f (x) = + .
2 2
9.6. Para cada una de las siguientes funciones definidas en [0, L], definir y graficar (a
mano) su extensión par y su extensión impar a [−L, L], la función arreglada de la extensión
par y de la impar en [−L, L] y graficar su extensión periódica en [−3L, 3L] (considerar
L = 1).
9.7. Consideremos L = 1 y la función f (x) definida sobre el intervalo [0, 1] cuyo gráfico
es el de la figura.
www.elsolucionario.net
(a) 1 (b) 1 (c) 1
SFs10 f (x), SFs f (x). Es decir, la suma de diez términos de la serie de Fourier de
senos de f y el lı́mite de las sumas parciales, o la serie total de senos.
SFc10 f (x), SFc f (x). Es decir, la suma de diez términos de la serie de Fourier de
cosenos de f y el lı́mite de las sumas parciales, o la serie total de cosenos.
Usando el principio del máximo determinar una solución aproximada uN (x, t) que cumpla
167
www.elsolucionario.net
(a) Usando el principio del máximo determinar una solución aproximada uN (x, t) que
cumpla
|u(x, t) − uN (x, t)| ≤ 0,005, 0 ≤ x ≤ L, t > 0.
(b) Usando (la computadora y) la solución aproximada del item anterior, determinar
(aproximadamente) el valor de t0 a partir del cual u(x, t) ≤ 2 para todo t ≥ t0 .
www.elsolucionario.net
(a) Usando el principio del máximo (enunciado en el Problema 8.17) y el Teorema 9.23
determinar una solución aproximada uN (x, t) que cumpla
(b) Usando (la computadora y) la solución aproximada del item anterior, determinar
(aproximadamente) el valor de t0 a partir del cual u(x, t) ≥ 0,5 para todo t ≥ t0 .
Bibliografı́a complementaria
[Bleecker-Csordas 1996] Bleecker, D., Csordas, G., Basic Partial Differential Equations,
International Press, Cambridge, Massachusetts, 1996.
[Haberman 1998] Haberman, R., Elementary Applied Partial Differential Equations, Pren-
tice Hall, Upper Saddle River, NJ, 1998.
168
www.elsolucionario.net
Capı́tulo 10
La ecuación de Laplace
www.elsolucionario.net
de Laplace resulta uxx = 0 cuyas soluciones son las funciones lineales u(x) = c1 x + c2 .
Imponer condiciones de borde es sencillo en esta situación, y salvo el caso en que ambos
extremos tengan condiciones de borde de flujo prescripto, pueden hallarse siempre las
constantes c1 y c2 para que las condiciones de borde se cumplan, ver Sección 8.5. En
este capı́tulo estudiaremos cómo resolver la ecuación de Laplace en diferentes regiones del
plano y del espacio, y también veremos algunas propiedades cualitativas de las soluciones,
como el principio del máximo y la propiedad del valor medio.
∂ 2u ∂ 2u
EDP: + =0 0 < x < L, 0 < y < H.
∂x2 ∂y 2
donde g1 (y), g2 (y), f1 (x) y f2 (x) son funciones dadas. Otras condiciones de borde serán
presentadas en los ejercicios. En este caso, la EDP es lineal y homogénea pero las condi-
ciones de borde, aunque lineales, no son homogéneas. No podremos aplicar el método de
169
www.elsolucionario.net
separación de variables a este problema tal como está. La razón de esto es que al separar
variables, el problema a valores de borde que determina las constantes de separación po-
sibles o autovalores debe tener condiciones de borde homogéneas. En este ejemplo, todas
las condiciones de borde son no-homogéneas.
Podemos superar esta dificultad dándonos cuenta de que el problema original es no-
homogéneo debido a las cuatro condiciones de borde no-homogéneas. Utilizaremos el
principio de superposición para desglosar nuestro problema en cuatro subproblemas, cada
uno de los cuales tiene sólo una condición no-homogénea. Escribimos
donde cada una de las funciones ui (x, y) (i = 1, 2, 3, 4) satisface la ecuación de Laplace con
una condición de borde no-homogénea y las otras tres condiciones de borde homogéneas,
como se muestra en la Figura 10.1.
u = f2 (x) u1 = 0 u2 = 0 u3 = f2 (x) u4 = 0
www.elsolucionario.net
u = g2 (y) u1 = 0 u2 = g2 (y) u3 = 0 u4 = 0
∇2 u = 0 = ∇2 u 1 = 0 + ∇2 u 2 = 0 + ∇2 u 3 = 0 + ∇2 u 4 = 0
u = g1 (y) u1 = 0 u2 = 0 u3 = 0 u4 = g1 (y)
u = f1 (x) u1 = f1 (x) u2 = 0 u3 = 0 u4 = 0
El método para hallar cualquiera de las ui (x, y) es el mismo; sólo algunos detalles
difieren. Haremos la resolución para encontrar u4 (x, y), y dejaremos los otros casos para
los ejercicios:
∂ 2 u4 ∂ 2 u4
EDP: + =0 0 < x < L, 0<y<H
∂x2 ∂y 2
CB1: u4 (0, y) = g1 (y) 0<y<H
CB2: u4 (L, y) = 0 0<y<H
CB3: u4 (x, 0) = 0 0<x<L
CB4: u4 (x, H) = 0 0 < x < L.
170
www.elsolucionario.net
Nos proponemos resolver este problema por el método de separación de variables. Comen-
zamos ignorando la condición no-homogénea u4 (0, y) = g1 (y). Eventualmente sumaremos
soluciones producto para obtener esta condición, pero por ahora la ignoraremos. Busca-
mos soluciones producto
u4 (x, y) = h(x)φ(y),
que satisfagan las condiciones de borde homogéneas CB2, CB3 y CB4. Éstas nos dicen
que
h(L) = 0, φ(0) = 0, φ(H) = 0.
Luego, la parte de la solución dependiente de y, φ(y) tiene dos condiciones de borde
homogéneas, mientras que la parte dependiente de x sólo tiene una. Si sustituimos la
solución producto en la ecuación de Laplace obtenemos
www.elsolucionario.net
h(x) φ(y)
El lado izquierdo es sólo función de la variable x, mientras que el lado derecho depende
sólo de y. Por lo tanto, ambos deben ser iguales a una constante de separación. ¿Qué usa-
mos? ¿λ ó −λ? Una de ellas es más conveniente. Si la constante de separación fuera
negativa (como lo era antes), la ecuación (10.1) implicarı́a que h(x) oscila (trigonométri-
ca) y φ(y) es una combinación de exponenciales. Pero las condiciones de borde para φ(y)
indican que este caso no es posible. Por otro lado, si la constante de separación fuera po-
sitiva, (10.1) implicarı́a que h(x) es combinación de exponenciales y φ(y) es combinación
de senos y cosenos. Esto parece más razonable dadas las condiciones de borde y entonces
introducimos la variable de separación λ (aunque no suponemos λ ≥ 0):
1 ′′ 1 ′′
h (x) = − φ (y) = λ.
h(x) φ(y)
Obtenemos entonces dos ecuaciones diferenciales ordinarias:
171
www.elsolucionario.net
L. Todos los autovalores son positivos λ > 0. Las autofunciones son senos pues φ(0) = 0.
Más aún, la condición φ(H) = 0 implica que
nπ 2 nπy
λ= , φ(y) = sen , n = 1, 2, 3, . . . .
H H
Para obtener soluciones producto, debemos ahora resolver la EDO para h(x) con la con-
dición de borde homogénea correspondiente y con λ un autovalor de los ya encontrados.
Como λ = (nπ/H)2 , la ecuación resulta
nπ 2
h′′ (x) = h(x). (10.2)
H
La solución general es una combinación lineal de exponenciales, que también puede ex-
presarse como una combinación lineal de funciones hiperbólicas. Podrı́amos utilizar cual-
quiera de estas dos expresiones de la solución general, pero ninguna de estas es la más
apropiada para la condición de borde h(L) = 0. Podemos encontrar la solución que bus-
camos de un modo más rápido y sencillo si observamos que las funciones trasladadas
www.elsolucionario.net
nπ(x − L) nπ(x − L)
cosh y senh
H H
son soluciones linealmente independientes de (10.2). Por lo que la solución general también
puede escribirse como (ver Ejercicios 6.8–6.9)
nπ(x − L) nπ(x − L)
h(x) = a1 cosh + a2 senh .
H H
Resulta natural preguntarse ahora por qué es más conveniente esta expresión de la solu-
ción general. Bien, si queremos forzar la condición de borde h(L) = 0 obtenemos
nπ(L − L) nπ(L − L)
0 = h(L) = a1 cosh + a2 senh = a1 cosh 0 + a2 senh 0 = a1 .
H H
Una vez que tenemos esta solución general, imponemos la condición de borde no-homogénea
u4 (0, y) = g1 (y):
∞
X nπy nπ(−L)
u4 (0, y) = bn sen senh = g1 (y).
n=1
H H
172
www.elsolucionario.net
Este es el mismo tipo de serie de senos que ya hemos discutido con una sola diferencia:
Antes tenı́amos los coeficientes bn multiplicando a las funciones seno, pero ahora tenemos
que bn senh nπ(−L)
H
multiplica a las funciones seno (notar que bn senh nπ(−L)
H
no depende
nπ(−L)
de la variable y). Por ende, si llamamos Bn = bn senh H nuestro problema consiste
en hallar Bn tal que
∞
X nπy
Bn sen = g1 (y).
n=1
H
Luego, por lo que ya hemos visto en secciones anteriores
Z
2 H nπy
Bn = g1 (y) sen dy,
H 0 H
y despejando bn obtenemos
Z H
Bn 2 nπy
bn = nπ(−L)
= nπ(−L)
g1 (y) sen dy.
senh H senh 0 H
H H
www.elsolucionario.net
n.
Por lo tanto, la solución u4 (x, y) es la serie
∞
X nπy nπ(x − L)
u4 (x, y) = bn sen senh .
n=1
H H
con sus coeficientes bn dados por
Z H
2 nπy
bn = nπ(−L)
g1 (y) sen dy.
H senh 0 H
H
173
www.elsolucionario.net
u(a, θ) = f (θ)
∂ 2 u 1 ∂u 1 ∂ 2u
EDP: + + =0 0<r<a −π <θ <π
∂r2 r ∂r r2 ∂θ2
CB: u(a, θ) = f (θ), −π < θ < π.
www.elsolucionario.net
A primera vista parecerı́a que no podemos utilizar el método de separación de variables
porque no tenemos ninguna condición de borde homogénea. Sin embargo, el utilizar coor-
denadas polares requiere una discusión acerca del método de separación de variables que
nos ayudará a comprenderlo un poco mejor. Si resolvemos la ecuación de Laplace en un
rectángulo, entonces se necesitan condiciones en los extremos de definición de las varia-
bles, x = 0, x = L, y = 0, y = H. Estos extremos coinciden con los bordes fı́sicos de la
placa rectangular que se considera. Por otro lado, para coordenadas polares, 0 ≤ r ≤ a,
−π ≤ θ ≤ π. Matemáticamente, se necesitan dos condiciones en los extremos r = 0,
r = a, y θ = −π, θ = π. La única condición que corresponde a una frontera fı́sica es en
r = a. Por lo tanto todavı́a nos hace falta determinar condiciones de borde en r = 0 y
en θ = ±π. Para hallarlas utilizaremos consideraciones del problema fı́sico. Las coorde-
nadas polares son singulares en r = 0; por razones fı́sicas, pediremos que la incógnita u
(temperatura, densidad) sea finita, o acotada allı́:
Aún nos falta definir condiciones en θ = ±π. Esto es similar a lo que ocurre en la situación
del anillo circular. El extremo θ = −π corresponde al mismo punto de la placa circular que
θ = π. Aunque allı́ no hay realmente un borde de la placa, el hecho de que la temperatura
o densidad sea continua allı́ y que el flujo de calor o de masa en la dirección de θ también
sea continuo, implica las siguientes condiciones de borde:
∂u ∂u
u(r, −π) = u(r, π) (r, −π) = (r, π).
∂θ ∂θ
Notamos que estas condiciones “de borde” son también lineales y homogéneas. De este
modo, el problema luce similar a la ecuación de Laplace en un rectángulo: Hay cuatro
condiciones, de las cuales sólo una es no-homogénea: u(a, θ) = f (θ). Aplicaremos entonces
el método de separación de variables.
174
www.elsolucionario.net
u(r, θ) = φ(θ)G(r),
que satisfagan la EDP y las tres condiciones de borde homogéneas (nuevamente ignoramos
momentáneamente la condición no-homogénea). Las condiciones periódicas implican
www.elsolucionario.net
Hemos introducido la constante de separación como λ (en lugar de −λ) porque hay dos
condiciones homogéneas para la función φ de la variable θ. El problema que determi-
nará los autovalores λ es
Los autovalores λ se determinan del modo usual. En efecto, este es uno de los problemas
estándar que ya hemos resuelto, el del alambre circular, con L = π. Por lo tanto los
autovalores son nπ 2
λ= = n2 ,
L
con las autofunciones correspondientes
ó
r2 G′′ (r) + rG′ (r) − n2 G(r) = 0. (10.3)
La ecuación (10.3) es lineal y homogénea pero tiene coeficientes no-constantes. Sin
embargo se puede resolver fácilmente. La forma más sencilla de resolver esta ecuación,
es dándose cuenta de que para el operador diferencial en (10.3) la función G(r) = rp se
reproduce a sı́ misma. Esto quiere decir que si reemplazamos rp en el lugar de G(r) en la
ecuación (10.3) obtenemos un múltiplo de rp . En efecto
d2 p d
r2 2
(r ) + r (rp ) − n2 rp = r2 p (p − 1) rp−2 + r p rp−1 − n2 rp
dr dr
= p (p − 1) + p − n2 rp .
175
www.elsolucionario.net
Luego nuestro problema se transforma en encontrar las potencias p para las cuales
p(p − 1) + p − n2 = 0.
Pero
p(p − 1) + p − n2 = p2 − p + p − n2 = p2 − n2
y esta última expresión es igual a cero cuando p = n ó p = −n. Si n 6= 0 estas raı́ces son
distintas y las funciones rn y r−n son linealmente independientes, por lo que la solución
general resulta
G(r) = c1 rn + c2 r−n .
El caso n = 0 es diferente, e igualmente importante porque λ = 0 también es un autovalor.
Si n = 0 la ecuación (10.3) resulta
r2 G′′ (r) + rG′ (r) = 0
que es equivalente a
rG′′ (r) + G′ (r) = 0
y también a
www.elsolucionario.net
d
(rG′ (r)) = 0
dr
lo que a su vez equivale a decir que rG′ (r) es constante. Luego G′ (r) = const/r, integrando
obtenemos
G(r) = c1 + c2 ln r.
Ya hemos discutido la condición en r = 0 y hemos dicho que requeriremos |u(0, θ)| <
∞, lo que implica para la solución producto que
|G(0)| < ∞.
Si ahora observamos las soluciones generales obtenidas, vemos que la parte multiplicada
por c2 tiende a ∞ cuando r → 0. Esto nos dice que c2 = 0 y luego la solución general de
la parte dependiente de r es
G(r) = c1 rn , n ≥ 0,
donde para n = 0 esta expresión se reduce a la constante c1 .
Concluyendo, las soluciones producto que satisfacen las condiciones homogéneas son
rn cos nθ (n ≥ 1) y rn sen nθ (n ≥ 0).
Por el principio de superposición, la solución general de la ecuación de Laplace en el
cı́rculo es ∞ ∞
X X
n
u(r, θ) = A0 + An r cos nθ + Bn rn sen nθ.
n=1 n=1
176
www.elsolucionario.net
Este problema es similar al que resolvimos en la Sección 8.3. La única diferencia es que
aquı́ el papel de los coeficientes an y bn es jugado por An an y por Bn an , respectivamente.
Usando las mismas fórmulas integrales de esa sección, obtenemos
Z π
1
A0 = f (θ) dθ
2π −π
Z
n 1 π
An a = f (θ) cos nθ dθ, n≥1
π −π
Z
n 1 π
Bn a = f (θ) sen nθ dθ, n ≥ 1.
π −π
www.elsolucionario.net
donde los coeficientes An y Bn están dados por
Z π
1
A0 = f (θ) dθ
2π −π
Z π
1
An = n f (θ) cos nθ dθ, n≥1
πa −π
Z π
1
Bn = n f (θ) sen nθ dθ, n ≥ 1.
πa −π
177
www.elsolucionario.net
Esto nos dice que si u es solución de la ecuación de Laplace (ecuación de difusión estacio-
naria), el flujo de calor neto a través de la frontera debe ser cero. Esto también está claro
desde el punto de vista fı́sico si pensamos en difusión de calor. Si el flujo de calor neto
a través de la frontera no fuera cero, habrı́a cambio (en el tiempo) de la energı́a térmica
total dentro de la región, violando la condición de estado estacionario. La condición (10.4)
se denomina condición de compatibilidad para la ecuación de Laplace.
2
Ejemplo 10.1. Consideremos la ecuación de Laplace ∇ u = 0 en el rectángulo [0, L] ×
[0, H] con las siguientes condiciones de borde
ux (0, y) = f (y), ux (L, y) = 0, uy (x, 0) = 0, uy (x, H) = g(x),
para 0 < x < L, 0 < y < H. Es decir, el flujo está prescripto en el lado izquierdo (x = 0)
y en el lado superior (y = H). La condición de compatibilidad (10.4) nos dice que
Z
0= ∇u · n ds
∂Ω
Z L Z H
= ∇u(x, 0) · (−j) dx + ∇u(L, y) · i dy
0 0
www.elsolucionario.net
Z L Z H
+ ∇u(x, H) · j dx + ∇u(0, y) · (−i) dy
0 0
Z L Z H Z L Z H
= −uy (x, 0) dx + ux (L, y) dy + uy (x, H) dx + −ux (0, y) dy.
0 0 0 0
178
www.elsolucionario.net
www.elsolucionario.net
R
(x0 , y0 )
Figura 10.3: Cı́rculo dentro de una región
general. La circunferencia no toca el borde
de la región.
179
www.elsolucionario.net
Entonces,
Z 2π
d 1 ∂
V (r) = u(x0 + r cos θ, y0 + r sen θ) dθ
dr 2π 0 ∂r
Z 2π
1
= ux (x0 + r cos θ, y0 + r sen θ) cos θ + uy (x0 + r cos θ, y0 + r sen θ) sen θ dθ
2π 0
Z 2π
1
= ∇u(x0 + r cos θ, y0 + r sen θ) · (cos θi + sen θj ) dθ
2π 0 | {z }
n
Z
1
= ∇u · n ds = 0,
2πr Cr
por (10.5). Por lo tanto V (r) es constante para 0 < r ≤ R y resulta que
Z Z
1 1
V (R) = u ds = u ds
2πR CR longitud(CR ) CR
Z Z
1 1
= lı́m+ V (r) = lı́m+ u ds = lı́m+ u ds.
www.elsolucionario.net
r→0 r→0 2πr Cr r→0 longitud(Cr ) Cr
Por otro lado, como los promedios sobre circunferencias que se encogen a un punto con-
vergen al valor de la función en ese punto
Z
1
lı́m u ds = u(x0 , y0 ),
r→0+ longitud(Cr ) Cr
y por lo tanto Z
1
u ds = u(x0 , y0 ).
2πR CR
Teorema 10.4 (Principio del máximo fuerte). Sea u una función C 2 y armónica en un
dominio conexo Ω. Si alcanza su máximo en un punto interior (que no está en el borde),
entonces u es constante.
180
www.elsolucionario.net
Demostración. Sea M = máxΩ u y sea (x0 , y0 ) ∈ Ω tal que u(x0 , y0 ) = M . Sea ahora R
tal que el disco DR con centro (x0 , y0 ) y radio R está completamente contenido en Ω,
entonces Z
1
u(x0 , y0 ) = M = u ds = promedio de u sobre CR .
2πR CR
Si existiera un punto sobre CR donde u < M , entonces para que el promedio sea M ,
deberı́a haber otros puntos sobre CR donde u > M , pero esto es imposible pues M es el
máximo de u. Por lo tanto u = M en todos los puntos de CR y esto se cumple para todo
radio R tal que DR ⊂ Ω. En otras palabras, u = M en todo disco centrado en (x0 , y0 )
contenido en Ω.
Ahora, si tomamos (x1 , y1 ) contenido en un disco de estos, resulta que también
u(x1 , y1 ) = M y por lo tanto, u = M en todo disco contenido centrado en (x1 , y1 ) y
contenido en Ω. Repitiendo este argumento cubrimos Ω con discos donde u = M y resul-
ta que u ≡ M en Ω.
Observación 10.5. La prueba matemática rigurosa del principio del máximo fuerte para
funciones armónicas puede hacerse demostrando que el conjunto de puntos de Ω donde
u = M es abierto y cerrado en Ω, y no vacı́o. Usando que Ω es conexo resulta que el
www.elsolucionario.net
conjunto es todo Ω.
Como consecuencia inmediata del principio del máximo fuerte obtenemos el principio
del máximo débil.
Teorema 10.6 (Principio del máximo débil). Sea Ω una región acotada de R2 . Si u es
2
una función C 2 en Ω y armónica (∇ u = 0) en Ω y es continua en Ω := Ω ∪ ∂Ω entonces
máx u = máx u,
Ω ∂Ω
181
www.elsolucionario.net
www.elsolucionario.net
Entonces
máx |u1 − u2 | ≤ máx |g1 − g2 |
Ω ∂Ω
Antes de proceder a la demostración vale la pena interpretar lo que este principio del
máximo significa en la práctica. Supongamos que g1 es el valor medido de la temperatura
en el borde de una región y g2 es el valor real o exacto (g1 es una aproximación de g2 ).
Sea ahora u1 la solución de la ecuación de Poisson con dato g1 en el borde. Si u2 denota
la temperatura real en la región, entonces u2 es la solución de la ecuación de Laplace
con dato g2 en el borde. Pregunta: ¿Cuán bien aproxima u1 (solución calculada) a u2
(temperatura real desconocida)? Respuesta: tan bien como g1 aproxime a g2 .
Demostración del Teorema 10.9. Definamos (como en la demostración de la proposición
anterior) v = u1 − u2 . Entonces v es solución de la ecuación de Laplace con dato de borde
g1 − g2 . Más precisamente
2
∇ v = 0 en Ω, v = g1 − g2 en ∂Ω.
Por el principio del máximo
máx v = máx v = máx g1 − g2 ≤ máx |g1 − g2 |,
Ω ∂Ω ∂Ω ∂Ω
y por lo tanto
v(x) ≤ máx |g1 − g2 |, para todo x ∈ Ω.
∂Ω
Análogamente,
−v(x) ≤ máx |g1 − g2 |, para todo x ∈ Ω.
∂Ω
Luego
− máx |g1 − g2 | ≤ v(x) = u1 (x) − u2 (x) ≤ máx |g1 − g2 |, para todo x ∈ Ω.
∂Ω ∂Ω
182
www.elsolucionario.net
10.5. Ejercicios
10.1. Hallar la solución de la ecuación de Laplace dentro de un rectángulo 0 ≤ x ≤ L,
0 ≤ y ≤ H con las siguientes condiciones de borde (para 0 < x < L, 0 < y < H):
∂u ∂u
(a) (0, y) = 0, (L, y) = 0, u(x, 0) = 0, u(x, H) = f (x).
∂x ∂x
∂u ∂u
(b) (0, y) = g(y), (L, y) = 0, u(x, 0) = 0, u(x, H) = 0.
∂x ∂x
10.2. Supongamos que u(x, y) es la solución de la ecuación de Laplace dentro del rectángu-
lo 0 ≤ x ≤ L, 0 ≤ y ≤ H con las condiciones de borde:
www.elsolucionario.net
sólo funciona bajo la condición de la parte (a). Notar que queda una constante sin
determinar, es decir que puede tomar cualquier valor.
(c) Considerar la ecuación del calor dependiente del tiempo vt = k(vxx + vyy ) con condi-
ción inicial v(x, y, 0) = g(x, y) y condiciones de borde
para 0 < x < L, 0 < y < H, t > 0, (las mismas que u). Demostrar que si v es una
solución C 2 , entonces la energı́a
Z LZ H
E(t) = v(x, y, t) dy dx
0 0
RLRH
se mantiene constantemente igual a 0 0
g(x, y) dy dx para todo t > 0.
(d) Como dijimos antes, la solución del ı́tem (b) tiene una constante arbitraria que queda
sin determinar. Determinarla considerando que u(x, y) es el estado estacionario de
la ecuación (del calor) dependiente del tiempo del ı́tem (c). Es decir, suponer que
(Ayuda: la constante a determinar tiene que ver con la condición inicial g(x, y).)
10.3. Resolver la ecuación de Laplace dentro de un cuarto de cı́rculo de radio 1 (0 < θ <
π/2, 0 < r < 1), sujeto a las condiciones de borde
∂u
(a) (r, 0) = 0, u(r, π/2) = 0, u(1, θ) = f (θ)
∂θ
∂u ∂u ∂u
(b) (r, 0) = 0, (r, π/2) = 0, (1, θ) = g(θ).
∂θ ∂θ ∂r
R π/2
Mostrar que esta última solución existe sólo si 0 g(θ) dθ = 0. ¿Hay solución única?
183
www.elsolucionario.net
www.elsolucionario.net
* 10.7. Demostrar rigurosamente el principio del máximo fuerte del Teorema 10.4 si-
guiendo la sugerencia de la Observación 10.5.
* 10.8. Demostrar con un contraejemplo que si Ω no es acotado, entonces no vale el
principio del máximo débil (Teorema 10.6). Más precisamente, hallar un dominio Ω y
una función armónica que valga cero en ∂Ω y resulte no acotada en Ω.
Bibliografı́a complementaria
[Haberman 1998] Haberman, R., Elementary Applied Partial Differential Equations, Pren-
tice Hall, Upper Saddle River, NJ, 1998.
[Pinchover-Rubinstein 2005] Pinchover, Y., Rubinstein, J., An Introduction to Partial Dif-
ferential Equations, Cambridge Univbersity Press, New Yotk, 2005.
184
www.elsolucionario.net
Capı́tulo 11
Dada una cuerda elástica de densidad lineal constante ρ a tensión constante T0 , ali-
neada con el eje x, el desplazamiento vertical u(x, t) del punto a distancia x del origen,
satisface la ecuación
www.elsolucionario.net
ρ utt (x, t) = T0 uxx (x, t).
Esta ecuación se conoce como ecuación de ondas unidimensional, y se obtiene después
de algunas simplificaciones que resultan de suponer pequeños desplazamientos y que
no hay desplazamiento horizontal. Los detalles pueden encontrarse en [Haberman 1998,
T0
Weinberger 1995]. Usualmente se reemplaza por c2 con c > 0 y la ecuación resulta:
ρ
Dado que la ecuación es de segundo orden en el tiempo, deben imponerse dos condi-
ciones iniciales: la posición inicial y la velocidad inicial:
185
www.elsolucionario.net
ux = v w w x + v z z x
ux = v w 1 + v z 1
ux = v w + v z
uxx = vww wx + vwz zx + vzw wx + vzz zx
uxx = vww 1 + vwz 1 + vzw 1 + vzz 1
uxx = vww + 2vwz + vzz .
ut = v w w t + v z z t
www.elsolucionario.net
ut = vw c + vz (−c)
ut = c v w − c v z
utt = c vww wt + c vwz zt − c vzw wt − c vzz zt
utt = c vww c + c vwz (−c) − c vzw c − c vzz (−c)
utt = c2 vww − c2 vwz − c2 vzw + c2 vzz
utt = c2 vww − 2 c2 vwz + c2 vzz .
vwz = 0.
vw (w, z) = C1 (w),
con C1 (w) una función de la variable w. Integrando esta última expresión respecto de w
obtenemos
v(w, z) = C2 (z) + C1 (w),
donde C1 es una primitiva o antiderivada de C1 . En sı́ntesis, v(w, z) es la suma de dos
funciones, una que depende de w = x + c t y otra que depende de z = x − c t. Podemos
decir entonces que
u(x, t) = F (x + c t) + G(x − c t).
Imponemos ahora las condiciones iniciales: La primera u(x, 0) = f (x) implica que
186
www.elsolucionario.net
Por lo tanto, F (x) y G(x) deben cumplir el siguiente sistema de dos ecuaciones con dos
incógnitas (f (x) y g(x) son datos conocidos):
www.elsolucionario.net
2 c 0
Z x+c t
f (x + c t) + f (x − c t) 1
u(x, t) = + g(s) ds..
2 2c x−c t
187
www.elsolucionario.net
(x0 , t0 )
ct
x=
b+
a−
x=
ct
x0 − c t 0 x0 + c t 0 a b
Figura 11.1: Dominio de dependencia y de influencia para la ecuación de ondas.
Por otro lado, los valores de los datos iniciales f (x) y g(x) sobre el intervalo [a, b] sólo
www.elsolucionario.net
influirán sobre los puntos (x, t) que estén sobre el cono de influencia:
que es equivalente a
1 h′′ (t) φ′′ (x)
= = −λ,
c2 h(t) φ(x)
y se cumple siempre que ambos lados sean iguales a la misma constante −λ. Las condicio-
nes de borde u(0, t) = 0 y u(L, t) = 0 implican que φ(0) = 0, y φ(L) = 0 pues queremos
soluciones no triviales.
188
www.elsolucionario.net
www.elsolucionario.net
hn (t) = an cos + bn sen .
L L
Las soluciones de variables separables son entonces
nπx nπct nπct
un (x, t) = sen an cos + bn sen , n = 1, 2, 3 . . . .
L L L
Superponiendo soluciones obtenemos la solución más general
∞
X nπx nπct nπct
u(x, t) = sen an cos + bn sen .
n=1
L L L
Observamos ahora que a diferencia de la ecuación de difusión, las soluciones son ahora
oscilantes, y dados los coeficientes an , bn , nunca existirá lı́mt→∞ u(x, t).
Al imponer las condiciones iniciales quedarán determinados los coeficientes an , bn ,
dependiendo del desplazamiento inicial f (x) y de la velocidad inicial g(x). Por un lado,
∞ X ∞
X nπx nπc0 nπc0 nπx
u(x, 0) = sen an cos +bn sen = an sen .
n=1
L | {z L } | {zL } n=1
L
1 0
189
www.elsolucionario.net
y entonces
∞ ∞
X nπx nπc nπc0 nπc nπc0 X nπc nπx
ut (x, 0) = sen − an sen +bn cos = bn sen .
n=1
L L | {zL } L | {zL } n=1
L L
0 1
www.elsolucionario.net
La solución de (11.2) es
∞
X nπx nπct nπct
u(x, t) = sen an cos + bn sen .
n=1
L L L
con
Z L Z L
2 nπx 2 nπx
an = f (x) sen , bn = g(x) sen dx, n = 1, 2, 3, . . . .
L 0 L nπc 0 L
190
www.elsolucionario.net
11.3. Ejercicios
11.1. Considerar la ecuación de ondas
es solución. Se dice entonces que u(x, t) es una superposición de ondas viajeras. ¿Por
qué?
(b) Hallar fórmulas para F y G para que la solución u(x, t) cumpla con las condiciones
iniciales
2
u(x, 0) = e−x ut (x, 0) = 0, x ∈ R.
www.elsolucionario.net
(d) ¿Qué diferencia se observarı́a en los gráficos si c fuera 6 en lugar de 3?
11.2. Consideremos una cuerda vibrante levemente amortiguada que satisface
Las autofunciones φn (x) y los autovalores λn son los mismos que antes, y las soluciones h deben
satisfacer nπ 2
ρh′′ (t) + βh′ (t) = − T0 h(t).
L
Resolver esta ecuación ordinaria proponiendo h(t) = ert . La condición dada sobre β en el
enunciado implicará que las soluciones r1 , r2 de la ecuación caracterı́stica son complejas.
Bibliografı́a complementaria
[Haberman 1998] Haberman, R., Elementary Applied Partial Differential Equations, Pren-
tice Hall, Upper Saddle River, NJ, 1998.
[Weinberger 1995] Weinberger, H.F., A First Course in Partial Differential Equations:
with Complex Variables and Transform Methods, Dover Publications Inc., 1995.
191
www.elsolucionario.net
Capı́tulo 12
www.elsolucionario.net
La ecuación de difusión en un dominio bidimensional Ω (de forma arbitraria) es
2
(ED) ut (x, y; t) = k∇ u(x, y; t) (x, y) ∈ Ω, t > 0
(CB) u(x, y; t) = 0 (x, y) ∈ ∂Ω, t > 0 (12.1)
(CI) u(x, y; 0) = f (x, y) (x, y) ∈ Ω, (t = 0).
Hemos considerado aquı́ condiciones de borde (CB) de tipo Dirichlet o de valor prescripto,
2
y recordemos que ∇ u(x, y; t) = uxx (x, y; t) + uyy (x, y; t).
Si esta u(x, y; t) es solución de la ecuación de difusión (12.1), entonces debe ocurrir que
2
′ 2 1 h′ (t) ∇ φ(x, y)
h (t) φ(x, y) = k h(t) ∇ φ(x, y) =⇒ = .
k h(t) φ(x, y)
Como los miembros izquierdo y derecho de la última igualdad dependen de diferentes
variables independientes, deben igualar a una misma constante −λ. Es decir
2
h′ (t) = −k λ h(t), ∇ φ(x, y) = −λ φ(x, y).
Para que u(x, y; t) cumpla la condición de borde (CB) y no sea la solución trivial, debe
cumplirse que
192
www.elsolucionario.net
0 < λ1 < λ2 ≤ λ3 ≤ λ4 ≤ · · · → ∞,
y correspondientes autofunciones
www.elsolucionario.net
φ1 , φ2 , φ3 , φ4 , . . . ,
ZZ
2
En el teorema anterior hemos usado la notación kφn k = [φn (x)]2 dA.
Ω
Observación 12.2. En el caso unidimensional Ω = (0, L), las autofunciones φn (x) y los
autovalores λn son
nπ 2 nπx
λn = , φn (x) = sen , n = 1, 2, . . . .
L L
L
En este caso sencillo, sabemos que kφn k2 = , y por eso la fórmula para el cálculo de los
Z2
2 L nπL
coeficientes de Fourier de senos es an = f (x) sen dx.
L 0 x
El teorema garantiza la existencia de autovalores y autofunciones para el Laplaciano.
Además, cualquier función continua y C 1 a trozos se puede aproximar por combinaciones
193
www.elsolucionario.net
www.elsolucionario.net
∞
X
an φn (x, y) = f (x, y), (x, y) ∈ Ω,
n=1
ZZ
1
que acorde al Teorema 12.1 se cumplirá si elegimos an = f (x, y)φn (x, y) dA,
kφn k2 Ω
n = 1, 2, 3, . . . .
Resumiendo:
La solución de (12.1) es
∞
X
u(x, y; t) = an φn (x, y)e−kλn t .
n=1
ZZ
1
con an = f (x, y)φn (x, y) dA, n = 1, 2, 3, . . . .
kφn k2 Ω
194
www.elsolucionario.net
Para que u(x, y; t) cumpla la condición de borde (CB) y no sea la solución trivial, debe
cumplirse que
www.elsolucionario.net
Arribamos entonces al siguiente problema de autovalores para el Laplaciano:
Como los autovalores λn son todos positivos, las soluciones de este problema son
p p
hn (t) = an cos(c λn t) + bn sen(c λn t).
Ası́, aunque no conozcamos los autovalores λn y sus autofunciones φn sabemos que las
soluciones de variables separables son
p p
un (x, y; t) = φn (x, y) an cos(c λn t) + bn sen(c λn t) , n = 1, 2, 3, . . . .
Obtenemos una solución más general sumando las de variables separables, es decir
∞
X p p
u(x, y; t) = φn (x, y) an cos(c λn t) + bn sen(c λn t) .
n=1
195
www.elsolucionario.net
www.elsolucionario.net
ZZ
1 1
que acorde al Teorema 12.1 se cumplirá si elegimos bn = √ g(x, y)φn (x, y) dA,
kφn k2 c λn Ω
n = 1, 2, 3, . . . .
Resumiendo:
La solución de (12.2) es
∞
X p p
u(x, y; t) = φn (x, y) an cos(c λn t) + bn sen(c λn t) ,
n=1
con
ZZ ZZ
1 n 1 1
an = f (x, y)φ (x, y) dA y bn = √ g(x, y)φn (x, y) dA,
kφn k2 Ω
2
kφn k c λn Ω
n = 1, 2, 3, . . . .
196
www.elsolucionario.net
Recordemos que tanto λ como φ(x, y) son incógnitas de este problema, y dada la
simplicidad de las condiciones de borde, tratemos de encontrar soluciones con φ(x, y) de
variables separables. Es decir, buscamos soluciones de la forma
φ(x, y) = f (x) g(y).
2
Como ∇ φ(x, y) = φxx (x, y)+φyy (x, y) = f ′′ (x)g(y)+f (x)g ′′ (y), resulta que queremos
hallar f (x), g(y) no triviales y λ ∈ R que satisfagan
f ′′ (x)g(y) + f (x)g ′′ (y) = −λf (x)g(y),
es decir,
f ′′ (x) g ′′ (y) f ′′ (x) g ′′ (y)
+ = −λ, =⇒ = −λ − .
f (x) g(y) f (x) g(y)
Los miembros derecho e izquierdo de la última igualdad son funciones de diferentes va-
riables, y por lo tanto, ambos miembros deben ser iguales a una misma constante de
separación que llamaremos −µ. Por otro lado, las condiciones de borde implican que
f (0) = f (L) = 0 y g(0) = g(H) = 0. El problema de hallar autovalores y autofunciones
www.elsolucionario.net
de variables separables para el rectángulo (0, L) × (0, H) resulta entonces:
f ′′ (x)
= −µ, 0 < x < L, f (0) = f (L) = 0,
f (x)
(12.4)
g ′′ (y)
= −(λ − µ), 0 < y < H, g(0) = g(H) = 0.
g(y)
A diferencia del Teorema 12.1 los autovalores están indexados por dos ı́ndices n y m,
en lugar de uno. Podrı́amos ordenarlos y numerarlos con un solo ı́ndice de ser necesario,
pero esto no hace falta en la práctica, además el orden dependerá de cuánto valgan L y
H. Veamos lo que ocurre en el problema de difusión y en el de ondas.
www.elsolucionario.net
(0, L) × (0, H), lo visto en la sección anterior nos permite decir que la solución es
∞ X
∞
nπx mπy −k ( nπ 2 2
X h i
L )
+( mπ
H )
t
u(x, y; t) = an m sen sen e
n=1 m=1
L H
con
Z L Z H
4 nπx mπy
an m = f (x, y) sen sen dy dx, n, m = 1, 2, . . . .
LH 0 0 L H
√ q
nπ 2 mπ 2
con ωn m = λn m = L
+ H
,
Z L Z H
4 nπx mπy
an m = f (x, y) sen sen dy dx,
LH 0 0 L H
y
Z L Z H
1 4 nπx mπy
bn m = g(x, y) sen sen dy dx,
c ωn m LH 0 0 L H
para n, m = 1, 2, . . . .
198
www.elsolucionario.net
∂ 2 φ 1 ∂φ 1 ∂ 2φ
+ + = −λφ, 0 < r < R, −π < θ < π,
∂r2 r ∂r r2 ∂θ2
φ(r, −π) = φ(r, π) 0 < r < R,
∂φ ∂φ (12.5)
(r, −π) = (r, π) 0 < r < R,
∂θ ∂θ
φ(R, θ) = 0 −π < θ < π,
|φ(0, θ)| < ∞ −π < θ < π.
www.elsolucionario.net
Busquemos nuevamente autofunciones de variables separables φ(r, θ) = f (r)g(θ).
Deseamos entonces hallar funciones no triviales f (r), g(θ) y números λ que satisfagan:
1 1
f ′′ (r)g(θ) + f ′ (r)g(θ) + 2 f (r)g ′′ (θ) = −λf (r)g(θ).
r r
r2
Multiplicando por , obtenemos
f (r)g(θ)
199
www.elsolucionario.net
r2 f ′′ (r) + rf ′ (r) + λr2 f (r) = µf (r), 0 < r < R, |f (0)| < ∞, f (R) = 0,
r2 f ′′ (r) + rf ′ (r) + (λr2 − n2 )f (r) = 0, 0 < r < R, |f (0)| < ∞, f (R) = 0. (12.6)
√
Para resolver este problema a valores de borde hacemos
√ el cambio de√variables z = λr
y la identificación F (z) = f (r), o sea F (z) = F ( λr) y f (r) = f (a/ λ), entonces
√
f ′ (r) = fr (r) = Fz (z)zr = F ′ (z) λ, y análogamente f ′′ (r) = λF ′′ (z).
www.elsolucionario.net
que se conoce como Ecuación Diferencial de Bessel de orden n. Las soluciones de esta
ecuación se llaman funciones de Bessel, y se describen a continuación.
Su definición precisa no es tan importante para este curso, dado que vienen implementadas
en diferentes bibliotecas de software1 . Lo importante es saber algunas de sus propiedades:
(
1, si n = 0,
lı́m+ Jn (z) =
z→0 0, si n = 1, 2, . . . .
200
www.elsolucionario.net
1 1
J0 Y0
J1 Y1
J2 Y2
0.5 0.5
J3 Y3
0 0
−0.5 −0.5
−1 −1
0 2 4 6 8 10 0 2 4 6 8 10
Figura 12.1: Gráfico de las funciones de Bessel de primera especie Jn (z) (izquierda) y de
segunda especie Yn (z) (derecha) para n = 0, 1, 2, 3.
www.elsolucionario.net
F (z) = a Jn (z) + b Yn (z), 0 < z < ∞,
1
En MATLAB y en Mathematica se llaman besselj y bessely; en wxMaxima se llaman bessel j y
bessel y.
201
www.elsolucionario.net
1 1
J0 J1
0.5 0.5
z0 4
0 z0 1 z0 2 z0 3 0
z1 1 z1 2 z1 3
−0.5 −0.5
−1 −1
0 2 4 6 8 10 12 0 2 4 6 8 10 12
Figura 12.2: Gráfico de las funciones de Bessel de primera especie J0 (z) (izquierda) y J1 (z)
www.elsolucionario.net
(derecha) y de los ceros z0 m , z1 m , m = 1, 2, . . . .
1 1
J0 (z01 r) J1 (z11 r)
J0 (z02 r) J1 (z12 r)
J0 (z03 r) J1 (z13 r)
0.5 0.5
0 0
−0.5 −0.5
−1 −1
0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1
202
www.elsolucionario.net
Z R
fn m1 (r)fn m2 (r) r dr = 0, si m1 6= m2 . (12.10)
0
n = 0:
z 2 z0 m
0m
λ0 m = , φ 0 m = J0 ( r), m = 1, 2, 3 . . . .
R R
n = 1, 2, 3, . . . :
www.elsolucionario.net
z 2
nm
λn m = ,
R
zn m
φ1n m (r, θ) = Jn ( r) cos nθ,
R
zn m
φ2n m (r, θ) = Jn ( r) sen nθ, m = 1, 2, 3, . . . ,
R
ZZ Z π Z R
φin m φjk ℓ dA = φin m (r, θ)φjk ℓ (r, θ) r dr dθ = 0, si i 6= j ó n 6= k ó m 6= ℓ.
DR −π 0
(12.11)
∞
X z0 m −( z0Rm )2 kt
u(r, θ; t) = a0 m J 0 ( r)e
m=1
R
∞ X
∞
X zn m zn m 2
+ Jn ( r) (an m cos nθ + bn m sen nθ) e−( R ) kt ,
n=0 m=1
R
203
www.elsolucionario.net
con RR
f (r, θ)Jn ( z0Rm r) dA
DR
a0 m = RR 2 , m = 1, 2, . . . ,
z0 m
J n ( R
r) dA
RR DR zn m
DR
f (r, θ)J n ( R r) cos nθ dA
an m = RR 2 , n = 1, 2, . . . , m = 1, 2, . . . ,
zn m
DR
J n ( R
r) cos nθ dA
RR
f (r, θ)Jn ( znRm r) sen nθ dA
DR
bn m = RR zn m 2 n = 1, 2, . . . , m = 1, 2, . . . .
DR
Jn ( R r) sen nθ dA
Para simplificar, suele denotarse con λ = (n, m, i) en el conjunto de ı́ndices
y para λ = (n, m, i)
zλ
J0 ( R r),
si n = 0,
i
φλ (r, θ) = φn m (r, θ) = Jn ( zRλ r) cos nθ, si n > 0 e i = 1,
www.elsolucionario.net
zλ
Jn ( R r) sen nθ, si n > 0 e i = 2,
con zλ = zn m .
Luego, la solución del problema de difusión puede escribirse simplificada de la siguiente
manera: X zλ 2
u(r, θ; t) = aλ φλ (r, θ)e−( R ) kt ,
λ∈Λ
con RR
DR
f φλ dA
aλ = RR .
DR
[φλ ]2 dA
con RR
D
f φλ dA
aλ = RR R ,
DR
φ2λ dA
y RR
R g φλ dA
bλ = RRDR 2 .
c zλ φ dA
DR λ
204
www.elsolucionario.net
1, 2, . . . dan cero, y sólo sobreviven los coeficientes con ı́ndice n = 0. Por lo tanto, la
solución de la ecuación de difusión resulta
∞
X z0 m −( z0Rm )2 kt
u(r, t) = a0m J0 ( r)e ,
m=1
R
con RR
f (r)J0 ( z0Rm r) r dr
a0m = R0 R 2 , m = 1, 2, 3, . . . .
z0 m
0
J 0 ( R
r) r dr
Por otro lado, la solución de la ecuación de ondas, bajo la hipótesis de simetrı́a radial,
resulta ∞
X z0 m h z0 m z0 m i
u(r, t) = J0 ( r) a0m cos( ct) + b0m sen( ct) ,
m=1
R R R
con
RR RR
0
f (r)J0 ( z0Rm r) r dr R 0 g(r)J0 ( z0Rm r) r dr
a0m = RR 2 , y b0m = R , m = 1, 2, 3, . . . .
J 0 ( z0 m
r) r dr cz0 m R J0 ( z0 m r) 2 r dr
0 R 0 R
www.elsolucionario.net
12.5. Ecuación de Laplace en un cilindro circular.
Funciones de Bessel modificadas
2
Como ya hemos visto, la ecuación de Laplace ∇ u = 0 modela el estado estacionario
de la ecuación de difusión en ausencia de términos fuente. Hemos resuelto la ecuación de
Laplace en un rectángulo y en un cı́rculo. La ecuación de Laplace en un prisma rectangular
también puede resolverse por el método de separación de variables (ver Ejercicio 12.7).
Como puede verse en los ejercicios al final del capı́tulo, las tres variables independientes
dan dos problemas de autovalores con soluciones oscilatorias, y el restante problema tiene
soluciones no-oscilatorias.
Un problema más desafiante presenta la ecuación de Laplace en un cilindro de radio
R y altura H. Utilizando las coordenadas cilı́ndricas
x = r cos θ, y = r sen θ, z = z,
205
www.elsolucionario.net
u(r, −π, z) = u(r, π, z), uθ (r, −π, z) = uθ (r, π, z), 0 < r < R, 0 < z < H. (12.12)
Nos encontramos frente a tres condiciones de borde no homogéneas. Como antes, podemos
dividir este problema en tres subproblemas, cada uno de los cuales involucra resolver
la ecuación de Laplace, pero tiene condición de borde no-homogénea en una cara, y
homogénea en las otras dos. Más precisamente,
u = u 1 + u 2 + u3 ,
con
2 2 2
∇ u1 = 0
∇ u2 = 0
∇ u3 = 0
u1 (r, θ, 0) = α(r, θ), u2 (r, θ, 0) = 0, u3 (r, θ, 0) = 0,
u1 (r, θ, H) = 0,
u2 (r, θ, H) = β(r, θ),
u3 (r, θ, H) = 0,
u1 (R, θ, z) = 0, u2 (R, θ, z) = 0, u3 (R, θ, z) = γ(θ, z),
www.elsolucionario.net
además de las condiciones de periodicidad (12.12).
A continuación, comenzamos a resolver por el método de separación de variables una
sola vez, y luego tratamos los diferentes casos individualmente.
2
Separación de variables. Buscamos soluciones de la ecuación de Laplace ∇ u = 0 de
la forma u(r, θ, z) = f (r)g(θ)h(z). Entonces f (r), g(θ) y h(z) deben satisfacer
1 1
f ′′ (r)g(θ)h(z) + f ′ (r)g(θ)h(z) + 2 f (r)g ′′ (θ)h(z) + f (r)g(θ)h′′ (z) = 0.
r r
Dividiendo por f (r)g(θ)h(z) obtenemos
La suma de los tres primeros términos depende de las variables r, θ y el cuarto término
depende sólo de z, por lo que concluimos que para alguna constante de separación λ,
206
www.elsolucionario.net
Las condiciones de borde periódicas (12.12) deben cumplirse en los tres casos, y se tra-
ducen en g(−π) = g(π), g ′ (−π) = g ′ (π). Esto implica que los autovalores µ son n2 ,
n = 0, 1, 2, . . . , con autofunciones:
n=0 g(θ) = 1,
N = 1, 2, 3 . . . , g(θ) = cos nθ, sen nθ.
Para cada n = 0, 1, 2, 3 . . . deberemos hallar los valores de λ para los cuales las siguientes
funciones tienen solución no trivial, y satisfacen las CB correspondientes (en cada caso):
h′′ (z)
= λ, r2 f ′′ (r) + r f ′ (r) + (λr2 − n2 )f (r) = 0.
h(z)
Para saber cuál de estas dos ecuaciones determina los autovalores λ, hace falta considerar
las condiciones de borde homogéneas. Éstas dependen de si estamos buscando u1 , u2 ó u3 .
www.elsolucionario.net
2
∇ u1 = 0
u1 (r, θ, 0) = α(r, θ),
u1 (r, θ, H) = 0,
u1 (R, θ, z) = 0.
h(H) = 0,
f (R) = 0.
Una condición de borde que debemos considerar se debe a la singularidad creada por
el cambio de variables en r = 0, que como en el caso de coordenadas polares debemos
imponer
|f (0)| < ∞.
Los autovalores λ estarán determinados por el problema que tiene dos condiciones de
borde homogéneas, que en este caso es el siguiente:
Las autofunciones de este problema están dadas por las funciones de Bessel de primera
especie, y resultan, para cada n = 0, 1, 2, . . . , las funciones (ver Figuras 12.2–12.3)
zn m
fn m (r) = Jn ( r), donde zn m es el m-ésimo cero positivo de Jn , m = 1, 2, 3 . . . ,
R
2
y los autovalores son λn m = znRm . Ahora, para estos valores de λn m debemos ha-
2
llar las funciones h(z) correspondientes. Como λn m = znRm > 0 resulta que h(z) =
207
www.elsolucionario.net
zn m zn m
C1 cosh z + C2 senh z. Para imponer la condición inicial conviene trasladar las
R R
funciones en H, escribiendo
zn m zn m
h(z) = C1 cosh (z − H) + C2 senh (z − H).
R R
zn m
Ası́ h(H) = 0 implica C1 = 0 y luego h(z) = C2 senh (z − H).
R
Resumiendo, las soluciones de variables separables que cumplen las mismas condicio-
nes de borde homogéneas que u1 son:
z0 m
n = 0, a0 m f0 m (r) senh( (z − H)), m = 1, 2, . . . ,
R
zn m
n = 1, 2, . . . , fn m (r) senh( (z − H)) (an m cos nθ + bn m sen nθ) , m = 1, 2, . . . ,
R
y una solución más general se obtiene sumándolas:
∞
X z0 m
www.elsolucionario.net
u1 (r, θ, z) = a0 m f0 m (r) senh (z − H)
m=1
R
∞ X
∞
X zn m
+ fn m (r) senh (z − H) (an m cos nθ + bn m sen nθ) .
n=1 m=1
R
Z π Z π
2
cos nθ dθ = sen2 nθ dθ = π, n = 1, 2, . . .
−π −π
208
www.elsolucionario.net
www.elsolucionario.net
es similar al de u1 y se deja como ejercicio (ver Ejercicio 12.9).
209
www.elsolucionario.net
lı́m Kn (w) = ∞, n = 0, 1, 2, . . . .
w→0+
lı́m Kn (z) = 0, n = 0, 1, 2, . . . .
w→+∞
www.elsolucionario.net
In (w) y Kn (w) son funciones que no oscilan y más aún, son positivas para todo
w > 0; ver Figura 12.4.
10 10
I0 K0
I1 K1
8 I2 8 K2
I3 K3
6 6
4 4
2 2
0 0
0 2 4 6 8 10 0 2 4 6 8 10
Figura 12.4: Gráfico de las funciones de Bessel modificadas de primera especie In (w)
(izquierda) y de segunda especie Kn (w) (derecha) para n = 0, 1, 2, 3.
210
www.elsolucionario.net
mπ
La condición de borde |f (0)| < ∞ implica b = 0 pues lı́m+ Kn ( r) = +∞. Luego la
r→0 H
solución de (12.13) que cumple la condición de borde |f (0)| < ∞ es
mπ
f (r) = a In ( r).
H
Resumiendo, las soluciones de variables separables que cumplen las mismas condicio-
nes de borde homogéneas que u3 son:
mπ mπz
n = 0, a0 m I 0 ( r) sen , m = 1, 2, . . . ,
H H
mπ mπz
n = 1, 2, . . . , In ( r) sen (an m cos nθ + bn m sen nθ) , m = 1, 2, . . . ,
H H
y una solución más general se obtiene sumándolas:
∞
X mπ mπz
u3 (r, θ, z) = a0 m I 0 ( r) sen
m=1
H H
∞ X
∞
www.elsolucionario.net
X mπ mπz
+ In ( r) sen (an m cos nθ + bn m sen nθ) .
n=1 m=1
H H
211
www.elsolucionario.net
ZZ Z π Z H
Usando que v(θ, z) dA = v(θ, z) R dz dθ, y que entonces
CR −π 0
ZZ h Z π Z H
mπz i2 mπz H
sen dA = R dθ sen2 dz = R2π = RπH,
CR H −π 0 H 2
ZZ Z π Z H
mπz 2 mπz RπH
sen cos nθ dA = R cos2 nθ dθ sen2 dz = ,
CR H −π 0 H 2
ZZ Z π Z H
mπz 2 2 mπz RπH
sen sen nθ dA = sen nθ dθ sen2 dz = ,
CR H −π 0 H 2
www.elsolucionario.net
Rπ RH
2 −π 0 γ(θ, z) sen mπz sen nθ dz dθ
bn m = mπ
H , n = 1, 2, . . . , m = 1, 2, . . . .
In H R πH
12.6. Ejercicios
12.1. Consideremos una membrana vibrante de forma genérica Ω levemente amortiguada
que satisface
212
www.elsolucionario.net
12.4. Decir cuáles son (aproximadamente) los primeros 6 autovalores del Laplaciano y sus
correspondientes autofunciones en Ω = D1 es decir, el disco de radio 1. Notar que algunos
autovalores se repiten con autofunciones ortogonales. Decir cuáles son las 6 frecuencias
principales de vibración de una membrana elástica (con T0 = 2 y ρ = 1) que ocupa el disco
unitario Ω = D1 . Para este ejercicio deberá hallar de manera aproximada el valor de los
ceros de las primeras funciones de Bessel de primera especie. Puede hacerlo, por ejemplo,
mediante la aplicación de zoom en los gráficos de las funciones de Bessel, o buscando en
alguna tabla. Puede resultarle útil saber que zn1 1 < zn2 1 siempre que n1 < n2 .
zn m
* 12.5. Demostrar que si n = 0, 1, 2, . . . , m = 1, 2, . . . y fn m (r) = Jn ( r) entonces
R
vale (12.10), es decir
Z R
fn m1 (r)fn m2 (r) r dr = 0, si m1 6= m2 .
0
2
Ayuda: Usar que r2 fn′′ m (r) + rfn′ m (r) + (λn m r2 − n2 )fn m = 0, para λn m = znRm , y demostrar
que
n2
www.elsolucionario.net
′
λn m fn m (r)r = fn m (r) − rfn′ m (r) .
r
Multiplicar esta igualdad para m = m1 por fn m2 , integrar por partes dos veces, usar la condición
de borde fn m (R) = 0 para arribar a
Z R Z R
λn m 1 fn m1 (r)fn m2 (r) r dr = λn m2 fn m1 (r)fn m2 (r) r dr.
0 0
(a)
u(x, y, 0) = 0, u(x, y, W ) = 0, 0 < x < L, 0 < y < H,
u(x, 0, z) = f (x, z), u(x, H, z) = 0, 0 < x < L, 0 < z < W,
u(0, y, z) = 0, u(L, y, z) = 0 0 < y < H, 0 < z < W.
(b)
u(x, y, 0) = 0, u(x, y, W ) = 0, 0 < x < L, 0 < y < H,
uy (x, 0, z) = 0, uy (x, H, z) = 0, 0 < x < L, 0 < z < W,
u(0, y, z) = 0, u(L, y, z) = f (y, z) 0 < y < H, 0 < z < W.
(c) ¿Cómo resulta la solución del item anterior si f (y, z) sólo depende de z, es decir
f (y, z) = f (z)?
213
www.elsolucionario.net
(a)
u(r, θ, 0) = 0, 0 < r < R, −π < θ < π,
u(r, θ, H) = β(r, θ), 0 < r < R, −π < θ < π,
u(R, θ, z) = 0, −π < θ < π, 0 < z < H.
(b) (Notar que las condiciones de borde son de tipo Neumann sobre las caras planas)
(c) ¿Cómo resulta la solución del ı́tem anterior si f (θ, z) sólo depende de θ, es decir
f (θ, z) = f (θ)? ¿Qué relación guarda con la solución del problema de Laplace en el
disco?
Bibliografı́a complementaria
www.elsolucionario.net
[Bleecker-Csordas 1996] Bleecker, D., Csordas, G., Basic Partial Differential Equations,
International Press, Cambridge, Massachusetts, 1996.
[Haberman 1998] Haberman, R., Elementary Applied Partial Differential Equations, Pren-
tice Hall, Upper Saddle River, NJ, 1998.
214