An automatic mesh generation method employing a dynamic bubble system can provide a high-quality mesh for electromagnetic finite element analysis. However, it takes very long time to compute bubbles' movement when a mesh with a huge number of elements is generated. The bubbles move according to a force among them like the van der Waals force. However, the force received from bubbles of very far distance can be ignored. Therefore, it is unnecessary to calculate the force from all the bubbles. In this paper, we propose the subdivided dynamic bubble system in order to shorten the computation time. The proposed method considers an analysis region to divide into some tetrahedral subdivisions, and then the bubbles moves independently in the tetrahedral subdivisions. As the result, the acceleration of the computation time is achieved.Index Terms-Dynamic bubble system, finite element analysis, mesh generation, parallel computing.