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

Ver Teoria

Enunciado

Questão 2: Resolva o problema de condução de calor em uma barra metálica:

25 u x x = u t ,                 0 < x < 5 ,   t > 0 u 0 ,   t = 0 ,           u 5 ,   t = 0 ,   t > 0 u x ,   0 = 2 - x ,     0 < x < 5

Passo 1

Então, talvez a gente consiga resolvê-la por separação:

u x ,   t = X x ⋅ T t

Se isso é solução, então...

u x x = X ' ' x ⋅ T t

u t = X x ⋅ T ' t

Logo...

25 u x x = u t   →

25 X ' ' x ⋅ T t = X x ⋅ T ' t →

25 X ' ' x X x = T ' t T t →

Essas duas funções diferentes variáveis só é igual se for uma constante λ :

25 X ' ' x X x = T ' t T t = λ

Então...

25 X ' ' x X x = λ

E

T ' t T t = λ

Portanto, ganhamos duas EDOs:

X ' ' x = λ 25 X x   T ' t = λ T t      

Passo 2

Antes de resolvermos essas EDO’s temos que encontrar as condições de contorno e condições iniciais para sabermos por qual EDO começarmos. Vamos nessa!

Sabemos que pelas condições de contorno...

u 0 ,   t = 0 ,           u 5 ,   t = 0

Como

u x ,   t = X x ⋅ T t

Então...

u 0 ,   t = 0   →       X 0 T t = 0   →       l o g o …

X 0 = 0

E

u 5 ,   t = 0   →       X 5 T t = 0   →       l o g o …

X 5 = 0

Passo 3

Beleza, então veja que a EDO que temos mais informações é a de X ...

X ' ' x = λ 25 X x   X 0 = 0 X 5 = 0      

Resolvendo essa EDO por Equação característica...

r 2 = λ 25       →           r = ± λ 25 = ± 1 5 λ

Então se λ > 0 ...

X = c 1 e 1 5 λ x + c 2 e - 1 5 λ x

Substituindo as condições de contorno:

X 0 = c 1 e 0 + c 2 e 0 = c 1 + c 2 = 0   →

c 1 = - c 2

E...

X 5 = - c 2 e λ + c 2 e - λ = 0

Como as exponenciais não são nulas...

c 2 = - c 1 = 0

Então essa solução é ruim! Vamos ver considerando λ = 0

X = c 1 + c 2 x   →

X 0 = 0     →         c 2 = 0

X 5 = 0     →           c 1 = 0

Também não queremos a solução constante!

Logo, a única opção que sobrou pra gente é...

λ < 0     →

X x = e 0 x ⋅ c 1 cos ⁡ 1 5 λ x + c 2 sin ⁡ 1 5 λ x =

c 1 cos ⁡ 1 5 λ x + c 2 sin ⁡ 1 5 λ x

Aplicando as condições agora...

X 0 = 0         →             X 0 = c 1 cos ⁡ 0 + c 2 sin ⁡ 0 =

c 1 = 0

E

X 5 = 0         →             X 5 = c 1 cos ⁡ λ + c 2 sin ⁡ λ   →

X 5 = c 2 sin ⁡ λ

Considerando c 2 diferente de 0 , para o seno ser 0 , esse termo dele tem que ser múltiplo de π ...

λ = m π   →

λ = m 2 π 2

Como λ tem que ser negativo...

λ = - m 2 π 2

Esses são os autovalores, e as autofunções são:

X x = sin ⁡ 1 5 m π x

Passo 4

Com esses autovalores podemos procurar a solução dessa outra EDO aqui:

T ' t = λ T t   →

T m ' t = - m 2 π 2   T m t

Perceba que a cada m , temos uma EDO diferente! Por equação característica podemos resolvê-la...

r = - m 2 π 2

Logo, a solução é...

T m = b m e - m 2 π 2 t

Onde b m é uma constante de solução para cada m .

Passo 5

A solução para cada m então é:

u m x ,   t = X x ⋅ T m t =

u m x ,   t = b m e - m 2 π 2 t sin ⁡ 1 5 m π x

Assim, o resultado do problema homogêneo é a soma de todas as soluções lineares:

u x ,   t = ∑ m = 1 ∞ b m e - m 2 π 2 t sin ⁡ 1 5 m π x

Passo 6

