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

Ver Teoria

Enunciado

Questão 4: Utilize o método de Separação de Variáveis para encontrar a solução do problema

u t - u x x = 0                                                                                       e m               0 ,   1 × 0 ,   + ∞   u 0 ,   t = u 1 ,   t = 0                                                         p a r a       t ∈ 0 ,   + ∞                   u x ,   0 = s e n ( 2 π x ) - s e n ( 5 π x )                 p a r a       x ∈ 0 ,   1                          

Passo 1

Primeiro vamos estabelecer aqui o que é o nosso objetivo principal. Esse objetivo é no que a gente quer chegar, saca? Ora, se é o problema nos pediu uma solução de uma EDP, o que ele quer encontrar é a função u x ,   t , e esse é o nosso objetivo!

Como vamos fazer isso? O problema já nos ordenou ali em cima: por separação de variáveis.

Só uma curiosidade mesmo, essa EDP é a famosa equação do calor, pois a derivada com relação t é uma derivada de primeira ordem . Agora vamos nessa?

Passo 2

O método de separação de variáveis, para você que tem muitas matérias para colocar em dia e não se recorda, é o método que você chuta a solução, como sendo um produto de soluções de variáveis isoladas. Tipo, no nosso caso, a solução é u x ,   t , e vamos falar que é a mesma coisa que um produto de uma função só de x por outra só de y :

u x ,   t = X x ⋅ T t

Saca? Agora que vem a graça do problema, se essa é uma solução então ela satisfaz a equação da EDP...

u t - u x x = 0           →             X x ⋅ T t t - X x ⋅ T t x x = 0         →

O legal disso é que, por exemplo, em X x ⋅ T t t , a derivada é com relação a t , como X x , independe de t , então...

X x ⋅ T t t - T t ⋅ X x x x = 0         →

Nas duas parcelas, como estamos derivando em relação a uma única variável, essa é uma derivada comum, entende? Tipo f x x = f ' x ...

X x ⋅ T ' t - T t ⋅ X ' ' x = 0       →

Repare que como temos coisas independentes, podemos separar esses caras, ou seja, jogar tudo que tem x por um lado e tudo que tem y por outro. Veja:

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

T ' t T t = X ' ' x X x

Quase pronto! As variáveis estão até separadas, mas falta mais uma coisinha: como são possíveis duas funções de duas variáveis diferentes serem iguais? A única solução é que elas sejam ambas iguais uma constante, que não depende nem de uma e nem de outra, saca?

T ' t T t = X ' ' x X x = - λ

O menos na frente da constante serve apenas para simplificar um pouco as continhas, mas você pode chamar simplesmente de λ , ou a , ou qualquer outro nome de constante que tiver em mente . Assim...

T ' t = - λ ⋅ T t           →                 T ' t + λ T t   = 0

e

X ' ' x = - λ ⋅ X x           →                 X ' ' x + λ ⋅ X x = 0

Passo 3

Separamos as variáveis na EDP, mas no nosso problema, existem outras condições que envolvem a nossa u ...

        u 0 ,   t = u 1 ,   t = 0                     u x ,   0 = sin ⁡ 2 π x - sin ⁡ 5 π x                          

Olhando para as condições de contorno:

u 0 ,   t = u 1 ,   t = 0                 →                   u 0 ,   t = 0 u 1 ,   t = 0

Substituindo a nossa solução...

X 0 ⋅ T t = 0                           e                           X 1 ⋅ T t = 0

Nos dois casos, vamos considerar que T t ≠ 0 , pois se não a gente recai na solução trivial da EDP, ou seja...

u x ,   t = X x ⋅ T t = 0                     ( e   n ã o   é   e s s a   s o l u ç ã o   q u e   a     g e n t e   q u e r )

Então, pela condição de contorno...

X 0 = 0

X 1 = 0

Já pela condição de temo inicial...

u x ,   0 = sin ⁡ 2 π x - sin ⁡ 5 π x         →

f x = sin ⁡ 2 π x - sin ⁡ 5 π x

Assim, podemos separar agora a nossa EDP, em duas EDO’s:

u t - u x x = 0       u 0 ,   t = u 1 ,   t = 0                     u x ,   0 = sin ⁡ 2 π x - sin ⁡ 5 π x                                 →                 X ' ' x + λ ⋅ X x = 0 X 0 = 0 X 1 = 0           e

T ' t + λ T t   = 0 u x ,   0 = f x

Veja que a primeiro podemos resolver a primeira EDO, já que dela temos muitas informações de condição de contorno ! Vamos nessa?

Passo 4

X ' ' x + λ ⋅ X x = 0 X 0 = 0 X 1 = 0

A gente pode resolver essa EDO, por equação característica, afinal λ é uma constante...

X ' ' x + λ ⋅ X x = 0             →               E q u a ç ã o   C a r a c t .           →               b 2 + λ = 0       →

b 2 = - λ         →             b = - λ o u b = - - λ

