The aeroelasticity and the flutter of the wind turbine blade have been emphasized by related fields. The flutter of the wind turbine blade airfoil and its condition will be focused on. The eigenvalue method and the time domain analysis method will be used to solve the flutter of the wind turbine blade airfoil respectively. The flutter problem will be firstly solved using eigenvalue approach. The flutter region, where the flutter will occur and anti-flutter region, where the flutter will not occur, will be obtained directly by judging the sign of the real part of the characteristic roots of the blade system. Then the time domain analysis of flutter of wind turbine blade will be carried out through the use of the four-order Runge-Kutta numerical methods, the flutter region and the anti-flutter region will be gotten in another way. The time domain analysis can give the changing treads of the aeroelastic responses in great detail than those of the eigenvalue method. The flap displacement of wind turbine blade airfoil will change from convergence to divergence, and change from divergence to convergence extremely suddenly. During the flutter region, the flutter of wind turbine blade will occur extremely dramatically. The flutter region provided by the time domain analysis of the flutter of the blade airfoil accurately coincides with the results of eigenvalue approach, therefore the simulation results are reliable and credible.