Ecuaciones Diferenciales Parciales

Saul Ancari Villca

21 de septiembre de 2021

Saul Ancari Villca EDP 21 de septiembre de 2021 1 / 34

Introducción
Una ecuación diferencial parcial (EDP) es una ecuación que involucra una función
desconocida de dos o más variables y algunas de sus derivadas parciales.

Definición 1
Sea U un subconjunto abierto de Rn y k ≥ 1 un entero. Una expresión de la forma
 
F D k u(x ), D k −1 u(x ), . . . , Du(x ), u(x ), x = 0 (x ∈ U ) (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.

Saul Ancari Villca EDP 21 de septiembre de 2021 2 / 34

Introducción
Una ecuación diferencial parcial (EDP) es una ecuación que involucra una función
desconocida de dos o más variables y algunas de sus derivadas parciales.

Definición 1
Sea U un subconjunto abierto de Rn y k ≥ 1 un entero. Una expresión de la forma
 
F D k u(x ), D k −1 u(x ), . . . , Du(x ), u(x ), x = 0 (x ∈ U ) (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.

Saul Ancari Villca EDP 21 de septiembre de 2021 2 / 34

Clasificación
(i) Lineal si tiene la forma
aα (x )D α u = f (x )
X

|α|≤k

para funciones dadas aα (|α| ≤ k ), f . Si f ≡ 0 es llamada homogénea.

(ii) Semilineal si tiene la forma

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

|α|=k

(iii) Casi lineal si tiene la forma

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

|α|=k

(iv) La EDP (1) es totalmente no lineal si depende de forma no lineal de las derivadas de
orden superior.
Saul Ancari Villca EDP 21 de septiembre de 2021 3 / 34
Clasificación
(i) Lineal si tiene la forma
aα (x )D α u = f (x )
X

|α|≤k

para funciones dadas aα (|α| ≤ k ), f . Si f ≡ 0 es llamada homogénea.

(ii) Semilineal si tiene la forma

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

|α|=k

(iii) Casi lineal si tiene la forma

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

|α|=k

(iv) La EDP (1) es totalmente no lineal si depende de forma no lineal de las derivadas de
orden superior.
Saul Ancari Villca EDP 21 de septiembre de 2021 3 / 34
Clasificación
(i) Lineal si tiene la forma
aα (x )D α u = f (x )
X

|α|≤k

para funciones dadas aα (|α| ≤ k ), f . Si f ≡ 0 es llamada homogénea.

(ii) Semilineal si tiene la forma

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

|α|=k

(iii) Casi lineal si tiene la forma

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

|α|=k

(iv) La EDP (1) es totalmente no lineal si depende de forma no lineal de las derivadas de
orden superior.
Saul Ancari Villca EDP 21 de septiembre de 2021 3 / 34
Clasificación
(i) Lineal si tiene la forma
aα (x )D α u = f (x )
X

|α|≤k

para funciones dadas aα (|α| ≤ k ), f . Si f ≡ 0 es llamada homogénea.

(ii) Semilineal si tiene la forma

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

|α|=k

(iii) Casi lineal si tiene la forma

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

|α|=k

(iv) La EDP (1) es totalmente no lineal si depende de forma no lineal de las derivadas de
orden superior.
Saul Ancari Villca EDP 21 de septiembre de 2021 3 / 34
Ejemplos

1. Ecuación del Laplace

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

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

ut − ∆u = 0.

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

Saul Ancari Villca EDP 21 de septiembre de 2021 4 / 34

Ejemplos

1. Ecuación del Laplace

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

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

ut − ∆u = 0.

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

Saul Ancari Villca EDP 21 de septiembre de 2021 4 / 34

Ejemplos

1. Ecuación del Laplace

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

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

ut − ∆u = 0.

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

Saul Ancari Villca EDP 21 de septiembre de 2021 4 / 34

Ecuación de Laplace y Poisson

En esta sección estudiaremos la ecuación de Laplace

∆u = 0 (2)
y la ecuación de Poisson
− ∆u = f (3)
En ambos, (2) and (3), U ⊂ Rn es abierto y u : U → R es la incógnita. En (3) la función
f : U → R esta dada.
Definición 2
Una función u ∈ C 2 que satisface (2) sera llamada función armónica.

Saul Ancari Villca EDP 21 de septiembre de 2021 5 / 34

Ecuación de Laplace y Poisson

En esta sección estudiaremos la ecuación de Laplace

∆u = 0 (2)
y la ecuación de Poisson
− ∆u = f (3)
En ambos, (2) and (3), U ⊂ Rn es abierto y u : U → R es la incógnita. En (3) la función
f : U → R esta dada.
Definición 2
Una función u ∈ C 2 que satisface (2) sera llamada función armónica.

Saul Ancari Villca EDP 21 de septiembre de 2021 5 / 34

Solución fundamental
Intentemos encontrar una solución u de la ecuación de Laplace (2) en U = Rn , que tiene la
forma
u(x ) = v (r ),
 1/2
donde r = |x | = x12 + ... + xn2 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 + · · · + xn2 2 xi = (x , 0).
∂xi 2 r
Luego

xi
uxi = v 0 (r ) ,
r
xi2
xi2 xi2  1 xi2 
 
00 0
r− r 00 0
uxi xi = v (r ) + v (r ) = v (r ) + v (r )  − 3 
r2 r2 r2 r r

Saul Ancari Villca EDP 21 de septiembre de 2021 6 / 34

Solución fundamental
Intentemos encontrar una solución u de la ecuación de Laplace (2) en U = Rn , que tiene la
forma
u(x ) = v (r ),
 1/2
donde r = |x | = x12 + ... + xn2 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 + · · · + xn2 2 xi = (x , 0).
∂xi 2 r
Luego

xi
uxi = v 0 (r ) ,
r
xi2
xi2 xi2  1 xi2 
 
00 0
r− r 00 0
uxi xi = v (r ) + v (r ) = v (r ) + v (r )  − 3 
r2 r2 r2 r r

Saul Ancari Villca EDP 21 de septiembre de 2021 6 / 34

Solución fundamental
Intentemos encontrar una solución u de la ecuación de Laplace (2) en U = Rn , que tiene la
forma
u(x ) = v (r ),
 1/2
donde r = |x | = x12 + ... + xn2 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 + · · · + xn2 2 xi = (x , 0).
∂xi 2 r
Luego

xi
uxi = v 0 (r ) ,
r
xi2
xi2 xi2  1 xi2 
 
00 0
r− r 00 0
uxi xi = v (r ) + v (r ) = v (r ) + v (r )  − 3 
r2 r2 r2 r r

Saul Ancari Villca EDP 21 de septiembre de 2021 6 / 34

Solución fundamental
Intentemos encontrar una solución u de la ecuación de Laplace (2) en U = Rn , que tiene la
forma
u(x ) = v (r ),
 1/2
donde r = |x | = x12 + ... + xn2 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 + · · · + xn2 2 xi = (x , 0).
∂xi 2 r
Luego

xi
uxi = v 0 (r ) ,
r
xi2
xi2 xi2  1 xi2 
 
00 0
r− r 00 0
uxi xi = v (r ) + v (r ) = v (r ) + v (r )  − 3 
r2 r2 r2 r r

Saul Ancari Villca EDP 21 de septiembre de 2021 6 / 34

Solución fundamental
Intentemos encontrar una solución u de la ecuación de Laplace (2) en U = Rn , que tiene la
forma
u(x ) = v (r ),
 1/2
donde r = |x | = x12 + ... + xn2 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 + · · · + xn2 2 xi = (x , 0).
∂xi 2 r
Luego

xi
uxi = v 0 (r ) ,
r
xi2
xi2 xi2  1 xi2 
 
00 0
r− r 00 0
uxi xi = v (r ) + v (r ) = v (r ) + v (r )  − 3 
r2 r2 r2 r r

Saul Ancari Villca EDP 21 de septiembre de 2021 6 / 34

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
Si v 0 , 0, obtenemos
v 00 1−n
log (|v 0 |)0 = 0
= ,
v r
y por lo tanto v 0 (r ) = a
r n−1
para alguna constante a. De esta forma, si r > 0, tenemos

b log r + c

 (n = 2)
v (r ) = 
 nb−2 + c
 (n ≥ 3),
r

donde b y c son constantes.

Saul Ancari Villca EDP 21 de septiembre de 2021 7 / 34

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
Si v 0 , 0, obtenemos
v 00 1−n
log (|v 0 |)0 = 0
= ,
v r
y por lo tanto v 0 (r ) = a
r n−1
para alguna constante a. De esta forma, si r > 0, tenemos

b log r + c

 (n = 2)
v (r ) = 
 nb−2 + c
 (n ≥ 3),
r

donde b y c son constantes.

Saul Ancari Villca EDP 21 de septiembre de 2021 7 / 34

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
Si v 0 , 0, obtenemos
v 00 1−n
log (|v 0 |)0 = 0
= ,
v r
y por lo tanto v 0 (r ) = a
r n−1
para alguna constante a. De esta forma, si r > 0, tenemos

b log r + c

 (n = 2)
v (r ) = 
 nb−2 + c
 (n ≥ 3),
r

donde b y c son constantes.

Saul Ancari Villca EDP 21 de septiembre de 2021 7 / 34

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
Si v 0 , 0, obtenemos
v 00 1−n
log (|v 0 |)0 = 0
= ,
v r
y por lo tanto v 0 (r ) = a
r n−1
para alguna constante a. De esta forma, si r > 0, tenemos

b log r + c

 (n = 2)
v (r ) = 
 nb−2 + c
 (n ≥ 3),
r

donde b y c son constantes.

Saul Ancari Villca EDP 21 de septiembre de 2021 7 / 34

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
Si v 0 , 0, obtenemos
v 00 1−n
log (|v 0 |)0 = 0
= ,
v r
y por lo tanto v 0 (r ) = a
r n−1
para alguna constante a. De esta forma, si r > 0, tenemos

b log r + c

 (n = 2)
v (r ) = 
 nb−2 + c
 (n ≥ 3),
r

donde b y c son constantes.

Saul Ancari Villca EDP 21 de septiembre de 2021 7 / 34

Definición 3
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).

Saul Ancari Villca EDP 21 de septiembre de 2021 8 / 34

Ahora usando la solución fundamental de la ecuación de Laplace, construiremos una solución
de la ecuación de Poisson.
Teorema 1 (Solución de la ecuación de Poisson)
Sea f ∈ Cc2 (Rn ) y Z
u(x ) = Φ(x − y )f (y )dy
Rn
 1 R (4)
− 2π R2 log(|

 x − y |)f (y )dy (n = 2)
= 1
R f (y )

 n n−2 dy
n(n−2)α(n) R |x −y |
(n ≥ 3)
Entonces
(i) u ∈ C 2 (Rn )
(ii) −∆u = f en Rn .

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

Saul Ancari Villca EDP 21 de septiembre de 2021 9 / 34

Ahora usando la solución fundamental de la ecuación de Laplace, construiremos una solución
de la ecuación de Poisson.
Teorema 1 (Solución de la ecuación de Poisson)
Sea f ∈ Cc2 (Rn ) y Z
u(x ) = Φ(x − y )f (y )dy
Rn
 1 R (4)
− 2π R2 log(|

 x − y |)f (y )dy (n = 2)
= 1
R f (y )

 n n−2 dy
n(n−2)α(n) R |x −y |
(n ≥ 3)
Entonces
(i) u ∈ C 2 (Rn )
(ii) −∆u = f en Rn .

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

Saul Ancari Villca EDP 21 de septiembre de 2021 9 / 34

Demostración del Teorema 1

1. Tenemos Z Z
u (x ) = Φ(x − y )f (y )dy = Φ(y )f (x − y )dy ;
Rn Rn
luego Z " #
u (x + hei ) − u(x ) f (x + hei − y ) − f (x − y )
= Φ(y ) dy ,
h Rn h
donde h , 0 y ei = (0, . . . , 1, . . . , 0), el 1 esta en el i −ésimo lugar.Ya que el soporte de f es
compacto , existe r suficientemente grande tal que spt(f (x + hei − y ) − f (x − y )) ⊂ B (0, r )
para todo |h | ≤ 1. Además por el teorema del valor medio tenemos que

 f (x + hei − y ) − f (x − y )   ∂f
   
 =  (x + θhei − y ) ≤ M ,

h ∂xi


para algún θ ∈ (0,1), donde M = máx | ∂∂xfi |.

Saul Ancari Villca EDP 21 de septiembre de 2021 10 / 34

Demostración del Teorema 1

1. Tenemos Z Z
u (x ) = Φ(x − y )f (y )dy = Φ(y )f (x − y )dy ;
Rn Rn
luego Z " #
u (x + hei ) − u(x ) f (x + hei − y ) − f (x − y )
= Φ(y ) dy ,
h Rn h
donde h , 0 y ei = (0, . . . , 1, . . . , 0), el 1 esta en el i −ésimo lugar.Ya que el soporte de f es
compacto , existe r suficientemente grande tal que spt(f (x + hei − y ) − f (x − y )) ⊂ B (0, r )
para todo |h | ≤ 1. Además por el teorema del valor medio tenemos que

 f (x + hei − y ) − f (x − y )   ∂f
   
 =  (x + θhei − y ) ≤ M ,

h ∂xi


para algún θ ∈ (0,1), donde M = máx | ∂∂xfi |.

Saul Ancari Villca EDP 21 de septiembre de 2021 10 / 34

Demostración del Teorema 1

1. Tenemos Z Z
u (x ) = Φ(x − y )f (y )dy = Φ(y )f (x − y )dy ;
Rn Rn
luego Z " #
u (x + hei ) − u(x ) f (x + hei − y ) − f (x − y )
= Φ(y ) dy ,
h Rn h
donde h , 0 y ei = (0, . . . , 1, . . . , 0), el 1 esta en el i −ésimo lugar.Ya que el soporte de f es
compacto , existe r suficientemente grande tal que spt(f (x + hei − y ) − f (x − y )) ⊂ B (0, r )
para todo |h | ≤ 1. Además por el teorema del valor medio tenemos que

 f (x + hei − y ) − f (x − y )   ∂f
   
 =  (x + θhei − y ) ≤ M ,

h ∂xi


para algún θ ∈ (0,1), donde M = máx | ∂∂xfi |.

Saul Ancari Villca EDP 21 de septiembre de 2021 10 / 34

Demostración del Teorema 1

1. Tenemos Z Z
u (x ) = Φ(x − y )f (y )dy = Φ(y )f (x − y )dy ;
Rn Rn
luego Z " #
u (x + hei ) − u(x ) f (x + hei − y ) − f (x − y )
= Φ(y ) dy ,
h Rn h
donde h , 0 y ei = (0, . . . , 1, . . . , 0), el 1 esta en el i −ésimo lugar.Ya que el soporte de f es
compacto , existe r suficientemente grande tal que spt(f (x + hei − y ) − f (x − y )) ⊂ B (0, r )
para todo |h | ≤ 1. Además por el teorema del valor medio tenemos que

 f (x + hei − y ) − f (x − y )   ∂f
   
 =  (x + θhei − y ) ≤ M ,

h ∂xi


para algún θ ∈ (0,1), donde M = máx | ∂∂xfi |.

Saul Ancari Villca EDP 21 de septiembre de 2021 10 / 34

Por otro lado, como f ∈ C 1 (Rn ) tenemos que

f (x + hei − y ) − f (x − y ) ∂f
→ (x − y )
h ∂xi
sobre Rn cuando h → 0. Aplicando el teorema de convergencia dominada, obtenemos

∂u ∂f
Z
(x ) = Φ(y ) (x − y )dy (i = 1, . . . , n).
∂ xi Rn ∂ xi
Similarmente
∂2 u ∂2 f
Z
(x ) = Φ(y ) (x − y )dy (i , j = 1, . . . , n). (5)
∂xi ∂xj Rn ∂xi ∂xj
Como la expresión del lado derecho de (5) es continua en la variable x, vemos que
u ∈ C 2 (Rn ).

Saul Ancari Villca EDP 21 de septiembre de 2021 11 / 34

Por otro lado, como f ∈ C 1 (Rn ) tenemos que

f (x + hei − y ) − f (x − y ) ∂f
→ (x − y )
h ∂xi
sobre Rn cuando h → 0. Aplicando el teorema de convergencia dominada, obtenemos

∂u ∂f
Z
(x ) = Φ(y ) (x − y )dy (i = 1, . . . , n).
∂ xi Rn ∂ xi
Similarmente
∂2 u ∂2 f
Z
(x ) = Φ(y ) (x − y )dy (i , j = 1, . . . , n). (5)
∂xi ∂xj Rn ∂xi ∂xj
Como la expresión del lado derecho de (5) es continua en la variable x, vemos que
u ∈ C 2 (Rn ).

Saul Ancari Villca EDP 21 de septiembre de 2021 11 / 34

Por otro lado, como f ∈ C 1 (Rn ) tenemos que

f (x + hei − y ) − f (x − y ) ∂f
→ (x − y )
h ∂xi
sobre Rn cuando h → 0. Aplicando el teorema de convergencia dominada, obtenemos

∂u ∂f
Z
(x ) = Φ(y ) (x − y )dy (i = 1, . . . , n).
∂ xi Rn ∂ xi
Similarmente
∂2 u ∂2 f
Z
(x ) = Φ(y ) (x − y )dy (i , j = 1, . . . , n). (5)
∂xi ∂xj Rn ∂xi ∂xj
Como la expresión del lado derecho de (5) es continua en la variable x, vemos que
u ∈ C 2 (Rn ).

Saul Ancari Villca EDP 21 de septiembre de 2021 11 / 34

2. Dado que Φ(x ) → ∞ cuando x → 0, necesitaremos cálculos posteriores para aislar esta
singularidad dentro de una bola pequeña. Fijemos  > 0, entonces
Z Z
∆u(x ) = Φ(y )∆x f (x − y )dy + Φ(y )∆x f (x − y )dy
por (5) B (0,) Rn −B (0,) (6)
=: I + J

Luego 
C  2 | log | (n = 2)
2 Z 

|I | ≤ C D (f ) L ∞ (Rn ) |Φ(y )|dy ≤  (7)
B (0,) C  2
 (n ≥ 3).

Saul Ancari Villca EDP 21 de septiembre de 2021 12 / 34

2. Dado que Φ(x ) → ∞ cuando x → 0, necesitaremos cálculos posteriores para aislar esta
singularidad dentro de una bola pequeña. Fijemos  > 0, entonces
Z Z
∆u(x ) = Φ(y )∆x f (x − y )dy + Φ(y )∆x f (x − y )dy
por (5) B (0,) Rn −B (0,) (6)
=: I + J

Luego 
C  2 | log | (n = 2)
2 Z 

|I | ≤ C D (f ) L ∞ (Rn ) |Φ(y )|dy ≤  (7)
B (0,) C  2
 (n ≥ 3).

Saul Ancari Villca EDP 21 de septiembre de 2021 12 / 34

2. Dado que Φ(x ) → ∞ cuando x → 0, necesitaremos cálculos posteriores para aislar esta
singularidad dentro de una bola pequeña. Fijemos  > 0, entonces
Z Z
∆u(x ) = Φ(y )∆x f (x − y )dy + Φ(y )∆x f (x − y )dy
por (5) B (0,) Rn −B (0,) (6)
=: I + J

Luego 
C  2 | log | (n = 2)
2 Z 

|I | ≤ C D (f ) L ∞ (Rn ) |Φ(y )|dy ≤  (7)
B (0,) C  2
 (n ≥ 3).

Saul Ancari Villca EDP 21 de septiembre de 2021 12 / 34

Integrando por partes obtenemos
Z
J = Φ(y )∆y f (x − y )dy
Rn −B (0,)
Z D E
=− ∇Φ(y ), ∇y f (x − y ) dy
Rn −B (0,) (8)
∂f
Z
+ Φ(y ) (x − y )dS (y )
∂B (0,) ∂ν
=: K + L ,

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

Saul Ancari Villca EDP 21 de septiembre de 2021 13 / 34

Integrando por partes obtenemos
Z
J = Φ(y )∆y f (x − y )dy
Rn −B (0,)
Z D E
=− ∇Φ(y ), ∇y f (x − y ) dy
Rn −B (0,) (8)
∂f
Z
+ Φ(y ) (x − y )dS (y )
∂B (0,) ∂ν
=: K + L ,

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

Saul Ancari Villca EDP 21 de septiembre de 2021 13 / 34

3. Integramos por partes el término K , para obtener
∂Φ
Z Z
K = ∆Φ(y )f (x − y )dy − (y )f (x − y )dS (y )
Rn −B (0,) ∂B (0,) ∂ν
∂Φ
Z
=− (y )f (x − y )dS (y )
∂B (0,) ∂ν
−1 y −y
ya que Φ es armónica fuera del origen. Ahora ∇Φ(y ) = nα(n) |y |n
(y , 0) y ν = |y | = − y sobre
∂Φ
∂B (0, ). Como consecuencia ∂ν (y ) = hν, ∇Φ(y )i = 1
nα(n) n−1
sobre ∂B (0, ). Puesto que
nα(n)ε n−1
es el área de la superficie de la esfera ∂B (0, ), tenemos
Z
1
K = − f (x − y )dS (y )
nα(n) n−1 ∂B (0,)
? (10)
=− f (y )dS (y ) → −f (x ) cuando →0
∂B (x ,)

4. Combinando (6)-(10) y haciendo  → 0, concluimos que −∆u(x ) = f (x ), como habı́amos

afirmado

Saul Ancari Villca EDP 21 de septiembre de 2021 14 / 34
3. Integramos por partes el término K , para obtener
∂Φ
Z Z
K = ∆Φ(y )f (x − y )dy − (y )f (x − y )dS (y )
Rn −B (0,) ∂B (0,) ∂ν
∂Φ
Z
=− (y )f (x − y )dS (y )
∂B (0,) ∂ν
−1 y −y
ya que Φ es armónica fuera del origen. Ahora ∇Φ(y ) = nα(n) |y |n
(y , 0) y ν = |y | = − y sobre
∂Φ
∂B (0, ). Como consecuencia ∂ν (y ) = hν, ∇Φ(y )i = 1
nα(n) n−1
sobre ∂B (0, ). Puesto que
nα(n)ε n−1
es el área de la superficie de la esfera ∂B (0, ), tenemos
Z
1
K = − f (x − y )dS (y )
nα(n) n−1 ∂B (0,)
? (10)
=− f (y )dS (y ) → −f (x ) cuando →0
∂B (x ,)

4. Combinando (6)-(10) y haciendo  → 0, concluimos que −∆u(x ) = f (x ), como habı́amos

afirmado

Saul Ancari Villca EDP 21 de septiembre de 2021 14 / 34
3. Integramos por partes el término K , para obtener
∂Φ
Z Z
K = ∆Φ(y )f (x − y )dy − (y )f (x − y )dS (y )
Rn −B (0,) ∂B (0,) ∂ν
∂Φ
Z
=− (y )f (x − y )dS (y )
∂B (0,) ∂ν
−1 y −y
ya que Φ es armónica fuera del origen. Ahora ∇Φ(y ) = nα(n) |y |n
(y , 0) y ν = |y | = − y sobre
∂Φ
∂B (0, ). Como consecuencia ∂ν (y ) = hν, ∇Φ(y )i = 1
nα(n) n−1
sobre ∂B (0, ). Puesto que
nα(n)ε n−1
es el área de la superficie de la esfera ∂B (0, ), tenemos
Z
1
K = − f (x − y )dS (y )
nα(n) n−1 ∂B (0,)
? (10)
=− f (y )dS (y ) → −f (x ) cuando →0
∂B (x ,)

4. Combinando (6)-(10) y haciendo  → 0, concluimos que −∆u(x ) = f (x ), como habı́amos

afirmado

Saul Ancari Villca EDP 21 de septiembre de 2021 14 / 34
3. Integramos por partes el término K , para obtener
∂Φ
Z Z
K = ∆Φ(y )f (x − y )dy − (y )f (x − y )dS (y )
Rn −B (0,) ∂B (0,) ∂ν
∂Φ
Z
=− (y )f (x − y )dS (y )
∂B (0,) ∂ν
−1 y −y
ya que Φ es armónica fuera del origen. Ahora ∇Φ(y ) = nα(n) |y |n
(y , 0) y ν = |y | = − y sobre
∂Φ
∂B (0, ). Como consecuencia ∂ν (y ) = hν, ∇Φ(y )i = 1
nα(n) n−1
sobre ∂B (0, ). Puesto que
nα(n)ε n−1
es el área de la superficie de la esfera ∂B (0, ), tenemos
Z
1
K = − f (x − y )dS (y )
nα(n) n−1 ∂B (0,)
? (10)
=− f (y )dS (y ) → −f (x ) cuando →0
∂B (x ,)

4. Combinando (6)-(10) y haciendo  → 0, concluimos que −∆u(x ) = f (x ), como habı́amos

afirmado

Saul Ancari Villca EDP 21 de septiembre de 2021 14 / 34
3. Integramos por partes el término K , para obtener
∂Φ
Z Z
K = ∆Φ(y )f (x − y )dy − (y )f (x − y )dS (y )
Rn −B (0,) ∂B (0,) ∂ν
∂Φ
Z
=− (y )f (x − y )dS (y )
∂B (0,) ∂ν
−1 y −y
ya que Φ es armónica fuera del origen. Ahora ∇Φ(y ) = nα(n) |y |n
(y , 0) y ν = |y | = − y sobre
∂Φ
∂B (0, ). Como consecuencia ∂ν (y ) = hν, ∇Φ(y )i = 1
nα(n) n−1
sobre ∂B (0, ). Puesto que
nα(n)ε n−1
es el área de la superficie de la esfera ∂B (0, ), tenemos
Z
1
K = − f (x − y )dS (y )
nα(n) n−1 ∂B (0,)
? (10)
=− f (y )dS (y ) → −f (x ) cuando →0
∂B (x ,)

4. Combinando (6)-(10) y haciendo  → 0, concluimos que −∆u(x ) = f (x ), como habı́amos

afirmado

Saul Ancari Villca EDP 21 de septiembre de 2021 14 / 34
Formula del valor medio
Consideraremos un conjunto abierto U ⊂ Rn . En lo que sigue usaremos las siguientes
notaciones:
? Z
1
f (y )dy := f (y )dy
B (x ,r ) α(n)r n B (x ,r )

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

? Z
1
f (y )dS (y ) := fdS
∂B (x ,r ) nα(n)r n−1 ∂B (x ,r )
= promedio de f sobre la esfera ∂B (x , r ),
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 .
Saul Ancari Villca EDP 21 de septiembre de 2021 15 / 34
Teorema 2 (Formula del valor medio)
Sea U un subconjunto abierto de Rn . Si u ∈ C 2 (U ) es armónica, entonces
? ?
u(x ) = u(y )dS (y ) = u(y )dy (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 ?
φ (r ) =
0
h∇u(x + rz ), z idS (z ).
∂B (0,1)

Saul Ancari Villca EDP 21 de septiembre de 2021 16 / 34

Teorema 2 (Formula del valor medio)
Sea U un subconjunto abierto de Rn . Si u ∈ C 2 (U ) es armónica, entonces
? ?
u(x ) = u(y )dS (y ) = u(y )dy (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 ?
φ (r ) =
0
h∇u(x + rz ), z idS (z ).
∂B (0,1)

Saul Ancari Villca EDP 21 de septiembre de 2021 16 / 34

Teorema 2 (Formula del valor medio)
Sea U un subconjunto abierto de Rn . Si u ∈ C 2 (U ) es armónica, entonces
? ?
u(x ) = u(y )dS (y ) = u(y )dy (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 ?
φ (r ) =
0
h∇u(x + rz ), z idS (z ).
∂B (0,1)

Saul Ancari Villca EDP 21 de septiembre de 2021 16 / 34

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

Entonces
?
r
φ (r ) =
0
∆u(y )dy (12)
n B (x ,r )
= 0 (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 )
Saul Ancari Villca EDP 21 de septiembre de 2021 17 / 34
Luego, ?
y −x
 
φ0 (r ) = ∇u(y ), dS (y )
∂B (x ,r ) r
∂u
Z
1
= (y )dS (y )
y −x nα(n)r n−1 ∂B (x ,r ) ∂ν
ν=
r
Z
1
= ∆u(y )dy .
formulas de Green nα(n)r n−1 B (x ,r )

Entonces
?
r
φ (r ) =
0
∆u(y )dy (12)
n B (x ,r )
= 0 (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 )
Saul Ancari Villca EDP 21 de septiembre de 2021 17 / 34
Luego, ?
y −x
 
φ0 (r ) = ∇u(y ), dS (y )
∂B (x ,r ) r
∂u
Z
1
= (y )dS (y )
y −x nα(n)r n−1 ∂B (x ,r ) ∂ν
ν=
r
Z
1
= ∆u(y )dy .
formulas de Green nα(n)r n−1 B (x ,r )

Entonces
?
r
φ (r ) =
0
∆u(y )dy (12)
n B (x ,r )
= 0 (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 )
Saul Ancari Villca EDP 21 de septiembre de 2021 17 / 34
Por tanto ?
u(y )dS (y ) = u(x ). (14)
∂B (x ,r )

2. Por otro lado, usando coordenadas polares tenemos

Z Z r Z !
u(y )dy = u(y )dS (y ) ds
B (x ,r ) 0 ∂B (x ,s )
Z r
= u (x ) nα(n)s n−1 ds = α(n)r n u(x )
por (14) 0

Saul Ancari Villca EDP 21 de septiembre de 2021 18 / 34

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

2. Por otro lado, usando coordenadas polares tenemos

Z Z r Z !
u(y )dy = u(y )dS (y ) ds
B (x ,r ) 0 ∂B (x ,s )
Z r
= u (x ) nα(n)s n−1 ds = α(n)r n u(x )
por (14) 0

Saul Ancari Villca EDP 21 de septiembre de 2021 18 / 34

Teorema 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 . 0 tenemos que existe x0 ∈ U y r0 > 0

tal que B (x0 , r0 )) ⊂ U y
∆u > 0 sobre B (x0 , r0 ). (15)
>
Considerando que φ(r ) = ∂B (x0 ,r )
u(y )dS (y ) = constante, tenemos
?
r
0 = φ (r )
0
= ∆u(y )dy > 0.
por (12) y (15) n B (x ,r )

Lo que es una contradicción.


Saul Ancari Villca EDP 21 de septiembre de 2021 19 / 34
Propiedades de las funciones armónicas
A continuación, presentamos una secuencia de resultados interesantes sobre funciones
armónicas, todas estas basadas en la formula del valor medio. En lo que sigue, U ⊂ Rn es
abierto y limitado.
Teorema 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 .
Saul Ancari Villca EDP 21 de septiembre de 2021 20 / 34
Propiedades de las funciones armónicas
A continuación, presentamos una secuencia de resultados interesantes sobre funciones
armónicas, todas estas basadas en la formula del valor medio. En lo que sigue, U ⊂ Rn es
abierto y limitado.
Teorema 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 .
Saul Ancari Villca EDP 21 de septiembre de 2021 20 / 34
Propiedades de las funciones armónicas
A continuación, presentamos una secuencia de resultados interesantes sobre funciones
armónicas, todas estas basadas en la formula del valor medio. En lo que sigue, U ⊂ Rn es
abierto y limitado.
Teorema 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 .
Saul Ancari Villca EDP 21 de septiembre de 2021 20 / 34
Propiedades de las funciones armónicas
A continuación, presentamos una secuencia de resultados interesantes sobre funciones
armónicas, todas estas basadas en la formula del valor medio. En lo que sigue, U ⊂ Rn es
abierto y limitado.
Teorema 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 .
Saul Ancari Villca EDP 21 de septiembre de 2021 20 / 34
Demostración del Teorema 4
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 ).

Saul Ancari Villca EDP 21 de septiembre de 2021 21 / 34

Demostración del Teorema 4
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 ).

Saul Ancari Villca EDP 21 de septiembre de 2021 21 / 34

Demostración del Teorema 4
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 ).

