In this paper, an efficient ghost-cell based immersed boundary method (IBM) is used to perform direct numerical simulation (DNS) of reactive fluid-particle systems. With an exothermic first order reaction proceeding at the exterior particle surface, the solid temperature rises and consequently increases the reaction rate via an Arrhenius temperature dependence. In other words, the heat and mass transport is dependently coupled through the particle thermal energy equation and the Arrhenius equation, and they offer dynamic boundary conditions for the fluid phase thermal energy equation and species equation respectively. The fluid-solid coupling is enforced at the exact position of the particle surface by implicit incorporation of the boundary conditions into the discretized momentum, species and thermal energy conservation equations of the fluid phase. Different fluid-particle systems are studied with increasing complexity: a single sphere, three spheres and a dense array consisting of hundreds of randomly generated particles. In these systems the mutual impacts between heat and mass transport processes are investigated. (C) 2019 Elsevier Ltd. All rights reserved.