Application of the numerical density-enthalpy method to the multi-phase flow through a porous medium