Saul Ancari Villca EDP 21 de septiembre de 2021 21 / 34

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 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


 u = g sobre ∂U .
 (16)


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

Demostración. Si u y e u), j = 0, 1,
son soluciones de 
−∆wj = 0 en U


sobre ∂U .


 wj = 0
Luego, aplicando el Teorema 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 = 0.
u −e
Saul Ancari Villca EDP 21 de septiembre de 2021 22 / 34
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 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


 u = g sobre ∂U .
 (16)


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

Demostración. Si u y e u), j = 0, 1,
son soluciones de 
−∆wj = 0 en U


sobre ∂U .


 wj = 0
Luego, aplicando el Teorema 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 = 0.
u −e
Saul Ancari Villca EDP 21 de septiembre de 2021 22 / 34
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 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


 u = g sobre ∂U .
 (16)


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

Demostración. Si u y e u), j = 0, 1,
son soluciones de 
−∆wj = 0 en U


sobre ∂U .


 wj = 0
Luego, aplicando el Teorema 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 = 0.
u −e
Saul Ancari Villca EDP 21 de septiembre de 2021 22 / 34
Regularidad
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. Para
 > 0 escribamos
U := {x ∈ U | dist(x , ∂U ) > }.

Definición 4
(i) Definamos η ∈ C ∞ (Rn ) por
  
