ECUACIONES DIFERENCIALES PARCIALES

(22 de septiembre de 2021)

Lectura 1

Saul Ancari Villca

Índice
Introducción 1

1. La ecuación de Laplace 4
1.1. Solución fundamental . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
1.2. Ecuación de Poisson . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
1.3. Fórmula del valor medio . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.4. Propiedades de las funciones armónicas . . . . . . . . . . . . . . . . . . . . 9
1.4.1. Principio del máximo fuerte, unicidad . . . . . . . . . . . . . . . . 9
1.4.2. Regularidad . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
1.4.3. Teorema de Liouville . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
1.5. Función de Green . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13

2. La ecuación del calor 15

2.1. Solución fundamental . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
2.2. Problema de valor inicial . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
2.3. Problema no-homogéneo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
2.3.1. Solución del problema no homogéneo con condiciones iniciales
generales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
2.4. Propiedades de la solución . . . . . . . . . . . . . . . . . . . . . . . . . . . 21

Ejercicios 24

Referencias 26

Introducción
Este curso está destinado a cubrir algunos hechos básicos del análisis de ecuaciones
diferenciales parciales (EDP). Una ecuación diferencial parcial es una ecuación que involu-
cra una función desconocida de dos o más variables y algunas de sus derivadas parciales.
Fijemos un entero k ≥ 1 y sea U un subconjunto abierto de Rn .

1
Saul Ancari Villca Ecuaciones diferenciales parciales

Definición 0.1. Una expresión de la forma

F Dk u(x), Dk−1 u(x), . . . , Du(x), u(x), x = 0

(x ∈ U ) (0.1)

sera llamado una ecuación diferencial parcial de orden k, donde

k k−1
F : Rn × Rn × · · · × Rn × R × U → R

esta dado, y
u:U →R
es desconocido.

Por lo general, la EDP viene sujeto a condiciones de frontera y condiciones iniciales.

Estaremos interesados en responder algunas de las siguientes preguntas:

(i) Existencia de soluciones: ¿Podemos encontrar soluciones explícitas? Si no, ¿ po-

demos probar existencia en un espacio adecuado? Si la solución se define en un
espacio "más débil "¿Podemos demostrar que en realidad es una solución clásica
(propiedades de regularidad)?

(ii) Unicidad de soluciones: ¿Cómo probamos que una determinada solución definida
en un el espacio adecuado es único? Para ello utilizaremos técnicas como el principio
del máximo, métodos de energía, etc.

(iii) Dependencia continua de las condiciones iniciales: Dado la condición inicial nos
gustaría que nuestras soluciones (nuevamente en un espacio apropiado) dependa
de forma continua de estos datos.

Por encontrar las soluciones u de (0.1) queremos decir, idealmente, obtener soluciones
simples y explícitas o, en su defecto, deducir la existencia y otras propiedades de las
soluciones.

Definición 0.2. (i) La ecuación diferencial parcial (0.1) es lineal si tiene la forma
X
aα (x)Dα u = f (x)
|α|≤k

para funciones dadas aα (|α| ≤ k), f . Esta EDP lineal es homogénea si f ≡ 0.

(ii) La EDP (0.1) es semilineal si tiene la forma

X
aα (x)Dα u + a0 Dk−1 u, . . . , Du, u, x = 0.

|α|=k

(iii) La EDP (0.1) es cuasilineal si tiene la forma

X
aα Dk−1 u, . . . , Du, u, x Dα u + a0 Dk−1 u, . . . , Du, u, x = 0

|α|=k

2
Saul Ancari Villca Ecuaciones diferenciales parciales

(iv) La EDP (0.1) es no lineal si depende de forma no lineal de las derivadas de orden superior.

No se conoce una teoría general sobre la solubilidad de todas las ecuaciones diferencia-
les parciales. Es extremadamente improbable que exista tal teoría, dada la rica variedad de
fenómenos físicos, geométricos y probabilísticos que pueden ser modelados por EDP. En
cambio, la investigación se centra en varias ecuaciones diferenciales parciales particulares
que son importantes para aplicaciones dentro y fuera de las matemáticas, con la esperanza
de que la comprensión de los orígenes de estas EDP pueda dar pistas sobre sus soluciones.
A continuación x ∈ U , donde U es un subconjunto abierto de Rn , y t ≥ 0. Además,
∇u = ∇x u = (ux1 , . . . , uxn ) es el gradiente u respecto a la variable x = (x1 , . . . , xn ).

Ecuaciones diferenciales parciales simples

Ecuaciones lineales

1. Ecuación del Laplace

n
X
∆u = uxi xi = 0.
i=1

2. Ecuación de transporte lineal

n
X
ut + bi uxi = 0.
i=1

3. Ecuación del calor (o difusión)

ut − ∆u = 0.

4. Ecuación de la onda
utt − ∆u = 0.

Ecuaciones no lineales

1. Ecuación de la Eikonal
|∇u| = 1.

2. Ecuación de Poisson no lineal

−∆u = f (u).

3. Ecuación de p−Laplace
div |∇u|p−2 ∇u = 0.

3
Saul Ancari Villca Ecuaciones diferenciales parciales

1. La ecuación de Laplace
Entre las más importantes de las ecuaciones diferenciales parciales se encuentran la
ecuación de Laplace

∆u = 0 (1.2)
y la ecuación de Poisson
−∆u = f (1.3)
En ambos, (1.2) and (1.3), U ⊂ Rn es abierto y u : U → R es la incógnita.
Pn En (1.3) la función
f : U → R esta dada. Recordemos que el Laplaciano de u es ∆u = i=1 uxi xi .

Definición 1.1. Una función u ∈ C 2 que satisface (1.2) sera llamada función armónica.

1.1. Solución fundamental

Motivación de la solución fundamental
Una buena estrategia para buscar soluciones de una EDP consiste en identificar prime-
ro algunas soluciones explícitas y luego, suponiendo que la ecuación es lineal, superponer
estas soluciones para obtener otras más complicadas. Además, cuando buscamos solucio-
nes explícitas, también es usual poner atención en alguna clase específica de soluciones
con ciertas propiedades de simetría. Ya que el operador Laplaciano es invariante bajo
rotaciones, en consecuencia, parece natural buscar soluciones radiales, esto es, funciones
de r = |x|. Por lo tanto, intentemos encontrar una solución u de la ecuación de Laplace
(1.2) en U = Rn , que tiene la forma

u(x) = v(r),
1/2
donde r = |x| = (x21 + · · · | x2n ) y v es seleccionado (si fuera posible) de tal forma que
∆u = 0. Primero notemos que para i = 1, . . . , n se verifica

