We perform one-dimensional fluid simulation of ion acoustic (IA) solitons propagating parallel to the magnetic field in electron-ion plasmas by assuming a large system length. To model the initial density perturbations (IDP), we employ a KdV soliton type solution. Our simulation demonstrates that the generation mechanism of IA solitons depends on the wavelength of the IDP. The short wavelength IDP evolve into two oppositely propagating identical IA solitons, whereas the long wavelength IDP develop into two indistinguishable chains of multiple IA solitons through a wave breaking process. The wave breaking occurs close to the time when electrostatic energy exceeds half of the kinetic energy of the electron fluid. The wave breaking amplitude and time of its initiation are found to be dependent on characteristics of the IDP. The strength of the IDP controls the number of IA solitons in the solitary chains. The speed, width, and amplitude of IA solitons estimated during their stable propagation in the simulation are in good agreement with the nonlinear fluid theory. This fluid simulation is the first to confirm the validity of the general nonlinear fluid theory, which is widely used in the study of solitary waves in laboratory and space plasmas. © 2013 AIP Publishing LLC.