Equação do Calor – Homogênea
Teoria
Condução de calor
Antes de entrar no Cálculo, vamos entender um pouquinho da Física por trás desse assunto, para você se situar um pouco, ok? O nosso objeto de estudo vai ser a seguinte barra:
QUE, uma barra? Sim! O que nós queremos estudar é a variação de temperatura na barra quando transmitimos calor a ela. Vamos chamar de a temperatura de um ponto no instante . Temos então que a temperatura é função de duas variáveis e !
Os problemas de calor que vamos ver aqui têm sempre a mesma carinha. Nós vamos ter que resolver a seguinte equação diferencial parcial:
Onde é a derivada parcial da função temperatura em relação ao tempo e é a derivada parcial de segunda ordem da temperatura em relação a . Esse é uma constante de proporcionalidade que é sempre dada no enunciado, não precisamos nos preocupar com ela.
Repare que temos derivadas parciais exatamente porque é uma função de duas variáveis.
Resolução pelo método da separação de variáveis
Beleza, e na prática? Vamos ver como resolver essa equação com um passo a passo. Exercícios de equação de calor homogênea têm sempre essa cara:
Exemplo: Encontre a função da temperatura da barra que satisfaz as seguintes condições:
Esse sistema vai ser dado pelo enunciado sempre, tá? A diferença é que às vezes ele pode te dar números mesmo como valor de e , e te dá a expressão de .
Nosso objetivo é resolver a equação diferencial parcial tendo essas três condições que o problema deu, ou seja, achar a expressão de . Essa é a primeira condição
Esse termo representa a temperatura quando . Ou seja, é a função que dá a temperatura ao longo da barra no instante inicial. Portanto, temos um problema de valor inicial em relação à variável . Isso porque temos uma informação de no inicial.
Agora as outras duas condições. Elas querem dizer que, tanto para como para , independente do tempo , a temperatura é zero.
Bom, como temos informações sobre para dois pontos distintos da barra (ponto inicial e ponto final), vemos que em relação à , nós temos um problema de valores contorno.
Passo 1: Separação de variáveis
Cara, vimos que em relação à temos um problema de valor inicial e em relação a temos um problema de valores de contorno. A ideia aqui é tentar separar esses dois problemas! Esse método que vamos usar se chama Separação de variáveis.
Vamos imaginar que a função que procuramos possa ser escrita como um produto entre uma função apenas da variável e outra apenas da variável , dessa forma:
E aí a gente vai supor que a nossa solução tem esse formato. Ok, tem que satisfazer a equação de calor, não é?
Então, vamos substituí-la ai! Temos:
Repare que, derivando em relação a , é vista como constante, o mesmo ocorre para quando derivamos em relação a ! Beleza, agora vamos separar as variáveis nessa equação, de um lado deixamos tudo que tem e do outro tudo que tem :
Vamos deixar do lado de para facilitar as contas, tá? Agora vamos analisar essa relação que achamos com cuidado. Temos duas variáveis independentes e .
Digamos que nós escolhemos um valor fixo para . Assim a fração em vira uma constante. Aí, nós teríamos uma expressão fixa no lado esquerdo da equação (um valor constante), que deve ser igual ao lado direito da equação para qualquer , pois não fixamos .
Por termos uma igualdade, o lado direito precisa ser constante independentemente do valor de . E isso só acontece se essa fração
for igual a uma constante para qualquer valor de . Acontece o mesmo se a gente fixar e olhar para a fração em . Com fixo, a fração em precisa ser também constante independentemente de . Por isso podemos igualar essas duas frações a uma constante
É padrão colocar esse menos, mas não é necessário. Cê não precisa colocar. Mas fica atento que você precisa chegar na mesma solução que a gente aqui. Se não chegou, há algo errado. Vamos agora separar essa igualdade em duas equações. Olhando só para :
E para :
Logo, temos as duas equações
Ok, em que isso melhora nossa vida? Lembra que eu te falei que tínhamos dois problemas juntos? Agora eles estão separadinhos em duas equações diferentes.
Precisamos entender agora como vão ficar as condições que ele deu. Começando pela condição do problema de valores de contorno. A gente tinha . Como a gente supôs que , a gente escreve:
Repara que precisa ser diferente de zero porque se ele for igual a zero, toda a nossa solução zera, ou seja se , isso implica que . Mas a gente não quer . Isso é a solução trivial. Ela não serve para a gente.
Então, nós precisamos ter . O mesmo ocorre aqui:
A condição nós vamos guardar para o final.
Passo 2: Resolver o problema de contorno
Beleza, nós transformamos tudo o que o problema deu em dois problemas separados. Supomos que:
Aí temos o problema de valor inicial em relação a
E o problema de valores de contorno em relação a :
E agora? Nós vamos encontrar o e o um de cada vez e, assim, encontrar . Temos duas EDO’s para resolver, uma na variável e a outra na , qual resolver primeiro? A que tiver 2 condições de contorno, sempre!
Temos duas condições para , então começamos por ela.
Essa é uma EDO com problemas de valores de contorno que vai ser resolvida pelo método de problemas de autovalor, porque temos aquele ali. Buscando uma solução não-trivial , vamos avaliar a equação característica:
Vamos ter que analisar cada caso de . Quando a gente substitui as condições iniciais que temos nesse problema nas soluções para e :
só vamos encontrar soluções triviais, onde . Elas não nos interessam. Então vamos usar aqui só o caso , cuja a solução é
Para esse caso, a solução da equação característica por bhaskara dá:
Daí e .
Substituindo a primeira condição de contorno :
Agora substituindo a segunda , já sabendo que :
Como não pode zerar para essa solução não ser trivial também, precisamos que o seno seja igual a zero. E isso acontece para valores de ângulos múltiplos de , ou seja, . Daí, precisamos que
Substituindo em e esquecendo do .
Esses são os autovalores a as autofunções para o problema de contorno em . Como a gente está usando autovalores aqui, aquela constante vai ser desconsiderada mesmo, tá? Como se a gente fizesse . Então esse é a nossa solução para a EDO.
Passo 3: Tendo o valor de , podemos resolver a outra EDO, em
Temos a seguinte EDO de primeira ordem
Aplicando o método da separação de variáveis, temos:
Integrando ambos os lados
Substituindo :
Sabendo que , podemos escrever:
Calma, estamos quase lá.
Passo 4: Escrever como série e achar os coeficientes da série
Bom, esse que encontramos é uma fórmula geral de soluções que atendem à nossa EDO em e à EDO em . Sabemos também, pela própria definição de autofunções, que múltiplos dessa solução ou combinações lineares de soluções desse tipo também satisfazem as EDO’s.
Então, podemos escrever o seguinte:
Uma série infinita nada mais é do que uma combinação linear de infinitos termos, não é? Esse é o coeficiente dos termos da série, você já vai entender porque ele está ali!
Só temos um pequeno problema: e a condição ? Não usamos ela em momento algum, portanto, essa solução não atende essa condição. Bom, então vamos fazer ela passar a atender!
Veja que, substituindo na série, temos :
Temos, então, que encontrar coeficientes de forma que a soma dessa série seja igual à função Reconhece isso de algum lugar? Isso é uma Série de Fourier de uma função ímpar.
Mas não é uma função periódica, pelo contrário, ela só está definida em , ao longo do comprimento da barra. Então temos que fazer uma extensão ímpar de , para depois representar a função como uma série de senos.
Como a extensão é ímpar, teremos a seguinte fórmula para :
Achando o coeficiente , é só substituir na série que encontramos:
E temos a nossa solução com representação em série. Finalmente!
Então, a notícia ruim é que esses problemas de equação de calor são enormes ☹ A parte boa é que eles são todos iguais! 😊 Ou muito parecidos, pelos menos. 😉
Casos especiais
Vamos te dar agora algumas dicas para você ganhar tempo na hora de resolver esses problemas. Se você conseguir reconhecer um desses dois casos especiais de condições, isso pode facilitar muito sua vida. Vamos lá!
Temperatura nas extremidades fixada em zero
O primeiro caso é o que acabamos de ver, quando temos a temperatura das extremidades sempre zero, ou seja, . Isso nos leva a um problema com condições de Dirichlet homogêneas, que é o seguinte:
Sua solução é a seguinte:
Dependendo do seu tipo de prova, você pode mencionar esse caso especial e escrever direto essa solução sem precisar perder tempo provando. Se você conseguir, pode memorizar isso aqui então.
Mesmo que você precise resolver passo a passo, é uma boa conhecer essa solução, porque esse problema aparece em muitas questões!
Extremidades isoladas
Esse caso é parecido com o anterior. Em vez de a temperatura ser fixada em zero nas extremidades da barra, sua derivada em relação a é zero: .
Esse caso nos leva a um problema com condições de Neumann homogêneas:
Suas soluções são do tipo:
Há também a seguinte solução:
Novamente, você não precisa decorar isso aqui, mas é uma boa ter alguma noção dessa solução porque esse problema também aparece em diversas questões.
Passo a Passo
Vamos colocar aquele passo a passo esperto então? Partiu!
- Falar que a solução é do tipo
- Colocar a forma apresentada no passo 1 dentro da equação do calor
- Separar tudo que for função de de um lado e o que for de do outro
- Igualar tudo a
- Obter as EDO’s separadas: uma para e outra para
- Substituir a forma do passo 1 nas condições de contorno, lembrando que não queremos solução trivial!
- Resolver a EDO que tiver duas condições de contorno!
- Estudar os casos de , não podemos ter solução trivial! (aqui obtemos os autovalores, , e autofunções)
- Resolver a EDO que falta, introduzindo o valor de que você obteve no passo 8
- Fazer a combinação linear das soluções obtidas em 9 (vai aparecer a série de Fourier)
- Para obter o resolver a série de Fourier da condição inicial!
Obs: Viu que esses problemas usam muitos conhecimentos de outros tópicos de cálculo, né? Mas fica tranquilo, temos tudo aqui no site bem explicadinho. Para tentar te ajudar, marquei o nome desses tópicos em azul. Aí se ficar alguma dúvida nesses assuntos, pesquisa aqui no site, tá?
Casos especiais
Passo a Passo
Exercícios Resolvidos
Exercício Resolvido #1
UNICAMP, Cálculo 3, P3 2011.1, Questão 4 item b
Dado que
e solução da equação do calor
Determine a solução do problema de valor inicial e de contorno
Passo 1
Bom, agora temos a equação de calor
e uma solução
podemos afirmar que uma combinação linear de também será solução:
Passo 2
Beleza, agora precisamos verificar se essa solução atende as condições de contorno e inicial.
Vamos começar com as de contorno. O que a gente precisa verificar é se
Se a gente colocar em , teremos . Logo, independente do valor de na exponencial a gente vai atender a condição de contorno!
Da mesma forma, em , temos . De novo, atenderemos a condição para qualquer valor de !
Passo 3
Vamos pra condição inicial
Okay, então façamos :
Então está deve ser a Série de Fourier de , com uma extensão ímpar (pois a série é em senos). Então, mãos à obra! Essa extensão será dada por:
Passo 4
Essa função tem período , logo . Beleza, vamos achar os coeficientes da série . Ok, você pode chama-los de ...
Mas calma aí pequeno gafanhoto, e são funções ímpares, logo seu produto é par. Então vamos simplificar a integral...
Agora, vamos separar em cada intervalo que temos a expressão de bem definida!
Portanto, vamos ficar apenas com
Passo 5
Mas vamos olhar pra esses cossenos: se for par (), então
E essa diferença vai ser sempre . Viu?
Agora se é ímpar, , .
E tudo se reduz a
Passo 6
Assim:
Resposta
Exercício Resolvido #2
UNICAMP, Cálculo 3, P2 2014.1, Questão 4 item c
Usando separação de variáveis, encontrar a solução da seguinte equação do calor. Explique detalhadamente como se resolve o problema
onde é definida como a extensão par e periódica de período da função onde e esboce o gráfico no intervalo .
Passo 1
Para resolver o problema por separação de variáveis, vamos escrever
Note que se trata de uma equação de calor.
Com a separação, teremos:
Portanto, nossa equação se torna
E isso origina duas EDO’s (uma em e outra em ):
Passo 2
Agora, vamos olhar para as condições de contorno:
De e , se , então ambas as condições seriam atendidas. Mas como isso implicaria em o que daria a solução trivial.
Passo 3
Neste caso, nos resta admitir o seguinte:
Assim, formamos o seguinte problema:
Buscando uma solução não-trivial , vamos avaliar a equação característica:
Temos um problema de contorno e vamos ter que analisar cada caso de .
Passo 4
Caso 1: . Ou seja, o que está dentro da raiz é positivo. Nesse caso, . As soluções são do tipo:
Calculando .
Agora vamos usar as condições de contorno:
Isso nos dá
Substituindo esses valores na nossa solução geral, temos
Mas isso não nos ajuda em nada! Essa é a solução trivial, que não nos interessa. Portanto, esse caso, , não nos dá soluções para a EDO.
Passo 5
Caso 2: . Nesse caso, . As soluções são do tipo:
Calculando
Agora vamos usar as condições de contorno:
Isso nos deixa com a solução
Vamos calcular para esse
Nossa solução fica
Temos que analisar agora a condição e isso vai nos dar que.
Passo 6
Caso 3: . Ou seja, o que está dentro da raiz é negativo. Nesse caso, . As soluções são do tipo:
Substituindo
Calculando
Aplicando as condições de contorno
Note que na última equação não podemos ter (se não teríamos ). Então, precisamos de . Como , sobra
Passo 7
Como queremos , basta fazer , de modo que obtemos . Estes são nossos autovalores!
Portanto, ficamos com
que são nada mais nada menos que nossas autofunções!
Passo 8
Agora que já conhecemos os autovalores , vamos voltar pra equação em :
Fácil, né? Vamos substituir por :
E, assim, uma solução deve ser
Passo 9
Agora vem o pulo do gato. Se qualquer é solução para , então uma combinação linear dessas soluções (para diferentes valores de ) também é solução:
Assim:
Sabendo que
Para calcularmos os coeficientes da nossa função temos que tratar a expansão em série de Fourier da nossa função .
Passo 10
Vamos primeiro achar essa nossa função. Se a função é par, então . Bom, assumimos então que . E como o período é , vemos que . 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.
Onde é o período; portanto, .
Passo 12
Vamos avaliar os coeficientes da série:
Note que as descontinuidades de não afetam a integral.
Agora, temos:
Além disso:
Passo 13
Combinando os resultados, a série de Fourier é simplesmente
Passo 14
No nosso problema de equação de calor os coeficientes são justamente os coeficientes da série de Fourier de : .
Lembre que e , então e . Logo, sobra um termo apenas:
Na realidade, a solução é exatamente .
Resposta
Exercício Resolvido #3
William E. Boyce e Richard C. DiPrima, Equações diferenciais elementares e problemas de contorno, 9ª ed., Rio de Janeiro: LTC, 2010, pp. 326. Exercício 7
Encontre a solução do problema de condução de calor
Passo 1
Vamos fazer a separação de variáveis. Vamos supor que a função que procuramos possa ser escrita como a multiplicação de uma função apenas da variável por outra apenas da variável , dessa forma:
Vamos procurar uma solução dessa forma, com as variáveis separadas. Ok, tem que satisfazer a equação de calor, não é?
(Nesse caso, temos .) Então, vamos substituí-la ai! Temos:
Repare que, derivando em relação a , é vista como constante, o mesmo ocorre para quando derivamos em relação a . Agora vamos separar as variáveis nessa equação, de um lado deixamos tudo que tem e do outro tudo que tem :
Como essa igualdade deve valer para todo e para todo e essas variáveis são independentes, sabemos que essas frações são iguais a uma constante:
Vamos agora separar essa igualdade em duas equações:
Repara nos dados que o problema deu, vamos reescrevê-los usando a separação de variáveis :
Veja que se , isso implica que , encontramos a solução trivial. Então, nós precisamos ter . O mesmo ocorre aqui:
Passo 2
Resolver o problema de contorno.
Nós vamos encontrar o e o um de cada vez e, assim, encontrar . Temos duas EDO’s para resolver, uma na variável e a outra na , qual resolver primeiro? A que tiver condições de contorno, sempre!
Bom esse é um caso especial, as condições de Dirichlet homogêneas, não é? Nós já falamos sobre a solução desse problema na teoria desse capítulo, mas vamos fazer aqui passo a passo.
Temos a seguinte equação característica:
Caso 1: . Ou seja, o que está dentro da raiz é positivo. Nesse caso, . As soluções são do tipo:
Agora vamos usar as condições de contorno:
Isso nos dá
Substituindo esses valores na nossa solução geral, temos
Mas isso não nos ajuda em nada! Essa é a solução trivial, que não nos interessa. Portanto, esse caso, , não nos dá soluções para a EDO.
Caso 2: . Nesse caso, . As soluções são do tipo:
Agora vamos usar as condições de contorno:
Isso nos dá
Substituindo esses valores na nossa solução geral, temos
Novamente, a solução trivial, que não nos ajuda em nada. Beleza, vamos ver o 3º caso agora!
Caso 3: . Ou seja, o que está dentro da raiz é negativo. Nesse caso, . As soluções são do tipo:
Substituindo
Agora vamos usar as condições de contorno:
Isso nos dá
Agora vamos com calma! Se tivermos , a solução é a trivial , que não nos interessa. Então, para essa igualdade ai em cima ser verdadeira, precisamos que o seno seja zero.
O seno se anula em , , ... , múltiplos inteiros de . Logo, podemos dizer que:
As soluções do 3º caso são então do tipo:
Agora vamos resolver o problema em .
Passo 3
Tendo o valor de , resolver a outra EDO e achar a solução fundamental
Temos a seguinte EDO do 1º grau:
Como nós sabemos, a solução dessa EDO é:
Substituindo :
Sabendo que , podemos escrever:
Passo 4
Escrever como série e achar os coeficientes da série
Bom, esse que encontramos é uma fórmula geral de soluções que atendem à nossa EDO, sabemos que combinações lineares dessas soluções também satisfazem a EDO.
Então, podemos escrever o seguinte:
Uma série infinita nada mais é do que uma combinação linear de infinitos termos, não é? Beleza, agora temos que achar esse de forma que a nosssa solução atenda à seguinte condição do problema, que ainda não usamos:
Veja que, substituindo na série, temos :
Bom, se você reparar, aqui não precisamos nem calcular os ’s pela integral que dá os coeficientes de Fourier.
Chamando o termo geral da série de , vemos que é a soma de com , ou seja, para valores diferentes de e , os termos da série são nulos. Nós podemos fazer essa comparação porque a é função de senos.
Beleza, então temos que:
Os outros ’s são zero, então a série que procuramos é finita e só possui dois termos, com e . Logo, substituindo por 2 e 5:
Se você calculasse os ’s por integração, teria muito mais trabalho!