Within a current research project, an Euler-Euler approach is used in a commercial CFD simulation program to describe an ice slurry flow in plate heat exchangers. The Eulerian model is based on the kinetic theory of granular flow. The ice flow patterns in different structures of channels of plate heat exchangers are investigated. In addition, it provides predictions for the pressure drop and the concentration profiles that represent the beginning of agglomeration. The ice agglomeration varies with concentration and particle size. The CFD simulation provides information about these parameters in combination with the channel geometry and the flow path. The simulation results show a comparison of pressure drop and volume fraction of the ice slurry flow in commercially available plate heat exchangers. The variation of the ice slurry parameters (concentration, particle size) at different locations in the plate heat exchangers (including gravity) is very crucial for the explanation of the surface shear stress and the particle interaction. The numerical investigations have focused mainly on illustrating the shear stress profiles, and mean velocity profiles in the channels by increasing the ice volume fraction. The numerical results provide force coefficients at different inlet velocities, ice concentrations and particle sizes. The wall shear stress profiles show the effect of gravity on the flow profile of ice slurry.