The paper studies a simulation of a one-dimensional nonlinear free electron laser that includes a square hybrid wiggler, prebunched electron beam, and ion-channel guiding mechanism. The Lorentz equations of motion and the Maxwell equations were combined to derive a set of non-linear differential equations solved numerically within the framework of slowly varying amplitude and wavenumber. The relativistic electron beam is assumed to be cold, and self-fields of the electron beam and the slippage of the radiation wave are ignored. The appropriate parameters are chosen in such a way that the Raman regime is being considered, where there is high density and low energy. The effect of ion-channel densities on the saturation of the device is studied for both groups Ι and ΙΙ orbits. It is found that as the ion channel density increases, the range of saturation radiation for group Ι circuits increases, too. In contrast, the saturated radiation amplitude for group ΙΙ orbits is decreased with increasing the ion-channel density. Moreover, the results of the numerical calculations show that the saturation length for an FEL with a prebunched beam is considerably shorter than that of a uniform electron beam. A free electron laser (FEL) model called the two-stream FEL, which uses two relativistic electron beams of different energy, is also studied. Comparing the radiation amplitude between FEL and TSFEL shows that TSFEL is more efficient than FEL.