This paper considers similarity solutions of the multi-dimensional transport equation for the unsteady flow of two viscous incompressible fluids. We show that in plane, cylindrical and spherical geometries, the flow equation can be reduced to a weakly-coupled system of two first-order nonlinear ordinary differential equations. This occurs when the two phase diffusivity £>(#) satisfies {D/D')' = I/or and the fractional flow function / (6) satisfies df/d6 = KD"
11, where n is a geometry index (1, 2 or 3), a and K are constants and primes denote differentiation with respect to the water content 6. Solutions are obtained for time dependent flux boundary conditions. Unlike single-phase flow, for two-phase flow with n = 2 or 3, a saturated zone around the injection point will only occur provided the two conditions / 0 ' £)/(l -f)d9 < oo and/'(I) ^ 0 are satisfied. The latter condition is important due to the prevalence of functional forms of / (6) in oil/water flow literature having the property/'(I) = 0.