# Inviscid Incompressible Flow in FreeFEM

**URL:** <https://community.freefem.org/t/inviscid-incompressible-flow-in-freefem/3284>\
**Category:** General Discussion\
**Created:** [June 3, 2024, 5:28pm UTC](https://community.freefem.org/t/inviscid-incompressible-flow-in-freefem/3284 "2024-06-03T17:28:10Z")\
**Posts on this page:** 1\
**Showing post:** 15

<div class="post-metadata">

**Author:** ![cmd](https://yyz2.discourse-cdn.com/flex030/user_avatar/community.freefem.org/cmd/32/67_2.png) [@cmd](https://community.freefem.org/u/cmd)\
**Post date:** [June 4, 2024, 10:56am UTC](https://community.freefem.org/t/inviscid-incompressible-flow-in-freefem/3284/15 "2024-06-04T10:56:38Z")

</div>

See the notes of @fb77 here:

> [@Boundary condition](https://community.freefem.org/t/boundary-condition/2945/6):
>
> If you want to avoid using a Lagrange multiplier (because it enlarges the matrix to invert), some other ways to penalize the normal boundary conditions are possible. In 3D you can use for example the boundary vertices. You have to define first a continuous approximate normal on the boundary (label 1 is for boundary) fespace Vh1(Th,[P2,P2,P2]); varf bdryn([w1,w2,w3],[v1,v2,v3])=int2d(Th,1)(v1\*N.x+v2\*N.y+v3\*N.z); Vh1 [normalappx,normalappy,normalappz]; normalappx[]=bdryn(0,Vh1); Then you can ad…

---

_[View the full topic](https://community.freefem.org/t/inviscid-incompressible-flow-in-freefem/3284)._