The lecture will present the application of a locally-conservative Eulerian- Lagrangian method that has been discussed by Dr. Son-Young Yi for semi-linear parabolic equations to multiphase flow in porous media. After a brief reminder of her presentation and of the history of the derivation of the method, I shall apply a version of the method to two-phase flow in fractally inhomogeneous porous media. Several experimental results will be shown to indicate the convergence of the method under mesh refinement and its overall validity.