Algebraic Dynamic Multilevel (ADM) Method For Simulations Of Multiphase Flow With An Adaptive Saturation Interpolator