A positivity-preserving and free energy dissipative hybrid scheme for the Poisson-Nernst-Planck equations on polygonal and polyhedral meshes. (15th February 2022)