Numerical Solution of 3D Poisson-Nernst-Planck Equations Coupled with Classical Density Functional Theory for Modeling Ion and Electron Transport in a Confined Environment. Commun. Comput. Phys. 2014;16(5):1298-1322. doi:10.4208/cicp.040913.120514a