A discrete bubble model has been developed taking into account multiple bubble-bubble interactions and a delayed coalescence method. The obtained simulation results were compared with experimental data reported in literature. The simulation results predicted by the developed model indicate clearly that the multiple interactions of bubbles lead to more reasonable results than those predicted by a binary interaction model. In addition, two types of interaction models were applied and predicted results were compared. The frequency of gas bubbles passing through the bed cross section versus bed height follows the same trend as the experimental data.