Equações de Laplace e Poisson - Eletrostática: mudanças entre as edições

De Física Computacional
Ir para navegação Ir para pesquisar
 
(29 revisões intermediárias pelo mesmo usuário não estão sendo mostradas)
Linha 91: Linha 91:


A melhoria na performance do algoritmo de Gauss-Siedel em comparação com o de Jacobi é apenas um fator 2.
A melhoria na performance do algoritmo de Gauss-Siedel em comparação com o de Jacobi é apenas um fator 2.
Conforme poderá ser visualizado pelos gráficos, não conseguimos distinguir uma diferença entre os resultados de ambos os métodos, pois são muito semelhantes. A principal vantagem do algoritmo de Gauss-Siedel reside no menor número de chamadas (atualizações) necessárias para que o sistema convirja para a tolerância escolhida. Com isso, ganhamos muito em tempo e gasto computacional.
=== Cálculo do campo elétrico ===
Já vimos como podemos calcular o potencial utilizando as equaçãoes de Laplace e de Poisson. Vamos agora calcular o campo elétrico a partir do potencial.
Vimos anteriormente que para encontrarmos o campo elétrio a partir do potencial, era necessário calcular um gradiente. Vamos então calcular a componente <math>x</math> do campo elétrico, uma vez que para as outras componentes, o cálculo é análogo. Podemos esolher entre as derivadas centrada, para frente ou para trás
<center><math> E_x(i,j) \approx - \frac{V(i+1,j) - V(i-1,j)}{2 \Delta x} </math></center>
<center><math> E_x(i,j) \approx - \frac{V(i+1,j) - V(i,j)}{\Delta x} </math></center>
<center><math> E_x(i,j) \approx - \frac{V(i,j) - V(i-1,j)}{\Delta x} </math>.</center>
A escolha do tipo de derivada numérica será devido à localização do potencial em que estivermos. Se estivermos nas bordas do contorno, então precisamos de uma derivada para a direita ou para esquerda, enquanto que se estivermos dentro do contorno, uma derivada centrada seria a melhor opção por ser mais simétrica.


== Exemplos ==
== Exemplos ==
Linha 96: Linha 114:
=== Caixa com bordas carregadas ===
=== Caixa com bordas carregadas ===


A caixa deste problema é bidimensional e ilustrada na figura ao lado. As laterais em x = 1 e x = -1 estão fixadas a um potencial elétrico de V = 1 e V = -1, respectivamente,pois são superfícies condutoras. As bordas superiores, em y = 1 e y = -1 são superfícies isolantes, logo o potencial varia de forma linear. No interior da caixa, o potencial é inicialmente zero pois são há cargas. Com a utilizando a equação de Laplace e as condições de contorno dadas, podemos resolver este problema analiticamente e numericamente.
[[Arquivo:placas_laplace1.png|thumb|right|Desenho esquemático da geometria do problema. Giordano e Nakanishi.|300px]]


[[Arquivo:placas_laplace1.png|thumb|right|Desenho esquemático da geometria do problema. Giordano e Nakanishi.|300px]]
A caixa deste problema é bidimensional e ilustrada na figura ao lado. As laterais em x = 1 e x = -1 estão fixadas a um potencial elétrico de V = 1 e V = -1, respectivamente, pois são superfícies condutoras. As bordas superiores, em y = 1 e y = -1 são superfícies isolantes, logo o potencial varia de forma linear. No interior da caixa, o potencial é inicialmente zero pois são há cargas. Com a utilizando a equação de Laplace e as condições de contorno dadas, podemos resolver este problema analiticamente e numericamente.


