Simulation of incompressible fluid flows using Galerkin weighted residual finite element method (C++)

To solve incompressible Navier-Stokes equations using Finite Element (FE) method, the simple lagging method is used to reform non-linear terms into linear ones. The pressure-velocity coupling problem has been also solved using inequal order method wherein velocities are approximated using second-degree polynomials, however, it would be first-degree for the pressure. This solver is written in C++. Matlab is also used to derive symbolic formulations needed in the FE discretization method. Then, Galerkin weighted residual method would be used in the context of the FE method. After applying boundary conditions, the system of linear equations has been solved using GMRES solver.