The Modified Galerkin Method (MGM) has been proposed as one of the most efficient methods for two-dimensional convection-diffusion equations. In the MGM, the non-symmetric matrices, which are derived from the convection term in the Galerkin formulation, are not used, and an artificial diffusion is introduced through an error analysis approach to improve its discretization accuracy in both time and space directions. In this study, the MGM is applied for two-dimensional viscous fluid flow analysis, and the driven cavity flow problems are solved up to Reynolds number of 10,000 using the vorticity-stream function formulation and non-uniform meshes. The results show the effectiveness of MGM.