Como não há condutores no interior da caixa, esperamos que quando o tempo evoluir, o potencial no interior da caixa seja igual ao potencial nas superfícies isolantes  y = 1 e y = -1. Na direção <math>y</math> o potencial é constante, enquanto que na direção <math>x</math>, ele varia de V = -1 à V = 1.
Como não há condutores no interior da caixa, esperamos que quando o tempo evoluir, o potencial no interior da caixa seja igual ao potencial nas superfícies isolantes  y = 1 e y = -1. Na direção <math>y</math> o potencial é constante, enquanto que na direção <math>x</math>, ele varia de V = -1 à V = 1.
Linha 110: Linha 128:
Com a condição de contorno inicial, aplicamos o algoritmo de Jacobi e o algoritmo de Gauss-Siedel e escolhemos uma tolerância de <math>10^{-5} \times L^2</math>. Toda a iteração é feita somente dentro da caixa, sem modificarmos as bordas. Um chute inicial para a inicialização da matriz dentro das bordas é que todas os sítios estejam no potencial nulo.
Com a condição de contorno inicial, aplicamos o algoritmo de Jacobi e o algoritmo de Gauss-Siedel e escolhemos uma tolerância de <math>10^{-5} \times L^2</math>. Toda a iteração é feita somente dentro da caixa, sem modificarmos as bordas. Um chute inicial para a inicialização da matriz dentro das bordas é que todas os sítios estejam no potencial nulo.


A seguir apresentamos o potencial elétrico em um gráfico 3D, as linhas equipotenciais e o campo elétrico obtido através do potencial.
A seguir apresentamos o potencial elétrico em um gráfico 3D e o campo elétrico obtido através do potencial.


{|
{|
|[[Arquivo:potencial_laplace1.png|thumb|upright=1.1|center|Potencial elétrico para a caixa com bordas carregadas.|330px]]
|[[Arquivo:potencial_laplace1.png|thumb|upright=1.1|center|Potencial elétrico para a caixa com bordas carregadas.|500px]]
|[[Arquivo:equipotencial_laplace1.png|thumb|upright=1.1|center|Linhas equipotenciais para o interior do prisma externo.|330px]]
|[[Arquivo:electric_laplace1.png|thumb|upright=1.1|center|Campo elétrico para o interior do prisma externo.|500px]]
|[[Arquivo:electric_laplace1.png|thumb|upright=1.1|center|Campo elétrico para o interior do prisma externo.|330px]]
|}
|}


No gráfico do potencial elétrico, plotamos apenas metade dos pontos para que a imagem não ficasse extremamente poluída.
No gráfico do potencial elétrico, plotamos apenas metade dos pontos para que a imagem não ficasse extremamente poluída. Enquanto que no gráfico do campo elétrico, plotamos um terço dos vetores com o dobro do tamanho. A justificativa foi que se plotássemos tudo em sua forma completa, não seria possível ver o sentido dos vetores. A escolha utilizada foi a que mais trouxe uma boa visibilidade para o gráfico.
 
O número de atualizações para o algoritmo de Jacobi foi '''684''' enquanto que para o método de Gauss-Siedel foi de '''51'''.


=== Prisma com condutor interno ===
=== Prisma com condutor interno ===
Linha 141: Linha 160:
|[[Arquivo:electric_laplace2.png|thumb|upright=1.1|center|Campo elétrico para o interior do prisma externo.|330px]]
|[[Arquivo:electric_laplace2.png|thumb|upright=1.1|center|Campo elétrico para o interior do prisma externo.|330px]]
|}
|}
O número de atualizações para o algoritmo de Jacobi foi '''2721''' enquanto que para o método de Gauss-Siedel foi de '''1760'''.


=== Capacitores paralelos ===
=== Capacitores paralelos ===


Neste problema, vamos encontrar o potencial e o campo elétrico para duas placas paralelas submetidas a um potencial de V = -1 e V = 1. Na Figura $$ apresentamos um desenho esquemático do problema. Novamente, o problema é em duas  dimensões por possuir simetria em relação ao eixo z. A caixa na qual as placas estão inseridas estão fixadas a um potencial nulo.  
[[Arquivo:capacitores_laplace3.png|thumb|right|Desenho esquemático da geometria do problema. Giordano e Nakanishi.|300px]]
 
Neste problema, vamos encontrar o potencial e o campo elétrico para duas placas paralelas submetidas a um potencial de V = -1 e V = 1. Ao lado, apresentamos um desenho esquemático do problema. Novamente, o problema é em duas  dimensões por possuir simetria em relação ao eixo <math>z</math>. A caixa na qual as placas estão inseridas estão fixadas a um potencial nulo.  


