In this work two-phase flows in porous media are simulated numerically with Discontinuous Galerkin methods. The three methods Symmetrical Interior Penalty Galerkin method (SIPG), Non-symmetrical Interior Penalty Galerkin method (NIPG) and the scheme from Oden, Babuška and Baumann (OBB) are considered. The terminology and the examples are taken from soil science. First the Richards equation is solved using these methods. Then a twophase flows problem in the saturation/pressure formation is solved with OBB and NIPG. The numerical methods are implemented using the software toolkit PDELab. They are tested with examples from other publications. Weighted averages for absolute and relative permeabilities are examined.