1
C exp |x |2 −1
 si |x | < 1
η(x ) := 

0
 si |x | ≥ 1,

elegimos la constante C > 0 de tal forma que ηdx = 1.

R
Rn
(ii) Para cada  > 0, sea
1 x 
η (x ) := η
n 
Saul Ancari Villca EDP 21 de septiembre de 2021 23 / 34
Regularidad
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. Para
 > 0 escribamos
U := {x ∈ U | dist(x , ∂U ) > }.

Definición 4
(i) Definamos η ∈ C ∞ (Rn ) por
  
1
C exp |x |2 −1
 si |x | < 1
η(x ) := 

0
 si |x | ≥ 1,

elegimos la constante C > 0 de tal forma que ηdx = 1.

R
Rn
(ii) Para cada  > 0, sea
1 x 
η (x ) := η
n 
Saul Ancari Villca EDP 21 de septiembre de 2021 23 / 34
Regularidad
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. Para
 > 0 escribamos
U := {x ∈ U | dist(x , ∂U ) > }.

Definición 4
(i) Definamos η ∈ C ∞ (Rn ) por
  
1
C exp |x |2 −1
 si |x | < 1
η(x ) := 

0
 si |x | ≥ 1,

elegimos la constante C > 0 de tal forma que ηdx = 1.

