In order to get the details of flow fields in steam turbines, three-dimensional turbulent flow calculations are useful. However in a design procedure, three-dimensional flow calculations are only possible in the last design stage, because they need in-depth boundary conditions of both geometries and flows. At such a late time in the procedure, it is difficult to go back and change main design parameters, such as flow area and stage load. Both three-dimensional flow patterns and non-equilibrium condensation caused by rapid expansions of steam have important roles with respect to steam turbine performance particularly in low-pressure sections. The steam turbine internal efficiency can be improved by taking account of these effects in the early design stage, especially in flow pattern design. This paper describes a multi-stage through-flow calculation technique including both three-dimensional flow efffects and phase changes from vapour to small droplets. To compute the high-speed two phase steam flow, a flux-splitting procedure including non-equilibrium homogeneously condensation is introduced. Three-dimensional blade forces are calculated by using angles of both blade camber and radial lean. The blade camber lines can be decided without in-depth blade geometries. Therefore this computational technique is applicable in the flow pattern design. The calculation results agree well with fully three-dimensional flow calculation and the calculation can predict supersaturating states and Wilson lines which are defined as the maximum supercooling.