Hello everyone, I am trying to apply DG for pressure correction scheme of navier stokes equation.

In the code the bilinear form,:-

problem U ([uu1,uu2], [v1,v2])= int2d(Th) (uu1*v1+uu2*v2) + int2d(Th) (phi*dt*(dx(v1) + dy(v2))) // velocity Correction

-intalledges(Th) (dt * mean(phi) * (jump(v1)*N.x+jump(v2)*N.y)/nTonEdge)
-int2d(Th)(u1 * v1+u2 * v2)
;
isn’t working. Here phi is a known value, calculated from a previous bilinear form. I think the problem occurs when i am trying to integrate the product "mean(phi)*( jump(v1) * N.x + jump(v2) * N.y)".

So is there any other way of doing this integration or the problem is in somewhere else.

Please enlighten me.