Numerical modeling of convection of powder was done. Convection of powder occurs when a strong vertical vibration is applied. When the control parameter, the acceleration amplitude Γ, exceeds some critical value Γc, convection and heaping start. In our simulation, the powder is modelled with a discrete elementary method but without rotation and with isotropic linear elastic interaction. The convection is reproduced, and Γc turns out to be about the gravity acceleration g, which agrees with the experiment, an other remarkable phenomena, surface fluidization, is also reproduced, and the critical value Γc again agrees with the experiment. It is suggestde that this convection is caused by the instability of the powder induced by the elastic interaction.