12 La ecuación del calor
La ecuación,
@u @ 2u
=k 2 (12.1)
@t @x
donde k > 0, se llama la ecuación del calor. A menudo se la conoce también como
ecuación de difusión. Por diferentes motivos esta, y sus variaciones, son extremada-
mente relevantes. En general se la piensa como una ecuación de evolución con
parámetro ‘temporal’ t.
12.1 El problema de Cauchy en la recta.
12.1.1 Existencia de soluciones con dato inicial de soporte compacto.
Denotamos por ⌦ = ( 1, 1) ⇥ [0, 1) y de coordenadas (x, t). En ⌦ consideramos
el problema de Cauchy, (
ut = kuxx ,
(12.2)
u(x, 0) = h(x),
donde el dato inicial h es C 1 y de soporte compacto. Buscamos una solución u, C 1
en ⌦. La manera usual de resolver el problema es a través del núcleo del calor en la
recta,
1 2
H(x, y, t) = p e (x y) /4kt (12.3)
4⇡kt
Fijo y, u(x, t) = H(x, y, t) es solución de (12.1) en ⌦ . Además, H(x, y, t) posee la
siguiente invariancia de escala,
1 x µ y µ
p H( p , p , t) = H(x, y, t), (12.4)
para todo > 0 y µ. De esta invariancia y de la forma Gaussiana normal de
H(x, 0, 1/2k),
1 2
H(x, 0, 1/2k) = p e x /2 (12.5)
2⇡
se deduce que, fijado y, Dt (x) = H(x, y, t) es un núcleo de Dirac en y cuando t tiende
a cero, esto es,
1. Dt (x) > 0,
R1
2. 1 Dt (x)dx = 1,
3. para todo ✏ > 0, lim max{Dt (x) : x 2 ( 1, y ✏) [ (y + ✏, 1)} = 0 cuando
t ! 0.
90
La integral, Z 1
u(x, t) = H(x, y, t)h(y)dy (12.6)
1
define una solución de (12.1) en ⌦ . Debido a que Dt (x) es un núcleo de Dirac,
dicha solución se extiende a una función continua en ⌦ con dato inicial h(x). Para
probar que es la solución que buscamos debemos mostrar que u es C 1 en ⌦. Para
eso escribimos u como,
Z 1
u(x, t) = H(z, 0, t)h(x z)dz (12.7)
1
Esta expresión muestra que para todo l 1, @x(l) u, como función en ⌦ , se extiende
a una función continua en ⌦ con dato inicial @x(l) h. Esto es suficiente para probar
que u es C 1 en ⌦. De hecho, debido a que u satisface la ecuación del calor en ⌦ ,
tenemos que, en ⌦ y para todo l 0 y j 1,
(j)
@t @x(l) u = k l @x(2j+l) u (12.8)
Se sigue que u es C 1 en ⌦ y que satisface la ecuación del calor.
12.1.2 Unicidad.
La solución u hallada en la sección enterior es única en el espacio de soluciones v tal
que lim v(x, t) = lim vx (x, t) = 0 cuando x ! ±1. La prueba de esto es la siguiente.
Sean u1 y u2 soluciones C 1 con el mismo dato inicial. Entonces u = u1 u2 tiene
dato inicial cero y es solución. Multiplicando (12.1) por u e integrando por partes
en x, obtenemos,
✓Z 1 ◆ Z 1
1d 2
u (x, t)dx = k u2x (x, t)dx (12.9)
2 dt 1 1
Por lo tanto la función de t, Z 1
u2 (x, t)dx (12.10)
1
es no-negativa, es monotónicamente decreciente y vale cero en t = 0. Se sigue que u
es idénticamente cero y por tanto u1 = u2 .
12.1.3 Dato inicial contı́nuo y ’smoothing’.
Si h es solamente continua, la fórmula (12.3) define una solución C 1 en ⌦ que se
extiende a una función continua en ⌦ con dato inicial h(x). La solución es única en
el espacio de las soluciones v tal que para todo t > 0, lim v(x, t) = lim vx (x, t) = 0
cuando x ! ±1. La prueba de esto se realiza de la misma manera que en las
secciones anteriores.
Independientemente de la irregularidad del dato inicial contı́nuo h, para todo t
91
u(x, t) es una función C 1 de x. Este hecho se conoce como ’suavizado’ o ’smoothing’
y es una propiedad fundamental del flujo del calor.
12.2 El problema de Cauchy en un intervalo con condiciones de frontera
de Dirichlet.
Denotamos ⌦ = [0, L] ⇥ [0, 1). En ⌦ consideramos el problema de Cauchy con
condición de frontera de Dirichlet,
8
>
< ut = kuxx ,
u(x, 0) = h(x), (12.11)
>
: u(0, t) = u(⇡, t) = 0
donde el dato inicial h es C 1 y, por compatibilidad con la condición de borde y la
ecuación del calor, debe satisfacer
0000 000000
h(0) = 0, h00 (0) = 0, h (0) = 0, h (0) = 0, . . . (12.12)
0000 000000
h(⇡) = 0, h00 (⇡) = 0, h (⇡) = 0, h (⇡) = 0, . . . (12.13)
Buscamos una solución u, C 1 en ⌦.
Observe que tomamos por simplicidad L = ⇡. La solución con L arbitrario se
obtiene de esta con la transformación...
12.2.1 Soluciones por separación de variables.
Al igual que con la ecuación de Laplace, buscamos soluciones de la forma,
u(x, t) = X(x)T (t) (12.14)
con u función C 1 . Reemplazando esta expresión de u en (12.1) obtenemos,
XTt = kXxx T (12.15)
Al igual que cómo analizamos separación de variables en en la ecuación de Laplace,
si u no es idénticamente cero entonces existe una constante tal que,
Xxx = X, Tt = kT (12.16)
Para satisfacer las condiciones de borde es necesario que = n2 = 1, 4, 9, . . ., de
aquı́ que las soluciones para cada n se escriben de la forma,
kn2 t
un = bn e sin nx (12.17)
Al menos formalmente, la solución más general que puede construirse con las fun-
ciones un es,
1
X
n2 kt
u= bn e sin nx (12.18)
n=1
92
Para satisfacer la condición inicial u(x, 0) = h(x), los coeficientes an deben ser los
coeficientes de Fourier de la extensión impar de h a [ ⇡, ⇡]. Por lo tanto,
Z ⇡
2
bn = h(y) sin nydy (12.19)
⇡ 0
Formalmente entonces, la solución es (12.18) y los coeficientes son (12.19).
De las condiciones (12.12) se deduce, luego de integrar por partes al igual que en
la Sección 10.0.1, que
1
bn (h) = k bn (h(k) ) (12.20)
n
Por lo tanto bn decae más rápido que toda potencia de n. Se sigue que u es C 1 en
[0, ⇡] ⇥ [0, 1) y que satisface la ecuación del calor, las condiciones de frontera y el
dato inicial. Por lo tanto es una solución C 1 al problema de Cauchy.
12.2.2 Unicidad.
Para ver que la solución es única procedemos igual que en la sección 12.1.2. Si u1 y
u2 son dos soluciones C 1 del problema de Dirichlet (12.11) y hacemos u = u1 u2 ,
entonces obtenemos,
Z ⇡ Z ⇡
d 2
u(x, t) dx = k ux (x, t)2 dx (12.21)
dt 0 0
R⇡
Por tanto 0 u2 (x, t)dx es no-negativa, es monotónicamente decreciente y a tiempo
cero vale cero. Esto implica que u es idénticamente cero y por lo tanto u1 = u2 .
12.2.3 El principio del máximo.
Existe un principio del máximo para las soluciones de la ecuación del calor del que
puede deducirse también la unicidad de soluciones. Dada su relevancia, en esta
sección damos una prueba.
Teorema 12.1. Principio del máximo. El máximo y el mı́nimo de toda solución
C 1 de la ecuación del calor en
⌦ = (a, b) ⇥ (0, T ] (12.22)
que es continua en [a, b] ⇥ [0, T ] se alcanzan en el siguiente borde de ⌦,
S = {0} ⇥ [0, T ] [ [a, b] ⇥ {0} [ {b} ⇥ [0, T ] (12.23)
Prueba. Sea M el máximo de u en ⌦ y sea m el máximo de u en el borde S. Supong-
amos que M > m y sea (x0 , t0 ) un punto en ⌦ donde se alcanza M . Sea l = (b a)/2.
Definimos la función,
M m
v(x, t) = u(x, t) + (x x0 ) 2 (12.24)
4l2
93
Observamos que en S tenemos,
M m
v(x, t) < m + l2 < M (12.25)
4l2
y que en (x0 , t0 ) tenemos,
v(x0 , t0 ) = u(x0 , t0 ) = M (12.26)
Por lo tanto el máximo de v en ⌦ se alcanza en ⌦. Sea (x1 , t1 ) un punto en ⌦ donde
se alcanza el máximo de v. Naturalmente debemos tener
vxx (x1 , t1 ) 0. (12.27)
Por otro lado si t1 < T entonces vt (x1 , t1 ) = 0 mientras que si t1 = T tenemos,
vt (x1 , t1 ) 0, (12.28)
de lo contrario M no se alcanza en (x1 , t1 ).
Ahora de la definición de v se desprende que,
vt (x1 , t1 ) = ut (x1 , t1 ), (12.29)
M m
vxx (x1 , t1 ) = uxx (x1 , t1 ) + (12.30)
l2
Como u satisface la ecuación del calor tenemos,
M m
vt (x1 , t1 ) = kvxx (x1 , t1 ) k < kvxx (x1 , t1 ) (12.31)
l2
lo que contradice (12.27) y (12.28).
Ejercicio 12.2. Probar la unicidad de la sección 12.2.2 usando el principio del
máximo.
En la siguiente sección brindamos otra applicación del principio del máximo.
12.2.4 Dato inicial contı́nuo y condición de frontera de Dirichlet. Unici-
dad.
Si el dato inicial h es solamente contı́nuo, la solución u en (12.18) es C 1 en [0, ⇡] ⇥
2
(0, 1), ya que, debido a la presencia de los coeficienes e kn t , la serie (12.18) converge
uniformemente en C 1 en todo rectángulo cerrado que no toque al intervalo inicial
[0, ⇡] ⇥ {0}. La serie (12.18) puede diverger en t = 0 ya que h es solamente continua.
Un problema similar ocurrı́a con la fórmula de Poisson. Sin embargo también para
la ecuación del calor, si h es continua entonces el lı́mite de u en t = 0 es h.
Una manera de probar esto es ver que u(x, t) es igual a la restricción a ⌦ = [0, ⇡] ⇥
[0, 1) de la solución û(x, t) en la recta, con dato inicial ĥ igual a la extensión impar
de h a R. Discutimos esto en lo que sigue. Para comenzar, observamos que, como ĥ es
94
impar y periódica, la fórmula (12.6) muestra que û(0, t) = û(⇡, t) = 0 para todo t 0.
Por otro lado û es continua y û(x, 0) = h(x). Consideramos ahora una secuencia ĥi
de funciones impares C 1 , con ĥi (0) = ĥi (⇡) = 0 que convergen uniformemente a
ĥ. Sea hi la restricción de ĥi a [0, ⇡]. Sean ui las soluciones en [0, ⇡] ⇥ [0, 1) con
dato inicial hi y condición de frontera de Dirichlet ui (0, t) = ui (⇡, t) = 0. De la serie
(11.13) se deduce que para todo t > 0,
max{|ui (x, t) u(x, t)| : x 2 [0, ⇡]} ! 0 (12.32)
cuando i ! 1. Por otro lado, del principio del máximo aplicado a û ui en
[0, ⇡] ⇥ [0, t] deducimos,
max{|û(x, t) ui (x, t)| : x 2 [0, ⇡]} ! 0 (12.33)
cuando i ! 1. Se sigue que u(x, t) = û(x, t) para todo t > 0. Por lo tanto u se
extiende continuamente a [0, ⇡] ⇥ [0, 1) con dato inicial h.
Otra forma de probar esto es escribiendo la solución a la ecuación del calor en la
forma, Z ⇡
u(x, t) = Q(x, y, t)h(y)dy (12.34)
0
donde Q(x, y, t) es el núcleo del calor (para este problema en particular) y es igual
a,
1
2X 2
Q(x, y, t) = e kn t sin nx sin ny (12.35)
⇡ n=1
y luego estudiar el lı́mite de Q(x, y, t) cuando t tiende a cero, de la misma forma que
se hizo con el núcleo de Poisson. La diferencia con el núcleo de Poisson es que Q no
posee una expresión cerrada y su análisis directo no es sencillo. La figura 14 muestra
los gráficos de Q(x, y, t) para y = ⇡/2 y tres valores de t pequeños.
Se observa que, fijo y, el comportamiento de Q(x, y, t) para t ! 0, es reminiscente
a un núcleo de Dirac. Esto es lo que realmente sucede. De hecho, fijado y, Q(x, y, t)
es aproximado por el núcleo de la ecuación del calor en la recta, H(x, y, t), cuando t
tiende a cero. El resultado se resume en el siguiente ejercicio.
Teorema 12.3. Sea y 2 (0, ⇡) y sea Di (x) un núncleo de Dirac en y con soporte
en (0, ⇡). Cuando i ! 1, la soluciónes ûi a la ecuación del calor en la recta con
dato inicial Di convergen a H(x, y, t), y las soluciónes ui a la ecuación del calor en
el intervalo [0, ⇡] con condición de frontera de Dirichlet nula, convergen a Q(x, y, t).
Usar este hecho y el principio del máximo para probar que,
max{|H(x, y, t) Q(x, y, t)| : x 2 [0, ⇡]} ! 0 (12.36)
cuando t tiende a cero.
La unicidad de soluciones puede probarse al igual que en la sección 12.1 (evitando
95
100
80
60
40
20
-20
-2 -1 0 1 2 3 4 5
Figure 14: El núcleo de calor t = 0.01, 0.001, 0.0001.
t = 0) o usando el principio del máximo.
12.2.5 El problema de Cauchy hacia el pasado.
Puede parece que el problema de Cauchy hacia el pasado (por ejemplo, con condi-
ciones de frontera de Dirichlet nulas) se resuelve de la misma forma que hacia el
futuro, pero no es el caso. El problema de Cauchy hacia el pasado no está bien
puesto. En general, dado un dato inicial, no existe una solución hacia el pasado.
Veamos esto en un ejemplo. Consideremos el dato inicial h dado por la serie,
1
X
n
h= e sin nx (12.37)
n=1
El dato h es C 1 por el rápido decaimiento de los coeficientes. Sin embargo h no puede
ser el dato inicial de una solución u(x, t) para t < 0. Supongamos que u(x, t) existe
con t 2 [t0 , 0], para algún t0 < 0. Entonces u(x, 0) es el resultado de la evolución a
tiempo t0 de u(x, t0 ). Por lo tanto los coeficientes de Fourier deben decaer tan o
2
más rápido que en /kt0 que no es el caso.
’Genéricamente’, los datos iniciales no general soluciones hacia el pasado. Este es
un hecho caraterı́stico de la ecuación del calor. Notar sin embargo que algunos datos
iniciales, como por ejemplo h = sin nx sı́ poseen una evolució hacia el pasado.
96
13 La ecuación de ondas.
La ecuación,
@ 2u 2
2@ u
c =0 (13.1)
@t2 @x2
donde c > 0 es la velocidad de propagación, se llama la ecuación de ondas. En general
se la piensa como una ecuación de evolución con parámetro ’temporal’ t. La ecuación
es simétrica hacia el futuro (t > 0) o el pasado (t < 0) y por lo tanto puede resolverse
en la dos direcciones temporales.
13.1 El problema de Cauchy en la recta.
Denotamos ⌦ = ( 1, 1) ⇥ [0, 1). En ⌦ consideramos el problema de Cauchy,
8
>
< utt c2 uxx = 0,
u(x, 0) = h(x), (13.2)
>
: u (x, 0) = w(x)
t
donde h y w son funciones C 1 sin más restricciones. Buscamos una solución C 1 en
⌦.
13.1.1 La fórmula de d’Alambert.
Una manera explı́cita de resolver la ecuación de ondas en la recta con la fórmula de
d’Alambert,
Z x+ct
1 1
u(x, t) = h(x + ct) + h(x ct) + w(y)dy (13.3)
2 2c x ct
Es fácil interpretar el primer paréntesis como el dato o posición inicial h desplazan-
dose hacia la derecha y la izquierda con velocidad c (y con un coeficiente 1/2). El
segundo término solo involucra a la velocidad inicial. Volveremos sobre esto más
adelante.
La fórmula de d’Alambert se prueba del siguiente modo. Observamos primero que
(13.1) se puede escribir como,
✓ ◆✓ ◆
@ 1@ @ 1@
+ u=0 (13.4)
@x c @t @x c @t
Hacemos el cambio de variables ⇠ = x + ct, ⌘ = x ct. Por lo tanto,
✓ ◆ ✓ ◆
@ 1 @ 1@ @ 1 @ 1@
= + , = (13.5)
@⇠ 2 @x c @t @⌘ 2 @x c @t
97
De aquı́ escribimos la ecuación (13.4) como,
@ @
u=0 (13.6)
@⇠ @⌘
Ejercicio 13.1. Si u C 1 satisface (13.6) entonces u es de la forma,
u = p(⇠) + q(⌘) (13.7)
con p y q funciones C 1 .
Para satisfaces los datos iniciales debemos tener,
h(y) = p(y) + q(y), w(y) = c(p0 (y) q 0 (y)) (13.8)
Derivando la primera identidad concluimos,
h 0 = p0 + q 0 , w = c(p0 q0) (13.9)
De aquı́ que,
1 1 1 1
p0 = (h0 + w), q 0 = (h0 w) (13.10)
2 c 2 c
Integrando entre 0 y y, y reemplazando el resultado en (13.7) obtenemos la fórmula
de d’Alambert (13.3).
13.1.2 Energı́a y unicidad.
La solución de la sección anterior se halló constructivamente y por lo tanto es única.
La solucón al problema de Cauchy debe ser necesariamente la dada por la fórmula
de d’Alambert.
Hay otra forma de demostrar unicidad que es utilizando la energı́a E. Esto re-
quiere sin embargo de datos iniciales con energı́a finita. Para simplificar el análisis
asumamos que tanto h como w son de soporte compacto. La fórmula de d’Alamabert
muestra que u(x, t) = 0 también es de soporte compacto para todo t 0. La energı́a,
Z 1
1 2
E(t) = 2
ut + u2x dx (13.11)
1 c
es conservada ya que derivando con respecto a t obtenemos,
Z 1 Z 1
dE 1
= 2 2 ut utt + utx ux dx = 2 ut uxx + utx ux dx = 0, (13.12)
dt 1 c 1
donde la última identidad se sigue de integrar por partes.
Por lo tanto E(t) = E(0) para todo t. Si E(0) = 0 entones u(x, t) = 0. La
unicidad es ahora directa. Si u1 y u2 son dos soluciones al problema de Cauchy con
dato inicial de soporte compacto entonces u1 u2 tiene energı́a inicial cero y por lo
tanto u1 = u2 .
98
13.2 El problema de Cauchy en un intervalo con condiciónes de frontera
de Dirichlet.
Denotamos ⌦ = [0, ⇡] ⇥ [0, 1). Consideramos el siguiente problema de Cauchy con
condiciones de frontera de Dirichlet nulas,
8
>
> utt c2 uxx = 0,
>
< u(x, 0) = h(x),
(13.13)
>
> ut (x, 0) = w(x),
>
:
u(0, t) = u(⇡, t) = 0
donde los datos iniciales h y w son C 1 , y, por compatibilidad con la condición de
frontera y la ecuación de ondas, debe satisfacer,
0000 000000
f (z) = 0, f 00 (z) = 0, f (z) = 0, f (z) = 0, . . . (13.14)
para f = h, f = w, z = 0 y z = ⇡. Buscamos una solución C 1 en ⌦.
La solución con L arbitrario se obtiene de esta haciendo...
13.2.1 Soluciones por separación de variables.
Al igual que con la ecuación de Laplace y la ecuación del calor, buscamos soluciones
de la forma,
u(x, t) = X(x)T (t) (13.15)
con u función C 1 . Reemplazando esta expresión de u en (12.1) obtenemos,
c2 XTtt = Xxx T (13.16)
Al igual que el Ejercicio ??, si u no es idénticamente cero entonces existe una con-
stante tal que,
Xxx = X, Ttt = c2 T (13.17)
Para satisfacer las condiciones de borde es necesario que = n2 = 1, 4, 9, . . ., de
aquı́ que las soluciones para cada n se escriben de la forma,
un = (an sin nct + bn cos nct) sin nx (13.18)
Al menos formalmente, la solución más general que puede construirse con las fun-
ciones un es,
1
X
u= (an cos nct + bn sin nct) sin nx (13.19)
n=1
Para satisfacer la condición inicial u(x, 0) = h(x), los coeficientes an deben ser los
coeficientes de Fourier de la extensión impar de h a [ ⇡, ⇡]. Por lo tanto,
Z ⇡
2
an = h(y) sin nydy (13.20)
⇡ 0
99
De igual manera, los coeficientes ncbn deben ser los coeficientes de Fourier de la
extensión impar de w a [ ⇡, ⇡]. Es decir,
Z ⇡
2
bn = w(y) sin nydy (13.21)
nc⇡ 0
Formalmente entonces, la solución es (13.19) y los coeficientes son (13.20) y
(13.21).
De las condiciones (12.12) se deduce, luego de integrar por partes al igual que en
la Sección 10.0.1, que
1 1
an (h) = k
an (h(k) ), bn (w) = k
bn (w(k) ) (13.22)
n n
Por lo tanto bn decae más rápido que toda potencia de n. Se sigue que u es C 1 en
[0, ⇡] ⇥ [0, 1) y que satisface la ecuación de ondas, las condiciones de frontera y el
dato inicial. Por lo tanto es una solución C 1 al problema de Cauchy.
100
