We present a numerical analysis of the stochastic population balance (SPB) theory for foam motion in porous media. The theory condenses into a set of non-linear partial differential equations in the saturation, pressure, and bubble density. We solve the equations using the IMPES method and perform sensitivity and parametric analysis. Finally, we compare the saturation profiles obtained numerically with those obtained from CT scan foam experiments. The agreement between the theory and experiments confirms that the stochastic population balance model describes adequately foam dynamics in porous media.