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

Ver Teoria

Enunciado

Resolva o problema

∂ t 2 u   =   c 2 ∂ x 2 u , 0   ≤   x   ≤   L ,   t   >   0 ;

u ( 0 , t ) = 0 ,   u ( L , t ) = 0 ;

u ( x ,   0 )   =   f ( x ) ∂ t u ( x ,   0 )   =   g ( x ) .

(d)   f ( x )   =     s e n 3 ( π L x ) ,   g ( x ) = 0 .  

Passo 1

Eaew, belezinha?

Temos que resolver o problema encontrando a função de onda u ( x ,   t ) que obedece às condições:

∂ t 2 u   =   c 2 ∂ x 2 u , 0   ≤   x   ≤   L ,   t   >   0 ;

u ( 0 , t ) = 0 ,   u ( L , t ) = 0 ;

u ( x ,   0 )   =   f ( x ) ∂ t u ( x ,   0 )   =   g ( x ) .

Sendo,

f ( x )   = s e n 3 ( π L x ) ,     g ( x ) = 0 .

Vamos utilizar algumas relações já conhecidas que foram explicadas no Responde Aí na seção de Cálculo: Série de Fourier e EDP, Tópico 5: EDP - Equações da Onda. Também vamos utilizar alguns resultados de integrais já conhecidos e algumas relações trigonométricas. Vamos lá!

Passo 2

Primeiro vamos lembrar que para o problema:

∂ t 2 u   =   c 2 ∂ x 2 u , 0   ≤   x   ≤   L ,   t   >   0 ;

u ( 0 , t ) = 0 ,   u ( L , t ) = 0 ;

u ( x ,   0 )   =   f ( x ) ∂ t u ( x ,   0 ) �� =   g ( x ) .

Com g ( x ) = 0 , temos, pelo material do Responde Aí na seção de Cálculo: Série de Fourier e EDP, Tópico 5: EDP - Equações da Onda, que a solução é dada pela série de Fourier:

u ( x , t ) = Σ ∞ n = 1 [   a n   s e n (   n π L x )   c o s (   n π   c L t ) ]

com o coeficiente sendo obtido pela relação:

a n = 1 L ∫ - L L f ( x )   s e n (   n π L x )   d x   ,   n   = 1 ,   2 ,   .   .   .  

Então, para obtermos a solução temos apenas que calcular esses coeficientes.

Passo 3

Vamos manipular a função para ela ser soma de senos e cossenos, pois dessa forma no momento que calcularmos os coeficientes da série as integrais já conhecidas irão aparecer. Então, utilizando as identidades trigonométricas:

s e n 2   ( α )   =   1   -   c o s ( 2 α ) 2 ,  

s e n ( α )   c o s ( β )   =   s e n ( α   +   β )   +   s e n ( α   −   β )   2  

Temos então que:

f ( x )   =   s e n 3 ( x ) = s i n ( x ) s i n 2 ( x ) = s i n ( x )   1   -   c o s ( 2 x ) 2

f ( x ) =   s i n ( x )   -   s i n ( x ) c o s ( 2 x ) 2 =   s i n ( x ) 2 -   s i n ( 3 x ) - s i n ( x ) 4

f ( x ) =   s i n ( x ) 2 -   s i n ( 3 x ) 4 +   s i n ( x ) 4 =   3   s i n ( x ) 4 -   s i n ( 3 x ) 4

Passo 4

Temos que:

f ( x )   =   s e n 3 ( π L x ) =   3   s i n ( x ) 4 -   s i n ( 3 x ) 4

Substituindo em a n :

a n = 1 L ∫ - L L (   3   s i n ( x ) 4 -   s i n ( 3 x ) 4 )   s e n (   n π L x )   d x

a n = 1 4 L ∫ - L L ( 3   s e n ( π L x )   s e n (   n π L x )     −   s e n ( 3 π L x )   s e n (   n π L x )   ) d x

Agora temos que lembrar das seguintes relações:

∫ - L L s e n (   n π L x )     s e n   (   m π L x )   d x   =   0 ,   n   =   1 ,   2 ,   .   .   .   ,   m   =   1 ,   2 ,   .   .   .     n   ≠ m   ;

∫ - L L s e n 2 (   n π L x )     d x   =   L   .

Então temos que as integrais de a n só serão diferentes de zero para   n = 1 e   n = 3 , então teremos os coeficientes:

a 1 = 3 4 L ∫ - L L   s e n 2 ( π L x ) d x = 3 4 L L = 3 4 ;

a 4 = - 1 4 L ∫ - L L   s e n 2 ( 3 π L x ) d x = - 1 4 L L = - 1 4 .

Passo 5

Assim podemos escrever a solução:

u ( x , t ) = Σ ∞ n = 1 [   a n   s e n (   n π L x )   c o s (   n π   c L t ) ] = 3   s e n (   π L x )   c o s ( π   c L t ) - s e n (   4 π L x )   c o s (   4 π   c L t )

u ( x , t ) = 3 4   s e n (   π L x )   c o s ( π   c L t ) - 1 4 s e n (   3 π L x )   c o s (   3 π   c L t ) .

A questão foi isso! Te vejo na próxima! Tchau Tchau!

Resposta

u ( x , t ) = 3 4   s e n (   π L x )   c o s ( π   c L t ) - 1 4 s e n (   3 π L x )   c o s (   3 π   c L t ) .

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.