Entretanto, diferente do problema anterior, agora nós temos uma dificuldade extra, caso quiséssemos atacar este problema de forma analítica. Temos efeitos de distorção devido às bordas finitas das placas, logo o campo elétrico não será, de forma alguma, simples nas vizinhanças superiores e inferiores das placas. Uma solução numérica é a única alternativa.
Entretanto, diferente do problema anterior, agora nós temos uma dificuldade extra, caso quiséssemos atacar este problema de forma analítica. Temos efeitos de distorção devido às bordas finitas das placas, logo o campo elétrico não será, de forma alguma, simples nas vizinhanças superiores e inferiores das placas. Uma solução numérica é a única alternativa.
Linha 150: Linha 173:
==== Resultado numérico ====
==== Resultado numérico ====


Vamos novamente escrever uma matriz de lado L = 201 e assumir que a caixa tem lado 1. Vamos inicializar todos os sítios em 0, menos nas placas, que estão a uma distância de 0.3 da origem. Escolhemos novamente a tolerância de <math>10^{-5} \times L^2</math>.  
Vamos novamente escrever uma matriz de lado L = 201 e assumir que a caixa tem lado 2. Vamos inicializar todos os sítios em 0, menos nas placas, que estão a uma distância de 0.3 da origem. Escolhemos novamente a tolerância de <math>10^{-5} \times L^2</math>. Aplicamos o algoritmo de Gauss-Siedel utilizando a mesma condição de contorno inicial.


Nas figuras posteriores apresentamos o potencial elétrico em um gráfico tridimensional, as linhas equipotenciais e o campo elétrico.
Nas figuras posteriores apresentamos o potencial elétrico em um gráfico tridimensional, as linhas equipotenciais e o campo elétrico.


Aplicamos o algoritmo de Gauss-Siedel utilizando a mesma condição de contorno inicial. Nas próximas figuras, mostramos o potencial em um gráfico em 3D, as linhas equipotenciais e o campo elétrico.
{|
|[[Arquivo:potencial_laplace3.png|thumb|upright=1.1|center|Potencial elétrico ao redor dos capacitores.|330px]]
|[[Arquivo:equipotencial_laplace3.png|thumb|upright=1.1|center|Linhas equipotenciais.|330px]]
|[[Arquivo:electric_laplace3.png|thumb|upright=1.1|center|Campo elétrico ao redor dos capacitores.|330px]]
|}
 
No gráfico do campo elétrico, plotamos somente um oitavo dos pontos e duplicamos os vetores para que ficasse uma boa visibilidade.
 
O número de atualizações para o algoritmo de Jacobi foi '''2146''' enquanto que para o método de Gauss-Siedel foi de '''1440'''.


=== Caixa com carga pontual ===
=== Caixa com carga pontual ===
[[Arquivo:caixa_poisson1.png|thumb|right|Campo elétrico para a carga pontual. Giordano e Nakanishi.|300px]]


Este é o primeiro problema em que cargas elétricas estão envolvidas, ou seja, temos que resolver numericamente a equação de Poisson ao invés da equação de Laplace. Vamos considerar uma caixa tridimensional (um cubo) de lado 2 e colocar a carga elétrica no centro desse cubo, conforme a Figura ao lado. Vamos fixar o potencial da caixa em V = 0 e calcular o potencial elétrico e o campo elétrico dessa carga. Devemos nos lembrar que a solução analítica para este problema, ou seja, o potencial elétrico para uma carga pontual é
Este é o primeiro problema em que cargas elétricas estão envolvidas, ou seja, temos que resolver numericamente a equação de Poisson ao invés da equação de Laplace. Vamos considerar uma caixa tridimensional (um cubo) de lado 2 e colocar a carga elétrica no centro desse cubo, conforme a Figura ao lado. Vamos fixar o potencial da caixa em V = 0 e calcular o potencial elétrico e o campo elétrico dessa carga. Devemos nos lembrar que a solução analítica para este problema, ou seja, o potencial elétrico para uma carga pontual é
[[Arquivo:caixa_poisson1.png|thumb|right|Campo elétrico para a carga pontual.|300px]]


