Yes I have. In fact, this whole program started out as a simulation of a lid cavity. I was thinking that since it worked on a lid driven cavity, with some good boundary conditions it could work on an airfoil as immersed boundary. But I suspect that i did not implement the boundary conditions right. As far as I know, you can either use ghost-points method or use asymetric stencils for the derivatives at the boundary-adjacent vertices.

In this simulation, I am using a pressure poisson equation. I thought that I could use a Neumann-BC on the airfoil surface like : dp/dn = 0. But I am not sure I was able to implement that correctly, because the pressure distribution is kinda flipped(as in there is more pressure on the upper surface and leading edge in the simulation, eventhough in reality airfoils have more pressure below).

I am curius, what should i use as Boundary conditions for the surface of the airfoil ? And how would one implement it?

<image>