Equação De Laplace No Disco
Teoria
Fala aí! Hoje vamos aprender sobre a equação de Laplace e sua aplicação em um disco.
Tá preparado? Então vem comigo!
⭐ Veja também ⭐
- EDPs gerais - Introdução e Método de Resolução.
- Série de Fourier.
- Métodos de resolução de EDOs de ordem superior.
O que é a Equação de Laplace?
A Equação de Laplace é uma equação diferencial parcial usada em vários ramos da ciência. A fórmula para a equação de Laplace é essa aqui:
Mas como aqui estamos interessados no estudo somente para o disco, vamos ter que pensar num jeito mais adequado de descrever essa equação aí.
A solução vem na forma das coordenadas polares!
Se liga, a região de interesse (o disco) será delimitada pela equação da circunferência.
Suponha que
Não convém demonstrar como chegar nessa equação, porque iria dar um trabalhão, e nosso foco aqui é saber como utilizá-la, tranquilo?
Bora ver como resolver essa parada aí?!
Como resolver a equação de Laplace?
Um tópico que volta a aparecer aqui é a resolução de equações diferenciais por separação de variáveis. Esse conteúdo nós já temos bem explicadinho aqui, então, se tiver precisando de uma revisada, corre lá que eu te espero 😁
Para resolver problemas da equação de Laplace, vamos sempre separar nossas variáveis na função.
Essas duas funções
Pelo método de separação de variáveis, chegamos em:
Agora, passamos tudo que tem R para um lado e tudo o que tem
Aí, a gente pode igualar essa expressão a uma constante arbitrária. Aqui, eu escolhi
Com isso, ficamos com duas EDOs que podemos resolver separadamente:
Análise de
Vamos supor que nossa função
Concorda que a temperatura irá se repetir à cada volta no disco? Então, a função deve ser
Quer dizer que, quando resolvermos uma questão e acharmos várias soluções para
Bora avaliar as faixas de valores de
Na equação:
Podemos ter
Como a derivada segunda é zero, se integrarmos isso aí duas vezes, chegamos em:
Entretanto, como
Agora vamos testar
Essa relação aí vai simplificar nossas contas.
E essa aí é periódica para todo
Pra finalizar, testamos
Essa solução não é periódica. Então, teremos:
Nossa função não pode ser zero. Então,
Análise de R(r)
Quando analisarmos a função R(r), teremos algumas soluções com fatores que explodem para infinito dentro dos limites do círculo. Esses fatores têm que ser desconsiderados.
Bora voltar pra equação:
Para
A solução
Vamos integrar essa equação para achar R(r). No final das contas, chegamos no seguinte:
E, com algumas manipulações algébricas:
Integra novamente:
Repare que essa equação não é limitada dentro do círculo que estamos estudando, pois
Quando testarmos
A nossa solução final é dada pelas soluções que encontramos em cada um dos casos, beleza?
Então bora dar sequência com mais alguns exercícios.
Como resolver a equação de Laplace?
Análise de
Análise de R(r)
Exercícios Resolvidos
Exercício Resolvido #1
William E. Boyce e Richard C. DiPrima, Equações diferenciais elementares e problemas de contorno, 9ª ed., Rio de Janeiro: LTC, 2010, pp. 334. Exercício 11
Encontre uma solução da equação de Laplace no interior do círculo , que satisfaça a condição de contorno sobre o círculo:
Passo 1
Enunciando a Equação de Laplace em coordenadas polares:
Fazendo a separação de variáveis:
Substituindo na equação e utilizando notação conveniente:
Passo 2
Agrupando termos:
E isolando:
Igualando ambas as expressões a uma constante arbitrária:
Passo 3
Agora vamos avaliar quais faixas de valores de são possíveis para a solução. Estudando a equação em , tentamos :
Entretanto, como deve ser -periódica, essa solução só funciona para . Portanto temos:
Passo 4
Testando , fazemos:
Supondo por simplicidade, de forma que . Seguindo:
Que tem equação característica:
Que resulta na solução:
Como essa solução é -periódica para todo , está valendo.
Passo 5
Por fim, testamos . Fazendo:
Plugando na equação:
Que tem equação característica:
Resultando na solução:
Esta solução não é -periódica para quaisquer valores de , portanto concluímos que:
Passo 6
Como para todo a solução trigonométrica funciona, fazendo (por elegância), temos então a autofunção de :
E a solução geral seria:
Substituindo os valores válidos de na equação em , temos para :
Arrumando:
Como buscamos uma solução para todo , a solução não interessa. Segue:
Passo 7
Arrumando mais para podermos integrar:
Integrando ambos os lados:
Pra quem lembra dos macetes de Cálculo 1:
Caso você não tenha sacado, pode resolver de outros jeitos.
Passo 8
Integrando novamente:
Agora analisando a função encontrada, temos que ela não é limitada na região que estamos estudando, pois estamos no interior do círculo , e temos:
Portanto nos resta concluir que:
A solução fica:
Passo 9
Testando a outra faixa possível , fazemos , :
Que é a famosa equação de Euler. Fazendo a solução genérica:
Plugamos:
Como para todos valores de :
As soluções então são:
Passo 10
Como estamos no interior do círculo , precisamos que a função não exploda para infinito em nenhum ponto desta região, logo como:
Concluímos que , portanto temos a autofunção:
Plugando na autofunção geral:
Sem esquecer que temos também a solução constante:
A solução geral portanto fica:
Passo 11
Precisamos agora determinar os coeficientes e . Utilizando a última condição que sobrou do enunciado:
Basta derivar com respeito à e depois fazer e temos:
Como à esquerda temos uma série de Fourier da função , os coeficientes podem ser encontrados pelas fórmulas:
Como a função é -periódica, temos que . Descobrimos os coeficientes:
Resposta
Exercício Resolvido #2
William E. Boyce e Richard C. DiPrima, Equações diferenciais elementares e problemas de contorno, 9ª ed., Rio de Janeiro: LTC, 2010, pp. 334. Exercício 6
Encontre a solução da equação de Laplace na região semicircular , , que satisfaz as condições de contorno:
Suponha que está bem definida e limitada na região dada.
Passo 1
Enunciando a Equação de Laplace em coordenadas polares:
Fazendo a separação de variáveis:
Substituindo na equação e utilizando notação conveniente:
Passo 2
Agrupando termos:
E isolando:
Igualando ambas as expressões a uma constante arbitrária:
Passo 3
E separando também as condições iniciais:
Como a solução é trivial, nós a descartamos. Portanto segue que:
Juntando estes valores com a equação em obtemos:
Passo 4
Para resolver a equação em precisamos considerar as possíveis faixas de valores de . Testando:
Substituindo os valores das condições:
Portanto temos a solução trivial:
Passo 5
Tentando , fazemos:
Para garantir que lambda seja positivo. Além disso, vamos substituir uma solução da forma:
E achar o . Plugando na EDO:
Passo 6
Ou seja, como a equação característica tem raízes complexas, temos que as soluções são da forma:
Naturalmente vamos utilizar a forma da direita por ser mais simples. Substituindo os valores vamos tentar descobrir as constantes:
Show de bola. A função fica:
Agora substituindo a outra:
Passo 7
Agora vejamos, para não cair na solução trivial de novo (), podemos buscar quais resolvem a equação acima. O que equivale à resolver:
Que tem solução:
Portanto, como só nos interessa :
Passo 8
Portanto a nossa solução geral fica:
Testando a gora a última faixa de valores possíveis , fazemos:
Temos que a nossa equação característica fica:
Que resulta em:
Como essa função não é periódica em podemos descartar. Daria no mesmo substituir os valores e chegar na solução trivial:
Resolvendo obtemos . Portanto só existem soluções não triviais para .
Passo 9
Seguindo, podemos voltar à resolver a equação em . Fazendo :
Supomos uma solução da forma;
Substituindo:
Temos:
Pondo em evidência:
Como para todo , temos que:
Ou seja, a solução geral é:
Passo 10
Agora o papo é analisar essa função na região estudada. Como a região é tal que:
Então precisamos garantir que dentro dela a função seja limitada, ou seja, não exploda para infinito. Como temos que:
Concluímos que:
Portanto a solução é:
Plugando na solução mais geral, sabendo que :
Passo 11
Falta nós determinarmos os coeficientes . Sabemos que, pelo enunciado, a informação que ficou sobrando é:
Que significa:
Podemos ver que à esquerda temos uma série de Fourier para a função , portanto os coeficientes são determinados pela fórmula:
Onde é o período da função. Como sabemos que, por estar num semicírculo, ela deve ser -periódica, temos: