Abstract: A computer-implemented simulation method of predicting local concentrations of constituents in coolant water anywhere along fuel rods within any fuel assembly mechanical design of a Boiling Water Reactor (BWR) potentially resulting in crud deposits on said fuel rods. The method is based on a sub-channel approach of predicting local mass fluxes of vapor and liquid in coolant water anywhere along fuel rods within any fuel assembly mechanical design of a Boiling Water Reactor (BWR) for given steady-state or transient boundary conditions. The sub-channel approach is based on the solution of mass, momentum and energy conservation equations for the vapor phase and the liquid phase, the liquid phase is represented by more than one field variable, and is specifically represented by three fields, with the vapor phase as a fourth field, consisting of droplets, a liquid base film, and disturbance waves.