The paper is about the development of numerical methods for the simulation of multiphase flow using the Eulerian approach. Different discretizations for the convective terms and solution methods for the source terms are considered. Also, a splitting concept was developed to deal with particle reflection at a rigid wall.