R
Rn
(ii) Para cada  > 0, sea
1 x 
η (x ) := η
n 
Saul Ancari Villca EDP 21 de septiembre de 2021 23 / 34
Llamamos a η el regularizador estándar. La función η es suave y satisface
Z
η dx = 1, spt (η ) = B (0, ).
Rn

Definición 5
Si f : U → R es localmente integrable, definimos su regularizante
Z Z
f  (x ) := η (x − y )f (y )dy = η (y )f (x − y )dy
U B (0,)

para x ∈ U

Saul Ancari Villca EDP 21 de septiembre de 2021 24 / 34

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

Definición 5
Si f : U → R es localmente integrable, definimos su regularizante
Z Z
f  (x ) := η (x − y )f (y )dy = η (y )f (x − y )dy
U B (0,)

para x ∈ U

Saul Ancari Villca EDP 21 de septiembre de 2021 24 / 34

Teorema 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.
Si 1 ≤ p < ∞ y f ∈ Lloc (U ), entonces f  → f en Lloc (U ).
(iv)
p p

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

u ∈ C ∞ (U ).

Saul Ancari Villca EDP 21 de septiembre de 2021 25 / 34

Demostración del Teorema de regularidad
Sea u el regularizante de u. Probaremos que u es suave al demostrar que u ≡ uε sobre Uε .
En efecto, si x ∈ Uε , entonces
Z
ε
u (x ) = ηε (x − y )u(y )dy
U
Z !
1 |x − y |
= n η u(y )dy
ε B (x ,ε) ε
Z ε   Z !
1 r
= n η udS dr
ε 0 ε ∂B (x ,r )
Z ε  
1 r
= u(x ) η nα(n)r n−1 dr
por (11) ε n
0 ε
Z
= u(x ) ηε dy = u(x ).
B (0,ε)

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

