A numerical technique for the analysis of three-dimensional compressible turbulent flows in a turbine stage is presented .To calculate the steady interaction flow fields in a nozzle and bucket simultaneously, the nozzle outlet elements and the bucket inlet elements are overlapped in the axial direction and are used for connecting boundary elements. To calculated the flows in arbitrarily shaped geometries, a control volume method combined with a body-fitted curvilinear coordinate system is used to obtain spatially discretized governing equations. In the present analysis, a two-equation model of turbulence is introduced to estimate the turbulence effect. In order to assure the effectiveness of the present method, a computation is carried out for the flow in a model turbine stage. Experimental data is also obtained by the use of a 5-hole pitot tube for the purpose of comparison with computational results. It is shown that three-dimensional flow phenomena due to the viscous effect are well predicted and the comparisons with experimental data give encouraging results for three-dimensional flow prediction in a turbine stage.