Tarea 5

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

Nombre: Rivaldo Mercado Código: 2019119032 Grupo: 2

1. Resolver

𝜕𝑢 𝜕 2𝑢
=𝑘 2 , 0 < 𝑥 < 𝐿, 𝑡>0
𝜕𝑡 𝜕𝑥
𝑢(0, 𝑡) = 𝑢(𝐿, 𝑡) = 0, 𝑡≥0
𝑢(𝑥, 0) = 𝑥(𝐿 − 𝑥), 0≤𝑥≤𝐿

Usando separación de variables tenemos que.

𝑢(𝑥, 𝑡) = 𝑋(𝑥 )𝑇(𝑡) (1)


Entonces

𝜕𝑢 𝜕2𝑢
= 𝑋𝑇 ′ , = 𝑋 ′′ 𝑇
𝜕𝑡 𝜕𝑥 2

Sustituimos y tenemos
′ ′′
𝑇 ′ 𝑋′′
𝑋𝑇 = 𝑘𝑋 𝑇 ⟹ =
𝑘𝑇 𝑋

𝑇 ′ 𝑋′′
= =𝛼
𝑘𝑇 𝑋

𝑇 ′ − 𝑘𝛼𝑇 = 0 , 𝑋 ′′ − 𝛼𝑋 = 0
.

𝑑𝑇
= 0 ⟹ 𝑇(𝑡) = 𝑐1
𝑑𝑡
𝑚2 = 0 ⟹ 𝑚 = 0

⟹ 𝑋(𝑥 ) = 𝑐2 + 𝑐3𝑥

Luego, sustituyendo en (1) se tiene que

𝑢(𝑥, 𝑡) = 𝑐1 (𝑐2 + 𝑐3 𝑥 )

Considerando la condición de frontera 𝑢(0, 𝑡) = 0. Tenemos que,


𝑢(𝑥, 𝑡) = 𝑐1(𝑐3 𝑥 ) o bien 𝑢(𝑥, 𝑡) = 𝐴𝑥 ∴ 𝐴 = 𝑐1𝑐3
Ahora, para 𝑢(𝐿, 𝑡) = 0 ⟹ 𝑢(𝐿, 𝑡) = 𝐴𝐿 = 0,

Caso II. 𝛼 > 0.

𝑑𝑇 𝑑𝑇
= 𝑘𝛼 𝑑𝑡 ⟹ ∫ = 𝑘𝛼 ∫ 𝑑𝑡
𝑇 𝑇

𝑇(𝑡) = 𝑐1 𝑒 𝑘𝛼𝑡

Para la segunda tenemos

𝑋 ′′ − 𝛼𝑋 = 0
La ecuación característica es
𝑚2 − 𝛼 = 0

𝑚 = ±√𝛼

Luego,
𝑋(𝑥 ) = 𝑐2𝑒 √𝛼 𝑥 + 𝑐3 𝑒 −√𝛼 𝑥
Entonces
𝑢(𝑥, 𝑡) = 𝑐1 𝑒𝑘𝛼𝑡 (𝑐2 𝑒√𝛼 𝑥 + 𝑐3 𝑒−√𝛼 𝑥 )

Usando la primera condición, tenemos que 𝑢(0, 𝑡) = 𝑐1 𝑒𝑘𝛼𝑡 (𝑐2 + 𝑐3 ) = 0 entonces

𝑢(𝑥, 𝑡) = 𝑐1 𝑒𝑘𝛼𝑡 (𝑐2 𝑒√𝛼 𝑥 − 𝑐2 𝑒−√𝛼 𝑥 ) ⟹ 𝑢(𝑥, 𝑡) = 𝑐1 𝑐2 𝑒𝑘𝛼𝑡 (𝑒√𝛼 𝑥 − 𝑒−√𝛼 𝑥 )

𝑢(𝑥, 𝑡) = 𝐴𝑒𝑘𝛼𝑡 (𝑒√𝛼 𝑥 − 𝑒−√𝛼 𝑥 ) ∴ 𝐴 = 𝑐1 𝑐2