Só que veja que aqui temos que ter cuidado! Não conhecemos lambda, então temos que supor cada caso, se lembra?

  1. Para Lambda menor que 0 ( λ < 0 ):
  2. A solução será real e vai ser do tipo...

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

    Aplicando as condições de contorno...

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

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

    Logo, isolando em cima e substituindo em baixo, ganhamos...

    - c 2 e - λ + c 2 e - - λ = 0     →

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

    Nesse caso, como a exponencial nunca é zero, então...

    c 2 = 0               o   q u e   i m p l i c a       c 1 = 0             →

    X x = c 1 e - λ ⋅ x + c 2 e - - λ   x = 0

    Essa é a solução trivial que a gente não quer, logo λ não pode ser menor que 0 .

  3. Para λ = 0 ...
  4. X ' ' x + λ ⋅ X x = 0               →               X ' ' x = 0                   →

    X ' x = c 1         →

    X x = c 2 + c 1 x

    Aplicando as condições de contorno...

    X 0 = 0         →               X 0 = c 2 = 0                  

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

    Novamente temos a solução trivial!

  5. Já para λ > 0 ...

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

O que está dentro da raiz é negativo, então a gente tem uma solução aparentemente imaginária! Mas se você se recordar um pouquinho de EDO’s, vai ver que essa resposta é real dessa forma aqui:

X x = c 1 e α x sin ⁡ β x + c 2 e α x cos ⁡ β x

Onde

b = α + β i   →         b = 0 + λ i    

Logo, a solução será do tipo...

X x = c 1 sin ⁡ λ x + c 2 cos ⁡ λ x

Aplicando as condições de contorno...

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

X 1 = 0         →               X 1 = c 1 sin ⁡ λ = 0

Se c 1 = 0 , temos a solução trivial de novo, mas dessa vez o c 1 pode ser diferente de 0 , tudo que precisa acontecer para isso é...

sin ⁡ λ = 0           →             λ = n π                 p a r a   n = 1,2 …

Logo...

λ n = n 2 π 2

Esse é o autovalor da nossa E D O ! E a autofunção é

X x = sin ⁡ n π x

Encontramos uma solução que não é trivial! . Segura ela aí, então e vamos embora!

Passo 5

Bem, agora que temos um conjunto de valores de λ que não dá uma solução trivial, a gente pode seguir com a solução indo agora para a segunda EDO:

T ' t + λ T t   = 0                                   →

T ' t + n 2 π 2 T t   = 0                                   →

Resolvendo por separação de variáveis, essa EDO...

d T d t = - n 2 π 2 T           →               d T T = - n 2 π 2 d t       →

∫ d T T = ∫ - n 2 π 2 d t         →

ln ⁡ T = - n 2 π 2 t + c 3       →

T = e - n 2 π 2 t + c 3       →

Lembrando que e - n 2 π 2 t + c 3 = e - n 2 π 2 t ⋅ e c 3 , e chamando e c 3 de constante 3 (é só um nome)...

T t = c 3 e - n 2 π 2 t

Passo 6

Veja que agora temos um conjunto de soluções desse tipo aqui:

u x ,   t = X x ⋅ T t           ⇒

u n x ,   t = c 1 sin ⁡ n π x ⋅ c 3 e - n 2 π 2 t               p a r a       n = 1 ,   2 …

Considerando c 1 ⋅ c 2 = b n

u n x ,   t = b n sin ⁡ n π x ⋅ e - n 2 π 2 t               p a r a       n = 1 ,   2 …

Logo, a nossa solução geral vai ser a soma de todas essas soluções e múltiplos b n delas...

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

Quem é b n ? Bem, para isso podemos olhar a última condição que não utilizamos para nada: a condição inicial!

u x ,   0 = sin ⁡ 2 π x - sin ⁡ 5 π x = f x

Então venha comigo mais um pouco (tá quase acabando! ).

Passo 7

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

u x ,   0 = ∑ n = 1 ∞ b n sin ⁡ n π x = f x     →

∑ n = 1 ∞ b n sin ⁡ n π x = sin ⁡ 2 π x - sin ⁡ 5 π x

Logo, desenvolvendo a série um pouco até n = 5 ...

b 1 sin ⁡ π x + b 2 sin ⁡ 2 π x + b 3 sin ⁡ 3 π x + b 4 sin ⁡ 4 π x + b 5 sin ⁡ 5 π x + ∑ n = 6 ∞ b n sin ⁡ n π x = sin ⁡ 2 π x - sin ⁡ 5 π x

Logo, b 2 = 1 e b 5 = - 1 , todos os demais tem que zerar! Assim, a solução fica sendo...

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

u x ,   t = b 2 sin ⁡ 2 π x e - 2 2 π 2 t + b 5 sin ⁡ n π x e - 5 2 π 2 t       →

u x ,   t = sin ⁡ 2 π x e - 4 π 2 t - sin ⁡ 5 π x e - 25 π 2 t

E essa é a solução final!!!

Resposta

u x ,   t = sin ⁡ 2 π x e - 4 π 2 t - sin ⁡ 5 π x e - 25 π 2 t

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