Dedução Leapfrog: mudanças entre as edições

De Física Computacional
Ir para navegação Ir para pesquisar
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...'
 
Asorander (discussão | contribs)
Sem resumo de edição
 
Linha 41: Linha 41:
</math>
</math>


onde <math>r = \frac{v\Delta v}{\Delta x}</math>.
onde <math>r = \frac{v\Delta t}{\Delta x}</math>.

Edição atual tal como às 17h10min de 25 de janeiro de 2018

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.