Based on the boundary integral method, a 3D bubble breakup model in a narrow flow field is established, and a corresponding computation program is developed to simulate the symmetrical and asymmetrical bubble breakup. The calculated results are compared with the experimental results and agree with them very well, indicating that the numerical model is valid. Based on the basic behavior of bubbles in a narrow flow field, the symmetrical and asymmetrical bubble breakup is studied systematically using the developed program. A feasibility rule of 3D bubble breakup is presented. The dynamics of sub-bubbles after splitting is studied. The influences of characteristic parameters on bubble breakup and sub-bubble dynamics are analyzed.