Collocation with ‘adaptive’ finite elements in Buckley-Leverett problem

In solving the Buckley-Leverett problem, the use of higher order Hermitian polynomials at the frontal location and lower order ones elsewhere in the solution domain yields an enhanced frontal definition. This improvement is achieved without a significant increase in computational effort.