charged particle trajectories, calculation, trajectory, electrostatic
Version 5.0 introduces the possibility to study trajectories of charged particles movement in the plane-parallel and axisymmetric electrostatic fields. It utilizes our original approaches based on our implementation of Finite-Element technology1 and modern computational algorithms2.
Trajectory calculation uses the following data:
Viewing calculation results, you see:
Calculating trajectories QuickField uses following assumptions:
According to these assumptions, we can describe the trajectory (x(t), y(t), z(t)) of a charged particle in two-dimensional electrostatic field E(x, y) with Newton's system of differential equations:
We reorganize this system of three second degree equations into six first-degree equations and append the following additional equation:
Defining the length l(t) of the trajectory covered by the particle in time t. We integrate the resulting system using the Runge-Kutta-Merson method with automatically defined integration step. Numerical integration stops immediately before the finite element's boundary, the step leading outside of the element being excluded. At the last point in the element, we extrapolate the trajectory with cubical segment of its Taylor series relative to time and solve the resulting equation using Tartaglia-Cardano formula and taking into account possible decrease of the equation's degree in homogeneous or zero fields.