Guiding Center 4d Deeply Trapped with Lobatto Runge-Kutta Methods

The time step and number of time steps are those of this orbit, set in src/guiding-center-4d.jl; the toroidal momentum is shown alongside the energy, as the second invariant of the guiding centre dynamics.

LobattoIIIA(2)

Plots

Guiding Center 4d Deeply Trapped

LobattoIIIA(3)

Plots

Guiding Center 4d Deeply Trapped

LobattoIIIA(4)

Plots

Guiding Center 4d Deeply Trapped

LobattoIIIB(2)

Plots

Guiding Center 4d Deeply Trapped

LobattoIIIB(3)

Plots

Guiding Center 4d Deeply Trapped

LobattoIIIB(4)

Plots

Guiding Center 4d Deeply Trapped

LobattoIIIC(2)

Plots

Guiding Center 4d Deeply Trapped

LobattoIIIC(3)

Plots

Guiding Center 4d Deeply Trapped

LobattoIIIC(4)

Plots

Guiding Center 4d Deeply Trapped

LobattoIIID(2)

Plots

Guiding Center 4d Deeply Trapped

LobattoIIID(3)

Plots

Guiding Center 4d Deeply Trapped

LobattoIIID(4)

Plots

Guiding Center 4d Deeply Trapped

LobattoIIIE(2)

Plots

Guiding Center 4d Deeply Trapped

LobattoIIIE(3)

Plots

Guiding Center 4d Deeply Trapped

LobattoIIIE(4)

Plots

Guiding Center 4d Deeply Trapped

LobattoIIIF(2)

Plots

Guiding Center 4d Deeply Trapped

LobattoIIIF(3)

Plots

Guiding Center 4d Deeply Trapped

LobattoIIIF(4)

Plots

Guiding Center 4d Deeply Trapped

LobattoIIIG(2)

Plots

Guiding Center 4d Deeply Trapped

LobattoIIIG(3)

Plots

Guiding Center 4d Deeply Trapped

LobattoIIIG(4)

Plots

Guiding Center 4d Deeply Trapped