A conservative discretization of the Poisson–Nernst–Planck equations on adaptive Cartesian grids