Study of stationary plasma thrusters using two-dimensional fully kinetic simulations