Modeling a sediment transport phenomena using Volume of Fluid (VOF) method for Multiphase problem and Eulerian-Larangian particle transport method for Discreet Phase Modeling (DPM). My mesh maximum skewness is around 0.5. For solver settings pressure-velocity coupling was introduced with Green-Gauss Cell Based gradient, PRESTO! pressure and Second Order Upwind momentum. I tightened all URF (Under-Relaxation Factors) to 0.1~0.2 and used 0.001 time step size.
After running the solution for some time everything looks fine in the animation but the solver fails showing "Error: Divergence detected in AMG solver: epsilon". Anyone has any idea on how it can be solved or the error can be removed?