<center><math> V(r) = \frac{1}{4 \pi \epsilon_0}\frac{q}{r} </math></center>
<center><math> V(r) = \frac{1}{4 \pi \epsilon_0}\frac{q}{r} </math></center>
Linha 174: Linha 205:
|[[Arquivo:potencial_poisson.png|thumb|upright=1.1|center|Campo elétrico para a carga pontual.|500px]]
|[[Arquivo:potencial_poisson.png|thumb|upright=1.1|center|Campo elétrico para a carga pontual.|500px]]
|}
|}
No gráfico do campo elétrico, os vetores foram diminuídos em 10 vezes para que a visibilidade ficasse boa. Vamos lembrar que a direção e sentido dos vetores e a diferença de tamanho entre eles que é  o fator relevante, não necessariamente seu módulo verdadeiro.
O número de atualizações para o algoritmo de Jacobi foi '''225''' enquanto que para o método de Gauss-Siedel foi de '''127'''.


== Códigos ==
== Códigos ==


Todos os códigos foram escritos em FORTRAN e podem ser acessados aqui.
Todos os códigos foram escritos em FORTRAN e podem ser acessados aqui. Cada link é direcionado à dois links diferentes, um com o código e outro com os comandos utilizandos para construir os gráficos no gnuplot.


[[ Exemplo 1 - Caixa com bordas carregadas ]]
[[ Exemplo 1 - Caixa com bordas carregadas ]]

Edição atual tal como às 00h16min de 6 de fevereiro de 2023

Nesta página resolvemos alguns problemas de eletrostática, através da resolução das equações de Laplace e de Poisson com condições de contorno. Apresentamos a solução analítica (se houver) e a solução numérica utilizando dois métodos. Os métodos numéricos aqui utlizados são denomindos de métodos de relaxação: método de Jacobi e método de Gauss-Seiedel.

Eletrostática

As equações de Laplace e Poisson surgem muito naturalmente no contexto da eletrostática. Para obtermos tais equações, comecemos com a Lei de Gauss para a eletrostática

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \nabla \cdot E = \frac{\rho}{\epsilon_0} }

onde Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle E} é o campo elétrico, Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \rho} é a densidade de carga elétrica e Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \epsilon_0} é a permissividade elétrica do vácuo.

Podemos relacionar o campo elétrico com o potencial elétrico pela seguinte relação

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle E = - \nabla V }

Substituindo na expressão anterior, encontramos a equação de Poisson

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \nabla^2 V = - \frac{\rho}{\epsilon_0}}

Se tomarmos o caso em que a densidade de cargas é nula, ou seja, não há cargas no meio, então obtemos a equação de Laplace

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \nabla^2 V = 0 }

Em problemas de eletrostática, o objetivo geralmente reside em encontrar o campo elétrico. Em muitos casos é mais fácil encontrarmos primeiro o potencial através das equações de Laplace ou Poisson e então, através de um gradiente, obtermos o campo elétrico. Isso ocorre devido à necessidade de termos condições de contorno para que as equações de Laplace e Poisson sejam devidamente resolvidas.

Nas próximas seções, resolveremos alguns problemas em que serão dados as condições de contorno para o potencial elétrico e através dessa informação será possível obter finalemente o campo elétrico.

Métodos númericos

Neste trabalho, resolvemos os problemas de eletrostática através dos chamados métodos de relaxação. Utilizamos dois algoritmos diferentes: o algoritmo de Jacobi e o algoritmo de Gauss-Siedel. Antes de apresentarmos estes algoritmos, vamos fazer algumas considerações gerais sobre o método de relaxação.

Vamos iniciar com a equação da difusão, uma equação que envolve derivadas parciais em relação oa tempo e espaço

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \frac{\partial \bar{V}(x,y,z,t)}{\partial t} = D \left( \frac{\partial^2 \bar{V}}{\partial x^2} + \frac{\partial^2 \bar{V}}{\partial y^2} + \frac{\partial^2 \bar{V}}{\partial z^2} \right) } .

Como as equações de Laplace e Poisson são equações independentes do tempo, as soluções estacionárias da equação da difusão serão também solução das equações de Laplace e Poisson. O processo de difusão tem um limite de equilíbrio, isto é, em algum momento o estado não vai depender do tempo. Conforme veremos a seguir, o algoritmo de Jacobi faz com que o sistema evolua de acordo com a equação de difusão. Mas como o sistema deve ser estacionário (não depender mais do tempo), a derivada de Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \bar{V}} com relação ao tempo deve se anular e recaímos na equação de Laplace.

A seguir apresentamos os dois algoritmos de relaxação.

Algoritmo de Jacobi

