SUMMARYMass transfer in the presence of chemical reactions for flows through porous media is of interest to many disciplines. The Lattice Boltzmann method (LBM) is particularly attractive in such cases due to the ease with which it handles complicated boundary conditions. However, useful Lagrangian information (such as solute survival distance, effective diffusivity, collision frequency) is challenging to obtain from the LBM. In this paper, we present a straightforward and efficient Lagrangian methodology (Lagrangian scalar tracking, LST) for performing solute transport simulations in the presence of heterogeneous, first-order, irreversible reactions, based on a velocity field obtained from LBM. The hybrid LST/LBM technique tracks passive mass markers that have two contributions to their movement: convective (obtained through interpolation of a previously obtained velocity field) and Brownian. Various Schmidt number solutes and different solute release modes can be modeled with a single solvent flow field using this method. Moreover, the mass markers can have a range of reaction rate coefficients. This allows for the exploration of the whole spectrum of first-order heterogeneous reaction rates with just a single simulation. In order to show the applicability of the LST/LBM scheme, results from a case study are presented in which the consumption of oxygen and/or nutrients within a porous bone tissue engineering scaffold is modeled under flow perfusion culturing conditions. Although the reactive LST methodology described in this paper compliments the LBM, it can also be used with any other flow simulation that can generate the velocity field.