Preconditioned GMRES method is widely used in solving linear systems in the coupled multi-physics problems. Block triangular preconditioners usually give satisfactory performance. We first analyze such preconditioned GMRES method for an abstract problem, and then apply this technique to linear systems arising from the finite element discretization.