Abstract: A new numerical method of integrating the nonlinear evolution equations, namely the Taylor expansion method, was presented. The standard Galerkin method can be viewed as the 0-th order Taylor expansion method; while the nonlinear Galerkin method can be viewed as the 1-st order modified Taylor expansion method. Moreover, the existence of the numerical solution and its convergence rate were proven. Finally, a concrete example, namely, the two-dimensional Navier-Stokes equations with a non slip boundary condition,was provided. The result is that the higher order Taylor expansion method is of the higher convergence rate under some assumptions about the regularity of the solution.

Key words: nonlinear evolution equation, Navier-Stokes equation, Taylor expansion method, convergence rate