Problems with the same model
 

The same model but with different boundary conditions is running. Boundaries had been calculated with 1D code for two different rotating velocities but, I emphasize, for the same model. One of two models is running with no problems. The other, with the smallest omega, crashed because of negative densities. I'm sure that boundaries are ok and I tryed to mesh with smaller cells wher I've problems but with no positive results. Have you got anything to sudgest me, please?