Mudanças entre as edições de "Análise espectral e estabilidade"
m (→O problema) |
m (→O problema) |
||
Linha 6: | Linha 6: | ||
Queremos aproximar a solução da equação de Navier Stokes em um duto. Para isso devemos resolver a cada passo de tempo uma equação de Poisson como | Queremos aproximar a solução da equação de Navier Stokes em um duto. Para isso devemos resolver a cada passo de tempo uma equação de Poisson como | ||
<math>\frac{\partial^2 p}{\partial x^2}+\frac{\partial^2 p}{\partial y^2}= f(u)\,\!</math> | <math>\frac{\partial^2 p}{\partial x^2}+\frac{\partial^2 p}{\partial y^2}= f(u)\,\!</math> | ||
+ | |||
+ | ===Condição de Neumann em todos os lados=== | ||
+ | O espectro é real com autovalores entre (-1,1] incluindo 0 e 1. | ||
+ | |||
+ | ===Condição de Neumann em todos os lados fixando um ponto, P(1,2)=1 === | ||
+ | O espectro tem o maior autovalor λ=1.6537 e os outros autovalores próximos do eixo real (talvez devido a erros de ponto flutuante) entre (-1,1) incluindo 0. | ||
+ | |||
+ | ===Condição de Neumann em 3 lados fixando P=1 na entrada (ou na saída)=== | ||
+ | O espectro tem o maior autovalor λ=11.23 e os outros autovalores próximos do eixo real (com um pequeno circulo de valores complexos) entre (-1,1) incluindo 0. | ||
+ | ===Teste 4=== | ||
+ | |||
+ | ===Teste 5=== |
Edição das 23h04min de 21 de julho de 2009
A análise dos autovalores de uma matriz de iteração pode ser usada para estudar a estabilidade de um método iterativo.
Vamos relatar um estudo para um problema específico.
Índice
O problema
Queremos aproximar a solução da equação de Navier Stokes em um duto. Para isso devemos resolver a cada passo de tempo uma equação de Poisson como <math>\frac{\partial^2 p}{\partial x^2}+\frac{\partial^2 p}{\partial y^2}= f(u)\,\!</math>
Condição de Neumann em todos os lados
O espectro é real com autovalores entre (-1,1] incluindo 0 e 1.
Condição de Neumann em todos os lados fixando um ponto, P(1,2)=1
O espectro tem o maior autovalor λ=1.6537 e os outros autovalores próximos do eixo real (talvez devido a erros de ponto flutuante) entre (-1,1) incluindo 0.
Condição de Neumann em 3 lados fixando P=1 na entrada (ou na saída)
O espectro tem o maior autovalor λ=11.23 e os outros autovalores próximos do eixo real (com um pequeno circulo de valores complexos) entre (-1,1) incluindo 0.