Abstract: The effect of different definitions of turbulence production terms in SST turbulence model on the transonic flow simulation is investigated numerically through solving the Reynolds-Averaged Navier-Stokes (RANS) equations with the use of fifth-order weighted compact nonlinear scheme (WCNS) and multi-block structured grid technology. The main purpose of present work is to provide technical support for engineering applications of high-order difference schemes on complex configurations. RAE2822 subcritical airfoil and DLR-F6 wing body configuration are adopted as the computational models in the paper. The numerical investigation focuses mainly on the influence of different definitions of turbulence production terms in SST turbulence model on residual convergence process, turbulent viscosity and velocity distribution in boundary layer, pressure coefficients and aerodynamic characteristics. The flow solutions with different definitions of turbulence production terms are also compared with some wind tunnel test results. The investigation shows that the influence of different definitions of turbulence production terms on high-order simulation is so trivial that no further consideration is needed for transonic flows without large separation bubble under a small angle of attack.