Probability Density Function Approach For Modeling Multi-Phase Flow In Porous Media