A numerical model for density-and-viscosity-dependent flows in two-dimensional variably saturated porous media