A Numerical solver of 3D Poisson Nernst Planck equations for functional studies of ion channels