A multi-fluid Eulerian model incorporating the kinetic theory of granular flow is used for the simulation of bubbling fluidized beds containing a binary mixture of Geldart B particles at low gas velocities. The cases of density, size and combined density/size segregation are investigated using computational fluid dynamic simulations. Various expressions for the drag force are evaluated for predicting different segregations. The simulation results show that summation of the particle-particle drag force, i.e., the "hindrance effect" term, and the Stokes drag of particles, which is modified based on the Wen-Yu drag model can be used for accurate simulation of a binary mixture of particles differing in size, density, or both. Bed expansion and dimensionless axial segregation profiles of CFD results are compared with the experimental data and good agreement is found.