Tanto na equação de Laplace quanto na equação de Poisson, há derivadas segundas de uma função (no nosso caso aqui será sempre o potencial elétrico). Com isso, para discretizar as equações de Laplace e Poisson, precisamos discretizar as derivadas. Podemos utilizar tanto derivadas à direita, quanto à esquerda, quanto centradas

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \frac{\partial V}{\partial x} \approx \frac{V(i+1,j,k) - V(i,j,k)}{\Delta x} }


Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \frac{\partial V}{\partial x} \approx \frac{V(i,j,k) - V(i-1,j,k)}{\Delta x} }


Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \frac{\partial V}{\partial x} \approx \frac{V(i+1,j,k) - V(i-1,j,k)}{2 \Delta x} }

Como estamos interessados em derivadas segundas, podemos escrever estas como uma combinação da derivada centrada, de forma que

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \frac{\partial^2 V}{\partial x^2} \approx \frac{1}{\Delta x} \left[ \frac{\partial V}{\partial x} \left(i+\frac{1}{2}\right) - \frac{\partial V}{\partial x} \left(i-\frac{1}{2}\right) \right] }

onde Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \frac{\partial V}{\partial x} \left(i+\frac{1}{2}\right)} pode ser entendida como uma derivada centrada em um ponto imaginário Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \left(i+\frac{1}{2}\right)} que é idêntica a uma derivada à direita. O mesmo raciocínio pode ser aplicado para Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \frac{\partial V}{\partial x} \left(i-\frac{1}{2}\right)} . Dessa forma, a derivada segunda pode ser escrita como

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \frac{\partial^2 V}{\partial x^2} \approx \frac{1}{\Delta x} \left[ \frac{V(i+1,j,k) - V(i,j,k)}{\Delta x} - \frac{V(i,j,k) - V(i-1,j,k)}{\Delta x} \right] }

E rearanjando os termos, encontramos

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \frac{\partial^2 V}{\partial x^2} \approx \frac{V(i+1,j,k) + V(i-1,j,k) - 2 V(i,j,k)}{(\Delta x)^2} }

As derivadas parciais (primeiras e segungas) para Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle y} e Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle z} são completamente análogas. Se nós considerarmos que os passos ao longo de Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle x} , Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle y} e Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle z} são iguais, então Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \Delta x = \Delta y = \Delta z} , e podemos escrever a equação de Laplace como

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle V(i,j,k) = \frac{1}{6} \left[ V(i+1,j,k) + V(i-1,j,k) + V(i,j+1,k) + V(i,j-1,k) + V(i,j,k+1) + V(i,j,k-1) \right] } .

Isso significa que o valor do potencial em qualquer ponto é a média dos valores de Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle V} nas vizinhanças do ponto. Para um dado Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle V_0} inicial, nós calculamos um Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle V_1} através desta última expressão. Com Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle V_1} , é possível calcularmos Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle V_2} e assim por diante, até que a solução convirja para um dado valor. É necessário estabelecer uma tolerância para saber até onde precisamos aplicar o algoritmo. Com o intuito de enfatizar essa evolução temporal de Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle V} , vamos reescrever a última expressão como

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle V_{i,j,k}^{n+1} = \frac{1}{6} \left[ V_{i+1,j,k}^{n} + V_{i-1,j,k}^{n} + V_{i,j+1,k}^{n} + V_{i,j-1,k}^{n} + V_{i,j,k+1}^{n} + V_{i,j,k-1}^{n} \right] } .

onde os índices Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle i,j,k} se referem às coordenadas espaciais Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle x,y,z} e Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle n} se refere ao instante de tempo.

De forma completamente análoga, podemos escrever a equação de Poisson discretizada da seguinte forma

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle V(i,j,k) = \frac{1}{6} \left[ V(i+1,j,k) + V(i-1,j,k) + V(i,j+1,k) + V(i,j-1,k) + V(i,j,k+1) + V(i,j,k-1) \right] + \frac{\rho (i,j,k) \Delta x^2}{6} }

onde assumimos novamente que Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle \Delta x = \Delta y = \Delta z} .

