Integrador simplético

Em matemática, um integrador simplético (ou simpléctico) é um métodos numéricos para equações diferenciais ordinárias para sistemas hamiltonianos.

Introdução

Um sistema hamiltoniano é localmente descrito por um aberto $𝒫\subset {ℝ}^N\times {ℝ}^N$ e uma função $\mathscr{H}:𝒫\to ℝ$, onde $N$ é o número de graus de liberdade do sistema, o espaço de fase é localmente homeomorfo ao conjunto $𝒫$ e a função $\mathscr{H}$ é chamada de hamiltoniano do sistema. A evolução do sistema satisfaz às equações de Hamilton$${\dot{q}}^n=\frac{\partial }{\partial p_n}\mathscr{H},{\dot{p}}_n=-\frac{\partial }{\partial q^n}\mathscr{H}.$$Neste artigo a solução para essas equações será denotada por $(q(t),p(t))$.

Definição

De acordo com o teorema de Liouville, dada uma região $\mathrm{\Omega }$ do espaço de fase, se cada ponto dessa região evolui de acordo com as equações de Hamilton então o volume dessa região é constante, essa é uma propriedade importante dos sistemas hamiltonianos e portanto temos interesse em encontrar métodos numéricos de integração que preservem essa propriedade, um integrador simplético respeita essa propriedade. Em outras palavras, a evolução do sistema é um simplectomorfismo e um integrador simplético é, por definição, um método cuja discretização é também um simplectomorfismo.

Dado que o integrador simplético respeita as propriedades geométricas do espaço de fase e das equações de Hamilton ele é indicado para o problema de muitos corpos e para integração de tempo longo em sistemas hamiltonianos, obtendo resultados consideravelmente satisfatórios se comparado com outros métodos como RK4.

Exemplos

Métodos de primeira ordem

Os integradores simpléticos de primeira ordem são $𝒪(\delta t)$, onde $\delta t$ é o tamanho do passo para discretização do tempo. Os métodos de primeira ordem são variantes do método de Euler ajustados de acordo com a geometria simplética do espaço de fase.

Euler simplético I:

• $\mathscr{H}:𝒫\to ℝ$
• $p(t+\delta t)$ previamente estimado utilizando outro método (na primeira iteração).

Euler simplético II:

• $\mathscr{H}:𝒫\to ℝ$
• $q(t+\delta t)$ previamente estimado utilizando outro método (na primeira iteração).

Predefinição:NumBlkEsses métodos são, em geral, implícitos, entretanto em alguns casos particulares eles dão lugar a métodos explícitos, por exemplo:

Euler semi-implícito:

• $\mathscr{H}(q,p)=T(p)+V(q)$

Teorema I: Os métodos (Predefinição:EquationNote), (Predefinição:EquationNote) e (Predefinição:EquationNote) são simpléticos de ordem 1.

Demonstração: Cf. referências e . Note que o método (Predefinição:EquationNote) é um caso particular do método (Predefinição:EquationNote).

Métodos de segunda ordem

Os integradores simpléticos de segunda ordem são $𝒪(\delta t^2)$. Em geral, métodos de segunda ordem são variantes do método de Verlet, bem como o próprio método de Verlet.

Störmer-Verlet I:

• $\mathscr{H}:𝒫\to ℝ$
• $p(t+\delta t/2)$ e $q(t+\delta t)$ previamente estimados utilizando outro método (na primeira iteração).

Predefinição:NumBlkStörmer-Verlet II:

• $\mathscr{H}:𝒫\to ℝ$
• $q(t+\delta t/2)$ e $p(t+\delta t)$ previamente estimados utilizando outro método (na primeira iteração).

Predefinição:NumBlkComo acontece para os métodos de Euler, o caso do hamiltoniano separável permite desenvolver métodos explícitos derivados desses, que são implícitos.

Verlet:

• $\mathscr{H}(q,p)=T(p)+V(q)$

Teorema II: Os métodos (Predefinição:EquationNote), (Predefinição:EquationNote) e (Predefinição:EquationNote) são simpléticos de ordem 2.

Demonstração: Cf. referência e . Note que o método (Predefinição:EquationNote) é um caso particular do método (Predefinição:EquationNote).

Métodos de terceira ordem

Os integradores simpléticos de terceira ordem são $𝒪(\delta t^3)$.

Ronald D. Ruth:

• $\mathscr{H}(q,p)=T(p)+V(q)$
• $c_1+c_2+c_3=1$, $d_1+d_2+d_3=1$, $c_2{d}_1+c_3(d_1+d_2)=1/2$, $c_2+d_1^2+c_3(d_1+d_2{)}^2=1/3$ e $d_3+d_2(c_1+c_2{)}^2+d_1{c}_1^2=1/3$.

Predefinição:NumBlkPor exemplo: $c_1=7/24$, $c_2=3/4$, $c_3=-1/24$, $d_1=2/3$, $d_2=-2/3$ e $d_3=1$.

Teorema III: O método (Predefinição:EquationNote) é simplético de ordem 3.

Demonstração: Cf. referência .