Why does the dispersion force become non-zero for volume fraction below 0.1 and over 0.9 ? It's also a very steep gradient, so maybe that's the cause of your problems: small change in volume fraction gives a very different force repsonse. With that you could trigger a feed back look that bounces the force between -10 and -60 for a negligble change in flow field.