Bubbling fluidized beds find application mainly in power conversion industries. For design, dimensioning, and operation of fluidized bed equipment, the understanding of multiphase gas-solid flows is of great importance. The use of computational fluid dynamics in the simulation of gas-solid systems is limited by the complexity of mathematical models, which rely on a series of empirical or theoretical correlations. In the present work, the code Multiphase Flow with Interphase eXchanges (MFIX) was employed to simulate flows in a bubbling fluidized bed and to compare results predicted using different gas-solid drag models. A two-fluid model with kinetic theory of granular flows (TFM-KTGF) was employed, in which gas-solid drag correlations, such as Gidaspow, Hill-Koch-Ladd, or Syamlal and O'Brien, were applied to model momentum transfer between phases. The results predicted were compared with each other and with experimental results from the literature. It was found that the results predicted using each model differ much. The Gidaspow and Hill-Koch-Ladd models yielded bubbles with shapes more similar to the experiments.