Lie symmetry method is applied to a system of nonlinear PDEs describing two-phase mass flows. Optimal systems of Lie subalgebras are constructed. Similarity forms are constructed using optimal systems of Lie subalgebras. Solutions are presented for symmetry reduced systems. Solutions describe interesting observable phenomena of two phase mass flow.