Unsymmetric conjugate gradient methods and sparse direct methods in finite element flow simulation