A numerical algorithm for simulating the electric field distribution and particle trajectories in tribo-charge powder coating systems is presented. It combines the finite element method for the electric field analysis and the numerical integration of the equation for the particle trajectories. The particle-in-cell technique is used to couple both methods, which are applied iteratively until a self-consistent solution is obtained.