Turbulent flow widely exists in the aerospace field, and it is still challenging to realise the accurate prediction in the numerical simulation. To realise the high-fidelity numerical simulation of compressible turbulent flow, a high-order accurate self-adaptive turbulence eddy simulation (SATES) method is developed on the PHengLEI-HyOrder open-source solver, combining with the high-order accurate weighted compact nonlinear schemes (WCNS). The compressible flow in the subsonic and transonic is numerically simulated, including some typical cases, such as subsonic flow past a circular cylinder and flow past a square cylinder, high-lift configuration DLR-F11, transonic flow around a circular cylinder. The results predicted by the current high-order accurate SATES are in good agreement with the available experimental and numerical data. The present numerical method can also accurately capture the interactions between shock waves and turbulence while accurately simulating flow separation, shear layer instability and large-scale vortex shedding. The results obtained show that the current high-order accurate SATES simulations based on PHengLEI-HyOrder solver can accurately simulate complex turbulent flows with high reliability.