A least-squares mixed formulation is developed for simulation of two-phase flow in porous media. Such problems arise in petroleum applications and ground-water flow. An adaptive strategy based on the element residual as an error indicator is developed in conjunction with unstructured remeshing and tested for the two-phase flow of oil and water. An element-by-element conjugate-gradient scheme (EBE-CG) is compared to a band solution algorithm.