Abstract: In order to investigate the spatial distribution characteristics of breaking waves and the entrainment characteristics of bubbles in the flow field during submarine surface sailing. A self-developed adaptive mesh two-phase flow solver based on curvature criterion was used to simulate the flow field of submarine surface sailing. By using the self-developed bubble recognition program, the simulation data were statistically analyzed, and the characteristics of bubble size distribution and space distribution in the flow field were obtained. By using the curvature criterion, the number of grids required for calculation is effectively reduced. When the submarine is sailing on the surface, the volume of the bubbles in the wake region accounts for more than 80% of the total volume of bubbles in the flow field. With the increase of submarine Froud number from 0.28 to 0.38, the intensity of wake breaking significantly rises, the total number of bubbles increases by about 60%, the total entrainment volume of bubbles increases by about 50%, and the main distribution position of bubbles gradually tends to the outside and rear area of the submarine. This study can provide a reference for exploring the hydrodynamic characteristics of submarine surface sailing.