Invariant Regions for the Nernst-Plank Equations