A Computational Study of the Influence of Drag Models and Heat Transfer Correlations on the Simulations of Reactive Polydisperse Flows in Bubbling Fluidized Beds

In this work, the influence of gas–solid drag and heat transfer coefficient models on the prediction capacity of the Euler–Euler approach to simulate reactive bubbling fluidized bed flows is studied. Three different cases are considered, a non-reactive bidisperse bubbling fluidized bed flow (Case 1), and two reactive polydisperse flows in bubbling fluidized beds, one for biomass gasification (Case 2), and the other for biomass pyrolysis (Case 3). The Gidaspow, Syamlal–O’Brien, and BVK gas–solid drag models and the Gunn, Ranz–Marshall, and Li–Mason gas–solid heat transfer correlations are investigated. A Eulerian multiphase approach in a two-dimensional Cartesian domain is employed for the simulations. Computational results for the three cases are compared with experimental data from the literature. The results obtained here contribute to a better understanding of the impacts of such closure models on the prediction ability of the Euler–Euler approach to simulate reactive flows. The results indicate that, for the simulation of reactive flows in bubbling fluidized bed reactors, the kinetic modeling of the reactions has a global effect, which superposes with the influence of the drag and heat transfer coefficient models. Nevertheless, local parameters can be noticeably affected by the choice of the interface closure models. Finally, this work also identifies the models that lead to the best results for the cases analyzed here, and thus proposes the use of such selected models for gasification and pyrolysis processes occurring in bubbling fluidized bed reactors.