We present a discrete-domain approach to three-phase displacements and hysteresis in porous media. In this method, constrained energy minimization leads to evolution equations for local saturations that describe a wide ...