Reactive particulate systems are of prime importance in varieties of practical applications in process engineering. As an example this study considers extraction of phosphorous from waste water by calcium silicate hydrate particles in the P-RoC process. For such systems modeling has a large potential to help to optimize process conditions, e.g., particle-size distributions or volume flows. The goal of this study is to present a new generic modeling framework to capture relevant aspects of reactive particle fluid flows using combined lattice Boltzmann method and discrete-element method. The model developed is Euler-Lagrange scheme which consist of three-components viz., a fluid phase, a dissolved reactive substance, and suspended particles. The fluid flow and reactive mass transport are described in a continuum framework using volume-averaged Navier-Stokes and volume-averaged advection-diffusion-reaction equations, respectively, and solved using lattice Boltzmann methods. The volume-averaging procedure ensures correctness in coupling between fluid flow, reactive mass transport, and particle motion. The developed model is validated through series of well-defined benchmarks. The benchmarks include the validation of the model with experimental data for the settling of a single particle in a cavity filled with water. The benchmark to validate the multi-scale reactive transport involves comparing the results with a resolved numerical simulation. These benchmarks also prove that the proposed model is grid convergent which has previously not been established for such coupled models. Finally, we demonstrate the applicability of our model by simulating a suspension of multiple particles in fluid with a dissolved reactive substance. Comparison of this coupled model is made with a one-way coupled simulation where the influence of particles on the fluid flow and the reactive solution transport is not considered. This elucidates the need for the two-way coupled model.