Saul Ancari Villca EDP 21 de septiembre de 2021 
26 / 34
Demostración del Teorema de regularidad
Sea u el regularizante de u. Probaremos que u es suave al demostrar que u ≡ uε sobre Uε .
En efecto, si x ∈ Uε , entonces
Z
ε
u (x ) = ηε (x − y )u(y )dy
U
Z !
1 |x − y |
= n η u(y )dy
ε B (x ,ε) ε
Z ε   Z !
1 r
= n η udS dr
ε 0 ε ∂B (x ,r )
Z ε  
1 r
= u(x ) η nα(n)r n−1 dr
por (11) ε n
0 ε
Z
= u(x ) ηε dy = u(x ).
B (0,ε)

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

Saul Ancari Villca EDP 21 de septiembre de 2021 
26 / 34
Demostración del Teorema de regularidad
Sea u el regularizante de u. Probaremos que u es suave al demostrar que u ≡ uε sobre Uε .
En efecto, si x ∈ Uε , entonces
Z
ε
u (x ) = ηε (x − y )u(y )dy
U
Z !
1 |x − y |
= n η u(y )dy
ε B (x ,ε) ε
Z ε   Z !
1 r
= n η udS dr
ε 0 ε ∂B (x ,r )
Z ε  
1 r
= u(x ) η nα(n)r n−1 dr
por (11) ε n
0 ε
Z
= u(x ) ηε dy = u(x ).
B (0,ε)

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

