Abstract: Four kinds of two-equation turbulence model (k-ε and k-ω model) with compressibility modification,including the dilation-dissipation term and pressure-dilation term respectively,coupling with a full Reynolds-averaged Navier-Stokes code were used to simulate numerically the shock/turbulent boundary layer interactions,which occur in three different two-dimensional transonic nozzles with different Mach numbers in front of the shock.The computational results of both the time-averaged parameters,such as the velocity (u,the velocity along the x-coordinate) and the pressure (P,along the lower wall surface),and the pulsating parameters like the turbulent kinetic energy (k) were compared with the experimental results.The computational density contour of large separation case was compared with the experimental density contour.It is shown that when the two-equation turbulence models with compressibility modification are used to simulate the shock/turbulent boundary layer interaction,they predicate not only the distribution of the time-averaged parameters and the pulsating parameters but also the basic feature of the flow with or without separation.