De la segunda condición de frontera, se tiene que 𝑢(𝐿, 𝑡) = 𝐴𝑒𝑘𝛼𝑡 (𝑒3√𝛼 𝐿 − 𝑒−3√𝛼 𝐿 ) = 0


entoces,

𝑒3√𝛼 𝐿 − 𝑒−3√𝛼 𝐿 = 0 ⟹ 𝑒3√𝛼 𝐿 = 𝑒−3√𝛼 𝐿


⟹ 𝑒 6√𝛼 𝐿
=1

Puesto que 𝛼 > 0, 𝑒 6√𝛼 𝐿 > 1. Esto muestra que no podemos tener 𝛼 > 0.

Caso III. 𝛼 < 0.

𝑇 ′ 𝑋′′
= = −𝜆2
𝑘𝑇 𝑋
de donde se tiene que
𝑇′ 𝑋′′
= −𝜆2 (3) ; = −𝜆2 (4)
𝑘𝑇 𝑋
𝑑𝑇 𝑑𝑇
= −𝑘𝜆2 𝑑𝑡 ⟹ ∫ = −𝑘𝜆2 ∫ 𝑑𝑡
𝑇 𝑇
2𝑡
𝑇(𝑡) = 𝑐1 𝑒 −𝑘𝜆

𝑋 ′′ + 𝜆2 𝑋 = 0
La ecuación característica es
𝑚2 + 𝜆2 = 0

𝑚 = ±𝜆𝑖
Luego,
𝑋(𝑥 ) = 𝑐2 cos 𝜆𝑥 + 𝑐3 𝑠𝑒𝑛 𝜆𝑥
Luego,
2
𝑢(𝑥, 𝑡) = 𝑐1 𝑒 −𝑘𝜆 𝑡 (𝑐2 cos 𝜆𝑥 + 𝑐3 𝑠𝑒𝑛 𝜆𝑥 )

Tenemos que.

2
𝑢(𝑥, 𝑡) = 𝑐1𝑒 −𝑘𝜆 𝑡 (𝑐3 𝑠𝑒𝑛 𝜆𝑥 )
2
𝑢(𝑥, 𝑡) = 𝐴𝑒 −𝑘𝜆 𝑡 𝑠𝑒𝑛 𝜆𝑥 ∴ 𝐴 = 𝑐1 𝑐3

Entonces tenemos que.

𝑠𝑒𝑛 𝜆𝐿 = 0

𝜆𝐿 = 𝑛𝜋 𝑐𝑜𝑛 𝑛 = 0, ±1, ±2, ⋯

𝑛𝜋
𝜆=
𝐿
Luego,
𝑛 2 𝜋2 𝑛𝜋
−𝑘 𝑡
𝑢(𝑥, 𝑡) = 𝐴 𝑒 𝐿2 𝑠𝑒𝑛 ( 𝑥)
𝐿

Satisfaciendo la condición inicial.



𝑛 2 𝜋2 𝑛𝜋
−𝑘 𝑡
𝑢(𝑥, 𝑡) = ∑ 𝐴𝑛 𝑒 𝐿2 𝑠𝑒𝑛 ( 𝑥)
𝐿
𝑛=1

entonces,

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

Reconocemos que.

2 𝐿 𝑛𝜋
𝐴𝑛 = ∫ 𝑓(𝑥 ) 𝑠𝑒𝑛 ( 𝑥 ) 𝑑𝑥
𝐿 0 𝐿

Sustituyendo.
∞ 𝐿
2 𝑛𝜋 𝑛𝜋 𝑛 2 𝜋2
−𝑘 2 𝑡
𝑢(𝑥, 𝑡) = ∑ (∫ 𝑓 (𝑥 ) 𝑠𝑒𝑛 ( 𝑥 ) 𝑑𝑥 ) 𝑠𝑒𝑛 ( 𝑥 ) 𝑒 𝐿
𝐿 0 𝐿 𝐿
𝑛=1

También podría gustarte