Dedução Leapfrog

De Física Computacional
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ΔtΔx.