I propose a numerical model which describes the dynamical features of vibrated beds. It succeeds in reproducing convective motion in vibrated beds which was observed by Faraday for the first time in 1831. In addition to this, this modeling can explain threshold value of instability and surface fluidization. Moreover, I numerically show vibrated bed without side wall can exhibit strong non-linear feature like turbulence or anomalous diffusion.