Tarea #7

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 16

UNIVERSIDAD NACIONAL DE TRUJILLO

FACULTAD DE INGENIERÍA
ESCUELA DE INGENIERÍA MECÁNICA

TAREA N° 7

CURSO:
ECUACIONES DIFERENCIALES

ALUMNOS:
ESQUIVEL CASTAÑEDA HAMNER NIKOLAE
GUTIERREZ RISCO EVERTH LEONEL
PAREDES MURGA ALEXIS PAUL
PINTADO CASEUX AARON WILFREDO
QUIROZ ALVA CRISTIAN DANIEL
SOLANO HERRERA JOE
VEGA BALTODANO DIEGO ALBERTO

DOCENTE:
VARGAS VERA, LUIS FERNANDO

CICLO Y SECCIÓN:
IV CICLO – SECCIÓN A

Trujillo – Perú
2022
Sección 11.2
Ejercicio 5: Encuentre la serie de Fourier de 𝑓 en el intervalo dado.
0, −𝜋 < 𝑥 < 0
𝑓(𝑥) = { 2
𝑥 , 0≤𝑥<𝜋
Solución:
Sabemos que si tenemos la función 𝑓(𝑥) con periodo 2𝐿 y definida en el intervalo
−𝐿 < 𝑥 < 𝐿 entonces la serie de Fourier viene dada por

𝑎0 𝑚𝜋𝑥 𝑚𝜋𝑥
𝑓(𝑥) = + ∑ (𝑎𝑚 cos + 𝑏𝑚 sin )
2 𝐿 𝐿
𝑚=1

Donde
1 𝐿
𝑎0 = ∫ 𝑓(𝑥)𝑑𝑥
𝐿 −𝐿
1 𝐿 𝑚𝜋𝑥
𝑎𝑚 = ∫ 𝑓(𝑥) cos 𝑑𝑥 , 𝑚 = 1, 2, …
𝐿 −𝐿 𝐿
1 𝐿 𝑚𝜋𝑥
𝑏𝑚 = ∫ 𝑓(𝑥) sin 𝑑𝑥 , 𝑚 = 1, 2, …
𝐿 −𝐿 𝐿
Ahora, evaluamos los coeficientes 𝑎0 , 𝑎𝑚 y 𝑏𝑚 para la función dada. Ahora vamos a
evaluar 𝑎0 reemplazando 𝐿 = 𝜋.
1 𝜋
𝑎0 = ∫ 𝑓(𝑥)𝑑𝑥
𝜋 −𝜋
1 0 1 𝜋
𝑎0 = ∫ 0𝑑𝑥 + ∫ 𝑥 2 𝑑𝑥
𝜋 −𝜋 𝜋 0
1 3𝜋
𝑎0 = [𝑥 ]0
3𝜋
1 3
𝑎0 = 𝜋
3𝜋
𝜋2
𝑎0 =
3

Ahora vamos a evaluar 𝑎𝑚


1 𝜋 𝑚𝜋𝑥
𝑎𝑚 = ∫ 𝑓(𝑥) cos 𝑑𝑥
𝜋 −𝜋 𝜋
1 𝜋 2
𝑎𝑚 = ∫ 𝑥 cos 𝑚𝑥 𝑑𝑥
𝜋 0
𝜋
1 𝑥2 2 𝜋
𝑎𝑚 = {[ sin 𝑚𝑥] − ∫ 𝑥 sin 𝑚𝑥 𝑑𝑥}
𝜋 𝑚 0
𝑚 0
2 𝑥 𝜋 1
𝑎𝑚 = − {[− cos 𝑚𝑥] + 2 [sin 𝑚𝑥]𝜋0 }
𝑚𝜋 𝑚 0 𝑚
2
𝑎𝑚 = cos 𝑚𝜋
𝑚2
2
𝑎𝑚 = 2 (−1)𝑚
𝑚
Ahora vamos a evaluar 𝑏𝑚
1 𝜋 𝑚𝜋𝑥
𝑏𝑚 = ∫ 𝑓(𝑥) sin 𝑑𝑥
𝜋 −𝜋 𝜋
1 𝜋 2
𝑏𝑚 = ∫ 𝑥 sin 𝑚𝑥 𝑑𝑥
𝜋 0
𝜋
1 𝑥2 2 𝜋
𝑏𝑚 = {[− cos 𝑚𝑥] + ∫ 𝑥 cos 𝑚𝑥 𝑑𝑥}
𝜋 𝑚 0
𝑚 0
𝜋 2 𝑥 𝜋 1
𝑏𝑚 = − cos 𝑚𝜋 + {[ sin 𝑚𝑥] + 2 [cos 𝑚𝑥]𝜋0 }
𝑚 𝑚𝜋 𝑚 0 𝑚
𝜋 2
𝑏𝑚 = − cos 𝑚𝜋 + 3 (cos 𝑚𝜋 − 1)
𝑚 𝑚 𝜋
𝜋 2
𝑏𝑚 = − (−1)𝑚 + 3 ((−1)𝑚 − 1)
𝑚 𝑚 𝜋
Por lo tanto, la serie de Fourier para la función dada es

𝑎0 𝑚𝜋𝑥 𝑚𝜋𝑥
𝑓(𝑥) = + ∑ (𝑎𝑚 cos + 𝑏𝑚 sin )
2 𝜋 𝜋
𝑚=1
∞ ∞
𝜋2 2 𝜋 2
𝑓(𝑥) = + ∑ ( 2 (−1)𝑚 ) cos 𝑚𝑥 + ∑ (− (−1)𝑚 + 3 ((−1)𝑚 − 1)) sin 𝑚𝑥
6 𝑚 𝑚 𝑚 𝜋
𝑚=1 𝑚=1

Ahora, tenemos que proporcionar el número al que converge la serie de Fourier en un


punto de discontinuidad de 𝑓. De la definición de 𝑓(𝑥), podemos concluir que 𝑓(𝑥) es
continua en el intervalo (−𝜋, 𝜋) y en los extremos del intervalo, converge al intervalo
𝑓(𝜋 − ) + 𝑓(−𝜋 + )
2
Además
𝑓(𝜋 − ) = lim− 𝑥 2 = 𝜋 2
𝑥→𝜋

Y
𝑓(−𝜋 + ) = lim + 0 = 0
𝑥→−𝜋

Por lo tanto, la serie de Fourier de 𝑓 en los extremos converge al promedio


𝑓(𝜋 − ) + 𝑓(−𝜋 + ) 𝜋 2 + 0 𝝅𝟐
= =
2 2 𝟐
Sección 12.2
Ejercicio 3: Una varilla de longitud 𝐿 coincide con el intervalo [0, 𝐿] en el eje x.
Establezca el problema con valores en la frontera para la temperatura 𝑢(𝑥, 𝑡). El extremo
izquierdo se mantiene a una temperatura de 100° y hay una transferencia de calor del
extremo derecho al medio que lo rodea a temperatura cero. La temperatura inicial es 𝑓(𝑥)
en toda la varilla.
Solución:
El problema pertenece a la ecuación de calor para la barra de longitud “L”, de este modo
la ecuación diferencial parcial que definirá el problema planteado es:
𝜕 2𝑢 𝜕𝑢
𝐾 2
= , 0 < 𝑥 < 𝐿, 𝑡 > 0 … (1)
𝜕𝑥 𝜕𝑡
Tenemos que la ecuación (1) está sujeta a los siguientes datos dados por el ejercicio, que
como ya sabemos se les llama condiciones de contorno.
Ahora, como el extremo izquierdo se mantiene a una temperatura de 100°, ese es el final
𝑥 = 0 tiene la temperatura constante 100° → 𝑢(0, 𝑡) = 100 … (2)
Además, el calor se transfiere desde el extremo derecho, es decir 𝑥 = 𝐿 al medio
circundante a la temperatura 0. La transferencia de calor viene dada por ℎ[𝑢(𝐿, 𝑡) −
𝑢𝑚 ] = ℎ𝑢(𝐿, 𝑡), donde 𝑢𝑚 es la temperatura del entorno. Aquí 𝑢𝑚 = 0. Entonces la
transferencia de calor es ℎ[𝑢(𝐿, 𝑡) − 0] = ℎ𝑢(𝐿, 𝑡),
𝜕𝑢
Por lo tanto, conduce a | = −ℎ 𝑢(𝐿, 𝑡) … (3)
𝜕𝑡 𝑥=𝐿

Siendo la temperatura inicial es 𝑓(𝑥) a lo largo de toda la varilla, esto implica


𝑢(𝑥, 0) = 𝑓(𝑥), 0 < 𝑥 < 𝐿 … (4)
Por lo tanto, el problema del valor en la frontera es:
𝜕 2𝑢 𝜕𝑢
𝐾 = , 0 < 𝑥 < 𝐿, 𝑡 > 0
𝜕𝑥 2 𝜕𝑡
𝜕𝑢
𝑢(0, 𝑡) = 100, | = −ℎ 𝑢(𝐿, 𝑡), 𝑡 > 0
𝜕𝑥 𝑥=𝐿
𝑢(𝑥, 0) = 𝑓(𝑥), 0 < 𝑥 < 𝐿
Sección 12.3
Ejercicio 5: Suponga que se pierde calor desde la superficie lateral de una varilla
delgada de longitud 𝐿 dentro del medio circundante a temperatura cero. Si se aplica la
ley lineal de transferencia de calor, entonces la ecuación de calor toma la forma

𝜕 2𝑢 𝜕𝑢
𝑘 2
− ℎ𝑢 =
𝜕𝑥 𝜕𝑡

0 < 𝑥 < 𝐿, 𝑡 > 0, ℎ una constante. Encuentre la temperatura 𝑢(𝑥, 𝑡) si la temperatura


inicial es 𝑓(𝑥) en toda la varilla y los extremos 𝑥 = 0 y 𝑥 = 𝐿 están aislados. Vea la
figura 12.3.3.

Solución:
Tendremos:
𝑢(0, 𝑡) = 0 , 𝑢(𝐿, 𝑡) = 0 𝑡>0
𝑢(𝑥, 0) = 𝑓(𝑥) , 0<𝑥<𝐿
Utilizando el método de separación de variables obtenemos:
𝑢(𝑥, 𝑡) = 𝑋(𝑥)𝑇(𝑡)
𝜕 2𝑢
= 𝑋 ′′ (𝑥)𝑇(𝑡)
𝜕𝑥 2
𝜕𝑢
= 𝑋(𝑥)𝑇 ′ (𝑡)
𝜕𝑡
𝜕𝑢
= 𝑋′(𝑥)𝑇(𝑡)
𝜕𝑡
Donde al sustituir en la ecuación inicial:
𝑘𝑋 ′′ (𝑥)𝑇(𝑡) − ℎ𝑋(𝑥)𝑇(𝑡) = 𝑋(𝑥)𝑇 ′ (𝑡)
Realizaremos una división que nos permita darle forma a esta nueva ecuación:
𝑋 ′′ (𝑥)𝑇(𝑡) − ℎ𝑋(𝑥)𝑇(𝑡) 𝑋(𝑥)𝑇 ′ (𝑡)
=
𝑘𝑇(𝑡) 𝑘𝑇(𝑡)
ℎ𝑋(𝑥) 𝑋(𝑥)𝑇 ′ (𝑡)
𝑋 ′′ (𝑥) − =
𝑘 𝑘𝑇(𝑡)
1 ℎ𝑋(𝑥) 𝑇 ′ (𝑡)
(𝑋 ′′ (𝑥) − )=
𝑋(𝑥) 𝑘 𝑘𝑇(𝑡)
𝑋 ′′ (𝑥) ℎ𝑋(𝑥) 𝑇 ′ (𝑡)
− =
𝑋(𝑥) 𝑋(𝑥)𝑘 𝑘𝑇(𝑡)
𝑋 ′′ (𝑥) ℎ 𝑇 ′ (𝑡)
− =
𝑋(𝑥) 𝑘 𝑘𝑇(𝑡)
𝑋 ′′ (𝑥) 𝑇 ′ (𝑡) ℎ
= +
𝑋(𝑥) 𝑘𝑇(𝑡) 𝑘
𝑋 ′′ (𝑥) 1 𝑇 ′ (𝑡)
= ( + ℎ)
𝑋(𝑥) 𝑘 𝑇(𝑡)
Asumimos que:
𝑋 ′′ (𝑥) 1 𝑇 ′ (𝑡)
= ( + ℎ) = −𝜆2
𝑋(𝑥) 𝑘 𝑇(𝑡)
Tendremos entonces dos términos con el mismo valor a los cuales llamaremos:
𝑋 ′′ (𝑥)
= −𝜆2 (𝑖)
𝑋(𝑥)
1 𝑇 ′ (𝑡)
( + ℎ) = −𝜆2 (𝑖𝑖)
𝑘 𝑇(𝑡)
Para 𝑖:
𝑋 ′′ (𝑥)
= −𝜆2
𝑋(𝑥)
𝑋 ′′ (𝑥) + 𝜆2 𝑋(𝑥) = 0
Considerando un 𝜆 ≠ 0, la ecuación auxiliar será:
𝑚2 + 𝜆2 = 0
𝑚2 = −𝜆2
𝑚 = ±𝑖𝜆
Tendremos la solución:
𝑋(𝑥) = 𝐶1 cos(𝜆𝑥) + 𝐶2 𝑠𝑒𝑛(𝜆𝑥)
Ahora tomando 𝑖𝑖:
1 𝑇 ′ (𝑡)
( + ℎ) = −𝜆2
𝑘 𝑇(𝑡)
𝑇 ′ (𝑡) + (𝑘𝜆2 + ℎ)𝑇(𝑡) = 0
Considerando un 𝜆 ≠ 0:
𝑇 ′ (𝑡) = −(𝑘𝜆2 + ℎ)𝑇(𝑡)
𝑇 ′ (𝑡)
= −(𝑘𝜆2 + ℎ)
𝑇(𝑡)
Integrando esta ecuación tenemos:
𝑇 ′ (𝑡)
∫ 𝑑𝑡 = ∫ −(𝑘𝜆2 + ℎ) 𝑑𝑡
𝑇(𝑡)
ln(𝑇(𝑡)) = −(𝑘𝜆2 + ℎ)𝑡 + ln(𝑐)
ln(𝑇(𝑡)) − ln(𝑐) = −(𝑘𝜆2 + ℎ)𝑡
𝑇(𝑡)
ln ( ) = −(𝑘𝜆2 + ℎ)𝑡
𝑐
Tendremos la solución:
2 +ℎ)𝑡
𝑇(𝑡) = 𝐶𝑒 −(𝑘𝜆
Si aplicamos las condiciones iniciales del problema.
Para 𝑢(0, 𝑡) = 0:
𝑢(𝑥, 𝑡) = 𝑋(𝑥)𝑇(𝑡)
𝑢(0, 𝑡) = 𝑋(0)𝑇(𝑡) = 0
𝑋(0) = 0
𝑋 ′ (0) = 0
Sustituyendo en la solución de 𝑖:
𝑋(𝑥) = 𝐶1 cos(𝜆𝑥) + 𝐶2 𝑠𝑒𝑛(𝜆𝑥)
𝑋′(𝑥) = −𝜆𝐶1 sen(𝜆𝑥) + 𝜆𝐶2 𝑐𝑜𝑠(𝜆𝑥)
𝑋′(0) = −𝜆𝐶1 sen(0) + 𝜆𝐶2 𝑐𝑜𝑠(0)
𝑋 ′(0) = 0 + 𝜆𝐶2
Ya que 𝜆 ≠ 0:
𝐶2 = 0
Tenemos entonces:
𝑋(𝑥) = 𝐶1 cos(𝜆𝑥)
Para 𝑢(𝐿, 𝑡) = 0:
𝑢(𝑥, 𝑡) = 𝑋(𝑥)𝑇(𝑡)
𝑢(𝐿, 𝑡) = 𝑋(𝐿)𝑇(𝑡) = 0
𝑋(𝐿) = 0
𝑋 ′ (𝐿) = 0
Sustituyendo en la solución de 𝑖:
𝑋(𝑥) = 𝐶1 cos(𝜆𝑥) + 𝐶2 𝑠𝑒𝑛(𝜆𝑥)
𝑋′(𝑥) = −𝜆𝐶1 sen(𝜆𝑥) + 𝜆𝐶2 𝑐𝑜𝑠(𝜆𝑥)
𝑋′(𝐿) = −𝜆𝐶1 sen(𝜆𝐿) + 𝜆𝐶2 𝑐𝑜𝑠(𝐿)
Sabiendo que 𝐶2 = 0:
𝑋′(𝐿) = −𝜆𝐶1 sen(𝜆𝐿)
−𝜆𝐶1 sen(𝜆𝐿) = 0
Para una solución no trivial, se tiene que 𝐶1 ≠ 0
Por lo tanto:
sen(𝜆𝐿) = sen(𝑛𝜋)
𝜆𝐿 = 𝑛𝜋
𝑛𝜋
𝜆=
𝐿
𝑛𝜋
Para 𝜆, 𝐶1 ≠ 0, los autovalores son 𝜆𝑛 = , donde 𝑛 = 1,2,3, …
𝐿

La solución para 𝜆 ≠ 0 es de la forma:


𝑛𝜋
𝑋(𝑥) = 𝐶1 cos 𝑥
𝐿
Ahora veremos que la solución general está dada por:
𝑢(𝑥, 𝑡) = 𝑋(𝑥)𝑇(𝑡)
𝑛𝜋 2
𝑢(𝑥, 𝑡) = (𝐶1 cos 𝑥) (𝐶𝑒 −(𝑘𝜆 +ℎ)𝑡 )
𝐿
2 𝑛𝜋
𝑢(𝑥, 𝑡) = 𝑘𝑒 −(𝑘𝜆 +ℎ)𝑡 𝑐𝑜𝑠 𝑥
𝐿
Donde:
𝑘 = 𝐶1 . 𝐶
Para 𝜆 = 0, las ecuaciones son 𝑋 ′′ (𝑥) = 0 y 𝑇 ′ (𝑡) + ℎ𝑇(𝑡) = 0
Al integrar la ecuación anterior tenemos:
𝑋 ′ (𝑥) = 𝐷
Si volvemos a integrar:
𝑋(𝑥) = 𝐷𝑥 + 𝐸
Por lo tanto, tendremos una solución de la forma:
𝑋(𝑥) = 𝐷𝑥 + 𝐸
Al resolver 𝑇 ′ (𝑡) + ℎ𝑇(𝑡) = 0:
𝑇 ′ (𝑡) + ℎ𝑇(𝑡) = 0
𝑇 ′ (𝑡) = −ℎ𝑇(𝑡)
𝑇 ′ (𝑡)
= −ℎ
𝑇(𝑡)
Al integrar la ecuación obtenemos:
𝑇 ′ (𝑡)
∫ 𝑑𝑡 = ∫ −ℎ 𝑑𝑡
𝑇(𝑡)
𝑙𝑛𝑇(𝑡) = −ℎ𝑡 + 𝑙𝑛𝑓
𝑙𝑛𝑇(𝑡) − 𝑙𝑛𝑓 = −ℎ𝑡
𝑇(𝑡)
𝑙𝑛 = −ℎ𝑡
𝑓
𝑇(𝑡)
= 𝑒 −ℎ𝑡
𝑓
𝑇(𝑡) = 𝑓𝑒 −ℎ𝑡
Tendremos entonces una solución de la forma:
𝑇(𝑡) = 𝑓𝑒 −ℎ𝑡
Si aplicamos las condiciones iniciales del problema.
Para 𝑢(0, 𝑡) = 0:
𝑢(𝑥, 𝑡) = 𝑋(𝑥)𝑇(𝑡)
𝑢(0, 𝑡) = 𝑋(0)𝑇(𝑡)
𝑋(0) = 0
𝑋 ′ (0) = 0
Donde:
𝑋(𝑥) = 𝐷𝑥 + 𝐸
𝑋′(𝑥) = 𝐷
𝑋′(0) = 𝐷
Podemos decir que:
𝐷=0
Por lo tanto:
𝑋(𝑥) = 𝐸
Finalmente tendremos la solución general de la ecuación:
𝑢(𝑥, 𝑡) = 𝑋(𝑥)𝑇(𝑡)
2 +ℎ)𝑡 𝑛𝜋
𝑢(𝑥, 𝑡) = 𝑘𝑒 −(𝑘𝜆 𝑐𝑜𝑠 𝑥
𝐿

𝑘𝑛2 𝜋 2 𝑛𝜋𝑥
−( +ℎ)𝑡
𝑢(𝑥, 𝑡) = ∑ 𝑎𝑛 𝑒 𝐿2 cos ( )
𝐿
𝑛=1

𝑘𝑛2 𝜋 2 𝑛𝜋𝑥
−( )
𝑢(𝑥, 𝑡) = ∑ 𝑎𝑛 𝑒 𝐿 𝑒 −ℎ𝑡 cos ( )
𝐿
𝑛=1

𝑘𝑛2 𝜋 2 𝑛𝜋𝑥
−ℎ𝑡 −( )
𝑢(𝑥, 𝑡) = 𝑒 ∑ 𝑎𝑛 𝑒 𝐿 cos ( )
𝐿
𝑛=1

Donde:
2 𝐿 𝑛𝜋𝑥
𝑎𝑛 = ∫ 𝑓(𝑥)𝑠𝑒𝑛 ( ) 𝑑𝑥
𝐿 0 𝐿
Ejercicio 6: Resuelva el problema 5, si los extremos 𝑥 = 0 𝑦 𝑥 = 𝐿, se mantienen a
temperatura cero.

Solución:
Del ejercicio 5 tenemos la siguiente ecuación:
𝐾𝜕 2 𝑢 𝜕𝑢
2
− ℎ𝑢 = , 0 < 𝑥 < 𝐿, 𝑡 > 0 𝑑𝑜𝑛𝑑𝑒 ℎ 𝑒𝑠 𝑢𝑛𝑎 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒
𝜕𝑥 𝜕𝑡
Como nos dicen que los extremos 𝑥 = 0 𝑦 𝑥 = 𝐿, se mantienen a temperatura cero,
entonces, tenemos las condiciones:
𝑢(0, 𝑡) = 0, 𝑢(𝐿, 𝑡) = 0, 𝑝𝑎𝑟𝑎 𝑡 > 0
Y, además:
𝑢(𝑥, 0) = 𝑓(𝑥)
Por el método de separación de variables, denotamos:
𝑢(𝑥, 𝑡) = 𝑀(𝑥)𝑁(𝑡)
Entonces tendríamos:
𝜕 2𝑢
= 𝑀′′(𝑥)𝑁(𝑡)
𝜕𝑥 2
𝜕𝑢
= 𝑀(𝑥)𝑁′(𝑡)
𝜕𝑡
Reemplazando y operando:
𝐾𝑀′′𝑁 − ℎ𝑀𝑁 = 𝑀𝑁′
𝐾𝑀′′𝑁 ℎ𝑀𝑁 𝑀𝑁′
− =
𝑀 𝑀 𝑀
𝐾𝑀′′𝑁
− ℎ𝑁 = 𝑁′
𝑀
𝐾𝑀′′𝑁
= 𝑁 ′ + ℎ𝑁
𝑀
𝑀′′ 𝑁 ′ + ℎ𝑁
=
𝑀 𝐾𝑁
𝑀′′ 𝑁 ′ ℎ𝑁
= +
𝑀 𝐾𝑁 𝐾𝑁
𝑀′′ 𝑁 ′ ℎ
= +
𝑀 𝐾𝑁 𝐾
𝑀′′ 1 𝑁′
= ( ℎ + ) = −𝜆
𝑀 𝐾 𝑁
Por lo que tendríamos:
𝑀′′ + 𝜆𝑀 = 0 y 𝑁 ′ + (ℎ + 𝐾𝜆)𝑁 = 0
Ahora, resolveremos estas ecuaciones para cuando: 𝜆 = 0, 𝜆 = −𝛼 2 < 0 𝑦 𝜆 = 𝛼 2 > 0
Para 𝜆 = 0:
𝑀′′ + (0)𝑀 = 0
𝑀′′ = 0
𝑀 = 𝑐1 𝑥 + 𝑐2
Y
𝑁 ′ + (ℎ + 𝐾(0))𝑁 = 0
𝑁 ′ + ℎ𝑁 = 0
𝑁 = 𝑐2 𝑒 −ℎ𝑡
Entonces, 𝑢(𝑥, 𝑡) sería:
𝑢(𝑥, 𝑡) = 𝑀(𝑥)𝑁(𝑡) = (𝑐1 𝑥 + 𝑐2 ) (𝑐3 𝑒 −ℎ𝑡 )
Para 𝜆 = −𝛼 2 :
𝑀′′ − 𝛼 2 𝑀 = 0
𝑀 = 𝑐4 cosh(𝛼𝑥) + 𝑐5 sinh(𝛼𝑥)
Y
𝑁 ′ + (ℎ − 𝛼 2 𝐾)𝑁 = 0
2 𝐾)𝑡
𝑁 = 𝑐6 𝑒 −(ℎ−𝛼
Entonces, 𝑢(𝑥, 𝑡) sería:
2 𝐾)𝑡
𝑢(𝑥, 𝑡) = 𝑀(𝑥)𝑁(𝑡) = (𝑐4 cosh(𝛼𝑥) + 𝑐5 sinh(𝛼𝑥)) (𝑐6 𝑒 −(ℎ−𝛼 )
Para 𝜆 = 𝛼 2 :
𝑀′′ + 𝛼 2 𝑀 = 0
𝑀 = 𝑐7 cos(𝛼𝑥) + 𝑐8 sin(𝛼𝑥)
Y
𝑁 ′ + (ℎ + 𝛼 2 𝐾)𝑁 = 0
2 𝐾)𝑡
𝑁 = 𝑐9 𝑒 −(ℎ+𝛼
Entonces, 𝑢(𝑥, 𝑡) sería:
2 𝐾)𝑡
𝑢(𝑥, 𝑡) = 𝑀(𝑥)𝑁(𝑡) = (𝑐7 cos(𝛼𝑥) + 𝑐8 sin(𝛼𝑥)) (𝑐9 𝑒 −(ℎ+𝛼 )
Analizamos por casos:
Reemplazamos las condiciones de frontera 𝑢(0, 𝑡) = 0, 𝑢(𝐿, 𝑡) = 0 en (𝑐1𝑥 +
𝑐2 ) (𝑐3 𝑒 −ℎ𝑡 )
𝑢(0, 𝑡) = (𝑐1 (0) + 𝑐2 ) (𝑐3 𝑒 −ℎ𝑡 ) = 0
𝑐2 = 0
𝑢(𝐿, 𝑡) = (𝑐1 (𝐿) + 0) (𝑐3 𝑒 −ℎ𝑡 ) = 0
𝑐1 = 0
Con esto, tendríamos la solución trivial, así que este caso no cumple.
Ahora reemplazamos las condiciones de frontera 𝑢(0, 𝑡) = 0, 𝑢(𝐿, 𝑡) = 0 en 𝑢(𝑥, 𝑡) =
2 𝐾)𝑡
𝑀(𝑥)𝑁(𝑡) = (𝑐4 cosh(𝛼𝑥) + 𝑐5 sinh(𝛼𝑥)) (𝑐6 𝑒 −(ℎ−𝛼 ):
2 𝐾)𝑡
𝑢(0, 𝑡) = (𝑐4 cosh(0) + 𝑐5 sinh(𝑥)) (𝑐6 𝑒 −(ℎ−𝛼 )=0
2 𝐾)𝑡
𝑢(0, 𝑡) = (𝑐4 (1) + 𝑐5 (0)) (𝑐6 𝑒 −(ℎ−𝛼 )=0
𝑐4 = 0
2 𝐾)𝑡
𝑢(𝐿, 𝑡) = ((0) + 𝑐5 sinh(𝛼𝐿) (𝑐6 𝑒 −(ℎ−𝛼 )=0
𝑐5 = 0 𝑜 𝑐6 = 0
Con lo que tenemos la solución trivial, por lo que este caso tampoco cumple.
Ahora reemplazamos las condiciones de frontera para (𝑐7 cos(𝛼𝑥) +
2 𝐾)𝑡
𝑐8 sin(𝛼𝑥)) (𝑐9 𝑒 −(ℎ+𝛼 ):
2 𝐾)𝑡
𝑢(0, 𝑡) = (𝑐7 cos(0) + 𝑐8 sin(0)) (𝑐9 𝑒 −(ℎ+𝛼 )=0
2 𝐾)𝑡
𝑢(0, 𝑡) = (𝑐7 ) (𝑐9 𝑒 −(ℎ+𝛼 )=0
𝑐7 = 0
2 𝐾)𝑡
𝑢(𝐿, 𝑡) = (0 cos(𝐿) + 𝑐8 sin(𝐿)) (𝑐9 𝑒 −(ℎ+𝛼 )=0
2 𝐾)𝑡
𝑢(𝐿, 𝑡) = (𝑐8 sin(𝛼𝐿)) (𝑐9 𝑒 −(ℎ+𝛼 )=0
𝑐8 = 0 𝑜 sin(𝛼𝐿) = 0
Como 𝑐8 = 0 da la solución trivial, entoncessin(𝛼𝐿) = 0.
Entonces:
sin(𝛼𝐿) = 0.
𝛼𝐿 = 𝑛𝜋
𝑛𝜋
𝛼=
𝐿
Reemplazamos y tenemos:
𝑛𝜋 𝑛𝜋 2
−(ℎ+( ) 𝐾)𝑡
𝑢(𝑥, 𝑡) = (𝑐8 sin ( 𝐿)) (𝑐9 𝑒 𝐿 )
𝐿
𝑛𝜋 𝑛𝜋 2
−(ℎ+( ) 𝐾)𝑡
𝑢(𝑥, 𝑡) = 𝑐𝑛 sin ( 𝐿) 𝑒 𝐿
𝐿
Por el principio de superposición:
𝑛𝜋 2 𝑛𝜋
−(ℎ+( ) 𝐾)𝑡
𝑢𝑛 = 𝐴𝑛 𝑒 𝐿 sin ( 𝑥)
𝐿

𝑛𝜋 2 𝑛𝜋
−(ℎ+( ) 𝐾)𝑡
𝑢(𝑥, 𝑡) = ∑ 𝐴𝑛 𝑒 𝐿 sin ( 𝑥)
𝐿
𝑛=1

Los 𝐴𝑛 tienen la forma de:


2 𝐿 𝑛𝜋
𝐴𝑛 = ∫ 𝑓(𝑥) sin ( 𝑥)
𝐿 0 𝐿
Tenemos la condición inicial:
𝑢(𝑥, 0) = 𝑓(𝑥)

𝑛𝜋
𝑢(𝑥, 0) = 𝑓(𝑥) = ∑ 𝐴𝑛 sin( 𝑥)
𝐿
𝑛=1

Por lo tanto, 𝑢(𝑥, 𝑡) es:



2 𝐿 𝑛𝜋 𝑛𝜋 2 𝑛𝜋
−(ℎ+( ) 𝐾)𝑡
𝑢(𝑥, 𝑡) = ∑ ( ∫ 𝑓(𝑥) sin ( 𝑥)) 𝑒 𝐿 sin ( 𝑥)
𝐿 0 𝐿 𝐿
𝑛=1
𝑛𝜋
Con: 𝑓(𝑥) = ∑∞
𝑛=1 𝐴𝑛 sin( 𝑥)
𝐿
Imagen Adjunta
Ejercicio 5: Resolver el problema de valor de frontera.
𝜕 2𝑢 𝜕 2𝑢
= + 𝑢, 𝑢(0, 𝑦) = 2𝑒 𝑦
𝜕𝑦𝜕𝑥 𝜕𝑦 2
Solución:
Sustituimos 𝑢(𝑥, 𝑦) = 𝑋(𝑥)𝑌(𝑦) en la EDP.
𝑋 ′ 𝑌 ′ = 𝑋𝑌 ′′ + 𝑋𝑌
𝑋 ′ 𝑌 ′ = 𝑋(𝑌 ′′ + 𝑌)
Separamos las variables y la igualamos a una constante ya que ambos lados son
independientes de 𝑥 como de 𝑦.
𝑋′ 𝑌 ′′ + 𝑌
= =𝜆
𝑋 𝑌′
𝑋 ′ − 𝜆𝑋 = 0 𝑌 ′′ + 𝑌 − 𝜆𝑌 ′ = 0
Para 𝑋
𝑋 ′ − 𝜆𝑋 = 0
Cambiamos a la ecuación característica.
𝑟−𝜆 =0
𝑟=𝜆
Su solución general.
𝑋 = 𝐶1 𝑒 𝜆𝑥
Para 𝑌
𝑌 ′′ + 𝑌 − 𝜆𝑌 ′ = 0
Cambiamos a la ecuación característica.
𝑟 2 − 𝜆𝑟 + 1 = 0
𝜆 ± √𝜆2 − 4
𝑟=
2
Su solución general.
𝜆+√𝜆2 −4 𝜆−√𝜆2 −4
𝑦 𝑦
𝑌 = 𝐶2 𝑒 2 + 𝐶3 𝑒 2

Ahora para 𝑢.
⇒ 𝑢(𝑥, 𝑦) = 𝑋𝑌
𝜆+√𝜆2 −4 𝜆−√𝜆2 −4
𝑢(𝑥, 𝑦) = 𝐶1 𝑒 𝜆𝑥
(𝐶2 𝑒 2 𝑦 + 𝐶3 𝑒 2 𝑦 )

Aplicamos la condición inicial.


𝑢(0, 𝑦) = 𝑋(0)𝑌(𝑦) = 2𝑒 𝑦
Reemplazamos.
𝜆+√𝜆2 −4 𝜆−√𝜆2 −4
𝐶1 𝑒 𝜆(0)
(𝐶2 𝑒 2 𝑦 + 𝐶3 𝑒 2 𝑦 ) = 2𝑒 𝑦

𝜆+√𝜆2 −4 𝜆−√𝜆2 −4
𝐶1 𝐶2 𝑒 2 𝑦 + 𝐶1 𝐶3 𝑒 2 𝑦 = 2𝑒 𝑦
Simplificamos constantes.
𝐴1 = 𝐶1 𝐶2 𝐵1 = 𝐶1 𝐶3
𝜆+√𝜆2 −4 𝜆−√𝜆2 −4
𝐴1 𝑒 2 𝑦 + 𝐵1 𝑒 2 𝑦 = 𝑎𝑒 𝑦 + 𝑏𝑒 𝑦
Donde:
𝑎+𝑏 = 2
Igualamos términos de cada lado.
𝜆+√𝜆2 −4
𝑦
⇒ 𝐴1 𝑒 2 = 𝑎𝑒 𝑦
𝜆−√𝜆2 −4
⇒ 𝐵1 𝑒 2 𝑦 = 𝑏𝑒 𝑦
Las constantes nos quedan:
𝐴1 + 𝐵2 = 2
Para el valor de 𝜆.
𝜆+√𝜆2 −4 𝜆−√𝜆2 −4
𝑦 𝑦
𝑒 2 = 𝑒𝑦 = 𝑒 2

𝜆=2
Finalmente 𝑢(𝑥, 𝑦) nos queda:
2+√22 −4 2−√22 −4
𝑦 𝑦
𝑢(𝑥, 𝑦) = 𝑒 2𝑥 (𝐴1 𝑒 2 + 𝐵1 𝑒 2 )

𝑢(𝑥, 𝑦) = 𝑒 2𝑥 (𝐴1 𝑒 𝑦 + 𝐵1 𝑒 𝑦 )
𝒖(𝒙, 𝒚) = 𝟐𝒆𝟐𝒙+𝒚

También podría gustarte