Implements the Buckley-Leverett equation
The initial condition is
and boundary conditions are
Release Notes: