With the purpose of researching the internal flow stability of the vertical mixed-flow pump device (VMFPD), the entire flow field was solved through the very large-eddy simulation (VLES) k- ω model and the internal flow characteristics of the pump device were quantitatively analyzed with the use of the energy gradient theory (EGT) along with the physical parameters of the flow field. The research elucidates the variety rule of the pressure along the VMFPD and the energy gradient function K as well as the time-frequency features of the pressure pulsation (PP) in the impeller and guide vane. The results demonstrate that the solid wall restriction of the guide vane and the impeller rotation cause a significant shift in the velocity and direction of the water flow. This causes a significant pressure difference between the water bodies, which makes it easier to create an unstable flow. The suction surface of the impeller blade possesses the great concentration of the area with energy gradient function K > 107, which is the critical area influencing the internal flow stability of the pump device.