I am trying to solve the stabilized stokes equation using the DG method with P2dc-P1dc elements to give an exact solution and verify the order of spatial convergence. This is my code, and the code run shows that the convergence order is guaranteed to be correct after adding the stabilization term.
macro interstab(p, pp)(jump(p)*jump(pp)) //
+ intalledges(Th)( (1-nTonEdge)*interstab(p, pp) )
However I also tried programming in Matlab and found that the optimal spatial convergence order can be achieved without the stabilization term.
steady_stokes.edp (4.5 KB)