Dear all,
I implemented a discontinous Galerkin scheme to solve the stokes equations. The case P2dc-P1dc for velocity-pressure is good but the case P1dc-P0 does not work and I do not why.
Does anyone have an idea what the problem is ?
Dear all,
I implemented a discontinous Galerkin scheme to solve the stokes equations. The case P2dc-P1dc for velocity-pressure is good but the case P1dc-P0 does not work and I do not why.
Does anyone have an idea what the problem is ?
I can look at it if you send a code.
Thank you for you reply. Please find enclosed the code.
stokesdg.edp (1.8 KB)
Your right-hand side f
was not well computed. Computing it by using dx()
and dy()
in the fespaces induces some interpolation errors and inconsistencies. In particular dx(pex)
and dy(pex)
with pex
in P0
gives simply zero.
I have defined f explicitly and it seems to work.
stokesdg.edp (2.2 KB)
Yes, thank you for your replay !