Novamente reescrevendo a expressão para a equação de Poisson, obtemos

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle V_{i,j,k}^{n+1} = \frac{1}{6} \left[ V_{i+1,j,k}^{n} + V_{i-1,j,k}^{n} + V_{i,j+1,k}^{n} + V_{i,j-1,k}^{n} + V_{i,j,k+1}^{n} + V_{i,j,k-1}^{n} \right] + \frac{\rho_{i,j,k}^{n} (\Delta x)^2}{6} } .

Como este algoritmo é um algoritmo iterativo, precisamos escolher uma tolerância na qual o algoritmo pare de calcular a solução, pois ela já é boa o suficiente para o nosso objetivo. Devemos ter em mente, que uma tolerância cada vez menor, implica em uma solução cada vez mais próxima da solução correta, entretanto, isto significa um tempo e gasto computacional maior. É necessário analisar as vantagens e/ou necessidades de uma solução mais, ou menos, exata.

Algoritmo de Gauss-Siedel

O algoritmo de Gauss-Siedel é uma melhoria no algoritmo de Jacobi. Enquanto o algoritmo de Jacobi utiliza somente os valores antigos para construir o novo, o algoritmo de Gauss-Siedel utiliza os novos valores, assim que eles ficam disponíveis. Consideremos o caso bidimensional, vamos supor que você vai percorrer as linhas (indo do menor Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle i} para o maior Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle i} ) da sua cela, mantendo a coluna (Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle j} ) constante. Nesse caso, os valores Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle V_{i-1,j}^{n+1}} e Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle V_{i,j-1}^{n+1}} já serão conhecidos, de modo que podemos escrever o algoritmo de Gauss-Siedel como

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle V_{i,j}^{n+1} = \frac{1}{4} \left[ V_{i+1,j}^{n} + V_{i-1,j}^{n+1} + V_{i,j+1}^{n} + V_{i,j-1}^{n+1} \right] } .

A melhoria na performance do algoritmo de Gauss-Siedel em comparação com o de Jacobi é apenas um fator 2.

Conforme poderá ser visualizado pelos gráficos, não conseguimos distinguir uma diferença entre os resultados de ambos os métodos, pois são muito semelhantes. A principal vantagem do algoritmo de Gauss-Siedel reside no menor número de chamadas (atualizações) necessárias para que o sistema convirja para a tolerância escolhida. Com isso, ganhamos muito em tempo e gasto computacional.

Cálculo do campo elétrico

Já vimos como podemos calcular o potencial utilizando as equaçãoes de Laplace e de Poisson. Vamos agora calcular o campo elétrico a partir do potencial.

Vimos anteriormente que para encontrarmos o campo elétrio a partir do potencial, era necessário calcular um gradiente. Vamos então calcular a componente Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle x} do campo elétrico, uma vez que para as outras componentes, o cálculo é análogo. Podemos esolher entre as derivadas centrada, para frente ou para trás

Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle E_x(i,j) \approx - \frac{V(i+1,j) - V(i-1,j)}{2 \Delta x} }


Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle E_x(i,j) \approx - \frac{V(i+1,j) - V(i,j)}{\Delta x} }


Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle E_x(i,j) \approx - \frac{V(i,j) - V(i-1,j)}{\Delta x} } .

A escolha do tipo de derivada numérica será devido à localização do potencial em que estivermos. Se estivermos nas bordas do contorno, então precisamos de uma derivada para a direita ou para esquerda, enquanto que se estivermos dentro do contorno, uma derivada centrada seria a melhor opção por ser mais simétrica.

Exemplos

Caixa com bordas carregadas

Desenho esquemático da geometria do problema. Giordano e Nakanishi.

A caixa deste problema é bidimensional e ilustrada na figura ao lado. As laterais em x = 1 e x = -1 estão fixadas a um potencial elétrico de V = 1 e V = -1, respectivamente, pois são superfícies condutoras. As bordas superiores, em y = 1 e y = -1 são superfícies isolantes, logo o potencial varia de forma linear. No interior da caixa, o potencial é inicialmente zero pois são há cargas. Com a utilizando a equação de Laplace e as condições de contorno dadas, podemos resolver este problema analiticamente e numericamente.

