[gerardosrez]

Hi guys
Im developing an unstructured CVFEM code for the NS equations using SIVA algorithm proposed by Baliga.

Im solving the induced flow in a square cavity for Re = 400, the mesh elements are triangular. For the first run my mesh had 100 nodes and I got convergence but when I refined the mesh I started to had divergence.

Im using MATLAB and I have tried with Gauss -Seidel, Bicgstab, Bicgstab with preconditioning without results, so what can you tell me about it

The code is using 1st order upwind and Im following the ideas of the following paper: Alexandre Lamoureux & B. Rabi Baliga (2011): Improved Formulations of the Discretized Pressure Equation and Boundary Treatments in Co-Located Equal-Order Control-Volume Finite-Element Methods for Incompressible Fluid Flow, Numerical Heat Transfer, Part B: Fundamentals: An International Journal of Computation and Methodology, 59:6, 442-472

Any help would be great, this is my master thesis and Im out of time so thanks in advance....