Numerical simulations of a stationary bubbling fluidized bed (BFB), based on an Eulerian approach, are performed. The closure of the Eulerian model is treated with two-equation models. Emphasis is put on the importance of three dimensionality in the simulations. The numerical results are compared to local instantaneous pressure measurements and time-averaged measurements (bed height and probability density function (pdf) of the spatial distribution of particles). The results of the Simulations show that two-dimensional simulations should be used with caution and only for sensitivity analysis, whereas three-dimensional simulations are able to reproduce both the statics (bed height and pdf of the spatial distribution of particles) and the dynamics (power spectrum of pressure fluctuations) of the bed. The non-sphericity of the particles is accounted for in the models and its influence on the results is presented.