Como não há condutores no interior da caixa, esperamos que quando o tempo evoluir, o potencial no interior da caixa seja igual ao potencial nas superfícies isolantes y = 1 e y = -1. Na direção Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle y} o potencial é constante, enquanto que na direção Falhou ao verificar gramática (MathML com retorno SVG ou PNG (recomendado para navegadores modernos e ferramentas de acessibilidade): Resposta inválida ("Math extension cannot connect to Restbase.") do servidor "https://wikimedia.org/api/rest_v1/":): {\displaystyle x} , ele varia de V = -1 à V = 1.

Como o campo elétrico é perpendicular às equipotenciais, o campo elétrico será inteiramente na direção e orientado na direção em que o potencial decresce, ou seja, para a esquerda.

Resultado numérico

Numericamente, construímos uma matriz quadrada de L = 61 linhas e colunas, com as devidas condições de contorno nas bordas. Utilizamos um lado de 61 para que o problema ficasse simétrico em relação à origem, logo temos 30 colunas à esquerda e 30 colunas à direita da coluna central. Esta abordagem simétrica foi utilizada em todos os exemplos resolvidos aqui.

Com a condição de contorno inicial, aplicamos o algoritmo de Jacobi e o algoritmo de Gauss-Siedel e escolhemos uma tolerância de . Toda a iteração é feita somente dentro da caixa, sem modificarmos as bordas. Um chute inicial para a inicialização da matriz dentro das bordas é que todas os sítios estejam no potencial nulo.

A seguir apresentamos o potencial elétrico em um gráfico 3D e o campo elétrico obtido através do potencial.

Potencial elétrico para a caixa com bordas carregadas.
Campo elétrico para o interior do prisma externo.

No gráfico do potencial elétrico, plotamos apenas metade dos pontos para que a imagem não ficasse extremamente poluída. Enquanto que no gráfico do campo elétrico, plotamos um terço dos vetores com o dobro do tamanho. A justificativa foi que se plotássemos tudo em sua forma completa, não seria possível ver o sentido dos vetores. A escolha utilizada foi a que mais trouxe uma boa visibilidade para o gráfico.

O número de atualizações para o algoritmo de Jacobi foi 684 enquanto que para o método de Gauss-Siedel foi de 51.

Prisma com condutor interno

Neste problema assumimos um prisma retangular infinito na direção . Este prisma é metálico e está sujeito a um potencial nulo fixo. O prisma condutor interno está sujeito a um potencial V = 1. Nosso objetivo é descobrir o potencial e o campo elétrico entre os dois prismas. Abaixo apresentamos um corte esquemático do prisma.

Desenho esquemático da geometria do problema. Giordano e Nakanishi.

É importante mencionar que este problema não possui solução analítica.Com isso, você pode se perguntar, como podemos testar o método, se não conhecemos a solução analítica do problema? A resposta para esta pergunta reside no fato de que nós podemos inferir qual deve ser o resultado esperado baseado na simetria do problema. Da mesma forma que no exemplo anterior, onde o potencial decresceu de V = 1 até V = -1, aqui nós esperamos que ele decresça de V = 1 no prisma central até V = 0 nas bordas. Qualquer resultado diferente disso, já nos mostra problemas com o método. Mas e como saberemos que o resultado mostrado, mesmo respeitando esta simetria, é o correto? Bom, aplicando o método a um número suficientes de problemas com resposta analítica conhecida nos dá confiança no método, de forma que podemos "confiar", com cautela, nos seus resultados fornecidos. Cm isso, vamos aos resultados numéricos.

Resultado numérico

Considerando o prisma exterior tendo um lado 2 e o condutor interno tendo um lado 0.6. Foi considerado um espaçamento espacial de 0.01, logo utilizamos uma matriz quadrada de tamanho L = 201.

Novamente, com a condição de contorno inicial, aplicamos o algoritmo de Jacobi e escolhemos uma tolerância de . Aplicamos o algoritmo de Gauss-Siedel utilizando a mesma condição de contorno inicial.

Abaixo apresentamos o potencial elétrico em um gráfico 3D, as linhas equipotenciais eo campo elétrico obtido através do potencial.

Potencial elétrico para o interior do prisma externo.
Linhas equipotenciais para o interior do prisma externo.
Campo elétrico para o interior do prisma externo.

O número de atualizações para o algoritmo de Jacobi foi 2721 enquanto que para o método de Gauss-Siedel foi de 1760.

Capacitores paralelos

Desenho esquemático da geometria do problema. Giordano e Nakanishi.

