We investigate and classify possible analytical solutions for a simplified version of the foam bubble population model, by varying injection conditions and kinetic foam generation parameter. We prove that the behavior of the analytical solutions changes at the transition between two regions, similar to rarefaction-shock solutions for the Buckley-Leverett equation. In one region (region I), the solutions are in the form of traveling waves, in rather good agreement with CT scanned foam experiments and numerical simulations reported in Simjoo and Zitha (2015). In region II, however, corresponds to solutions as a sequence of waves: one spreading wave and one traveling wave. These corresponding flow profiles are different from those found so far in the experiments.