Saul Ancari Villca EDP 21 de septiembre de 2021 
26 / 34
Demostración del Teorema de regularidad
Sea u el regularizante de u. Probaremos que u es suave al demostrar que u ≡ uε sobre Uε .
En efecto, si x ∈ Uε , entonces
Z
ε
u (x ) = ηε (x − y )u(y )dy
U
Z !
1 |x − y |
= n η u(y )dy
ε B (x ,ε) ε
Z ε   Z !
1 r
= n η udS dr
ε 0 ε ∂B (x ,r )
Z ε  
1 r
= u(x ) η nα(n)r n−1 dr
por (11) ε n
0 ε
Z
= u(x ) ηε dy = u(x ).
B (0,ε)

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

Saul Ancari Villca EDP 21 de septiembre de 2021 
26 / 34
Demostración del Teorema de regularidad
Sea u el regularizante de u. Probaremos que u es suave al demostrar que u ≡ uε sobre Uε .
En efecto, si x ∈ Uε , entonces
Z
ε
u (x ) = ηε (x − y )u(y )dy
U
Z !
1 |x − y |
= n η u(y )dy
ε B (x ,ε) ε
Z ε   Z !
1 r
= n η udS dr
ε 0 ε ∂B (x ,r )
Z ε  
1 r
= u(x ) η nα(n)r n−1 dr
por (11) ε n
0 ε
Z
= u(x ) ηε dy = u(x ).
B (0,ε)

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

