This package provides test particle tracing in an analytic or numerical electromagnetic field via DifferentialEquations.jl and native solvers that are compatible with the DifferentialEquations interface.
In the Julia REPL,
julia> ]
pkg> add TestParticleVisualization via Makie and Plots are supported via recipes. Please refer to each visualization library's documentation for installations.
TestParticle.jl is designed to work together with OrdinaryDiffEq.jl. For example, a proton in a static magnetic field can be traced in SI units via
using TestParticle, OrdinaryDiffEq, StaticArrays
# Magnetic field
B(x) = SA[0, 0, 1e-8]
# Electric field
E(x) = SA[0,0, 0.0, 0.0]
# Initial conditions
stateinit = let x0 = [1.0, 0.0, 0.0], v0 = [0.0, 1.0, 0.1]
[x0..., v0...]
end
# Time span
tspan = (0, 20)
# Assemble particle + fields
param = prepare(E, B, species=Proton)
prob = ODEProblem(trace!, stateinit, tspan, param)
# Trace trajectory and save positions & velocities
sol = solve(prob, Vern9())Native Boris particle pusher also follows a similar interface:
dt = 3e-11 # fixed time step
savestepinterval = 10
prob = TraceProblem(stateinit, tspan, param)
sol = TestParticle.solve(prob; dt, savestepinterval)[1]For plotting with Makie,
using GLMakie
plot(sol, idxs=(1, 2, 3))More tutorials and examples can be found in the doc.