The two fluids Eulerian model based on the kinetic theory of granular flow has been applied to predict the ice slurry flow in a vertical pipe without considering ice melting process. It can be seen that the ice particle velocity profile is symmetrical along the central axis in the vertical pipe. When ice slurry moves at low speed

the velocity gradient becomes small

and the concentration profile approaches even in the pipe cross section. When the moving speed increases

the ice particle velocity and concentration stratify. The ice particle concentration is relatively small near the pipe wall

but larger close to the central axis with a concentration peak nearby. The flow direction imposes minimal effects to velocity profiles but makes a great impact on particle concentration distribution. In addition

there are meaningful differences in pressure drop when flowing with different directions.