[jhabriel]

Hello, I'm having some problems validating a numerical simulation (ImPES [Implicit Pressure Explicit Saturation]) with the Buckley-Leverett (analytical) model. As you know, the Buckley-Leverett model only needs a constant inlet flow of a phase, on the other hand the ImPES method only deals with pressures values. How should I treat the boundary conditions in the ImPES method in order to compare the both cases?. I tried setting a fixed total Darcy velocity, but, in order to conservate the mass, in theory, the total velocity through each faces should be the same, so I would have the same velocity on all the time levels, and the saturation front won't be able to propagate.

Thanks for any idea.
Jhabriel Varela
National University of Asuncion
jhabriel@gmail.com