O movimento das cargas é baseada na segunda lei de Newton e na lei de Coulomb da eletrostática.
A força resultante na carga $q_i$ é dada por
$$
\vec{F}_i = \dfrac{1}{4\pi\epsilon_0}\sum_{j \neq i}^N \dfrac{q_iq_j}{r_{ij}^2}\hat{r}_{ij},
$$
onde
- $N$: número de cargas livres
- $\epsilon_0$: constante eletrostática do vácuo
O movimento das cargas é amortecido por uma força proporcional à velocidade:
$$
\vec{F}_\gamma = -\gamma \vec{v}
$$
O sistema de EDO's é resolvido pelo método de Runge-Kutta de quarta ordem.
Outras considerações:
- Os valores de $q$, $\gamma$ e da massa das cargas são dados em unidades arbitrárias.
- As colisões das cargas com a fronteira do condutor são parcialmente elásticas.
Leia mais: