Equação de Liouville-Bratu-Gelfand

De Física Computacional
Ir para navegação Ir para pesquisar

Equação de Liouville-bratu-Gelfand

Na matemática, a Equação Liouville–Bratu–Gelfand ou Equação de Liouville é uma equação de Poisson não linear, nomeada em homenagem aos matemáticos Joseph Liouville, Gheorghe Bratu e Israel Gelfand, que é descrita da seguinte forma

Essa equação aparece em problemas de fuga térmica, como na teoria de Frank-Kamenetskii, e na astrofísica, por exemplo, na equação Emden–Chandrasekhar. Esta equação pode descrever a carga espacial de eletricidade em torno de um fio brilhante ou até mesmo uma nebulosa planetária.

A solução de Liouville

Em duas dimensões, com coordenadas cartesianas (x,y), Joseph Liouville propôs uma solução em 1853 como

onde é uma função analítica arbitrária com . Em 1915, G.W. Walker[1] encontrou uma solução assumindo uma forma para . Se , então a solução de Walker é

onde é algum raio finito. Essa solução vai ao infinito para qualquer , mas vai ao infinito na origem , finito na origem para e vai a zero na origem para . Walker também propôs mais duas soluções em seu artigo de 1915.


Formas radialmente simétricas

Se o sistema a ser estudado for radialmente simétrico, então a equação em dimensões torna-se

onde é a distância a partir da origem. Com as condições de contorno

e para , uma solução real existe apenas para , onde é o parâmetro crítico chamado de parâmetro de Frank-Kamenetskii. O parâmetro crítico é para , para e para . Para , existem duas soluções e para existem infinitas soluções com as soluções oscilando em torno do ponto . Para , a solução é única e, nesses casos, o parâmetro crítico é dado por . A multiplicidade de soluções para foi descoberta por Israel Gelfand em 1963 e, posteriormente, em 1973, generalizada para todos os por Daniel D. Joseph e Thomas S. Lundgren.[2]

A solução para que é válida no intervalo é dada por

onde está relacionada a como

A solução para que é válida no intervalo é dada por

onde está relacionada a como


Método Crank-Nicolson

Para resolver a equação diferencial parcial de Liouville-Bratu-Gelfand foi utilizado o método de Crank-Nicolson. Sendo a média ponderada dos métodos explícito e implícito de FTCS (Forward Time Central Space) da forma:

Falhou ao verificar gramática (erro de sintaxe): {\displaystyle \frac{\psi_j^{n+1}-\psi_j^{n}}{\Delta t} = \theta . implícito + (1-\theta) . explícito }

Método de Relaxação

Como a equação de Liouville-Bratu-Gelfand não depende do tempo, é necessária uma aproximação para uma equação similar dependente do tempo para resolvê-la numericamente aplicando o método de FTCS (Foward Time Central Space). Fazendo então a solução dessa equação convergir ao estado estacionário diante de uma evolução temporal longa o suficiente ().

Utilizando uma equação da difusão genérica chegamos na seguinte forma:

.

Onde é a constante de difusão. Assim é possível aplicar essa equação no método de Jacobi, método numérico de relaxação.

Referências

  1. https://en.wikipedia.org/wiki/Liouville%E2%80%93Bratu%E2%80%93Gelfand_equation
  2. Scherer, CLÁUDIO. Métodos Computacionais da Física. 2010.
  1. Walker, George W. "Some problems illustrating the forms of nebulae." Proceedings of the Royal Society of London. Series A, Containing Papers of a Mathematical and Physical Character 91.631 (1915): 410-420.https://www.jstor.org/stable/pdf/93512.pdf?refreqid=excelsior%3Af4a4cc9656b8bbd9266f9d32587d02b1
  2. Joseph, D. D., e T. S. Lundgren. "Quasilinear Dirichlet problems driven by positive sources." Archive for Rational Mechanics and Analysis 49.4 (1973): 241-269.