Uncertainty Assessment of Subsurface Flow by a Probability Density Function Method