# DG Pressure Correction for Navier Stokes

**URL:** <https://community.freefem.org/t/dg-pressure-correction-for-navier-stokes/3386>\
**Category:** General Discussion\
**Created:** [July 5, 2024, 1:36pm UTC](https://community.freefem.org/t/dg-pressure-correction-for-navier-stokes/3386 "2024-07-05T13:36:42Z")\
**Posts on this page:** 1\
**Page:** 1

<div class="post-metadata">

**Author:** ![Ravi](https://yyz2.discourse-cdn.com/flex030/user_avatar/community.freefem.org/ravi/32/2468_2.png) [@Ravi](https://community.freefem.org/u/Ravi)\
**Post date:** [July 5, 2024, 1:36pm UTC](https://community.freefem.org/t/dg-pressure-correction-for-navier-stokes/3386/1 "2024-07-05T13:36:42Z")

</div>

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.
