Turbulent flow at a wide range of Mach numbers may occur during design based accidents (DBA) or severe accidents in nuclear power plants (NPPs), for instinct, incompressible gas mixing flow (subsonic), gas mixture critical flow (transonic) or detonation (supersonic). All of them need to be simulated in safety analysis of NPP containments, therefore we need the computer code to have the capability to simulate all-speed flow. GASFLOW-MPI is a well-developed parallel all-speed CFD software used to predict gas turbulent mixing, chemical combustion and other related thermal-hydraulic phenomena in containments of NPPs. In order to validate the all-speed turbulent flow capability of GASFLOW-MPI, four turbulent jet cases over a broad range of Mach number have been performed in this paper. The large eddy simulation (LES) turbulent model is employed to capture more details of turbulence and flow features in this study. The standard Smagorinsky sub-grid scale model is utilized to model the unresolved turbulent behaviors at sub-grid scale. Both the instantaneous and time-averaged flow fields are analyzed as well as turbulent statistics to compare with the experimental data sets. The simulation results are in good agreement with the experimental data sets available in the literature. The shear layer instability is observed in the instantaneous flow field which leads to the end of the potential core and then the centerline velocity begins to decay. The size of the potential core and the decay rate of the centerline velocity are consistent well with the experiments. The turbulence intensity in the fully developed region agree well with the experimental data and the velocity power spectrum reveals the -5/3 law. Simulation results show that the compressibility has a significant effect on a free shear layer growing rate which is consistent with the experiments. The complex shock wave structure is also captured exactly by the second-order total variation diminishing (TVD) scheme. GASFLOW-MPI code has been successfully validated using the experimental data of turbulent jet flow at the wide regime of Mach numbers.