Abstract: The Taylor bubble rising in a vertical tube filled with liquid lead-bismuth eutectic was simulated based on MPS method. Several performed results including shape development, bubble rising terminal speed fitting line, and the thickness and axial velocity profile of the falling film were presented. The simulation results were compared with some experimental results from literature and reveal the corresponding laws obeyed by the terminal speed, velocity profile and film thickness. The computational results agree well with both theoretical analysis and experimental results, which demonstrates the reasonable selection of model as well as the accuracy and reliability of MPS method.