Achou que tinha acabado, jovem gafanhoto? Ainda falta saber quem é b m pela outra condição que não utilizamos:

u x ,   0 = 2 - x ,     0 < x < 5

Então...

u x ,   0 = ∑ m = 1 ∞ b m e 0 sin ⁡ 1 5 m π x →

u x ,   0 = ∑ m = 1 ∞ b m sin ⁡ 1 5 m π x = 2 - x

Perceba, que esse b m são os índices de uma série de Fourier em senos, e ela pode ser encontrada de acordo com Fourier para série apenas de senos:

b m = 2 5 ∫ 0 5 2 - x sin ⁡ 1 5 m π x   d x →

Integrando por partes...

u = 2 - x                     e                   d v = sin ⁡ 1 5 m π x d x

Logo...

d u = - d x                     e                   v = - cos ⁡ 1 5 m π x 1 5 m π

Como...

∫ u d v = u v - ∫ v d u   →

2 5 ∫ 0 5 2 - x sin ⁡ 1 5 m π x   d x = 2 5 - 2 - x ⋅ cos ⁡ 1 5 m π x 1 5 m π - ∫ 0 5 cos ⁡ 1 5 m π x 1 5 m π   d x →

2 5 ∫ 0 5 2 - x sin ⁡ 1 5 m π x   d x = 2 5 - 5 2 - x ⋅ cos ⁡ 1 5 m π x m π - 5 m π ∫ 0 5 cos ⁡ 1 5 m π x   d x →

2 5 ∫ 0 5 2 - x sin ⁡ 1 5 m π x   d x = 2 5 - 5 2 - x ⋅ cos ⁡ 1 5 m π x m π - 25 m 2 π 2 sin ⁡ 1 5 m π x 0 5 →

2 5 ∫ 0 5 2 - x sin ⁡ 1 5 m π x   d x = 2 5 - 5 2 - 5 ⋅ cos ⁡ m π m π - 25 m 2 π 2 sin ⁡ m π - - 5 2 ⋅ cos ⁡ 0 m π - 25 m 2 π 2 sin ⁡ 0 →

2 5 ∫ 0 5 2 - x sin ⁡ 1 5 m π x   d x = 2 5 - 5 2 - 5 ⋅ cos ⁡ m π m π - 25 m 2 π 2 sin ⁡ m π - - 5 2 ⋅ cos ⁡ 0 m π - 25 m 2 π 2 sin ⁡ 0 →

2 5 ∫ 0 5 2 - x sin ⁡ 1 5 m π x   d x = 2 5 15 ⋅ - 1 m + 10 m π →

2 5 ∫ 0 5 2 - x sin ⁡ 1 5 m π x   d x = 6 ⋅ - 1 m + 10 m π →

SOOO! Finished at least!

u x ,   t = ∑ m = 1 ∞ b m e - m 2 π 2 t sin ⁡ 1 5 m π x →

u x ,   t = ∑ m = 1 ∞ 6 ⋅ - 1 m + 10 m π e - m 2 π 2 t sin ⁡ 1 5 m π x

Resposta

u x ,   t = ∑ m = 1 ∞ 6 ⋅ - 1 m + 10 m π e - m 2 π 2 t sin ⁡ 1 5 m π x

Exercícios de Livros Relacionados

Em cada um dos problemas de 1 a 6, determine se o método de separação de variáveis pode ser usado para substituir a equação diferencial parcial dada por um par de equações diferenciais ordinárias. Se
Ao resolver equações diferenciais, quase sempre os cálculos podem ser simplificados por meio da utilização de variáveis adimensionais.Mostre que, se introduzirmos a variável adimensional ζ = x L , equ
Considere a condução de calor em uma barra com 40 cm de comprimento cujas extremidades sãomantidas à temperatura 0°C para todo t0 . Em cada um dos problemas de 9 a 12, encontre umaexpressão para a tem
Considere a condução de calor em uma barra com 40 cm de comprimento cujas extremidades sãomantidas à temperatura 0°C para todo t0 . Em cada um dos problemas de 9 a 12, encontre umaexpressão para a tem
Encontre a solução do problema de condução de calor 100 u x x = u t ,0x1 ,t0 ; u ( 0 , t ) = 0 ,u ( 1 , t ) = 0 ,t0 ; u ( x , 0 ) = s e n 2 π x - s e n 5 π x ,0 ≤ x ≤ 1