Dedução Leapfrog

De Física Computacional
Revisão de 16h28min de 24 de janeiro de 2018 por Asorander (discussão | contribs) (Criou página com 'Queremos resolver as equações que temos para <math>k_{i+\frac{1}{2}}^{n+1}</math>: <math> k_{i+\frac{1}{2}}^{n+1} = v \frac{U_{i+1}^{n+1}-U_i^{n+1}}{\Delta x} </math> Sabe...')
(dif) ← Edição anterior | Revisão atual (dif) | Versão posterior → (dif)
Ir para navegação Ir para pesquisar

Queremos resolver as equações que temos para ki+12n+1:

ki+12n+1=vUi+1n+1Uin+1Δx

Sabendo que

ki+12n=vUi+1nUinΔx,

e

sin+12=Uin+1UinΔt,

podemos escrever equações para Ui+1n+1, Uin+1 e Ui+1n:

Ui+1n+1=si+1n+12Δt+Ui+1n,(1)

Uin+1=sin+12Δt+Uin(2)

e

Ui+1n=Δxvki+12n+Uin.(3)

Substituindo as equações (1) e (2) na equação para ki+12n+1, obtemos:

Δxvki+12n+1=(si+1n+12Δt+Ui+1n)(sin+12Δt+Uin).

Ao substituirmos a equação (3) nessa última equação obtida, obtemos a equação citada no desenvolvimento do Método de Leapfrog, dada por

Δxvki+12n+1=si+1n+12Δt+Δxvki+12n+Uin(sin+12Δt+Uin).

E, por fim, dela obtemos a equação para k:

ki+12n+1=ki+12n+r(si+1n+12sin+12),

onde r=vΔvΔx.