Neste problema, vamos encontrar o potencial e o campo elétrico para duas placas paralelas submetidas a um potencial de V = -1 e V = 1. Ao lado, apresentamos um desenho esquemático do problema. Novamente, o problema é em duas dimensões por possuir simetria em relação ao eixo . A caixa na qual as placas estão inseridas estão fixadas a um potencial nulo.

Entretanto, diferente do problema anterior, agora nós temos uma dificuldade extra, caso quiséssemos atacar este problema de forma analítica. Temos efeitos de distorção devido às bordas finitas das placas, logo o campo elétrico não será, de forma alguma, simples nas vizinhanças superiores e inferiores das placas. Uma solução numérica é a única alternativa.

Resultado numérico

Vamos novamente escrever uma matriz de lado L = 201 e assumir que a caixa tem lado 2. Vamos inicializar todos os sítios em 0, menos nas placas, que estão a uma distância de 0.3 da origem. Escolhemos novamente a tolerância de . Aplicamos o algoritmo de Gauss-Siedel utilizando a mesma condição de contorno inicial.

Nas figuras posteriores apresentamos o potencial elétrico em um gráfico tridimensional, as linhas equipotenciais e o campo elétrico.

Potencial elétrico ao redor dos capacitores.
Linhas equipotenciais.
Campo elétrico ao redor dos capacitores.

No gráfico do campo elétrico, plotamos somente um oitavo dos pontos e duplicamos os vetores para que ficasse uma boa visibilidade.

O número de atualizações para o algoritmo de Jacobi foi 2146 enquanto que para o método de Gauss-Siedel foi de 1440.

Caixa com carga pontual

Campo elétrico para a carga pontual. Giordano e Nakanishi.

Este é o primeiro problema em que cargas elétricas estão envolvidas, ou seja, temos que resolver numericamente a equação de Poisson ao invés da equação de Laplace. Vamos considerar uma caixa tridimensional (um cubo) de lado 2 e colocar a carga elétrica no centro desse cubo, conforme a Figura ao lado. Vamos fixar o potencial da caixa em V = 0 e calcular o potencial elétrico e o campo elétrico dessa carga. Devemos nos lembrar que a solução analítica para este problema, ou seja, o potencial elétrico para uma carga pontual é

onde é a carga da partícula e é a distância entre o ponto em que estamos calculando o potencial e a carga.

Resultado numérico

Para resolver com o método de Jacobi, inicializamos uma matriz com lado L = 21, ou seja, o tamanho dos sítios foi de 0.1. Utilizamos unidades em que e onde a constante dielétrica no vácuo é a unidade. Inicializamos uma matriz de mesmo tamanho com todos os sítios nulos exceto o sítio central e seus primeiros vizinhos (esquerda, direita, cima, baixo, frente e atrás). A seguir mostramos os gráficos do potencial e do campo elétrico para este problema.

Campo elétrico para a carga pontual.
Campo elétrico para a carga pontual.

No gráfico do campo elétrico, os vetores foram diminuídos em 10 vezes para que a visibilidade ficasse boa. Vamos lembrar que a direção e sentido dos vetores e a diferença de tamanho entre eles que é o fator relevante, não necessariamente seu módulo verdadeiro.

O número de atualizações para o algoritmo de Jacobi foi 225 enquanto que para o método de Gauss-Siedel foi de 127.

Códigos

Todos os códigos foram escritos em FORTRAN e podem ser acessados aqui. Cada link é direcionado à dois links diferentes, um com o código e outro com os comandos utilizandos para construir os gráficos no gnuplot.

Exemplo 1 - Caixa com bordas carregadas

Exemplo 2 - Prisma com condutor interno

Exemplo 3 - Capacitores paralelos

Exemplo 4 - Caixa com carga pontual

Referências

Giordano, N. J. Nakanishi, H. Computational Physics. Pearson Prentice Hall, 2ª ed. 2006.

Griffiths, D. Eletrodinâmica. Pearson 3ª ed. São Paulo.

Scherer, C. Métodos Computacionais da Física. Livraria da Física, 2ª ed. 2010.

Disponível em “http://fiscomp.if.ufrgs.br/index.php?title=Equações_de_Laplace_e_Poisson_-_Eletrostática&oldid=9375”