We examine the qualitative properties of ionic flows through membrane channels via Poisson–Nernst–Planck (PNP) type models with steric effects under relaxed electroneutrality boundary conditions, and more realistic setups in the study of ion channel problems. Of particular interest are the vital roles played by some critical potentials identified for both individual fluxes and current–voltage relations. These critical potentials split the whole electric potential interval into different subintervals, over which distinct dynamics of ionic flows are observed. The discussion provides an efficient way to control the boundary conditions to observe distinct dynamics of ionic flows through membrane channels. This is important for future analytical studies and critical for future numerical and even experimental studies of ion channel problems.