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

Ver Teoria

Enunciado

Usando separação de variáveis, encontrar a solução da seguinte equação do calor. Explique detalhadamente como se resolve o problema

u t = u x x ,     0 < x < 1 ,     t > 0 u x 0 , t = 0 ,     u x 1 , t = 0 ,     t ≥ 0 u x , 0 = f x ,     0 ≤ x ≤ 1

onde f é definida como a extensão par e periódica de período P = 2 da função f x =   1   onde     0 < x < 1 e esboce o gráfico no intervalo - 1 < x < 1 .

Passo 1

Para resolver o problema por separação de variáveis, vamos escrever

u x , t = X x T t

Note que se trata de uma equação de calor.

Com a separação, teremos:

u t x , t = X x T ' ( t )

u x x x , t = X ' ' x T ( t )

Portanto, nossa equação se torna

u t = u x x

X x T ' t = X ' ' x T ( t )

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

E isso origina duas EDO’s (uma em x e outra em t ):

T ' t + λ T t = 0 X ' ' x + λ X x = 0

Passo 2

Agora, vamos olhar para as condições de contorno:

u x 0 , t = 0 → X ' 0 T t = 0       ( I )

u x 1 , t = 0 → X ' 1 T t = 0       ( I I )

De ( I ) e ( I I ) , se T t = 0 , então ambas as condições seriam atendidas. Mas como isso implicaria em u x , t = X x T t = 0 o que daria a solução trivial.

Passo 3

Neste caso, nos resta admitir o seguinte:

X ' 0 = 0

X ' 1 = 0

Assim, formamos o seguinte problema:

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

Buscando uma solução não-trivial ( X x ≠ 0 ) , vamos avaliar a equação característica:

X ' ' x + λ X x = 0

r 2 + λ = 0

r = ± - λ

Temos um problema de contorno e vamos ter que analisar cada caso de Δ .

Passo 4

Caso 1: ∆ > 0 . Ou seja, o que está dentro da raiz é positivo. Nesse caso, λ < 0 . As soluções são do tipo:

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

Calculando X ' ( x ) .

X ' x = - λ c 1 e - λ x - - λ c 2 e - - λ x

Agora vamos usar as condições de contorno: X ' 0 = 0 X ' 1 = 0

X ' ( 0 ) = - λ c 1 - - λ c 2 = 0

X ' 1 = - λ c 1 e - λ - - λ c 2 e - - λ = 0

Isso nos dá

c 2 = c 1

c 1 e - λ - c 1 e - - λ = 0

c 1 ( e - λ - e - - λ ⏟ ≠ 0 ) = 0

c 1 = 0

c 2 = c 1 = 0

Substituindo esses valores na nossa solução geral, temos

X x = 0

Mas isso não nos ajuda em nada! Essa é a solução trivial, que não nos interessa. Portanto, esse caso, ∆ > 0 , não nos dá soluções para a EDO.

Passo 5

Caso 2: ∆ = 0 . Nesse caso, λ = 0 . As soluções são do tipo:

X x = c 1 + c 2 x

Calculando X ' ( x )

X ' x = c 2

Agora vamos usar as condições de contorno: X ' 0 = 0 X ' 1 = 0

X ' 0 = c 2 = 0

Isso nos deixa com a solução

X x = c 1

Vamos calcular T ( t ) para esse λ

T ' t = 0

T t = c 3

Nossa solução fica

u x , t = c 1 c 3

Temos que analisar agora a condição u x , 0 = f ( x ) e isso vai nos dar que.

c 1 c 3 = f x ⇒ u x , t = f ( x )

Passo 6

Caso 3: ∆ < 0 . Ou seja, o que está dentro da raiz é negativo. Nesse caso, λ > 0 . As soluções são do tipo:

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

Substituindo r = ± - λ = 0 ± i λ

X x = C 1 cos ⁡ λ x + C 2 sen ⁡ λ x

Calculando X ' ( x )

