A computationally efficient technique for solving two point boundary value problems in porous media

A numerical technique for solving two point boundary value problems in the porous media is presented. The technique is the combination of orthogonal collocation and finite element methods. One provides the accuracy whereas the other provides the stability to the numerical results. The technique is much simpler and the computer codes can be easily developed. The numerical results are compared with the analytic results for verification. The numerical results are found to match the analytic ones for more than four decimal places. The method can be applied to the stiff problems without any limit on the range of parameters.