Exercícios Resolvidos de Modelagem com Equação do Calor

Ver Teoria

Enunciado

Resolva o problema

E D     ∂ t u = 3 ∂ x 2 u ,       0 ≤ x ≤ 3 ,   t ≥ 0 ;

                    C F     ∂ x u 0 , t = 0 ,       ∂ x u 3 , t = 0 ,       t ≥ 0 ;                     C I   u x , 0 = f x ,         0 ≤ x ≤ 3

no caso seguinte:

c )   f x =   - s e n 2 2 π 3 x

Passo 1

Olá pessoal, como vocês estão?

Temos aqui um probleminha das nossas famosas EDP: equação diferencial parcial, mais especificamente, a equação do calor.

Como foi dado nas condições de fronteira que:

∂ x u 0 , t = ∂ x u 3 , t = 0 ,   t ≥ 0

trata-se de um problema em que as extremidades estão isoladas termicamente. Pois, a taxa de variação da temperatura em x nas extremidades da barra é nula.

Também foi dado a condição inicial:

u x , 0 = - s e n 2 2 π 3 x  

Massa! Sabendo de tudo isso, vamos usar o método de separação das variáveis pra resolver essa EDP, ou seja, usaremos que:

u x , t = f x g ( t )

Vamo nessa!

Passo 2

Pra resolver esse tipo de Equação, usamos o método da separação de variáveis. Daí, fazendo:

u x , t = f x g ( t )

E, substituindo na EDP, vem:

f x g ' t - 3 g t f ' ' x = 0

∴ f ' ' x f x = g ' t 3 g t

A única forma da razão dessas funções serem sempre iguais, independentemente do valor de x e de t , é se tivermos:

f ' ' x f x = g ' t 3 g t = λ

onde λ é uma constante.

Assim, reduzimos o nosso problema a duas EDO’s:

f ' ' x - λ f x = 0 f ' 0 = f ' 3 = 0

E:

g ' t - 3 λ g t = 0

Beleza!

Daí precisamos analisar os possíveis valores de λ :

- Se λ > 0 , então:

r = ± λ

f x = c 1 e - λ x + c 2 e λ x

Testando as condições iniciais:

f ' x = c 1 - λ ⋅ e - λ x + c 2 λ ⋅ e λ x

f ' 0 = - c 1 + c 2 λ = 0

Como não queremos a solução trivial:

c 2 - c 1 = 0

∴ c 2 = c 1

Logo:

f ' 3 = c 1 - λ ⋅ e - λ 3 + c 1 λ ⋅ e λ 3 = 0

Daí, temos:

c 1 = c 2 = 0

Ou seja, nos deu a solução trivial, que não estamos interessados...

- Se λ = 0 , então:

r = 0

∴ f x = c 1 + c 2 x

Colocando as condições iniciais:

f ' 0 = c 2 = 0

E:

f ' 3 = c 2 = 0

Daí, temos uma solução:

f x = c 1

Massa!

- Se λ < 0 , então:

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

Como:

f ' x = - - λ ⋅ c 1 ⋅ sen ⁡ - λ x + - λ ⋅ c 2 c o s ⁡ ( - λ x )

Pela condição inicial:

f ' ( 0 ) = c 2 - λ = 0

∴ c 2 = 0

Logo:

f ' ( 3 ) = - λ c 1 s e n 3 - λ = 0

Como não queremos a solução trivial, devemos ter:

s e n 3 - λ = 0

∴ 3 - λ = n π

⇒ λ = - n 2 π 2 3 2

Daí, desconsiderando a constante c 1   ( c 1 = 1 ) , temos que as autofunções são:

f n x = c o s n π x 3 ,     n = 1,2 , 3 , …

Perfeito.

Passo 3

Substituindo agora o valor que encontramos para λ na EDO em g ( t ) , teremos que:

g ' t - 3 - n 2 π 2 9 g t = 0

Logo:

g t = e - 3 n 2 π 2 9 t

= e - n 2 π 2 3 t

Daí, juntando todas as soluções de f e g , e lembrando que:

u x , t = f x g t

Chegamos que:

u x , t = c 0 2 + ∑ n = 1 ∞ c n e - n 2 π 2 3 t ⋅ cos ⁡ n π x 3

Pra matar a questão precisamos apenas descobrir quanto vale c n , e é justamente aí que entra nossa condição inicial, C I , dada no problema.

u x , 0 = c 0 2 + ∑ n = 1 ∞ c n ⋅ cos ⁡ n π x 3 = - s e n 2 2 π 3 x

Ou seja gente, a nossa solução da EDP para t = 0 é justamente a extensão par e periódica de f ( x ) .

Assim, o coeficiente c n , para n = 0,1 , 2 , … são os coeficientes da série de Fourier de f ( x ) , que podem ser calculados como:

c n = 2 3 ⋅ ∫ 0 3 - s e n 2 2 π 3 x ⋅ cos ⁡ n π x 3 d x

Massa! Nossa missão então se resume a resolver essa integral!

Passo 4

Para calcular c n , precisamos resolver a integral:

∫ 0 3 s e n 2 2 π 3 x ⋅ cos ⁡ n π x 3 d x

Vamos usar a seguinte identidade trigonométrica:

s e n 2 x = 1 - cos ⁡ 2 x 2

Assim, vamos ter:

∫ 0 3 - s e n 2 2 π 3 x ⋅ cos ⁡ n π x 3 d x = 1 2 ∫ 0 2 cos ⁡ n π x 3 d x - 1 2 ∫ 0 3 cos ⁡ 4 π 3 x ⋅ cos ⁡ n π x 3 d x

Vamos começar resolvendo a primeira integral:

∫ 0 3 cos ⁡ n π x 3 d x = - 3 n π ⋅ s e n n π x 3 0 3

=   - 3 n π ⋅ s e n n π = 0

Esse resultado é válido para qualquer n ≠ 0 . Mas, se n = 0 , temos:

∫ 0 3 cos ⁡ 0 ⋅ π x 3 d x = ∫ 0 3 d x = 3

Conclusão:

∫ 0 2 cos ⁡ n π x 3 d x = 3 ,     n = 0 0 ,     n ≠ 0

Indo agora para a segunda integral, precisamos calcular:

∫ 0 3 cos ⁡ 4 π 3 x ⋅ cos ⁡ n π x 3 d x

Vamos usar a seguinte transformação trigonométrica:

cos ⁡ p + cos ⁡ q = 2 ⋅ cos ⁡ p + q 2 ⋅ cos ⁡ p - q 2

Daí, vem que:

cos ⁡ 2 x ⋅ cos ⁡ n π x 3 = 1 2 cos ⁡ n + 4 π x 3 - cos ⁡ n - 4 π x 3

Substituindo na integral:

∫ 0 3 cos ⁡ 4 π 3 x ⋅ cos ⁡ n π x 3 d x = 1 2 ⋅ ∫ 0 3 cos ⁡ n + 4 π x 3 d x + 1 2 ∫ 0 3 cos ⁡ n - 4 π x 3 d x

Se n ≠ 4 , temos como resultado das integrais:

= 1 2 ⋅ - 3 n + 4 π ⋅ s e n n + 4 π x 3 0 3 + 1 2 ⋅ - 3 n - 4 π ⋅ s e n n - 4 π x 3 0 3

- 3 2 n + 5 π ⋅ s e n ( n + 3 π ) - 3 2 n - 5 π ⋅ s e n ( n - 3 π )  

O problema é que, como n ∈ N , n ± 4 sempre será um número inteiro. Logo:

s e n n + 4 π = s e n n - 4 π = 0

Perfeito! Mas, se n = 4 , vamos ter:

∫ 0 3 cos ⁡ 4 - 4 π x 3 d x = ∫ 0 3 1 ⋅ d x

= 3

Portanto:

∫ 0 3 cos ⁡ 4 π x 3 ⋅ cos ⁡ n π x 3 d x = 3 ,     n = 4 0 ,     n ≠ 4

Passo 5

Juntando os resultados que encontramos em c n , temos:

c n = 2 3 ⋅ ∫ 0 3 - s e n 2 2 π 3 x ⋅ cos ⁡ n π x 3 d x

Então:

c n =   - 2 3 ⋅ 3 ,     n = 0 - 2 3 ⋅ 3 ,     n = 4 0 ,     n ≠ { 0,4 }

Assim, temos que:

u x , t = c 0 2 + ∑ n = 1 ∞ c n e - n 2 π 2 3 t ⋅ cos ⁡ n π x 3

Vale:

u x , t =   - 1 - 2 ⋅ e - 16 π 2 3 t ⋅ cos ⁡ 4 π x 3

E é isso galera! Espero que vocês tenham gostado.

Valeeeu 😊

Resposta

u x , t =   - 1 - 2 ⋅ e - 16 π 2 3 t ⋅ cos ⁡ 4 π x 3

Exercícios de Livros Relacionados

Encontre a temperatura no estado estacionário u(ρ, φ) na esfera de raio unitário se a temperatura for independente de θ e satisfizer a condição de contorno:
Encontre a temperatura no estado estacionário em uma barra semi-infinita 0
Considere o fluxo de calor em um cilindro de comprimento infinito e raio 1: 0 ≤ r e as seguintes condições inicial e de contorno:Mostre que em que J0(λn) = 0. Encontre uma fórmula para cn.