A CFD-DEM (computational fluid dynamics discrete element method) model, with heat transfer, and homogeneous and heterogeneous chemical reactions accounted for, is implemented here to numerically study the biomass gasification in a three-dimensional bubbling fluidized bed. The results show that the biomass inlet on the right sidewall results in laterally asymmetrical profiles of important bed hydrodynamics and gas concentrations. Increasing the temperature, steam-biomass ratio and equivalence ratio all enlarge the gas voidage, fluid force, collision force and dispersion coefficient of the solid phase, but have different effects on the different homogeneous reaction rates. Vigorous horizontal and vertical solid dispersions occur in the lateral central and sidewalls of the bed, respectively. The results of the current work provide valuable insights regarding the hydrodynamics of the bubbling fluidized bed for biomass gasification, and are expected to be useful for the operation, scale-up and optimization of such systems for sustainable energy production. (C) 2019 Elsevier Ltd. All rights reserved.