We report the bubbling fluidized bed reactor is widely used in fast pyrolysis of biomass. Discrete simulation of this reactor is challenging due to many sand particles and lack of accurate drag corrections accounting for the interaction of two different solid particles with different properties. In this research, the computational cost is reduced by using the coarse-grained computational fluid dynamic-discrete element method, where many sand particles are lumped into a larger numerical parcel. The Syamlal?O'Brien drag model is used for sand, while Ganser correction coupled with Gidaspow model is used for the nonspherical biomass particles. This hybrid approach shows superior behavior over other drag models using pressure drops as a benchmark. The predicted bed height and pressure fluctuating frequencies compare well with experiment. Finally, the mixing of biomass is close to perfect if the superficial velocity is larger than four times the minimum fluidization velocity.