In this paper, we present an inclusive tracking algorithm for water droplets in a wet steam flow through a multi-stage steam turbine. This algorism is based on the Eulerian-Lagrangian coupled solver. The solver continuously computes water droplet growth, kinematic non-equilibrium between vapor and droplets, capture and kinetics of droplets on turbine blades, departure of large droplets from the trailing edge of blades, acceleration and atomization of large droplets, and recollisions between blades and droplets. Our Eulerian-Lagrangian coupled solver is used to predict wetness in unsteady three-dimensional (3D) wet steam flows through three-stage stator rotor cascade channels in a low pressure (LP) steam turbine model which is developed by Mitsubishi Heavy Industries (MHI). Droplet groups tracked by the discrete droplet model (DDM) are placed in the computational domain according to the predicted wetness. Interference from the gas phase on the droplets is considered, to track their kinetic and behavior, until they reach the outlet of the computational domain. The aim of this research is to investigate those multi-physics phenomena that trigger all forms of loss in steam turbines. In addition, this method will also be applied to multi-physics problems such as erosion in future work. This paper is presented as a first step in the research. Overviews of model of current coupling solver and several test calculations are presented.