Fully discrete pseudo spectral numerical schemes to 2-D and 3-D incompressible Euler equation are considered in this talk. To ensure the numerical stability for this inviscid equation, an artificial viscosity is added with a required numerical accuracy. A local in time convergence analysis is established for smooth solution. A global in time energy stability is assured with a suitable choice of the artificial viscosity term. Some numerical simulation results are also presented.