In this paper, we present a multi-fluid ionization model. We prove that this stationary, mono-dimensional model has a maximal solution which is not global at variance with the mono-species case and we present a numerical method for solving this highly singular system of ordinary differential equations. Numerical results are compared with those obtained for other models.