Abstract
In analyzing infiltration into ground by using numerical models, favorable models should be characterized as a method with good convergence and good mass-balance as well as no oscillation at seepage fronts. This study compares the results obtained by finite element models for pressure-based equations and mixed-type, pressure-and-saturation-based equations with and without saturation parameterization to solve water-NAPL, two-phase flow equations with a special focus on fluid mass-balance and oscillation behavior near seepage fronts. Investigation of the accuracy of numerical results for both types of equations showed that the mixed-type, pressure-and-saturation-based equations give reliable results compared to the pressure-based equation, and the mixed-type equations with saturation parameterization give better results than those without saturation parameterization.
