A equação de Laplace em coordenadas polares planas, equação 8.11, é complicada pelo raio variável. Considere a malha de diferenças finitas na figura abaixo, com nós (i,j) igualmente espaçados de Δ θ e Δ r . Deduza um modelo de diferenças finitas para a equação 8.11, semelhante à expressão cartesiana 8.109.

Passo 1

E aí, gente!! Tudo bem? Então vou mostrar pra vocês que mecânica dos fluidos não é nem de longe um bicho de 7 cabeças (nem de perto também)!

Foi pedido para modelar

1 r ∂ ∂ r   r ∂ ψ ∂ r + 1 r 2 ∂ 2 ψ ∂ θ 2 = 0

Passo 2

Nesse caso, usando diferenças finitas há duas possibilidades de prosseguirmos, dependendo de onde você queira truncar o primeiro termo. Eu sugiro fazer assim:

1 r i j Δ r r i j + Δ r 2 ψ i , j + 1 - ψ i , j Δ r - r i j - Δ r 2 ψ i , j - ψ i , j - 1 Δ r + ψ i + 1 , j - 2 ψ i j + ψ i - 1 , j     r i j 2 Δ θ 2 ≈ 0

Passo 3

Fazendo um ajuste bonitinho, ficamos com:

2 + 2 ζ ψ i j ≈ ψ i + 1 , j + ψ i - 1 , j + ζ 1 + η ψ i , j + 1 + ζ 1 - η ψ i , j - 1  

Onde ζ = r i j Δ θ Δ r 2 e η = Δ r / 2 r i j

Resposta

2 + 2 ζ ψ i j ≈ ψ i + 1 , j + ψ i - 1 , j + ζ 1 + η ψ i , j + 1 + ζ 1 - η ψ i , j - 1

Ver Outros Exercícios desse livro

Exercícios de Livros Relacionados