The methods of multiphase flow dynamics modelling are known. The theoretical calculation of the flow dynamic parameters of multiphase medium is based on the assumption of the mutual interaction between the phases, which becomes essentially manifested if mass concentration of the dispersed phase exceeds the value of 20 %. In many cases this limit seems to be even too high. The results of the experimental investigation of moderate and dense suspension flow confirm that the influence of solid particles on the suspension flow behaviour and the rheological parameters can be significant even for the mass concentration less than 20 %. In this paper the model of multiphase flow, which in based on the joint solution of vector equations of dispesed particle motion in Lagrange coordinates and the equation of continuum in Euilerian coordinates is presented. The method is applied to axially symmetrical flow of incompressible fluid in the rectangular channel and to flow in the zigzag channel. In the case of zigzag channel the numerical scheme is proposed. and Obsahuje seznam literatury