The CFD/DEM approach was used for investigating the forces playing a role in a furnace raceway formation and stability. The configuration is an actual pilot-scale hot blast furnace filled only with coke particles. In such a system, the raceway was unstable, with successively a growing phase and a collapse. The CFD/DEM numerical simulations were coupled with a core-shrinking model to mimic coke particle combustion. However, the kinetic reactions and heat transfers were not numerically predicted. Instead, the characteristic combustion timescale of one coke particle was imposed, and the combustion zone was adjusted to match the global combustion measured in the pilot-scale experiment. The results showed that the standard contact model was not enough to resist the pressure exerted by the granular weight on the raceway. However, the addition of a cohesive force, through the Johnson-Kendall-Roberts (JKR) model, allowed the qualitative reproduction of the gas pressure fluctuations and the collapse cycles in accordance with the experiment. A sensitivity analysis of the flow rate showed that CFD/DEM is able to reproduce quantitatively the time between two collapses, as observed in the experiment. Predicted raceway size and shapes are also in agreement with the experimental observations in the range of investigated parameters.