Invariant regions for the Nernst-Planck equations