The paper discusses a reuse of matrix factorization as building block in the Augmented Lagrangian (AL) and modified AL preconditioners for non-symmetric saddle point linear algebraic systems. strategy is applied to solve two-dimensional incompressible fluid problems with efficiency rates independent Reynolds number. solver then tested simulate motion surface fluid, an example 2D flow motivated ...