Computational Fluid Dynamics–Discrete Element Method Numerical Investigation of Binary Particle Mixing in Gas–Solid Fluidized Bed with Different Drag Models

The fluidized bed is a critical reactor in the energy and chemical industries, where the mixing and agglomeration behaviors of binary particles significantly influence both the efficiency of reaction processes and the uniformity of final products. However, the selection of appropriate drag force models remains a subject of debate due to the variability in particle properties and operating conditions. In this study, we investigated the fluidization behavior of binary mixtures composed of two different sizes of Geldart-D particles within a fluidized bed, evaluating nine distinct drag force models, including Wen and Yu; Schiller and Naumann; Ergun; Gidaspow, Bezburuah, and Ding; Huilin and Gidaspow; De Felice; Syamlal and O’Brien; and Hill, Koch, and Ladd. We focused on four key parameters: particle mixing degree, migration characteristics, temperature variation, and mean pressure drop. Simulation results revealed that the choice of drag model markedly affected mixing behavior, migration dynamics, and temperature distribution. Notably, the Ergun; Gidaspow, Bezburuah, and Ding; and Hill, Koch, and Ladd models exhibited superior particle mixing uniformity. While the drag model had a relatively minor impact on particle temperature changes, its selection became critical in simulations requiring high-temperature precision. Regarding pressure drop, the Huilin and Gidaspow and Gidaspow, Bezburuah, and Ding models demonstrated smaller and more stable pressure drop fluctuations. These findings offer valuable theoretical insights into gas–solid two-phase flow under binary particle mixing and provide practical guidance for the design and operation of fluidized bed reactors.