Saul Ancari Villca EDP 21 de septiembre de 2021 
26 / 34
Demostración del Teorema de regularidad
Sea u el regularizante de u. Probaremos que u es suave al demostrar que u ≡ uε sobre Uε .
En efecto, si x ∈ Uε , entonces
Z
ε
u (x ) = ηε (x − y )u(y )dy
U
Z !
1 |x − y |
= n η u(y )dy
ε B (x ,ε) ε
Z ε   Z !
1 r
= n η udS dr
ε 0 ε ∂B (x ,r )
Z ε  
1 r
= u(x ) η nα(n)r n−1 dr
por (11) ε n
0 ε
Z
= u(x ) ηε dy = u(x ).
B (0,ε)

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

Saul Ancari Villca EDP 21 de septiembre de 2021 
26 / 34
Demostración del Teorema de regularidad
Sea u el regularizante de u. Probaremos que u es suave al demostrar que u ≡ uε sobre Uε .
En efecto, si x ∈ Uε , entonces
Z
ε
u (x ) = ηε (x − y )u(y )dy
U
Z !
1 |x − y |
= n η u(y )dy
ε B (x ,ε) ε
Z ε   Z !
1 r
= n η udS dr
ε 0 ε ∂B (x ,r )
Z ε  
1 r
= u(x ) η nα(n)r n−1 dr
por (11) ε n
0 ε
Z
= u(x ) ηε dy = u(x ).
B (0,ε)

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

Saul Ancari Villca EDP 21 de septiembre de 2021 
26 / 34
Teorema de Liouville
Aquı́ probaremos que no existen funciones armónicas no triviales, limitadas y definidas en
todo Rn .
Teorema 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 formula del valor medio a uxi con r > 0,
tenemos
?
2n
Z
|uxi (x0 )| = | uxi dx | = | uνi dS |
B (x0 , 2r ) α(n)r n ∂B (x0 , 2r )
2n 2n
≤ kukL ∞ (∂B (x0 , 2r )) ≤ M.
r r
Haciendo r → ∞, concluimos que uxi (x0 ) = 0 para todo x0 ∈ Rn , i = 1, 2, ...n. Esto prueba que
u es constante.
Saul Ancari Villca EDP 21 de septiembre de 2021 27 / 34
Teorema de Liouville
Aquı́ probaremos que no existen funciones armónicas no triviales, limitadas y definidas en
todo Rn .
Teorema 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 formula del valor medio a uxi con r > 0,
tenemos
?
2n
Z
|uxi (x0 )| = | uxi dx | = | uνi dS |
B (x0 , 2r ) α(n)r n ∂B (x0 , 2r )
2n 2n
≤ kukL ∞ (∂B (x0 , 2r )) ≤ M.
r r
Haciendo r → ∞, concluimos que uxi (x0 ) = 0 para todo x0 ∈ Rn , i = 1, 2, ...n. Esto prueba que
u es constante.
Saul Ancari Villca EDP 21 de septiembre de 2021 27 / 34
Teorema de Liouville
Aquı́ probaremos que no existen funciones armónicas no triviales, limitadas y definidas en
todo Rn .
Teorema 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 formula del valor medio a uxi con r > 0,
tenemos
?
2n
Z
|uxi (x0 )| = | uxi dx | = | uνi dS |
B (x0 , 2r ) α(n)r n ∂B (x0 , 2r )
2n 2n
≤ kukL ∞ (∂B (x0 , 2r )) ≤ M.
r r
Haciendo r → ∞, concluimos que uxi (x0 ) = 0 para todo x0 ∈ Rn , i = 1, 2, ...n. Esto prueba que
u es constante.
Saul Ancari Villca EDP 21 de septiembre de 2021 27 / 34
Teorema de Liouville
Aquı́ probaremos que no existen funciones armónicas no triviales, limitadas y definidas en
todo Rn .
Teorema 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 formula del valor medio a uxi con r > 0,
tenemos
?
2n
Z
|uxi (x0 )| = | uxi dx | = | uνi dS |
B (x0 , 2r ) α(n)r n ∂B (x0 , 2r )
2n 2n
≤ kukL ∞ (∂B (x0 , 2r )) ≤ M.
r r
Haciendo r → ∞, concluimos que uxi (x0 ) = 0 para todo x0 ∈ Rn , i = 1, 2, ...n. Esto prueba que
u es constante.
Saul Ancari Villca EDP 21 de septiembre de 2021 27 / 34
Teorema de Liouville
Aquı́ probaremos que no existen funciones armónicas no triviales, limitadas y definidas en
todo Rn .
Teorema 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 formula del valor medio a uxi con r > 0,
tenemos
?
2n
Z
|uxi (x0 )| = | uxi dx | = | uνi dS |
B (x0 , 2r ) α(n)r n ∂B (x0 , 2r )
2n 2n
≤ kukL ∞ (∂B (x0 , 2r )) ≤ M.
r r
Haciendo r → ∞, concluimos que uxi (x0 ) = 0 para todo x0 ∈ Rn , i = 1, 2, ...n. Esto prueba que
u es constante.
Saul Ancari Villca EDP 21 de septiembre de 2021 27 / 34
Teorema 9 (Formula de representación)
Sea f ∈ Cc2 (Rn ), n ≥ 3. Entonces cualquier solución limitada de

−∆u = f en Rn

tiene la forma Z
u(x ) = Φ(x − y )f (y )dy + C
Rn
para alguna constante C.
R
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.

Saul Ancari Villca EDP 21 de septiembre de 2021 28 / 34

Teorema 9 (Formula de representación)
Sea f ∈ Cc2 (Rn ), n ≥ 3. Entonces cualquier solución limitada de

−∆u = f en Rn

tiene la forma Z
u(x ) = Φ(x − y )f (y )dy + C
Rn
para alguna constante C.
R
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.

Saul Ancari Villca EDP 21 de septiembre de 2021 28 / 34

Función de Green
Sea U un subconjunto abierto, limitado de Rn con ∂U de clase C 1 . El propósito de esta parte
es presentar una formula 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
Z
[u(y )∆Φ(y − x ) − Φ(y − x )∆u(y )]dy
Vε
Z
= [u(y )h∇y Φ(y − x ), νi − Φ(y − x )h∇u, νi]dS (y ),
∂Vε

donde ν es el vector unitario exterior sobre ∂Vε .

Saul Ancari Villca EDP 21 de septiembre de 2021 29 / 34
Recordemos que ∆Φ(x − y ) = 0 para x , y. Entonces
Z Z
− Φ(y − x )∆u(y )dy = [u(y )h∇y Φ(y − x ), νi − Φ(y − x )h∇u, νi]dS (y ). (17)
Vε ∂Vε
Además observamos que
Z  Z
Φ(y − x )h∇u, νidS (y ) ≤
 
 Φ(y − x )|∇u|dS (y )
∂B (x ,ε)  ∂B (x ,ε)
Z
= Φ(z )|∇u|dS (z )
∂B (0,ε)
≤ C máx |Φ|εn−1 → 0, cuando ε → 0.
∂B (0,ε)

Por otro lado

Z Z * +
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 |
Z
1
= u(y )dS (y ) → u(x ), cuando ε → 0.
nα(n)εn−1 ∂B (x ,ε)

Saul Ancari Villca EDP 21 de septiembre de 2021 30 / 34

Observe que ∂Vε = ∂U ∪ ∂B (x , ε). Haciendo ε → 0 en (17) obtenemos
Z
u(x ) = [Φ(y − x )h∇u, νi − u(y )h∇y Φ(y − x ), νi]dS (y ) (18)
∂U
Z
− Φ(y − x )∆u(y )dy .
U

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

Remarca. 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 .

Saul Ancari Villca EDP 21 de septiembre de 2021 31 / 34

Observe que ∂Vε = ∂U ∪ ∂B (x , ε). Haciendo ε → 0 en (17) obtenemos
Z
u(x ) = [Φ(y − x )h∇u, νi − u(y )h∇y Φ(y − x ), νi]dS (y ) (18)
∂U
Z
− Φ(y − x )∆u(y )dy .
U

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

Remarca. 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 .

Saul Ancari Villca EDP 21 de septiembre de 2021 31 / 34

Observe que ∂Vε = ∂U ∪ ∂B (x , ε). Haciendo ε → 0 en (17) obtenemos
Z
u(x ) = [Φ(y − x )h∇u, νi − u(y )h∇y Φ(y − x ), νi]dS (y ) (18)
∂U
Z
− Φ(y − x )∆u(y )dy .
U

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

Remarca. 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 .

Saul Ancari Villca EDP 21 de septiembre de 2021 31 / 34

Usando la formula de Green para integrales, obtenemos
Z Z
− φ ∆udy =
x
[uh∇φx , νi − φx (y )h∇u, νi]dS (y ) (19)
U ∂U
Z
= [uh∇φx , νi − Φ(y − x )h∇u, νi]dS (y )
∂U

Reemplazando (19) en (18) conseguimos

Z
u(x ) = − u(y )[h∇y Φ(y − x ), νi − h∇φx , νi]dS (y ) (20)
∂U
Z
− [Φ(y − x ) − φx (y )]∆u(y )dy .
U

Ahora supongamos que u ∈ C 2 (Ū ) resuelve

(
−∆u = f en U
(21)
u=g sobre ∂U ,

para funciones f , g continuas.

Saul Ancari Villca EDP 21 de septiembre de 2021 32 / 34
Usando la formula de Green para integrales, obtenemos
Z Z
− φ ∆udy =
x
[uh∇φx , νi − φx (y )h∇u, νi]dS (y ) (19)
U ∂U
Z
= [uh∇φx , νi − Φ(y − x )h∇u, νi]dS (y )
∂U

Reemplazando (19) en (18) conseguimos

Z
u (x ) = − u(y )[h∇y Φ(y − x ), νi − h∇φx , νi]dS (y ) (20)
∂U
Z
− [Φ(y − x ) − φx (y )]∆u(y )dy .
U

Ahora supongamos que u ∈ C 2 (Ū ) resuelve

(
−∆u = f en U
(21)
u=g sobre ∂U ,

para funciones f , g continuas.

Saul Ancari Villca EDP 21 de septiembre de 2021 32 / 34
Usando la formula de Green para integrales, obtenemos
Z Z
− φ ∆udy =
x
[uh∇φx , νi − φx (y )h∇u, νi]dS (y ) (19)
U ∂U
Z
= [uh∇φx , νi − Φ(y − x )h∇u, νi]dS (y )
∂U

Reemplazando (19) en (18) conseguimos

Z
u (x ) = − u(y )[h∇y Φ(y − x ), νi − h∇φx , νi]dS (y ) (20)
∂U
Z
− [Φ(y − x ) − φx (y )]∆u(y )dy .
U

Ahora supongamos que u ∈ C 2 (Ū ) resuelve

(
−∆u = f en U
(21)
u=g sobre ∂U ,

para funciones f , g continuas.

Saul Ancari Villca EDP 21 de septiembre de 2021 32 / 34
Entonces, usando (20), u tiene la siguiente expresión
Z
u (x ) = − g (y )[h∇y Φ(y − x ) − ∇φx (y ), νi]dS (y ) (22)
∂U
Z
+ [Φ(y − x ) − φx (y )]f (y )dy .
U
Motivado a esto, definimos la función de Green.
Definición 6
La función de Green para la región U es definida por

G (x , y ) := Φ(y − x ) − φx (y ) (x , y , x , y ).

Teorema 10 (Fórmula de representación usando la función de Green)

Si u ∈ C 2 (Ū ) resuelve el problema (21), entonces
Z Z
u(x ) = − g (y )h∇y G (x , y ), ν(y ) + f (y )G (x , y )dy .
∂U U

Saul Ancari Villca EDP 21 de septiembre de 2021 33 / 34

Entonces, usando (20), u tiene la siguiente expresión
Z
u (x ) = − g (y )[h∇y Φ(y − x ) − ∇φx (y ), νi]dS (y ) (22)
∂U
Z
+ [Φ(y − x ) − φx (y )]f (y )dy .
U
Motivado a esto, definimos la función de Green.
Definición 6
La función de Green para la región U es definida por

G (x , y ) := Φ(y − x ) − φx (y ) (x , y , x , y ).

Teorema 10 (Fórmula de representación usando la función de Green)

Si u ∈ C 2 (Ū ) resuelve el problema (21), entonces
Z Z
u(x ) = − g (y )h∇y G (x , y ), ν(y ) + f (y )G (x , y )dy .
∂U U

Saul Ancari Villca EDP 21 de septiembre de 2021 33 / 34

Entonces, usando (20), u tiene la siguiente expresión
Z
u (x ) = − g (y )[h∇y Φ(y − x ) − ∇φx (y ), νi]dS (y ) (22)
∂U
Z
+ [Φ(y − x ) − φx (y )]f (y )dy .
U
Motivado a esto, definimos la función de Green.
Definición 6
La función de Green para la región U es definida por

G (x , y ) := Φ(y − x ) − φx (y ) (x , y , x , y ).

Teorema 10 (Fórmula de representación usando la función de Green)

Si u ∈ C 2 (Ū ) resuelve el problema (21), entonces
Z Z
u(x ) = − g (y )h∇y G (x , y ), ν(y ) + f (y )G (x , y )dy .
∂U U

Saul Ancari Villca EDP 21 de septiembre de 2021 33 / 34

Bibliografı́a

Lawrence C Evans. “Partial differential equations”. En: Graduate studies in mathematics

19.4 (1998), pág. 7.

Saul Ancari Villca EDP 21 de septiembre de 2021 34 / 34