I am struggling with DG methods for 2 phase flow through porous media. I managed to get the global pressure, total flux equations to work with an unconditionally stable mixed FE DG formulation as described here: http://math.aalto.fi/reports/a531.ps
but the saturation equation is giving me problems. Basically I tried LDG and IP and both gave pretty weird unphysical results.
Does anyone have any good references for DG with flow through porous media? I found at least one good compositional flow paper, but not much about blackoil models.
Thanks in advance.
PS: I probably should add that I am writing the code with FEniCS and CBCPDESys. So far it has saved me years in development time.