Abstract: By adopting VOF method and a momentum equation coupled with surface tension model and in combination with the function of UDF, the interactions between two horizontal distribution bubbles and the effects of swinging motion on bubbles’ motion characteristics in a vertical channel were numerically simulated. The attraction and repulsion of the two bubbles in process of rising were analyzed in detail. The change in form of the two rising bubbles given by simulation agrees well with the experimental results. The results show that in the process of rising, the two horizontal distribution bubbles will swing up and down, and the position relation between them cyclically changes. Swinging motion makes the bubble become more flat, and the minimum distance between the two bubbles can be smaller. In addition, the swinging motion will destroy the cyclicity of position relation between the two bubbles in process of rising.