Resolva o seguinte problema de valor de contorno usando o método de separação de variáveis justificando detalhadamente TODA a análise

9 u x x   =   u t ,     0   ≤   x   <   4 ,   t > 0 u 0 ,   t = u 4 ,   t =   0 ,     t > 0 u x ,   0 =   5 s e n 3 π x 4 + 8 s e n 2 π x ,     0   ≤   x   ≤   4

Passo 1

Buscamos soluções por separação de variáveis :

u x ,   t = X ( x ) T ( t )

Substituindo   u na equação

9 X " x T ( t ) = X x T ' ( t )

Passando tudo que tem x para um lado e tudo que tem t para o outro lado

X " ( x ) X ( x )   = T ' ( t ) 9 T ( t )

Como essa igualdade vale para todo x e todo y , vamos igualá-la a uma constante

X " ( x ) X ( x )   = T ' ( t ) 9 T ( t ) =   - λ  

Separando a igualdade em duas equações, você vai encontrar isso aqui...

X " ( x ) + λ X ( x ) =  0

T ' t + λ T ( t ) = 0

Substituindo as condições de contorno...

X ( 0 ) T   ( t )   =   0   ⇒   X   ( 0 )   =   0  

X ( 4 ) T   ( t )   =   0   ⇒   X   ( 4 )   =   0  

X x T   0 =     5 s e n 3 π x 4 + 8 s e n 2 π x ⇒   T   ( 0 )   =   5 s e n 3 π x 4 + 8 s e n 2 π x

Passo 2

Isso nos dá um PVC em   X x .

Esqueça, por enquanto, a condição que tem uma função: u x ,   0 = 5 s e n 3 π x 4 + 8 s e n 2 π x .

Precisamos agora achar X x e   T   ( t ) u x ,   t = X ( x )   T   ( t )

Passo 3

Como temos um problema de autovalor, vamos montar e resolver a equação característica da equação de   X x

r 2 + λ = 0

r = ±   - λ  

Agora vamos analisar cada caso.

Se ∆ > 0   e   λ < 0 ,   teremos duas raízes r 1 = - λ   e     r 2 = - - λ    

X ( x ) = c 1 e x r 1 + c 2 e x r 2 = c 1 e x - λ   + c 2 e x - - λ

Nas condições de contorno...

X 0 = c 1 + c 2 =   0 ∴   c 1 = - c 2

X ( 4 ) = c 1 e 4 λ   - e - 4 λ =   0 ∴ c 1 = - c 2 = 0

Neste caso, só teremos soluções triviais.

Se ∆ = 0   e   λ = 0 ,   teremos duas raízes iguais a r = λ a solução da seguinte forma:

Que no nosso caso vai ficar...

X x = c 1 + c 2 x

Aplicando as condições iniciais

X 0 = c 1 + c 2 x = c 1 = 0

X 4 = c 1 + 4 c 2 = 4 c 2 = 0   ∴   c 2 = 0  

Também só há soluções triviais aqui!

Se ∆ < 0   e   λ > 0 , r = ± i   λ   e a solução tem essa cara aqui

X x = c 1 s e n λ   x + c 2 cos ⁡ λ   x

X 0 = c 2 = 0

X 4 = c 1 s e n 4 λ   = 0

Como queremos soluções não triviais

λ = n π 4 ,   n = 1,2 , 3 …

λ n = n ² π ² 16 ,     n = 1,2 , 3 …

X n = s e n ( n π 4 ) x ,     n = 1,2 , 3 …

Passo 4

Agora vamos trabalhar nossa função T ( t ) usando a condição inicial

T ' t + 9 λ T t = T ' t + 9 n ² π ² 16 T t = 0

Resolvendo a EDO de primeira ordem vamos ter

d T ( t ) T ( t ) = - 9 n ² π ² 16 d t

∫ d T ( t ) T ( t ) = l n T ( t )  

- 9 n ² π ² 16 ∫ d t = - 9 n ² π ² 16 t

l n T t = - 9 n ² π ² 16 t + C  

T t =   A e - 9 n ² π ² 16 t

Agora sim,

u n x ,   t = X x T t = s e n ( n π x 4 ) e - 9 n ² π ² 16 t

    n = 1,2 , 3 …

Passo 5

Pelo princípio da superposição, podemos escrever a série

u n x ,   y = ∑ n = 1 ∞ A n s e n ( n π x 4 ) e - 9 n ² π ² 16 t

Vamos aplicar a condição: u x ,   0 =   f t = 5 s e n 3 π x 4 + 8 s e n 2 π x

w n x ,   0 = ∑ n = 1 ∞ A n s e n ( n π x 4 ) = 5 s e n 3 π x 4 + 8 s e n 2 π x

De cara, percebemos que para n = 3 , A 3 = 5 ,   para n = 8 , A 8 = 8 . Para n ≠ 3 ,   8 ,   A n = 0   !

Assim, a série vai ser finita e sua soma:

u x ,   y = ∑ n = 1 ∞ A n s e n n π x 4 e - 9 n ² π ² 16 t = 5 e - 9 π ² 16 t s e n 3 π x 4 + 8 e - 64 π 2 16 t s e n 2 π x

Resposta

u n x ,   t = 5 e - 9 π ² 16 t s e n 3 π x 4 + 8 e - 64 π 2 16 t s e n 2 π x

Exercícios Resolvidos de Equação do Calor – Homogênea

Ver Teoria

Exercícios de Livros Relacionados