A quasi-one-dimensional steady-state Poisson-Nernst-Planck model with Bikerman's local hard-sphere potential for ionic flows of two oppositely charged ion species through a membrane channel is analyzed. Of particular interest is the qualitative properties of ionic flows in terms of individual fluxes without the assumption of electroneutrality conditions, which is more realistic to study ionic flow properties of interest. This is the novelty of this work. Our result shows that i) boundary concentrations and relative size of ion species play critical roles in characterizing ion size effects on individual fluxes; ii) the first order approximation J k1 = D k J k1 in ion volume of individual fluxes J k = D k J k is linear in boundary potential, furthermore, the signs of ∂ V J k1 and ∂ 2 V λ J k1 , which play key roles in characterizing ion size effects on ionic flows can be both negative depending further on boundary concentrations while they are always positive and independent of boundary concentrations under electroneutrality conditions (see Corollaries 3.2-3.3, Theorems 3.4-3.5 and Proposition 3.7). Numerical simulations are performed to identify some critical potentials defined in (2). We believe our results will provide useful insights for numerical and even experimental studies of ionic flows through membrane channels.