This paper presents a comprehensive validation and verification study of turbomachinery Reynolds-averaged Navier-Stokes flow solvers on the transonic axial compressor TUDa-GLR-OpenStage. Two commercial solvers namely Ansys CFX and Numeca FineTurbo are adopted to provide the benchmark solutions, which can be used for verification of other RANS solvers in the future. Based on these solvers, five sets of grids, two advection schemes (i.e., central difference and second-order upwind), four turbulence models (i.e., SA, SA-RC, SST and EARSM) and two rotor-stator interface models (i.e., mixing plane and sliding plane) are investigated to quantify their effects on predicting the performance and the flow field of the compressor stage. Results show that the choices of grid density and turbulence model are most sensitive to the prediction, leading to 5% and 7% variation in compressor performance characteristics, respectively. Regarding the choice of grid density, a method to estimate the grid discretization error is demonstrated, which is transferrable to other cases. Regarding the choice of turbulence model, the EARSM model is found overall most accurate among the investigated models, and the limitations and deficiencies of the rest models are discussed in detail based on the analysis of the mean flow fields and the eddy viscosity fields. The grids and the major CFD results presented in this work are open-accessed to the community for further research. The results and discussions presented in this paper provide a useful reference for future practices of RANS simulations for compressors.