Exercícios Resolvidos de Equação da Onda

Ver Teoria

Enunciado

Determina as soluções de

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

u x , 0 = f x ,           ∂ t u x , 0 = g x ,

Nos casos seguintes:

b )   f x = e - x 2 ,   g x = 2 c x e - x 2 .

Passo 1

E aí gente, como vocês estão?

Nessa questão recebemos a EDP conhecida como equação da onda, e precisamos determinar as soluções para cada um dos problemas de valor inicial dados.

Temos a EDP:

∂ t 2 u = c 2 ∂ x 2 u

Com as condições iniciais:

u x , 0 = e - x 2

∂ t u x , 0 = 2 c x e - x 2

Pra resolver esse problema tem uma receitinha de bolo que precisamos seguir. Vamos dar uma revisada rápida:

1 ) Usar o método da separação de variáveis para obter as EDO’s separadas em X ( x ) e T ( t ) , onde:

u x , t = X x T ( t )

2 ) Resolver a EDO que tiver duas condições de contorno e substituir o valor de λ encontrado na outra EDO.

3 ) Fazer a combinação linear das soluções obtidas

4 ) Obter os coeficientes resolvendo as séries de Fourier das condições iniciais

É galera... dá um trabalho bom mesmo! Mas vamos dar o gás que acabamos rapidinho

Passo 2

Começando pelo método de separação de variáveis, façamos:

u x , t = X x T ( t )

Então:

∂ x 2 u = T x ⋅ X ' ' x

∂ t 2 u = X x ⋅ T ' ' ( t )

Substituindo na EDP:

X x ⋅ T ' ' t = c 2 ⋅ X ' ' x ⋅ T t

Logo:

T ' ' t c 2 T t = X ' ' x X x = λ

Daí, temos duas EDO’s:

X ' ' x - λ X x = 0 T ' ' x - λ c 2 T t = 0

Massa! Vamos analisar algumas condições de contorno do nosso problema.

Como se trata da equação da onda, vamos considerar que as extremidades estão fixas (embora não tenha sido dado no enunciado, essa informação é fundamental pra continuarmos resolvendo nosso probleminha).

Assim:

u 0 , t = u L , t = 0

Como separamos a função u , temos:

X 0 ⋅ T t = 0

Para fugir da solução trivial (isto é, T ( t ) constante e igual a 0 ), ficamos com:

X 0 = 0

Analogamente, teremos também:

X L = 0

Daí, podemos montar o seguinte PVI em x :

X ' ' x - λ X x = 0 X 0 = X L = 0

E:

T ' ' t - λ c 2 T t = 0

Vamos resolver primeiro o PVI, depois que descobrirmos o valor de λ substituímos na EDO em T para encontrar a sua solução.

Passo 3

X ' ' x - λ X x = 0 X 0 = X L = 0

A EDO nos dá a equação característica:

r 2 - λ = 0

∴ r = λ

Os casos λ ≥ 0 nos dará a solução trivial, que não estamos interessados. Daí, vamos considerar o caso λ < 0 . Ficamos com:

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

Como:

X 0 = 0

Temos:

X 0 = c 2 = 0

Como X L = 0 , temos:

X L = c 1 s e n - λ L = 0

Logo:

- λ L = n π

∴ λ n =   - n 2 π 2 L ,     n = 1,2 , 3 , …

Muito bom! Temos então a solução particular:

X x = s e n n π x L ,     n = 1,2 , 3 , …

Substituindo esse resultado na EDO em T ( t ) , vamos ter que:

T n t = a n cos ⁡ n π c t L + b n s e n n π c t L  

Excelente!

Aí agora, pra escrevermos a solução formal da equação da onda, basta considerarmos a combinação linear das soluções particulares que encontramos.

Ou seja:

u x , t = ∑ n = 1 ∞ a n cos ⁡ n π c t L + b n s e n n π c t L s e n n π x L ,       n = 1,2 , 3 , …

Muito bom!

Passo 4

Mas quem é a n e b n ?

Então, pra determinarmos esses coeficientes vamos finalmente usar as condições iniciais dadas no problema. Como:

u x , 0 = e - x 2

Substituindo na nossa solução formal, temos que:

u x , 0 = ∑ n = 1 ∞ a n s e n n π x L = e - x 2

Conseguiu perceber alguma coincidência com a série de Fourrier? No somatório estamos desenvolvendo justamente a extensão ímpar da função e - x 2 . Assim, o seu coeficiente pode ser calculado como:

a n = 2 L ∫ 0 L e - x 2 s e n n π x L d x

Perfeito!

Substituindo agora a outra condição de contorno, vamos ter que:

∂ t u x , 0 = 2 c x e - x 2

Derivando a solução formal da EDP com relação a t , temos:

∂ t u = ∑ n = 1 ∞ - a n n π c L ⋅ s e n n π c t L + b n n π c L ⋅ cos ⁡ n π c t L s e n n π x L

Fazendo t = 0 , ficamos com:

∂ t u x , 0 = ∑ n = 1 ∞ b n n π c L s e n n π x L = 2 c x e - x 2

Novamente, percebendo a semelhança com a série de Fourrier, vamos ter que o coeficiente que multiplica o seno é dado por:

b n n π c L = 2 L ∫ 0 L 2 c x e - x 2 s e n n π x L d x

Logo:

b n = 2 n π c ∫ 0 L 2 c x e - x 2 s e n n π x L d x

b n = 4 n π ∫ 0 L x e - x 2 s e n n π x L d x

Excelente!

Essas integrais que definem a n e b n são calculadas computacionalmente, por isso vamos parar por aqui, beleza?

Espero que vocês tenham curtido a solução e até a próxima 😊

Resposta

u x , t = ∑ n = 1 ∞ a n cos ⁡ n π c t L + b n s e n n π c t L s e n n π x L ,       n = 1,2 , 3 , …

a n = 2 L ∫ 0 L e - x 2 s e n n π x L d x

b n = 4 n π ∫ 0 L x e - x 2 s e n n π x L d x

Exercícios de Livros Relacionados

Os Problemas 19 e 20 indicam como é possível mostrar que a solução formal (20), (22) das Eqs. (1), (3) e (9) é, de fato, a solução desse problema.19- Usando a identidade trigonométrica s e nAc o sB=½
Seja h ( ξ ) o deslocamento inicial em [ 0 ,L ] estendido a ( - L ,0 ) como uma função ímpar e estendido ao resto da reta como uma função periódica de período 2 L . Supondo a continuidade de h e suas
Combinando os resultados dos Problemas 16 e 17, mostre que a solução do problemaé dada por
O movimento de uma membrana circular elástica, como a membrana de um tambor, é determinado pela equação de onda bidimensional em coordenadas polaresSupondo que u ( r ,θ ,t )=R r Θ ( θ ) T ( t ) , enco
Considere a equação de onda a 2 u x x = u t t em um meio unidimensional infinito, sujeita às condições iniciaisUsando a forma da solução obtida no Problema 13, mostre queUse a primeira das equações no