∂r 1 2
−1/2 xi
= x1 + · · · + x2n 2xi = (x 6= 0).
∂xi 2 r
Luego
xi
uxi = v 0 (r) ,
r
x2i 1 x2i
 
00 0
uxi xi = v (r) 2 + v (r) − 3
r r r
para i = 1, . . . , n, y así
n−1 0
∆u = v 00 (r) + v (r).
r
Luego, ∆u = 0 sí y sólo sí
n−1 0
v 00 + v = 0.
r

4
Saul Ancari Villca Ecuaciones diferenciales parciales

Si v 0 6= 0, obtenemos
0 v 00 1−n
log (|v 0 |) = 0
= ,
v r
y por lo tanto v 0 (r) = a
rn−1
para alguna constante a. De esta forma, si r > 0, tenemos
(
b log r + c (n = 2)
v(r) = b
rn−2
+c (n ≥ 3),

donde b y c son constantes.

Estas consideraciones motivan la siguiente

Definición 1.2. La función

(
1
− 2π log |x| (n = 2)
Φ(x) := 1 1
n(n−2)α(n) |x|n−2
(n ≥ 3),

definida para x ∈ Rn \{0}, es la solución fundamental de la ecuacion de Laplace. Aquí α(n)

denota el volumen de la bola unitaria B(0, 1)..

1.2. Ecuación de Poisson

Por construcción la función x 7→ Φ(x) es armónica para x 6= 0. Si cambiamos el origen a
un nuevo punto y, tenemos que x 7→ Φ(x − y) es una función armónica de x, x 6= y. Ahora
tomemos f : Rn → R y notemos que la aplicación x 7→ Φ(x − y)f (y) (x 6= y) es armónica
para cada punto y ∈ Rn . Este razonamiento nos sugiere pensar que la convolución
ˆ
u(x) = Φ(x − y)f (y)dy
n
(R
1
´ (1.4)
− 2π log(|x − y|)f (y)dy (n = 2)
R2 ´
= 1 f (y)
n(n−2)α(n) Rn |x−y|n−2
dy (n ≥ 3)

resuelve la ecuación de Laplace (1.2). Sin embargo, este razonamiento es erróneo debido a
la singularidad en x = y. Debemos proceder con más cuidado al calcular ∆u.

Teorema 1.1. (Solución de la ecuación de Poisson). Supongamos que f ∈ Cc2 (Rn ) y definamos u
por (1.4). Entonces

(i) u ∈ C 2 (Rn )

(ii) −∆u = f en Rn .

En consecuencia, vemos que (1.4) nos proporciona una fórmula para una solución de
la ecuación de Poisson (1.3) en Rn .

5
Saul Ancari Villca Ecuaciones diferenciales parciales

Demostración. 1. Tenemos
ˆ ˆ
u(x) = Φ(x − y)f (y)dy = Φ(y)f (x − y)dy;
Rn Rn

luego ˆ  
u (x + hei ) − u(x) f (x + hei − y) − f (x − y)
= Φ(y) dy,
h Rn h
donde h 6= 0 y ei = (0, . . . , 1, . . . , 0), el 1 esta en el i−ésimo lugar. Ya que el soporte de f es
compacto y f ∈ C 1 (Rn ) tenemos que

f (x + hei − y) − f (x − y) ∂f
→ (x − y)
h ∂xi
uniformemente sobre Rn cuando h → 0, y de esta forma
ˆ
∂u ∂f
(x) = Φ(y) (x − y)dy (i = 1, . . . , n).
∂xi Rn ∂xi

Similarmente
ˆ
∂ 2u ∂ 2f
(x) = Φ(y) (x − y)dy (i, j = 1, . . . , n). (1.5)
∂xi ∂xj Rn ∂xi ∂xj

Como la expresión del lado derecho de (1.5) es continua en la variable x, vemos que
u ∈ C 2 (Rn ).
2. Dado que Φ(x) → ∞ cuando x → 0, necesitaremos cálculos posteriores para aislar
esta singularidad dentro de una bola pequeña. Fijemos  > 0, entonces
ˆ ˆ
∆u(x) = Φ(y)∆x f (x − y)dy + Φ(y)∆x f (x − y)dy
por (1.5) B(0,) Rn −B(0,) (1.6)
=: I + J

Luego
ˆ (
C2 | log | (n = 2)
|I | ≤ C D2 (f ) L∞ (Rn ) |Φ(y)|dy ≤ (1.7)
B(0,) C2 (n ≥ 3).
Integrando por partes obtenemos
ˆ
J = Φ(y)∆y f (x − y)dy
Rn −B(0,)
ˆ
=− h∇Φ(y), ∇y f (x − y)i dy
Rn −B(0,) (1.8)
ˆ
∂f
+ Φ(y) (x − y)dS(y)
∂B(0,) ∂ν
=: K + L ,

6
Saul Ancari Villca Ecuaciones diferenciales parciales

ν denota el vector normal apuntando hacia adentro a lo largo de ∂B(0, ). Podemos ver
fácilmente que
ˆ (
C| log | (n = 2)
|L | ≤ k∇f kL∞ (Rn ) |Φ(y)|dS(y) ≤ (1.9)
∂B(0,) C (n ≥ 3)

3. Integramos por partes el término K , para obtener

ˆ ˆ
∂Φ
K = ∆Φ(y)f (x − y)dy − (y)f (x − y)dS(y)
Rn −B(0,) ∂B(0,) ∂ν
ˆ
∂Φ
=− (y)f (x − y)dS(y)
∂B(0,) ∂ν

ya que Φ es armónica fuera del origen. Ahora ∇Φ(y) = nα(n) −1 y

|y|n
(y 6= 0) y ν = −y
|y|
= − y
sobre ∂B(0, ). Como consecuencia ∂Φ ∂ν
(y) = hν, ∇Φ(y)i = nα(n)1
n−1 sobre ∂B(0, ). Puesto

que nα(n)εn−1 es el área de la superficie de la esfera ∂B(0, ), tenemos

ˆ
1
K = − f (x − y)dS(y)
nα(n)n−1 ∂B(0,)
(1.10)
=− f (y)dS(y) → −f (x) cuando  → 0
∂B(x,)

4. Combinando (1.6)-(1.10) y haciendo  → 0, concluimos que −∆u(x) = f (x), como

habíamos afirmado

1.3. Fórmula del valor medio

Ahora consideraremos un conjunto abierto U ⊂ Rn y supondremos que u es función
armónica en U . Nuestro primer objetivo, y uno de los principales en esta sección, es el de
obtener la propiedad del valor medio para funciones armónicas. Esta propiedad establece
que u(x) es igual al promedio de u sobre la esfera ∂B(x, r) y también al promedio de u
sobre la bola B(x, r), sobre cualquier bola B(x, r) ⊂ U . Aunque las fórmulas obtenidas
no son explícitas para u, a partir de ellas podremos obtener otros resultados importantes,
como veremos más adelante.
Usaremos la siguiente notación:
ˆ
1
f (y)dy := f (y)dy
B(x,r) α(n)rn B(x,r)
= promedio de f sobre la bola B(x, r),

ˆ
1
f (y)dS(y) := f dS
∂B(x,r) nα(n)rn−1 ∂B(x,r)

= promedio de f sobre la esfera ∂B(x, r),

7
Saul Ancari Villca Ecuaciones diferenciales parciales

donde
π n/2
α(n) = volumen de la bola unitaria B(0, 1) en Rn =
Γ n2 + 1

y
nα(n) = área de la esfera unitaria ∂B(0, 1) en Rn .

Teorema 1.2. (Fórmula del valor medio). Si u ∈ C 2 (U ) es armónica, entonces

u(x) = u(y)dS(y) = u(y)dy (1.11)

∂B(x,r) B(x,r)

para toda bola B(x, r) ⊂ U .

Demostración. 1. Definamos

φ(r) := u(y)dS(y)
∂B(x,r)

= u(x + rz)dS(z).
y = x + rz ∂B(0,1)

Entonces
φ0 (r) = h∇u(x + rz), zidS(z).
∂B(0,1)

Luego,  
0 y−x
φ (r) = ∇u(y), dS(y)
∂B(x,r) r
ˆ
1 ∂u
= n−1
(y)dS(y)
y − x nα(n)r ∂B(x,r) ∂ν
ν=
r
ˆ
1
= ∆u(y)dy.
formulas de Green nα(n)rn−1 B(x,r)

Entonces
ˆ
0 r
φ (r) = ∆u(y)dy (1.12)
n B(x,r)
= 0 (1.13)
u es armónica

De donde φ es constante y

φ(r) = lı́m φ(t) = lı́m u(y)dS(y) = u(x).

t→0 t→0 ∂B(x,t)

Por tanto
u(y)dS(y) = u(x). (1.14)
∂B(x,r)

8
Saul Ancari Villca Ecuaciones diferenciales parciales

2. Por otro lado, usando coordenadas polares tenemos

ˆ ˆ r ˆ 
u(y)dy = u(y)dS(y) ds
B(x,r) 0 ∂B(x,s)
ˆ r
= u(x) nα(n)sn−1 ds = α(n)rn u(x)
por (1.14) 0

Teorema 1.3. Si u ∈ C 2 (U ) satisface

u(x) = u(y)dS(y)
∂B(x,r)

para toda bola B(x, r) ⊂ U , entonces u es armónica.

Demostración. La prueba es por contradicción: Si ∆u 6≡ 0 tenemos que existe x0 ∈ U y
r0 > 0 tal que B(x0 , r0 )) ⊂ U y

∆u > 0 sobre B(x0 , r0 ). (1.15)

ffl
Considerando que φ(r) = ∂B(x,r) u(y)dS(y) = constante, tenemos

r
0 = φ0 (r) = ∆u(y)dy > 0.
por (1.12) y (1.15) n B(x,r)

Lo que es una contradicción.

1.4. Propiedades de las funciones armónicas

A continuación, presentamos una secuencia de resultados interesantes sobre funciones
armónicas, todas estas basadas en la fórmula del valor medio. En lo que sigue, U ⊂ Rn es
abierto y limitado.

1.4.1. Principio del máximo fuerte, unicidad

Teorema 1.4. (Principio máximo fuerte). Supongamos que u ∈ C 2 (U ) ∩ C(U ) es armónica en U .
(i) Entonces
máx u = máx u.
U ∂U

(ii) Además, si U es conexo y existe un punto x0 ∈ U tal que

u (x0 ) = máx u
U

entonces
u es constante en U.

9
Saul Ancari Villca Ecuaciones diferenciales parciales

La afirmación (i) es el principio del máximo para la ecuación de Laplace y (ii) es el

principio del máximo fuerte. Reemplazando u por −u, obtenemos afirmaciones similares con
“ mı́n ” reemplazada por “ máx ”.
Demostración. Supongamos que existe x0 ∈ U con u (x0 ) = M := máxU u. Entonces, para
0 < r < dist(x0 , ∂U ), la propiedad del valor medio nos da

M = u (x0 ) = u(y)dy ≤ M,
B(x0 ,r)
1=1
B(x0 ,r)

luego
(M − u(y)) dy = 0.
B(x0 ,r)

De M − u ≥ 0 en U , podemos ver que u(y) = M para todo y ∈ B(x, r). Por tanto,
{x ∈ U | u(x) = M } es abierto y cerrado en U , luego igual a U ya que U es conexo. Esto
prueba la afirmación (ii), como consecuencia tenemos (i).

Ejercicio. Suponga que U es conexo y que u ∈ C 2 (U ) ∩ C(U ) satisface

(
∆u = 0 en U
u = g sobre ∂U,

donde g ≥ 0. Entonces u es positivo en U si g no se anula sobre ∂U .

Una aplicación importante del principio del máximo es la unicidad de las soluciones
de la ecuación de Poisson con ciertas condiciones en la frontera.
Teorema 1.5. (Unicidad). Sea g ∈ C(∂U ), f ∈ C(U ). Existe a lo más una solución u ∈ C 2 (U ) ∩
C(U ) del problema con valores en la frontera
(
−∆u = f en U
(1.16)
u = g sobre ∂U.

e son dos soluciones de (1.16) tenemos que wj := (−1)j (u−e

Demostración. Si u y u u), j = 0, 1,
son soluciones de (
−∆wj = 0 en U
wj = 0 sobre ∂U.
Luego, aplicando el Teorema 1.4(i) a w0 y w1 tenemos

máx w0 = máx w0 = 0 = máx w1 = máx w1 = − mı́n w0 ,

U ∂U ∂U U U

por tanto
u−u
e = 0.

10
Saul Ancari Villca Ecuaciones diferenciales parciales

1.4.2. Regularidad
Comenzaremos introduciendo algunas definiciones y resultados sobre el regularizador.
A continuación demostramos que si u ∈ C 2 es armónica, entonces necesariamente u ∈
C ∞ . Por tanto, las funciones armónicas son automáticamente infinitamente diferenciables. Esta
afirmación se llama teorema de regularidad.
Para  > 0 escribamos

U := {x ∈ U | dist(x, ∂U ) > }.

Definición 1.3. (i) Definamos η ∈ C ∞ (Rn ) por

(  
C exp |x|21−1 si |x| < 1
η(x) :=
0 si |x| ≥ 1,
´
elegimos la constante C > 0 de tal forma que Rn
ηdx = 1.

(ii) Para cada  > 0, sea

1 x
η (x) := n η
 
Llamamos a η el regularizador estándar. La función η es suave y satisface
ˆ
η dx = 1, spt (η ) = B(0, ).
Rn

Definición 1.4. Si f : U → R es localmente integrable, definimos su regularizante

ˆ ˆ

f (x) := η (x − y)f (y)dy = η (y)f (x − y)dy
U B(0,)

para x ∈ U

Teorema 1.6 (Propiedades de los regularizantes). (i) f  ∈ C ∞ (U ).

(ii) f  → f c.t.p. cuando  → 0.

(iii) Si f ∈ C(U ), entonces f  → f uniformemente sobre subconjuntos compactos de U .

(iv) Si 1 ≤ p < ∞ y f ∈ Lploc (U ), entonces f  → f en Lploc (U ).

Teorema 1.7. (Regularidad). Si u ∈ C(U ) satisface la propiedad del valor medio (1.11) para cada
bola B(x, r) ⊂ U , entonces
u ∈ C ∞ (U ).

Demostración. Sea η el regularizador estándar, observemos que η es una función radial.

Sea u el regularizante de u.

11
Saul Ancari Villca Ecuaciones diferenciales parciales

Probaremos que u es suave al demostrar que u ≡ uε sobre Uε . En efecto, si x ∈ Uε ,

entonces ˆ
ε
u (x) = ηε (x − y)u(y)dy
U
ˆ  
1 |x − y|
= n η u(y)dy
ε B(x,ε) ε
ˆ ˆ
1 ε r

= n η udS dr
ε 0 ε ∂B(x,r)
ˆ ε  
1 r
= n
u(x) η nα(n)rn−1 dr
por (1.11) ε 0 ε
ˆ
= u(x) ηε dy = u(x).
B(0,ε)

De esta forma uε ≡ u en U , y así u ∈ C ∞ (Ue ) para cada ε > 0.

1.4.3. Teorema de Liouville

Aquí probaremos que no existen funciones armónicas no triviales, limitadas y definidas
en todo Rn .
Teorema 1.8 (Teorema de Liouville). Sea u : Rn → R una función armónica. Si u es limitada
(Existe M > 0 tal que |u(x)| ≤ M ) para todo x ∈ Rn , entonces u es constante.
Demostración. Sea x0 ∈ Rn un punto arbitrario. Como u es armónica, es claro que uxi , es
también armónica para i = 1, 2, ...n. Aplicando la fórmula del valor medio a uxi con r > 0,
tenemos

|uxi (x0 )| = | uxi dx|

B(x0 , r2 )
n ˆ
2
=| uνi dS|
α(n)rn ∂B(x0 , r2 )
2n 2n
≤ kukL∞ (∂B(x0 , r2 )) ≤ M.
r r
Haciendo r → ∞, concluimos que uxi (x0 ) = 0 para todo x0 ∈ Rn , i = 1, 2, ...n. Esto prueba
que u es constante.
Teorema 1.9 (Fórmula de representacíon). Sea f ∈ Cc2 (Rn ), n ≥ 3. Entonces cualquier
solución limitada de
−∆u = f en Rn
tiene la forma ˆ
u(x) = Φ(x − y)f (y)dy + C
Rn
para alguna constante C.

12
Saul Ancari Villca Ecuaciones diferenciales parciales

´
Demostración. Como para n ≥ 3, Φ(x) → 0 cuando |x| → ∞, ū(x) := Rn Φ(x − y)f (y)dy
es una solución limitada de −∆u = f en Rn . Si u es otra solución limitada, w = u − ū es
constante, de acuerdo al teorema de Liouville.

1.5. Función de Green

Suponga ahora que U es un subconjunto abierto, limitado de Rn con ∂U de clase C 1 . El
propósito de esta parte es presentar una fórmula general para la solución de la ecuación
de Poisson

−∆u = f en U,

sujeta a la condición de frontera

u = g sobre ∂U.

Derivación de la fórmula de Green

Sea u ∈ C 2 (Ū ) una función arbitraria. Sea x ∈ U y ε > 0 muy pequeño tal que B(x, ε) ⊂
U . Aplicando la fórmula de Green a u(y) y Φ(y − x), sobre la región Vε := U \ B(x, ε),
tenemos
ˆ
[u(y)∆Φ(y − x) − Φ(y − x)∆u(y)]dy
Vε
ˆ
= [u(y)h∇y Φ(y − x), νi − Φ(y − x)h∇u, νi]dS(y),
∂Vε

donde ν es el vector unitario exterior sobre ∂Vε . Recordemos que ∆Φ(x − y) = 0 para
x 6= y. Entonces
ˆ ˆ
− Φ(y − x)∆u(y)dy = [u(y)h∇y Φ(y − x), νi − Φ(y − x)h∇u, νi]dS(y). (1.17)
Vε ∂Vε

Además observamos que

ˆ  ˆ
 

 Φ(y − x)h∇u, νidS(y) ≤ Φ(y − x)|∇u|dS(y)
∂B(x,ε) ∂B(x,ε)
ˆ
= Φ(z)|∇u|dS(z)
∂B(0,ε)

≤ C máx |Φ|εn−1 → 0, cuando ε → 0.

∂B(0,ε)

Por otro lado

ˆ ˆ  
y−x y−x
u(y)h∇y Φ(y − x), νidS(y) = u(y) − ,− dS(y)
∂B(x,ε) ∂B(x,ε) nα(n)|y − x|n |y − x|
ˆ
1
= u(y)dS(y) → u(x), cuando ε → 0.
nα(n)εn−1 ∂B(x,ε)

13
Saul Ancari Villca Ecuaciones diferenciales parciales

Observe que ∂Vε = ∂U ∪ ∂B(x, ε). Haciendo ε → 0 en (1.17) obtenemos

ˆ
u(x) = [Φ(y − x)h∇u, νi − u(y)h∇y Φ(y − x), νi]dS(y) (1.18)
ˆ
∂U

− Φ(y − x)∆u(y)dy.
U

Esta identidad es valido para cualquier punto x ∈ U y cualquier función u ∈ C 2 (Ū ).

Remarca 1. Si los valores de ∆u son conocidos en U y los valores de u, ∇u son conocidos en ∂U ,

entonces las función u esta completamente determinada en Ū .
La idea ahora es remover el termino h∇u, νi. Para x ∈ U fijo, introducimos una función
corrector ϕx = φx (y), que es una solución del siguiente problema

∆φx = 0 en U
φx = Φ(y − x) sobre ∂U.
Usando la fórmula de Green para integrales, obtenemos
ˆ ˆ
x
− φ ∆udy = [uh∇φx , νi − φx (y)h∇u, νi]dS(y) (1.19)
U ˆ∂U
= [uh∇φx , νi − Φ(y − x)h∇u, νi]dS(y)
∂U
Reemplazando (1.19) en (1.18) conseguimos
ˆ
u(x) = − u(y)[h∇y Φ(y − x), νi − h∇φx , νi]dS(y) (1.20)
ˆ∂U
− [Φ(y − x) − φx (y)]∆u(y)dy.
U

Ahora supongamos que u ∈ C 2 (Ū ) resuelve


−∆u = f en U
(1.21)
u=g sobre ∂U,
para funciones f , g continuas. Entonces, usando (1.20), u tiene la siguiente expresión
ˆ
u(x) = − g(y)[h∇y Φ(y − x) − ∇φx (y), νi]dS(y) (1.22)
ˆ∂U
+ [Φ(y − x) − φx (y)]f (y)dy.
U
Motivado a esto, definimos la función de Green.
Definición 1.5. La función de Green para la región U es definida por
G(x, y) := Φ(y − x) − φx (y) (x, y, x 6= y).
Teorema 1.10 (Fórmula de representación usando la función de Green). Si u ∈ C 2 (Ū )
resuelve el problema (1.21), entonces
ˆ ˆ
u(x) = − g(y)h∇y G(x, y), ν(y) + f (y)G(x, y)dy.
∂U U

14
Saul Ancari Villca Ecuaciones diferenciales parciales

2. La ecuación del calor

A continuación estudiaremos la ecuación del calor.

ut − ∆u = 0 (1)

y la ecuación no-homogénea del calor

ut − ∆u = f, (2)

sujeto a condiciones iniciales y de borde adecuadas. Aquí t > 0 y x ∈ U , donde U ⊂ Rn es

abierto. La incógnita es u : U × [0, ∞) → R, u = u(x, t), y el Laplaciano
Pn ∆ es considerado
con respecto a la variable x = (x1 , . . . , xn ): ∆u = ∆x u = i=1 uxi xi . En (2) la función
f : U × [0, ∞) → R esta dada.
Una gran cantidad de resultados sobre funciones armónicas tienen un resultado análo-
go (pero más complicado) en el contexto de la ecuación del calor. Por ello, organizaremos
nuestro estudio de la ecuación del calor siguiendo las ideas surgidas en el estudio de la
ecuación de Laplace.

2.1. Solución fundamental

Motivación de la solución fundamental
En la Sección 1.1 pudimos observar que un paso importante en el estudio de una EDP
es a menudo encontrar algunas soluciones específicas.
Observemos que la ecuación del calor involucra una derivada con respecto al tiempo
t, pero dos derivadas con respecto a las variables xi (i = 1, . . . , n). En consecuencia, si u
resuelve 1, entonces también lo hace u(λx, λ2 t) para λ ∈ R. Este escalamiento indica que
2
el radio rt (r = |x|) es importante para la ecuación de calor y sugiere que busquemos una
2
solución de (1) que tenga la forma u(x, t) = v( rt ) = v ( |x|t ) (t > 0, x ∈ Rn ), para alguna
2

función v aún por determinar.

Aunque este enfoque eventualmente conduce a lo que queremos (ver el Ejercicio 8), es
más rápido buscar una solución u teniendo la estructura especial

1 x
u(x, t) = v β (x ∈ Rn , t > 0) , (3)
tα t
donde las constantes α, β and the function v : Rn → R serán encontradas. Llegamos a (3)
si buscamos una solución u de la ecuación de calor invariante bajo la dilatación

u(x, t) 7→ λα u λβ x, λt .

Es decir, pedimos que

u(x, t) = λα u λβ x, λt

para todo λ > 0, x ∈ Rn , t > 0. Escribamos λ = t−1 , de esta forma obtenemos (3) para
v(y) := u(y, 1).

15
Saul Ancari Villca Ecuaciones diferenciales parciales

Reemplacemos (3) en (1) y luego calculamos

αt−(α+1) v(y) + βt−(α+1) hy, ∇v(y)i + t−(α+2β) ∆v(y) = 0 (4)

para y := t−β x. Para transformar (4) en una expresión que involucre solamente la variable
y tomamos β = 12 . Entonces los términos con t son idénticos, por lo que (4) se reduce a
1
αv + hy, ∇vi + ∆v = 0. (5)
2
Podemos simplificar aún más si suponemos que v es radial; es decir, v(y) = w(|y|) para
algún w : R → R. Entonces (5) se convierte en
1 n−1 0
αw + rw0 + w00 + w = 0,
2 r
para r = |y|, 0 = d
dr
. Ahora, si tomamos α = n2 , simplificamos

0 1
rn−1 w0 + (rn w)0 = 0.
2
De esta forma
1
rn−1 w0 + rn w = a
2
para alguna constante a. Asumiendo que lı́mr→∞ w = w0 = 0, tenemos a = 0, de donde
1
w0 = − rw.
2
Luego, para alguna constante
r2
w = be− 4 . (6)
|x|2
b − 4t
Combinando (3), (6) y nuestras elecciones para α, β, concluimos que tn/2
e resuelve la
ecuación del calor (1). Este cálculo motiva la siguiente
Definición 2.1. La función
( |x|2
1
(4πt)n/2
e− 4t (x ∈ Rn , t > 0)
Φ(x, t) :=
0 (x ∈ Rn , t < 0)

es llamada la solución fundamental de la ecuación del calor.

Podemos notar que Φ es singular en el punto (0, 0). Algunas veces escribiremos
Φ(x, t) = Φ(|x|, t) para enfatizar que la solución fundamental es radial en la variable
x. El motivo por el que elegimos la constante (4π)−n/2 se debe al siguiente
Lema 2.1. (Integral de la solución fundamental). Para cada tiempo t > 0,
ˆ
Φ(x, t)dx = 1.
Rn

16
Saul Ancari Villca Ecuaciones diferenciales parciales

Demostración. Calculemos
ˆ ˆ
1 |x|2
− 4t
Φ(x, t)dx = e dx
(4πt)n/2 Rn
Rn
ˆ
1 2
= n/2 e−|z| dz
π Rn
n ˆ
1 Y ∞ −zi2
= n/2 e dzi = 1.
π i=1 −∞

2.2. Problema de valor inicial

Ahora emplearemos Φ para encontrar una solución del problema de valor inicial (o
Cauchy) (
ut − ∆u = 0 en Rn × (0, ∞)
(7)
u = g sobre Rn × {t = 0}.
Observemos que la función (x, t) 7→ Φ(x, t) resuelve la ecuación de calor menos en la
singularidad (0, 0), y de igual forma (x, t) 7→ Φ(x − y, t) para cada y ∈ Rn fijo. Luego la
convolución ˆ
u(x, t) = Φ(x − y, t)g(y)dy
Rn ˆ (8)
1 |x−y|2
− 4t n
= e g(y)dy (x ∈ R , t > 0)
(4πt)n/2 Rn
también es una solución.

Teorema 2.1. (Solución del problema de valor inicial). Asumiremos que g ∈ C (Rn ) ∩ L∞ (Rn ),
y definamos u por (8). Entonces

(i) u ∈ C ∞ (Rn × (0, ∞))

(ii) ut (x, t) − ∆u(x, t) = 0 (x ∈ Rn , t > 0),

(iii) lı́m u(x, t) = g (x0 ) para todo punto x0 ∈ Rn .

(x,t)→(x0 ,0)
x∈Rn ,t>0

|x|2
Demostración. 1. Ya que la función tn/21
e− 4t es infinitamente diferenciable, con derivadas
de todos los ordenes uniformemente limitadas, sobre Rn ×[δ, ∞) para todo δ > 0, podemos
ver que u ∈ C ∞ (Rn × (0, ∞)). Es más
ˆ
ut (x, t) − ∆u(x, t) = [(Φt − ∆x Φ) (x − y, t)] g(y)dy
Rn
=0 (x ∈ Rn , t > 0) ,

17
Saul Ancari Villca Ecuaciones diferenciales parciales

ya que Φ resuelve la ecuación de calor.

2. Fijemos x0 ∈ Rn ,  > 0. Elijamos δ > 0 tal que
  
g(y) − g x0  <  si y − x0  < δ, y ∈ Rn .


(9)
Entonces, si |x − x0 | < 2δ , tenemos, según el lema 2.1,
ˆ 
 0

   0

 
u(x, t) − g x  = 
 n Φ(x − y, t) g(y) − g x dy 

ˆ R

Φ(x − y, t) g(y) − g x0  dy


≤
B(x0 ,δ)
ˆ

Φ(x − y, t) g(y) − g x0  dy


+
Rn −B(x0 ,δ)

=:I + J.
Ahora, ˆ
I≤ Φ(x − y, t)dy = ,
Rn
δ
debido a (9) y el lema 2.1. Además, si |x − x0 | ≤ 2
y |y − x0 | ≥ δ, entonces

y − x0  ≤ |y − x| + δ ≤ |y − x| + 1 y − x0  .
   
2 2
De esta forma |y − x| ≥ 12 |y − x0 |. Como consecuencia
ˆ
J ≤ 2kgkL∞ Φ(x − y, t)dy
Rn −B(x0 ,δ)
ˆ
C |x−y|2
≤ n/2 e− 4t dy
t n 0
ˆR −B(x ,δ)
C |y−x0 |2
≤ n/2 e− 16t dy
t Rn −B(x0 ,δ)
ˆ
|z|2
− 16
= C √
e dz → 0 cuando t → 0+ .
√
0
y−x =z t R −B (0,δ/ t)
n

δ
Por tanto, si |x − x0 | < 2
y t > 0 es suficientemente pequeño,

u(x, t) − g x0  < 2.



2.3. Problema no-homogéneo

Ahora veremos el problema de valor inicial no-homogéneo
(
ut − ∆u = f en Rn × (0, ∞)
(10)
u = 0 sobre Rn × {t = 0}.

18
Saul Ancari Villca Ecuaciones diferenciales parciales

¿Cómo podemos obtener una fórmula para la solución? Recordemos la motivación que
nos llevo a (8), podemos observar que la función (x, t) 7→ Φ(x − y, t − s) es la solución de
la ecuación del calor (para y ∈ Rn , 0 < s < t). Fijemos s, la función
ˆ
u = u(x, t; s) = Φ(x − y, t − s)f (y, s)dy
Rn

resuelve (
ut (·; s) − ∆u(·; s) = 0 en Rn × (s, ∞)
(11)
u(·; s) = f (·, s) sobre Rn × {t = s},
que es un problema de valor inicial de la forma (7), con tiempo inicial t = 0 reemplazado
por t = s y g reemplazado por f (·, s). Observemos que u(·; s) no es una solución de (10).
Sin embargo, el principio de Duhamel asegura que podemos construir una solución de
(10) a partir de las soluciones de (11), integrando con respecto a s. La idea es considerar
ˆ t
u(x, t) = u(x, t; s)ds (x ∈ Rn , t ≥ 0) .
0

Reescribiendo, tenemos
ˆ tˆ
u(x, t) = Φ(x − y, t − s)f (y, s)dyds
Rn
ˆ0
t ˆ (12)
1 |x−y|2
− 4(t−s)
= e f (y, s)dyds,
0 (4π(t − s))n/2 Rn
n
para x ∈ R , t > 0.
A continuación confirmaremos que la fórmula (12) funciona.
Teorema 2.2. (Solución del problema no-homogéneo). Supongamos que f ∈ C12 (Rn × [0, ∞)) y,
además, que tiene soporte compacto. Definamos u por (12). Entonces
(i) u ∈ C12 (Rn × (0, ∞)),
(ii) ut (x, t) − ∆u(x, t) = f (x, t) (x ∈ Rn , t > 0),
y
(iii) lı́m u(x, t) = 0 para todo punto x0 ∈ Rn .
(x,t)→(x0 ,0)
x∈Rn ,t>0

Demostración. 1. Usando cambio de variable, reescribimos

ˆ tˆ
u(x, t) = Φ(y, s)f (x − y, t − s)dyds.
0 Rn

Ya que f ∈ C12 (Rn × [0, ∞)) tiene soporte compacto y Φ = Φ(y, s) es suave cerca de
s = t > 0, tenemos ˆ tˆ
ut (x, t) = Φ(y, s)ft (x − y, t − s)dyds
0 ˆ Rn
(13)
+ Φ(y, t)f (x − y, 0)dy
Rn

19
Saul Ancari Villca Ecuaciones diferenciales parciales

y
ˆ tˆ
uxi xj (x, t) = Φ(y, s)fxi xj (x − y, t − s)dyds (i, j = 1, . . . , n) (14)
0 Rn
Así ut , uxi xj , y de igual forma u, ∇x u, pertenecen a C (Rn × (0, ∞)).
2. Ahora calculamos
ut (x, t) − ∆u(x, t)
ˆ tˆ   
∂
= Φ(y, s) − ∆x f (x − y, t − s) dyds
por (13) y (14) 0 Rn ∂t
ˆ
+ Φ(y, t)f (x − y, 0)dy
Rn
ˆ tˆ   
∂ (15)
= Φ(y, s) − − ∆y f (x − y, t − s) dyds
Rn ∂s
ε
ˆ εˆ   
∂
+ Φ(y, s) − − ∆y f (x − y, t − s) dyds
∂s
ˆ0 Rn

+ Φ(y, t)f (x − y, 0)dy

Rn
=: Iε + Jε + K
Ahora, ˆ εˆ


|Jε | ≤ kft kL∞ + D2 f L∞ Φ(y, s)dyds ≤ C, (16)
0 Rn
por el lema 2.1. Integrando por partes podemos encontrar
ˆ t ˆ   
∂
Iε = − ∆y Φ(y, s) f (x − y, t − s)dyds
∂s
ˆR
ε n

+ Φ(y, ε)f (x − y, t − ε)dy

ˆ
n
R (17)
− Φ(y, t)f (x − y, 0)dy
ˆ R n

= Φ(y, ε)f (x − y, t − ε)dy − K,

Rn

ya que Φ es resuelve la ecuación del calor. Combinando (15) - (17), afirmamos

ˆ
ut (x, t) − ∆u(x, t) = lı́m Φ(y, ε)f (x − y, t − ε)dy
ε→0 Rn
= f (x, t) (x ∈ Rn , t > 0) ,
el limite cuando  → 0 lo calculamos similarmente como hicimos en la demostración del
Teorema 2.1. Finalmente notemos que
ˆ tˆ
ku(·, t)kL∞ ≤ kf kL∞ Φ(y, t − s)dyds
0 Rn
= tkf kL∞ → 0 cuando t → 0.
por lema 2.1

20
Saul Ancari Villca Ecuaciones diferenciales parciales

2.3.1. Solución del problema no homogéneo con condiciones iniciales generales

Por supuesto que podemos combinar los Teoremas 2.1 y 2.2 para descubrir que
ˆ ˆ tˆ
u(x, t) = Φ(x − y, t)g(y)dy + Φ(x − y, t − s)f (y, s)dyds
Rn 0 Rn

es, bajo las hipótesis de g y f como arriba, una solución de

(
ut − ∆u = f en Rn × (0, ∞)
u = g sobre Rn × {t = 0}

2.4. Propiedades de la solución

Sea U ⊂ Rn un subconjunto abierto y limitado, y T > 0 fijo.
Definición 2.2.
(i) Definimos el cilindro parabólico

UT := U × (0, T ].

(ii) El borde parabólico de UT es

ΓT := ŪT − UT .

Notemos que UT incluye la tapa superior U × {t = T } y el borde parabólico ΓT

comprende la tapa inferior y los lados verticales de U × [0, T ], pero no la tapa superior.


Teorema 2.3. (Principio máximo). Supongamos que u ∈ C12 (UT ) ∩ C U T tal que

ut − ∆u ≤ 0,

entonces
máx u = máx u.
UT ΓT

Demostración. 1. Suponga que ut − ∆u < 0 en UT . Para ε > 0, defina

UT −ε = {(x, t) : x ∈ U, 0 < t ≤ T − ε}

Sea (x0 , t0 ) ∈ ŪT −ε tal que u(x, t0 ) = máxŪT −ε u.

Si (x0 , t0 ) ∈ UT −ε \ U × {T − ε}, entonces ut = 0 y ∆u ≤ 0, que induce una contradic-
ción.

Si (x0 , t0 ) ∈ U × {T − ε}, entonces ut ≥ 0 y ∆ ≤ 0. Esto también induce una

contradicción.

21
Saul Ancari Villca Ecuaciones diferenciales parciales

Por tanto (x0 , t0 ) ∈ ΓT −ε , así

máx u = máx u ≤ máx u,
ŪT −ε ΓT −ε ΓT

haciendo ε → 0, concluimos que

máx u ≤ máx u.
ŪT ΓT

2. En el caso general, para c > 0, consideremos

v(x, t) = u(x, t) − ct,
luego
vt − ∆v = ut − ∆u − c < 0.
De la primera parte de la prueba sabemos que
máx u = máx(v + ct)
ŪT ŪT

≤ máx(v) + cT
ŪT

≤ máx(v) + cT
ΓT

≤ máx(u) + 2cT.
ΓT

Haciendo c → 0, obtenemos que

máx u = máx u.
ŪT ΓT

Una aplicación importante del principio máximo es la siguiente afirmación.

Teorema 2.4. (Unicidad).

Sea g ∈ C (ΓT ), f ∈ C (UT ). Entonces existe a lo más una solución
2
u ∈ C1 (UT ) ∩ C U T del problema de valores iniciales y de frontera
(
ut − ∆u = f en UT
(18)
u = g sobre ΓT .

e son dos soluciones de (18) tenemos que wj := (−1)j (u − u

Demostración. Si u y u e), j = 0, 1,
son soluciones de (
(wj )t − ∆wj = 0 en UT
wj = 0 sobre ΓT .
Luego, aplicando el Teorema 2.3 a w0 y w1 tenemos
máx w0 = máx w0 = 0 = máx w1 = máx w1 = − mı́n w0 ,
UT ΓT ΓT UT UT

por tanto
u−u
e = 0.

22
Saul Ancari Villca Ecuaciones diferenciales parciales

Teorema 2.5 (Principio del máximo en U = Rn ). Sea u1 ∈ C12 (Rn × (0, T ]) ∩ C(Rn × [0, T ])
tal que

ut − ∆u ≤ 0 en Rn × (0, T ]
u=g sobre Rn × {t = 0}
2
y u(x, t) ≤ Aea|x| para todo x ∈ Rn y 0 ≤ t ≤ T , donde A, a > 0 son constantes. Entonces
sup u = sup g.
Rn ×[0,T ] Rn

Demostración. Primero supongamos que

4aT < 1.
Consideremos ε > 0, tal que
4a(T + ε) < 1. (19)
Para un punto fijo y ∈ Rn , µ > 0, definimos
µ |x−y|2
v(x, t) : = u(x, t) − n e 4(T +ε−t) (x ∈ Rn , t > 0)
(T + ε − t) 2

= u(x, t) − µK.
Recordemos que la solución fundamental es
1 − |x−y|2
Φ(x, y, t) = n e
4t
t2
y observe que
K(x, y, t) = Φ(ix, iy, T + ε − t),
así K es solución de la ecuación del calor. Entonces
vt − ∆v = ut − ∆u ≤ 0 en Rn × (0, T ].
Fijemos r > 0 y sea U := B 0 (y, r), UT = B 0 (y, r)×[0, T ]. Aplicando el principio de máximo
anterior, tenemos
máx v ≤ máx v.
ŪT ΓT

Si x ∈ Rn
µ |x−y|2
v(x, 0) = u(x, 0) − n e 4(T +ε) ≤ u(x, 0) = g(x).
(T + ε) 2
Si |x − y| = r, 0 ≤ t ≤ T , entonces
µ r2
v(x, t) = u(x, t) − n e
4(T +ε−t)
(T + ε − t) 2
2 µ r2
≤ Aea|x| − n e
4(T +ε−t)
(T + ε − t) 2
2 µ r2
≤ Aea(|y|+r) − n e
4(T +ε) .
(T + ε) 2

23
Saul Ancari Villca Ecuaciones diferenciales parciales

Ahora de acuerdo a (19), para algún γ > 0

1
= a + γ.
4(T + ε)
Así, continuando con los cálculos, tenemos
2 n 2 (a+γ)
v(x, t) ≤ Aea(|y|+r) − µ(4(a + γ)) 2 er → −∞ cuando r → ∞
≤ sup g,
Rn

para r suficientemente grande. Por lo tanto

v(x, t) ≤ sup g
Rn

para todo y ∈ Rn , [0 ≤ t ≤ T ]. Como

µ |x−y|2
v(x, t) = u(x, t) − n e
4(T +ε−t) ≤ sup g,
(T + ε − t) 2 Rn

haciendo µ → 0, concluimos que

u(x, t) ≤ sup g.
Rn

En el caso general donde 4aT < 1 falla, aplicamos el resultado anterior sobre los intervalos
1
[0, T1 ], [T1 , 2T1 ], etc., para T1 = 8a .

Teorema 2.6. Sea g ∈ C(Rn ), f ∈ C(Rn × [0, T ]). Entonces el problema


ut − ∆u = f enRn × (0, T )
u=g sobre Rn × {t = 0}

tiene a lo más una solución u ∈ C12 (Rn × (0, T ]) ∩ C(Rn × [0, T ]), satisfaciendo
2
|u(x, t)| ≤ Aea|x| (x ∈ Rn , 0 ≤ t ≤ T )

para constantes A, a>0.

Demostración. Ejercicio.

Ejercicios
1. En Rn , pruebe que la ecuación de Laplace ∆u = 0 es invariante bajo rotaciones; esto
es, si O es una matriz ortogonal de tamaño n × n y definimos

v(x) := u (Ox) (x ∈ Rn ) ,

entonces ∆v = 0.

24
Saul Ancari Villca Ecuaciones diferenciales parciales

2. Encuentre todas las funciones armónicas u : R2 → R de la forma

2
X
u(x) = aij xi xj .
i,j=1

3. (a) Sean u y v armónicas sobre el abierto U ⊂ Rn . Muestre que uv es armónica sobre

U sí y sólo sí ∇u · ∇v = 0.
(b) Sea U ⊂ Rn abierto y conexo. Si u y u2 son armónicas sobre U , muestre que u es
la función constante.

4. Definamos
P = {x ∈ Rn | |xn | = 1}.
Encuentre una función armónica h 6≡ 0 de tal forma que h = 0 sobre P .
Sugerencia: Pruebe con las funciones del tipo eaxn−1 cos bxn , donde a y b son constan-
tes.

5. Sea U un subconjunto abierto de Rn . Decimos que una función v ∈ C 2 (Ū ) es sub

armónica si ∆v ≥ 0 en U .

(a) Pruebe que una función sub armónica v, satisface

v(x) ≤ v(y)dy para todo B(x, r) ⊂ U.

B(x,r)

(b) Si U es limitado, pruebe que

(i) máxŪ v = máx∂U v.
(ii) Además, si U es conexo y existe un punto x0 ∈ U tal que

v(x0 ) = máx v,
Ū

entonces v es constante en U .

6. Sea U un subconjunto abierto, limitado de Rn . Pruebe que existe una constante C,

dependiendo solo de U , tal que

máx u ≤ (máx |g| + máx |f |)

Ū ∂ Ū Ū

siempre que u sea una solución suave de


−∆u = f en U
u = g sobre ∂U.

|x|2
Sugerencia: −∆(u + 2n
λ) ≤ 0, para λ := máxŪ |f |.

7. Suponga que u es suave y es la solución de ut − ∆u = 0 en Rn × (0, ∞).

25
Saul Ancari Villca Ecuaciones diferenciales parciales

(a) Muestre que uλ (x, t) := u (λx, λ2 t) también resuelve la ecuación del calor para
cada λ ∈ R.
(b) Utilice (a) para mostrar que v(x, t) := hx, ∇u(x, t)i + 2tut (x, t) también resuelve
la ecuación del calor.
 
8. Asuma n = 1 y u(x, t) = v √xt .

(a) Muestre que

ut = uxx
si y solo si
z
v 00 + v 0 = 0. (20)
2
Demuestre que la solución general de (20) es
ˆ z
2
v(z) = c e−s /4 ds + d.
0
 
x
(b) Derive u(x, t) = v √con respecto a x y seleccione la constante c de forma
t
apropiada para obtener la solución fundamental Φ para n = 1. Explique por qué
este procedimiento produce la solución fundamental. (Sugerencia: ¿Cuál es la
condición inicial para u? )

9. Asuma n = 1. El propósito de este ejercicio es mostrar que el principio máximo no

es cierto para la ecuación ut − xuxx = 0, que tiene un coeficiente variable.

(a) Verifique que u = −2xt − x2 es una solución.

(b) Encuentre la ubicación del máximo de u en el rectángulo cerrado [−2, 2] × [0, 1].

10. Sea U un subconjunto abierto, limitado de Rn y T > 0. Suponga que u1 , u2 ∈

C12 (UT ) ∩ C(ŪT ) satisfacen

uit − ∆ui = f en UT
ui = gi sobre ΓT ,

para i = 1, 2, donde f ∈ C(UT ) y g1 , g2 ∈ C(ΓT ). Entonces |u1 − u2 | ≤ máx |g1 − g2 |.

ΓT

Referencias
[1] Lawrence C Evans. “Partial differential equations”. En: Graduate studies in mathematics
19.4 (1998).

26