X ' x = - C 1 λ  sen ⁡ λ x + C 2 λ  cos ⁡ λ x

Aplicando as condições de contorno

X ' 0 = C 2 λ = 0 → C 2 = 0

X ' 1 = - C 1 λ  sen ⁡ λ = 0 →  sen ⁡ λ = 0

Note que na última equação não podemos ter C 1 = 0 (se não teríamos X x = 0 ). Então, precisamos de  sen ⁡ λ = 0 . Como C 2 = 0 , sobra

X x = C 1 cos ⁡ λ x

Passo 7

Como queremos  sen ⁡ λ = 0 , basta fazer λ = n π ,   n ∈ Z + , de modo que obtemos λ n = n π 2 . Estes são nossos autovalores!

Portanto, ficamos com

X x = C 1 cos ⁡ n π x

que são nada mais nada menos que nossas autofunções!

Passo 8

Agora que já conhecemos os autovalores λ n = n π 2 , vamos voltar pra equação em t :

T ' t + λ T t = 0

d T d t = - λ T

∫       d T T = - λ ∫     d t

ln ⁡ T = - λ t

T t = e - λ t

Fácil, né? Vamos substituir λ por λ n :

T t = e - n π 2 t

E, assim, uma solução deve ser

u x , t = X x T t

u x , t = C 1 cos ⁡ n π x e - n π 2 t

Passo 9

Agora vem o pulo do gato. Se qualquer C 1 cos ⁡ n π x é solução para X ( x ) , então uma combinação linear dessas soluções (para diferentes valores de n ) também é solução:

X x = d 0 2 + ∑ n = 1 ∞ d n   cos ⁡ n π x

Assim:

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

Sabendo que u x , 0 = f ( x )

f x = d 0 2 + ∑ n = 1 ∞ d n  cos ⁡ n π x

Para calcularmos os coeficientes da nossa função temos que tratar a expansão em série de Fourier da nossa função f ( x ) .

Passo 10

Vamos primeiro achar essa nossa função. Se a função f x = 1 ,   0 < x < 1 é par, então f x = 1 ,   - 1 < x < 0 . Bom, assumimos então que f - 1 = f 0 = f 1 = 0 . E como o período é P = 2 , vemos que f 1 = f - 1 + 2 = f ( - 1 ) . Um esboço é ilustrado abaixo.

Passo 11

Agora que nós já sabemos como é a nossa função podemos escrever ela como uma série de Fourier e achar os coeficientes que estamos procurando.

f x = a 0 2 + ∑ n = 1 ∞ a n cos ⁡ n π x L + b n s e n n π x L

Onde 2 L = P é o período; portanto, L = 1 .

Passo 12

Vamos avaliar os coeficientes da série:

a 0 = 1 L ∫ - L L f x d x

a 0 = ∫ - 1 1 f x d x

a 0 = 1 . 1 - - 1 = 2

Note que as descontinuidades de f ( x ) não afetam a integral.

Agora, temos:

a n = 1 L ∫ - L L f x cos ⁡ n π x L d x

a n = ∫ - 1 1 cos ⁡ n π x d x

a n = s e n n π x n π - 1 1

a n = 0

Além disso:

b n = 1 L ∫ - L L f x sen ⁡ n π x L d x

b n = ∫ - 1 1 sen ⁡ n π x d x

b n = - c o s n π x n π - 1 1

b n = 0

Passo 13

Combinando os resultados, a série de Fourier é simplesmente

f x = a 0 2

f x = 1

Passo 14

No nosso problema de equação de calor os coeficientes d n   são justamente os coeficientes da série de Fourier de f ( x ) : d n = a n .

Lembre que a 0 = 2 e a 1 = a 2 = … = 0 , então d 0 = 2 e d 1 = d 2 = … = 0 . Logo, sobra um termo apenas:

u x , t = d 0 2 = 1

Na realidade, a solução é exatamente u x , t = f ( x ) .

Resposta

u x , t = f ( 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