I am studying eulerian granular mullti-phase flow in 2-D fluidized bed. I have always divergence or floating point exception error after a few iterations. I tried different under relaxation factor and time step. Small time step is better for converging but it took too much time. Mesh qualities are acceptable limits for skewness and orthogonal quality. I will be appreciated for any advice for my problem.