In the end you need to solve a matrix vector system for ##x## of the form $$A \cdot x=b$$ for the unknown velocities, pressures etc. at every node in your mesh. At every iteration you have an approximation for ##x##, say ##\tilde x##. This leads to a nonzero vector $$R=b - a\cdot \tilde x $$. The problem is solved if ##R=0##. In the mean time, you can square all elements in ##R##, sum them and then take the square root to get the norm of the residual, which gives an indication for you how far you are from the actual solution and if you are making progress.