[mnroelofs]

I have been trying to implement the SIMPLE algorithm as described by Ferzinger and Peric in the book "Computational Methods for Fluid Dynamics".
I want to simulate thec2D airflow around an airfoil of which I can dynamically create the shape. I use a triangular mesh for the flow field and use collocated variables. To solve the linear systems of equations I implemented a SOR Gauss-Seidel algorithm (but actually underrelax).

I implemented the algorithm in a C# application and the good news is it runs. However, since I am getting wrong results, I want to pinpoint whether the solver diverges, oscillates or that I just have gotten the theory wrong and thus implemented it incorrectly.

So, here are two pictures in which I have written down my discretization for the momentum and pressure correction equations. Can you spot any mistakes? The code file where the matrices and RHS vectors are created is also included. I hope the code is self-explanatory, otherwise ask.

Furthermore, should the coefficient matrix for u and v momentum be the same?

Thanks in advance for your time.