A virtual element method for the steady-state Poisson-Nernst-Planck equations on polygonal meshes. (15th November 2021)