Quasi-implicit treatment of velocity-dependent mobilities in underground porous media flow simulation