ExF drift

This example demonstrates a single proton motion under uniform B and gravity fields.

using TestParticle
using TestParticle: get_gc
using OrdinaryDiffEq
using StaticArrays
using CairoMakie

# Analytic EM fields
B(x) = SA[0.0, 1e-8, 0.0]
E(x) = SA[0.0, 0.0, 0.0]

# Earth's gravity
F(x) = SA[0.0, 0.0, -TestParticle.mᵢ*9.8]

# Initial static particle
stateinit = let x0 = [1.0, 0.0, 0.0], v0 = [0.0, 0.0, 0.0]
   [x0..., v0...]
end
# Time span
tspan = (0, 1.0)

param = prepare(E, B, F, species=Proton)
prob = ODEProblem(trace!, stateinit, tspan, param)
sol = solve(prob, Vern9())
# drift in x-direction + free fall in z-direction
f = lines(sol, idxs=(3,1);
   figure = (; size = (800, 400), fontsize=18),
   axis = (; title="ExF Drift", xlabel="Z [m]", ylabel="X [m]", aspect = DataAspect())
)

This page was generated using DemoCards.jl and Literate.jl.