Modeling the saturation process of flows through rigid porous media by the solution of a nonlinear hyperbolic system with one constrained unknown