EDP
Piero Rodriguez Garcı́a
November 21, 2023
1 Problemas de Valor Inicial con Frontera en un
Dominio Infinito
En esta última sección del capı́tulo, consideraremos problemas de valor inicial-
frontera para la ecuación del calor en el dominio 0 < x < ∞, t > 0. Las
soluciones de nuestros problemas se obtendrán mediante la transformada de
Fourier seno o coseno, según lo apropiado para este tipo de problema. Primero
consideramos el problema
(7.1) ut − kuxx = 0 (x > 0, t > 0)
(7.2) u(x, 0) = f (x) (x ≥ 0)
(7.3) u(t, 0) = 0 (t ≥ 0)
que corresponde al problema del flujo de calor en una barra seminfinita x > 0,
con una distribución inicial de temperatura descrita por la función f (x) para
x ≥ 0, y cuyo extremo en x = 0 se mantiene a temperatura cero. Para resolver
este problema, utilizaremos la transformada de Fourier seno. Esta elección se
basa en el hecho de que el problema de autovalor singular relacionado con el
problema en cuestión es
X ′′ + λX = 0 (Xacotado cuando x → ∞)
que tiene como autofunciones Xs (x) = sin sx correspondientes a los auto-
valores λ = s2 , s > 0. En consecuencia, asumimos que el problema (7.1), (7.2),
(7.3) tiene una solución u que tiene la propiedad de que u y sus derivadas par-
ciales hasta segundo orden son todas continuas por tramos y absolutamente
integrables en 0 < x < ∞ para t > 0. Luego, la transformada de Fourier seno
de u,
21 Z ∞
2
(7.4) Us (s, t) = u(x, t) sin(sx) dx
π 0
existe y
1
21 Z ∞
2
(7.5) u(x, t) = Us (s, t) sin(sx) ds
π 0
Aquı́, Us no significa la derivada parcial de U con respecto a s. Siguiendo el
procedimiento descrito en la sección anterior, procedemos a encontrar la trans-
formada Us de nuestra solución. Una vez que se encuentra esta función, nuestra
solución se da por (7.5).
Al diferenciar (7.4) con respecto a t bajo el signo integral y usando la
ecuación diferencial (7.1), tenemos
21 Z ∞
∂Us (s, t) 2
= ut (x, t) sin(sx) dx
∂t π 0
12 Z ∞
2
= k uxx (x, t) sin(sx) dx
π 0
Después de integrar por partes dos veces, esto da
12
∂Us (s, t) 2 ∞
=k (ux sin(sx) − su cos(sx))|0
∂t π
12 Z ∞
2
(7.6) = −ks2 u(x, t) sin(sx) dx
π 0
Si asumimos, como antes, que u y ux tienden a cero cuando x → ∞, entonces,
debido a la condición frontera (7.3), los términos de frontera en el lado derecho
de (7.6) desaparecen, y mediante (7.4) llegamos a la ecuación diferencial
∂Us (s, t)
(7.7) + ks2 Us (s, t) = 0
∂t
para la función Us nosotros supongamos que la función f en (7.2) es suave por
partes y absolutamente integrable en (0, oo). Entonces, a partir de la condición
inicial (7.2) y (7.4), encontramos
21 Z ∞
2
Us (s, 0) = u(x, 0) sin(sx) dx
π 0
12 Z ∞
2
(7.8) Us (s, t) = f (x) sin(sx) dx
π 0
= Fs (s)
que es la transformada de Fourier seno de la función f . La solución de la
ecuación (7.7) que satisface la condición inicial (7.8) es
2
(7.9) Us (s, t) = Fs (s)e−ks t
2
Por lo tanto, mediante (7.5) y (7.8), obtenemos
12 Z ∞
2 2
u(x, t) = Fs (s)e−ks t sin(sx) ds
π 0
Z ∞Z ∞
2 2
(7.10) f (ξ)e−ks t sin(sξ) sin(sx) dξ ds
π 0 0
Esto proporciona la solución del problema (7.1), (7.2), (7.3).
Si reemplazamos la función sin s en la integranda de la integral en (7.10) por
1
[cos s(x − ξ) − cos s(x + ξ)]
2
y cambiamos formalmente el orden de integración, podemos escribir (7.10)
en la forma
Z ∞
(x − ξ)2 (x + ξ)2
1
(7.11) u(x, t) = √ exp(− ) − exp(− ) f (ξ) dξ
2πk 0 4kt 4kt
en la cual hemos utilizado el resultado (6.13). En esta forma, se puede
verificar que u es la solución de nuestro problema bajo la condición más débil
de que f sea continua y acotada para x ≥ 0 con f (0) = 0.
Observamos que la solución (7.11) de nuestro problema (7.1), (7.2), (7.3)
también se puede obtener extendiendo el dato inicial f (7.2) como una función
impar para −∞ < x < ∞ y luego resolviendo el problema de valor inicial resul-
tante mediante la fórmula (6.14). [Ver Problema 4, Ejercicios 6.5.(Demuestre
que si la función f(x) (6.2) [u(x,0)=f(x)] es par (o impar), entonces la solución
(6.14) también es par (o impar) con respecto a la variable x. Por lo tanto,
deducimos que u(0,t) = 0 (o u x(O,t) = 0))]
Usando el hecho de que f (−x) = −f (x), la solución dada por (6.14) se re-
duce entonces a (7.11). Observamos que este es precisamente el método que em-
pleamos en la solución de un problema de valor inicial-frontera para la ecuación
de la onda. El método funciona porque tanto la ecuación de la onda como la
ecuación del calor permanecen sin cambios en su forma cuando la variable x se
reemplaza por −x.
De la misma manera, podemos encontrar la solución del problema (7.1),
(7.2) con la condición frontera
(7.12) ux (x, t) = 0 (t ≥ 0)
extendiendo la función como una función par para −ω < x < ∞. En este
caso, la solución está dada por
Z ∞
(x − ξ)2 (x + ξ)2
1
(7.13) u(x, t) = √ exp(− ) + exp(− ) f (ξ) dξ
2πk 0 4kt 4kt
Destacamos aquı́ que para resolver el problema (7.1), (7.2), (7.12) mediante
el método de la transformada de Fourier, será necesario emplear la transformada
3
de coseno de Fourier de la solución u. Esto también es cierto incluso en el prob-
lema más difı́cil donde la condición frontera (7.12) es no homogénea. Ilustramos
el procedimiento en relación con el problema
(7.14) ut − kuxx = 0
u(x, 0) = 0
ux (0, t) = h(t)
Como es habitual, suponemos que este problema tiene una solución que
es continua por tramos y absolutamente integrable en (0, ∞), incluyendo sus
derivadas de primer y segundo orden. Consideramos la transformada de coseno
de Fourier de u,
21 Z ∞
2
(7.15) Uc (s, t) = u(x, t) cos(sx) dx
π 0
Diferenciando esto con respecto a t y procediendo como en casos anteriores,
tenemos
21 Z ∞
∂Uc (s, t) 2
= ut (x, t) cos(sx) dx
∂t π 0
12 Z ∞
2
= k uxx (x, t) cos(sx) dx
π 0
12
2 ∞
=k (ux cos(sx) + su sin(sx))|0
π
21 Z ∞
2
= −ks2 u(x, t) cos(sx) dx
π 0
Bajo la suposición de que u y ux tienden a cero cuando x → ∞, y teniendo
en cuenta la condición frontera en (7.14), esto lleva a la ecuación diferencial no
homogénea
21
∂Uc (s, t) 2 2
(7.16) + ks Uc (s, t) = −k h(t)
∂t π
La solución de esta ecuación que, por la condición inicial en (7.14), satisface
la condición Uc (s, 0) = 0
es
21 Z t
2 2
(7.17) Uc (s, t) = −k e−ks (t−τ ) h(τ ) dτ
π 0
Por lo tanto, tomando la transformada inversa, obtenemos
21 Z ∞
2
u(x, t) = Uc (s, t) cos(sx) ds
π 0
4
Z ∞ Z t
2k 2
(7.18) =− e−ks (t−τ )
h(τ ) cos(sx) dτ ds
π 0 0
Si intercambiamos formalmente el orden de integración y usamos el resultado
en (6.13), esta fórmula se simplifica a la forma
12 Z t −x2
k exp 4k(t−τ )
(7.19) u(x, t) = − √ h(τ ) dτ
π 0 t − τ
Se puede demostrar que esto proporciona la solución del problema (7.14)
bajo la suposición de que h es diferenciable para t ≥ 0.
Destacamos que no podemos verificar que (7.19) cumple con la condición
frontera de (7.14) simplemente diferenciando (7.19) bajo el signo integral y
estableciendo x = 0. De hecho, tal manipulación formal conducirı́a al resultado
falso ux (0, t) = 0. Ilustramos la situación en este caso mediante el siguiente
ejemplo.
1.1 Ejemplo 7.1
Halla la solucion del problema
ut − kuxx = 0 (x > 0, t > 0)
u(x, 0) = f (x) (x ≥ 0)
u(t, 0) = 1 (t ≥ 0)
Solution: Por la formula (7.19) tenemos que:
12 Z t −x2
1 exp 4k(t−τ )
u(x, t) = − √ dτ
π 0 t−τ
√
Sea z= x/ 4t − τ entonces:
∞ 2
e−(z )
Z
x
u(x, t) = − √ √ , dz
π x/ 4t z2
por lo que se ve fácilmente que
2 2
!
∞
e−(z )
e−x /4t
Z
1 4t
ux (x, t) = − √ √ , dz +
π x/ 4t z2 π x
2
! 2
! 2
!
e−x /4t e−x /4t e−x /4t
4t 4t
uxx (x, t) = − − √
π x π x πt
2
!
e−x /4t
ut (x, t) = √
πt
5
For x > 0 y t > 0. Ası́, ut − ux = 0, mostrando que u satisface la ecuación del
calor para x > 0, t > 0. A medida que t → 0, es claro que u(x, t) → 0. Ahora,
para verificar la condición frontera ux (0, t) = 1, realizamos una integración por
partes para poner u en la forma
Z ∞
x 1 −z2 ∞ √ 2x −(z 2 )
u(x, t) = − √ (− e )|x/ 4t + √ √ e dz
π z π x/ 4t
21 Z ∞
4t −x2 /4t 2x 2
=− e + √ √ e−(z ) dz
π π x/ 4t
Entonces Z ∞
2x 2
ux (x, t) = √ √ e−(z ) dz
π x/ 4t
R∞ 2 p
lo que inmediatamente da ux (0, t) = 1. Desde que 0
e−(z ) dz = π/2
1.2 Ejercicio 1
ut − kuxx = 0 (x > 0, t > 0)
u(x, 0) = xe−x (x ≥ 0)
u(0, t) = 0 (t ≥ 0)
Solución Paso a Paso: 1. Utilizamos la fórmula (7.11) del texto:
Z ∞
(x − ξ)2 (x + ξ)2
1
u(x, t) = √ exp − − exp − (ξe−ξ ) dξ
2πk 0 4kt 4kt
2. Sustituimos u(x, 0) = xe−x en la fórmula:
Z ∞
1
u(x, 0) = √ (ξe−ξ ) dξ
2πk 0
3. Simplificamos la integral y resolvemos:
Z ∞
1
u(x, 0) = √ ξe−ξ dξ
2πk 0
1
=√ · Γ(2)
2πk
4. La solución final es:
Γ(2)
u(x, t) = √
2πk
6
1.3 Ejercicio 2
ut − kuxx = 0 (x > 0, t > 0)
u(x, 0) = e−ax (a > 0, x ≥ 0)
ux (0, t) = 0 (t ≥ 0)
Solución Paso a Paso: 1. Utilizamos la fórmula (7.19) del texto:
12 Z t exp − x2
1 4k(t−τ )
u(x, t) = − √ · e−aτ dτ
π 0 t−τ
2. Sustituimos u(x, 0) = e−ax en la fórmula:
12 Z t exp − x2
1 4k(t−τ )
u(x, t) = − √ · e−aτ dτ
π 0 t−τ
q 3. Resolvemos la integral indefinida utilizando el cambio de variable z =
x2
4k(t−τ ) :
12 Z t
1 x
exp −z 2 · e−aτ ·
u(x, t) = − dτ
π 0 2kz
4. La solución final es:
∞ 2
e−z
Z
x
u(x, t) = − √ √ dz
π x/ 4t z2
1.4 Ejercicio 3
ut − kuxx = 0 (x > 0, t > 0)
u(x, 0) = 0 (x > 0)
u(0, t) = A (Aconstante, t > 0)
Solución Paso a Paso: 1. Utilizamos la fórmula (7.13) del texto:
Z ∞
(x − ξ)2 (x + ξ)2
1
u(x, t) = √ exp − + exp − f (ξ) dξ
2πk 0 4kt 4kt
2. Sustituimos u(x, 0) = 0 en la fórmula:
Z ∞ 2 2
1 x x
u(x, 0) = √ exp − + exp − · 0 dξ
2πk 0 4k 4k
3. La integral se simplifica a cero y la solución final es:
u(x, t) = 0
7
1.5 Ejercicio 4
ut − kuxx = 0 (x > 0, t > 0)
u(x, 0) = 0 (x ≥ 0)
u(0, t) = h(t) (t ≥ 0)
Solución Paso a Paso: 1. Utilizamos la fórmula (7.17) del texto:
12 Z t
2 2
Uc (s, t) = −k e−ks (t−τ ) h(τ ) dτ
π 0
2. Sustituimos u(x, 0) = 0 en la fórmula:
12 Z t
2 2
Uc (s, t) = −k e−ks (t−τ ) · 0 dτ
π 0
3. La integral se simplifica a cero y la solución final es:
Uc (s, t) = 0
